Search NASA⌕ Search

Engineering topics

Sulman, Benjamin N.

Publications and source records attributed to Sulman, Benjamin N..

At least 19 records

Effects of fire and fire-induced changes in soil properties on post-burn soil respiration

Boreal forests cover vast areas of land in the northern hemisphere and store large amounts of carbon (C) both aboveground and belowground. Wildfires, which are a primary ecosystem disturbance of boreal forests, affect soil C via combustion and transformation of organic matter during the fire itself and via changes in plant growth and microbial activity post-fire. Wildfire regimes in many areas of the boreal forests of North America are shifting towards more frequent and severe fires driven by changing climate. As wildfire regimes shift and the effects of fire on belowground microbial community composition are becoming clearer, there is a need to link fire-induced changes in soil properties to changes in microbial functions, such as respiration, in order to better predict the impact of future fires on C cycling. We used laboratory burns to simulate boreal crown fires on both organic-rich and sandy soil cores collected from Wood Buffalo National Park, Alberta, Canada, to measure the effects of burning on soil properties including pH, total C, and total nitrogen (N). We used 70-day soil incubations and two-pool exponential decay models to characterize the impacts of burning and its resulting changes in soil properties on soil respiration. Laboratory burns successfully captured a range of soil temperatures that were realistic for natural wildfire events. We found that burning increased pH and caused small decreases in C:N in organic soil. Overall, respiration per gram total (post-burn) C in burned soil cores was 16% lower than in corresponding unburned control cores, indicating that soil C lost during a burn may be partially offset by burn-induced decreases in respiration rates. Simultaneously, burning altered how remaining C cycled, causing an increase in the proportion of C represented in the modeled slow-cycling vs. fast-cycling C pool as well as an increase in fast-cycling C decomposition rates. Together, our findings imply that C storage in boreal forests following wildfires will be driven by the combination of C losses during the fire itself as well as fire-induced changes to the soil C pool that modulate post-fire respiration rates. Moving forward, we will pair these results with soil microbial community data to understand how fire-induced changes in microbial community composition may influence respiration.

54 ENVIRONMENTAL SCIENCES↗

Simulated plant-mediated oxygen input has strong impacts on fine-scale porewater biogeochemistry and weak impacts on integrated methane fluxes in coastal wetlands

Methane (CH 4 ) emissions from wetland ecosystems are controlled by redox conditions in the soil, which are currently underrepresented in Earth system models. Plant-mediated radial oxygen loss (ROL) can increase soil O 2 availability, affect local redox conditions, and cause heterogeneous distribution of redox-sensitive chemical species at the root scale, which would affect CH 4 emissions integrated over larger scales. In this study, we used a subsurface geochemical simulator (PFLOTRAN) to quantify the effects of incorporating either spatially homogeneous ROL or more complex heterogeneous ROL on model predictions of porewater solute concentration depth profiles (dissolved organic carbon, methane, sulfate, sulfide) and column integrated CH 4 fluxes for a tidal coastal wetland. From the heterogeneous ROL simulation, we obtained 18% higher column averaged CH 4 concentration at the rooting zone but 5% lower total CH 4 flux compared to simulations of the homogeneous ROL or without ROL. This difference is because lower CH 4 concentrations occurred in the same rhizosphere volume that was directly connected with plant-mediated transport of CH 4 from the rooting zone to the atmosphere. Sensitivity analysis indicated that the impacts of heterogeneous ROL on model predictions of porewater oxygen and sulfide concentrations will be more important under conditions of higher ROL fluxes or more heterogeneous root distribution (lower root densities). Despite the small impact on predicted CH 4 emissions, the simulated ROL drastically reduced porewater concentrations of sulfide, an effective phytotoxin, indicating that incorporating ROL combined with sulfur cycling into ecosystem models could potentially improve predictions of plant productivity in coastal wetland ecosystems.

54 ENVIRONMENTAL SCIENCES↗

Modeling Interactive Effects of Manganese Bioavailability, Nitrogen Deposition, and Warming on Soil Carbon Storage

