Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “natural soiling”

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 613 records · Page 34

Global Methane Budget 2000–2020

Abstract. Understanding and quantifying the global methane (CH4) budget is important for assessing realistic pathways to mitigate climate change. CH4 is the second most important human-influenced greenhouse gas in terms of climate forcing after carbon dioxide (CO2), and both emissions and atmospheric concentrations of CH4 have continued to increase since 2007 after a temporary pause. The relative importance of CH4 emissions compared to those of CO2 for temperature change is related to its shorter atmospheric lifetime, stronger radiative effect, and acceleration in atmospheric growth rate over the past decade, the causes of which are still debated. Two major challenges in quantifying the factors responsible for the observed atmospheric growth rate arise from diverse, geographically overlapping CH4 sources and from the uncertain magnitude and temporal change in the destruction of CH4 by short-lived and highly variable hydroxyl radicals (OH). To address these challenges, we have established a consortium of multidisciplinary scientists under the umbrella of the Global Carbon Project to improve, synthesise, and update the global CH4 budget regularly and to stimulate new research on the methane cycle. Following Saunois et al. (2016, 2020), we present here the third version of the living review paper dedicated to the decadal CH4 budget, integrating results of top-down CH4 emission estimates (based on in situ and Greenhouse Gases Observing SATellite (GOSAT) atmospheric observations and an ensemble of atmospheric inverse-model results) and bottom-up estimates (based on process-based models for estimating land surface emissions and atmospheric chemistry, inventories of anthropogenic emissions, and data-driven extrapolations). We present a budget for the most recent 2010–2019 calendar decade (the latest period for which full data sets are available), for the previous decade of 2000–2009 and for the year 2020. The revision of the bottom-up budget in this 2025 edition benefits from important progress in estimating inland freshwater emissions, with better counting of emissions from lakes and ponds, reservoirs, and streams and rivers. This budget also reduces double counting across freshwater and wetland emissions and, for the first time, includes an estimate of the potential double counting that may exist (average of 23 Tg CH4 yr−1). Bottom-up approaches show that the combined wetland and inland freshwater emissions average 248 [159–369] Tg CH4 yr−1 for the 2010–2019 decade. Natural fluxes are perturbed by human activities through climate, eutrophication, and land use. In this budget, we also estimate, for the first time, this anthropogenic component contributing to wetland and inland freshwater emissions. Newly available gridded products also allowed us to derive an almost complete latitudinal and regional budget based on bottom-up approaches. For the 2010–2019 decade, global CH4 emissions are estimated by atmospheric inversions (top-down) to be 575 Tg CH4 yr−1 (range 553–586, corresponding to the minimum and maximum estimates of the model ensemble). Of this amount, 369 Tg CH4 yr−1 or ∼ 65 % is attributed to direct anthropogenic sources in the fossil, agriculture, and waste and anthropogenic biomass burning (range 350–391 Tg CH4 yr−1 or 63 %–68 %). For the 2000–2009 period, the atmospheric inversions give a slightly lower total emission than for 2010–2019, by 32 Tg CH4 yr−1 (range 9–40). The 2020 emission rate is the highest of the period and reaches 608 Tg CH4 yr−1 (range 581–627), which is 12 % higher than the average emissions in the 2000s. Since 2012, global direct anthropogenic CH4 emission trends have been tracking scenarios that assume no or minimal climate mitigation policies proposed by the Intergovernmental Panel on Climate Change (shared socio-economic pathways SSP5 and SSP3). Bottom-up methods suggest 16 % (94 Tg CH4 yr−1) larger global emissions (669 Tg CH4 yr−1, range 512–849) than top-down inversion methods for the 2010–2019 period. The discrepancy between the bottom-up and the top-down budgets has been greatly reduced compared to the previous differences (167 and 156 Tg CH4 yr−1 in Saunois et al. (2016, 2020) respectively), and for the first time uncertainties in bottom-up and top-down budgets overlap. Although differences have been reduced between inversions and bottom-up, the most important source of uncertainty in the global CH4 budget is still attributable to natural emissions, especially those from wetlands and inland freshwaters. The tropospheric loss of methane, as the main contributor to methane lifetime, has been estimated at 563 [510–663] Tg CH4 yr−1 based on chemistry–climate models. These values are slightly larger than for 2000–2009 due to the impact of the rise in atmospheric methane and remaining large uncertainty (∼ 25 %). The total sink of CH4 is estimated at 633 [507–796] Tg CH4 yr−1 by the bottom-up approaches and at 554 [550–567] Tg CH4 yr−1 by top-down approaches. However, most of the top-down models use the same OH distribution, which introduces less uncertainty to the global budget than is likely justified. For 2010–2019, agriculture and waste contributed an estimated 228 [213–242] Tg CH4 yr−1 in the top-down budget and 211 [195–231] Tg CH4 yr−1 in the bottom-up budget. Fossil fuel emissions contributed 115 [100–124] Tg CH4 yr−1 in the top-down budget and 120 [117–125] Tg CH4 yr−1 in the bottom-up budget. Biomass and biofuel burning contributed 27 [26–27] Tg CH4 yr−1 in the top-down budget and 28 [21–39] Tg CH4 yr−1 in the bottom-up budget. We identify five major priorities for improving the CH4 budget: (i) producing a global, high-resolution map of water-saturated soils and inundated areas emitting CH4 based on a robust classification of different types of emitting ecosystems; (ii) further development of process-based models for inland-water emissions; (iii) intensification of CH4 observations at local (e.g. FLUXNET-CH4 measurements, urban-scale monitoring, satellite imagery with pointing capabilities) to regional scales (surface networks and global remote sensing measurements from satellites) to constrain both bottom-up models and atmospheric inversions; (iv) improvements of transport models and the representation of photochemical sinks in top-down inversions; and (v) integration of 3D variational inversion systems using isotopic and/or co-emitted species such as ethane as well as information in the bottom-up inventories on anthropogenic super-emitters detected by remote sensing (mainly oil and gas sector but also coal, agriculture, and landfills) to improve source partitioning. The data presented here can be downloaded from https://doi.org/10.18160/GKQ9-2RHT (Martinez et al., 2024).

