Spatial biases of information influence global estimates of soil respiration: How can we improve global predictions?

Soil respiration (Rs), the efflux of CO2 from soils to the atmosphere, is a major component of the terrestrial carbon cycle, but is poorly constrained from regional to global scales. The global soil respiration database (SRDB) is a compilation of in situ Rs observations from around the globe that has been consistently updated with new measurements over the past decade. It is unclear whether the addition of data to new versions has produced better-constrained global Rs estimates. We compared two versions of the SRDB (v3.0 n = 5173 and v5.0 n = 10,366) to determine how additional data influenced global Rs annual sum, spatial patterns and associated uncertainty (1 km spatial resolution) using a machine learning approach. A quantile regression forest model parameterized using SRDBv3 yielded a global Rs sum of 88.6 Pg C year-1 , and associated uncertainty of 29.9 (mean absolute error) and 57.9 (standard deviation) Pg C year-1 , whereas parameterization using SRDBv5 yielded 96.5 Pg C year-1 and associated uncertainty of 30.2 (mean average error) and 73.4 (standard deviation) Pg C year-1 . Empirically estimated global heterotrophic respiration (Rh) from v3 and v5 were 49.9-50.2 (mean 50.1) and 53.3-53.5 (mean 53.4) Pg C year-1 , respectively. SRDBv5's inclusion of new data from underrepresented regions (e.g., Asia, Africa, South America) resulted in overall higher model uncertainty. The largest differences between models parameterized with different SRDVB versions were in arid/semi-arid regions. The SRDBv5 is still biased toward northern latitudes and temperate zones, so we tested an optimized global distribution of Rs measurements, which resulted in a global sum of 96.4 ± 21.4 Pg C year-1 with an overall lower model uncertainty. These results support current global estimates of Rs but highlight spatial biases that influence model parameterization and interpretation and provide insights for design of environmental networks to improve global-scale Rs estimates.

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 7821 daily site-scale RH measurements, to evaluate the spatiotemporal performance of one first-order decay model (CASA-CNP) and two microbially explicit biogeochemical models (CORPSE and MIMICS) that were forced by two different input datasets. CORPSE and MIMICS did not provide any measurable performance improvement; instead, the models were highly sensitive to the input data used to drive them. Spatial variability in RH fluxes 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 spatiotemporal 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.

COSORE: A community database for continuous soil respiration and other soil-atmosphere greenhouse gas flux data.

Globally, soils store two to three times as much carbon as currently resides in the atmosphere, and it is critical to understand how soil greenhouse gas (GHG) emissions and uptake will respond to ongoing climate change. In particular, the soil-to-atmosphere CO2 flux, commonly though imprecisely termed soil respiration (RS ), is one of the largest carbon fluxes in the Earth system. An increasing number of high-frequency RS measurements (typically, from an automated system with hourly sampling) have been made over the last two decades; an increasing number of methane measurements are being made with such systems as well. Such high frequency data are an invaluable resource for understanding GHG fluxes, but lack a central database or repository. Here we describe the lightweight, open-source COSORE (COntinuous SOil REspiration) database and software, that focuses on automated, continuous and long-term GHG flux datasets, and is intended to serve as a community resource for earth sciences, climate change syntheses and model evaluation. Contributed datasets are mapped to a single, consistent standard, with metadata on contributors, geographic location, measurement conditions and ancillary data. The design emphasizes the importance of reproducibility, scientific transparency and open access to data. While being oriented towards continuously measured RS , the database design accommodates other soil-atmosphere measurements (e.g. ecosystem respiration, chamber-measured net ecosystem exchange, methane fluxes) as well as experimental treatments (heterotrophic only, etc.). We give brief examples of the types of analyses possible using this new community resource and describe its accompanying R software package.

Constraining estimates of global soil respiration by quantifying sources of variability.

