Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data depth”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Retrospective Analysis of Backwater Habitat Availability Using Remote Sensing

Low-velocity channel-margin habitats serve as important nursery habitats for the endangered Colorado pikeminnow (Ptychocheilus lucius) in the middle Green River between Jensen and Ouray, Utah. These habitats, known as backwaters, are associated with emergent sandbars, and are shaped and reformed annually by peak flows. Our recent knowledge about backwater characteristics and dynamics that is summarized in the synthesis report (Grippo et al. 2015) was based on detailed annual survey data from a relatively small sample of backwaters that were collected from 2003 to 2014 and reach-wide evaluations of backwater surface area based on manual interpretation of aerial and satellite imagery. Methods that bridge the gap between the detailed surveys from a small number of backwaters and the reach-wide assessment of their surface area would enable an assessment of the availability of backwater habitats that meet the minimum depth requirements for suitable habitat for Colorado pikeminnow. In 2015 Argonne National Laboratory (Argonne) tested three regression models—linear, multiple, and partial least square (PLS) regression models—for estimating backwater depth using National Agriculture Imagery Program (NAIP) imagery collected in July 2006 that covered the Jensen-Ouray reach of the Green River (Hamada and LaGory 2016). The results suggested that a PLS regression model showed high correlation with the reference depth (R 2 = 0.69) and had the most unbiased and consistent estimate of backwater depth. The results also indicated that the PLS model would provide reasonable estimates of the amount of habitat providing a minimum suitable depth of 30 cm for young-of-the-year Colorado pikeminnow, even though absolute depth estimates may be uncertain for backwater areas deeper than approximately 40 cm. The study also provided insights regarding the amount and the selection of calibration and validation (cal-val) data needed for improving the accuracy of depth prediction.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Sub-daily virus sampling at the Bermuda Atlantic Time Series reveals diel and depth-structured population dynamics without community-level shifts

Ocean microbes contribute to biogeochemical cycles and ecosystem function, but they do so under top-down pressure imposed by viruses. While viruses are increasingly understood spatially and beginning to be incorporated into predictive modeling, high-frequency ocean virus dynamics remain understudied due to methodological challenges. Here we sampled stratified Bermuda Atlantic Time Series (BATS) waters for 112 hours at sub-daily 4- (surface) or 12- (deep chlorophyll maximum) hour intervals, purified viral particles from these samples, sequenced their metagenomes, and used the resulting data to characterize high-frequency virus community dynamics. Aggregated community diversity metrics changed with depth, but were not statistically significant temporally at a fixed location. However, finer-scale population-level analyses revealed both depth and temporal change, including physicochemical depth-driven differences and, in surface waters, thousands of viral populations that exhibited statistically significant diel rhythms. Statistical analyses revealed three main archetypes of temporal dynamics that themselves differed in abundance patterns, host predictions, viral taxonomy, and gene functions. Among these, highlights include viruses resembling an archetype with a night peaking pattern in activity that include an over-representation of viruses that putatively infect Prochlorococcus, a phototrophic cyanobacteria. Together, these efforts provide baseline community- and population-scale short-time-frame observations relevant to future climate state modeling.

