Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “soil depth”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A continental scale analysis reveals widespread root bimodality

An improved understanding of root vertical distribution is crucial for assessing plant-soil-atmosphere interactions and their influence on the land carbon sink. Here, we analyze a continental-scale dataset of fine roots reaching 2 meters depth, spanning from Alaskan tundra to Puerto Rican forests. Contrary to the expectation that fine root abundance decays exponentially with depth, we found root bimodality at ~20% of 44 sites, with secondary biomass peaks often below 1m. Root bimodality was more likely in areas with low total fine root biomass and was more frequent in shrublands than grasslands. Notably, secondary peaks coincided with high soil nitrogen content at depth. Our analyses suggest that deep soil nutrients tend to be underexploited, while root bimodality offers plants a mechanism to tap into deep soil resources. Our findings add to the growing recognition that deep soil dynamics are systematically overlooked, and calls for more research attention to this deep frontier in the face of global environmental change.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Measured and Simulated Carbon Dynamics in Midwestern U.S. Corn‐Soybean Rotations

Abstract Corn ( Zea mays L.) and soybean ( Glycine max [L.] Merr.) production dominate Midwestern U.S. agriculture and impact the regional carbon and nitrogen cycles. Sustaining soil carbon is important for corn‐soybean production (CS); however, quantifying soil carbon changes requires long‐term field measurements and/or model simulations. In this study, changes in soil organic (SOC), inorganic (SIC), and total (TC) carbon; pH; total nitrogen (TN); and net ecosystem production (NEP) were measured in a conventional corn‐soybean rotation in central Iowa. Soil samples ( n = 42; 0–120 cm depth) were collected from two adjacent fields in 2005 and 2016. Eddy‐flux stations set up in the fields continuously monitored NEP from 2005–2016, and net biome production (NBP) was calculated using yield records. The DayCENT model was used to simulate the effects of conventional management practices on soil carbon and calibrated with field‐measured NEP and SOC. Measured soil TC (0–120 cm) decreased by −14.19 ± 6.25 Mg ha −1 , with highest reductions in SOC and SIC ( p < 0.05) at 0–15 and 90–120 cm, respectively. Measured TN decreased by −0.7 ± 0.29 Mg ha −1 with N‐accumulation at 60–90 cm ( p < 0.05). Eddy‐flux NBP decreased by −13.19 ± 0.05 Mg ha −1 . Soil and eddy‐flux records show a carbon reduction by −1.14 ± 0.63 and −1.20 ± 0.06 Mg ha −1 yr −1 , respectively. The validated DayCENT model suggests that all SOC pools declined. We postulate that conventional CS production has adverse effects on C and N dynamics in Midwestern United States.

Dold, C.↗
DOE OSTI2021

Effects of Slope Aspect and Rainfall on Belowground Deep Fine Root Traits and Aboveground Tree Height

The vertical root distribution and rooting depth are the main belowground plant functional traits used to indicate drought resistance in arid and semiarid regions. The effects of the slope aspect on the aboveground traits are visible but not the belowground deep root traits. We aimed to investigate the fine root traits of the locust tree ( Robinia pseudoacacia L.) planted on southerly and northerly aspects, and the variations in the rooting depth in regions with different rainfall, as well as assessing how deep rooting, might affect the response to drought in a loess region. We selected three study sites with different rainfall amounts, with six sampling plots at each site (three each with southerly and northerly aspects). Soil core samples were collected down to the depth where no roots were present. The locust trees tended to develop deeper fine roots rather than greater heights. The tree height and diameter were greater for locust trees on northerly aspects, whereas trees on southerly aspects had significantly deeper rooting depths. Fine root traits (root length, root area, and root dry weight density) were higher in the southerly aspect for both Changwu and Ansai, but lower in Suide. The ratio of the root front depth tree height ranged from 1.04 to 3.17, which was higher on southerly than northerly aspects, and it increased as the rainfall decreased. Locust tree growth traits (belowground fine root and aboveground tree height) were positively correlated with the mean annual rainfall. The soil moisture content of the topsoil decreased as the rainfall decreased, but the pattern varied in the deep layer. Our results suggest that the variations in the belowground rooting depth under different slope aspects may be related to plant survival strategies. The vertical extension of the rooting depth and tree height may be key functional traits that determine plant growth in drought-prone regions.

Lihui, Ma↗
DOE OSTI2021

Theoretical evidence that root penetration ability interacts with soil compaction regimes to affect nitrate capture

