Hydrogeologic Framework Research Articles

Geophysical study of the Salto–Arapey Aquifer System in Salto, Uruguay

The Salto–Arapey Aquifer System (SAAS) represents a complex hydrogeological framework composed of a detrital aquifer that lays over a fissured one in the Chaco-Paraná Basin. To enhance our understanding of this underground system, an initiative was undertaken to characterize a specific sector of the SAAS utilizing a geophysical methods known as Audiomagnetotelluric sounding. The primary outcomes of this endeavor include the geoelectrical characterizations of key formations and the development of a 2.5D geological-geophysical model. This endeavor significantly advances our geological and geophysical comprehension of the targeted sector, thereby offering substantial societal and economic benefits. Although the method did not allow differences between sedimentary formations, it did allow the recognition of altered basalts, which would facilitate the identification of areas with considerable potential for the productive development of water wells. Notably, this investigation holds significant local relevance as it presents a replicable methodology to characterize the rest of the system, or even aquifers sharing similar characteristics.

Hydrogeological conceptual model of Stampriet transboundary aquifer system in Southern Africa

Stampriet Transboundary Aquifer System (STAS) is shared between Botswana, Namibia and South Africa with massive utilization especially in Namibia. However, the understanding of the comprehensive hydrogeological framework of the aquifer system is limited. This study aimed to develop a comprehensive hydrogeological conceptual model of the aquifer system which mainly entailed 3D hydrostratigraphic modelling using Rockworks 3D geological modelling software, analysis of regional groundwater occurrence and flow system, sources and sinks, hydraulic parameters and definition of STAS boundaries. Six hydrostratigraphic units composed of four aquifers, including a new one (the Dwyka), and two aquitards were identified. Analysis of the regional groundwater flow system revealed that groundwater flows in two directions, indicating presence of a regional groundwater divide. 3D hydrostratigraphic geometry revealed that groundwater flow in STAS is influenced by topography but also by faults such as the Zoetfontein Fault, which displaces the pre-Kalahari rocks vertically creating a barrier to groundwater flow. Hydrostratigraphic and 3D geometry work also led to a 42% increase of the size of STAS compared to the boundary initially delineated by UNESCO. The proposed hydrogeological conceptual model forms a basis for development of numerical integrated hydrological model to simulate surface-groundwater interactions, regional groundwater flow and aquifer development potential.

Tracing Hydrological Processes: Insights from Hydrochemical and Isotopic Investigations in the Northern Part of Croatian Dinaric Karst

This study presents the findings derived from a three-year monitoring programme focusing on the hydrochemical composition and stable isotope signatures (δ18O, δ2H) within the vulnerable karst groundwater system in Croatia's northern Dinaric karst region. Covering an area of 1980 km2 across Kapela Mountain and its foothills, this groundwater system falls within the Kupa River basin and encompasses the catchment areas of the main springs across two spring levels of the Dobra, Mrežnica, and Slunjčica Rivers (namely the Vitunj, Zagorska Mrežnica, Dretulja, Veliko Vrelo, Gojak, Tounjčica, Mrežnica, and Slunjčica springs). Given the complex hydrogeological framework, prior studies have not extensively addressed the hydrochemical characteristics of this region, thereby necessitating a comprehensive investigation to elucidate system dynamics. This paper outlines the first investigation employing stable isotopic analyses within the karst groundwater system of Kapela Mt. The main results and conclusions of the research include: (i) the aquifers across Kapela Mt. drain mainly limestones, (ii) water–rock interaction and carbonate weathering are significant contributors to water geochemistry, (iii) there is a slight human impact on the Gojak and Tounjčica spring waters, (iv) the mean residence time of water in the observed aquifer is up to 1.5 years, (v) the new LMWL was calculated for Kapela Mt. based on a three-year dataset (June 2018 – May 2021), (vi) in the north Dinaric karst, the predominant origin of precipitation is from the Mediterranean air mass, (vii) Velebit Mt. has a strong influence on the precipitation isotopic composition of the study area.

Hydrogeologic Framework Model‐Based Numerical Simulation of Groundwater Flow and Salt Transport and Analytic Hierarchy Process‐Based Multi‐Criteria Evaluation of Optimal Pumping Location and Rate for Mitigation of Seawater Intrusion in a Complex Coastal Aquifer System

