Reply to Comment on “Impact of groundwater withdrawals on the interaction of multi-layered aquifers in the Viterbo geothermal area (central Italy)”: report published in Hydrogeology Journal (2013) 21:1339–1353, by Antonella Baiocchi,Francesca Lotti and Vincenzo Piscopo

Hydrogeology Journal -     

Structural evolution of a 2M 1 phengite mica up to 11 GPa: an in situ single-crystal X-ray diffraction study

The structural evolution at high pressure of a natural 2M ₁-phengite [(K_0.98Na_0.02)_Σ=1.00(Al_1.55Mg_0.24Fe_0.21Ti_0.02)_Σ=2.01(Si_3.38Al_0.62)O₁₀(OH)₂; a = 5.228(2), b = 9.057(3), c = 19.971(6)Å, β = 95.76(2)°; space group: C2/c] from the metamorphic complex of Cima Pal (Sesia Zone, Western Alps, Italy) was studied by single-crystal X-ray diffraction with a diamond anvil cell under hydrostatic conditions up to ~11 GPa. A series of 12 structure refinements were performed at selected pressures within the P range investigated. The compressional behaviour of the same phengite sample was previously studied up to ~25 GPa by synchrotron X-ray powder diffraction, showing an irreversible transformation with a drastic decrease of the crystallinity at P > 15–17 GPa. The elastic behaviour between 0.0001 and 17 GPa was modelled by a third-order Birch–Murnaghan Equation of State (BM-EoS), yielding to K _T0 = 57.3(10) GPa and K′ = ?K _T0/?P = 6.97(24). The single-crystal structure refinements showed that the significant elastic anisotropy of the 2M ₁-phengite (with β(a):β(b):β(c) = 1:1.17:4.60) is mainly controlled by the anisotropic compression of the K-polyhedra. The evolution of the volume of the inter-layer K-polyhedron as a function of P shows a negative slope, Fitting the P–V(K-polyhedron) data with a truncated second-order BM-EoS we obtain a bulk modulus value of K _T0(K-polyhedron) = 26(1) GPa. Tetrahedra and octahedra are significantly stiffer than the K-polyhedron. Tetrahedra behave as quasi-rigid units within the P range investigated. In contrast, a monotonic decrease is observed for the octahedron volume, with K _T0 = 120(10) GPa derived by a BM-EoS. The anisotropic response to pressure of the K-polyhedron affects the P-induced deformation mechanism on the tetrahedral sheet, consisting in a cooperative rotation of the tetrahedra and producing a significant ditrigonalization of the six-membered rings. The volume of the K-polyhedron and the value of the ditrigonal rotation parameter (α) show a high negative correlation (about 93%), though a slight discontinuity is observed at P >8 GPa. α increases linearly with P up to 7–8 GPa (with ?α/?P ≈ 0.7°/GPa), whereas at higher Ps a “saturation plateau” is visible. A comparison between the main deformation mechanisms as a function of pressure observed in 2M ₁- and 3T-phengite is discussed.     

High-pressure polymorphism and structural transitions of norsethite,BaMg(CO3)2

In situ high-pressure investigations on norsethite, BaMg(CO₃)₂, have been performed in sequence of diamond-anvil cell experiments by means of single-crystal X-ray and synchrotron diffraction and Raman spectroscopy. Isothermal hydrostatic compression at room temperature yields a high-pressure phase transition at P _c ≈ 2.32 ± 0.04 GPa, which is weakly first order in character and reveals significant elastic softening of the high-pressure form of norsethite. X-ray structure determination reveals C2/c symmetry (Z = 4; a = 8.6522(14) Å, b = 4.9774(13) Å, c = 11.1542(9) Å, β = 104.928(8)°, V = 464.20(12) Å³ at 3.00 GPa), and the structure refinement (R ₁ = 0.0763) confirms a distorted, but topologically similar crystal structure of the so-called γ-norsethite, with Ba in 12-fold and Mg in octahedral coordination. The CO₃ groups were found to get tilted off the ab-plane direction by ~16.5°. Positional shifts, in particular of the Ba atoms and the three crystallographically independent oxygen sites, give a higher flexibility for atomic displacements, from which both the relatively higher compressibility and the remarkable softening originate. The corresponding bulk moduli are K ₀ = 66.2 ± 2.3 GPa and dK/dP = 2.0 ± 1.8 for α-norsethite and K ₀ = 41.9 ± 0.4 GPa and dK/dP = 6.1 ± 0.3 for γ-norsethite, displaying a pronounced directional anisotropy (α: β _a ^?1  = 444(53) GPa, β _c ^?1  = 76(2) GPa; γ: β _a ^?1  = 5.1(1.3) × 10³ GPa, β _b ^?1  = 193(6) GPa β _c ^?1  = 53.4(0.4) GPa). High-pressure Raman spectra show a significant splitting of several modes, which were used to identify the transformation in high-pressure high-temperature experiments in the range up to 4 GPa and 542 K. Based on the experimental series of data points determined by XRD and Raman measurements, the phase boundary of the α-to-γ-transition was determined with a Clausius–Clapeyron slope of 9.8(7) × 10^?3 GPa K^?1. An in situ measurement of the X-ray intensities was taken at 1.5 GPa and 411 K in order to identify the nature of the structural variation on increased temperatures corresponding to the previously reported transformation from α- to β-norsethite at 343 K and 1 bar. The investigations revealed, in contrast to all X-ray diffraction data recorded at 298 K, the disappearance of the superstructure reflections and the observed reflection conditions confirm the anticipated \(R\bar{3}m\) space-group symmetry. The same superstructure reflections, which disappear as temperature increases, were found to gain in intensity due to the positional shift of the Ba atoms in the γ-phase.     