54 ENVIRONMENTAL SCIENCES↗

Improving CLM5.0 Biomass and Carbon Exchange Across the Western United States Using a Data Assimilation System

The Western United States is dominated by natural lands that play a critical role for carbon balance, water quality, and timber reserves. This region is also particularly vulnerable to forest mortality from drought, insect attack, and wildfires, thus requiring constant monitoring to assess ecosystem health. Carbon monitoring techniques are challenged by the complex mountainous terrain, thus there is an opportunity for data assimilation systems that combine land surface models and satellite-derived observations to provide improved carbon monitoring. Here, we use the Data Assimilation Research Testbed to adjust the Community Land Model (CLM5.0) with remotely sensed observations of leaf area and above-ground biomass. The adjusted simulation significantly reduced the above-ground biomass and leaf area, leading to a reduction in both photosynthesis and respiration fluxes. The reduction in the carbon fluxes mostly offset, thus both the adjusted and free simulation projected a weak carbon sink to the land. This result differed from a separate observation-constrained model (FLUXCOM) that projected strong carbon uptake to the land. Simulation diagnostics suggested water limitation had an important influence upon the magnitude and spatial pattern of carbon uptake through photosynthesis. We recommend that additional observations important for water cycling (e.g., snow water equivalent, land surface temperature) be included to improve the veracity of the spatial pattern in carbon uptake. Furthermore, the assimilation system should be enhanced to maximize the number of the simulated state variables that are adjusted, especially those related to the recommended observed quantities including water cycling and soil carbon.

54 ENVIRONMENTAL SCIENCES↗

Replacing liquid fossil fuels and hydrocarbon chemical feedstocks with liquid biofuels from large-scale nuclear biorefineries

Liquid fossil fuels (1) enable transportation and (2) provide energy for mobile work platforms and (3) supply dispatchable energy to highly variable demand (seasonal heating and peak electricity). We describe a system to replace liquid fossil fuels with drop-in biofuels including gasoline, diesel and jet fuel. Because growing biomass removes carbon dioxide from the air, there is no net addition of carbon dioxide to the atmosphere from burning biofuels. In addition, with proper management, biofuel systems can sequester large quantities of carbon as soil organic matter, improving soil fertility and providing other environmental services. In the United States liquid biofuels can potentially replace all liquid fossil fuels. The required system has two key features. First, the heat and hydrogen for conversion of biomass into high-quality liquid fuels is provided by external low-carbon energy sources--nuclear energy or fossil fuels with carbon capture and sequestration. Using external energy inputs can almost double the energy content of the liquid fuel per unit of biomass feedstock by fully converting the carbon in biomass into a hydrocarbon fuel. Second, competing effectively with fossil fuels requires very large biorefineries—the equivalent of a 250,000 barrel per day oil refinery. This requires commercializing methods for converting local biomass into high-density storable feedstocks that can be economically shipped to large-scale biorefineries.