AbstractA series of hydrogeologic framework model (HFM)‐based steady‐ and transient‐state numerical simulations is performed first using a coupled subsurface flow‐transport numerical model to analyze groundwater flow and salt transport in an actual three‐dimensional complex coastal aquifer system before and during groundwater pumping. A series of analytic hierarchy process (AHP)‐based multi‐criteria evaluations is then performed applying a multi‐criteria decision‐making approach to determine optimal pumping location and rate for a new pumping well in the complex coastal aquifer system during groundwater pumping. The complex coastal aquifer system is composed of six anisotropic fractured porous geologic media (five rock formations and one fault) and three isotropic porous geologic media (three soil formations) and shows high geometric irregularity and significant heterogeneity and anisotropy of the nine geologic media. Results of the steady‐state numerical simulations show successful model calibration with 26 measured groundwater levels and two observed seawater intrusion front lines. The latter two are determined by spatial interpolation and extrapolation of electrical conductivity logging data and electrical resistivity survey data, respectively. Based on the status and prospect of necessary water uses and available groundwater resources, the field observations of groundwater and seawater intrusion, and the analyses of the steady‐state numerical simulation after the model calibration, six candidate pumping locations are selected for the new pumping well. In addition, from six preliminary individual transient‐state numerical simulations, maximum pumping rates at the six candidate pumping locations are calculated first, and a set of six incremental candidate pumping rates is then assigned at each of the six candidate pumping locations. Results of the transients‐state numerical simulations show that groundwater flow and salt transport are spatially and temporally changed, and seawater intrusion is further intensified by groundwater pumping. In addition, the magnitudes of such spatial and temporal changes and intensification are significantly different depending on the candidate pumping locations and rates. Results of the steady‐ and transient‐state numerical simulations also show that both complexity (geometric irregularity, heterogeneity, and anisotropy including the fault) and topography have significant effects on the spatial distributions and temporal changes of groundwater flow and salt transport in the coastal aquifer system before and during groundwater pumping. In addition, results of statistical estimations of the mesh Peclet and Courant numbers confirm acceptabilities of minimizing numerical dispersion in the steady‐ and transient‐state numerical simulations. Based on the analyses of the transient‐state numerical simulations, eight multiple criteria are chosen to judge, prioritize, and rank the six candidate pumping locations and six candidate pumping rates for optimal pumping. Results of the multi‐criteria evaluations determine the optimal pumping location and rate for the new pumping well among the six candidate pumping locations and six candidate pumping rates. In addition, results of consistency checks confirm consistencies of judgments in the multi‐criteria evaluations.

Using an Open-Source Tool to Develop a Three-Dimensional Hydrogeologic Framework of the Kobo Valley, Ethiopia

Groundwater resource management requires understanding the groundwater basin’s hydrogeology and would be improved with the development of a three-dimensional hydrogeologic framework model (HFM). A wide range of methods and software exist to quantify the extent, structure, and properties of geologic systems. However, most geologic software is proprietary and cost-prohibitive for use in developing countries. GemPy is a Python-based, open-source (no-cost) tool for generating three-dimensional geological models. This study uses available data and GemPy to develop the Kobo Valley Hydrogeologic Framework Model (KV-HFM), a three-dimensional HFM for Kobo Valley in northern Ethiopia, which is part of the East African Rift System. The KV-HFM is a conceptual model that comprises the hydrostratigraphy, structural features, and hydraulic properties of the Kobo Valley groundwater system. The limited data described the extent and altitude of the hydrostratigraphic units using the GemPy implicit potential–field interpolation. The KV-HFM showed the existence of an east-to-west, structural-based groundwater divide composed of volcanic rock and clay. This divide splits the catchment into two groundwater systems with limited interconnected flow. This study illustrates the use of open-source software for developing an HFM using sparse, existing geologic data.

Drone Magnetic and Time Domain Electromagnetic Exploration in Metamorphic Formations: Tool for the Identification of Strategic Sites for Aquifer Exploitation

