Search NASA⌕ Search

SEARCH · Search NASA

Results for “Groundwater Modeling”

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.

At least 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Numerical Study of Coupled Water and Vapor Flow, Heat Transfer, and Solute Transport in Variably-Saturated Deformable Soil During Freeze-Thaw Cycles

As climate change intensifies, soil water flow, heat transfer, and solute transport in the active, unfrozen zones within permafrost and seasonally frozen ground exhibit progressively more complex interactions that are difficult to elucidate with measurements alone. For example, frozen conditions impede water flow and solute transport in soil, while heat and mass transfer are significantly affected by high thermal inertia generated from water-ice phase change during the freeze-thaw cycle. To assist in understanding these subsurface processes, the current study presents a coupled two-dimensional model, which examines heat conduction-convection with water-ice phase change, soil water (liquid water and vapor) and groundwater flow, advective-dispersive solute transport with sorption, and soil deformation (frost heave and thaw settlement) in variably saturated soils subjected to freeze-thaw actions. This coupled multiphysics problem is numerically solved using the finite element method. The model's performance is first verified by comparison to a well-documented freezing test on unsaturated soil in a laboratory environment obtained from the literature. Then based on the proposed model, we quantify the impacts of freeze-thaw cycles on the distribution of temperature, water content, displacement history, and solute concentration in three distinct soil types, including sand, silt and clay textures. The influence of fluctuations in the air temperature, groundwater level, hydraulic conductivity, and solute transport parameters was also comparatively studied. The results show that (a) there is a significant bidirectional exchange between groundwater in the saturated zone and soil water in the vadose zone during freeze-thaw periods, and its magnitude increases with the combined influence of higher hydraulic conductivity and higher capillarity; (b) the rapid dewatering ahead of the freezing front causes local volume shrinkage within the non-frozen region when the freezing front propagates downward during the freezing stage and this volume shrinkage reduces the impact of frost heave due to ice formation. This gradually recovers when the thawed water replenishes the water loss zone during the thawing stage; and (c) the profiles of soil moisture, temperature, displacement, and solute concentration during freeze-thaw cycles are sensitive to the changes in amplitude and freeze-thaw period of the sinusoidal varying air temperature near the ground surface, hydraulic conductivity of soil texture, and the initial groundwater levels. Our modeling framework and simulation results highlight the need to account for coupled thermal-hydraulic-mechanical-chemical behaviors to better understand soil water and groundwater dynamics during freeze-thaw cycles and further help explain the observed changes in water cycles and landscape evolution in cold regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Post-fire time series photos from five sites across the Oak Creek watershed, Washington

This dataset supports a broader study examining wildfire impacts on hydrologic connectivity across 5 sites within the Oak Creek watershed and the resulting biogeochemical impacts. Sites were selected using the Advanced Terrestrial Simulator (ATS) hydrologic model to identify locations with varying groundwater contributions and hydrologic responses across different burn severity scenarios. The Retreat Fire burned from July 23 to August 2, 2024, affecting all five sites. This dataset provides time series game camera photos, while the broader study includes continuous water quality monitoring, biogeochemical sampling of water and soils, precipitation data, and organic matter analysis. The other data types and additional metadata (include site environmental information) can be found at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/3018020. Because this study is ongoing, this data package will be updated regularly to include newly collected photos. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. In addition to a readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes all column/row headers and variable definitions. This dataset is comprised of (1) file-level metadata; (2) data dictionary; (3) field metadata; (4) readme; (5) field protocol; and (6) folders of game camera photos. The game camera photos are organized by site with subfolders by month of collection. The field metadata contains a subset of the information collected that is most relevant to photo-processing. The full set of field metadata can be found at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/3018020. All files are .csv, .pdf, or .jpg.

Burn severity↗
DOE OSTI2020

Impact of Lateral Groundwater Flow and Subsurface Lower Boundary Conditions on Atmospheric Boundary Layer Development over Complex Terrain