Carrillo, Alfonso [The Ohio State University, Colu↗
DOE OSTI2021

Geothermal Play-Fairway Analysis of Washington State Prospects: Final Report

The Washington State Geothermal Play-Fairway Analysis overcomes the exploration challenges posed by dense vegetation, glacial deposits, and extreme precipitation. The geothermal play-fairways we target are locations where heat, permeability, and saturated porosity are present in sufficient volume to provide adequate heat exchange at depths accessible by modern drilling technology. The three study areas lie along the Cascade Range magmatic arc and are near Mount Baker, Mount St. Helens, and the Wind River Valley. The seven-year project is divided into three phases. In Phase 1 we build on a previous statewide assessment of geothermal resources and develop an initial modeling approach. The results are a series of favorability, uncertainty, and risk maps for three targeted study areas. Based on these initial results, we collect new geologic and geophysical data to further refine our modeling and reduce exploration uncertainty in Phase 2. We improve the modeling method to handle the new data and update the favorability, uncertainty, and risk maps. We also update the conceptual geothermal resource models. In Phase 3 we validate our modeling approach by drilling two temperature-gradient holes and collecting and analyzing core, image logs, and new geochemistry. Our modeling approach improves on an earlier statewide method through a more-rigorous and detailed assessment of heat and permeability. Permeability potential is assessed through geomechanical modeling of the deformation that can generate and maintain reservoir porosity and permeability. Metrics to inform heat potential include temperature-gradient wells, which are sparse in Washington; proximity of Quaternary volcanic vents and young intrusive rock; spring temperature; and reservoir temperature inferred from geothermometry. We weight the individual components using an expert-guided approach known as the Analytical Hierarchy Process. During Phase 2 we also develop a fluid-filled fracture model, and an infrastructure model that helps to delineate areas which are more favorable for geothermal development based on proximity to transmission lines, elevation, land ownership and use restrictions, and availability of process water. New geologic and geophysical data is collected during Phase 2 in each of our three main study areas. At Mount Baker and north of Mount St. Helens we conduct 1:24,000-scale geologic mapping and lidar analysis to better constrain the location and character of surface faults; detailed mapping in the Wind River Valley was completed just prior to the start of this project. Ages of intrusive rocks are determined with 40 Ar/ 39 Ar geochronology, though all of our samples are Miocene or older. We collect ground based gravity observations (a total of 1,580 new stations) in all of our study areas and ground-based magnetic lines (a total of 93 km) at Mount Baker. These data are combined with existing gravity and aeromagnetic data and used to constrain fault locations and geometry. Two to three cross sections are constructed at each study area using the mapped surface geology and forward-modeling of the gravity and magnetic data; these cross sections form the basis for our updated conceptual models. We collect magnetotelluric surveys at Mount Baker and Mount St. Helens and these data are inverted to form a resistivity model from the surface to about 10 km depth; each model shows conductive zones that can be interpreted as upwelling geothermal fluids. At Mount St. Helens we deploy a passive seismic array and use the newly detected events to refine the location of the Saint Helens seismic zone. We also employ ambient-noise tomography to develop a detailed seismic-velocity model for the study area and use this model to help constrain our cross sections and conceptual model. Based on the new data collected during Phase 2—and our updated models—we develop a campaign of temperature-gradient holes and core analysis to validate our modeling in Phase 3. Drill hole MB76-31 is located near Little Park Creek, 11 km west-southwest of the summit of Mount Baker, and is 1,471 ft deep. About 410 ft of core from the lower portion of the hole—and image logs from ~175 ft below ground surface to the bottom—are collected and analyzed. Water samples are collected and processed for geothermometry. Drill hole MSH17-24 is located along upper Schultz Creek, 16 km north-northeast of Mount St. Helens and has core from 470 ft to the bottom at 1,053 ft. We did not collect image logs due to borehole stability concerns, but water samples are collected and analyzed for geothermometry. Repeat temperature-gradient measurements are made at both sites and thermal conductivity is measured from core samples. At MB76-31, the equilibrated temperature gradient of 64°C/km and calculated heat flow of 141–159 mW/m 2 is more than twice the regional average. Detailed mapping and analysis of the core, coupled with correlation to the image logs, indicates a history of permeability generation consistent with our predictions of high permeability. Because the site has high favorability in the Phase 2 model, we consider the results a positive validation of the modeling. At site MSH17-24, the equilibrated temperature gradient of ~15°C/km and calculated heat flow of 41–43 mW/m 2 are similar to regional. Geochemical analysis of the water samples indicates a meteoric source without any geothermal component. Detailed outcrop-based mapping of fault exposures near the drill site and analysis of image logs from nearby boreholes indicates a history of permeability generation consistent with our predictions. Because the site has low favorability in the Phase 2 model, we consider the results a positive validation of the modeling. Together, the two sites provide a reasonably positive validation of the Phase 2 modeling and should encourage future use of this modeling approach.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Identification of Galaxy–Galaxy Strong Lens Candidates in the DECam Local Volume Exploration Survey Using Machine Learning

We perform a search for galaxy–galaxy strong lens systems using a convolutional neural network (CNN) applied to imaging data from the first public data release of the DECam Local Volume Exploration Survey, which contains ~520 million astronomical sources covering ~4000 deg 2 of the southern sky to a 5σ point–source depth of g = 24.3, r= 23.9, i = 23.3, and z = 22.8 mag. Following the methodology of similar searches using Dark Energy Camera data, we apply color and magnitude cuts to select a catalog of ~11 million extended astronomical sources. After scoring with our CNN, the highest-scoring 50,000 images were visually inspected and assigned a score on a scale from 0 (not a lens) to 3 (very probable lens). We present a list of 581 strong lens candidates, 562 of which are previously unreported. We categorize our candidates using their human-assigned scores, resulting in 55 Grade A candidates, 149 Grade B candidates, and 377 Grade C candidates. We additionally highlight eight potential quadruply lensed quasars from this sample. Due to the location of our search footprint in the northern Galactic cap (b > 10 deg) and southern celestial hemisphere (decl. < 0 deg), our candidate list has little overlap with other existing ground-based searches. Where our search footprint does overlap with other searches, we find a significant number of high-quality candidates that were previously unidentified, indicating a degree of orthogonality in our methodology. We report properties of our candidates including apparent magnitude and Einstein radius estimated from the image separation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

TRACER-Coastal Urban Boundary-Layer Interactions with Convection (TRACER-CUBIC) Field Campaign Report

To better understand the complicated web of processes governing convective cloud life cycle and aerosol-convection interactions, the U.S. Department of Energy (DOE)’s Atmospheric Radiation Measurement (ARM) user facility supported deployment of a variety of advanced atmospheric measurement systems to the greater Houston, Texas, area from 1 October 2021 to 30 September 2022 as part of the Tracking Aerosol Convection Interactions Experiment (TRACER). Houston was selected as a study area because isolated convection and a variety of aerosol conditions are common in this region. This one-year ARM Mobile Facility (AMF) deployment featured a four-month intensive operational period (IOP) during summer 2022 (1 June–30 September). The ARM instrumentation was deployed at three sites along an east-west transect from La Porte, Texas to an ancillary site in a less-polluted rural region southwest of downtown Houston (Figure 1). At the La Porte Site, which is located near the Houston ship channel in an area that experiences significant pollution, the first ARM Mobile Facility (AMF1) was deployed. During the IOP, the ARM tethered balloon system (TBS) operated at the ancillary site. The second-generation C-Band Scanning ARM Precipitation Radar (CSAPR) operated near Pearland, Texas, roughly halfway between the Laporte and ancillary sites. As part of the TRACER- Coastal Urban Boundary-Layer Interactions with Convection (CUBIC) project, three boundary-layer profiling systems) were deployed along a north-south transect spanning from the University of Houston Coastal Center to the Aldine site north of downtown Houston (also blue dot in Figure 1) during the TRACER IOP. These systems included the National Oceanic and Atmospheric Administration (NOAA) National Severe Storms Laboratory CLAMPS2 (C2), which was deployed at the UHCC, the University of Wisconsin SPARC, which was deployed at the ARM CSAPR site near Pearland (orange diamond in middle of map in Figure 1), and the University of Oklahoma CLAMPS1 (C1), which was deployed at Aldine. These three systems have been successfully operated in various field campaigns, providing data sets that collectively offer new insights into atmospheric-boundary-layer (ABL) processes, sea-breeze (SB) circulations, and convection initiation (CI). For the TRACER IOP window, these systems ran continuously between 1 June and 26 September, 2022. Due to commitments to NOAA projects, the Doppler lidar at the UHCC site was not available until 24 June 2022. The CLAMPS and SPARC profiling systems are self-contained platforms that have benefited from several years of development and deployment. Instruments and data processing were maintained remotely, which made the 4-month deployment for the TRACER-CUBIC IOP period possible. The same basic instrument configuration comprises each system: a scanning Doppler wind lidar for flow characterization and passive profiler(s) for characterizing planetary-boundary-layer (PBL) thermodynamic properties. Each platform includes a Halo Streamline Doppler wind lidar, an Atmospheric Emitted Radiance Interferometer (AERI), and a surface meteorology station. CLAMPS1 and CLAMPS2 each also include a microwave radiometer (MWR, Figure 1d). The TRACER-CUBIC hypotheses included (i) Interactions of SB and urban circulations and how they affect the PBL structure in the Houston environment, causing spatially (horizontally and vertically) and temporally highly variable flow patterns, (ii) heat, moisture, and aerosol transport and mixing depend on these flow dynamics, and (iii) an improved understanding of the flow patterns and PBL structure are critical for investigating the processes leading to CI. To test these hypotheses, the project aimed at (i) characterizing SB circulations and their impacts on the diurnal evolution of the structure of the ABL, (ii) studying the evolution of Houston’s complex urban boundary layer, and (iii) identifying effects of urban-induced circulations on pre-convective environments. TRACER-CUBIC observations generally provide good coverage during the summer IOP. Initial screening of the data indicates a good number of cases with bay-breeze (BB) and/or SB signatures, local CI, and interesting boundary-layer features such as strong nocturnal low-level jets (LLJs, Table 1). The numbers listed in rows 3-5 in this table will be further updated as part of ongoing in-depth analyses and systematic identification of local circulations and CI events. More detailed information about the data availability and quality for each instrument is provided in the “readme” files that were submitted to the ARM Data Center along with each archived data sets. These “readme” files also provide instrument descriptions, information about the data collection and processing procedures, data formats, and any additional information relevant for further data analysis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Cation- and lattice-site-selective magnetic depth profiles of ultrathin Fe 3 O 4 (001) films