Abstract Background and Aims Although root penetration of strong soils has been intensively studied at the scale of individual root axes, interactions between soil physical properties and soil foraging by whole plants are less clear. Here we investigate how variation in the penetration ability of distinct root classes and bulk density profiles common to real-world soils interact to affect soil foraging strategies. Methods We utilize the functional–structural plant model ‘OpenSimRoot’ to simulate the growth of maize (Zea mays) root systems with variable penetration ability of axial and lateral roots in soils with (1) uniform bulk density, (2) plow pans and (3) increasing bulk density with depth. We also modify the availability and leaching of nitrate to uncover reciprocal interactions between these factors and the capture of mobile resources. Key Results Soils with plow pans and bulk density gradients affected overall size, distribution and carbon costs of the root system. Soils with high bulk density at depth impeded rooting depth and reduced leaching of nitrate, thereby improving the coincidence of nitrogen and root length. While increasing penetration ability of either axial or lateral root classes produced root systems of comparable net length, improved penetration of axial roots increased allocation of root length in deeper soil, thereby amplifying N acquisition and shoot biomass. Although enhanced penetration ability of both root classes was associated with greater root system carbon costs, the benefit to plant fitness from improved soil exploration and resource capture offset these. Conclusions While lateral roots comprise the bulk of root length, axial roots function as a scaffold determining the distribution of these laterals. In soils with high soil strength and leaching, root systems with enhanced penetration ability of axial roots have greater distribution of root length at depth, thereby improving capture of mobile resources.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

A Rapid Decision Support Tool for Estimating Impacts of a Vadose Zone Volatile Organic Compound Source on Groundwater and Soil Gas

Diminishing rates of subsurface volatile contaminant removal by soil vapor extraction (SVE) oftentimes warrants an in‐depth performance assessment to guide remedy decision‐making processes. Such a performance assessment must include quantitative approaches to better understand the impact of remaining vadose zone contamination on soil gas and groundwater concentrations. The spreadsheet‐based Soil Vapor Extraction Endstate Tool (SVEET) software functionality has recently been expanded to facilitate quantitative performance assessments. The updated version, referred to as SVEET2, includes expansion of the input parameter ranges for describing a site (site geometry, source characteristics, etc.), an expanded list of contaminants, and incorporation of elements of the Vapor Intrusion Estimation Tool for Unsaturated‐zone Sources software to provide soil gas concentration estimates for use in vapor intrusion evaluation. As part of the update, SVEET2 was used to estimate the impact of a tetrachloroethene (PCE) vadose zone source on groundwater concentrations, comparing SVEET2 results to field‐observed values at an undisclosed site where SVE was recently terminated. PCE concentrations from three separate monitoring wells were estimated by SVEET2 to be within the range of 6.0–6.7 μg/L, as compared to actual field concentrations that ranged from 3 to 11 μg/L PCE. These data demonstrate that SVEET2 can rapidly provide representative quantitative estimates of impacts from a vadose zone contaminant source at field sites. Finally, in the context of the SVE performance assessment, such quantitative estimates provide a basis to support remedial and/or regulatory decisions regarding the continued need for vadose zone volatile organic compound remediation or technical justification for SVE termination, which can significantly reduce the cost to complete for a site.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Unraveling the depth-dependent causal dynamics of methanogenesis and methanotrophy in a high-latitude fen peatland

The dynamics of methane (CH 4 ) cycling in high-latitude peatlands through different pathways of methanogenesis and methanotrophy are still poorly understood due to the spatiotemporal complexity of microbial activities and biogeochemical processes. Additionally, long-term in situ measurements within soil columns are limited and associated with large uncertainties in microbial substrates (e.g. dissolved organic carbon, acetate, hydrogen). To better understand CH 4 cycling dynamics, we first applied an advanced biogeochemical model, ecosys , to explicitly simulate methanogenesis, methanotrophy, and CH 4 transport in a high-latitude fen (within the Stordalen Mire, northern Sweden). Next, to explore the vertical heterogeneity in CH 4 cycling, we applied the PCMCI/PCMCI+ causal detection framework with a bootstrap aggregation method to the modeling results, characterizing causal relationships among regulating factors (e.g. temperature, microbial biomass, soil substrate concentrations) through acetoclastic methanogenesis, hydrogenotrophic methanogenesis, and methanotrophy, across three depth intervals (0–10 cm, 10–20 cm, 20–30 cm). Our results indicate that temperature, microbial biomass, and methanogenesis and methanotrophy substrates exhibit significant vertical variations within the soil column. Soil temperature demonstrates strong causal relationships with both biomass and substrate concentrations at the shallower depth (0–10 cm), while these causal relationships decrease significantly at the deeper depth within the two methanogenesis pathways. In contrast, soil substrate concentrations show significantly greater causal relationships with depth, suggesting the substantial influence of substrates on CH 4 cycling. CH 4 production is found to peak in August, while CH 4 oxidation peaks predominantly in October, showing a lag response between production and oxidation. Overall, this research provides important insights into the causal mechanisms modulating CH 4 cycling across different depths, which will improve carbon cycling predictions, and guide the future field measurement strategies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Permafrost thermal conditions are sensitive to shifts in snow timing