Credible soil moisture redistribution schemes are essential to meteorological models, as lower boundary moisture influences the balance of surface turbulent fluxes and atmospheric boundary layer (ABL) development. While land surface models (LSMs) have vastly improved in their hydrologic representation, several commonly held assumptions, such as free-draining lower boundary, one-dimensional moisture flux, and lack of groundwater representation, can bias the terrestrial water balance. This study explores the impact of LSM hydrology representation on ABL development in the Weather Research and Forecasting (WRF) meteorological model. The results of summertime WRF simulations with Noah LSM, characterized by 2-m-thick soil and one-dimensional flow, are shown for a domain in the Colorado Rocky Mountain headwaters region. A reference WRF simulation is compared to 1) the same model with soil moisture initialized by the hydrologic model ParFlow; 2) a deep, free-draining simulation; and 3) WRF coupled to ParFlow, a three-dimensional, integrated groundwater-surface water model. Results show that both lateral transport of groundwater and the rate of drainage from the lower soil layer can weaken or reverse the coupling strength between evaporative fraction and ABL over a 5-month summer period. The resulting shifts in low-level moist convection in river valleys and thermally driven airflows yield strengthened anabatic upslope winds and perturbations to regional precipitation.

Forrester, Mary M.↗
DOE OSTI2022

Physical Controls on the Hydrology of Perennially Ice-Covered Lakes, Taylor Valley, Antarctica (1996–2013)

The McMurdo Dry Valleys, Antarctica, are a polar desert populated with numerous closed-watershed, perennially ice-covered lakes primarily fed by glacial melt. Lake levels have varied by as much as 8 m since 1972 and are currently rising after a decade of decreasing. Precipitation falls as snow, so lake hydrology is dominated by energy available to melt glacier ice and to sublimate lake ice. To understand the energy and hydrologic controls on lake level changes and to explain the variability between neighboring lakes, only a few kilometers apart, we model the hydrology for the three largest lakes in Taylor Valley. Here, we apply a physically based hydrological model that includes a surface energy balance model to estimate glacial melt and lake sublimation to constrain mass fluxes to and from the lakes. Results show that lake levels are very sensitive to small changes in glacier albedo, air temperature, and wind speed. We were able to balance the hydrologic budget in two watersheds using meltwater inflow and sublimation loss from the ice-covered lake alone. A third watershed, closest to the coast, required additional inflow beyond model uncertainties. We hypothesize a shallow groundwater system within the active layer, fed by dispersed snow patches, contributes 23% of the inflow to this watershed. The lakes are out of equilibrium with the current climate. If the climate of our study period (1996–2013) persists into the future, the lakes will reach equilibrium starting in 2300, with levels 2–17 m higher, depending on the lake, relative to the 2020 level.

58 GEOSCIENCES↗
DOE OSTI2024

Mitigating drought mortality by incorporating topography into variable forest thinning strategies

Abstract Drought-induced productivity reductions and tree mortality have been increasing in recent decades in forests around the globe. Developing adaptation strategies hinges on an adequate understanding of the mechanisms governing the drought vulnerability of forest stands. Prescribed reduction in stand density has been used as a management tool to reduce water stress and wildfire risk, but the processes that modulate fine-scale variations in plant water supply and water demand are largely missing in ecosystem models. We used an ecohydrological model that couples plant hydraulics with groundwater hydrology to examine how within-stand variations in tree spatial arrangements and topography might mitigate forest vulnerability to drought at individual-tree and stand scales. Our results demonstrated thinning generally ameliorated plant hydraulic stress and improved carbon and water fluxes of the remaining trees, although the effectiveness varied by climate and topography. Variable thinning that adjusted thinning intensity based on topography-mediated water availability achieved higher stand productivity and lower mortality risk, compared to evenly-spaced thinning at comparable intensities. The results from numerical experiments provided mechanistic evidence that topography mediates the effectiveness of thinning and highlighted the need for an explicit consideration of within-stand heterogeneity in trees and abiotic environments when designing forest thinning to mitigate drought impacts.

Environmental Sciences & Ecology↗
NASA NTRS1982

The application of remote sensing to resource management and environmental quality programs in Kansas

Capabilities to process data for state agencies in Kansas were upgraded through the vehicle of a low cost data processing system. Short term projects in which agencies identified areas of immediate needs, and longer terms projects, continued from previous years, are described including studies of Arkansas River irrigation; evaluation of rangeland in the Cimarron National Grassland; a model of the Walnut Creek watershed groundwater; selection of a pronghorn antelope release site; the establishment of a geographical data base for tax reassessment in Finney County; a land use/land cover inventory and hazards assessment; and applied R & D in agricultural remote sensing. The topics discussed at a short course in remote sensing and publications are listed.

Barr, B. G.↗
NASA NTRS2021

