Search NASA⌕ Search

Engineering topics

Maher, Kate

Publications and source records attributed to Maher, Kate.

Quantifying Groundwater Response and Uncertainty in Beaver‐Influenced Mountainous Floodplains Using Machine Learning‐Based Model Calibration

Abstract Beavers ( Castor canadensis ) alter river corridor hydrology by creating ponds and inundating floodplains, and thereby improving surface water storage. However, the impact of inundation on groundwater, particularly in mountainous alluvial floodplains with permeable gravel/cobble layers overlain by a soil layer, remains uncertain. Numerical modeling across various floodplain structures considers topographic and sediment complexity and multidirectional flow, linking inundation to groundwater response. This study develops a model‐data integration workflow to address uncertainty in groundwater response to beaver‐induced inundations in a mountainous alluvial floodplain in the Upper Colorado River Basin. Uncertain factors include seasonal hydrologic dynamics, hydraulic conductivities, floodplain structures, and meteorological forcings. We employed an ensemble of groundwater models, based on geophysical and hydrologic data, with machine learning‐based calibration using a neural density estimator. This allowed us to quantify the vertical flux from the soil layer to the permeable gravel bed, the down‐valley underflow within the gravel bed, and their ratios. Results show a significant increase in the vertical flux relative to down‐valley underflow, from 2 during dry pond periods to 20 during wet periods, serving as an analogy for conditions without and with beaver ponds. The study highlights the influence of floodplain structure on groundwater storage, water balance, and water quality impacted by beaver ponds. A thick gravel bed layer, with a large down‐valley underflow, minimizes the effect of beaver‐induced inundation on water quality. We emphasize the need for field‐scale measurements of floodplain structure and improved characterization of evapotranspiration changes to reduce uncertainty in groundwater response. Plain Language Summary Beavers change the flow of water in river corridors by creating ponds, expanding wetlands, and flooding floodplains. This increases surface water area, promotes plant growth, and enhances biodiversity. However, the impact of this flooding on groundwater flow is not well understood, especially in mountainous areas with gravel layers where water moves easily beneath soil. In this study, we used numerical modeling to investigate how beaver ponds influence groundwater in a mountainous floodplain of the Upper Colorado River Basin. We adapted a machine learning method to validate our numerical models using multiple field data sets. Our findings show that beaver ponds significantly increase vertical water flow from the soil to the gravel during wet periods, compared to when the ponds are fully drained. The study also highlights the importance of floodplain structure in controlling both water flow in gravel layers along the river direction and vertical flow from the soil to the gravel with the presence of beavers. To reduce uncertainty in groundwater response, we emphasize the need for more field‐scale measurements of floodplain structure, hydraulic properties, and evapotranspiration changes. Key Points Floodplain structures and hydraulic conductivities are important for groundwater response with beaver ponds in mountainous floodplains Large down‐valley underflow in permeability‐stratified floodplains reduces beaver‐induced impacts on groundwater storage and water quality Machine learning‐based model calibration methods are effective for estimating posterior distributions of groundwater model parameters

Wang, Lijing↗

Mountainous Floodplain Connectivity in Response to Hydrological Transitions

In mountainous watersheds, floodplain sediments are typically characterized by gravel bed layers capped by an overlying soil unit that serves as a hotspot for biogeochemical reactivity. However, the influence of soil biogeochemistry on gravel bed underflow composition remains unclear, especially during hydrological transitions that alter the vertical connectivity between overlaying soils and the underlying gravel bed. This study investigates these dynamics by measuring hydraulic gradients and water compositions over three hydrological years in a typical mountainous, low-order stream floodplain in the Upper Colorado River Basin. Results indicate that the timing of hydrological conditions strongly influences the vertical exchanges that control water quality. Specifically, during flooding events such as beaver ponding, that induce downward flushing of the soil, anoxic conditions prevalent in the biogeochemically active soil are transferred downstream via gravel bed underflow. Conversely, snowmelt and drought conditions increase oxic conditions in the gravel bed due to diminished hydrological connectivity with the overlying soil. To compare water quality response to hydrological transitions across similar floodplain environments, we propose a conceptual model that quantifies the inundation-induced flushing of soil porewater to measure solute exchange efficiency with the gravel bed solute convergence efficiency (SCE). This model provides a framework for quantifying biogeochemical processes in hydrological underflow systems, which is critical for water and elemental budgets in these globally important mountainous ecosystems.