Abstract Changes in snow precipitation at high latitudes can significantly affect permafrost thermal conditions and thaw depth, potentially exposing more carbon-laden soil to microbial decomposition. A fully coupled process-based surface/subsurface thermal hydrology model with surface energy balance is used to analyze the impact of intra-annual variability in snow on permafrost thermal regime and the active layer thickness. In the four numerical scenarios considered, simulations were forced by the same meteorological data, except the snow precipitation, which was systematically altered to change timing of snowfall. The scenarios represent subtle shifts in snow timing, but the snow onset/melt days, the end of winter snowpack depth, and total annual snow precipitation are unchanged among scenarios. The simulations show that small shifts in the timing of snow accumulation can have significant effects on subsurface thermal conditions leading to active layer deepening and even talik formation when snowfall arrives earlier in the winter. The shifts in snow timing have a stronger impact on wetter regions, especially soil underneath small ponds, as compared to drained regions. This study highlights the importance of understanding potential changes in winter precipitation patterns for reliable projections of active-layer thickness in a changing Arctic climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Long‐term changes in soil carbon and nitrogen fractions in switchgrass, native grasses, and no‐till corn bioenergy production systems

Cellulosic bioenergy is a primary land-based climate mitigation strategy, with soil carbon (C) storage and nitrogen (N) conservation as important mitigation elements. Here, we present 13 years of soil C and N change under three cellulosic cropping systems: monoculture switchgrass (Panicum virgatum L.), a five native grasses polyculture, and no-till corn (Zea mays L.). Soil C and N fractions were measured four times over 12 years. Bulk soil C in the 0–25 cm depth at the end of the study period ranged from 28.4 (± 1.4 se) Mg C ha −1 in no-till corn, to 30.8 (± 1.4) Mg C ha −1 in switchgrass, and to 34.8 (± 1.4) Mg C ha −1 in native grasses. Mineral-associated organic matter (MAOM) ranged from 60% to 90% and particulate organic matter (POM) from 10% to 40% of total soil C. Over 12 years, total C as well as both C fractions persisted under no-till corn and switchgrass and increased under native grasses. In contrast, POM N stocks decreased 33% to 45% across systems, whereas MAOM N decreased only in no-till corn and by less than 13%. Declining POM N stocks likely reflect pre-establishment land use, which included alfalfa and manure in earlier rotations. Root production and large soil aggregate formation explained 69% (p < 0.001) and 36% (p = 0.024) of total soil C change, respectively, and 60% (p = 0.020) and 41% (p = 0.023) of soil N change, demonstrating the importance of belowground productivity and soil aggregates for producing and protecting soil C and conserving soil N. Differences between switchgrass and native grasses also indicate a dependence on plant diversity. Soil C and N benefits of bioenergy crops depend strongly on root productivity and pre-establishment land use.

09 BIOMASS FUELS↗
DOE OSTI2022

Dry Season Transpiration and Soil Water Dynamics in the Central Amazon