A detailed understanding of ultrathin film surface properties is crucial for the proper interpretation of spectroscopic, catalytic, and spin-transport data. We present x-ray magnetic circular dichroism (XMCD) and x-ray resonant magnetic reflectivity (XRMR) measurements on ultrathin Fe 3 O 4 films to obtain magnetic depth profiles for the three resonant energies corresponding to the different cation species F e o c t 2 + , F e t e t 3 + , and F e o c t 3 + located on octahedral and tetrahedral sites of the inverse spinel structure of Fe 3 O 4 . By analyzing the XMCD spectrum of Fe 3 O 4 using multiplet calculations, the resonance energy of each cation species can be isolated. Performing XRMR on these three resonant energies yields magnetic depth profiles that each correspond to one specific cation species. The depth profiles of both kinds of Fe 3+ cations reveal a (3.9 ± 1.0) -Å -thick surface layer of enhanced magnetization, which is likely due to an excess of these ions at the expense of the F e o c t 2 + species in the surface region. The magnetically enhanced F e t e t 3 + layer is additionally shifted about 2.9 ± 0.4 Å farther from the surface than the F e o c t 3 + layer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Resolution of Virtual Depth Sectioning from Four-Dimensional Scanning Transmission Electron Microscopy

Abstract One approach to three-dimensional structure determination using the wealth of scattering data in four-dimensional (4D) scanning transmission electron microscopy (STEM) is the parallax method proposed by Ophus et al. (2019. Advanced phase reconstruction methods enabled by 4D scanning transmission electron microscopy, Microsc Microanal25, 10–11), which determines the scattering matrix and uses it to synthesize a virtual depth-sectioning reconstruction of the sample structure. Drawing on an equivalence with a hypothetical confocal imaging mode, we derive contrast transfer and point spread functions for this parallax method applied to weakly scattering objects, showing them identical to earlier depth-sectioning STEM modes when only bright field signal is used, but that improved depth resolution is possible if dark field signal can be used. Through a simulation-based study of doped Si, we show that this depth resolution is preserved for thicker samples, explore the impact of shot noise on the parallax reconstructions, discuss challenges to making use of dark field signal, and identify cases where the interpretation of the parallax reconstruction breaks down.