Bias Correction of Hydrologic Projections Strongly Impacts Inferred Climate Vulnerabilities in Institutionally Complex Water Systems

Water-resources planners use regional water management models (WMMs) to identify vulnerabilities to climate change. Frequently, dynamically downscaled climate inputs are used in conjunction with land-surface models (LSMs) to provide hydrologic streamflow projections, which serve as critical inputs for WMMs. Here, we show how even modest projection errors can strongly affect assessments of water availability and financial stability for irrigation districts in California. Specifically, our results highlight that LSM errors in projections of flood and drought extremes are highly interactive across timescales, path-dependent, and can be amplified when modeling infrastructure systems (e.g., misrepresenting banked groundwater). Common strategies for reducing errors in deterministic LSM hydrologic projections (e.g., bias correction) can themselves strongly distort projected climate vulnerabilities and misrepresent their inferred financial consequences. Overall, our results indicate a need to move beyond standard deterministic climate projection and error management frameworks that are dependent on single simulated climate change scenario outcomes.

Keyvan Malek↗
DOE OSTI2021

Recovering the effects of subgrid heterogeneity in simulations of radionuclide transport through fractured media

Groundwater flow and contaminant transport through fractured media can be simulated using Discrete Fracture Network (DFN) models which provide a natural description of structural heterogeneity. However, this approach is computationally expensive, with the large number of intersecting fractures necessitated by many real-world applications requiring modeling simplifications to be made for calculations to be tractable. Upscaling methods commonly used for this purpose can result in some loss of local-scale variability in the groundwater flow velocity field, resulting in underestimation of particle travel times, transport resistance and retention in transport calculations. In this paper, a transport downscaling algorithm to recover the transport effects of heterogeneity is tested on a synthetic Brittle Fault Zone model, motivated by the problem of large safety assessment calculations for geological repositories of spent nuclear fuel. We show that the variability in the local-scale velocity field which is lost by upscaling can be recovered by sampling from a library of DFN transport paths, accurately reproducing DFN transport statistic distributions and radionuclide breakthrough curves in an upscaled model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Kilometer-Scale Hydrologic Exchange Flows in a Gravel Bed River Corridor and Their Implications to Solute Migration

In this work, a well-characterized field site along a major river corridor was used to investigate the dynamic pathways and impacts of subsurface hydrogeologic structures on kilometer-scale hydrologic exchange flows between river water and groundwater. A aqueous uranium (U aq ) plume exists at the site that discharges to the Columbia River. We performed temporally intensive monitoring of specific conductance (SpC) and U aq concentrations within the plume at increasing distances from the river shoreline and at locations both within and outside a presumed subsurface pathway of lateral hydrologic exchange. Specific conductance and U aq were utilized as in-situ tracers of hydrologic exchange and associated groundwater-surface water mixing. Monitoring data was collected for an approximate two year period during which seasonal river elevations varied by more than 2 m causing distinct events of river water intrusion and retreat from the near-shore aquifer, and large temporal variations in well-water analyte concentrations. Statistical clustering analysis was applied to identify well groupings exhibiting common time series behaviors. PFLOTRAN was implemented within a modeling domain of 3.3 x 1.9 km with dynamic flow boundary conditions and accurate subsurface topography for the basal aquitard underlying the aquifer. In-silico calculations of river water tracer behavior were performed to assess predominant flow directions and velocities influenced by river stage variations and variable aquitard topography, and their collective influences on SpC and U aq concentrations. We demonstrate complicated multidirectional flow behaviors at the kilometer scale resulting from hydrologic exchange and subsurface structure that influence plume dynamics and determine the overall system behavior.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Vadose and Saturated Zone Flow and Transport Calculations for the Active Trenches of the Low-Level Burial Grounds, Hanford Site, Washington