Manganese (Mn) is a redox–active micronutrient that has been shown to accelerate plant litter decomposition; however, the effect of Mn–promoted decomposition on soil C storage is unclear. We present a novel biogeochemical model simulating how Mn bioavailability influences soil organic C (SOC) stocks in a soil profile (<50 cm) within a temperate forest. In our model, foliar Mn increased in response to increasing soluble Mn released through Mn–oxide (birnessite) dissolution in mineral soil layers. The ensuing Mn enrichment in leaf litter redistributed Mn to the surface forest floor layer, promoted enzymatic oxidation of lignin, and decreased SOC stocks. Total SOC loss was partially mitigated by accumulation of lignin–oxidation products as mineral–associated organic C. We also explored how Mn–driven changes to C storage interacted with effects of N deposition and warming. Nitrogen enrichment inhibited Mn–dependent lignin degradation, increasing SOC stocks and weakening their dependence on Mn bioavailability. Warming stimulated decomposition and reduced C stocks but was less effective at low Mn bioavailability. Our model results suggest that SOC stocks are sensitive to Mn bioavailability because increased plant uptake redistributes Mn to surface soils where it can enhance litter decomposition. Finally, based on our simulations, we predict that Mn becomes limiting to litter decomposition where Mn is poorly soluble. Depletion of bioavailable Mn or other cofactors that are critical to decomposition could limit the response of organic C stocks to warming over time, but quantitative projections of the role of Mn bioavailability in regulating decomposition requires additional measurements to constrain model uncertainties.

58 GEOSCIENCES↗

Integrating Tide‐Driven Wetland Soil Redox and Biogeochemical Interactions Into a Land Surface Model

Abstract Redox processes, aqueous and solid‐phase chemistry, and pH dynamics are key drivers of subsurface biogeochemical cycling and methanogenesis in terrestrial and wetland ecosystems but are typically not included in terrestrial carbon cycle models. These omissions may introduce errors when simulating systems where redox interactions and pH fluctuations are important, such as wetlands where saturation of soils can produce anoxic conditions and coastal systems where sulfate inputs from seawater can influence biogeochemistry. Integrating cycling of redox‐sensitive elements could therefore allow models to better represent key elements of carbon cycling and greenhouse gas production. We describe a model framework that couples the Energy Exascale Earth System Model (E3SM) Land Model (ELM) with PFLOTRAN biogeochemistry, allowing geochemical processes and redox interactions to be integrated with land surface model simulations. We implemented a reaction network including aerobic decomposition, fermentation, sulfate reduction, sulfide oxidation, methanogenesis, and methanotrophy as well as pH dynamics along with iron oxide and iron sulfide mineral precipitation and dissolution. We simulated biogeochemical cycling in tidal wetlands subject to either saltwater or freshwater inputs driven by tidal hydrological dynamics. In simulations with saltwater tidal inputs, sulfate reduction led to accumulation of sulfide, higher dissolved inorganic carbon concentrations, lower dissolved organic carbon concentrations, and lower methane emissions than simulations with freshwater tidal inputs. Model simulations compared well with measured porewater concentrations and surface gas emissions from coastal wetlands in the Northeastern United States. These results demonstrate how simulating geochemical reaction networks can improve land surface model simulations of subsurface biogeochemistry and carbon cycling.

54 ENVIRONMENTAL SCIENCES↗

Reparameterizing Litter Decomposition Using a Simplified Monte Carlo Method Improves Litter Decay Simulated by a Microbial Model and Alters Bioenergy Soil Carbon Estimates

Litter decomposition determines soil organic matter (SOM) formation and plant-available nutrient cycles. Therefore, accurate model representation of litter decomposition is critical to improving soil carbon (C) projections of bioenergy feedstocks. Soil C models that simulate microbial physiology (i.e., microbial models) are new to bioenergy agriculture, and their parameterization is often based on small datasets or manual calibration to reach benchmarks. Here, we reparameterized litter decomposition in a microbial soil C model (CORPSE - Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment) using the continental-scale Long-term Inter-site Decomposition Experiment Team (LIDET) dataset which documents decomposition across a range of litter qualities over a decade. Here, we conducted a simplified Monte Carlo simulation that constrained parameter values to reduce computational costs. The LIDET-derived parameters improved modeled C and nitrogen (N) remaining, decomposition rates, and litter mean residence times as compared to Baseline parameters. We applied the LIDET litter decomposition parameters to a microbial bioenergy model (Fixation and Uptake of Nitrogen – Bioenergy Carbon, Rhizosphere, Organisms, and Protection) to examine soil C estimates generated by Baseline and LIDET parameters. LIDET parameters increased estimated soil C in bioenergy feedstocks, with even greater increases under elevated plant inputs (i.e., by increasing residue, N fertilization). This was due to the integrated effects of plant litter quantity, quality, and agricultural practices (tillage, fertilization). Collectively, we developed a simple framework for using large-scale datasets to inform the parameterization of microbial models that impacts projections of soil C for bioenergy feedstocks.

