Metabolic adaptations underpin high productivity rates in relict subsurface water.

Groundwater aquifers are ecological hotspots with diverse microbes essential for biogeochemical cycles. Their ecophysiology has seldom been studied on a basin scale. In particular, our knowledge of chemosynthesis in the deep aquifers where temperatures reach 60 °C, is limited. Here, we investigated the diversity, activity, and metabolic potential of microbial communities from nine wells reaching ancient groundwater beneath Israel's Negev Desert, spanning two significant, deep (up to 1.5 km) aquifers, the Judea Group carbonate and Kurnub Group Nubian sandstone that contain fresh to brackish, hypoxic to anoxic water. We estimated chemosynthetic productivity rates ranging from 0.55 ± 0.06 to 0.82 ± 0.07 µg C L-1 d-1 (mean ± SD), suggesting that aquifer productivity may be underestimated. We showed that 60% of MAGs harbored genes for autotrophic pathways, mainly the Calvin-Benson-Bassham cycle and the Wood-Ljungdahl pathway, indicating a substantial chemosynthetic capacity within these microbial communities. We emphasize the potential metabolic versatility in the deep subsurface, enabling efficient carbon and energy use. This study set a precedent for global aquifer exploration, like the Nubian Sandstone Aquifer System in the Arabian and Western Deserts, and reconsiders their role as carbon sinks.

Deep desert aquifers as an archive for Mid- to Late Pleistocene hydroclimate: An example from the southeastern Mediterranean.

Many efforts have been made to illuminate the nature of past hydroclimates in semi-arid and arid regions, where current and future shifts in water availability have enormous consequences on human subsistence. Deep desert aquifers, where groundwater is stored for prolonged periods, might serve as a direct record of major paleo-recharge events. To date, groundwater-based paleoclimate reconstructions have mainly focused on a relatively narrow timescale (up to ∼40 kyr), limited by the relatively short half-life of the widely used radiocarbon (5.73 kyr). Here we demonstrate the usage of deep regional aquifers in the arid southeastern Mediterranean as a hydroclimate archive for earlier Mid-to-Late Pleistocene epochs. State-of-the-art dating tools, primarily the 81Kr radioisotope (t1/2 = 229 kyr), were combined with other atmosphere-derived tracers to illuminate the impact of four distinguishable wetter episodes over the past 400 kyr, with differences in climatic conditions and paleo-recharge locations. Variations in stable water isotope composition suggest moisture transport from more proximal (Mediterranean) and distal (Atlantic) sources to different parts of the region at distinct times. Large variability in the computed noble gas-based recharge temperature (NGT), ranging ~15-30 °C, cannot be explained by climate variations solely, and points to different recharge pathways, including geothermal heating in the deep unsaturated zone and recharge from high-elevation (colder) regions. The obtained groundwater record complements and enhances the interpretation of other terrestrial archives in the arid region, including a contribution of valuable information regarding the moisture source origin as reflected in the deuterium-excess values, which is unattainable from the common practice analysis of calcitic cave deposits. We conclude that similar applications in other deep (hundred-m-order) regional groundwater systems (e.g., the Sahara desert aquifers) can significantly advance our understanding of long-term (up to 1 Myr) paleo-hydroclimate in arid regions, including places where no terrestrial remnants, such as cave, lake, and spring sediments, are available.

Controls on the 36Cl/Cl input ratio of paleo-groundwater in arid environments: New evidence from 81Kr/Kr data.

Measurements of the long-lived 81Kr and 36Cl radioisotopes in groundwater from the Negev Desert (Israel) were used to assess the 36Cl/Cl input ratios and Cl- contents for paleorecharge into the Nubian Sandstone Aquifer (NSA). The reconstructed Cl- content of the recharge flux was on the order of 300-400 mg/L. An initial 36Cl/Cl ratio of 50 × 10-15 was assessed for the groundwater replenishment in the Negev Desert since the late Pleistocene, in agreement with the 36Cl/Cl ratios in recent local rainwater. This is despite possible changes in the climatic conditions and the 36Cl production rates in the atmosphere over this timeframe. This similarity in values is explained by the major role played by the erosion and weathering of near-surface materials in the desert environment that dominate the hydrochemistry of rains, floods, and the consequent groundwater recharge. Spatial variation in the reconstructed initial 36Cl/Cl ratio is accounted for by the differences in the mineral aerosol sources for specific recharge areas of the NSA. Accordingly, regional variations in the 36Cl/Cl input in groundwater reservoirs surrounding the Mediterranean Sea indicate various processes that govern the 36Cl/Cl system. Finally, the results of this study highlight the great advantage of integrating 81Kr age information in evaluating the initial 36Cl/Cl and Cl- input, which is essential for the calibration of 36Cl radioisotope as an available long-term dating tool for a given basin.

Correlating fracture trends and hydraulic head using semivariogram cloud analysis.

Fracture trends (defined as kilometer-scale linear features interpolated between field observations of fractures along their strikes) often have a dominant orientation. Finding a correlation between this orientation and hydraulic data could shed light on their hydraulic influence. A significant correlation between head residuals from first-order regional drift and the orientation of 2- to 4-km-long fracture trends was found in a study site in the Negev, Israel, using the semivariogram cloud analysis. Correlation of head residuals rather than the head itself implies that the orientation of the fracture trends controls the anisotropy and heterogeneity at this scale, mainly because the fracture trends define the orientation of blocks, which differ in their hydraulic properties. Preferential transmissive pathways are probably shorter than the full extent of the fracture trends, causing a relatively high head difference along the trends on the 2- to 4-km scale. Fracture trend density and additional data from short-range hydraulic tests helped characterize two blocks separated by a fault zone. The identification of hydraulic features on a kilometer scale is necessary for better modeling of regional ground water flow and transport. Hydraulic tests at this scale are not feasible, thereby rendering combined analyses of head and structural data, such as the one presented here, essential.

Particle transport in unsaturated fractured chalk under arid conditions.

A series of field and laboratory experiments were conducted to study the mechanisms of particle detachment and transport from fractures in vadose chalk. Experiments of intermittent flow events along fracture surfaces were carried out in the laboratory. In the field, water was percolated from land surface via a discrete fracture into a compartmental sampler installed inside a horizontal corehole located I m below the surface. The mass, size distribution, and composition of the particles drained from the fracture voids were examined along with flow rates and salt dissolution. Two boreholes penetrating the underlying saturated zone were sampled and analyzed for colloidal concentration and composition. Most of the particle and solute release at the drained effluents occurred during the first several hours of flow, but erratic pulses of particles were still observed after long periods of time. Most of the detached particles had a mean diameter of >2 microm, while the mobile colloidal phase in the groundwater had a mean diameter of approximately 1 microm. Mineralogical composition of the groundwater colloids and the particles detached from the upper vadose fracture were similar. Laboratory observations demonstrated the importance of the existence of a coating layer, made of weathered particles and salts, on particle detachment. The results of this study suggest that: (1) particle detachment causes flow-rate variability in the unsaturated fracture; (2) the mechanisms of particle detachment and salt dissolution within the fracture are linked: and (3) although most of the detached particles are large and likely to accumulate inside fractures, some colloidal particles also eroded from the fracture void and are likely to be transported to the groundwater.