The purpose of the fate and transport modeling described in this environmental calculation file (ECF) is to evaluate the impacts to groundwater associated with waste disposal operations at Mixed Waste Low-Level Burial Ground (LLBG) Trenches 31 and 34 to satisfy requirements in DOE O 435.1, Radioactive Waste Management. The model integrates the flow and transport in the vadose zone beneath the active trenches with the saturated zone downgradient of the trenches to predict the radionuclide concentration at the point of assessment (POA). DOE M 435.1-1, Radioactive Waste Management Manual, defines the POA as the point of highest projected dose or concentration beyond a 100 m (328 ft) buffer zone surrounding the disposed waste. The modeling is conducted in accordance with the DOE G 435.1, Implementation Guide for Use with DOE M 435.1-1, performance assessment (PA) guidelines. The modeling involves evaluation of the groundwater concentrations and radionuclide arrival times during the 1,000-year compliance and 10,000-year sensitivity-uncertainty periods per DOE O 435.1 and DOE M 435.1. This analysis does not consider radionuclide release during facility operations, only the post-closure impacts of the radionuclides to the environment. The evaluation of potential radiological dose to groundwater receptors caused by releases from a closed facility containing radioactive waste typically includes the following: (1) Release of radionuclides from that facility (2) Transport of those radionuclides through the environment, and (3) Exposure to humans to environmental concentration levels of those radionuclides The fate and transport three-dimensional (3D) model analysis involves the post-closure impacts to the environment of the technetium-99, iodone-129, and uranium (all isotopes in the waste). The residual inventory estimates include several radionuclides, but technetium-99 is typically responsible for almost all of the beta-gamma dose equivalent associated with groundwater (water resources) protection per 40 CFR 141, “National Primary Drinking Water Regulations” (e.g., see the results in WCH-520, Performance Assessment of Environmental Restoration Disposal Facility, Hanford Site, Washington; hereinafter referred to as the ERDF PA), and iodone-129 can also be a significant dose contributor for some waste (e.g., RPP-RPT-59958, Performance Assessment for the Integrated Disposal Facility, Hanford Site, Washington; hereinafter referred to as the IDF PA). Uranium does not typically factor significantly into the impacts to groundwater, even during the 10,000-year sensitivity-uncertainty period, but always remains of interest as a contaminant. This ECF does not address vadose and saturated zone modeling for Trench 94 of the 200 East Area LLBG. Current information confirms the validity of the low corrosion rate of the naval reactor plant carbon steel (HY-80), and the even lower corrosion rate of the nickel-iron-chromium alloy reactor vessel (Inconel Alloy 600) presented in DOE/EIS-0259, Final Environmental Impact Statement on the Disposal of Decommissioned, Defueled Cruiser, Ohio Class, and Los Angeles Class Naval Reactor Plants. Based on these low corrosion rates, the time to breach the reactor vessel to allow release of radionuclides from the activated metal of the reactor vessel internal structure is at least 10,000 years. This time to breach precludes the need to evaluate the vadose and saturated zone transport of contaminants released from the reactor compartment disposal packages in the 200 East Area LLBGs PA (CP-63826, Waste Release Model Package Report for the Active Trenches of the Low-Level Burial Grounds, Hanford Site, Washington).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1979

Crustal structure, geophysical models and contemporary tectonism of the Colorado Plateau

A regional analysis of the crust and upper mantle of the Colorado Plateau is presented, using existing geophysical and geological data combined with new surface wave dispersion and groundwater geothermometry data; the tectonic implications of these models are also investigated. Surface wave and seismic refraction data indicate that the crust of the interior of the Colorado Plateau is 44 + or - 3 km thick, and its crustal structure is typical of stable continental areas. Pn velocities, however, appear to be lower (7.8 km/s) than would be expected in a stable region, while silica geothermometry indicates that the average heat flow for the plateau is 55 mW per sq m (1.3 HFU).

Keller, G. R.↗
DOE OSTI2020

Whose Gas is it anyway? Differentiating the Source of a Large Soil Vapor Plume beneath Two Adjacent Waste Sites - 20487