3D imaging↗
DOE OSTI2024

Surface elevation, snow depth, vegetation height, and color imagery from multiple UAV surveys from 2018 to 2023 across a watershed near Teller road mile marker 27, Seward Peninsula, Alaska

Multiple photogrammetric surveys using an unoccupied aerial system (UAS) were performed across a watershed along Teller road on the Seward Peninsula (Alaska) to investigate terrain, vegetation and snowpack properties. This work was led by the environmental geophysics team from Lawrence Berkeley National Laboratory as part of the Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic). The multiple photogrammetric surveys used mostly a DJI Matrice 210 UAS with a DJI X5S camera, Ground Control Points (GCPs) surveyed with a Real Time Kinematic (RTK)-GPS, and a structure from motion (SfM) technique for photogrammetric reconstruction. Digital Surface elevation Models (DSMs) were inferred at various times of year, including near peak in plant and leaf density (July 19 2017), close to peak in snow depth (April 1 2019 and April 11 2022), and at low plant and leaf density of tall after the first bare-ground date (June 9 2019). Snow and canopy height were obtained by subtracting DSMs from a Digital Terrain Model (DTM) proxy inferred from the June 9 2019 survey. The derived products include ortho-mosaiced RGB map, DSMs, a DTM proxy, snow depth (April 2019 and 2022) and canopy height (July 2017). Unprocessed and processed data products are included in this package (processing levels 0-2). Data and metadata are provided as text (*.csv), image (*.jpg), and GeoTiff (*.tif) formats. This metadata document contains a description of the survey and processing steps and the inferred products.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic) was a 15-year research effort (2012-2027) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy’s Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy’s Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Fine Root and Soil Organic Carbon Depth Distributions are Inversely Related Across Fertility and Rainfall Gradients in Lowland Tropical Forests

