Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “M CODES”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

First-collision source treatment for ray effect mitigation in discrete-ordinate radiation transport solutions

Deterministic transport codes play a fundamental role in the modelling and simulation of neutron transport. One of the most common deterministic methods is the method of discrete ordinates, also known as the S method. While offering significant advantages over other deterministic methods or stochastic methods like Monte Carlo, the method of discrete ordinates suffers from non-physical artifacts in its local solution due to its discretization of angle. These artifacts, referred to as ray effects because of their ray-like appearance, tend to be worse in problems with small sources in areas with little scattering. Significant effort has gone into developing methods to mitigate ray effects, such as the first-collision source treatment, which separates the angular flux into the uncollided and collided fluxes and solving them using non-traditional techniques such as ray tracing. One such code capable of doing this is Lawrence Livermore National Laboratory's deterministic transport code ARDRA. Current ray tracing methods typically trace to a set of points inside a zone to compute an overall flux. However, this approach has significant drawbacks, such as a low order of convergence and not being conservative. Therefore, a new method has been developed that traces instead to a set of points on each of a zone's surfaces and computing the currents, before using these to obtain the flux. A comparison between these two ray tracing methods showed significant advantages to the new surface method, including inherent conservation, a higher convergence rate, and an increase in calculable information like leakage. This work performed under the auspices of the U.S. Department of Energy by Lawrence Liver- more National Laboratory under Contract DE-AC52-07NA27344. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Radiation burnthrough measurements to infer opacity at conditions close to the solar radiative zone–convective zone boundary

Recent measurements at the Sandia National Laboratory of the x-ray transmission of iron plasma have inferred opacities much higher than predicted by theory, which casts doubt on modeling of iron x-ray radiative opacity at conditions close to the solar convective zone-radiative zone boundary. An increased radiative opacity of the solar mixture, in particular iron, is a possible explanation for the disagreement in the position of the solar convection zone-radiative zone boundary as measured by helioseismology and predicted by modeling using the most recent photosphere analysis of the elemental composition. In this report we present data from radiation burnthrough experiments, which do not support a large increase in the opacity of iron at conditions close to the base of the solar convection zone and provide a constraint on the possible values of both the mean opacity and the opacity in the x-ray range of the Sandia experiments. The data agree with opacity values from current state-of-the-art opacity modeling using the CASSANDRA opacity code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Ballooning theory for micro-tearing mode in tokamak

This paper aims to investigate the impact of magnetic drift on the linear micro-tearing mode by using a kinetic approach to derive a reduced two-field eigen system in real space. Here, since the magnetic drift in real space has derivatives, it is more efficient to solve the mode equations in a Fourier-ballooning representation using the two-dimensional (2D) ballooning transform. The lowest-order eigen system in the Fourier-ballooning representation consists of two integral equations, which are numerically solved using the finite difference method for both eigenvalues and wave functions. The main results will be presented through graphical eigenvalue scans for each parameter. Furthermore, we present a graphical comparison between the predictions of the ballooning theory and GENE gyrokinetic code simulation in the pedestal region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

2D radial-azimuthal particle-in-cell benchmark for E × B discharges

In this paper we propose a representative simulation test-case of E × B discharges accounting for plasma wall interactions with the presence of both the electron cyclotron drift instability and the modified-two-stream-instability. Seven independently developed particle-in-cell (PIC) codes have simulated this benchmark case, with the same specified conditions. The characteristics of the different codes and computing times are given. Here, results show that both instabilities were captured in a similar fashion and good agreement between the different PIC codes is reported as main plasma parameters were closely related within a 5% interval. The number of macroparticles per cell was also varied and statistical convergence was reached. Detailed outputs are given in the supplementary data, to be used by other similar groups in the perspective of code verification.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Progress in DIII-D towards validating divertor power exhaust predictions