In the contemporary era, the exploitation of aquifers in the agricultural sector has become increasingly important. In response, researchers have directed their efforts towards the formulation of effective methodologies, with geophysical prospecting emerging as a fundamental tool in locating the best underground deposits. The magnetic prospecting technique can discriminate between different categories of rocks, which facilitates the localisation of geological contacts—an essential factor in determining the strategic location of boreholes, while electromagnetic time-domain prospecting helps in the definition of sedimentary strata. In particular, this process reveals the important influence of tertiary and metamorphic formations on the regional hydrogeological framework of the studied area. The variable yields recorded in the wells in the area that have yielded good results are a clear indication of the presence of aquifers. However, it is important to note that numerous wells have been drilled in this region that have yielded negligible or even zero flow rates. Prudent selection of the location and depth of boreholes is essential to ensure proper management of this resource. The use of drones equipped with magnetometers is essential to speed up the spatial mapping process. Empirical results corroborate the accurate classification of lithological units, thus facilitating the selection of sites for groundwater abstraction. These studies serve to validate initial hypotheses and profoundly enrich our understanding of the hydrogeological dynamics of the site, thus providing avenues for optimal and sustainable exploitation and future academic research.

Water Balance Calculation for a Transboundary Aquifer System between Estonia and Latvia

Groundwater management in transboundary aquifers is complex due to the hidden nature of groundwater and its intricate interactions with surface water. The overarching issue lies in the need for improved tools and strategies to manage transboundary aquifers effectively, given their critical importance to regional water resources and ecosystems. In addressing these challenges, this study focuses on a specific transboundary aquifer system between Estonia and Latvia, situated within the sedimentary rocks of North-East Europe. Here, we propose a conceptual model approach as a valuable tool for comprehending the hydrogeological characteristics of transboundary aquifers and supporting more efficient management strategies. Utilizing open-source software MODFLOW-6, our developed conceptual model integrates the aquifer system’s geological and hydrogeological framework and simulates groundwater flow under transient conditions. Our findings highlight a profound hydrological connection between groundwater and surface water, particularly within the first 200 m of the aquifer. Despite varying cross-border groundwater flows, collaboration remains vital for sustainable transboundary aquifer management, aligning with international agreements such as the Water Convention and the EU Water Framework Directive, benefiting similar systems.

Integration of Geological, Geochemical Modelling and Hydrodynamic Condition for Understanding the Geometry and Flow Pattern of the Aquifer System, Southern Nyírség–Hajdúság, Hungary

Geological heterogeneity impacts groundwater flow patterns, necessitating a detailed hydrogeological framework for conceptualization process of aquifer systems. This research developed a new conceptual model of detailed geologic geometry by integrating 133 well-logs, 366 hydrodynamic data and 118 water samples. As new results, systematic 3D log correlation detected four distinct hydrostratigraphic units in the Southern Nyírség–Hajdúság Groundwater Body (East Hungary). The primary aquifer was identified as an incised valley 10–13 km wide and a NE–SW strike. Logan’s approach estimated the average hydraulic conductivity of the Incised Valley Unit (IVU) at 11 m/d, higher than the other three aquifers (3.2 m/d to 4.6 m/d). The average specific capacity of wells screening the IVU is 315.6 m3/d/m, in contrast with the remaining aquifers ranging from 31.6 m3/d/m to 92 m3/d/m. Pressure–depth profiles, dynamic pressure increment and hydraulic head maps revealed recharge–discharge zones and hydraulic windows between hydrostratigraphic units. The elongated pattern on the hydraulic head map at the depth of the IVU showed the existence of a preferential path along its axis within the mapped borders of the IVU. Hydrochemical analysis revealed Ca-Mg-HCO3 water type within the primary aquifer and Na-HCO3 water type in the laterally connected aquifer. The saturation index values indicated a transition from undersaturated to supersaturated state inside the main aquifer for calcite and dolomite minerals. The correlation matrix and PCA results demonstrated that the carbonate weathering process is the main factor controlling the groundwater chemistry. This integrated approach holds significance for future applications of the regional conceptual model in water management planning, sustainable aquifer development and contaminant transport modelling. It provides essential contributions to informed decision-making and the formulation of effective strategies, ensuring the long-term availability and utilization of groundwater resources.

A Spatiotemporal Characterization of Water Resource Conditions and Demands as Influenced by the Hydrogeologic Framework of the Willcox Groundwater Basin, Southeastern Arizona, USA