58 GEOSCIENCES↗

NEON AOP Survey of Upper East River CO Watersheds: Waveform LiDAR Binary Data

The waveform Light Detection and Ranging (LiDAR) data in this package were generated through a National Ecological Observatory Network Airborne Observation Platform (NEON AOP) acquisition over watersheds of interest surrounding Crested Butte, Colorado. The remote sensing imagery acquired by the NEON AOP supports an interdisciplinary project on hydrology, biogeochemistry, and ecosystem functioning in a snow-dominated headwater environment. These waveform LiDAR data enable spatially continuous estimation of vegetation structure parameters to facilitate analyses of the major environmental drivers of structural and compositional variability. The package contains 97 compressed file archives in 7-zip (.7z) format, each corresponding to one acquisition flightpath. Within each .7z archive is a set of constituent files describing properties of the LiDAR waveforms, such as return intensity, geolocation, outgoing pulse and other behavior of the sensor and signals. Once downloaded, the files must first be unzipped using the widely distributed command-line software utility 7z, using the command '7z x \[filename\].7z \[target_directory\]'. All files within the .7z archives can be opened in IDL, MatLab, or the open-source R statistical computing environment. Further details about the data package are in the attached user guide (neon_aop_crbu_waveformlidar_userguide.pdf). This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

2018 NEON and 2025 CHESS Campaigns↗

Soil water content, matric potential, carbon dioxide and oxygen concentrations, Oct 2018-Dec 2021, Slate River Floodplain, Crested Butte, Colorado

This data package includes a time series of soil sensor data (temperature, water content, bulk electrical conductivity, porewater dissolved oxygen and porewater dissolved carbon dioxide) in a vertical profile from the Slate River floodplain outside Crested Butte, Colorado, a focus field site for the SLAC Floodplain Hydro-Biogeochemistry SFA. The data was generated as part of the work targeting the overarching research question for the SLAC SFA: How do ubiquitous subsurface interfaces mediate molecular-scale biogeochemical processes and groundwater quality in floodplains and watersheds? The package includes: (1) soil temperature, volumetric water content and electrical conductivity at 40, 60 and 82.5 cm depth; (2) soil matric potential at 40, 60, 79 and 100 cm depth; (3) soil CO2 concentrations at 40, 60 and 82.5 cm depth; and (4) soil oxygen concentrations at 60, 82.5, 100, 135, 170 and 182 cm depth. Both the carbon dioxide and oxygen sensors are optical sensors that can measure the partial pressure of oxygen in both saturated and unsaturated conditions. Unfortunately, soil CO2 in the profile is unexpectedly high and above the sensor calibration range (0-25,000 ppm). In addition, soil CO2 sensors failed within a year of deployment, so we only report CO2 data from 2019-2020.Within the data package, "FLMD.csv" describes file-level metadata and "dd.csv" defines column headers and universal terms across the dataset. The data package includes 4 "*data.csv" files, one for each calendar year in the dataset. Each "*data.csv" file has a corresponding "*_InstallationMethods.csv" file that describes the location, sensor model, sensor serial number and other metadata corresponding for each measured parameter. Because sensors have been added over time, not every sensor has data dating back to Oct 2018. Note that there is a data gap over winter 2019-2020 due to a power outage. While this repository currently only contains data through December 2021, the dataset will be updated as additional years are collected and processed.

54 ENVIRONMENTAL SCIENCES↗

Meteorological, reference evapotranspiration and estimated transpiration data, July 2020-Dec 2021, Slate River Floodplain, Crested Butte, Colorado