In this study, UEDGE simulations highlight the role of cross-field drifts on the onset of detached conditions, and new calibrated divertor vacuum ultra violet (VUV) spectroscopy is used to challenge the predictions of radiative constituents in these simulations. UEDGE simulations for DIII-D H-mode plasmas with the open divertor with the ion ∇B-drift towards the X-point show a bifurcated onset of the low field side (LFS) divertor detachment, consistent with experimentally observed step-like detachment onset (Jaervinen A.E. et al 2018 Phys. Rev. Lett. 121 075001). The divertor plasma in the simulations exhibits hysteresis in upstream separatrix density between attached and detached solution branches. Reducing the drift magnitude by a factor of 3 eliminates the step-like detachment onset in the simulations, confirming the strong role of drifts in the bifurcated detachment onset. When measured local plasma densities and temperatures are within proximity of predicted values in the simulations, there is no shortfall of the local emission of the dominant resonant radiating lines. However, the simulations systematically predict a factor of two lower total integrated radiated power than measured by the bolometer with the difference lost through radial heat flow out of the computational domain. Even though there is no shortfall in the emission of the dominant lines, a shortfall of total radiated power can be caused by underpredicted spatial extent of the radiation front, indicating a potential upstream or divertor transport physics origin for the radiation shortfall, or shortfall of radiated power in the spectrum between the dominant lines. In addition to the underpredicted spatial extent, in detached conditions, the simulations overpredict the peak radiation and dominant carbon lines near the X-point, which can be alleviated by manually increasing divertor diffusivity in the simulations, highlighting the ad hoc cross-field transport as one of the key limitations of the predictive capability of these divertor fluid codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Prediction of tissue optical properties using the Monte Carlo modeling of photon transport in turbid media and integrating spheres

Monte Carlo methods are an established technique for simulating light transport in biological tissue. Integrating spheres make experimental measurements of the reflectance and transmittance of a sample straightforward and inexpensive. This work presents an extension to existing Monte Carlo photon transport methods to simulate integrating sphere experiments. Crosstalk between spheres in dual-sphere experiments is accounted for in the method. Analytical models, previous works on Monte Carlo photon transport, and experimental measurements of a synthetic tissue phantom validate this method. We present two approaches for using this method to back-calculate the optical properties of samples. Experimental and simulation uncertainties are propagated through both methods. Both back-calculation methods find the optical properties of a sample accurately and precisely. Our model is implemented in standard Python 3 and CUDA C++ [J. Nickolls, I. Buck, M. Garland, and K. Skadron, ACM Queue 6 , 40 ( 2008 ) ] and is publicly available in Code 1.

Cook, Patrick D. (ORCID:0000000279345428)↗
DOE OSTI2026

Meteorological and Soil Data from Ecohydrology Sensor Towers at Pump House and Snodgrass Mountain in East River Watershed, Colorado, 2019-2025