In the Willcox Groundwater Basin (WGB), increasing rates of agricultural groundwater withdrawal have led to significant regional groundwater level decline, threatening the basin’s long-term water resource security. Updated characterization of the basin’s water resource conditions and agricultural water demand is critically important for informing groundwater resource management efforts. We developed the hydrogeologic framework of the WGB and linked groundwater level data with land cover classification data to provide a spatiotemporal assessment of water resource conditions and agricultural development in the WGB. A correlation analysis evaluated the degree of association between the basin’s mean annual depth-to-groundwater and agricultural land cover extent. Results of this study indicate that between 2008 and 2021, agricultural land cover in the WGB increased by 29%. The average rate of groundwater level change in the basin’s measured wells was calculated at −13.8 m between 2006 and 2021. We found a strong correlation between the basin’s mean annual measured depth-to-groundwater and the annual agricultural land cover extent, further reinforcing the understanding of agricultural water use in the basin as a principal driver of groundwater level decline. The methodological framework employed proved a simple and effective way to assess groundwater resources as influenced by geology and land use. The outcomes of this study provide critical information toward improved water resources management by providing an integrated understanding of local hydrogeology, groundwater level variability, and changes in agricultural land cover in arid inland basins such as those found in Arizona, USA.

Simulation of Groundwater levels in River basin using MODFLOW

MODFLOW is a groundwater modelling program. It can be compiled according to the practical applications. Groundwater flow models are used to calculate the rate and direction of movement of groundwater through aquifers. The outputs from model simulation are the groundwater flow and hydraulic heads which are in equilibrium with the specified hydrogeological conditions (i.e. hydrologic boundaries, hydrogeological framework, initial and transient conditions, hydraulic properties and sources) demarcated for the modelled area. Groundwater models perform an important role in the management and development of groundwater resources. This paper presents the results of a mathematical groundwater model developed for the Banda districts which is part of Yamuna River basin. For this purpose, groundwater modelling software (GMS) MODFLOW was used. The model results display that the computed values are in goodfitness of the measure data, which indicate the model is reliable.

Stable isotope modeling of the groundwater discharge in complex watersheds of the state of São Paulo, Brazil

The increasing pressure on water resources demands integrated and multidisciplinary approaches to address environmental, scientific, and social issues related to water availability in watersheds. In large watersheds, many factors control hydrologic processes, requiring the application of a methodology able to describe water dynamics in complex situations. This study uses stable isotope (2H and 18O) modeling to characterize the groundwater contribution in an area of heterogeneous hydrogeological framework and advanced anthropization, located in the Piracicaba-Capivari-Jundiaí (PCJ) and Tietê-Jacaré (TJ) water management units – central-eastern portion of the state of São Paulo, Brazil. Groundwater dynamics are controlled by the direct interaction of the rainfall regime with hydrogeological domains. The modeling in this study indicates that some 60–80% of the total streamflow volume comes from groundwater. However, in the crystalline domain, the aquifer capacity for surplus water storage is smaller, increasing the importance of rainfall for generating direct runoff and maintaining the streamflow. For this characteristic, the PCJ unit is more dependent on surface water and vulnerable to prolonged droughts. On the other hand, in the sedimentary areas, the contribution of groundwater discharge is around 80%, showing that the aquifer storage capacity guarantees homogenization of the discharge throughout the year, and ensures water security in drought periods. It is concluded that the successful application of stable isotope modeling in large basins. The method, combined with others already established, can add a more dynamic and sensitive set of information for understanding hydrological processes in complex hydrographic basins, with heterogeneous hydrogeological structures and the effects of anthropic action.

Seven hydrogeological terrains characteristic of southern Ontario