Humid tropical forests contain some of the largest soil organic carbon (SOC) stocks on Earth. Much of this SOC is in subsoil, yet variation in the distribution of SOC through the soil profile remains poorly characterized across tropical forests. We used a correlative approach to quantify relationships among depth distributions of SOC, fine root biomass, nutrients and texture to 1 m depths across 43 lowland tropical forests in Panama. The sites span rainfall and soil fertility gradients, and these are largely uncorrelated for these sites. We used fitted b parameters to characterize depth distributions, where b is a numerical index based on an asymptotic relationship, such that larger b values indicate greater concentrations of root biomass or SOC at depth in the profile. Root b values ranged from 0.82 to 0.95 and were best predicted by soil pH and extractable potassium (K) stocks. For example, the three most acidic (pH < 4) and Kpoor (< 20 g K m -2 ) soils contained 76 ± 5% of fine root biomass from 0 to 10 cm depth, while the three least acidic (pH > 6.0) and most K-rich (> 50 g K m -2 ) soils contained only 41 ± 9% of fine root biomass at this depth. Root b and SOC b values were inversely related, such that a large fine root biomass in surface soils corresponded to large SOC stocks in subsoils (50–100 cm). SOC b values were best predicted by soil pH and base cation stocks, with the three most base-poor soils containing 34 ± 8% of SOC from 50 to 100 cm depth, and the three most base-rich soils containing just 9 ± 2% of SOC at this depth. Nutrient depth distributions were not related to Root b or SOC b values. Here, these data show that large surface root biomass stocks are associated with large subsoil C stocks in strongly weathered tropical soils. Further studies are required to evaluate why this occurs, and whether changes in surface root biomass, as may occur with global change, could in turn influence SOC storage in tropical forest subsoils.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

High redshift LBGs from deep broadband imaging for future spectroscopic surveys

Lyman break galaxies (LBGs) are promising probes for clustering measurements at high redshift, z > 2, a region only covered so far by Lyman-α forest measurements. Here, in this paper, we investigate the feasibility of selecting LBGs by exploiting the existence of a strong deficit of flux shortward of the Lyman limit, due to various absorption processes along the line of sight. The target selection relies on deep imaging data from the HSC and CLAUDS surveys in the g, r, z and u bands, respectively, with median depths reaching 27 AB in all bands. The selections were validated by several dedicated spectroscopic observation campaigns with DESI. Visual inspection of spectra has enabled us to develop an automated spectroscopic typing and redshift estimation algorithm specific to LBGs. Based on these data and tools, we assess the efficiency and purity of target selections optimised for different purposes. Selections providing a wide redshift coverage retain 57% of the observed targets after spectroscopic confirmation with DESI, and provide an efficiency for LBGs of 83 ± 3%, for a purity of the selected LBG sample of 90 ± 2%. This would deliver a confirmed LBG density of ~ 620 deg$^{-2}$ in the range 2.3 < z < 3.5 for a r-band limiting magnitude r < 24.2. Selections optimised for high redshift efficiency retain 73% of the observed targets after spectroscopic confirmation, with 89 ± 4% efficiency for 97 ± 2% purity. This would provide a confirmed LBG density of ~ 470 deg$^{-2}$ in the range 2.8 < z < 3.5 for a r-band limiting magnitude r < 24.5.A preliminary study of the LBG sample 3d-clustering properties is also presented and used to estimate the LBG linear bias. A value of b$_{LBG}$ = 3.3 ± 0.2 (stat.) is obtained for a mean redshift of 2.9 and a limiting magnitude in r of 24.2, in agreement with results reported in the literature.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Resolution dependence of the turbulent atmospheric boundary layer in global storm-resolving climate simulations