09 - BIOMASS FUELS↗

Data from TropiRoot 1.0 database: tropical root characteristics across environments

TropiRoot 1.0 is a new tropical root database with root characteristics across environment gradients. It has data extracted from 104 new sources, resulting in more than 8000 rows of data (either species or community data). Most of the data in TropiRoot 1.0 includes root characteristics such as root biomass, morphology, root dynamics, mass fraction, architecture, anatomy, physiology and root chemistry. This initiative represents an approximately 30% increase in the currently available data for tropical roots in the Fine Root Ecology Database (FRED). TropiRoot 1.0, contains root characteristics from 25 different countries where seven are located in Asia, six in South America, five in Central America and the Caribbean, four in Africa, two in North America, and 1 in Oceania. Due to the volume of data, when ancillary data was available, including soil data, these data was either extracted and included in the database or their availability was recorded in an additional column. Multiple contributors checked the entries for outliers during the collation process to ensure data quality. For text-based observations, we examined all cells to ensure that their content relates to their specific categories. For numerical observations, we ordered each numerical value from least to greatest and plotted the values, checking apparent outliers against the data in their respective sources and correcting or removing incorrect or impossible values. Some data (soil and aboveground) have different columns for the same variable presented in different units, including originally published units, but root characteristics data had units converted to match the ones reported in FRED. By filling a gap from global databases, TropiRoot 1.0 expands our knowledge of otherwise so far underrepresented regions, and our ability to assess global trends. This advancement can be used to improve tropical forest representation in vegetation models.

54 ENVIRONMENTAL SCIENCES↗

Mass spectrometry and spectroscopic characterization of a tetrameric photosystem I supercomplex from Leptolyngbya ohadii , a desiccation-tolerant cyanobacterium

Cyanobacteria inhabiting desert biological soil crusts face the harsh conditions of the desert. They evolved a suite of strategies toward desiccation-hydration cycles mixed with high light irradiations, etc. In this study we purified and characterized the structure and function of Photosystem I (PSI) from Leptolyngbya ohadii, a desiccation-tolerant desert cyanobacterium. We discovered that PSI forms tetrameric (PSI-Tet) aggregate. We investigated it by using sucrose density gradient centrifugation, clear native PAGE, high performance liquid chromatography, mass spectrometry (MS), time-resolved fluorescence (TRF) and time-resolved transient absorption (TA) spectroscopy. MS analysis identified the presence of two PsaB and two PsaL proteins in PSI-Tet and uniquely revealed that PsaLs are N-terminally acetylated in contrast to non-modified PsaL in the trimeric PSI from Synechocystis sp. PCC 6803. Chlorophyll (Chl) a fluorescence decay profiles of the PSI-Tet performed at 77 K revealed two emission bands at ∼690 nm and 725 nm with the former appearing only at early delay time. The main fluorescence emission peak, associated with emission from the low energy Chls a, decays within a few nanoseconds. TA studies demonstrated that the 725 nm emission band is associated with low energy Chls a with absorption band clearly resolved at ∼710 nm at 77 K. As a result, our work suggests that the heterogenous composition of PsaBs and PsaL in PSI-Tet is related with the adaptation mechanisms needed to cope with stressful conditions under which this bacterium naturally grows.

59 BASIC BIOLOGICAL SCIENCES↗

STILGAR End-of-Project Report