Quantifying global soil respiration (RSG ) and its response to temperature change are critical for predicting the turnover of terrestrial carbon stocks and their feedbacks to climate change. Currently, estimates of RSG range from 68 to 98 Pg C year-1 , causing considerable uncertainty in the global carbon budget. We argue the source of this variability lies in the upscaling assumptions regarding the model format, data timescales, and precipitation component. To quantify the variability and constrain RSG , we developed RSG models using Random Forest and exponential models, and used different timescales (daily, monthly, and annual) of soil respiration (RS ) and climate data to predict RSG . From the resulting RSG estimates (range = 66.62-100.72 Pg), we calculated variability associated with each assumption. Among model formats, using monthly RS data rather than annual data decreased RSG by 7.43-9.46 Pg; however, RSG calculated from daily RS data was only 1.83 Pg lower than the RSG from monthly data. Using mean annual precipitation and temperature data instead of monthly data caused +4.84 and -4.36 Pg C differences, respectively. If the timescale of RS data is constant, RSG estimated by the first-order exponential (93.2 Pg) was greater than the Random Forest (78.76 Pg) or second-order exponential (76.18 Pg) estimates. These results highlight the importance of variation at subannual timescales for upscaling to RSG. The results indicated RSG is lower than in recent papers and the current benchmark for land models (98 Pg C year-1 ), and thus may change the predicted rates of terrestrial carbon turnover and the carbon to climate feedback as global temperatures rise.

Vectorized dataset of silted land formed by check dams on the Chinese Loess Plateau.

Check dams on the Chinese Loess Plateau (CLP) have captured billions of tons of eroded sediment, substantially reducing sediment load in the Yellow River. However, uncertainties persist regarding the precise sediment capture and the role of these dams in Yellow River flow and sediment dynamics due to the lack of available spatial distribution datasets. We produced the first vectorized dataset of silted land formed by check dams on the CLP, combining high-resolution and easily accessible Google Earth images with object-based classification methods. The accuracy of the dataset was verified by 1947 collected test samples, and the producer's accuracy and user's accuracy of the dam lands were 88.9% and 99.5%, respectively. Our dataset not only provides fundamental information for accurately assessing the ecosystem service functions of check dams, but also helps to interpret current changes in sediment delivery of the Yellow River and plan future soil and water conservation projects.

Warming promotes accumulation of microbial- and plant-derived carbon in terrestrial ecosystems.

The impact of global warming on soil carbon pools has been extensively investigated, however, there is still a lack of understanding regarding the specific response of microbial- and plant-derived carbon to warming. To address this knowledge gap, we conducted a comprehensive meta-analysis of 142 studies and evaluated 986 observations comparisons of different carbon source responses to warming. Our results revealed several key insights. Firstly, climate warming resulted in an average increase of 5.46 % in the terrestrial soil carbon pool. Specifically, microbial-derived carbon showed an average increase of 6.32 %, while plant-derived carbon exhibited an average increase of 3.70 %. Secondly, while warming duration and magnitude do not significantly affect the response of microbial-derived carbon to warming, they did impact the response of plant-derived carbon. Lastly, we observed that the response of different carbon sources to warming was affected by the specific environmental backgrounds:ecosystem and climatic zone types affect the response of warming to microbial-derived carbon, while differences in climatic region affect response of warming to plant-derived carbon. The variations in the response of different soil carbon sources to warming can be attributed to the nature of the carbon source themselves, as well as the complex transformations that occur between them through microbial metabolic processes and their interactions with soil mineral particles. We suggest that interactions at the soil-plant-microbe interface should be considered more carefully, and the response of ecosystems to warming should be observed from the perspective of soil organic carbon sources, so as to better understand the response of terrestrial ecosystems carbon cycle to global warming.

Global cross-biome patterns of soil respiration responses to individual and interactive effects of nitrogen addition, altered precipitation, and warming.