The current generation of state-of-the-art global climate models are being run at increasingly higher horizontal resolutions, with the goal of resolving organised deep convection explicitly. How a kilometre-scale resolution impacts the representation of the atmospheric boundary layer is, however, not well known. Using statistical analysis on global fields as well as high-frequency data at selected locations, produced with the Integrated Forecasting System (IFS) model for the Next Generation Earth-system Models (nextGEMS) project, we investigate the horizontal resolution dependence of some boundary-layer processes. We find that a change in resolution from 9 to 2.8 km causes no substantial changes to boundary-layer properties and processes at most of the locations studied, although some global changes are detected that indicate circulation changes. Small changes to the boundary-layer depth and structure are found in the Tropics. The short simulation length and lack of data for optimal boundary-layer analysis limits the conclusions, especially in relation to the connection between the boundary layer and the atmosphere general circulation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

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

Constraining the Dynamo Layers in Jupiter and Saturn with Observations and Scaling Laws

The dipole-dominated magnetic fields of Jupiter and Saturn provide evidence for active dynamos operating within their deep interiors, yet the depth of the convecting dynamo layers remains poorly constrained. While magnetic field observations, gravity data, and interior models each provide partial insight, they have not been combined into a single, self-consistent picture of the internal structure. Here, we develop a framework that links observed magnetic field strength with intrinsic heat flux and gravity-constrained interior structure using energy-based dynamo scaling laws. By relating the axial magnetic field strength to the convective power, we infer the radial thickness of the dynamo-generating region for both Jupiter and Saturn. The constants of proportionality in the scaling relations are derived using independent constraints from Earth observations, Jupiter observations, and numerical dynamo simulations. Applied to Jupiter, this framework shows how the inferred dynamo layer thickness is coupled to the outer boundary of the dynamo region. Thinner dynamo layers are predicted when the outer boundary shifts to shallower depths, and no solutions are possible when the outer boundary is less than 73% of Jupiter’s radius. These results constrain plausible geometries for future numerical dynamo simulations. Extending the analysis to Saturn, we find a thick, deep-seated dynamo layer with an outer radius at 42% of the radius to be most plausible. An alternative solution with an inner radius of the dynamo region at 60% of the planetary radius, as suggested by ring seismology models, requires a very thin dynamo layer, occupying only 2%–3% of the total radius.

Geosciences↗
DOE OSTI2021

Integrating New Arctic Plant Functional Types in a Land Surface Model Using Above- and Belowground Field Observations: Modeling Archive