With current observations and future projections of more intense and frequent droughts in the tropics, understanding the impact that extensive dry periods may have on tree and ecosystem-level transpiration and concurrent carbon uptake has become increasingly important. Here, we investigate paired soil and tree water extraction dynamics in an old-growth upland forest in central Amazonia during the 2018 dry season. Tree water use was assessed via radial patterns of sap flow in eight dominant canopy trees, each a different species with a range in diameter, height, and wood density. Paired multi-sensor soil moisture probes used to quantify volumetric water content dynamics and soil water extraction within the upper 100 cm were installed adjacent to six of those trees. To link depth-specific water extraction patterns to root distribution, fine root biomass was assessed through the soil profile to 235 cm. To scale tree water use to the plot level (stand transpiration), basal area was measured for all trees within a 5 m radius around each soil moisture probe. The sensitivity of tree transpiration to reduced precipitation varied by tree, with some increasing and some decreasing in water use during the dry period. Tree-level water use scaled with sapwood area, from 11 to 190 L per day. Stand level water use, based on multiple plots encompassing sap flow and adjacent trees, varied from ∼1.7 to 3.3 mm per day, increasing linearly with plot basal area. Soil water extraction was dependent on root biomass, which was dense at the surface (i.e., 45% in the upper 5 cm) and declined dramatically with depth. As the dry season progressed and the upper soil dried, soil water extraction shifted to deeper levels and model projections suggest that much of the water used during the month-long dry-down could be extracted from the upper 2–3 m. Results indicate variation in rates of soil water extraction across the research area and, temporally, through the soil profile. These results provide key information on whole-tree contributions to transpiration by canopy trees as water availability changes. In addition, information on simultaneous stand level dynamics of soil water extraction that can inform mechanistic models that project tropical forest response to drought.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Soil Temperature and Moisture, Council Road Mile Marker 71, Seward Peninsula, Alaska, beginning 2016

Daily averages of soil temperature and moisture measured once every hour at different heights, as well as daily averages of hourly measured snow depths located at Intensive Monitoring Stations at Council Road Mile Marker 71 site. Deployed at each site is an Onset HOBO U30 data logger with five smart temperature sensors and three smart soil moisture sensors (10HS). Three sites (CN_IS_4, CN_IS_6A, CN_IS_7) are equipped with a snow depth sensor. Data are retrieved annually since 2016. Contains 53 *.CSV files including a file inventory list by year. Data files have header rows, NaN fields indicate invalid or missing data, and negative vertical offsets are above ground. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Detecting Sub-Micron Space Weathering Effects in Lunar Grains With Synchrotron Infrared Nanospectroscopy

Space weathering processes induce changes to the physical, chemical, and optical properties of space-exposed soil grains. For the Moon, space weathering causes reddening, darkening, and diminished contrast in reflectance spectra over visible and near-infrared wavelengths. The physical and chemical changes responsible for these optical effects occur on scales below the diffraction limit of traditional far-field spectroscopic techniques. Recently developed super-resolution spectroscopic techniques provide an opportunity to understand better the optical effects of space weathering on the sub-micrometer length scale. This paper uses synchrotron infrared nanospectroscopy to examine depth-profile samples from two mature lunar soils in the mid-infrared, 1,500 –700 cm -1 (6.7–14.3 μm). Overall, our findings are broadly consistent with prior bulk observations and theoretical models of space weathered spectra of lunar materials. These results provide a direct spatial link between the physical/chemical changes in space-exposed grain surfaces and spectral changes of space weathered bodies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Local-scale heterogeneity of soil thermal dynamics and controlling factors in a discontinuous permafrost region

In permafrost regions, the strong spatial and temporal variability in soil temperature cannot be explained by the weather forcing only. Understanding the local heterogeneity of soil thermal dynamics and their controls is essential to understand how permafrost systems respond to climate change and to develop process-based models or remote sensing products for predicting soil temperature. In this study, we analyzed soil temperature dynamics and their controls in a discontinuous permafrost region on the Seward Peninsula, Alaska. We acquired one-year temperature time series at multiple depths (at 5 or 10 cm intervals up to 85 cm depth) at 45 discrete locations across a 2.3 km 2 watershed. We observed a larger spatial variability in winter temperatures than that in summer temperatures at all depths, with the former controlling most of the spatial variability in mean annual temperatures. We also observed a strong correlation between mean annual ground temperature at a depth of 85 cm and mean annual or winter season ground surface temperature across the 45 locations. We demonstrate that soils classified as cold, intermediate, or warm using hierarchical clustering of full-year temperature data closely match their co-located vegetation (graminoid tundra, dwarf shrub tundra, and tall shrub tundra, respectively). We show that the spatial heterogeneity in soil temperature is primarily driven by spatial heterogeneity in snow cover, which induces variable winter insulation and soil thermal diffusivity. These effects further extend to the subsequent summer by causing variable latent heat exchanges. Finally, we discuss the challenges of predicting soil temperatures from snow depth and vegetation height alone by considering the complexity observed in the field data and reproduced in a model sensitivity analysis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Gap-filling eddy covariance methane fluxes: Comparison of machine learning model predictions and uncertainties at FLUXNET-CH4 wetlands