Groundwater is an important component of the hydrological cycle of southern Ontario, significant to water supply and discharge to surface water. To improve overall understanding of key hydrological mechanisms, a proof-of-concept framework is presented that consists of seven typical terrains or hydrogeological settings. Geologic and topographic controls influence each setting as recorded by climate, streamflow, and groundwater level data, which trace how water moves through various parts of the southern Ontario landscape. This hydrogeological framework considers the geology of surface sediment and bedrock, permeability and porosity (inferred), physiography, and topographic gradients to define seven distinctive hydrogeological terrains. The data and watershed selection process allow for concentration on patterns of hydrologic response across the landscape, with less focus on hydrologic accounting of the terrain. The analysis identifies events in climate, stream, and well-monitoring data indicative of influxes and temporal patterns of hydrological and hydrogeological response to local geology and topography, thus differentiating styles of groundwater movement in sediment/bedrock landscapes of the Paleozoic basin and Shield margin. A case study is provided for each setting: five in sediment (clay, sand, gravel, till upland, and thick till); and two in bedrock (crystalline and carbonate). The selected hydrogeological terrains characterize ∼90% of the landscape based on a simplified geological map of southern Ontario. The main hydrogeological characteristics and behaviour of these terrains are proposed as a set of conceptual models representative of near-surface groundwater regimes. Such models can inform water and land resource management as future climates change.

A dual-continuum model (TOUGH2) for characterizing flow and discharge in a mechanically disrupted sandstone overburden

Underground hard coal mining usually disrupts the mechanical equilibrium of rock sequences, creating fractures within minor permeable rocks. The present study employs a dual-continuum model to assess how both fractured and porous sandstone media influence the percolation process in postmining setups. To test the approach, the software TOUGH2 was employed to simulate laminar fluid flow in the unsaturated zone of the Ibbenbüren Westfield mining area. Compared to other coal mining districts in Germany, this area is delineated by the topography and local geology, leading to a well-defined hydrogeological framework. Results reveal good agreement between the calculated and measured mine water discharge for the years 2008 and 2017. The constructed model was capable of reproducing the bimodal flow behavior of the adit by coupling a permeable fractured continuum with a low-conductivity rock matrix. While flow from the fractured continuum results in intense discharge events during winter months, the rock matrix determines a smooth discharge limb in summer. The study also evaluates the influence of individual and combined model parameters affecting the simulated curve. A detailed sensitivity analysis displayed the absolute and relative permeability function parameters of both continua among the most susceptible variables. However, a strong a priori knowledge of the value ranges for the matrix continuum helps to reduce the model ambiguity. This allowed for calibration of some of the fractured medium parameters for which sparse or variable data were available. However, the inclusion of the transport component and acquisition of more site-specific data is recommended to reduce their uncertainty.

Characterization of Typical Complex Hydrogeological Condition of Structurally Disturbed Shear Zone Areas of Palakkad District in Western Ghats Region, Central Kerala

Aquifer characteristics and hydrogeological condition of structurally disturbed shear zone areas of Palakkad district in Western Ghats Region of Palakkad district have been studied. The study of sub-surface hydrogeological framework of aquifers up to 200 meters below ground level (m bgl) was carried out to quantify groundwater resources, water quality, hydraulic properties with an aim to evolve a sustainable aquifer management plan in consideration with the water level monitoring and exploratory drilling. The hard rock terrains have been structurally controlled by shear a zone, which facilitate to develop numerous fractures and becomes potential repositories of groundwater. The thickness of lateritization changes between 1.02 and 24.1 m. The premonsoon depth to water level varies from 1.58 to 14.20mbgl, but that of PSM in between 0.60 and 21.50 m mbgl with a fluctuation about 6m. The quantification of groundwater potential in the aquifers of consolidated formation within the shear zones and establishment of a futuristic model of aquifer management plan for the area are the highlights of the work. Keywords: Laterite, Aquifer, Potential, Management, GIS, SWOT

Hydrogeophysical investigation using remote sensing and geoelectrical data in southeast Hiw, Qena, Egypt

The current study focuses mainly on understanding the hydrogeological framework of southeast Hiw, Qena, Egypt, as one of the promising regions for land reclamation and future development. This includes defining the groundwater aquifers and delineating their lateral extents, depths, and thicknesses, as well as getting some information on the groundwater quality in the area. These objectives have been accomplished by the integrated interpretation of remote sensing imagery and vertical electrical sounding (VES) data. The digital elevation model (DEM) of the Shuttle Radar Topography Mission (SRTM) was processed to calculate the surface gradient and delineate the structural lineaments and drainage channels. Moreover, VES data were measured at 32 locations using Schlumberger configuration. The result of the VES survey has been presented in the form of 2D cross-sections, 3D fence diagram, and hydrogeophysical contour maps. They reveal that the Quaternary sand is the main groundwater aquifer in the region. It has a depth of 30 to 58 m, a thickness of 30 to 68 m, and a resistivity range of 48 to 93 Ω.m. Furthermore, an empirical resistivity–salinity relationship was developed. Accordingly, total dissolved solids (TDS) values vary from 600 to 1500 ppm. This study improves the understanding of the subsurface geology of the region and serves as a guide for further studies regarding digging new wells and groundwater management in such an arid environment.