This data package includes hourly meteorological and soil sensor data at eight ecohydrology monitoring sites in East River Watershed, Colorado as part of the Watershed Function Scientific Focus Area (WFSFA) research led by Lawrence Berkeley National Lab (LBNL). Four field sites were located on the hillslope of East River (ER) near Pump House (PH) at Mount Crested Butte (ER-PHS1 to 4), and the other four are in the Snodgrass Mountain (SG) area (SG-EHS5 to 8). In terms of vegetation cover, three sites are in montane grasslands (ER-PHS1, ER-PHS2, and SG-EHS5), three are below evergreen conifer canopy (ER-PHS3, SG-EHS6, and SG-EHS7), and two are below deciduous aspen canopy (ER-PHS4 and SG-EHS8). The monitoring period began in October 2019 at the East River sites, in October 2020 at SG-EHS5 and SG-EHS6, and in October 2021 at SG-EHS7 and SG-EHS8. In September 2024, all four East River sites were fully retired. The four Snodgrass Mountain sites remain active. Each site is equipped with a comprehensive suite of meteorological sensors on a tripod and soil sensors that measure weather, energy fluxes, and soil variables. This data package includes measurements from ten different types of sensors and up to thirteen individual sensors per site, including (1) a weather station (measurement height ranges from 2.8~3.8 meters (m) above ground), (2) a quantum sensor for photosynthetic active radiation (PAR) (2.4~3.3m), (3) a net radiometer (1.7~2.1m), (4) an infrared radiometer (1.6~2.2m), (5) a sonic distance sensor (1.5~1.9m), (6) a soil carbon dioxide (CO2) flux chamber (0m), (7) a soil heat flux plate (-0.05m below ground), (8) a soil oxygen sensor (-0.3m), (9) a soil water potential sensor (-0.3m), and (10) soil water content sensors at 3~4 depths (-1.15 ~ -0.1m). A total of twenty-three variables is reported in this data package, including (1) atmospheric variables: air temperature (TA), atmospheric pressure (PA), vapor pressure (VP), and vapor pressure deficit (VPD), (2) precipitation variables: rain precipitation (P) and snow depth (D_SNOW), (3) energy fluxes variables: four-component net radiation (NETRAD) (shortwave/longwave incoming/outgoing radiation, SW_IN, SW_OUT, LW_IN, LW_OUT), photosynthetic photon flux density (PPFD), and soil heat flux (G), (4) soil variables: soil water content (SWC), soil water potential (SWP), soil temperature (TS), soil bulk electrical conductivity (COND_SOIL), and soil gaseous oxygen concentration (O2_SOIL), (5) wind variables: two-dimensional wind speed (WS), gust speed (WS_MAX), and wind direction (WD), and (6) surface variables: surface infrared temperature (T_CANOPY) and soil CO2 flux (CO2_SOIL). Please see the Methods section for data processing and QA/QC steps taken to generate the hourly datasets. The following files are included in this data package (notes on version: v{x}-{y}, where x is the metadata version, and y is the data version, when applicable): (1) “metadata_site_v{x}-{y}.csv” - a site metadata file that summarizes location information of all sites, including site ID, description, coordinates, timeframe, elevation, and vegetation cover, (2) “metadata_instrument_v{x}-{y}.csv” - an instrument metadata file that summarizes sensor information of all sites, including sensor manufacturer and model, measurement height, and sampling and averaging interval of all variables, (3) "data_{SITE_ID}_v{x}-{y}.csv" - eight data files that contain hourly data of each site indicated by {SITE_ID} in the filename, (4) “/figure/data_{SITE_ID}_v{x}-{y}.png" - eight figures that help visualize data of each site indicated by {SITE_ID} in the filename, (5) “/photo/*” - photos of each site indicated by {SITE_ID} in the filename, and (6) four file level metadata (flmd.csv) and data dictionary (*_dd.csv) files that summarize file, header, column, and variable information of all files. Notes: (1) Measurement height: Each variable name is followed by conventional positional qualifiers “H_V_R”, where H indicates the relative horizontal positions of that specific variable, V the vertical positions, and R the replicates. In this data package, only the vertical qualifier V varies, and V increases from the highest vertical position (V=1) to the lowest. Variables with the same qualifier are not necessarily measured by the same sensor, and the same variable with the same qualifier across different sites are not necessarily measured at the same height. Please refer to “metadata_instrument.csv” for the sensor information and measurement heights, and whether a variable is measured below the canopy. (2) Variable availability: Snow depth is not available at ER-PHS3 and SG-EHS7. SWC, soil temperature, and soil bulk EC at the deepest depth (<-1m) are not available at SG-EHS6 and SG-EHS7. The missing value code for numeric variables is -9999, except for SWP. For SWP, the missing value code is +9999, because SWP values are negative. (3) Sampling frequency: Please refer to “metadata_instrument.csv” for the increase of sampling frequency of some variables from 30-min to 1-min at ER-PHS1 to 4 in July 2020. (4) Sensors: While the methods of each sensor are not detailed, all sensors are commercially available, and their methods can be found in their manuals. Please refer to “metadata_instrument.csv” for the sensor manufacturer and model information. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

pyJSPEC - A Python Module for IBS and Electron Cooling Simulation