This Modeling Archive is in support of an NGEE Arctic publication "Integrating New Arctic Plant Functional Types in a Land Surface Model Using Above- and Belowground Field Observations". We used field observations of biomass and leaf traits (Salmon et al., 2019b) across a gradient of plant communities on the Seward Peninsula in northwest Alaska (Breen et al., 2020; Iversen et al., 2017) to replace the original ELM configuration with nine, arctic-specific PFTs. Original arctic PFTs in the model were (1) broadleaf deciduous boreal shrub and (2) C3 arctic grass. New PFTs were (1) lichen, (2) bryophyte, (3) evergreen dwarf shrub, (4) deciduous dwarf shrub, (5) deciduous low shrub, (6) deciduous low to tall shrub, (7) alder shrub, (8) forb, and (9) graminoid. New PFTs were parameterized and evaluated using site-level trait and biomass measurements (Salmon et al., 2019b), root C:N ratios and rooting depth patterns (Salmon et al., 2019c), and root biomass (Salmon et al., 2020). This archive contains forcing data and model output from ELM simulations conducted at the Kougarok Hillslope site (Kougarok Road Mile Marker 64). Simulations were conducted for the Kougarok Hillslope site using meteorological driving data from the Scenarios Network for Alaska and Arctic Planning (SNAP) downscaled climate projection dataset using NCAR-CCSM forcing (Bieniek et al., 2020 dataset; Walsh et al., 2018). The archive includes three versions of model simulations representing a progression from the original model configuration to an Arctic-specific configuration (see Output section below). Simulation 1 used grid cell data from a global E3SM configuration, including the fractional area of E3SM PFTs assigned to the grid cell containing the Kougarok Hillslope site in global simulations. This simulation did not distinguish between different plant communities on the Kougarok Hillslope, but instead used a single point simulation to represent the entire area. Simulation 2 used default E3SM PFT definitions combined with adjusted depth to bedrock for plant communities with shallow rocky layers to represent the role of abiotic soil factors in driving site differences (based on depth measurements of Iversen et al., 2019a), and adjusted the relative areas of the broadleaf deciduous boreal shrub and C3 arctic grass PFTs to reflect the observed spatial coverage of shrub and non-shrub PFTs across Kougarok Hillslope plant communities. Simulation 3 used the new arctic PFT definitions and parameterizations based on the vegetation types present at the site, along with the adjusted soil depths used in Simulation 2. Included are *.pdf, *.nc (NetCDF), *tar.gz, *.py, and *.bash files. We would like to thank Mary's Igloo Native Corporation for allowing us to perform this research on their land and giving us the opportunity to share this research with their community. Thanks to Shawn Serbin and Alistair Rogers for helpful discussions. Holly Vander Stel, Joanne Childs, and Stan Wullschleger helped with collection of field data. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

DOE JGI Metagenome Workflow

