Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Earth surface”

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 181 records · Page 10

Harmonization of global land use change and management for the period 850–2100 (LUH2) for CMIP6

Abstract. Human land use activities have resulted in large changes to the biogeochemical and biophysical properties of the Earth's surface, with consequences for climate and other ecosystem services. In the future, land use activities are likely to expand and/or intensify further to meet growing demands for food, fiber, and energy. As part of the World Climate Research Program Coupled Model Intercomparison Project (CMIP6), the international community has developed the next generation of advanced Earth system models (ESMs) to estimate the combined effects of human activities (e.g., land use and fossil fuel emissions) on the carbon–climate system. A new set of historical data based on the History of the Global Environment database (HYDE), and multiple alternative scenarios of the future (2015–2100) from Integrated Assessment Model (IAM) teams, is required as input for these models. With most ESM simulations for CMIP6 now completed, it is important to document the land use patterns used by those simulations. Here we present results from the Land-Use Harmonization 2 (LUH2) project, which smoothly connects updated historical reconstructions of land use with eight new future projections in the format required for ESMs. The harmonization strategy estimates the fractional land use patterns, underlying land use transitions, key agricultural management information, and resulting secondary lands annually, while minimizing the differences between the end of the historical reconstruction and IAM initial conditions and preserving changes depicted by the IAMs in the future. The new approach builds on a similar effort from CMIP5 and is now provided at higher resolution (0.25∘×0.25∘) over a longer time domain (850–2100, with extensions to 2300) with more detail (including multiple crop and pasture types and associated management practices) using more input datasets (including Landsat remote sensing data) and updated algorithms (wood harvest and shifting cultivation); it is assessed via a new diagnostic package. The new LUH2 products contain > 50 times the information content of the datasets used in CMIP5 and are designed to enable new and improved estimates of the combined effects of land use on the global carbon–climate system.

54 ENVIRONMENTAL SCIENCES↗

NGEE Arctic Plant Traits: Fine Roots, Kougarok Road Mile Marker 64, Seward Peninsula, Alaska, 2016

Soil cores were collected from twelve vegetation biomass plots at the Kougarok hillslope in late July of 2016. The sampled plots were located across six ecotypes present at this site (n=2 replicates per ecotype). Soil cores were separated into depth intervals in the field and frozen for transport to Oak Ridge National Laboratory. Data package contains one *.csv and one *.pdf.Laboratory sample processing of these soils generated depth-specific data on soil properties as well as fine root biomass, length, %C, %N, delta13C, and delta15N. Only live fine roots were analyzed. Fine roots were not separated by species or plant functional type.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), 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↗

Tunneling away the relic neutrino asymmetry

The Earth acts as a matter potential for relic neutrinos which modifies their index of refraction from vacuum by δ ∼ 10 − 8 . It has been argued that the refractive effects from this potential should lead to a large O ( δ ) neutrino-antineutrino asymmetry at the surface of the Earth. This result was computed by treating the Earth as flat. In this work, we revisit this calculation in the context of a perfectly spherical Earth. We demonstrate, both numerically and through analytic arguments, that the flat-Earth result is only recovered under the condition δ 3 / 2 k R ≫ 1 , where k is the typical momentum of the relic neutrinos and R is the radius of the Earth. This condition is required to prevent antineutrinos from tunneling into classically inaccessible trajectories below the Earth’s surface and washing away the large asymmetry. As the physical parameters of the Earth do not satisfy this condition, we find that the asymmetry at the surface should only be O ( δ ) . While the asphericity of the Earth may serve as a loophole to our conclusions, we argue that it is still difficult to generate a large asymmetry even in the presence of local terrain. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Urban-Scale Processes in High-Spatial-Resolution Earth System Models

Cities and their associated urban areas have a much larger impact on the environment than their spatial footprint on the Earth's surface. This is undoubtedly true at local and regional scales, but also at global scales. At this time, more than half of the world's population lives in urban areas (80% in the United States), and the global proportion is projected to climb to 70% by 2050. Therefore, numerical models to study physics, chemistry, and biology affecting the Earth system at regional and global scales must represent the effects of urban areas on climate and the effects of a changing climate on urban areas. At the same time, it is essential to develop state of the art, simple and accurate urban models to better understand the relevant processes and also to address issues related to urban security against the spectra of chemical, biological and radiological (CBR) hazards. Towards bringing urban scientific communities together, a two-and-a-half-day international workshop was held at the Argonne National Laboratory (ANL) in the Chicago area on May 22-24, 2019. This workshop brought together more than 60 national and international experts to develop a roadmap for a better understanding of the issues associated with urban areas and at enhancing the capabilities of regional and global Earth System models (ESMs) in representing the atmospheric dynamics and chemistry, unique aspects of the biosphere and land use, and human dimensions of the urban environment.