The Subsurface Tunnel Imaging LeveraGed by Analysis of Rayleigh wave ellipticity (STILGAR) project demonstrated an integrated geophysical approach for detecting, locating, and characterizing underground structural changes using dense seismic arrays and advanced inversion techniques. Field campaigns were conducted at two operational mines—the Redmond salt mine (Utah) and Graymont Pleasant Gap limestone mine (Pennsylvania)—providing real-world testbeds for monitoring anthropogenic subsurface activity. At the Redmond salt mine, seismic interferometry combined with back-projection inversion successfully identified continuous, low-amplitude signals from mining operations. The approach differentiated stationary from migrating anthropogenic sources, captured daily operational cycles, and validated the potential of passive seismic monitoring for remote detection of underground activity. At the Graymont Pleasant Gap mine, two dense seismic deployments in the spring and fall of 2023 generated over 4 TB of high-resolution data. Key outcomes included the relocation of 199 underground and 8 surface explosions with accuracies within tens of meters and the development of a 3D P-wave velocity model using the triple-difference tomography algorithm (tomoTD) that resolved major structural features such as the mine entrance, low-velocity tunnels, and roof-collapse areas. Ambient noise cross-correlation and back-projection analyses revealed persistent sources linked to ongoing mining activity, whereas horizontal-to-vertical spectral ratio (HVSR) and ellipticity studies confirmed stable site responses across seasons and identified soil thickness trends consistent with regional erosional and depositional processes. Checkerboard and sensitivity tests further validated the robustness of the tomographic results. Overall, the findings emphasize that although significant progress has been made in subsurface imaging, further work is needed to enhance the detection and localization of underground structures. Accurate imaging requires higher frequencies, yet anthropogenic sources tend to dominate the seismic record at those frequencies, and high-frequency surface waves are affected by higher modes that complicate interpretation. The improved detection and localization of human-induced signals enabled detailed temporal and spatial mapping of daily mine operations, demonstrating the feasibility of continuous anthropogenic source monitoring. Sensitivity to signals from nontraditional sources, such as fan operations, highlights the broader applicability of this approach to other industrial environments in which continuous and impulsive signals are present. The field campaigns produced a substantial volume of high-quality seismic data, supporting the development and testing of new methods for seismic source characterization and subsurface imaging. Future deployments should include sensors capable of recording lower frequencies to probe deeper structures, increase bandwidth to enhance resolution and sensitivity to both shallow and deep targets, and collect additional large-scale datasets to refine imaging and source characterization techniques. Moreover, conducting 3D modeling studies of seismic wavefields at higher frequencies will provide a better understanding of wave scattering and cavity–wavefield interactions in complex underground environments. In conclusion, the STILGAR project demonstrated that integrated seismic monitoring can effectively characterize underground operations, capturing both natural and anthropogenic signals. The approaches developed provide a foundation for improved detection, localization, and imaging of subsurface structures and are directly transferable to broader industrial monitoring applications.

58 GEOSCIENCES↗

2017 Meander C sediment characterization from the East River, Colorado

This dataset includes characterization data collected on sediment samples from Meander C as part of the Watershed Function Scientific Focus Area (SFA) located in the Upper Colorado River Basin. The data were collected to investigate the nature of mineral-organic associations across the meander transect and with depth. Samples were collected in 2017, and analyzed using sequential extractions (water, pyrophosphate, HCl, and dithionite) and density fractionation. Data collected include: texture, pH, root biomass, and metals (Fe, Al, Ca, Si) of bulk sediment samples; C and N of heavy fraction; and extractable OC, Fe, and Al, and organic matter characterization using high-resolution mass spectrometry (FT-ICR-MS) of sequential extractions. Sample locations are included in the ‘sample_metadata.xlsx” file. Sediment characterization data is included in the accompanying csv files: bulk and density fractions (MCB2017_density-and-bulk-fractions.csv), sequential extractions (MCB2017_sequential-extractions.csv), and FT-ICR-MS data from sequential extractions in three files to be read by the R package ftmsRanalysis (MCB2017_FTICRMS-sequential-edata.csv, MCB2017_FTICRMS-sequential-emeta.csv, and MCB2017_FTICRMS-sequential-fdata.csv). Sediment was collected in triplicate; ID’s refer to individual samples.

54 ENVIRONMENTAL SCIENCES↗

Dominance of Particulate Mercury in Stream Transport and Rapid Watershed Recovery from Wildfires in Northern California, USA

Frequency and intensity of wildfires are expected to increase due to climate change, especially in areas with a long summer drought. Forests are a major sink for the global pollutant mercury (Hg), and fluvial transport of Hg from recently burned watersheds has not been widely investigated. Here, we examined two years of fluvial transport of Hg and its speciation (total Hg, methyl-Hg, particulate, and dissolved forms) under storm events and baseflow in two recently burned watersheds with different burned proportions and one nonburned reference watershed in the Coastal Ranges of northern California. We examined postfire storm-event transport of Hg and its methylated form (methyl-Hg), addressed the importance of the “initial runoff pulse” to postfire Hg fluvial transport and its predominant association with suspended solids, and elucidated potential sources of Hg exports from the burned landscapes using geochemical indicators, which suggested that ash materials were likely the significant sources of particulates in the first high-flow season postfire but not subsequently. The maximum total suspended solid and total Hg levels in the “first pulse” at the severely burned watershed were 442 and 46 times higher, respectively, than those at the reference watershed. Stream suspended solid and Hg levels declined substantially in the burned watersheds after just a few months of rainfall likely due to the rapid regrowth of vegetation commonly observed in postfire landscapes, implying that the wildfire effects on immediate Hg inputs from the burned landscape are at most transient in nature.