Anthropogenic activities have increased atmospheric N, precipitation, and temperature events in terrestrial ecosystems globally, with N deposition increasing by 3- to 5-fold during the previous century. Despite decades of scientific research, no consensus has been achieved on the impact of climate conditions on soil respiration (Rs). Here, we reconstructed 110 published studies across diverse biomes, magnitudes, and driving variables to evaluate how Rs responds to N addition, altered precipitation (both enhanced and reduced precipitation or precipitation changes), and warming. Our findings show that N addition significantly increased Rs by 44 % in forests and decreased it by 19 % and 26 % in croplands and grasslands, respectively (P < 0.05). In forests and croplands, altered precipitation significantly increased Rs by 51 % and 17 % (all, P < 0.05), respectively, whereas impacts on grassland were insignificant. In comparison, warming stimulated Rs by 62 % in forests but inhibited it by 10 % in croplands (all, P < 0.05), whereas impacts on grassland were again insignificant. In addition, across all biomes, the responses of Rs to altered precipitation and warming followed a Gaussian response, increasing up to a threshold of 1800 mm and 25 °C, respectively, above which respiration rates decreased with further increases in precipitation and temperature. Our work suggests that the dual interaction of warming × altered precipitation promotes belowground CO2 emission, thus enhancing global warming. In general, the interactive effect of N addition × altered precipitation decreases Rs. Soil moisture was identified as a primary driver of Rs. Given these findings, we recommend future research on warming vs. changed precipitation to better forecast and understand the interaction between Rs and climate change.

Controls and variability of soil respiration temperature sensitivity across China.

Understanding the temperature sensitivity (Q10) of soil respiration is critical for benchmarking the potential intensity of regional and global terrestrial soil carbon fluxes-climate feedbacks. Although field observations have demonstrated the strong spatial heterogeneity of Q10, a significant knowledge gap still exists regarding to the factors driving spatial and temporal variabilities of Q10 at regional scales. Therefore, we used a machine learning approach to predict Q10 from 1994 to 2016 with a spatial resolution of 1 km across China from 515 field observations at 5 cm soil depth using climate, soil and vegetation variables. Predicted Q10 varied from 1.54 to 4.17, with an area-weighted average of 2.52. There was no significant temporal trend for Q10 (p = 0.32), but annual vegetation production (indicated by normalized difference vegetation index, NDVI) was positively correlated to it (p < 0.01). Spatially, soil organic carbon (SOC) was the most important driving factor in 62 % of the land area across China, and varied greatly, demonstrating soil controls on the spatial pattern of Q10. These findings highlighted different environmental controls on the spatial and temporal pattern of soil respiration Q10, which should be considered to improve global biogeochemical models used to predict the spatial and temporal patterns of soil carbon fluxes to ongoing climate change.

Historically inconsistent productivity and respiration fluxes in the global terrestrial carbon cycle.

The terrestrial carbon cycle is a major source of uncertainty in climate projections. Its dominant fluxes, gross primary productivity (GPP), and respiration (in particular soil respiration, RS), are typically estimated from independent satellite-driven models and upscaled in situ measurements, respectively. We combine carbon-cycle flux estimates and partitioning coefficients to show that historical estimates of global GPP and RS are irreconcilable. When we estimate GPP based on RS measurements and some assumptions about RS:GPP ratios, we found the resulted global GPP values (bootstrap mean [Formula: see text] Pg C yr-1) are significantly higher than most GPP estimates reported in the literature ([Formula: see text] Pg C yr-1). Similarly, historical GPP estimates imply a soil respiration flux (RsGPP, bootstrap mean of [Formula: see text] Pg C yr-1) statistically inconsistent with most published RS values ([Formula: see text] Pg C yr-1), although recent, higher, GPP estimates are narrowing this gap. Furthermore, global RS:GPP ratios are inconsistent with spatial averages of this ratio calculated from individual sites as well as CMIP6 model results. This discrepancy has implications for our understanding of carbon turnover times and the terrestrial sensitivity to climate change. Future efforts should reconcile the discrepancies associated with calculations for GPP and Rs to improve estimates of the global carbon budget.

The interplay of labile organic carbon, enzyme activities and microbial communities of two forest soils across seasons.