54 ENVIRONMENTAL SCIENCES↗

Sulfur isotopes from the Paleoproterozoic Francevillian Basin record multigenerational pyrite formation, not depositional conditions

Bulk-rock sulfur isotope data from pyrite in the ~2.1 billion-year sedimentary rocks of the Francevillian Basin, Gabon, have underpinned ideas about initial oxygenation of Earth’s surface environments and eukaryote evolution. In this paper, we show, using micro-scale analytical methods, that the bulk sulfur isotope record represents progressive diagenetic modification. Our findings indicate no significant change in microbial sulfur cycling processes and seawater sulfate composition throughout that initial phase of atmosphere-ocean oxygenation of Paleoproterozoic time. This offers an alternative view of Earth system evolution during the transition from an anoxic to an oxic state and highlights the need for a judicious reappraisal of conceptual models using sulfur isotope data as primary depositional signals linked to global-scale biogeochemical processes.

58 GEOSCIENCES↗

Resolving glacial isostatic adjustment (GIA) in response to modern and future ice loss at marine grounding lines in West Antarctica

Accurate glacial isostatic adjustment (GIA) modelling in the cryosphere is required for interpreting satellite, geophysical and geological records and for assessing the feedbacks of Earth deformation and sea-level change on marine ice-sheet grounding lines. GIA modelling in areas of active ice loss in West Antarctica is particularly challenging because the ice is underlain by laterally varying mantle viscosities that are up to several orders of magnitude lower than the global average, leading to a faster and more localised response of the solid Earth to ongoing and future ice-sheet retreat and necessitating GIA models that incorporate 3-D viscoelastic Earth structure. Improvements to GIA models allow for computation of the viscoelastic response of the Earth to surface ice loading at sub-kilometre resolution, and ice-sheet models and observational products now provide the inputs to GIA models at comparably unprecedented detail. However, the resolution required to accurately capture GIA in models remains poorly understood, and high-resolution calculations come at heavy computational expense. We adopt a 3-D GIA model with a range of Earth structure models based on recent seismic tomography and geodetic data to perform a comprehensive analysis of the influence of grid resolution on predictions of GIA in the Amundsen Sea Embayment (ASE) in West Antarctica. Through idealised sensitivity testing down to sub-kilometre resolution with spatially isolated ice loading changes, we find that a grid resolution of ~$\frac{1}{3}$ of the radius of the load or higher is required to accurately capture the elastic response of the Earth. However, when we consider more realistic, spatially coherent ice loss scenarios based on modern observational records and future ice-sheet model projections and adopt a viscoelastic Earth, we find that predicted deformation and sea-level change along the grounding line converge to within 5 % with grid resolutions of 7.5 km or higher, and to within 2% for grid resolutions of 3.75 km and higher, even when the input ice model is on a 1 km grid. Furthermore, we show that low mantle viscosities beneath the ASE lead to viscous deformation that contributes to the instrumental record on decadal timescales and equals or dominates over elastic effects by the end of the 21st century. Our findings suggest that for the range of resolutions of 1.9–15 km that we considered, the error due to adopting a coarser grid in this region is negligible compared to the effect of neglecting viscous effects and the uncertainty in the adopted mantle viscosity structure.

54 ENVIRONMENTAL SCIENCES↗

Organic layer thickness and carbon concentration in burned and unburned sites, Seward Peninsula, AK, 2022

Measurements associated with organic layer samples collected from naturally burned (1971, 2002, 2015, 2019) and unburned sites at the Kougarok Fire Complex, Seward Peninsula, AK, 2022. Here, a discontinuous permafrost underlies an arctic tundra ecosystem. Measurements include elemental carbon and nitrogen concentrations and stocks, organic layer thickness, and thaw depth. There are five files in *.csv format with one data file and four data description files including data dictionary, methods, terminology, and file-level metadata. 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↗

Experimental and Theoretical Studies of the Surface Oxidation Process of Rare–Earth Tritellurides