This data package includes a time series of meteorological data (air temperature, relative humidity, barometric pressure, precipitation, wind speed, soil heat flux, incoming shortwave radiation, incoming longwave radiation, outgoing shortwave radiation and outgoing longwave radiation) from the Slate River floodplain outside Crested Butte, Colorado, a focus field site for the SLAC Floodplain Hydro-Biogeochemistry SFA. The data was generated as part of the work targeting the overarching research question for the SLAC SFA: How do ubiquitous subsurface interfaces mediate molecular-scale biogeochemical processes and groundwater quality in floodplains and watersheds? The package also includes reference evapotranspiration calculated via the Penman-Monteith equation and an estimate of willow (Salix spp.) evapotranspiration using the method of Irmak et al., J. Hydrology, 2013. Due to two malfunctioning sensors (incoming shortwave radiation from Oct 2020-May 2021 and incoming longwave radiation from May 2021-March 2022), observed net radiation and calculated evapotranspiration are only included from July-October 2020. However, users may be able to extend calculated evapotranspiration using estimates of these two parameters from other sources. Both sensors have been replaced and data beyond 2021 be uploaded to this repository when it is available.Within the data package, "FLMD.csv" describes file-level metadata and "dd.csv" defines column headers and universal terms across the dataset. The data package includes 3 "*data.csv" files, two for meteorological data for 2020 and 2021 and one for reference evapotranspiration and transpiration data from 2020 ("SLT_OBJ2_PenmanMonteith_2020_data.csv"). Each "*data.csv" file of meteorological data has a corresponding "*_InstallationMethods.csv" file that describes the location, sensor model, sensor serial number and other metadata corresponding to each measured parameter. We also include a Jupyter notebook ("PenmanMonteithET_Calculations.ipynb") that details the assumptions and equations used to calculate daily and hourly reference evapotranspiration and daily willow transpiration from raw meteorological data.

54 ENVIRONMENTAL SCIENCES↗

The circular nutrient economy of terrestrial ecosystems and the consequences for rock weathering

Earth’s biosphere is thought to exert a substantial influence on regolith evolution and chemical weathering rates. However, ecosystems are also highly efficient at retaining and recycling nutrients. Thus, when the ecological demand for rock-derived nutrients (e.g., P, Ca, K) exceeds the rates of regolith supply, ecological retention and recycling strategies can minimize nutrient limitations. To evaluate the balance between nutrient recycling and new nutrient input, we combined a plant model that drives growth according to foliar P levels with a weathering model that includes regolith rejuvenation via erosion and export via chemical weathering according to water flow, regolith thickness, mineral dissolution rates, secondary minerals, and nutrient storage in organic and mineral phases. We find that plant growth is strongly dependent on the total regolith nutrient inventory, resulting in a strong correlation between plant productivity and erosion. Increased water export or decreased regolith thickness diminish the total inventory of nutrient corresponding to lower rates of recycling and lower plant growth. In contrast, purported biogenic drivers of weathering, such as enhanced mineral dissolution, only support higher growth rates at high erosion rates. At erosion rates typical of the global land surface, more rapid mineral dissolution combined with enhanced formation of secondary minerals, depletes the inventory of mineral P, resulting in no benefit for plant growth. We also find that the increased chemical weathering export does not scale directly with plant growth. For example, accelerated mineral weathering does increase chemical weathering export but not potential plant growth. Conversely, thicker regolith is associated with a small increase in weathering export, but a large increase in potential plant growth. Collectively, when plant growth is coupled to regolith weathering our calculations suggest that plant productivity is not directly correlated with silicate weathering fluxes, and that biotic drivers of silicate weathering may only be effective at high erosion rates not typical at the Earth’s surface.

58 GEOSCIENCES↗

Simulation of anoxic lenses as exporters of reactivity in alluvial aquifer sediments

