Soil Respiration Phenology Improves Modeled Phase of Terrestrial net Ecosystem Exchange in Northern Hemisphere

The Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2) is the latest atmospheric reanalysis of the modern satellite era produced by NASA’s Global Modeling and Assimilation Office (GMAO). MERRA-2 assimilates observation types not available to its predecessor, MERRA, and includes updates to the Goddard Earth Observing System (GEOS) model and analysis scheme so as to provide a viable ongoing climate analysis beyond MERRA’s terminus. While addressing known limitations of MERRA, MERRA-2 is also intended to be a development milestone for a future integrated Earth system analysis (IESA) currently under development at GMAO. This paper provides an overview of the MERRA-2 system and various performance metrics. Among the advances in MERRA-2 relevant to IESA are the assimilation of aerosol observations, several improvements to the representation of the stratosphere including ozone, and improved representations of cryospheric processes. Other improvements in the quality of MERRA-2 compared with MERRA include the reduction of some spurious trends and jumps related to changes in the observing system, and reduced biases and imbalances in aspects of the water cycle. Remaining deficiencies are also identified. Production of MERRA-2 began in June 2014 in four processing streams, and converged to a single near-real time stream in mid 2015. MERRA-2 products are accessible online through the NASA Goddard Earth Sciences Data Information Services Center (GES DISC).

Understanding and predicting ecosystem functioning (e.g., carbon and water fluxes) and the role of soils in carbon storage requires an accurate assessment of plant rooting distributions. Here, in a comprehensive literature synthesis, we analyze rooting patterns for terrestrial biomes and compare distributions for various plant functional groups. We compiled a database of 250 root studies, subdividing suitable results into 11 biomes, and fitted the depth coefficient β to the data for each biome (Gale and Grigal 1987). β is a simple numerical index of rooting distribution based on the asymptotic equation Y=1-βd, where d = depth and Y = the proportion of roots from the surface to depth d. High values of β correspond to a greater proportion of roots with depth. Tundra, boreal forest, and temperate grasslands showed the shallowest rooting profiles (β=0.913, 0.943, and 0.943, respectively), with 80-90% of roots in the top 30 cm of soil; deserts and temperate coniferous forests showed the deepest profiles (β=0.975 and 0.976, respectively) and had only 50% of their roots in the upper 30 cm. Standing root biomass varied by over an order of magnitude across biomes, from approximately 0.2 to 5 kg m-2. Tropical evergreen forests had the highest root biomass (5 kg m-2), but other forest biomes and sclerophyllous shrublands were of similar magnitude. Root biomass for croplands, deserts, tundra and grasslands was below 1.5 kg m-2. Root/shoot (R/S) ratios were highest for tundra, grasslands, and cold deserts (ranging from 4 to 7); forest ecosystems and croplands had the lowest R/S ratios (approximately 0.1 to 0.5). Comparing data across biomes for plant functional groups, grasses had 44% of their roots in the top 10 cm of soil. (β=0.952), while shrubs had only 21% in the same depth increment (β=0.978). The rooting distribution of all temperate and tropical trees was β=0.970 with 26% of roots in the top 10 cm and 60% in the top 30 cm. Overall, the globally averaged root distribution for all ecosystems was β=0.966 (r 2=0.89) with approximately 30%, 50%, and 75% of roots in the top 10 cm, 20 cm, and 40 cm, respectively. We discuss the merits and possible shortcomings of our analysis in the context of root biomass and root functioning.

This paper describes the technical development and accuracy assessment of the most recent and improved version of the SoilGrids system at 250m resolution (June 2016 update). SoilGrids provides global predictions for standard numeric soil properties (organic carbon, bulk density, Cation Exchange Capacity (CEC), pH, soil texture fractions and coarse fragments) at seven standard depths (0, 5, 15, 30, 60, 100 and 200 cm), in addition to predictions of depth to bedrock and distribution of soil classes based on the World Reference Base (WRB) and USDA classification systems (ca. 280 raster layers in total). Predictions were based on ca. 150,000 soil profiles used for training and a stack of 158 remote sensing-based soil covariates (primarily derived from MODIS land products, SRTM DEM derivatives, climatic images and global landform and lithology maps), which were used to fit an ensemble of machine learning methods—random forest and gradient boosting and/or multinomial logistic regression—as implemented in the R packages ranger, xgboost, nnet and caret. The results of 10–fold cross-validation show that the ensemble models explain between 56% (coarse fragments) and 83% (pH) of variation with an overall average of 61%. Improvements in the relative accuracy considering the amount of variation explained, in comparison to the previous version of SoilGrids at 1 km spatial resolution, range from 60 to 230%. Improvements can be attributed to: (1) the use of machine learning instead of linear regression, (2) to considerable investments in preparing finer resolution covariate layers and (3) to insertion of additional soil profiles. Further development of SoilGrids could include refinement of methods to incorporate input uncertainties and derivation of posterior probability distributions (per pixel), and further automation of spatial modeling so that soil maps can be generated for potentially hundreds of soil variables. Another area of future research is the development of methods for multiscale merging of SoilGrids predictions with local and/or national gridded soil products (e.g. up to 50 m spatial resolution) so that increasingly more accurate, complete and consistent global soil information can be produced. SoilGrids are available under the Open Data Base License.