DOE contractor CH2M Hill Plateau Remediation Company is currently responsible for conducting groundwater contamination monitoring at several RCRA treatment, storage, and disposal units located on the Hanford Site in Richland, Washington State. The Nonradioactive Dangerous Waste Landfill treatment, storage, and disposal unit presents a distinct groundwater monitoring problem because of a large multi-contaminant soil vapor plume beneath it that is a likely source of low-level volatile organic compound groundwater contamination. Adjacent to Nonradioactive Dangerous Waste Landfill is the Solid Waste Landfill. Volatile organic compounds are inventory components of both the Nonradioactive Dangerous Waste Landfill and the Solid Waste Landfill. Therefore, it is possible that both sites could be contributing to the soil vapor plume. For regulatory purposes, it is important to differentiate which site is the primary contributor of volatile organic compounds to the plume. An approach was developed to identify the primary volatile organic compound source of the soil vapor plume beneath Nonradioactive Dangerous Waste Landfill and Solid Waste Landfill. The site conceptual model hypothesis of vapor-phase volatile organic compound transport to the dissolved phase in groundwater was tested by a simple mathematical model of vapor/liquid equilibrium concentrations at the groundwater/air interface. Once it was shown that vapor-phase volatile organic compound transport to groundwater was a valid conceptual model for Nonradioactive Dangerous Waste Landfill and Solid Waste Landfill, spatial and statistical methods were used to determine the primary site contributing to the majority of volatile organic compounds to the soil vapor plume. Average groundwater chloroform, tetrachloroethene, and trichloroethene concentrations from Nonradioactive Dangerous Waste Landfill and Solid Waste Landfill monitoring network wells were plotted on maps of the facilities and immediate vicinities and compared to soil vapor sampling probe locations. Principal component analysis and mixing ratios were used to identify source contributions of each treatment, storage, and disposal unit to the plume. Results of the vapor/liquid equilibrium concentrations mathematical model showed that transport phenomena outweigh steady-state equilibria. Estimated vapor/liquid equilibrium concentrations were considerably lower than soil vapor measurements. The results indicate that dynamic vadose zone and groundwater factors such as decreased vapor concentrations with depth, vapor dilution from dispersion in the vadose zone, and advective and diffusional volatile organic compound dilution in groundwater result in groundwater volatile organic compound concentrations much less than would be measured under steady-state equilibrium conditions. Site source contribution differentiation by principal component analysis and mixing ratios was inconclusive using actual soil gas data because of the similarity in concentration values in both datasets for Nonradioactive Dangerous Waste Landfill and Solid Waste Landfill. Similar data populations suggest mixing of the vapor contributions from both sites by dispersion through the soil matrix pore spaces. However, when groundwater volatile organic compound data were compared between the Nonradioactive Dangerous Waste Landfill and Solid Waste Landfill monitoring networks, Solid Waste Landfill mean concentrations were higher, suggesting more vapor-phase volatile organic compound transport to groundwater at those locations. Simulated volatile organic compound soil vapor and groundwater datasets created to test the methods developed for this study show that the method can be successful in source differentiation when significantly different datasets are compared. This paper will describe a method of testing a conceptual model for vapor-phase contaminant transport to groundwater and for differentiating site sources of contaminants comprising a mixed-constituent soil vapor plume. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

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↗
DOE OSTI2022

Summary Analysis Hanford Site Composite (Analysis Update)

This Summary Analysis has three primary purposes: (1) Define the overall technical scope and approach for conducting the Hanford Site CA update; (2) Define the approach for complying with modeling requirements in the Williams (2012) memorandum; and (3) Serve as the mechanism for the Groundwater Vadose Zone Executive Council to review and approve the overall modeling approach for the Hanford Site CA.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2020

Enhanced Groundwater Flow on and Below Vera Rubin Ridge, the Murray Formation, Gale Crater: Evidence from Thermochemical Modeling.

NASA’s Mars Science Laboratory Curiosity rover has been exploring Vera Rubin ridge (VRR), part of the Murray formation in Gale crater, Mars, between sol 1809 and 2302. Evidence for Fe-oxides and phyllosilicates in mineralogical and geochemical data for this region was returned by Curiosity [1-5]. We applied thermochemical modeling to con-strain the formation conditions of the phyllosilicate-hematite assemblage identified on and below VRR. Average alteration compositions for the Murray formation on and below VRR were derived using CheMin and APXS data. These compositions were reacted with Gale Portage Water (GPW) between 25–100 °C and for 10% and 50% Fe3+/Fetot of the host rock [6]. Here we summarize models run at 50 °C and 10% Fe3+/Fetot for alteration compositions derived from Murray host rock compositions.

Turner, S. M. R.↗
DOE OSTI2020

Vadose Zone Flow and Transport Parameters Data Package for the Hanford Site Composite Analysis and Cumulative Impact Evaluation