Time series of methane fluxes measured by eddy-covariance require gap-filling to estimate annual emissions. Gap-filling methane fluxes is challenging because of high variability and complex responses to multiple drivers. To date, there is no widely established gap-filling standard for methane, with regards both to the best model algorithms and predictors. In this study, we address the need for standardization by synthesizing results of gap-filling methods applied at 17 wetland sites spanning boreal to tropical regions including all major wetlands classes and two rice paddies. We introduce new procedures for: 1) creating realistic artificial gap scenarios, 2) training and evaluating gap-filling models without overstating performance, and 3) predicting half-hourly methane fluxes and annual emissions with robust uncertainty estimates. We tested a conventional method (marginal distribution sampling) and four machine learning algorithms - penalized linear regression, artificial neural networks, random forests, and boosted decision trees - and four predictor sets, including temporal, meteorological, ecosystem carbon and energy flux, and soil predictors. We find that the conventional method can achieve similar median performance to the machine learning models but is worse than the best machine learning models and relatively insensitive to predictor choices. Of the machine learning models, decision tree algorithms performed the best in cross-validation experiments, even with a baseline predictor set, and artificial neural networks showed comparable performance when using all predictors. Soil temperature was frequently the most important predictor whilst water table depth was important at sites with substantial water table fluctuations, highlighting the value of data on soil conditions. Raw gap-filling uncertainties from the machine learning models were underestimated and we propose a method to calibrate uncertainties to observations. Finally, we gap-fill and provide summary evaluation metrics for all 81 sites in the FLUXNET-CH4 community dataset and publicly release the python code for model development, evaluation, and uncertainty estimation.

42 ENGINEERING↗
DOE OSTI2021

Stability of tropical forest tree carbon-water relations in a rainfall exclusion treatment through shifts in effective water uptake depth

Increasing severity and frequency of drought is predicted for large portions of the terrestrial biosphere, with major impacts already documented in wet tropical forests. Using a 4-year rainfall exclusion experiment in the Daintree Rainforest in northeast Australia, we examined canopy tree responses to reduced precipitation and soil water availability by quantifying seasonal changes in plant hydraulic and carbon traits for 11 tree species between control and drought treatments. Even with reduced soil volumetric water content in the upper 1 m of soil in the drought treatment, we found no significant difference between treatments for predawn and midday leaf water potential, photosynthesis, stomatal conductance, foliar stable carbon isotope composition, leaf mass per area, turgor loss point, xylem vessel anatomy, or leaf and stem nonstructural carbohydrates. While empirical measurements of aboveground traits revealed homeostatic maintenance of plant water status and traits in response to reduced soil moisture, modeled belowground dynamics revealed that trees in the drought treatment shifted the depth from which water was acquired to deeper soil layers. Finally, these findings reveal that belowground acclimation of tree water uptake depth may buffer tropical rainforests from more severe droughts that may arise in future with climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Soil carbon stocks from soil organic carbon at the Teller 47, Teller 27, and Kougarok 86 field sites, Seward Peninsula, Alaska, 2021 – 2023

Carbon stored in permafrost soils represents the largest terrestrial storage and hence has key significance for the global carbon cycle, yet the site-specific estimates of permafrost soil carbon are limited. To better understand permafrost carbon storage, soil organic content (SOC) was measured from soils collected at the Teller 27, Teller 47, and Kougarok 86 field sites on the Seward Peninsula, Alaska in the summers of 2019 and 2021-2023. Soil samples were collected from locations with and without near-surface permafrost to varied depths within the permafrost, where present. In addition to sample collection for soil organic carbon analysis, bulk density samples were also collected within each soil profile to allow for calculation of carbon stocks. Carbon to nitrogen ratios are also reported. Location data were collected at each observation point using a handheld GPS unit with waypoint averaging to reduce location measurement error to within 3 m. SOC was measured using elemental analysis in the laboratory. This dataset contains a *.csv file of soil carbon stocks, a *.kml file of measurement locations, a *.pdf user guide, a *.csv data dictionary, and a *.csv file level metadata.

EARTH SCIENCE > CRYOSPHERE > FROZEN GROUND↗
DOE OSTI2024

Responses of Drying-rewetting (Transient Soil Moisture) and Steady State Soil Moisture Incubation on Soil Organic Carbon Dynamics in Three US Soils, 2017