The intrabeam scattering is an important collective effect that can deteriorate the property of a high-intensity beam and electron cooling is a method to mitigate the IBS effect. JSPEC (JLab Simulation Package on Electron Cooling) is an open-source C++ program developed at Jefferson Lab, which simulates the evolution of the ion beam under the IBS and/or the electron cooling effect. The Python wrapper of the C++ code, pyJSPEC, for Python 3.x environment has been recently developed and released. It allows the users to run JSPEC simulations in a Python environment. It also makes it possible for JSPEC to collaborate with other accelerator and beam modeling programs as well as plentiful python tools in data visualization, optimization, machine learning, etc. In this paper, we will introduce the features of pyJSPEC and demonstrate how to use it with sample codes and numerical results.

Zhang, H.↗
DOE OSTI2020

CTF Theory Manual (V.4.2)

Coolant-Boiling in Rod Arrays—Two Fluids (COBRA-TF) is a thermal/hydraulic (T/H) simulation code designed for light water reactor (LWR) vessel analysis. It uses a two-fluid, three-field (i.e. fluid film, fluid drops, and vapor) modeling approach. Both sub-channel and 3D Cartesian forms of 9 conservation equations are available for LWR modeling. The code was originally developed by Pacific Northwest Laboratory in 1980 and had been used and modified by several institutions over the last few decades. COBRA-TF also found use at the Pennsylvania State University (PSU) by the Reactor Dynamics and Fuel Modeling Group (RDFMG) and has been improved, updated, and subsequently re-branded as CTF. As part of the improvement process, it was necessary to generate sufficient documentation for the open-source code which had lacked such material upon being adopted by RDFMG. This document serves mainly as a theory manual for CTF, detailing the many two-phase heat transfer, drag, and important accident scenario models contained in the code as well as the numerical solution process utilized. Coding of the models is also discussed, all with consideration for updates that have been made when transitioning from COBRA-TF to CTF. Further documentation outside of this manual is also available at RDFMG which focus on code input deck generation and source code global variable and module listings

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Representing the Future Role of Hydropower and Pumped Storage Hydropower (PSH) in Electricity Planning Tools

Existing tools for long-term electric sector planning struggle to represent hydropower's nuanced site-specific technical and operating characteristics, which depend on technical specifications as well as water management practices and regulations. As a result, long-term planning models and tools insufficiently characterize hydropower value and incentives, and they cannot fully represent the role hydropower can play in a future electricity system that could include a high penetration of variable wind and solar generation, battery storage, and other low-carbon technologies. This presentation demonstrates the culmination of a multi-year effort to enhance hydropower representations in electricity planning models at the National Renewable Energy Laboratory (NREL), as part of the U.S. Department of Energy (USDOE) HydroWIRES Initiative. New modeling techniques are demonstrated using the NREL Regional Energy Deployment System (ReEDS), an open-access electric sector capacity expansion model used extensively in a wide range of technology deployment and integration analysis, including the 2016 USDOE Hydropower Vision. ReEDS uses a least-cost optimization approach to understand investment and operation of electricity generation, storage, and transmission technologies under future scenarios of electricity technology innovation, demand, policy, and other sectoral drivers. ReEDS was modified to better represent value and opportunities for both pumped storage hydropower (PSH) and hydropower systems without pumping. We incorporated a new national closed-loop PSH resource and cost assessment to explore new PSH deployment opportunities and added plant-level data to better represent the existing PSH fleet. New upgrade pathways enable opportunities for enhanced hydropower flexibility by adding pumps, upgrading dispatchability, increasing capacity, or increasing energy availability. The model was also modified to better represent the value of long-duration energy storage beyond diurnal time scales, allowing both hydropower and PSH to better balance energy supply and demand variations in high-renewable systems. These new features are demonstrated under reference and high-renewable futures and a range of sensitivity scenarios to understand which hydropower and PSH deployment and upgrade opportunities are the most attractive. These scenarios indicate potential for new closed-loop PSH deployment and for hydropower flexibility improvements to have important impacts on long-term electricity system emissions and economic outcomes. Increasing flexibility of the existing hydropower fleet can reduce the need to invest in new flexible grid technologies and help achieve decarbonization goals. Systems with sufficient energy storage could also be valuable for balancing seasonal differences in renewable energy availability, particularly from solar energy. The methods developed for ReEDS and subsequent scenario results reveal important considerations for future hydropower and grid system planning, and all data and code is freely available in a public code repository for use throughout the hydropower industry.