Sedimentary interfaces between contrasting hydrogeological facies in alluvial aquifers drive the development of biogeochemical interfaces that influence subsurface and surface water quality. Here, we calibrate a reactive transport model on a series of dual-domain column experiments, where centimeter-scale, low-permeability, organic-rich anoxic lenses are embedded in coarser-grained aquifer material. Simulations explicitly account for C, Fe, and S cycling at the interface between the lenses and the aquifer in water-saturated conditions. Our results highlight the role of fine-grained, organic-rich inclusions not only as sources or sinks for redox-sensitive species, but also as exporters of nutrients that stimulate downgradient biogeochemical cycling. By releasing large amounts of organic carbon into the surrounding aquifer, such lenses drive the development of proximal secondary reduction zones ("halos"), characterized by high microbial activity (e.g., sulfate reduction) and accumulation of reduced reaction products (e.g., iron sulfide). If similar secondary reduction zones develop within the hydraulically conductive domain of an aquifer, they would be highly susceptible to changes in hydrologic conditions, for instance oxygen pulses associated with seasonal snowmelt. Furthermore, our results also emphasize the limitations of relying solely on aqueous species measurements to inform reactivity in systems where fast redox cycling and/or sizeable particulate transport may limit the signature of reactivity in the dissolved phase.

54 ENVIRONMENTAL SCIENCES↗

Effective kinetics driven by dynamic concentration gradients under coupled transport and reaction

Biogeochemical reaction kinetics are generally established from batch reactors where concentrations are uniform. In natural systems, many biogeochemical processes are characterized by spatially and temporally variable concentration gradients that often occur at scales which are not resolved by field measurements or biogeochemical and reactive transport models. Yet, it is not clear how these sub-scale chemical gradients affect reaction kinetics compared to batch kinetics. Here we investigate this question by studying the paradigmatic case of localized pulses of solute reacting with a solid or a dissolved species in excess. Additionally, we consider non-linear biogeochemical reactions, representative of mineral dissolution, adsorption and redox reactions, which we quantify using simplified power-law kinetics. The combined effect of diffusion and reaction leads to effective kinetics that differ quantitatively and qualitatively from the batch kinetics. Depending on the nonlinearity (reaction order) of the local kinetics, these effects lead to either enhancement or decrease of the overall reaction rate, and result in a rich variety of reaction dynamics. We derive analytical results for the effective kinetics, which are validated by comparison to direct numerical simulations for a broad range of Damköhler numbers and reaction order. Our findings provide new insights into the interpretation of imperfectly mixed lab experiments, the effective kinetics of field systems characterized by intermittent reactant release and the integration of sub-scale concentration gradients in reactive transport models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Global Sensitivity Analysis of a Reactive Transport Model for Mineral Scale Formation During Hydraulic Fracturing

Injection of water-based hydraulic fracturing fluid (HFF) into tight shale gas/oil formations can increase formation permeability and enhance production rates, but this process frequently causes mineral scale formation that can occlude pore space and hinder flow. To identify the most important factors that control the formation of mineral scales, we applied a novel global sensitivity analysis method—distance-based generalized sensitivity analysis (DGSA)—to a reactive transport model (RTM) that was previously built and calibrated to simulate precipitation of barite [BaSO4] and iron (hydr)oxide [Fe(OH) 3 ] in shale matrices and on fracture surfaces. Reactive transport simulations were run with model parameters randomly sampled based on assigned uncertainties. Modeling results for barite and Fe(OH)3 formation were clustered using machine-learning algorithms. A list of ranked critical input parameters was obtained after statistical quantification of cumulative distribution functions of input parameters. We found that barite formation is most sensitive to the rate of sulfate ion generation, which is determined by the pyrite dissolution rate coefficient and oxidant availability. In addition, barite formation is sensitive to the initial amounts of barite in HFF and shale, followed by barite thermodynamics/kinetics. For Fe(OH) 3 formation, the ranked factors are Fe(OH)3 precipitation rate coefficients, initial HFF pH, initial Fe(OH) 3 amount in HFF, and oxidant availability. Overall, our results provide insights into managing mineral scale formation during hydraulic fracturing to enhance production. Meanwhile, this study serves as an example of global sensitivity analysis of RTMs using the efficient, straightforward, and open-source DGSA method.