Recent studies have established Van der Waals (vdW) layered and 2D rare-earth tritellurides (RTe 3 ) as superconductors and near room-temperature charge density wave (CDW) materials. Their environmental stability raises natural concern owing to aging/stability effects observed in other tellurium-based layered crystals. Here, the results establish the stability and environmental aging characteristics of these RTe 3 systems involving a variety of metals such as La, Nd, Sm, Gd, Dy, and Ho. The atomic force microscopy (AFM) and scanning electron microscopy (SEM) results show that all the RTe 3 sheets oxidize to form thin TeO x layers that are primarily confined to the surface, edges, and grain boundaries. Time-resolved in situ Raman spectroscopy measurements are used to understand the kinetics of the oxidization process for different lanthanide metal cations and establish their relative stability/resilience to oxidization. Overall results indicate that the vdW layers show higher air stability as the 4f electron number decreases going from Ho to La, resulting in the most stable LaTe 3 compared to the least stable HoTe 3 . Comprehensive quantum mechanical simulations reveal that environmental degradation originates from a strong oxidizing reaction with O 2 molecules, while humidity (H 2 O) plays a negligible role unless Te vacancies are present. Moreover, the simulations explain the effects of 4f electrons on the work function and Te vacancies formation, which directly impact the aging characteristics of RTe 3 layers. Interestingly, optical and electrical measurements show that the CDW response is still observed in aged RTe 3 layers owing to the presence of underlying pristine/nonoxidized RTe 3 layers, except CDW transition temperatures increase due to the thickness effect. Overall results offer the first in-depth environmental aging studies on these materials, which can be applied to engineer and design their chemical stability, surface properties, and overall CDW characteristics.

36 MATERIALS SCIENCE↗

CondiDiag1.0: a flexible online diagnostic tool for conditional sampling and budget analysis in the E3SM atmosphere model (EAM)

Abstract. Numerical models used in weather and climate prediction take into account a comprehensive set of atmospheric processes (i.e., phenomena) such as the resolved and unresolved fluid dynamics, radiative transfer, cloud and aerosol life cycles, and mass or energy exchanges with the Earth's surface. In order to identify model deficiencies and improve predictive skills, it is important to obtain process-level understanding of the interactions between different processes. Conditional sampling and budget analysis are powerful tools for process-oriented model evaluation, but they often require tedious ad hoc coding and large amounts of instantaneous model output, resulting in inefficient use of human and computing resources. This paper presents an online diagnostic tool that addresses this challenge by monitoring model variables in a generic manner as they evolve within the time integration cycle. The tool is convenient to use. It allows users to select sampling conditions and specify monitored variables at run time. Both the evolving values of the model variables and their increments caused by different atmospheric processes can be monitored and archived. Online calculation of vertical integrals is also supported. Multiple sampling conditions can be monitored in a single simulation in combination with unconditional sampling. The paper explains in detail the design and implementation of the tool in the Energy Exascale Earth System Model (E3SM) version 1. The usage is demonstrated through three examples: a global budget analysis of dust aerosol mass concentration, a composite analysis of sea salt emission and its dependency on surface wind speed, and a conditionally sampled relative humidity budget. The tool is expected to be easily portable to closely related atmospheric models that use the same or similar data structures and time integration methods.

58 GEOSCIENCES↗

Mineral surface-catalyzed oxidation of Mn(II) by bromate: Implications for the occurrence of Mn oxides on Mars