Soil labile organic carbon (LOC) responds rapidly to environmental changes and plays an important role in carbon cycle. In this study, the seasonal fluctuations in LOC, the activities of carbon-cycle related enzymes, and the bacterial and fungal communities were analyzed for soils collected from two forests, namely Betula albosinensis (Ba) and Picea asperata Mast. (Pa), in the Qinling Mountains of China. Results revealed that the seasonal average contents of microbial biomass carbon (MBC), easily oxidized organic carbon (EOC), and dissolved organic carbon (DOC) of Pa forest soil were 13.5%, 30.0% and 15.7% less than those in Ba soil. The seasonal average enzyme activities of ß-1,4-glucosidase (ßG), and ß-1,4-xylosidase (ßX) of Ba forest soils were 30.0% and 32.3% higher than those of Pa soil while the enzyme activity of cellobiohydrolase (CBH) was 19.7% lower. Furthermore, the relative abundance of Acidobacteria was significantly higher in summer than in winter, whereas the relative abundance of Bacteroidetes was higher in winter. Regarding the fungal communities, the relative abundance of Basidiomycota was lowest in winter, whereas Ascomycota predominated in the same season. In addition, the soil LOC was significantly positively correlated with the CBH, ßG and ßX activities. Changes in LOC were significantly correlated with Acidobacteria, Bacteroidetes and Basidiomycota. We conclude that the seasonal fluctuations in forest soil LOC fractions relied on carbon cycle-associated enzymatic activities and microorganisms, which in turn were affected by climatic conditions.

A decreasing carbon allocation to belowground autotrophic respiration in global forest ecosystems.

Belowground autotrophic respiration (RAsoil) depends on carbohydrates from photosynthesis flowing to roots and rhizospheres, and is one of the most important but least understood components in forest carbon cycling. Carbon allocation plays an important role in forest carbon cycling and reflects forest adaptation to changing environmental conditions. However, carbon allocation to RAsoil has not been fully examined at the global scale. To fill this knowledge gap, we first used a Random Forest algorithm to predict the spatio-temporal patterns of RAsoil from 1981 to 2017 based on the most updated Global Soil Respiration Database (v5) with global environmental variables; calculated carbon allocation from photosynthesis to RAsoil (CAB) as a fraction of gross primary production; and assessed its temporal and spatial patterns in global forest ecosystems. Globally, mean RAsoil from forests was 8.9 ± 0.08 Pg C yr-1 (mean ± standard deviation) from 1981 to 2017 and increased significantly at a rate of 0.006 Pg C yr-2, paralleling broader soil respiration changes and suggesting increasing carbon respired by roots. Mean CAB was 0.243 ± 0.016 and decreased over time. The temporal trend of CAB varied greatly in space, reflecting uneven responses of CAB to environmental changes. Combined with carbon use efficiency, our CAB results offer a completely independent approach to quantify global aboveground autotropic respiration spatially and temporally, and could provide crucial insights into carbon flux partitioning and global carbon cycling under climate change.

Quantifying cover crop effects on soil health and productivity.

The dataset presented here supports the research paper entitled "A calculator to quantify cover crop effects on soil health and productivity". Soil health (sometimes used synonymously with soil quality) is a concept that describes soil as a living system to sustain plants, animals, and human. Soil physical, chemical, and biological properties, along with their interactions, are required to quantify soil health. The use of cover crops in agricultural rotations may enhance soil health, yet there has been little progress in understanding how external factors such as climate, soil type, and agronomic practices affect soil and cash crop responses. In response, this dataset compiles measurements from 281 studies and provides an analysis of field-measured changes in 38 soil health indicators due to cover crop usage. Environmental and background indicators were also compiled to assess how climatic and management practices affect soil and cash crop responses to cover crops, with specific categories including climate type (tropical, arid, temperate, and continental), soil texture (coarse, medium, and fine), cover crop type (legume, grass, multi-species mixture, and other), and cash crop type (corn, soybean, wheat, vegetable, corn-soybean rotation, corn-soybean-wheat rotation, and other). An unbalanced analysis of variation was used to determine the hierarchy of most to least important factors that affected responsiveness of each soil health indicator. Based on the hierarchy structure, a soil health calculator was then developed to quantify the response of 13 parameters - erosion, runoff, weed suppression, soil aggregate stability, leaching, infiltration, microbial biomass carbon, soil bulk density, soil organic carbon, soil nitrogen, microbial biomass nitrogen, cash crop yield, and saturated hydraulic conductivity - to cover crops. The presented data in the calculator report the mean change in parameter values based on all combinations of climate, soil texture, cover crop type, and cash crop type.

A database for global soil health assessment.