The Cryogenic AntiCoincidence detector for ATHENA X-IFU: a scientific assessment of the observational capabilities in the hard X-ray band

ATHENA is a large X-ray observatory, planned to be launched by ESA in 2028 towards an L2 orbit. One of the two instruments of the payload is the X-IFU: a cryogenic spectrometer based on a large array of TES microcalorimeters, able to perform integral field spectrography in the 0.2–12 keV band (2.5 eV FWHM at 6 keV). The X-IFU sensitivity is highly degraded by the particle background expected in the L2 orbit, which is induced by primary protons of both galactic and solar origin, and mostly by secondary electrons. To reduce the particle background level and enable the mission science goals, the instrument incorporates a Cryogenic AntiCoincidence detector (CryoAC). It is a 4 pixel TES based detector, placed < 1 mm below the main array. In this paper we report a scientific assessment of the CryoAC observational capabilities in the hard X-ray band (E > 10 keV). The aim of the study has been to understand if the present detector design can be improved in order to enlarge the X-IFU scientific capability on an energy band wider than the TES array. This is beyond the CryoAC baseline, being this instrument aimed to operate as anticoincidence particle detector and not conceived to perform X-ray observations.     

The particle background of the <Emphasis Type="Italic">X-IFU</Emphasis> instrument

In this paper we are going to review the latest estimates for the particle background expected on the X-IFU instrument onboard of the ATHENA mission. The particle background is induced by two different particle populations: the so called “soft protons” and the Cosmic rays. The first component is composed of low energy particles (< 100s keV) that get funnelled by the mirrors towards the focal plane, losing part of their energy inside the filters and inducing background counts inside the instrument sensitivity band. The latter component is induced by high energy particles (> 100 MeV) that possess enough energy to cross the spacecraft and reach the detector from any direction, depositing a small fraction of their energy inside the instrument. Both these components are estimated using Monte Carlo simulations and the latest results are presented here.     

Micro-seismic monitoring after the shipwreck of the Costa Concordia at Giglio Island (Italy)

A micro-seismic network was used for monitoring the wreck of the Costa Concordia cruise ship, wrecked and run agrounded along the Giglio Island coasts during the night of 13 January 2012, until its removal. The seismic traces were processed by means of real-time and “a posteriori” procedures to detect transients that could be ascribed to wreck movements on the sea bed to integrate this information in an early warning system for assessing the wreck stability. After a first discrimination of the transients using amplitude criteria we proceeded to the localization of the detected signals to focus the attention only on the transients originated in the shipwreck resting area. The results showed that most of the events localized on the wreck were likely related to human work activities or sudden internal brittle failure but not to displacements on the seafloor. Instead, the displacements are associated to the impact on the vessel of great sea storms which approach were well correlated with the increasing seismic noise at low frequency. The carried out procedures based on this unique dataset represent an opportunity to test seismic monitoring techniques also in not usual engineering context to support emergency management activities.     

An XMM-Newton proton response matrix

Soft protons constitute an important source of background in focusing X-ray telescopes, as Chandra and XMM-Newton experience has shown. The optics in fact transmit them to the focal plane with efficiency similar to the X-ray photon one. This effect is a good opportunity to study the environment of the Earth magnetosphere crossed by the X-ray satellite orbits, provided that we can link the spectra detected by the instruments with the ones impacting on the optics. For X-ray photons this link has the form of the so-called response matrix that includes the optics effective area and the energy redistribution in the detectors. Here we present a first attempt to produce a proton response matrix exploiting ray-tracing and GEANT4 simulations with the final aim to be able to analyse XMM-Newton soft proton data and link them to the external environment. If the procedure is found to be reliable, it can be applied to any future X-ray missions to predict the soft particles spectra impacting on the focal plane instruments.     

Impact of groundwater withdrawals on the interaction of multi-layered aquifers in the Viterbo geothermal area (central Italy)

The impact of groundwater withdrawals on the interaction between multi-layered aquifers with different water qualities in the Viterbo geothermal area (central Italy) was studied. In this area, deep thermal waters are used to supply thermal spas and public pools. A shallow overlying aquifer carries cold and fresh water, used for irrigation and the local drinking-water supply. Starting with a conceptual hydrogeological model, two simplified numerical models were implemented: a steady-state flow model of the entire groundwater system, and a steady-state flow and heat transport model of a representative area, which included complex interactions between the aquifers. The impact of increased withdrawals associated with potential future development of the thermal aquifer must be considered in terms of the water temperature of the existing thermal sources. However, withdrawals from the shallow aquifer might also influence the discharge of thermal sources and quality of the water withdrawn from the shallow wells. The exploitation of the two aquifers is dependent on the hydraulic conductivity and thickness of the intervening aquitard, which maintains the delicate hydrogeological equilibrium. Effective methods to control this equilibrium include monitoring the vertical gradient between the two aquifers and the residual discharge of natural thermal springs.