The occurrence of manganese (Mn) oxides on Mars is believed to be an indicator of an O 2 -rich paleoenvironment of Mars because Mn oxides often form through the oxidation of Mn(II) by O 2 on the surface of Earth. An alternative formation pathway was recently proposed, in which Mn(II) is oxidized by bromate (BrO 3 - ), a common oxidant in contemporary Martian regolith. However, the oxidation of Mn(II) by bromate in solution is kinetically controlled and slow unless using very high concentrations (100 mM) of reactants that may be irrelevant to the conditions of Mars. We conducted laboratory simulations to determine whether iron (Fe) oxides (hematite and goethite) and a phyllosilicate (montmorillonite), abundant minerals on the surface of Mars, could catalyze the oxidation of Mn(II) by bromate. Hematite and goethite, but not montmorillonite, dramatically accelerated the oxidation with a low concentration (1 mM) of Mn(II) and bromate under various solution conditions. Further, the reaction system was autocatalytic with Fe oxides initiating the oxidation of Mn(II) at the early stage and the subsequent catalysis mainly provided by the Mn oxide products. In contrast to producing Mn(IV)O 2 only during the homogeneous oxidation of Mn(II) by bromate in solutions, the heterogeneous mineral-surface catalyzed oxidation resulted in a mixture of Mn(III)OOH and Mn(IV)O 2 phases. Mn(III)OOH was an intermediate product and can be further oxidized by bromate to Mn(IV)O 2 . The occurrence and accumulation of the intermediate product MnOOH can be attributed to its rapid formation due to surface-enhanced nucleation and growth on Fe oxide surfaces and to its higher resistance to oxidation by bromate than Mn(III) ions or clusters. Overall, mineralsurface catalyzed oxidation of Mn(II) by bromate is favorable from both thermodynamic and kinetic perspectives, and can be a major pathway for the occurrence of Mn oxides on Mars where microorganisms are lacking to catalyze the reaction. Our study further improves our understanding of the thermodynamic and kinetic controls on Mn(II) oxidation.

54 ENVIRONMENTAL SCIENCES↗

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↗

DeFault: DEep‐Learning‐Based FAULT Delineation Using the IBDP Passive Seismic Data at the Decatur CO2 Storage Site

Abstract The carbon capture, utilization, and storage (CCUS) framework is an essential component in reducing greenhouse gas emissions, with its success hinging on the comprehensive knowledge of subsurface geology and geomechanics. Passive seismic event relocation and fault detection offer vital insights into subsurface structures and the ability to monitor fluid migration pathways. Accurate identification and localization of seismic events, however, face significant challenges, including the necessity for high‐quality seismic data and advanced computational methods. To address these challenges, we introduce a novel deep learning method, , specifically designed for passive seismic source relocation and fault delineating for passive seismic monitoring projects. By leveraging data domain‐adaptation, allows us to train a neural network with labeled synthetic data and apply it directly to field data. Using , the passive seismic sources are automatically clustered based on their recording time and spatial locations, and subsequently, faults and fractures are delineated accordingly. We demonstrate the efficacy of on a field case study involving injection related microseismic data from Decatur, Illinois area. Our approach accurately and efficiently relocated passive seismic events, identified faults and could aid in potential damage induced by seismicity. Our results highlight the potential of as a valuable tool for passive seismic monitoring, emphasizing its role in ensuring CCUS project safety. This research bolsters the understanding of subsurface characterization in CCUS, illustrating machine learning’s capacity to refine these methods. Ultimately, our work has significant implications for CCUS technology deployment, an essential strategy in combating climate change. Plain Language Summary In our quest to tackle climate change, we use a strategy known as carbon capture, utilization, and storage (CCUS) to keep greenhouse gases out of the atmosphere. This strategy relies heavily on our ability to understand what's happening deep under the earth's surface. To make sure we store super critical safely, we need to accurately map out the geological structure, especially faults, but this is tough without high‐quality data and complex computer programs. We've developed a new tool called “DeFault,” which uses advanced machine learning to improve how we find and map these underground features. “DeFault” is smart enough to learn from numerically simulated data and then apply what it’s learned to real‐world situations. It groups together seismic activity—tiny tremors and shifts in the earth—based on when and where they happen, which helps us spot where there might be cracks or faults. We tested “DeFault” in Illinois, where CO 2 is injected underground, and it successfully pinpointed where these tremors occurred and mapped out the faults, helping to prevent accidents accurately in the future. Our study shows that “DeFault” will be a powerful ally in making CCUS safer and more effective, especially for the Illinois Basin Decatur Project. Key Points Faults and fractures introduced by carbon storage can be monitored by passive seismicity DeFault algorithm enables an automatic process for accurate and efficient passive seismic event locating and clustering

58 GEOSCIENCES↗

Robust historical evapotranspiration trends across climate regimes