54 ENVIRONMENTAL SCIENCES↗

Data for Aboveground rather than belowground productivity drives variability in Miscanthus x giganteus net primary productivity

This dataset contains the data used for the publication “Aboveground rather than belowground productivity drives variability in Miscanthus x giganteus net primary productivity”. This dataset contains Miscanthus x giganteus biomass, carbon, and nitrogen tissue data for aboveground and belowground plant parts collected in 2021 for three different sites in Iowa with three different nitrogen application rates. Data at the Iowa sites were collected via biometric hand harvesting, belowground excavations, and soil coring both in-clump and beside-clump. Data were collected at two collection timepoints to calculate the contributions of belowground parts to Miscanthus x giganteus net primary productivity. This dataset also includes Miscanthus x giganteus and Switchgrass soil coring and excavation data collected in 2012 at the University of Illinois Urbana Champaign Energy Farm.

Belowground Biomass↗

Air temperature and precipitation constraining the modelled wetland methane emissions in a boreal region in northern Europe

Wetland methane responses to temperature and precipitation are studied in a boreal wetland-rich region in northern Europe using ecosystem process models. Six ecosystem models (JSBACH-HIMMELI, LPX-Bern, LPJ-GUESS, JULES, CLM4.5, and CLM5) are compared to multi-model means of ecosystem models and atmospheric inversions from the Global Carbon Project and upscaled eddy covariance flux results for their temperature and precipitation responses and seasonal cycles of the regional fluxes. Two models with contrasting response patterns, LPX-Bern and JSBACH-HIMMELI, are used as priors in atmospheric inversions with Carbon Tracker Europe–CH4 (CTE-CH4) in order to find out how the assimilation of atmospheric concentration data changes the flux estimates and how this alters the interpretation of the flux responses to temperature and precipitation. Inversion moves wetland emissions of both models towards co-limitation by temperature and precipitation. Between 2000 and 2018, periods of high temperature and/or high precipitation often resulted in increased emissions. However, the dry summer of 2018 did not result in increased emissions despite the high temperatures. The process models show strong temperature and strong precipitation responses for the region (51 %–91 % of the variance explained by both). The month with the highest emissions varies from May to September among the models. However, multi-model means, inversions, and upscaled eddy covariance flux observations agree on the month of maximum emissions and are co-limited by temperature and precipitation. The setup of different emission components (peatland emissions, mineral land fluxes) has an important role in building up the response patterns. Considering the significant differences among the models, it is essential to pay more attention to the regional representation of wet and dry mineral soils and periodic flooding which contribute to the seasonality and magnitude of methane fluxes. The realistic representation of temperature dependence of the peat soil fluxes is also important. Furthermore, it is important to use process-based descriptions for both mineral and peat soil fluxes to simulate the flux responses to climate drivers.

54 ENVIRONMENTAL SCIENCES↗

Monthly averages of ED2 model simulations initialized with airborne lidar structure, Jan 1981-Dec 2018, Brazilian Amazon