Field studies have been performed for decades to analyze effects of different management practices on agricultural soils and crop yields, but these data have never been combined together in a way that can inform current and future cropland management. Here, we collected, extracted, and integrated a database of soil health measurements conducted in the field from sites across the globe. The database, named SoilHealthDB, currently focuses on four main conservation management methods: cover crops, no-tillage, agro-forestry systems, and organic farming. These studies represent 354 geographic sites (i.e., locations with unique latitudes and longitudes) in 42 countries around the world. The SoilHealthDB includes 42 soil health indicators and 46 background indicators that describe factors such as climate, elevation, and soil type. A primary goal of this effort is to enable the research community to perform comprehensive analyses, e.g., meta-analyses, of soil health changes related to cropland conservation management. The database also provides a common framework for sharing soil health, and the scientific research community is encouraged to contribute their own measurements.

Fecal indicator dynamics at the watershed scale: Variable relationships with land use, season, and water chemistry.

Tracking fecal contamination in surface waters is critical to remediating water quality; however, general and source-specific fecal indicators often provide conflicting results. To understand the spatial and temporal dynamics of multiple fecal indicators and the sources they represent, we measured weekly concentrations of two general fecal indicator bacteria (FIB), a genetic indicator of human-associated Bacteroides (HF183), and surface water chemistry in nine mixed land-use watersheds in southwest Virginia, USA. At the watershed scale, general and source-specific indicators were decoupled, with distinct spatial, temporal, and chemical patterns. Random Forest analysis of individual sample variability identified temperature, watershed, nutrients, and cations as top predictors of indicator concentrations. However, these patterns - and the specific nutrients and cations identified - varied by indicator type. Among watersheds, FIB increased with developed land cover and during the summer months, while HF183 increased during the winter and only in urban watersheds. Nutrients generally related poorly to FIB and HF183, except E. coli, which correlated with total nitrogen. In contrast, all fecal indicators showed strong correlations with cations. FIB were more strongly related to calcium, magnesium, and potassium concentrations, while HF183 was related to sodium. These results suggest that, even at the watershed scale, 1) HF183 detects mainly human fecal contamination, while FIB detect broader ecosystem fecal inputs, and 2) poor correlation between specific and generalist fecal indicators is caused by unique spatial, temporal, and transport dynamics of different fecal sources in watersheds.

An analytical approach to ascertain saturation-excess versus infiltration-excess overland flow in urban and reference landscapes.

Uncontrolled overland flow drives flooding, erosion, and contaminant transport, with the severity of these outcomes often amplified in urban areas. In pervious media such as urban soils, overland flow is initiated via either infiltration-excess (where precipitation rate exceeds infiltration capacity) or saturation-excess (when precipitation volume exceeds soil profile storage) mechanisms. These processes call for different management strategies, making it important for municipalities to discern between them. In this study, we derived a generalized one-dimensional model that distinguishes between infiltration-excess overland flow (IEOF) and saturation-excess overland flow (SEOF) using Green-Ampt infiltration concepts. Next, we applied this model to estimate overland flow generation from pervious areas in 11 U.S. cities. We used rainfall forcing that represented low- and high-intensity events and compared responses among measured urban versus predevelopment reference soil hydraulic properties. The derivation showed that the propensity for IEOF versus SEOF is related to the equivalence between two nondimensional ratios: (a) precipitation rate to depth-weighted hydraulic conductivity and (b) depth of soil profile restrictive layer to soil capillary potential. Across all cities, reference soil profiles were associated with greater IEOF for the high-intensity set of storms, and urbanized soil profiles tended towards production of SEOF during the lower intensity set of storms. Urban soils produced more cumulative overland flow as a fraction of cumulative precipitation than did reference soils, particularly under conditions associated with SEOF. These results will assist cities in identifying the type and extent of interventions needed to manage storm water produced from pervious areas.

Measurement on physical parameters of raindrop energy.

Rainfall erosivity factor (R) is one of the most commonly used factors in soil erosion models. While rainfall energy (E) is the most elementary physical parameter to predict R. Based on comparative analysis of previous soil erosion models and rainfall erosivity factor measuring methods, integrated application of modern photogrammetric techniques, image analytic methods and automatic control theories, this paper provided a new method based on image analytic to calculate the rainfall energy and R factor, which obtains raindrop's volume and velocity by means of modern photogrammetric technique. Results show that this method can improve both efficiency and accuracy of rainfall energy calculation and other rainfall physical parameters measurement.