This report provides a description of the basis for the development and implementation of a conceptual model for vadose zone flow and transport for the composite analysis (CA) groundwater pathway analysis and the cumulative impact evaluation (CIE). The parameterization for a numerical model is intimately linked to the conceptual model framework. The report describes the basis for the selection of hydraulic and transport parameters for the hydrostratigraphic units (HSUs) identified in the 200 East and 200 West Areas. Whenever data are sparse or unavailable, surrogate hydraulic properties are chosen based on samples collected within the 200 Areas and nearby locations that are representative of sediments characteristic of the HSUs identified elsewhere.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Globally Gridded Groundwater Extraction Volumes and Costs under Six Depletion and Ponded Depth Targets

This repository contains simulated outputs from superwell – a hydro-economic tool for long-term assessment of groundwater cost and supply – providing globally gridded groundwater extractable volumes and associated unit costs ($/km³) for accessible groundwater production, based on a variety of user-defined depletion and ponded depth scenarios. Key model documentation: 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 Find the source code of the superwell model on GitHub: https://github.com/JGCRI/superwell Repository Overview Main output: superwell_outputs.7z contains 6 files (4.5 GB) named as superwell_py_deep_all_0.*PD_0.*DL.csv. These files present superwell outputs of global groundwater extraction volumes and cost estimates on a 0.5° scale for six scenarios with different Ponded Depth (PD; 0.3 and 0.6 m) and Depletion Limit (DL; 5%, 25%, and 40% of available volume) targets over the entire pumping lifetime of a grid cell superwell_py_deep_all_0.3PD_0.25DL_sample_100.csv contains superwell outputs for 100 data points sampled to match the global inputs' distribution superwell_py_deep_all_0.3PD_0.25DL_Grid_72548.csv contains superwell output for a single grid cell concept_v5.png provides an overview of the superwell workflow Outputs Description year_number: year of pumping depletion_limit: set depletion limit (DL) as a volume fraction of total available groundwater Mappings: continent, country, gcam_basin_id, Basin_long_name, grid_id: geographic identifiers and basin information Inputs: grid_area (km²): area of the grid cell whyclass: hydrogeological classification of the aquifer permeability (m/day), porosity (%), total_thickness (m), depth_to_water (m): aquifer properties. The geo-processed input data has been published separately: https://doi.org/10.57931/2307831 Model outputs: orig_aqfr_sat_thickness (m), aqfr_sat_thickness (m): original and remaining/instantaneous saturated thickness of the aquifer hydraulic_conductivity (m/day), transmissivity (m²/day): hydraulic properties of the aquifer radius_of_influence (m), areal_extent (km²): well radius and area of influence from the center of the well number_of_wells (-): number of wells in a grid cell determined by a ratio of well area and grid area max_drawdown (m), drawdown (m), drawdown_interference (m): well and aquifer drawdown during extraction total_head (m): total lift for the groundwater (depth to water plus drawdown) total_well_length (m): total depth of wells drilled well_yield (m³/day): pumping rate or well yield power (kW), energy (kWh): power and energy required for pumping groundwater Volume Outputs: volume_produced_perwell (m³), cumulative_vol_produced_perwell (m³): production volume metrics per well volume_produced_allwells (m³), cumulative_vol_produced_allwells (m³): aggregate extraction volumes for all wells in a grid cell available_volume (m³): available groundwater in storage for the grid cell as determined by aquifer properties depleted_vol_fraction: fraction of total volume pumped over available volumes in a grid cell (same as depletion limit) Cost Outputs: well_installation_cost ($): well installation cost based on the hydrogeological complexity of the aquifer annual_capital_cost, maintenance_cost, nonenergy_cost ($): nonenergy costs energy_cost_rate ($/kWh): electricity rate energy_cost ($): energy cost of pumping groundwater total_cost_perwell ($), total_cost_allwells ($): total annual energy and non-energy cost for each and all wells in a grid cell a unit_cost ($/m³), unit_cost_per_km3 ($/km³), unit_cost_per_acreft ($/acre-ft): total cost of pumping a unit of groundwater, indicated for different spatial units Key Resources Model documentation: 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. Input data: 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 superwell source code: https://github.com/JGCRI/superwell Cite as Niazi, H., Ferencz, S., Yoon, J., Graham, N., Wild, T., Hejazi, M., Watson, D., & Vernon, C. (2024). Globally Gridded Groundwater Extraction Volumes and Costs under Six Depletion and Ponded Depth Targets. MSD-LIVE Data repository. https://doi.org/10.57931/2307832 Contact Reach out to Hassan Niazi or Stephen Ferencz or open an issue in the superwell repository for questions or suggestions.

Earth Systems↗