54 ENVIRONMENTAL SCIENCES↗

Developing a Redox Network for Coastal Saltmarsh Systems in the PFLOTRAN Reaction Model

Coastal ecosystems have been largely ignored in Earth system models but are important zones for carbon and nutrient processing. Interactions between water, microbes, soil, sediments, and vegetation are important for mechanistic representation of coastal processes and ecosystem function. To investigate the role of these feedbacks, we used a reactive transport model (PFLOTRAN) that has the capability to be connected to the Energy Exascale Earth System Model (E3SM). PFLOTRAN was used to incorporate redox reactions and track chemical species important for coastal ecosystems as well as define simple representations of vegetation dynamics. Here our goal was to incorporate oxygen flux, salinity, pH, sulfur cycling, and methane production along with plant-mediated transport of gases and tidal flux. Using porewater profile and incubation data for model calibration and evaluation, we were able to create depth-resolved biogeochemical soil profiles for saltmarsh habitat and use this updated representation to simulate direct and indirect effects of elevated CO 2 and temperature on subsurface biogeochemical cycling. We found that simply changing the partial pressure of CO 2 or increasing temperature in the model did not fully reproduce observed changes in the porewater profile, but the inclusion of plant or microbial responses to CO 2 and temperature manipulations was more accurate in representing porewater concentrations. This indicates the importance of characterizing tightly coupled vegetation-subsurface processes for developing predictive understanding and the need for measurement of plant-soil interactions on the same time scale to understand how hotspots or moments are generated.

54 ENVIRONMENTAL SCIENCES↗

Subsurface Redox Interactions Regulate Ebullitive Methane Flux in Heterogeneous Mississippi River Deltaic Wetland

Abstract As interfaces connecting terrestrial and ocean ecosystems, coastal wetlands develop temporally and spatially complex redox conditions, which drive uncertainties in greenhouse gas emission as well as the total carbon budget of the coastal ecosystem. To evaluate the role of complex redox reactions in methane emission from coastal wetlands, a coupled reactive‐transport model was configured to represent subsurface biogeochemical cycles of carbon, nitrogen, and sulfur, along with production and transport of multiple gas species through diffusion and ebullition. This model study was conducted at multiple sites along a salinity gradient in the Barataria Basin at the Mississippi River Deltaic Plain. Over a freshwater to saline gradient, simulated total flux of methane was primarily controlled by its subsurface production and consumption, which were determined by redox reactions directly (e.g., methanogenesis, methanotrophy) and indirectly (e.g., competition with sulfate reduction) under aerobic and/or anaerobic conditions. At fine spatiotemporal scales, surface methane fluxes were also strongly dependent on transport processes, with episodic ebullitive fluxes leading to higher spatial and temporal variability compared to the gradient‐driven diffusion flux. Ebullitive methane fluxes were determined by methane fraction in total ebullitive gas and the frequency of ebullitive events, both of which varied with subsurface methane concentrations and other gas species. Although ebullition thresholds are constrained by local physical factors, this study indicates that redox interactions not only determine gas composition in ebullitive fluxes but can also regulate ebullition frequency through gas production.

03 NATURAL GAS↗

Soil Respiration Responses to Throughfall Exclusion Are Decoupled From Changes in Soil Moisture for Four Tropical Forests, Suggesting Processes for Ecosystem Models

Abstract Climatic drying is predicted for many tropical forests yet models remain poorly parameterized for these ecosystems, hampering predictions of forest‐climate interactions. We applied an integrated model–experiment approach, parameterizing an ecosystem model with tropical forest observational data and comparing model predictions to a field drying manipulation. We hypothesized that drying suppresses soil CO 2 fluxes (i.e., respiration) in relatively dry tropical forests but increases CO 2 fluxes in wetter tropical forests by alleviating anaerobiosis. We measured soil CO 2 fluxes during wet‐dry cycles from 2015 to 2022 in four Panamanian forests that vary in rainfall and soil fertility. Measured soil CO 2 fluxes declined in the dry season and peaked in the early wet season ahead of peak soil moisture, resulting in lower soil moisture optima for respiration than previously modeled. We then parameterized the model using field data and the new moisture‐respiration response functions. The updated model predicted increased soil CO 2 fluxes with drying in wetter and fertile forests and suppressed fluxes in drier, infertile forests. In contrast to model predictions, a chronic throughfall exclusion experiment initially suppressed soil respiration across forests, with sustained suppression for four years in the wettest forest only (−28% ± 4% during the dry season). In the fertile forest, drying eventually elevated CO 2 fluxes over this period (+75% ± 28% during the late wet season). The unexpected negative drying effect in the wettest, infertile forest could have resulted from reduced vertical flushing of nutrients into soils. Including hydro‐nutrient interactions in ecosystem models could improve predictions of tropical forest‐climate feedbacks.