58 GEOSCIENCES↗

Stochastic simulations of temporally and spatially variable Fe cycling within a floodplain aquifer

This data package contains model input and output files for local and global sensitivity analysis (SA) of a reactive transport model simulating redox cycling within a floodplain aquifer. The reaction network is implemented in CrunchFlow and focuses on spatially and temporally heterogeneous Fe cycling, though it contains several other species and redox pathways (37 aqueous species and 7 minerals in total). The package contains data from 600 Monte Carlo simulations (for global SA) and 104 one-at-a-time simulations (for local SA). It includes all input files required to run the CrunchFlow model, as well as limited model output – .csv files containing time series of primary species concentrations at 21 grid cells within the domain. These time series files contain the following species: H+, CO32-, SO42-, Cl-, Ca2+, Mg2+, Na+, K+, Fe(II), Fe(III), HS-, CH3COO-, O2(aq), S(aq), and Br-.These data were used to demonstrate the application of a novel form of global SA, distance-based generalized sensitivity analysis (DGSA) to a reactive transport model. The 2D model simulates the export of reduced species from pockets of fine-grained, organic-rich sediments embedded in a coarse sand aquifer. Stochastic simulations jointly varied 17 key model input parameters across several orders of magnitude, including the spatial variability of each parameter. DGSA on the model results reveal that the amplitude and variability of exported Fe(II) are most sensitive to the interaction between sand permeability and individual reaction rates. By contrast, the propagation of reducing conditions downgradient from the fine-grained lenses depends most heavily on the rate of dissolved organic carbon and sulfate release from the fines.

54 ENVIRONMENTAL SCIENCES↗

Soil Respiration Response to Rainfall Modulated by Plant Phenology in a Montane Meadow, East River, Colorado, USA

Soil respiration is a primary component of the terrestrial carbon cycle. However, predicting the response of soil respiration to climate change remains a challenge due to the complex interactions between environmental drivers, especially plant phenology, temperature, and soil moisture. In this study, we use a 1-D diffusion-reaction model to calculate depth-resolved CO 2 production rates from soil CO 2 concentrations and surface efflux observations in a subalpine meadow in the East River watershed, CO. Modeled rates are compared to in situ soil temperature and moisture conditions and MODIS satellite enhanced vegetation index (EVI) representing plant phenology across three hydrologically distinct growing seasons from 2016-2018. While soil respiration correlated with temperature on diel timescales (p<0.05), seasonal variability was dominated by soil moisture and plant phenology (p<0.05). We observed significant respiration increases in response to precipitation events; however, magnitude and duration were significantly higher in 2017 than 2016 despite similar wetting characteristics. Based on MODIS EVI, we suggest the respiration response to rainfall is controlled by plant phenology, which in turn reflects the capacity of plants to respond to precipitation via increased photosynthesis and autotrophic respiration, behavior that is not captured in typical soil respiration pulse models. Projected changes in montane climate such as earlier snowmelt and prolonged fore-summer drought may decrease soil respiration fluxes by decreasing the overlap between peak productivity and the summer monsoon. Lastly, we observed significant late-season CO 2 fluxes from the deep sub-soil (>165 cm) that support growing evidence for the importance of sub-soil processes in driving integrated respiration fluxes.

54 ENVIRONMENTAL SCIENCES↗

Chromium isotope fractionation during reduction of Chromium(VI) by Iron(II/III)-bearing clay minerals