Evapotranspiration (ET) links the hydrological, energy and carbon cycles on the land surface. Quantifying ET and its spatio-temporal changes is also key to understanding climate extremes such as droughts, heatwaves and flooding. Regional ET estimates require reliable observation-based gridded ET datasets, and while many have been developed using physically based, empirically based and hybrid techniques, their efficacy, and particularly the efficacy of their uncertainty estimates, is difficult to verify. In this work, we extend the methodology used in Hobeichi et al. (2018) to derive two new versions of the Derived Optimal Linear Combination Evapotranspiration (DOLCE) product, with observationally constrained spatio-temporally varying uncertainty estimates, higher spatial resolution, more constituent products and extended temporal coverage (1980–2018). After demonstrating the efficacy of these uncertainty estimates with out-of-sample testing, we derive novel ET climatology clusters for the land surface, based on the magnitude and variability of ET at each location on land. The new clusters include three wet and three dry regimes and provide an approximation of Köppen–Geiger climate classes. The verified uncertainty estimates and extended time period then allow us to examine the robustness of historical trends spatially and in each of these six ET climatology clusters. We find that despite robust decreasing ET trends in some regions these do not correlate with behavioural ET clusters. Each cluster, and the majority of the Earth's surface, shows clear robust increases in ET over the recent historical period. The new datasets DOLCE V2.1 and DOLCE V3 can be used for benchmarking global ET estimates and for examining ET trends respectively.

54 ENVIRONMENTAL SCIENCES↗

rabpro: global watershed boundaries, river elevation profiles, and catchment statistics

River and Basin Profiler (rabpro) is a Python package to delineate watersheds, extract river flowlines and elevation profiles, and compute watershed statistics for any location on the Earth’s surface. As fundamental hydrologically-relevant units of surface area, watersheds are areas of land that drain via aboveground pathways to the same location, or outlet. Delineations of watershed boundaries are typically performed on digital elevation models (DEMs) that represent surface elevations as gridded rasters. Depending on the resolution of the DEM and the size of the watershed, delineation may be very computationally expensive. With this in mind, we designed rabpro to provide user-friendly workflows to manage the complexity and computational expense of watershed calculations given an arbitrary coordinate pair. In addition to basic watershed delineation, rabpro will extract the elevation profile for a watershed’s mainchannel flowline. This enables the computation of river slope, which is a critical parameter in many hydrologic and geomorphologic models. Finally, rabpro provides a user-friendly wrapper around Google Earth Engine’s (GEE) Python API to enable cloud-computing of zonal watershed statistics and/or time-varying forcing data from hundreds of available datasets. Altogether, rabpro provides the ability to automate or semi-automate complex watershed analysis workflows across broad spatial extents.

54 ENVIRONMENTAL SCIENCES↗

C fluxes from incubations of California upland soil with and without synthetic iron (hydr)oxide associations

There is a growing understanding that interactions with soil mineral phases contribute to the accumulation and retention of otherwise degradable organic matter in soils and sediments. However, the bioavailability of organic compounds in mineral-organic-associations, especially under varying environmental conditions is not well known. To assess the impact of mineral association and warming on the decomposition of an easily respirable organic substrate (glucose), we conducted a series of laboratory incubations at different temperatures with field-collected soils from 10 to 20 cm, 50–60 cm, and 80–90 cm depth at the Blodgett Forest Research Station. We added 13C-labeled glucose either directly to native soil or sorbed to one of two synthetic iron (hydr)oxide phases (goethite and ferrihydrite) that differ in crystallinity and affinity for sorbing glucose. We incubated the soils at 25 deg C. This dataset comprises the added glucose C and the cumulative fluxes over 72 and 1896 hours. We found that: (1) association with the Fe (hydr)oxide minerals reduced the decomposition rate of glucose by > 99.5% relative to rate of decomposition for free glucose in soil; (2) the respiration rate per gram carbon did not differ appreciably with depth, suggesting a similar degree of decomposability for native C across depths and that under the incubation conditions total carbon availability represents the principal limitation on respiration under these conditions as opposed to reduced abundance of decomposers or moisture and oxygen limitations; (3) addition of free glucose enhanced native carbon respiration at all soil depths with the largest effect at 50–60 cm; (4) in general respiration of the organo-mineral complex (glucose and iron-(hydr)oxide) was less temperature sensitive than was respiration of native carbon; (5) the addition of organic free mineral decreased the rate of soil respiration in the intermediate 50–60 cm depth soil. The results emphasize the key role of MOAs in regulating the fluxes of carbon from soils to the atmosphere and in turn the stocks of soil carbon.This dataset lists the glucose additions made to soils with and without added geothite or ferrihidrite and the C respired over the timespan of the incubations.

54 ENVIRONMENTAL SCIENCES↗