54 ENVIRONMENTAL SCIENCES↗

Simulated Hydrological Dynamics and Coupled Iron Redox Cycling Impact Methane Production in an Arctic Soil

The fate of organic carbon (C) in permafrost soils is important to the climate system due to the large global stocks of permafrost C. Thawing permafrost can be subject to dynamic hydrology, making redox processes an important factor controlling soil organic matter (SOM) decomposition rates and greenhouse gas production. In iron (Fe)-rich permafrost soils, Fe(III) can serve as a terminal electron acceptor, promoting anaerobic respiration of SOM and increasing pH. Current large-scale models of Arctic C cycling do not include Fe cycling or pH interactions. Here, a geochemical reaction model was developed by coupling Fe redox reactions and C cycling to simulate SOM decomposition, Fe(III) reduction, pH dynamics, and greenhouse gas production in permafrost soils subject to dynamic hydrology. In this study we parameterized the model using measured CO 2 and CH 4 fluxes as well as changes in pH, Fe(II), and dissolved organic C concentrations from oxic and anoxic incubations of permafrost soils from polygonal permafrost sites in northern Alaska, United States. In simulations of repeated oxic-anoxic cycles, Fe(III) reduction during anoxic periods enhanced CO 2 production, while the net effect of Fe(III) reduction on cumulative CH 4 fluxes depended on substrate C availability. With lower substrate availability, Fe(III) reduction decreased total CH 4 production by further limiting available substrate. With higher substrate availability, Fe(III) reduction enhanced CH 4 production by increasing pH. Our results suggest that interactions among Fe-redox reactions, pH and methanogenesis are important factors in predicting CH 4 and CO 2 production as well as SOM decomposition rates in Fe-rich, frequently waterlogged Arctic soils.

54 ENVIRONMENTAL SCIENCES↗

Confronting the water potential information gap

Water potential directly controls the function of leaves, roots and microbes, and gradients in water potential drive water flows throughout the soil–plant–atmosphere continuum. Notwithstanding its clear relevance for many ecosystem processes, soil water potential is rarely measured in situ, and plant water potential observations are generally discrete, sparse, and not yet aggregated into accessible databases. These gaps limit our conceptual understanding of biophysical responses to moisture stress and inject large uncertainty into hydrologic and land-surface models. Here, we outline the conceptual and predictive gains that could be made with more continuous and discoverable observations of water potential in soils and plants. We discuss improvements to sensor technologies that facilitate in situ characterization of water potential, as well as strategies for building new networks that aggregate water potential data across sites. Here, we end by highlighting novel opportunities for linking more representative site-level observations of water potential to remotely sensed proxies. Together, these considerations offer a road map for clearer links between ecohydrological processes and the water potential gradients that have the ‘potential’ to substantially reduce conceptual and modelling uncertainties.

58 GEOSCIENCES↗

A new bioenergy model that simulates the impacts of plant‐microbial interactions, soil carbon protection, and mechanistic tillage on soil carbon cycling

Abstract Advancing our predictive understanding of bioenergy systems is critical to design decision tools that can inform which feedstock to plant, where to plant it, and how to manage its production to provide both energy and ecosystem carbon (C) benefits. Here, we lay the foundation for that advancement by integrating recent developments in the science of belowground processes in shaping the C cycle into a new bioenergy model, FUN‐BioCROP (Fixation and Uptake of Nitrogen‐Bioenergy Carbon, Rhizosphere, Organisms, and Protection). We show that FUN‐BioCROP can approximate the historical trajectory of soil C dynamics as natural ecosystems were successively converted into intensive agriculture and bioenergy systems. This ability relies in part on a novel tillage representation that mechanistically models tillage as a process that increases microbial access to C. Importantly, the impacts of tillage and feedstock choice also influence FUN‐BioCROP simulations of warming responses with no‐till perennial feedstocks, miscanthus, and switchgrass, having more C that is unprotected and susceptible to warming than tilled annual feedstocks like corn–corn–soybean. However, this susceptibility to warming is balanced by a greater potential for increases in belowground C allocation to enhance soil C stocks in perennial systems. Collectively, our model results highlight the importance of belowground processes in evaluating the ecosystem C benefits of bioenergy production.