Deforestation and forest degradation (selective logging, fires, fragmentation) have impacted nearly 40% of the original extent of the Brazilian Amazon, and have markedly impacted forest structure across the region. To date, few studies analysed how shifts in forest structure from degradation influence the forest sensitivity to climate extremes, because of the complex interactions between forest structure and micro-environmental conditions. To address this knowledge gap, we carried out a series of simulations across the Brazilian Amazon using the Ecosystem Demography Model (ED2), using observed forest structure derived from 541 airborne lidar transects (375 ha each) and two scenarios representing forest recovery and expansion of degradation to investigate how shifts in forest structure impact ecosystem function under near-average and extreme climate conditions, as part of the manuscript Longo et al 2025 "Degradation and Deforestation Increase the Sensitivity of the Amazon Forest to Climate Extremes". This dataset provides the output results from the ED2 model simulations for the three simulations at monthly time scales, in NetCDF format. For all simulations, we used bias-corrected hourly reanalyses (WFDE5) for most meteorological drivers, except for precipitation, which was obtained from CHIRPS. The meteorological drivers used in the study span 38 years (Jan 1981–Dec 2018). The output results correspond to the last 38 years of simulation (one full cycle of meteorological drivers), in which ED2 simulations used static stand structure (i.e., the forest structure was held constant). The following files are provided:ED2_emean_Global_R004_BrAmaz_s1c0t0l0f0.nc. This corresponds to the Control simulation. The forest structure was obtained from the airborne lidar.ED2_emean_Global_R005_BrAmaz_s1c0t1l1f0.nc. This corresponds to the Degraded simulation. The forest structure was obtained from a spin-up simulation initialized with airborne lidar and a scenario that expanded deforestation and selective logging across the Amazon.ED2_emean_Global_R006_BrAmaz_s1c0t1l0f0.nc. This corresponds to the Recovery simulation. The forest structure was obtained from a spin-up simulation initialized with airborne lidar and a scenario that completely halted deforestation and degradation, allowing degraded forests to recover for 38 years.We also provide file ED2_zones_R004_BrAmaz_s1c0t0l0f0.nc, which classifies each grid cell into zones used in the reference manuscript: 1: Southeast. 2: South. 3: West. 4: Central. 5: Northeast. 6: North. 7: Northwest". Index 0 corresponds to grid cells excluded from sub-region analyses because they were dominated by flooded forests, deforestation, and naturally non-forest vegetation.

54 ENVIRONMENTAL SCIENCES↗

Interrogating organo-mineral interactions by Py-GC/MS and ToF-SIMS

Soil organic matter (SOM) is crucial to soil fertility and agricultural sustainability in the face of climate change, and organo-mineral interactions are a key mechanism for stabilizing SOM. The stability of mineral-associated organic matter (MAOM) can be approximated by the strength of bonds between biomolecules and mineral surfaces. We hypothesized that factors such as mineral reactivity and solution pH would impact bonding strength. To test these factors, sorption isotherms were performed on minerals differing in particle size, morphology, and reactivity (quartz, goethite, gibbsite, kaolinite, montmorillonite), under two pH conditions (4 and 6), and across a range of tannic acid concentration (0-1400 ppm). The quantity and quality of organo-mineral bonds was assessed by quantifying post-sorption solution compounds using an absorbance plate reader and assessing the thermal energy needed to break bonds using evolved gas analysis (EGA) on a pyrolysis gas chromatograph-mass spectrometer (Py-GC/MS). Organo-mineral surfaces were also interrogated using Time-of-flight secondary ion mass spectrometry (ToF-SIMS). Sorption experiments revealed an asymptotic relationship between sorbed tannic acid and increasing initial concentrations which suggests our concentration range saturated available mineral surfaces. Results indicated that mineral identity, chiefly surface area and native pH, impacted the amount of compound sorbed to mineral surfaces. Evolved gas analyses of organo-mineral samples uncovered increased thermal stability of sorbed compounds compared to our tannic acid standard particularly on the iron oxide, goethite. ToF-SIMS revealed that goethite surface ions were covered by increasing concentrations of adsorbed tannic acid with the replacement of chlorine and sulfate ions by tannic acid related ions. Our results indicate that Py-GC/MS and ToF-SIMS are useful techniques to determine the nature and strength of organo-mineral bonds and are steppingstones to understanding controls on carbon sequestration in terrestrial systems.

58 GEOSCIENCES↗

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↗

Morphological and crystallographic controls in the replacement of calcite and aragonite by cerussite and otavite

Mineral replacement reactions are essential solid-fluid interactions in natural and industrial processes including metasomatism, diagenesis, and metamorphism, and sequestration of metal toxins from polluted water. In this paper, we explore the morphological evolution of aragonite and calcite (polymorphs of CaCO 3 ) during their replacement by cerussite (PbCO 3 , an isomorph of aragonite) and otavite (CdCO 3 , an isomorph of calcite) in acidic aqueous solutions (with initial pH 3 to 4). Observations using scanning electron and synchrotron-based transmission X-ray microscopies reveal the formation of pseudomorphic shells (~5 to 10 µm thick) of cerussite and otavite when Pb 2+ is reacted with either calcite or aragonite and when Cd 2+ is reacted with aragonite. The formation of pore space at the substrate-precipitate interface and within the precipitate is found to be key in promoting chemical exchange between the dissolving phase and the solution. In contrast, minimal reactivity is observed when Cd 2+ is reacted with calcite. These results demonstrate that the lack of substrate-precipitate epitaxy (either because of inconsistent symmetries or due to large lattice constant mismatch) is a critical factor that enables replacement of primary CaCO 3 minerals by Pb- or Cd-bearing secondary minerals while the presence of epitaxial relationships inhibit those reactions. When calcite and aragonite coexist, the replacement of aragonite by cerussite and otavite occurs preferentially over calcite. We postulate that the relative reactivity of calcite and aragonite to Pb 2+ and Cd 2+ can be determined by interplay between morphological and crystallographic properties of CaCO 3 materials as well as solution chemistry. Our results demonstrate that coupled dissolution-precipitation processes of carbonate minerals play an important role in controlling mobility and bioavailability of lead and cadmium in environments where acidic metal-contaminated water interacts with carbonate-rich rocks and soils.