Chromium stable isotope ratios are used to trace the reduction of Cr(VI) to Cr(III) in both ancient and modern systems. However, quantitative interpretation of Cr isotopic signatures has been stymied by the large variability in isotopic fractionation factors for Cr(VI) reduction by different reductants. Here we determine Cr isotope fractionation factors during Cr(VI) reduction by Fe(II/III)-bearing clay minerals, which are abundant in subsurface environments. Several variables were tested: pH, total Fe content of the clay, and the fraction of reduced Fe within the clay (Fe(II)/Fe(total)). The latter controls the standard reduction potential of the clay. Our results demonstrate that neither pH nor total Fe content of the clay have major effects on isotopic fractionation. In contrast, as the effective standard reduction potential of the clay and thus the standard free energy of Cr(VI) reduction become more negative, Cr isotope fractionation factors decrease in magnitude from - 4.9 to -1.3‰ according to a linear free energy relationship. This linear free energy relationship can be predicted from Marcus electron transfer theory and allows first-order predictions of Cr isotope fractionation factors to be made from the standard reduction potential or Fe(II)/Fe(total) of a clay, potentially improving our ability to model Cr isotope signatures in geochemical systems. Chromium is the first isotope system to show such a linear free energy relationship over a diverse range of reductants, including both aqueous and solid-phase reductants, and may provide a model for determining other redox-driven kinetic isotope effects in environmentally important isotope systems.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Stability of Floodplain Subsurface Microbial Communities Through Seasonal Hydrological and Geochemical Cycles

Riparian floodplains represent an interaction zone between the terrestrial subsurface and rivers, where regional groundwater flows, infiltration, and evapotranspiration drive mixing of water and import/export of nutrients and contaminants. These dynamics create seasonally transient redox conditions that drive biogeochemical transformations, which strongly modify groundwater quality. Microbial responses to changing hydrological conditions are perhaps the critical step connecting hydrology to geochemical transformations and groundwater quality, yet are not well understood. We aimed to address this knowledge gap by monitoring seasonal transitions at the U.S. Department of Energy legacy uranium ore processing site in Riverton, WY, through spring-summer-fall hydrological transitions. Our goal was to characterize the microbial community throughout the soil profile, down to the saturated aquifer, and observe its response to wet-dry transitions across a full season and compare to changes in geochemistry and hydrology. Next-generation sequencing was employed to identify biogeochemically-relevant microbial taxa based on the 16S rRNA gene; we found a broad diversity of microbial clades including taxa involved in sulfur and metal cycling, as well as nitrification. These data were paired with measurements of soil moisture, major nutrients and cations, and trace elements. Overall microbial community composition was dependent on soil depth or type, with seasonal effects only observed in the topsoil or subsurface aquifer. This finding indicates that microbial communities in the transiently-reduced center of the soil profile at the Riverton, WY site are remarkably stable, despite moisture and redox inversions. In addition, these communities likely impact the communities in surrounding soil horizons through export of metabolites and solutes as the water table rises and falls throughout the season.

16S rRNA↗

Persistence of soil organic carbon caused by functional complexity

Soil organic carbon management has the potential to aid climate change mitigation through drawdown of atmospheric carbon dioxide. To be effective, such management must account for processes influencing carbon storage and re-emission at different space and time scales. Achieving this requires a conceptual advance in our understanding to link carbon dynamics from the scales at which processes occur to the scales at which decisions are made. Here, we propose that soil carbon persistence can be understood through the lens of decomposers as a result of functional complexity derived from the interplay between spatial and temporal variation of molecular diversity and composition. For example, co-location alone can determine whether a molecule is decomposed, with rapid changes in moisture leading to transport of organic matter and constraining the fitness of the microbial community, while greater molecular diversity may increase the metabolic demand of, and thus potentially limit, decomposition. This conceptual shift accounts for emergent behaviour of the microbial community and would enable soil carbon changes to be predicted without invoking recalcitrant carbon forms that have not been observed experimentally. Finally, functional complexity as a driver of soil carbon persistence suggests soil management should be based on constant care rather than one-time action to lock away carbon in soils.

54 ENVIRONMENTAL SCIENCES↗

Depth-resolved soil data from DOE-LM Riverton Wyoming site, April-September 2017