09 BIOMASS FUELS↗

Leveraging observed soil heterotrophic respiration fluxes as a novel constraint on global-scale models

Microbially-explicit models may improve understanding and projections of carbon dynamics in response to future climate change, but their fidelity in simulating global-scale soil heterotrophic respiration (RH), a stringent test for soil biogeochemical models, has never been evaluated. We used statistical global RH products, as well as 7,821 daily site-scale RH measurements, to evaluate the spatio-temporal performance of one first-order decay model (CASA-CNP) and two microbially-explicit biogeochemical models (CORPSE and MIMICS) that were forced by two different climate datasets. CORPSE and MIMICS did not provide any measurable performance improvement; instead, the models were highly sensitive to the meteorological input data used to drive them. Spatial RH variability was generally well simulated except in the northern middle latitudes (~50°N) and arid regions; models captured the seasonal variability of RH well, but showed more divergence in tropic and arctic regions. Our results demonstrate that the next generation of biogeochemical models shows promise, but also needs to be improved for realistic spatio-temporal variability of RH. Finally, we emphasize the importance of net primary production, soil moisture, and soil temperature inputs, and that jointly evaluating soil models for their spatial (global scale) and temporal (site scale) performance provides crucial benchmarks for improving biogeochemical models.

Jian, Jinshi↗

Mycorrhizal Distributions Impact Global Patterns of Carbon and Nutrient Cycling

Most tree species predominantly associate with a single type of mycorrhizal fungi, which can differentially affect plant nutrient acquisition and biogeochemical cycling. Uncertainties in mycorrhizal distributions are non-trivial, and current estimates disagree in up to 50% over 40% of the land area, including tropical forests. Remote sensing capabilities for mycorrhizal detection show promise for refining these estimates further. Here, we address for the first time the impact of mycorrhizal distributions on global carbon and nutrient cycling. Using the state-of-the-art carbon-nitrogen economics within the Community Land Model version 5, we found Net Primary Productivity (NPP) increased throughout the 21st century by 20%; however, as soil nitrogen has progressively become limiting, the costs to NPP for nitrogen acquisition—that is, to mycorrhizae—have increased at a faster rate by 60%. This suggests that nutrient acquisition will increasingly demand a higher portion of assimilated carbon to support the same productivity.

54 ENVIRONMENTAL SCIENCES↗

SoDaH: the SOils DAta Harmonization database, an open-source synthesis of soil data from research networks, version 1.0