58 GEOSCIENCES↗

Invasive Plant Species Management Plan for Los Alamos National Laboratory (Rev.3)

Native species are plants and animals that continually occupy a natural range without direct or indirect introduction and/or care by humans. They are adapted to the environmental conditions and processes of the ecosystem in which they reside. A species introduced into a novel ecosystem can either exploit that ecosystem and thrive or be unable to survive in that ecosystem (Hobbs et al. 2006). Alien or non-native species are species that are intentionally or accidentally introduced into a novel ecosystem and are capable of living and propagating within the physical parameters of that ecosystem. Invasive species is a species that is non-native (or alien) to the ecosystem under consideration and a noxious species are those whose introduction causes or is likely to cause economic or environmental harm. The term invasive species is applicable to plants and animals alike; however, this invasive species management plan currently focuses on invasive vegetation. Invasive plant species are usually capable of rapid colonization of disturbed ground, such as after a change in wildfire regime intensity and frequency (Reilly et al. 2020) or anthropogenic ground disturbance (Burke and Grime 1996; Hobbs and Huenneke 1992). Climate change facilitates the spread and establishment of many alien species and creates new opportunities for them to become invasive (Turbelin and Catford 2021). As climate change impacts increase in the coming decades, we may see in an increase in invasive species establishment. See a list of definitions of terms pertinent to this document in Appendix A: Definitions. Los Alamos National Laboratory (LANL) hosts populations of non-native and invasive species all typical of the northern New Mexico region (Martin 2004; NMDA 2020). By implementing an invasive species management plan, LANL will have readily accessible management strategies for invasive species that are found on-site. The benefits of managing invasive species include a decrease in wildland fire risks, an increase in soil productivity, an increase in (productive or beneficial) wildlife habitat, an increase in water quantity and quality, and the restoration of impacted areas (Burke and Grime 1996; Hobbs and Huenneke 1992; D’Antonio and Hobbie 2005; MacDougall et al. 2013; Reilly et al. 2020). The aforementioned benefits from invasive species management directly enable the LANL mission by ensuring compliance requirements are met and site-wide programs, such as the Vulnerability Assessment and Resilience Plan, are supported for a mutually beneficial outcome. An example for the LANL site specifically, controlling non-native annual plant species, for example, could reduce the costs associated with stabilizing soils during the stormwater pollution prevention compliance process. Roadway and utility right-of-way areas are another place where an integrated vegetation management strategy would promote low growing perennial plants in a way that is mutually beneficial to habitat and the institution through lowered maintenance costs. There is also an economic benefit to managing invasive species. In one nationwide study, invasive plants had an estimated impact cost of $190.45 billion (Fantel-Lepczyk et al. 2022). Investing in preventative measures and surveillance could help to offset future control and management costs of invasive species that have the potential to become established. The State of New Mexico has developed plant species lists and recommendations through the New Mexico Department of Agriculture’s Noxious Weeds Management Act, Article 7D (NM Statute § 76-7D-4 2021); however, effective invasive species management must rely on local knowledge of the site and region. Los Alamos County (LAC) has already compiled a target invasive plant species list and species-specific management objectives (Martin 2004). Therefore, the New Mexico Noxious Weed List and the LAC invasive plant species list, as well as management objectives from those documents, are integrated into LANL’s invasive plant species management plan.

54 ENVIRONMENTAL SCIENCES↗

Supporting data for "Frozen flume experiments indicate rapid permafrost riverbank erosion depends on bank roughness"