Development and Description of a Composite Hydrogeologic Framework for Inclusion in a Geoenvironmental Assessment of Undiscovered Uranium Resources in Pliocene- to Pleistocene-Age Geologic Units of the Texas Coastal Plain

A previously completed mineral resources assessment of the Texas Coastal Plain indicated the potential for the future discovery of uranium resources. Geoenvironmental assessments that include the hydrogeologic framework can be used as a tool to understand the potential effects of mining operations. The hydrogeologic framework for this study focused on the composite hydrogeologic unit of the tract permissive for the occurrence of uranium consisting of the upper part of the Miocene-age Fleming Formation/Lagarto Clay, Pliocene-age Goliad and Pleistocene-age Willis Sands, Pleistocene-age Lissie and Beaumont Formations, and Holocene-age alluvial sediments (fluvial alluvium and eolian sand deposits). This composite hydrogeologic unit, which contains the Chicot and Evangeline aquifers of the Gulf Coast aquifer system, is intended for inclusion in a regional-scale geoenvironmental assessment of as yet undiscovered uranium resources. This article provides (1) a brief literature review describing the geologic and hydrogeologic settings, (2) the methodology used to develop a composite hydrogeologic framework, and (3) descriptions and maps of the land-surface altitude, composite hydrogeologic unit base and midpoint depth, water-level altitude, depth of water, unsaturated and saturated zone thickness, and transmissivity and hydraulic conductivity. A composite hydrogeologic unit, created by combining geologic and hydrogeologic data and maps for individual geologic and hydrogeologic units, is intended for use as a tool in a geoenvironmental assessment to evaluate potential contaminant migration through various avenues. Potential applications include using the hydrogeologic framework as an input into a geoenvironmental assessment to help estimate the potential for (1) runoff of contaminants into surface water, (2) infiltration of contaminants into the groundwater (aquifers), or (3) movement of contaminants from the mining area through wind, groundwater-flow, or streamflow in a given permissive tract. The procedures outlined in this paper also provide a method for developing hydrogeologic frameworks that can be applied in other areas where mining may occur.

A Hydrogeologic Framework for Understanding Surface Water and Groundwater Interactions in a Watershed System in the Willamette Basin in Western Oregon, USA

A broad understanding of local geology and hydrologic processes is important for effective water resources management. The objectives of this project were to characterize the hydrogeologic framework of the Oak Creek Watershed (OCW) geographical area and examine the connections between surface water and groundwater at selected locations along the main stem of Oak Creek. The OCW area comprises the Siletz River Volcanic (SRV) Formation in the upper portion of the watershed and sedimentary rock formations in the valley. Past hydrologic and geologic studies and our field measurement data were synthesized to create a hydrogeologic framework of the watershed, including a geologic interpretation and a conceptual model of shallow, deep, and lateral groundwater flow throughout the OCW. The highly permeable geology of the SRV formation juxtaposed against the Willamette Basin’s sedimentary geology creates areas of opposing groundwater flow characteristics (e.g., hydraulic conductivity) in the watershed. The Corvallis Fault is the primary interface between these two zones and generally acts as a hydraulic barrier, deflecting groundwater flow just upstream of the fault interface. The extreme angle of the Corvallis Fault and adjacent less permeable sedimentary geology might facilitate subsurface bulk water storage in selected locations. The stream-aquifer relationships investigated showed gaining conditions are prominent in the upper watershed’s northern volcanic region and transition into neutral and losing conditions in the downstream southern sedimentary region in the valley. Agriculture irrigation seepage in the valley appeared to contribute to streamflow gaining conditions. Results from this case study contribute critical information toward enhancing understanding of local hydrogeologic features and potential for improved SW-GW resources management in areas near coastal ranges such as those found in the Pacific Northwest, USA.

