AQUIFER Review of AQUIFER Technology Feasibility June 15 -- 19, 2020
Explore the source record for details and available documents.
nano-electric fuel cells↗
SEARCH · Engineering Papers
Results for “Aquifer”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins
This repository of global hydrogeologic datasets contains aquifer properties on 0.5° scale, including depth to groundwater (Fan et al., 2013), aquifer thickness (de Graaf et al., 2015), WHYMap aquifer classes (Richts et al., 2011), porosity and permeability (Gleeson et al., 2014), digitized and geo-processed from their respective sources. Globally gridded aquifer properties could be used independently to estimate global groundwater availability or used as critical inputs to the superwell model to simulate groundwater extraction and provide estimates of pumped volumes and unit costs under user-specific scenarios. Key resources related to this data are: Niazi, H., Ferencz, S., Graham, N., Yoon, J., Wild, T., Hejazi, M., Watson, D., & Vernon, C. (2024; In-prep). Long-term Hydro-economic Assessment Tool for Evaluating Global Groundwater Cost and Supply: Superwell v1. Geoscientific Model Development. superwell model repository which uses this data to simulate groundwater extraction and provides estimates of the global extractable volumes and unit-costs ($/km3) of accessible groundwater production under user-specified extraction scenarios. Repository Overview Main output: aquifer_properties.csv contains all processed outputs, including aquifer properties like porosity, permeability, aquifer thickness, and depth to groundwater. shapefiles.zip: contains all digitized GIS databases and shapefile for all aquifer properties prep_inputs.R: R script that processes the shapefiles to produce the aquifer_properties file plot_inputs.R: R script for plotting the maps and conducting preliminary analysis on the available groundwater volume basin_to_country_mapping.csv, basin_country_region_mapping.csv and continent_county_mapping.csv provide the mapping between continents, 32 energy-economic macro regions, countries, and water basins for post-processing aquifer_properties.csv Maps: Each map visualizes the spatial distribution of one of the aquifer properties across the globe map_in_Porosity.png map_in_Permeability.png map_in_Aquifer_thickness.png map_in_Depth_to_water.png map_in_Grid_area_km.png map_in_WHYClass.png Sample inputs sample_inputs.py: this script samples inputs from the aquifer_properties dataset, ensuring the sampled and original inputs maintain the same distributions sampled_data_100.csv contains 100 sampled data points and sampled_data_100.png compares their distributions Dataset Overview The main outputs are consolidated in a comprehensive aquifer_properties.csv file and include the following fields: GridCellID: Unique identifier for each (roughly 0.5°) grid cell Continent: Continent name Country: Country name GCAM_basin_ID: Identifier for GCAM hydrologic basin Basin_long_name: Full name of the basin WHYClass: Hydrogeologic classification based on WHYMap aquifer classes (Richts et al., 2011) Porosity: Soil porosity (%) (Gleeson et al., 2014) Permeability: Soil permeability (in square meters; Gleeson et al., 2014) Aquifer_thickness: Thickness of the aquifer (in meters; de Graaf et al., 2015) Depth_to_water: Depth to groundwater (in meters; Fan et al., 2013) Grid_area: Area of the grid cell (in square meters) Key References The datasets are digitized versions of global hydrogeologic properties from the following key literature sources: Depth to Groundwater: Fan, Y., Li, H., & Miguez-Macho, G. (2013). Global Patterns of Groundwater Table Depth. Science, 339(6122), 940-943. https://doi.org/10.1126/science.1229881 Aquifer Thickness: de Graaf, I. E. M., Sutanudjaja, E. H., van Beek, L. P. H., & Bierkens, M. F. P. (2015). A high-resolution global-scale groundwater model. Hydrol. Earth Syst. Sci., 19(2), 823-837. https://doi.org/10.5194/hess-19-823-2015 Porosity and Permeability: Gleeson, T., Moosdorf, N., Hartmann, J., & van Beek, L. P. H. (2014). A glimpse beneath earth's surface: GLobal HYdrogeology MaPS (GLHYMPS) of permeability and porosity. Geophysical Research Letters, 41(11), 3891-3898. https://doi.org/10.1002/2014GL059856 Aquifer classes: Richts, A., Struckmeier, W. F., & Zaepke, M. (2011). WHYMAP and the Groundwater Resources Map of the World 1:25,000,000. In J. A. A. Jones (Ed.), Sustaining Groundwater Resources: A Critical Element in the Global Water Crisis (pp. 159-173). Springer Netherlands. https://doi.org/10.1007/978-90-481-3426-7_10 Cite as Niazi, H., Watson, D., Hejazi, M., Yonkofski, C., Ferencz, S., Vernon, C., Graham, N., Wild, T., & Yoon, J. (2024). Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins. MSD-LIVE Data Repository. https://doi.org/10.57931/2307831 Contact Reach out to Hassan Niazi or open an issue in superwell repository for questions or suggestions.
aquifer thickness↗
Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins
This repository of global hydrogeologic datasets contains aquifer properties on 0.5° scale, including depth to groundwater (Fan et al., 2013), aquifer thickness (de Graaf et al., 2015), WHYMap aquifer classes (Richts et al., 2011), recharge (Döll and Fiedler, 2008; Gleeson et al., 2016), lakes (Messager et al., 2016), porosity and permeability (Gleeson et al., 2014), digitized and geo-processed from their respective sources. Globally gridded aquifer properties could be used independently to estimate global groundwater availability or used as critical inputs to the superwell model to simulate groundwater extraction and provide estimates of pumped volumes and unit costs under user-specific scenarios. Key resources related to this data are: Niazi, H., Ferencz, S. B., Graham, N. T., Yoon, J., Wild, T. B., Hejazi, M., Watson, D. J., & Vernon, C. R. (2025). Long-term hydro-economic analysis tool for evaluating global groundwater cost and supply: Superwell v1.1. Geoscientific Model Development, 18(5), 1737-1767. https://doi.org/10.5194/gmd-18-1737-2025 superwell model repository which uses this data to simulate groundwater extraction and provides estimates of the global extractable volumes and unit-costs ($/km3) of accessible groundwater production under user-specified extraction scenarios. Repository Overview Main output: aquifer_properties_rec.csv contains all processed outputs, including aquifer properties like porosity, permeability, recharge, lake areas, aquifer thickness, and depth to groundwater. shapefiles.zip: contains all digitized GIS databases and shapefile for all aquifer properties prep_inputs.R and prep_inputs_recharge_lakes.R: R scripts that process the shapefiles to produce the aquifer_properties_rec.csv file plot_inputs.R: R script for plotting the maps and conducting preliminary analysis on the available groundwater volume basin_to_country_mapping.csv, basin_country_region_mapping.csv and continent_county_mapping.csv provide the mapping between continents, 32 energy-economic macro regions, countries, and water basins for post-processing aquifer_properties_rec.csv Maps: Each map visualizes the spatial distribution of one of the aquifer properties across the globe map_in_Porosity.png map_in_Permeability.png map_in_Aquifer_thickness.png map_in_Depth_to_water.png map_in_Recharge.png map_in_Grid_area_km.png map_in_Lake_area_km.png map_in_WHYClass.png Sample inputs sample_inputs.py: this script samples inputs from the aquifer_properties_rec dataset, ensuring the sampled and original inputs maintain the same distributions sampled_data_100.csv contains 100 sampled data points and sampled_data_100.png compares their distributions Dataset Overview The main outputs are consolidated in a comprehensive aquifer_properties_rec.csv file and include the following fields: GridCellID: Unique identifier for each (roughly 0.5°) grid cell Continent: Continent name Country: Country name GCAM_basin_ID: Identifier for GCAM hydrologic basin Basin_long_name: Full name of the basin WHYClass: Hydrogeologic classification based on WHYMap aquifer classes (Richts et al., 2011) Porosity: Soil porosity (%) (Gleeson et al., 2014) Permeability: Soil permeability (in square meters; Gleeson et al., 2014) Aquifer_thickness: Thickness of the aquifer (in meters; de Graaf et al., 2015) Depth_to_water: Depth to groundwater (in meters; Fan et al., 2013) Recharge: long-term annual averaged recharge rates (in m/yr; Döll and Fiedler, 2008; Gleeson et al., 2016) Grid_area: Area of the grid cell (in square meters) Lakes_area: Area of inland lakes (in square meters; Messager et al., 2016) Key References The datasets are digitized versions of global hydrogeologic properties from the following key literature sources: Depth to Groundwater: Fan, Y., Li, H., & Miguez-Macho, G. (2013). Global Patterns of Groundwater Table Depth. Science, 339(6122), 940-943. https://doi.org/10.1126/science.1229881 Aquifer Thickness: de Graaf, I. E. M., Sutanudjaja, E. H., van Beek, L. P. H., & Bierkens, M. F. P. (2015). A high-resolution global-scale groundwater model. Hydrol. Earth Syst. Sci., 19(2), 823-837. https://doi.org/10.5194/hess-19-823-2015 Porosity and Permeability: Gleeson, T., Moosdorf, N., Hartmann, J., & van Beek, L. P. H. (2014). A glimpse beneath earth's surface: GLobal HYdrogeology MaPS (GLHYMPS) of permeability and porosity. Geophysical Research Letters, 41(11), 3891-3898. https://doi.org/10.1002/2014GL059856 Aquifer classes: Richts, A., Struckmeier, W. F., & Zaepke, M. (2011). WHYMAP and the Groundwater Resources Map of the World 1:25,000,000. In J. A. A. Jones (Ed.), Sustaining Groundwater Resources: A Critical Element in the Global Water Crisis (pp. 159-173). Springer Netherlands. https://doi.org/10.1007/978-90-481-3426-7_10 Recharge: Döll, P., & Fiedler, K. (2008). Global-scale modeling of groundwater recharge. Hydrol. Earth Syst. Sci., 12(3), 863-885. https://doi.org/10.5194/hess-12-863-2008; Gleeson, T., Befus, K. M., Jasechko, S., Luijendijk, E., & Cardenas, M. B. (2016). The global volume and distribution of modern groundwater. Nature Geoscience, 9(2), 161-167. https://doi.org/10.1038/ngeo2590 Inland Lakes: Messager, M. L., Lehner, B., Grill, G., Nedeva, I., & Schmitt, O. (2016). Estimating the volume and age of water stored in global lakes using a geo-statistical approach. Nature Communications, 7(1), 13603. https://doi.org/10.1038/ncomms13603 Cite as Niazi, H., Watson, D., Hejazi, M., Yonkofski, C., Ferencz, S., Vernon, C., Graham, N., Wild, T., & Yoon, J. (2024). Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins. MultiSector Dynamics-Living, Intuitive, Value-adding, Environment. https://doi.org/10.57931/2484226 Contact Reach out to Hassan Niazi or open an issue in the superwell repository for questions or suggestions.
aquifer thickness↗
Pantex Plant Ogallala Aquifer and Perched Groundwater Contingency Plan
The Pantex Plant Ogallala Aquifer and Perched Groundwater Contingency Plan has been developed in accordance with the requirements identified in the: • Interagency Agreement for the Pantex Superfund Site, Article 8.5 Work to be Performed, • Compliance Plan Provision of Hazardous Waste Permit No. 50284, and • Record of Decision for Groundwater, Soil, and Associated Media, Pantex Plant. A Long‐Term Monitoring System Design has been designed to monitor conditions in the perched groundwater including changes in the perched aquifer as a result of implementing the response actions. Monitoring is required for verifying the effectiveness of perched groundwater response actions (i.e., conditions in the perched aquifer are being affected as intended) and for confirming that the perched aquifer and Ogallala Aquifer characterization as defined in the Resource Conservation and Recovery Act Facility Investigation Report and the Corrective Measure Studies/Feasibility Study remains accurate. If monitoring results obtained through the monitoring network identify an unexpected condition or deviation, contingent actions will be considered and implemented as necessary to ensure continued protection of the Ogallala Aquifer and human health and the environment. Potential deviations to expected technology performance may be encountered for each of the four primary response actions that compose the selected remedy for perched groundwater; Playa 1 Pump and Treat System, Southeast Area Pump and Treat System, Southeast Area In‐Situ Bioremediation System (comprised of the Southeast In‐Situ Bioremediation System Original System, Southeast Area In‐Situ Bioremediation System Extension System, Offsite In‐Situ Bioremediation System, Perchlorate/Chromium ISB, Northeast ISB and County Road 8 ISB), and Zone 11 In‐Situ Bioremediation System. Monitoring will also be conducted to determine if there are deviations to the expected characterization, e.g., contaminants not expected as a result of the RCRA Facility Investigation characterization. Deviations to expected conditions in the Ogallala Aquifer could also be encountered if the response actions in the perched groundwater are not performing as expected, i.e., preventing contaminants from migrating to the Ogallala Aquifer. Currently, Pantex has begun investigation of detections of high explosives above groundwater protection standards in wells on the Texas Tech University property and a plume that is moving to the northeast from that area. Due to those detections, this Plan recognizes the fact that future detections in the Ogallala will be focused on first‐ time detections of analytes. After a remedy is determined, this Plan will require modification to address-deviations and contingent actions. This Plan was developed to identify the contingent actions necessary to mitigate impacts resulting from deviations to site conditions or response action performance. The Plan defines the environmental problem being addressed by the response actions, clarifies the expected conditions and objectives of the response actions, and identifies the potential deviations to the response actions (due to site conditions or technology performance) that could be encountered. The deviations were evaluated to determine the likelihood of occurrence, potential impact, and time to respond to avoid impact. The Plan also identifies the monitoring outlined in the Long‐Term Monitoring System Design Report (Consolidated Nuclear Security, 2024) and Sampling Analysis Plan (PanTeXas Deterrence, 2024) that will be used to detect the deviations. Lastly, the Plan specifies the contingent actions that could be implemented in response to the deviations. Because each response focuses on a discrete portion of the perched aquifer and contaminant plume, each response action has a different set of expected conditions, and therefore differing impacts from deviations to the site and technology expectations. As a result, the contingent actions are identified for each response action and potential deviation including specific constituents, location, and conditions. If deviations are encountered that impact the ability of the response action to meet performance objectives, the contingent actions will be focused on ensuring the response action can meet the performance objective. Contingent actions may be implemented as interim actions (ISMs/removal actions) in accordance with the Record of Decision, Interagency Agreement, and Hazardous Waste Permit‐50284, if warranted by the specific circumstances. For deviations to site characterization expected conditions, the contingent action will focus on determination of the source of the deviation, determination of the appropriate response, and evaluation of additional work to be completed. However, if the deviation to characterization impacts the performance of the response action, the contingent action will again focus on ensuring performance objectives can be met. Early source term removals and cleanup actions have been implemented to protect the Ogallala Aquifer. Because of these actions and based on modeling results, the expected conditions in the Ogallala Aquifer are that constituents of concern will not be detected above the Groundwater Protection Standards (GWPSs) nor will they reach potential points of exposure above the GWPS. The primary deviation of concern for the Ogallala is if constituents are detected in the Ogallala Aquifer near or above GWPSs. If it occurs, this change in expected conditions would require further evaluation of site and contaminant characteristics to determine an appropriate course of action. The evaluation would include additional monitoring, source identification, implementation of interim protective measures (if necessary), and delineation of extent. These evaluations are necessary to determine an appropriate response action for the Ogallala. The primary goal of the Plan is to provide for the continued protection of the Ogallala Aquifer and the health of its consumers. In recognition, this Plan presents a flexible and rational approach for making future decisions associated with confirming the change in perched and Ogallala aquifer conditions and identifying a response (technical activities, changes to response actions, regulatory oversight, and public involvement).
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
Identifying recharge sources and their impacts on a North Central New Mexico shallow aquifer using unsupervised machine learning
In this article, shallow aquifers are important but highly variable resources in arid to semi-arid regions. Limited shallow aquifer volume results in high sensitivity to recharge fluctuations, which can impact the local fauna and flora, and transport of contaminants in the aquifer or vadose zone. Aquifer response to external forcing (e.g., precipitation) is usually solved by estimating aquifer parameters and running physics-based models to match known fluctuations of hydraulic head. However, this technique is time and computationally expensive. Furthermore, high aquifer complexity decreases precision in physics-based models. Alternatively supervised machine learning is used to predict aquifer dynamics. However, these techniques rely on input data and struggle to interpret aquifer response for missing sources (i.e., snowpack data). To counter these problems, we propose an unsupervised machine learning technique (NMFk) to estimate the impact of different sources on aquifer recharge. NMFk is used to understand the influence of external forcing on shallow aquifer recharge in the Pajarito Plateau (Los Alamos, NM, USA). The results show how NMFk can be used to reduce the data dimension in a complex field dataset to three recharge signals that cause fluctuations within the field data. Here, the source signals are interpreted as rainfall, snowmelt, and a delayed aquifer response to the previous two signals. These results evidence how heterogeneous aquifers delimited by canyons incised into the Pajarito Plateau respond in similar ways to the source signals identified by NMFk. Furthermore, results show the importance of the local geology where faults act as sinks, and anthropogenic disturbances can facilitate infiltration amplifying the interpreted signal.
54 ENVIRONMENTAL SCIENCES↗
River–aquifer interactions enhancing evapotranspiration in a semiarid riparian zone: A modelling study
The hydrologic flows across the river–aquifer interface play an important role in groundwater dynamics and biogeochemical reactions within the subsurface; however, little is known about the effects of river–aquifer interactions on land surface processes. In this study, we developed a fully coupled three-dimensional (3D) land surface and subsurface model at a high resolution (~1 km) that accounts for high-frequency hydrologic exchange flow conditions to investigate how river–aquifer interactions modulate surface water budgets in the Upper Columbia-Priest Rapids watershed, a typical semiarid watershed located in the northwestern United States where river stage fluctuates in response to reservoir releases changing. Our results show that the spatiotemporal dynamics of river–aquifer interactions are highly heterogeneous, driven mainly by river-stage fluctuations. Adding 6.64 × 10 6 m 3 year –1 of water over the watershed from the river to groundwater owing to the lateral flow, river–aquifer interactions led to an increase in soil evaporation and transpiration supplied by higher soil moisture content, particularly in deeper subsurface. In a hypothetic future scenarios where a 5-m rise in river stage was assumed, the hydrologic flow exchange rates were intensified, resulting in higher surface water over the entire watershed. Overall, lateral flow induced by river–aquifer exchanges leads to an increase in evapotranspiration of ~75% in the historical period and of ~83% in the hypothetical future scenario. Finally, our study demonstrates the potential of coupled model as an effective tool for understanding river–aquifer–land surface interactions, and indicates that river–aquifer interactions fundamentally alter the water balance of the riparian zone for the semiarid watershed and will likely become more frequent and intense in the future under the effects of climate change.
54 ENVIRONMENTAL SCIENCES↗
NRAP-Open-IAM: Generic Aquifer Component Development and Testing
The Generic Aquifer Model calculates the concentrations of dissolved salt and dissolved CO 2 surrounding a leaking legacy well. The Generic Aquifer model can also estimate the size of an “impact plume” where concentration changes exceed user-specified thresholds. The model is a component of NRAP-Open-IAM, an open-source Integrated Assessment Model (IAM) developed by the National Risk Assessment Partnership (NRAP) to perform risk assessment for geologic CO 2 storage. The input parameters were selected to cover a wide range of groundwater aquifers and leakage rates. The generic aquifer model was developed using a generative adversarial deep learning network, trained using a large synthetic dataset of STOMP multiphase flow simulations. The deep learning model predictions of dissolved salt and dissolved CO 2 in the aquifer compare well to the original STOMP simulation results. The extent of aquifer impacted by leaking CO 2 or brine is calculated using a user-defined mass fraction threshold. The aquifer impact volumes calculated based on STOMP simulation results compare well to those calculated based on the deep learning model. In a provided python script, gridded observation results from the generic aquifer component of NRAP-Open-IAM are converted to HDF5 format files for monitoring design with the DREAM code.
54 ENVIRONMENTAL SCIENCES↗
Sulfate reduction accelerates groundwater arsenic contamination even in aquifers with abundant iron oxides
Groundwater contamination by geogenic arsenic is a global problem affecting nearly 200 million people. In South and Southeast Asia, a cost-effective mitigation strategy is to use oxidized low-arsenic aquifers rather than reduced high-arsenic aquifers. Aquifers with abundant oxidized iron minerals are presumably safeguarded against immediate arsenic contamination, due to strong sorption of arsenic onto iron minerals. However, preferential pumping of low-arsenic aquifers can destabilize the boundaries between these aquifers, pulling high-arsenic water into low-arsenic aquifers. We investigate this scenario in a hybrid field-column experiment in Bangladesh where naturally high-arsenic groundwater is pumped through sediment cores from a low-arsenic aquifer, and detailed aqueous and solid-phase measurements are used to constrain reactive transport modelling. Here we show that elevated groundwater arsenic concentrations are induced by sulfate reduction and the predicted formation of highly mobile, poorly sorbing thioarsenic species. Furthermore, this process suggests that contamination of currently pristine aquifers with arsenic can occur up to over 1.5 times faster than previously thought, leading to a deterioration of urgently needed water resources.
54 ENVIRONMENTAL SCIENCES↗
Impact of inland boundary conditions on seawater intrusion modeling in stratified coastal aquifers under sea-level rise
This study systematically examines the impact of inland boundary conditions on seawater intrusion modeling in general stratified coastal aquifers under sea-level rise (SLR). We derived analytical solutions for seawater wedge toe positions in stratified aquifers with general head boundary (GHB) conditions and compiled solutions for constant head boundary (CHB), constant flux boundary (CFB), and GHB conditions in both homogeneous and stratified aquifers. Additionally, we introduced the concept of equivalent boundary conditions and developed equations to ensure consistent steady-state toe positions across different boundary types. Unlike previous studies focusing on specific layering configurations, we conducted a comprehensive sensitivity analysis of seawater intrusion in general stratified aquifers using transmissivity centroid elevation (TCE). Our results show that SLR generally exacerbates seawater intrusion, with unconfined aquifers exhibiting higher sensitivity across all boundary conditions. Among boundary types, CHB conditions were the most sensitive to SLR, CFB the least, and GHB showed an intermediate response depending on its distance from the model boundary. Furthermore, seawater intrusion increases with the TCE of the stratified aquifer under SLR, suggesting that high-permeability or preferential flow layers near the surface amplify the impact of SLR. These findings emphasize the importance of selecting appropriate boundary conditions when evaluating SLR-driven seawater intrusion in stratified coastal aquifers.
Boundary condition↗
Manganese exposure from spring and well waters in the Shenandoah Valley: interplay of aquifer lithology, soil composition, and redox conditions
Manganese (Mn) is of particular concern in groundwater, as low-level chronic exposure to aqueous Mn concentrations in drinking water can result in a variety of health and neurodevelopmental effects. Much of the global population relies on drinking water sourced from karst aquifers. Thus, we seek to assess the relative risk of Mn contamination in karst by investigating the Shenandoah Valley, VA region, as it is underlain by both karst and non-karst aquifers and much of the population relies on water wells and spring water. Water and soil samples were collected throughout the Shenandoah Valley, to supplement pre-existing well water and spring data from the National Water Information System and the Virginia Household Water Quality Program, totaling 1815 wells and 119 springs. Soils were analyzed using X-ray fluorescence and Mn K-Edge X-ray absorption near-edge structure spectroscopy. Factors such as soil type, soil geochemistry, and aquifer lithology were linked with each location to determine if correlations exist with aqueous Mn concentrations. Analyzing the distribution of Mn in drinking water sources suggests that water wells and springs within karst aquifers are preferable with respect to chronic Mn exposure, with < 4.9% of wells and springs in dolostone and limestone aquifers exceeding 100 ppb Mn, while sandstone and shale aquifers have a heightened risk, with > 20% of wells exceeding 100 ppb Mn. The geochemistry of associated soils and spatial relationships to various hydrologic and geologic features indicates that water interactions with aquifer lithology and soils contribute to aqueous Mn concentrations. Relationships between aqueous Mn in spring waters and Mn in soils indicate that increasing aqueous Mn is correlated with decreasing soil Mn(IV). These results point to redox conditions exerting a dominant control on Mn in this region.
54 ENVIRONMENTAL SCIENCES↗
Linkages Between GRACE Water Storage, Hydrologic Extremes, and Climate Teleconnections in Major African Aquifers
Water resources management is a critical issue in Africa where many regions are subjected to sequential droughts and floods. The objective of our work was to assess spatiotemporal variability in water storage and related controls (climate, human intervention) in major African aquifers and consider approaches toward more sustainable development. Different approaches were used to track water storage, including GRACE/GRACE Follow On satellites for Total Water Storage (TWS); satellite altimetry for reservoir storage, MODIS satellites for vegetation indices, and limited ground-based monitoring. Results show that declining trends in TWS (60–73 km3 over the 18 yr GRACE record) were restricted to aquifers in northern Africa, controlled primarily by irrigation water use in the Nubian and NW Saharan aquifers. Rising TWS trends were found in aquifers in western Africa (23–49 km3), attributed to increased recharge from land use change and cropland expansion. Interannual variability dominated TWS variability in eastern and southern Africa, controlled primarily by climate extremes. Climate teleconnections, particularly El Nino Southern Oscillation and Indian Ocean Dipole, strongly controlled droughts and floods in eastern and southern Africa. Huge aquifer storage in northern Africa suggests that the recent decadal storage declines should not impact the regional aquifers but may affect local conditions. Increasing groundwater levels in western Africa will need to be managed because of locally rising groundwater flooding. More climate resilient water management can be accomplished in eastern and southern Africa by storing water from wet to dry climate cycles. Accessing the natural water storage provided by aquifers in Africa is the obvious way to manage the variability between droughts and floods.
Africa water storage↗
An Empirical Algorithm to Map Perennial Firn Aquifers and Ice Slabs Within the Greenland Ice Sheet Using Satellite L-band Microwave Radiometry
Perennial firn aquifers are subsurface meltwater reservoirs consisting of a meters-thick water-saturated firn layer that can form on spatial scales as large as tens of kilometers. They have been observed within the percolation facies of glaciated regions experiencing intense seasonal surface melting and high snow accumulation. Widespread perennial firn aquifers have been identified within the Greenland Ice Sheet (GrIS) via field expeditions, airborne ice-penetrating radar surveys, and satellite microwave sensors. In contrast, ice slabs are nearly continuous ice layers that can also form on spatial scales as large as tens of kilometers as a result of surface and subsurface water-saturated snow and firn layers sequentially refreezing following multiple melting seasons. They have been observed within the percolation facies of glaciated regions experiencing intense seasonal surface melting but in areas where snow accumulation is at least 25 % lower as compared to perennial firn aquifer areas. Widespread ice slabs have recently been identified within the GrIS via field expeditions and airborne ice-penetrating radar surveys, specifically in areas where perennial firn aquifers typically do not form. However, ice slabs have yet to be identified from space. Together, these two ice sheet features represent distinct, but related, sub-facies within the broader percolation facies of the GrIS that can be defined primarily by differences in snow accumulation, which influences the englacial hydrology and thermal characteristics of firn layers at depth. Here, for the first time, we use enhanced-resolution vertically polarized L-band brightness temperature () imagery (2015–2019) generated using observations collected over the GrIS by NASA's Soil Moisture Active Passive (SMAP) satellite to map perennial firn aquifer and ice slab areas together as a continuous englacial hydrological system. We use an empirical algorithm previously developed to map the extent of Greenland's perennial firn aquifers via fitting exponentially decreasing temporal L-band signatures to a set of sigmoidal curves. This algorithm is recalibrated to also map the extent of ice slab areas using airborne ice-penetrating radar surveys collected by NASA's Operation IceBridge (OIB) campaigns (2010–2017). Our SMAP-derived maps show that between 2015 and 2019, perennial firn aquifer areas extended over 64 000 km2, and ice slab areas extended over 76 000 km2. Combined together, these sub-facies are the equivalent of 24 % of the percolation facies of the GrIS. As Greenland's climate continues to warm, seasonal surface melting will increase in extent, intensity, and duration. Quantifying the possible rapid expansion of these sub-facies using satellite L-band microwave radiometry has significant implications for understanding ice-sheet-wide variability in englacial hydrology that may drive meltwater-induced hydrofracturing and accelerated ice flow as well as high-elevation meltwater runoff that can impact the mass balance and stability of the GrIS.
Julie Z Miller↗
Reactive transport modeling of the Aquifer Thermal Energy Storage (ATES) system at Stockton University, New Jersey during seasonal operations
Hydrogeochemical processes associated with Aquifer Thermal Energy Storage (ATES) operations can often impact the system performance owing to mineral precipitation either at the wellbore or in the aquifer owing to changes in temperature and fluid disequilibria. Although failure of ATES systems due to mineral precipitation ("fouling") is common, predictive reactive-transport models have rarely been applied to plan their design and operation. Here, the objective of this study is to develop a reactive-transport model by coupling thermal, hydrological, and chemical (THC) processes to evaluate effects of introduced atmospheric oxygen on water chemistry, mineral precipitation/dissolution, porosity, and permeability changes associated with an ATES system at Stockton University (New Jersey, USA). The THC model builds on a Thermal-Hydrological-Mechanical (THM) model of the site that evaluated system failure owing to possible fracturing in the caprock or around the wellbore. The causes of the system failure are not known – potential causes include hydraulic fracturing owing to elevated pump pressures that took place, a flow pathway created by one of the boreholes, or a pre-existing natural hydrologic connection between the upper unconfined aquifer and the ATES aquifer, any of which could have led to oxygenated water entering the reservoir and causing the observed Fe-oxide fouling on well screens. The THC model is used to evaluate some of the hypotheses and observations regarding system failure owing to geochemical processes. The reactive-transport code TOUGHREACT V4 was used to model the THC processes during seasonal heating and cooling operations at the Stockton ATES site over 6 years of operation. In the THC simulations, the primary effects on geochemistry were observed when the injection water is saturated with atmospheric oxygen. Simulations show greater precipitation of goethite near the cold wells as compared to the warm wells. Although volume fractions of Fe-hydroxides were relatively small, the model was aimed at processes in the aquifer at the scale of meters and larger rather than at the scale of mm or cm (i.e., a well screen). Kaolinite is the dominant precipitating phase, also around the cold wells. Illite dissolves near the cold wells and precipitates near the warm wells. There is a net decrease in the porosity near the cold wells and increase near the warm wells, although a slight amount of thermal contraction near the cold wells and expansion near the warm wells is responsible for a significant proportion of the porosity change. Owing to the coarse discretization of the numerical grid near the wells (compared to the screen thickness) the magnitude of permeability changes at the wellbore are likely underestimated. The reactive transport model in this study can be used for characterization of aquifers, optimizing the operational parameters (temperature, pressure, pH etc.), and planning of mitigation strategies for ATES systems.
15 GEOTHERMAL ENERGY↗
Impacts of Soil-aquifer Heat and Water Fluxes on Simulated Global Climate
Climate models have traditionally only represented heat and water fluxes within relatively shallow soil layers, but there is increasing interest in the possible role of heat and water exchanges with the deeper subsurface. Here, we integrate an idealized 50m deep aquifer into the land surface module of the GISS ModelE general circulation model to test the influence of aquifer-soil moisture and heat exchanges on climate variables. We evaluate the impact on the modeled climate of aquifer-soil heat and water fluxes separately, as well as in combination. The addition of the aquifer to ModelE has limited impact on annual-mean climate, with little change in global mean land temperature, precipitation, or evaporation. The seasonal amplitude of deep soil temperature is strongly damped by the soil-aquifer heat flux. This not only improves the model representation of permafrost area but propagates to the surface, resulting in an increase in the seasonal amplitude of surface air temperature of >1K in the Arctic. The soil-aquifer water and heat fluxes both slightly decrease interannual variability in soil moisture and in landsurface temperature, and decrease the soil moisture memory of the land surface on seasonal to annual timescales. The results of this experiment suggest that deepening the modeled land surface, compared to modeling only a shallower soil column with a no-flux bottom boundary condition, has limited impact on mean climate but does affect seasonality and interannual persistence.
soil↗
Impact of aquifer properties on the extent and timeline of CO 2 trapping
Abstract Geologic CO 2 sequestration in porous saline aquifers is a promising approach to reducing atmospheric concentrations of CO 2 . Reactive transport simulations provide the opportunity to analyze which factors influence geochemical reactivity in the reservoir, understand those most important for promoting CO 2 trapping, and assess individual sites. Field‐scale aquifer characterization is time and resource intensive such that here, reactive transport simulations are leveraged to enhance understanding of selected aquifer properties including porosity, permeability, depth of storage, and carbonate mineralogy on the overall CO 2 trapping potential to better select sites promoting geochemical reactivity for CO 2 trapping. There are different mechanisms for sequestrating CO 2 . Once injected, CO 2 will dissolve into the brine to create an acidic environment, resulting in the dissolution of pre‐injection formation minerals. Released ions can reprecipitate as secondary minerals. The dissolved CO 2 and mineralized CO 2 are considered as a more secure form of CO 2 trapping in this study compared to the free supercritical CO 2 . Here, a framework leveraging a controlled set of field scale simulations is developed to facilitate rapid, optimized site selection. Simulations vary aquifer properties to understand the impact of each unique property on CO 2 trapping, tracking, and comparing the amount of supercritical, aqueous, and mineralized CO 2 . The rate at which the CO 2 injected into the aquifer is converted to aqueous or mineralized CO 2 is newly defined here as the sequestration efficiency and used to compare simulation results. The reservoir depth and fraction of carbonate minerals in the formation are shown to be more important factors than reservoir porosity and permeability in affecting CO 2 trapping. However, the impact of aquifer properties on the evolution of injected CO 2 depends on the stage of the sequestration project. © 2023 Society of Chemical Industry and John Wiley & Sons, Ltd.
Energy & Fuels↗
Hydraulic Testing of the Lost Lake Aquifer Near Recovery Wells RWM 6 and RWM 7
Hydrologic tests were conducted on the Lost Lake Aquifer Zone (LLAZ) at recovery wells RWM 6 and RWM 7 in accordance with the approved test plan (Dixon, 2023). The objective of the testing was to determine well performance parameters of RWM 6 and aquifer hydraulic properties near RWM 6 and RWM 7. Tests at RWM 6 consisted of a step-drawdown test to determine well performance parameters, a constant rate pump test to determine aquifer hydraulic properties, and an aquifer recovery test also used to determine aquifer hydraulic properties. Well performance parameters determined included specific capacity, well efficiency, and head loss coefficients. Prior to testing at RWM 6, a constant rate pump test was conducted at RWM 7. Data from this testing was used to determine aquifer hydraulic properties near RWM 7.
42 ENGINEERING↗
Hydraulic Testing of the Lost Lake Aquifer Near Recovery Well RWM001R
An aquifer pumping test was conducted on the Lost Lake Aquifer Zone (LLAZ) at the recently installed recovery well RWM001R in accordance with the approved test plan (Dixon, 2024). The objective of the testing was to determine baseline well performance parameters and aquifer hydraulic conductivity. This testing consisted of two step-drawdown tests to determine well performance properties and two constant pumping rate aquifer tests with recovery periods to determine aquifer hydraulic properties. Well performance parameters determined included specific capacity, well efficiency, and head loss coefficients. The results of the step-drawdown testing are presented in Tables ES1 and ES2 whereas the results of the constant pumping rate testing are presented in Table ES3. The average specific capacity of RWM001R operating at pump capacity (~54 gpm) was 2.2 gpm/ft using data from the final step of both tests. The average well efficiency determined from both step-drawdown tests was 56% using data from the final step of both tests. For both step-drawdown tests, well efficiency was inversely related to the pumping rate with the average efficiency decreasing from approximately 79% to 56%. Aquifer head loss coefficients from both tests were comparable (Table ES2).
54 ENVIRONMENTAL SCIENCES↗
Leveraging High Resolution Remotely Sensed Evapotranspiration Data to Assess Water Availability and Vulnerability in the Memphis Aquifer Area in West Tennessee
The Memphis aquifer is an indispensable source of water for municipal, agricultural, and industrial use across west Tennessee. The aquifer's water is thousands of years old and in pristine condition due to the natural geological filtering by the aquifer’s confining clay layers. To sustain the aquifer’s water supply, recharge must happen at a rate efficient enough to keep up with pumping from the aquifer. In partnership with Protect Our Aquifer and the University of Memphis Center for Applied Earth Science and Engineering Research, this project examined multiple factors to discover thriving areas of recharge in west Tennessee. The team used Integrated Multi-satellite Retrievals for GPM (IMERG) precipitation data, ECOsystem Spaceborne Thermal Radiometer on Space Station (ECOSTRESS) evapotranspiration data, Landsat 8 and National Land Cover Database (NLCD) to track landcover change, and the North American Land Data Assimilation System (NLDAS) Noah land surface model to estimate monthly runoff. These factors were averaged seasonally from 2019 to 2022 to understand seasonal relationships within the study area. Precipitation and evapotranspiration data were used to create water balance maps and with the evaporative stress index (ESI), landcover, and runoff maps, the team pinpointed locations that are most conducive to recharge. The team determined that areas where the water balance is low, ESI is high, and where there is sufficient pervious landcover are locations where recharge happens most effectively. Identification of effective recharge locations may be used by the partners to prioritize areas in need of protection before they become impacted by industrialization and urbanization.
Lauren Webster↗