The DOE Joint Genome Institute (JGI) Metagenome Workflow performs metagenome data processing, including assembly; structural, functional, and taxonomic annotation; and binning of metagenomic data sets that are subsequently included into the Integrated Microbial Genomes and Microbiomes (IMG/M) (I.-M. A. Chen, K. Chu, K. Palaniappan, A. Ratner, et al., Nucleic Acids Res, 49:D751–D763, 2021, https://doi.org/10.1093/nar/gkaa939) comparative analysis system and provided for download via the JGI data portal (https://genome.jgi.doe.gov/portal/). This workflow scales to run on thousands of metagenome samples per year, which can vary by the complexity of microbial communities and sequencing depth. Here, we describe the different tools, databases, and parameters used at different steps of the workflow to help with the interpretation of metagenome data available in IMG and to enable researchers to apply this workflow to their own data. We use 20 publicly available sediment metagenomes to illustrate the computing requirements for the different steps and highlight the typical results of data processing. The workflow modules for read filtering and metagenome assembly are available as a workflow description language (WDL) file (https://code.jgi.doe.gov/BFoster/jgi_meta_wdl). The workflow modules for annotation and binning are provided as a service to the user community at https://img.jgi.doe.gov/submit and require filling out the project and associated metadata descriptions in the Genomes OnLine Database (GOLD) (S. Mukherjee, D. Stamatis, J. Bertsch, G. Ovchinnikova, et al., Nucleic Acids Res, 49:D723–D733, 2021, https://doi.org/10.1093/nar/gkaa983).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Improving precision and accuracy of genetic mapping with genotyping‐by‐sequencing data in outcrossing species

Abstract Genotyping‐by‐sequencing (GBS) is a widely used strategy for obtaining large numbers of genetic markers in model and non‐model organisms. In crop plants, GBS‐derived marker datasets are frequently used to perform quantitative trait locus (QTL) mapping. In some plant species, however, high heterozygosity and complex genome structure mean that researchers must use care in handling GBS data to conduct QTL mapping most effectively. Such outbred crops include most of the perennial grass and tree species used for bioenergy. To identify strategies for increasing accuracy and precision of QTL mapping using GBS data in outbred crops, we conducted an empirical study of SNP‐calling and genetic map‐building pipeline parameters in a Miscanthus sinensis population, and a complementary simulation study to estimate the relationship between genome‐wide error rate, read depth, and marker number. The bioenergy grass Miscanthus is an obligate outcrossing species with a recent (diploidized) whole‐genome duplication. For the study of empirical M. sinensis data, we compared two SNP‐calling methods (one non‐reference‐based and one reference‐based), a series of depth filters (12×, 20×, 30×, and 40×) and two map‐construction methods (i.e., marker ordering: linkage‐only and order‐corrected based on a reference genome). We found that correcting the order of markers on a linkage map by using a high‐quality reference genome improved QTL precision (shorter confidence intervals). For typical GBS datasets of between 1000 and 5000 markers to build a genetic map for biparental populations, a depth filter set at 30× to 40× applied to outbred populations provided a genome‐wide genotype‐calling error rate of less than 1%, improved accuracy of QTL point estimates and minimized type I errors for identifying QTL. Based on these results, we recommend using a reference genome to correct the marker order of genetic maps and a robust genotype depth filter to improve QTL mapping for outbred crops.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Data from Stewart et al. 2026 "Organic Colloid Composition in Variable-Redox Porewaters within a Mountainous Floodplain"

Redox gradients, often driven by changes in sediment moisture levels in porous, heterogeneous groundwater systems, create dynamic conditions that may promote the production and transport of colloids within natural waters. While much research has focused on the inorganic composition of colloids, the organic composition remains less well understood. Organic matter (OM) in colloids may associate with minerals, complex metal ions, and serve as an electron donor for microbial respiration; therefore, its composition is of high interest. We examined the composition of porewater OM along a redox gradient in a riparian soil located along the Slate River in Crested Butte, Colorado, USA as a function of depth (90, 130, 200, and 350 cm below ground surface). All depths were oxic to suboxic, except 200 cm, where the products of iron and sulfate reduction were observed concomitant with an increase in dissolved and/or colloidal OM, pH, alkalinity, and conductivity. We investigated the composition of porewater using correlated scanning transmission X-ray microscopy and transmission electron microscopy. We observed a change in the OM chemistry from carboxylate-rich at the 200 cm depth (where it was also enmeshed with non-crystalline iron) to phenol- and substituted-aromatic-rich at other depths. Radiocarbon dating revealed carbon in porewater at 200 cm was younger than depths above and below. Soil porewater can flow down into the underlaying gravel bed during baseflow conditions, thus we speculate whether riparian porewater could transport iron- and carboxylate-rich organic colloids into surrounding surface waters through the gravel bed conduit. This data package includes geochemical data (.csv) as a function of depth from porewaters in OBJ1. Measurements included are pH, conductivity, TOC, Fe concentration, sulfate concentration, dissolved oxygen. Additionally radiocarbon dating age of carbon in samples is provided. Additionally, this dataset contains metadata files (.csv) describing the methods, file contents, and column headers.

alpine floodplain↗
DOE OSTI2024

Validation Exercise of a Coarse Finite Element Model of Laser Welds

The objective of this project is to validate low-fidelity models of 304L to 304L stainless steel partial-penetration laser welds for thin sheets. Low-fidelity means that the weld is represented by coarsely meshed element blocks. Here, the hexahedral element size is approx imately half the weld penetration depth. The material behavior of the block is represented by a J2 plasticity model with a Voce hardening function. The source of the data used in this work is an extensive experimental study conducted by Sharlotte Kramer (1528) and published in 2015. Figure 1 shows a cross-section of the weld of interest. The nominal thickness of the sheets is 0.063 in. while the target penetration depth of the weld is in the range of 0.028 to 0.032 in., extending about half the sheet thickness. Uniaxial tension tests provided data for calibration of base material and weld models. Results of two validation geometries were also provided. The principal validation geometry is shown in Fig. 2. It consists of a plate specimen with in-plane dimensions 6 in × 2.875 in loaded in tension. A circular plug with a 1.5 in. diameter was cut from the center of the plate and then welded in place. The details of the welding schedule are given. An important assumption is that the welds in the calibration and validation specimens have similar geometric and material properties as those in the validation tests. The task was to first calibrate models for the base material and the welds and then simulate the validation tests until the point of weld first failure.

36 MATERIALS SCIENCE↗