Data collected from research networks present opportunities to test theories and develop models about factors responsible for the long-term persistence and vulnerability of soil organic matter (SOM). Synthesizing datasets collected by different research networks presents opportunities to expand the ecological gradients and scientific breadth of information available for inquiry. Synthesizing these data is challenging, especially considering the legacy of soil data that have already been collected and an expansion of new network science initiatives. To facilitate this effort, here we present the SOils DAta Harmonization database (SoDaH; https://lter.github.io/som-website, last access: 22 December 2020), a flexible database designed to harmonize diverse SOM datasets from multiple research networks. SoDaH is built on several network science efforts in the United States, but the tools built for SoDaH aim to provide an open-access resource to facilitate synthesis of soil carbon data. Moreover, SoDaH allows for individual locations to contribute results from experimental manipulations, repeated measurements from long-term studies, and local- to regional-scale gradients across ecosystems or landscapes. Finally, we also provide data visualization and analysis tools that can be used to query and analyze the aggregated database. The SoDaH v1.0 dataset is archived and available at https://doi.org/10.6073/pasta/9733f6b6d2ffd12bf126dc36a763e0b4 (Wieder et al., 2020).

54 ENVIRONMENTAL SCIENCES↗

Editorial: Forest Rhizosphere Interactions: Cascading Consequences for Ecosystem-Level Carbon and Nutrient Cycling

Mycorrhizal fungimine organicmatter and subsequently transfer plant-available nutrients to roots. While there is disagreement considering the role of mycorrhizal fungi in ecosystem carbon (C) and nutrient cycling relative to factors such as climate (Cotton, 2018; Bennett and Classen, 2020), soil nutrient availability (Brzostek et al., 2015; Frey, 2019), or host plant community composition (Read and Perez-Moreno, 2003), it is now clear that mycorrhizal fungal symbionts control a substantial proportion of C and nutrient flow belowground (Steidinger et al., 2019). Yet, determining the relative importance of arbuscular mycorrhizal (AM) and ectomycorrhizal (EM) fungi to these processes is still an open field of inquiry in ecosystemecology (Averill et al., 2019a) and earth system modeling (Sulman et al., 2019). At large scales, ecosystems dominated by EM fungi are thought to store more soil C per unit nitrogen (N) (Averill et al., 2014), have slower N cycles (Phillips et al., 2013), and respond more strongly to global change drivers such as N-deposition (Jo et al., 2019). However, it is unclear how generalizable these patterns are across sites, and plant and mycorrhizal taxa. If site or species-dependent relationships exist between mycorrhizal plants, fungi, and soil biogeochemical cycles, forecasting hot spots and hot moments of mycorrhizal fungal activities and their consistency across space and time would be evenmore challenging (Averill et al., 2019b; Kivlin and Hawkes, 2020).

54 ENVIRONMENTAL SCIENCES↗

Integrating Arctic Plant Functional Types in a Land Surface Model Using Above‐ and Belowground Field Observations

Abstract Accurate simulations of high‐latitude ecosystems are critical for confident Earth system model (ESM) projections of carbon cycle feedbacks to global climate change. Land surface model components of ESMs, including the E3SM Land Model (ELM), simulate vegetation growth and ecosystem responses to changing climate and atmospheric CO 2 concentrations by grouping heterogeneous vegetation into like sets of plant functional types (PFTs). Many such models represent high‐latitude vegetation using only two PFTs (shrub and grass), thereby missing the diversity of vegetation growth forms and functional traits in the Arctic. Here, we use field observations of biomass and leaf traits across a gradient of plant communities on the Seward Peninsula in northwest Alaska to replace the original ELM configuration for the first time with nine Arctic‐specific PFTs. The newly developed PFTs include: (1) nonvascular mosses and lichens, (2) deciduous and evergreen shrubs of various height classes, including an alder PFT, (3) graminoids, and (4) forbs. Improvements relative to the original model configuration included greater belowground biomass allocation, persistent fine roots and rhizomes of nonwoody plants, and better representation of variability in total plant biomass across sites with varying plant communities and depth to bedrock. Simulations through 2100 using the RCP8.5 climate scenario and constant PFT fractional areas showed alder‐dominated plant communities gaining more biomass and lichen‐dominated communities gaining less biomass compared to default PFTs. Our results highlight how representing the diversity of arctic vegetation and confronting models with measurements from varied plant communities improves the representation of arctic vegetation in terrestrial ecosystem models.

54 ENVIRONMENTAL SCIENCES↗

A trade-off between plant and soil carbon storage under elevated CO 2

Terrestrial ecosystems remove about 30 per cent of the carbon dioxide (CO 2 ) emitted by human activities each year, yet the persistence of this carbon sink depends partly on how plant biomass and soil organic carbon (SOC) stocks respond to future increases in atmospheric CO 2 . Although plant biomass often increases in elevated CO 2 (eCO 2 ) experiments, SOC has been observed to increase, remain unchanged or even decline. The mechanisms that drive this variation across experiments remain poorly understood, creating uncertainty in climate projections. In this study we synthesized data from 108 eCO 2 experiments and found that the effect of eCO 2 on SOC stocks is best explained by a negative relationship with plant biomass: when plant biomass is strongly stimulated by eCO 2 , SOC storage declines; conversely, when biomass is weakly stimulated, SOC storage increases. This trade-off appears to be related to plant nutrient acquisition, in which plants increase their biomass by mining the soil for nutrients, which decreases SOC storage. We found that, overall, SOC stocks increase with eCO 2 in grasslands (8 ± 2 per cent) but not in forests (0 ± 2 per cent), even though plant biomass in grasslands increase less (9 ± 3 per cent) than in forests (23 ± 2 per cent). Ecosystem models do not reproduce this trade-off, which implies that projections of SOC may need to be revised.

54 ENVIRONMENTAL SCIENCES↗