This dataset contains supporting files detailing five frozen flume experiments conducted at the Caltech Earth Surface Dynamics Laboratory to investigate rates of ablation-limited permafrost riverbank erosion under controlled conditions. Water flowed past a bank of saturated, frozen sand and ice and gradually eroded the bank by thawing pore ice and immediately entraining sand and washing it downstream. Experiments were scaled for flow hydraulics and heat transfer allowing comparisons between our results and natural permafrost riverbanks. For each experiment, we measured the initial and final sand bank topography using a Keyence laser scanner, water surface slope at 3-min intervals throughout the experiment using a Massa sonar scanner, bank erosion using 10-sec overhead timelapse imagery taken by an overhead camera, water and bank temperature using thermistors frozen into the sand bank and sampling at 2 Hz, and water discharge using an in-line flow meter. We include calibration data for the carriage (engineered by the Saint Anthony Falls Laboratory) used to make sonar and laser topography measurements. We also include calibration data for temperature sensors, water discharge measurements, and images of a regular grid placed in the flume to align overhead camera images with the carriage datum. Grain size analysis for the channel bed (gravel) was produced using a pebble count and bank sand was measured using a Camsizer X2. In addition to the five frozen experiments, we include sonar scans of water surface slope and Keyence scans of bed and bank topography for calibration experiments ran with an immobile gravel bank and bed.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of the Large-Scale and Regional Climatic Response Across North Africa to Natural Variability in Oceanic Modes and Terrestrial Vegetation Among the CMIP5 Models

Hydrologic variability is a serious threat to the poverty-stricken regions of North Africa. Meanwhile, the scientific community struggles to attribute these extreme climatic episodes to specific oceanic and land drivers. Prior modeling studies have not assessed simulated feedbacks over North Africa against an established observational benchmark, so their results are considered as untested, model-specific findings. The Coupled Model Intercomparison Project Phase Five (CMIP5) archive represents the state-of-the-art in climate projections, with most CMIP5 models now containing interactive vegetation phenology; however, there have been no studies to date of their simulated vegetation feedbacks. The capability of CMIP5 models at accurately simulating the forcing of both sea-surface temperature (SST) and regional leaf area index (LAI) anomalies on North African climate needs to be a key consideration in assessing the models’ overall credibility and determining appropriate model weighting for developing climate projections. The Generalized Equilibrium Feedback Assessment (GEFA) is a promising statistical method, which can be applied to either model output or observations, for isolating the local and remote impacts of individual oceanic or terrestrial forcings on regional climate. We performed a combined observational and modeling assessment of land-ocean-atmosphere interactions across the distinct ecological and moisture gradients of North Africa. We evaluated and demonstrated the reliability of the GEFA statistical method over North Africa using the Community Earth System Model (CESM), applied GEFA to observational data to quantify the observed forcing of ocean basin SST anomalies and regional LAI anomalies on North African climate, evaluated the CMIP5 models’ performance in terms of representing these key observed feedbacks, and formulated CMIP5 feedback performance metrics for weighting North African climate projections. Our study represented the first attempt to separate the observed roles of oceanic and vegetation feedbacks across North Africa, the first systematic assessment and intercomparison of land-ocean-atmosphere feedbacks in CMIP5, and the first exploration of vegetation feedbacks among CMIP5 models. The following work was accomplished. (1) The team generated 14 scientific publications and 30 presentations based on the DOE-funded research. (2) A stepwise version of the GEFA statistical method was developed in which unimportant forcings were dropped to increase the reliability of the results. (3) SGEFA was successfully validated through experiments with CESM in terms of its ability to isolate the atmospheric responses to individual oceanic or land forcings. (4) Observational evidence was revealed for the Sahel’s positive vegetation-rainfall feedback on the seasonal to interannual time scale, and it was attributed to a moisture recycling mechanism rather than an albedo mechanism. (5) An approach with SGEFA was developed in which the individual contributions of soil moisture versus vegetation forcings could be separated, revealing that the former forcing outweighed the latter for sub-Saharan Africa. (6) Tropical ocean temperatures were found to be key regulators of pan-tropical vegetation variability, especially for arid and semi-arid regions, including sub-Saharan Africa. (7) The CMIP5 models largely underestimated the importance of land feedbacks across the Sahel. (8) The general consensus among CMIP5 models indicates, for the late 21st century, a diminished seasonal predictability of sub-Saharan African regional climate and an elevated role of the land surface compared to oceanic drivers in regulating regional climate variability. (9) Land-ocean-atmosphere interactions were demonstrated to be key contributors to the seasonal predictability of African wildfire activity. (10) The El Djouf was determined to be of greater importance than the Bodélé depression in terms of providing trans-Atlantic dust transport to the Americas.

54 ENVIRONMENTAL SCIENCES↗