capacity expansion↗
DOE OSTI2026

Validation of SPH code Spheral to model interacting solid bodies in a supersonic flow

Contemporary discussions of planetary defense involve analyzing the risks posed by smaller sized, 20 to 200 m diameter, asteroids which are capable of breaking up in the atmosphere and generating a blast wave. Consequence assessments for this size class of asteroids are performed through fast-running analytic or semi-analytic models which are informed by high-fidelity hydrocode simulations of asteroid entry and breakup. However, insufficient historical data necessitates validating the independent physical processes which dominate airburst events. Here, the Fluid Solid Interface Smoothed Particle Hydrodynamics solver was previously used by Pearl et al. in 2023 to model the Chelyabinsk airburst and is used here to perform a series of validation simulations. The first effort involves modeling a cylinder in a hypersonic flow and comparing the bow shock geometry to that predicted by analytic theory. The second effort involves modeling the separation of two spherical bodies in supersonic flow and validating against experimental footage. Combined, these exercises demonstrate the ability of the code to model the flight-path of interacting solid bodies in a hypersonic flow.

Airburst↗
DOE OSTI2022

Simulation of Stark-broadened Hydrogen Balmer-line Shapes for DA White Dwarf Synthetic Spectra

Abstract White dwarfs (WDs) are useful across a wide range of astrophysical contexts. The appropriate interpretation of their spectra relies on the accuracy of WD atmosphere models. One essential ingredient of atmosphere models is the theory used for the broadening of spectral lines. To date, the models have relied on Vidal et al., known as the unified theory of line broadening (VCS). There have since been advancements in the theory; however, the calculations used in model atmosphere codes have only received minor updates. Meanwhile, advances in instrumentation and data have uncovered indications of inaccuracies: spectroscopic temperatures are roughly 10% higher and spectroscopic masses are roughly 0.1 M ⊙ higher than their photometric counterparts. The evidence suggests that VCS-based treatments of line profiles may be at least partly responsible. Gomez et al. developed a simulation-based line-profile code Xenomorph using an improved theoretical treatment that can be used to inform questions around the discrepancy. However, the code required revisions to sufficiently decrease noise for use in model spectra and to make it computationally tractable and physically realistic. In particular, we investigate three additional physical effects that are not captured in the VCS calculations: ion dynamics, higher-order multipole expansion, and an expanded basis set. We also implement a simulation-based approach to occupation probability. The present study limits the scope to the first three hydrogen Balmer transitions (H α , H β , and H γ ). We find that screening effects and occupation probability have the largest effects on the line shapes and will likely have important consequences in stellar synthetic spectra.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Accelerating x-ray tracing for exascale systems using Kokkos

The upcoming exascale computing systems Frontier and Aurora will draw much of their computing power from GPU accelerators. The hardware for these systems will be provided by AMD and Intel, respectively, each supporting their own GPU programming model. The challenge for applications that harness one of these exascale systems will be to avoid lock-in and to preserve performance portability. We report here on our results of using Kokkos to accelerate a real-world application on NERSC's Perlmutter Phase 1 (using NVIDIA A100 accelerators) and Crusher, the testbed system for OLCF's Frontier (using AMD MI250X). By porting to Kokkos, we successfully ran the same X-ray tracing code on both systems and achieved speed-ups between 13 % and 66 % compared to the original CUDA code. Finally, these results are a highly encouraging demonstration of using Kokkos to accelerate production science code.

97 MATHEMATICS AND COMPUTING↗
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↗