Appraisal of hydrogeological framework of shallow aquifer system of Dimapur alluvial plain of Nagaland in Northeastern India

Appraisal of hydrogeological framework of shallow aquifer system of Dimapur alluvial plain of Nagaland in Northeastern India

N2 in deep subsurface fracture fluids of the Canadian Shield: Source and possible recycling processes

In addition to high concentrations of CH4 and H2, abundant dissolved N2 is found in subsurface fracture fluids in Precambrian cratons around the world. These fracture fluids have hydrogeological isolation times on order of thousands to millions and even billions of years. Assessing the sources and sinks of N2 and related (bio)geochemical processes that drive the nitrogen cycle in these long isolated systems can shed insights into the nitrogen cycle on early Earth with implications for other planets and moons. In this study, we collected dissolved gas samples from deep subsurface fracture fluids at seven sites (Kidd Creek, LaRonde, Nickel Rim, Fraser, Copper Cliff South, Thompson, and Birchtree) in the Canadian Shield. Multiple gas components (e.g., H2, O2 and Ar) were integrated with δ15NN2 values to characterize the N2 signatures. Results show that the dissolved N2 in deep subsurface fracture fluids from the Canadian Shield sites are more 15N-enriched than those from the Fennoscandian Shield and the Witwatersrand Basin in the Kaapvaal Craton. The nitrogen isotopic signatures of the Canadian Shield samples coupled with their hydrogeological framework indicate the N2 was sourced from fixed ammonium in silicate minerals in host rocks and was generated by metamorphic devolatilization. Modeling of nitrogen devolatilization from host rocks supports this interpretation, but also suggests that a second process, likely abiotic N2 reduction, is required to account for the observed 15N enrichment in the N2 samples from the Canadian Shield. A 10-year monitoring study for one of the boreholes, at 2.4 km of the Kidd Creek Observatory, shows a steady decrease in δ15NN2 values with time, which coincides with the temporal isotopic evolution of some other gas components in this borehole. Although it cannot be confirmed at this time, this isotopic shift in N2 may be potentially attributed to microbial processes (e.g., anaerobic oxidation of ammonium). Nevertheless, the large 15N enrichments for the majority of the samples in this study suggest that the nitrogen cycle in the deep saline fracture fluids in the Canadian Shield is dominated by abiotic processes. This is in contrast to the nitrogen cycles in the subsurface fracture fluids in the Fennoscandian Shield and the Witwatersrand Basin, which have been shown to be strongly affected by extant microbial ecosystems discovered in those fracture waters.

Predicting adverse scenarios for a transboundary coastal aquifer system in the Atacama Desert (Peru/Chile)

The Caplina/Concordia transboundary coastal aquifer system, located in the Atacama Desert, is the primary source of water supply for domestic use and irrigation for La Yarada-Los Palos (Peru) and Concordia (Chile) agriculture districts, and to a lesser extent, for Tacna province public supply use (Peru). Despite the scarce amount of rainfall (<20 mm/year) in the area and the limited recharge coming from the Andean highlands, this transboundary aquifer system has been overexploited mainly for agriculture since before the 2000s on the Peruvian side. Consequently, this has caused groundwater depletion and seawater intrusion. In this study, comprehensive hydrogeological information was integrated to understand the aquifer system's behavior and the effects to which it has been subjected to groundwater overexploitation. To that end, a 3D hydrogeological framework was developed using the LEAPFROG software and a constant-density groundwater flow model with equivalent heads was generated in FEFLOW software, which was adjusted with Monte Carlo analysis and conventional automated calibration. Finally, eight scenarios, considering various water resource management options proposed by the authority and potential climatic trends (CMIP6), were simulated from 2020 to 2040. The results showed that between 2002 and 2020, the increase in the seawater wedge and the average groundwater level decline were 216 hm3/year and 7 m, respectively. It is expected that the depletion will continue with a groundwater level decline between 5 and 8 m and an increase in the seawater wedge between 1120 hm3/year and 1175 hm3/year for the forecast period. The study concludes that the aquifer system will remain unsustainable for the next 20 years, regardless of the selected scenarios, and suggests that any mitigation measure requires the participation of stakeholders from Peru, Chile, and Bolivia.