This data set contains measurements of soil characteristics (aggregate size distribution and mean size, total aggregate associated carbon, extractable organic C, and microbial biomass C), microbial respiration, and soil metabolite concentrations from a transient and steady soil moisture incubation experiment using soils of different textures (sandy, loamy, and clayey). The study investigated mechanisms driving the Birch effect (increased carbon mineralization pulses with wetting following a drying period) in differing soil textures. Three different soils of distinctly different textures were collected from 0-15cm depth in Georgia (sandy, 2017-05-01), Missouri (loamy, 2017-06-14), and Texas (clayey, December 2017). Soils were incubated for 140 days with destructive harvests done on days 1, 29, 33, 56, 112, 116, and 140 in transient soil moisture incubation and on days 1, 33, 116, and 140 in steady state soil moisture incubation. This dataset contains six data files in comma separate (.csv) format. Additional metadata are provided: six data dictionaries and a file-level metadata file in comma separate (*.csv) format and a user guide in PDF (*.pdf) format.

1-Methyladenosine concentration↗
DOE OSTI2026

Soil Water Retention and Hydraulic Conductivity Data and Model at Trail Creek in Taylor River Watershed, Colorado 2024-2025

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors near Trail Creek. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format TR-X-Y, where TR refers to Trail Creek, X is the treatment block identifier, and Y is the location identifier. Specifically, TR-ASCC1 is the control treatment block under the Adaptive Silviculture for Climate Change (ASCC) project, and TR-ASCC2 is the clear-cut treatment block. TR-ASCC-EHSn is associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and TR-ASCC-ERTn (upslope n=1) are ecohydrology sites along the electrical resistivity tomography transects. The sample and location information can be found in metadata.csv, and the data from the soil sensors will be included in a future data version when the observation period becomes sufficiently long for data analysis. Sampling and Measurements Each sample falls into one of the two sampling methods – (1) intact cores or (2) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. The intact cores were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. 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.

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2021

New insights into the drainage of inundated ice-wedge polygons using fundamental hydrologic principles

Abstract. The pathways and timing of drainage from the inundated centers of ice-wedge polygons in a warming climate have important implications for carbon flushing, advective heat transport, and transitions from methane to carbon dioxide dominated emissions. Here, we expand on previous research using a recently developed analytical model of drainage from a low-centered polygon. Specifically, we perform (1) a calibration to field data identifying necessary model refinements and (2) a rigorous model sensitivity analysis that expands on previously published indications of polygon drainage characteristics. This research provides intuition on inundated polygon drainage by presenting the first in-depth analysis of drainage within a polygon based on hydrogeological first principles. We verify a recently developed analytical solution of polygon drainage through a calibration to a season of field measurements. Due to the parsimony of the model, providing the potential that it could fail, we identify the minimum necessary refinements that allow the model to match water levels measured in a low-centered polygon. We find that (1) the measured precipitation must be increased by a factor of around 2.2, and (2) the vertical soil hydraulic conductivity must decrease with increasing thaw depth. Model refinement (1) accounts for runoff from rims into the ice-wedge polygon pond during precipitation events and possible rain gauge undercatch, while refinement (2) accounts for the decreasing permeability of deeper soil layers. The calibration to field measurements supports the validity of the model, indicating that it is able to represent ice-wedge polygon drainage dynamics. We then use the analytical solution in non-dimensional form to provide a baseline for the effects of polygon aspect ratios (radius to thaw depth) and coefficient of hydraulic conductivity anisotropy (horizontal to vertical hydraulic conductivity) on drainage pathways and temporal depletion of ponded water from inundated ice-wedge polygon centers. By varying the polygon aspect ratio, we evaluate the relative effect of polygon size (width), inter-annual increases in active-layer thickness, and seasonal increases in thaw depth on drainage. The results of our sensitivity analysis rigorously confirm a previous analysis indicating that most drainage through the active layer occurs along an annular region of the polygon center near the rims. This has important implications for transport of nutrients (such as dissolved organic carbon) and advection of heat towards ice-wedge tops. We also provide a comprehensive investigation of the effect of polygon aspect ratio and anisotropy on drainage timing and patterns, expanding on previously published research. Our results indicate that polygons with large aspect ratios and high anisotropy will have the most distributed drainage, while polygons with large aspect ratios and low anisotropy will have their drainage most focused near their periphery and will drain most slowly. Polygons with small aspect ratios and high anisotropy will drain most quickly. These results, based on parametric investigation of idealized scenarios, provide a baseline for further research considering the geometric and hydraulic complexities of ice-wedge polygons.

Harp, Dylan R. (ORCID:0000000197778000)↗