SEARCH · Search NASA
Results for “soil aggregation”
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.
Microbial responses to long-term warming differ across soil microenvironments
Soil carbon loss is likely to increase due to climate warming, but microbiomes and microenvironments may dampen this effect. In a 30-year warming experiment, physical protection within soil aggregates affected the thermal responses of soil microbiomes and carbon dynamics. In this study, we combined metagenomic analysis with physical characterization of soil aggregates to explore mechanisms by which microbial communities respond to climate warming across different soil microenvironments. Long-term warming decreased the relative abundances of genes involved in degrading labile compounds (e.g. cellulose), but increased those genes involved in degrading recalcitrant compounds (e.g. lignin) across aggregate sizes. These changes were observed in most phyla of bacteria, especially for Acidobacteria, Actinobacteria, Bacteroidetes, Chloroflexi, and Planctomycetes. Microbial community composition was considerably altered by warming, leading to declined diversity for bacteria and fungi but not for archaea. Microbial functional genes, diversity, and community composition differed between macroaggregates and microaggregates, indicating the essential role of physical protection in controlling microbial community dynamics. Our findings suggest that microbes have the capacity to employ various strategies to acclimate or adapt to climate change (e.g. warming, heat stress) by shifting functional gene abundances and community structures in varying microenvironments, as regulated by soil physical protection.
Cover crops and poultry litter impact on soil structural stability in dryland soybean production in southeastern United States
Abstract This study explored the efficacy of soil aggregate indices in quantifying soil structural development, utilizing 5‐year field experiment data from the Southeastern United States. The experiment utilized a split‐plot design with cover crops (native vegetation as control, cereal rye (Secale cerealeL.), winter wheat (Triticum aestivum), hairy vetch (Vicia villosa), and mustard (Brassica rapa) plus cereal rye as the main factor and fertilizer source (no fertilizer as control, inorganic fertilizer with phosphorus, potassium, and elemental sulfur, and poultry litter) as the secondary factor. Aggregate size fractions were determined using the wet‐sieving method, and aggregate stability index (ASI), mean weight diameter (MWD), geometric mean diameter (GMD), and fractal dimension (FD) were calculated to assess soil structural stability. Main effects results indicated that cereal rye (55.11%) and poultry litter (50.97%) exhibited the highest ASI values. The highest MWD, GMD, and FD were observed under mustard plus cereal rye (1.187 mm), cereal rye (0.462 mm), and hairy vetch (2.573), respectively. Principal component analysis revealed that cover crops significantly improved soil aggregate structure and stability, overcoming limitations of sole fertilization practices. Regression analysis suggested that ASI, MWD, and GWD positively correlated with soil organic carbon, whereas FD negatively correlated with MWD, GMD, and ASI. Principal component analysis exhibited that FD decreased with increasing soil organic carbon, ASI, MWD, and GMD, demonstrating that lower FD values indicate enhanced soil aggregation and structure. Assessed indices, FD included, effectively gauged soil structural stability. These metrics should be prioritized in managerial decisions to support soil productivity and health in agricultural systems.
Emerging multiscale insights on microbial carbon use efficiency in the land carbon cycle
Microbial carbon use efficiency (CUE) affects the fate and storage of carbon in terrestrial ecosystems, but its global importance remains uncertain. Accurately modeling and predicting CUE on a global scale is challenging due to inconsistencies in measurement techniques and the complex interactions of climatic, edaphic, and biological factors across scales. The link between microbial CUE and soil organic carbon relies on the stabilization of microbial necromass within soil aggregates or its association with minerals, necessitating an integration of microbial and stabilization processes in modeling approaches. In this perspective, we propose a comprehensive framework that integrates diverse data sources, ranging from genomic information to traditional soil carbon assessments, to refine carbon cycle models by incorporating variations in CUE, thereby enhancing our understanding of the microbial contribution to carbon cycling.
Getting to the root of the problem: Soil carbon and microbial responses to root inputs within a buried paleosol along an eroding hillslope in southwestern Nebraska, USA
Large quantities of soil carbon (C) can persist within paleosols for millennia due to burial and subsequent isolation from plant-derived inputs, atmospheric conditions, and microbial activity at the modern surface. Erosion exposes buried soils to modern root-derived C influx via root exudation and root turnover, thus stimulating microbial activity leading to SOC decomposition and accumulation through organo-mineral stabilization of modern C. With this study we aim to quantify how modern root-derived C inputs impact paleosol C decomposition and stabilization across varying degrees of isolation from modern surface conditions in southwestern Nebraska, USA, where hillslope erosion is bringing a buried Late-Pleistocene-early Holocene paleosol (the “Brady Soil”) closer to the modern surface. We collected Brady Soil samples from 0.2m, 0.4m, and 1.2m below the modern surface and conducted two lab-based incubations. Soils were amended with either (1) a lab-synthesized mixture of low molecular weight compounds (12 atom% 13 C), or (2) 13 C enriched root residues (92 atom% 13 C), in 30-day and 240-day incubation experiments, respectively. Here we determined microbial responses to synthetic root exudates and residues by partitioning the 13 C label from Brady Soil C, including measurements of total, root, and primed C respiration, microbial biomass C (MBC), microbial C use efficiency (CUE). To assess the capacity of isolated paleosols to accrue modern plant C, we used Nano-scale Secondary Ion Mass Spectrometry imaging. We found that: (1) adding root-derived C inputs primed Brady Soil C across all depths, and was mediated by depth and composition of root additions; (2) root-derived C inputs stimulated microbial biomass C (MBC) growth similarly across depths, but the magnitude of CUE and MBC varied by chemistry of root-derived additions; (3) new particulate organic matter was incorporated into mineral-associated pools over time; (4) material from the added root residues was found in association with bacterial cells and fungal hyphae as well as with soil aggregate and mineral surfaces. Our study shows that paleosols defy expectations of C content and reactivity with depth, and changes in land cover and climate will expose buried paleosols to modern surface conditions, increasing respired C. This work highlights the importance of evaluating the role resurfacing buried soils through landscape change plays in C cycle feedbacks to the climate system.
Data for: Miscanthus × giganteus changes soil structure and increases maximum water holding capacity
The data provided include results from a comparative study evaluating the impact of Miscanthus × giganteus (miscanthus) versus maize on soil structural properties and maximum water holding capacity (MWHC) across two Iowa sites. The dataset includes MWHC values determined using the Funnel Filter Paper Drainage (FFPD) method, as well as additional measurements of MWHC following structural disruption of the soil to isolate the effect of aggregation. It also contains three-dimensional micro-computed tomography (microCT) data used to quantify total porosity and pore size distribution (PSD) of soil aggregates at a 5 µm resolution. All data are provided as raw replicate-level measurements, organized by site, crop, and depth, along with processed summary files in table form in CSV (.csv) format to support reproducibility and downstream analysis.
WHONDRS laboratory time series moisture manipulative experiment from soil core layers across eastern contiguous US: time series aerobic respiration, geochemistry, and aggregates
This dataset supports a broader study examining the effects of wetting and drying on soil layers across the eastern contiguous United States (CONUS). The dataset provides data generated from a laboratory moisture manipulation experiment. The contents include time series aerobic respiration and moisture; dissolved oxygen; sediment geochemistry data; and field metadata. Samples were collected as part of a collaboration between WHONDRS (Worldwide Hydrobiogeochemistry Observation Network for Dynamic River Systems; https://whondrs.pnnl.gov) and MONet (Molecular Observation Network; https://www.emsl.pnnl.gov/monet). The field samples (soil cores) were labeled as MEL_##_COR and subsequent subsamples begin with MEL_##. Additional subsamples were taken for the laboratory experiment and were labeled as EL_##. The labels from the MEL field samples and the EL subsamples can be mapped directly based on the digits following the prefix and underscore (i.e., EL_01 is a subsample from MEL_01). See the critical details section below for more details on sample naming and experimental design.For details on how to navigate this data package, see this infographic from the River Corridor SFA https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About.In addition to this readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes all column/row headers and variable definitions.This dataset is comprised of (1) a folder containing environmental context photos; (2) file-level metadata; (3) data dictionary; (4) field metadata; (5) readme; (6) international generic sample number (IGSN) mapping file; and (7) a subfolder with soil sample data from field samples and the incubation experiment. The sample data subfolder contains (1) effect size; (2) gravimetric moisture from field samples and incubation experiment; (3) respiration rates, raw dissolved oxygen values, and plots; (4) specific conductance, pH, and temperature from the incubation; (5) soil aggregates; (6) a summary containing median values of each data type for each treatment (wet and dry) in the incubation; (7) a summary containing averages for each data type of each soil layer; and (8) methods codes. All files are .csv, .pdf, .jpeg, or .jpg.
Dominant Controls on Preferential Flow and Their Implications for Future Soil Water Fluxes
Abstract Soil water flow, particularly preferential flow (PF), is a critical control on hydrological and biogeochemical processes, including groundwater recharge, contaminant transport, and carbon cycling. However, it remains challenging to predict PF occurrence across large environmental gradients. Here, we developed a deep learning (DL) model to estimate event‐scale soil water flow velocity and the probability of PF occurrence using high‐frequency soil moisture and precipitation data from 33 sites across the National Ecological Observatory Network. The model demonstrated high skill in predicting the binary occurrence of PF (91% F1‐score; 85% accuracy) but the performance was limited in predicting soil water velocity ( R 2 = 0.31). We found that precipitation characteristics (duration, volume, and intensity) were the most important predictors for soil water velocity. Among the non‐precipitation event variables, sand content showed relatively high predictive skill, though differences among non‐event climate variables were generally modest. Lower sand content was associated with increased predicted soil water velocity, a finding that highlights the role of soil structure in producing more non‐uniform flow, which contrasts with traditional uniform flow models. Projecting a reduced DL model under both moderate and high‐emissions future climate scenarios (2060–2099 Representative Concentration Pathways 4.5 and 8.5), we found ∼7.3% increase under RCP4.5 and ∼15% under RCP8.5 of soil water velocities compared to the historical simulation, while modeled likelihood of PF changed little. These findings suggest climate change is not making PF more frequent, but it is making existing PF pathways more efficient with important consequences for associated nutrient and contaminant transport under climate change. Plain Language Summary Water movement in soil is critical for water quality. While often modeled as a uniform flow process, in reality water moves rapidly through cracks and burrows in what is called “preferential flow” (PF), which limits natural filtration and can transport pollutants. We developed a deep learning model, trained on data from 33 U.S. sites, to predict when and how fast this PF occurs based on precipitation, soil, and climate data. The model showed that precipitation characteristics (duration, intensity, volume) were the most important predictors of PF. Lower soil sand content/higher clay content was associated with faster water flow, likely due to clay soils forming aggregates and cracks that water moves through rather than infiltrating uniformly. Further analyses based on climate projections suggest that the speed at which PF occurs will become more rapid under future climate scenarios compared to historical simulation. This highlights the need to represent PF in soil water models when assessing future water quality. Key Points The effect of precipitation peak intensity on soil water velocities declined with increasing precipitation intensity Antecedent soil moisture failed to predict preferential flow (PF), contrasting the high predictive power of sand content Climate predictions suggest that soil water velocities through PF paths will increase ∼15% by 2099
A framework for testing soil carbon dynamics post land-use transition in a multisector dynamics model
Soil carbon plays a crucial role in the global carbon cycle. Changes in land use can determine whether carbon is stored or is emitted into the atmosphere as carbon dioxide, which has broad implications for the human and Earth systems. These feedbacks to the carbon cycle and their socio-economic drivers are modelled by many global multisector dynamics models to project future possibilities for the human-Earth system. One notable model of this class is the Global Change Analysis Model (GCAM), which uses a simplified process to model soil organic carbon (SOC) content after land-use transition across 384 land units. While the current GCAM soil carbon framework is based on scientific principles, it has not been tested against experimental data. This work examines rates of SOC change from GCAM input data. Specifically, first order rate constants derived from model inputs were compared to values from two syntheses to assess GCAM’s accuracy. Welch’s t-tests and linear models were used to determine if rate constants were consistent across all tested geographical areas and land-use transition types. While we found that there was general agreement on the direction and magnitude (i.e., rate) of SOC change, the rate constant derived from GCAM and empirical values differed strongly in a subset of specific instances. These results indicate that GCAM’s current SOC dynamics during land use transition successfully capture broad patterns of change in this critical carbon pool, but should be interpreted with caution at finer spatial scales. One potential cause of these discrepancies is our highly aggregated variable, soil timescale, which could be made more granular to improve accuracy. When using economically rooted multisector dynamics models, such as GCAM, it is critical to understand such model limitations for representing specific Earth system processes.
Comparison of automated chemical-guided segmentation and human annotation of soil organic matter in X-ray microcomputed tomography imaging in contrasted soil types
Soil organic matter (OM) formation and persistence is strongly influenced by the spatial distribution of organic substrates and microscale soil heterogeneity by dictating OM accessibility to microorganisms. However, traditional size and/or density fractionation techniques disrupt aggregate architecture, eliminating spatial information needed to fully understand intra-aggregate OM distribution. To quantify three-dimensional OM spatial distribution and automate segmentation in X-ray microcomputed tomography (µCT) imaging without human annotation bias, we developed an iodine gas vapor (I2) based staining workflow that eliminates labor-intensive manual annotation while maintaining segmentation accuracy, using aggregates from four taxonomically diverse soils (Xerofluvent, Haploxeroll Sphagnofibrist, Palehumult) with an 8-fold range of soil organic carbon. Human annotation of 10 µCT slices by the experienced and inexperienced annotators resulted in variations up to 3% in the Dice similarity coefficient (DSC), reflecting a degree of inherent subjectivity of manual labeling. Such inconsistencies are expected to compound as the number of manually annotated slices increases. Dual-energy µCT imaging at 33.1 keV (below the iodine (I) K-edge) and 33.2 keV (above the I K-edge) was used to resolve aggregate microstructure following I2 staining. The automated image subtraction pipeline identified OM regions by the I Kedge induced brightness increases, achieving DSC values of 0.58–0.83 relative to an experienced annotator. Sensitivity analyses revealed that the reconstruction alpha value—optimized via the open-source tool TomocuPy—and the 3D registration slice count were the primary determinants of accuracy, providing a novel benchmark for dual-energy soil imaging. The pipeline without GPU acceleration achieved 9.6 to 43.2 times faster than manual annotation. Using GPU-accelerated image post-processing and affine transformation matrices, the pipeline successfully segmented OM elements for large-scale datasets (3232×3232 pixel, 2048 slices) within ~5200 s from raw file acquisition to segmented output. The high-throughput approach enables the quantification of OM spatial distribution across diverse and heterogeneous soil.
Composition and persistence of soil organic matter along eroding and depositional transects in buried vs. modern soil layers: A case of the Brady paleosol at Wauneta, Nebraska
Paleosols form when soils are buried through deposition by aeolian, colluvial, alluvial or other processes. Burial of former topsoil isolates soil organic matter (SOM) from surface conditions, allowing carbon to accumulate and potentially remain stable for millennia. In this study, SOM composition, distribution, and persistence were analyzed in the Brady Soil of Nebraska, USA to compare SOM spatial variability in modern and buried soils, as well as the impact of erosional exposure on SOM stability. The Brady Soil, formed as a surface soil during the Pleistocene-Holocene transition and now a paleosol buried up to 6 m deep (or more) by loess deposition during the Holocene, was sampled along burial (up to 5.8 m depth) and erosional (up to 1.8 m depth) transects to compare SOM dynamics in different geomorphic settings. Fourier Transform Infrared Spectroscopy (FTIR) and Fourier Transform ion cyclotron resonance mass spectrometry (FTICR-MS) were used to analyze SOM composition, while δ 13 C isotope analyses identified SOM sources and radiocarbon values were used to estimate turnover rates. Results confirmed a vegetation shift from C3 to C4 plants after Brady Soil formation, reflecting warming climatic conditions. Increasing SOM age and decreasing δ 13 C and δ 15 N values with depth indicated slowing of decomposition rate in buried soils. Higher pH in the Brady Soil suggested greater base cation content, supporting SOM stabilization through organo-mineral associations and aggregate formation. However, exposure of the Brady Soil due to surface erosion caused faster SOM turnover. This result suggested susceptibility of buried SOM to losses via decomposition upon erosional exposure, possibly accelerated by priming in response to modern SOM inputs. These findings highlight the potential loss of carbon stocks in buried soils under future climate change, as shifts in soil physicochemical properties may destabilize long-preserved SOM.
Linking Plant and Microbial Traits to Soil Carbon for Reliable and Resilient Bioenergy Systems
Bioenergy systems in the United States offer a dual opportunity to supply renewable feedstocks while enhancing ecosystem services such as hydrologic regulation, erosion control, and soil carbon (C) storage. National assessments highlight the potential to grow perennial energy crops to improve soil function and ecosystem resilience. Realizing this potential requires understanding the ecological mechanisms that govern how C is added, transformed, and stabilized in soils. Plant traits determine the quantity, depth, and chemistry of organic inputs, while microbial processes—including carbon use efficiency, necromass formation, and trophic interactions—mediate their transformation and partitioning among soil carbon pools. These biological pathways are shaped by soil physical and chemical properties, including aggregation, texture, and mineralogy, and by environmental drivers such as temperature, moisture, and disturbance, leading to context-dependent outcomes across landscapes. Management practices that diversify feedstocks, minimize disturbance, and maintain soil cover can promote both biomass production and C retention, while microbial amendments and rhizosphere engineering offer emerging, but often context-dependent, tools to optimize plant–microbe interactions. Trade-offs between biomass yield and soil carbon storage may arise when systems favor rapid aboveground productivity at the expense of belowground inputs and microbial processing, underscoring the importance of trait combinations that support both functions. Advances in monitoring, reporting, and verification—spanning precision agriculture, remote sensing, and biosensing—are improving predictive capacity through microbial-explicit process models and model–experiment (ModEx) frameworks. By connecting soil, plant, and microbial processes with advances in modeling and biosensing, this review outlines research priorities focused on trait-based parameterization and ModEx integration. These priorities will support the design of bioenergy systems that are both reliable and resilient, enhancing renewable energy production and ecosystem sustainability.
Non-Destructive, Three-Dimensional Imaging of Processes in the Rhizosphere Utilizing High Energy Photons
Soil structure, which can be described as the aggregation and distribution of pore spaces, regulates carbon, nutrient, and water cycling across the Earth system. Yet the inability to make quantitative, dynamic, in situ measurements of soil structure and rhizosphere carbon flow has prevented meaningful incorporation of soil structural processes into Earth System Models (ESMs). The key limitations are: (i) Lack of scale integration between micron-scale soil structure and ecosystem-scale models, (ii) Poor functional linkage between soil structural properties and biogeochemical processes, and (iii) Absence of dynamic 3D measurements of rhizosphere structural changes and carbon transformations. To address these challenges, we developed a new integrated positron emission tomography (PET) - microcomputed tomography (CT) imaging platform that enables the first 4D (spatial 3D and time), non-invasive, quantitative imaging of carbon allocation and rhizosphere structural dynamics in living plants and intact soils. The system combines: • Rhizo-PET (R-PET): a high-resolution positron emission tomography scanner optimized for 11 CO 2 tracing in plant roots. • a-Se micro-CT: a high-contrast, high-spatial resolution CT system based on amorphous selenium (a-Se) direct conversion technology, enabling micron-scale visualization of soil structure and root–soil interfaces. Together, these advances allow us to quantify how carbon exudates move, transform, and stabilize within the rhizosphere, directly informing missing processes in BER-relevant carbon cycle models.
Coarse-grained simulation of colloidal self-assembly, cation exchange, and rheology in Na/Ca smectite clay gels
Knowledge Gap: The aggregation of clay minerals—layered silicate nanoparticles—strongly impacts fluid flow, solute migration, and solid mechanics in soils, sediments, and sedimentary rocks. Experimental and computational characterization of clay aggregation is inhibited by the delicate water-mediated nature of clay colloidal interactions and by the range of spatial scales involved, from 1 nm thick platelets to flocs with dimensions up to micrometers or more. Simulations: Using a new coarse-grained molecular dynamics (CGMD) approach, we predicted the microstructure, dynamics, and rheology of hydrated smectite (more precisely, montmorillonite) clay gels containing up to 2,000 clay platelets on length scales up to 0.1 μm. Further, simulations investigated the impact of simulation time, platelet diameters (6 to 25nm), and the ratio of Na to Ca exchangeable cations on the assembly of tactoids (i.e., stacks of parallel clay platelets) and larger aggregates (i.e., assemblages of tactoids). We analyzed structural features including tactoid size and size distribution, basal spacing, counterion distribution in the electrical double layer, clay association modes, and the rheological properties of smectite gels. Findings: Our results demonstrate new potential to characterize and understand clay aggregation in dilute suspensions and gels on a scale of thousands of particles with explicit representation of counterion clouds and with accuracy approaching that of all-atom molecular dynamics (MD) simulations. For example, our simulations predict the strong impact of Na/Ca ratio on clay tactoid formation and the shear-thinning rheology of clay gels.
Horizontal transfer of chromosomal DNA mediated by an integrative and conjugative element generates frequent localized recombination in Novosphingobium aromaticivorans
Horizontal gene transfer is an important evolutionary process by which DNA is exchanged between cells that are physically co-located but not direct evolutionary descendants. Horizontal transfer of highly divergent DNA is relatively easy to detect and can produce major phenotypic changes, exemplified by the acquisition of antibiotic resistance determinants. However, transfer of high-identity DNA, for example, between strains of the same species, is likely to be more frequent, harder to detect, and highly impactful in aggregate. In this work, we demonstrate that soil isolates of the alphaproteobacterium Novosphingobium aromaticivorans can exchange chromosomal DNA, leading to multiple unselected recombination events spanning approximately 10% of the chromosome. Chromosomal recombination was directional and more efficient near an integrative and conjugative element (ICE), and required a relaxase found in the ICE. Recombination could not be observed in strains from closely related Novosphingobium species. In combination, these results suggest that ICE-mediated recombination can efficiently recombine DNA within N. aromaticivorans, increasing the adaptive potential of the species while also enforcing species boundaries through preferential intraspecific recombination.
Particle dynamics of nanoplastics suspended in water with soil microparticles: insights from small angle neutron scattering (SANS) and ultra-SANS
Small-angle neutron scattering (SANS) and ultra-SANS (USANS) were employed to understand the aggregation behavior and observe the size reduction for nanoplastics (NPs) formed from a biodegradable mulch film, and microparticles of vermiculite (V), an artificial soil, suspended in water in the presence of low convective shear (ex situ stirring) prior to measurements. Neutron contrast matching was employed to minimize the signal of V (by 100-fold) and thereby isolate the signal due to NPs in the neutron beam, as the contrast match point (CMP) for V (67 vol% deuteration of water) differed from that of NPs by more than 20%. The original NPs' size distribution was bimodal: <200 nm and 500–1200 nm, referred to as small and large NPs, i.e., SNPs and LNPs, respectively. In the absence of V, SNPs formed homoaggregates at higher concentrations that decreased with stirring time, while the size of LNPs remained unchanged. The presence of V at 2-fold lower concentration than NPs did not change the size of SNPs but reduced the size of LNPs by nearly 2-fold as stirring time increased. Because the size of SNPs and LNPs did not differ substantially between CMP and 100% D 2 O solvents, it is evident that SNPs and LNPs are mainly composed of NPs and not V. In conclusion, the results suggest that LNPs are susceptible to size reduction through collisions with soil microparticles via convection, yielding SNPs near soil–water interfaces within vadose zones.
Divide and conquer: using RhizoVision Explorer to aggregate data from multiple root scans using image concatenation and statistical methods
Roots are important in agricultural and natural systems for determining plant productivity and soil carbon inputs. Sometimes, the amount of roots in a sample is too much to fit into a single scanned image, so the sample is divided among several scans, and there is no standard method to aggregate the data. Here, we describe and validate two methods for standardizing measurements across multiple scans: image concatenation and statistical aggregation. We developed a Python script that identifies which images belong to the same sample and returns a single, larger concatenated image. These concatenated images and the original images were processed with RhizoVision Explorer, a free and open-source software. An R script was developed, which identifies rows of data belonging to the same sample and applies correct statistical methods to return a single data row for each sample. These two methods were compared using example images from switchgrass, poplar, and various tree and ericaceous shrub species from a northern peatland and the Arctic. Most root measurements were nearly identical between the two methods except median diameter, which cannot be accurately computed by statistical aggregation. We believe the availability of these methods will be useful to the root biology community.
Spatially explicit terrestrial carbon densities for calibrating the carbon cycle in human-Earth system Models
Soil and vegetation carbon stocks play a critical role in human-Earth system models. These stocks (denominated as densities in MgC/ha) affect variables such as land use change emissions and also influence land use change pathways under climate forcing scenarios where terrestrial carbon is assigned a carbon price. Here we present reharmonized soil and vegetation carbon densities both at the 5-arcmin resolution grid cell level and also aggregated to 235 water sheds for 4 land use types (Cropland, Grazed land, Urban land and unmanaged vegetation) and 15 unmanaged land cover types. Moreover, we use the distribution of carbon within and across pixels to define statistical "states" of carbon, once again differentiated by land type. These statistical states are used to define a range of possible carbon values that can be used for defining initial conditions of soil and vegetation carbon in human-Earth system models. We implement these data in a state-of-the-art multi sector dynamics model, namely the Global Change Analysis Model (GCAM), and show that these new data improve several land use responses, especially when terrestrial carbon is assigned a carbon price.