NGEE Arctic Soil Pit and Core Inventory for Samples Collected Between 2016-2019: Locations, Sampling Characteristics, Collection, and Contacts, Seward Peninsula, Alaska, USA

This dataset includes associated information about individual pit and core samples (e.g., location, contact, date, data collected) sampled between 2016 and 2019 from NGEE Arctic sites at Teller (MM 27 and MM47), Kougarok (MM47 and MM56, Hillslope, and vegetation plot areas), and Council (MM71). The inventory is qualitative and excludes those datasets on soil pit cores that included time-series information. The package includes one each *.csv data file, *.kml and *.pdf. All quantitative data is available in the related datasets https://doi.org/10.5440/1417652, https://doi.org/10.5440/1423892, https://doi.org/10.5440/1342956, https://doi.org/10.5440/1854940, https://doi.org/10.5440/1544760, https://doi.org/10.5440/1346200, and https://doi.org/10.5440/1856042. The goal of this dataset is to help modelers understand when, where, and what data was collected with respect to physical soil properties on the Seward Peninsula, AK, USA over the course of NGEE Arctic’s campaigns in the region from 2016-2019. 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↗

Reconciling evidence of oxidative weathering and atmospheric anoxia on Archean Earth

Evidence continues to emerge for the production and low-level accumulation of molecular oxygen (O 2 ) at Earth’s surface before the Great Oxidation Event. Quantifying this early O 2 has proven difficult. Here, we use the distribution and isotopic composition of molybdenum in the ancient sedimentary record to quantify Archean Mo cycling, which allows us to calculate lower limits for atmospheric O 2 partial pressures ( P O 2 ) and O 2 production fluxes during the Archean. We consider two end-member scenarios. First, if O 2 was evenly distributed throughout the atmosphere, then P O 2 > 10 –6.9 present atmospheric level was required for large periods of time during the Archean eon. Alternatively, if O 2 accumulation was instead spatially restricted (e.g., occurring only near the sites of O 2 production), then O 2 production fluxes >0.01 Tmol O 2 /year were required. Archean O 2 levels were vanishingly low according to our calculations but substantially above those predicted for an abiotic Earth system.

58 GEOSCIENCES↗

Surface heat flux and its association with the MJO in the tropical western Pacific using ARM observations

The Madden-Julian oscillation (MJO) is a major atmospheric phenomenon in the tropics that moves eastward every 20 to 100 days. It brings heavy rain and strong winds, influencing extreme weather events far beyond the tropics – including flooding, hurricanes, tornadoes and heavy snow in the United States. Therefore, an improved understanding of the MJO is critical for enhancing weather forecasts and supporting better decision-making for communities, emergency managers, and the private sector. Past studies, based on short-term observations or long-term model data, have emphasized the dominant role of atmospheric humidity in driving the MJO. Our research, using long-term observations (2000-2014) from three U.S. Department of Energy’s Atmospheric Radiation Measurement (ARM) sites in the tropical western Pacific (Manus, Nauru, and Darwin) confirms those past findings. However, we also show that atmospheric temperature, especially in the mid and upper troposphere, also plays a key role in the MJO’s evolution from its quiet (suppressed) to convective (active) phase. These findings provide valuable insights into how the MJO evolves and may help evaluate and improve weather and climate models. Our study also examined how the exchange of heat between the Earth's surface and the atmosphere, called surface heat flux, interacts with the MJO. This flux includes components such as latent heat (related to evaporation), sensible heat (related to temperature difference between surface and atmosphere), and radiation. We found that modulation in MJO convection is well connected to surface heat flux over the tropical western Pacific. One surface heat flux component often overlooked is the sensible heat flux due to precipitation (Q P ). Since falling raindrops are typically cooler than the land or ocean surface, they can cool the surface and affect local weather. To test the impact of Q P on convection during MJO, we incorporated it into a weather model, and ran simulations for two periods: April 2009 (when the MJO was active) and June 2006 (when it was not). Including Q P reduced errors in simulating the daily peak of rainfall—by 83% during the strong MJO and 23% during the inactive phase. It also improved the timing of peak daily rainfall and reduced the overall precipitation error by about 10%. These improvements were especially noticeable during periods of heavy rainfall. Our results suggest that even small heat exchanges from rainfall can play a significant role in shaping local weather. Accounting for these effects, particularly in regions like the tropical islands, can lead to more accurate simulations of the MJO-associated precipitation.

54 ENVIRONMENTAL SCIENCES↗