Averaged (from 3 to 5 replicates) and QA/QC:d data for soil cores taken roughly every month in the period April 5 to September 13 in 2017 at a location close to DOE Legacy Management well 855 at the Riverton, Wyoming floodplain site. The groundwater at this site exhibits persistent U, Mo, and sulfate plumes and is one of the field sites in focus for the SLAC Groundwater Quality SFA program. Cores were separated into 5-20 cm segments based on soil horizonation and each segment was subsampled for microbial analyses (data available at the NCBI Single Read Archive (SRA) Database, https://www.ncbi.nlm.nih.gov/sra/, BioProject ID# PRJNA626616) in addition to the soil analyses included in this dataset. The soil analyses dataset includes aerobic and anaerobic respiration rate estimates (MicroResp method), water content estimates (on collected samples, not in situ), total C, N, S, Ca, K, Mg, Na, Fe, As, Mo, U, and organic-C concentrations, water extractable concentrations of the same elements, and HCl-extractable Fe(II) concentrations.There are four files included in this data package: metadataFile.csv, avgDataFile.csv, methodFile.csv, and RVT2017soilData.xlsx. The excel file contains the same information as the individual csv-files (in separate tabs).

54 ENVIRONMENTAL SCIENCES↗

Soil Environmental Conditions in the East River Watershed, CO

This dataset contains depth-resolved soil environmental variables collected between June 2016 and October 2018 from the Bradley Meadow site in the East River watershed. Datalogger output from 7 Campbell CS655 12 cm Soil Moisture and Temperature Sensors including depth, temperature, volumetric water content, and electrical conductivity are included in ‘BM_Merged_T_VWC_0616_10_18.csv’. Sensors were installed in June 2016. Data have been processed to merge raw sensor data and to remove additional variables. No additional QA/QC was performed. Measured pCO2 and depth are contained in ‘BM_EGM_WELL_CO2.csv’. These data represent gas well samples measured via injection using a PP Systems EGM-5 infra-red gas analyzer as described in the methods section. Measured CO2 surface flux values are contained in 'BM_EGM_Surface_Flux.csv'. These were measured using the PP Systems EGM-5 with the Soil Respiration Chamber attachment. Multiple reps are reported for each measurement including CO2 flux rate calculations using built-in linear and quadratic function methods, and average and standard deviations are taken across all reps and calculation methods.

54 ENVIRONMENTAL SCIENCES↗

Sample Identifiers and Metadata Reporting Format for Environmental Systems Science

The ESS-DIVE sample identifiers and metadata reporting format primarily follows the System for Earth Sample Registration (SESAR) Global Sample Number (IGSN) guide and template, with modifications to address Environmental Systems Science (ESS) sample needs and practicalities (IGSN-ESS). IGSNs are associated with standardized metadata to characterize a variety of different sample types (e.g. object type, material) and describe sample collection details (e.g. latitude, longitude, environmental context, date, collection method). Globally unique sample identifiers, particularly IGSNs, facilitate sample discovery, tracking, and reuse; they are especially useful when sample data is shared with collaborators, sent to different laboratories or user facilities for analyses, or distributed in different data files, datasets, and/or publications. To develop recommendations for multidisciplinary ecosystem and environmental sciences, we first conducted research on related sample standards and templates. We provide a comparison of existing sample reporting conventions, which includes mapping metadata elements across existing standards and Environment Ontology (ENVO) terms for sample object types and environmental materials. We worked with eight U.S. Department of Energy (DOE) funded projects, including those from Terrestrial Ecosystem Science and Subsurface Biogeochemical Research Scientific Focus Areas. Project scientists tested the process of registering samples for IGSNs and associated metadata in workflows for multidisciplinary ecosystem sciences. We provide modified IGSN metadata guidelines to account for needs of a variety of related biological and environmental samples. While generally following the IGSN core descriptive metadata schema, we provide recommendations for extending sample type terms, and connecting to related templates geared towards biodiversity (Darwin Core) and genomic (Minimum Information about any Sequence, MIxS) samples and specimens. ESS-DIVE recommends registering samples for IGSNs through SESAR, and we include instructions for registration using the IGSN-ESS guidelines. Our resulting sample reporting guidelines, template (IGSN-ESS), and identifier approach can be used by any researcher with sample data for ecosystem sciences.

54 ENVIRONMENTAL SCIENCES↗