Abstract. Understanding and quantifying the global methane (CH4) budget is important for assessing realistic pathways to mitigate climate change. Atmospheric emissions and concentrations of CH4 are continuing to increase, making CH4 the second most important human-influenced greenhouse gas in terms of climate forcing, after carbon dioxide (CO2). Assessing the relative importance of CH4 in comparison to CO2 is complicated by its shorter atmospheric lifetime, stronger warming potential, and atmospheric growth rate variations over the past decade, the causes of which are still debated. Two major difficulties in reducing uncertainties arise from the variety of geographically overlapping CH4 sources and from the destruction of CH4 by short-lived hydroxyl radicals (OH). To address these difficulties, we have established a consortium of multi-disciplinary scientists under the umbrella of the Global Carbon Project to synthesize and stimulate new research aimed at improving and regularly updating the global methane budget. Following Saunois et al. (2016), we present here the second version of the living review paper dedicated to the decadal methane budget, integrating results of top-down studies (atmospheric observations within an atmospheric inverse-modelling framework) and bottom-up estimates (including process-based models for estimating land surface emissions and atmospheric chemistry, inventories of anthropogenic emissions, and data-driven extrapolations). For the 2008–2017 decade, global methane emissions are estimated by atmospheric inversions (top-down approach) to be 572 Tg CH4 yr−1 (range 538–593, corresponding to the minimum and maximum estimates of the ensemble), of which 357 Tg CH4 yr−1 or ~ 60 % are attributed to anthropogenic sources (range 50–65 %). This total emission is 27 Tg CH4 yr−1 larger than the value estimated for the period 2000–2009 and 24 Tg CH4 yr−1 larger than the one reported in the previous budget for the period 2003–2012 (Saunois et al. 2016). Since 2012, global CH4 emissions have been tracking the carbon intensive scenarios developed by the Intergovernmental Panel on Climate Change (Gidden et al., 2019). Bottom-up methods suggest larger global emissions (737 Tg CH4 yr−1, range 583–880) than top-down inversion methods, mostly because of larger estimated natural emissions from sources such as natural wetlands, other inland water systems, and geological sources. However the strength of the atmospheric constraints on the top-down budget, suggest that these bottom-up emissions are overestimated. The latitudinal distribution of atmospheric-based emissions indicates a predominance of tropical emissions (~ 65 % of the global budget,

Abstract. Terrestrial biosphere models (TBMs) have become an integral tool for extrapolating local observations and understanding of land-atmosphere carbon exchange to larger regions. The North American Carbon Program (NACP) Multi-scale synthesis and Terrestrial Model Intercomparison Project (MsTMIP) is a formal model intercomparison and evaluation effort focused on improving the diagnosis and attribution of carbon exchange at regional and global scales. MsTMIP builds upon current and past synthesis activities, and has a unique framework designed to isolate, interpret, and inform understanding of how model structural differences impact estimates of carbon uptake and release. Here we provide an overview of the MsTMIP effort and describe how the MsTMIP experimental design enables the assessment and quantification of TBM structural uncertainty. Model structure refers to the types of processes considered (e.g. nutrient cycling, disturbance, lateral transport of carbon), and how these processes are represented (e.g. photosynthetic formulation, temperature sensitivity, respiration) in the models. By prescribing a common experimental protocol with standard spin-up procedures and driver data sets, we isolate any biases and variability in TBM estimates of regional and global carbon budgets resulting from differences in the models themselves (i.e. model structure) and model-specific parameter values. An initial intercomparison of model structural differences is represented using hierarchical cluster diagrams (a.k.a. dendrograms), which highlight similarities and differences in how models account for carbon cycle, vegetation, energy, and nitrogen cycle dynamics. We show that, despite the standardized protocol used to derive initial conditions, models show a high degree of variation for GPP, total living biomass, and total soil carbon, underscoring the influence of differences in model structure and parameterization on model estimates.

Clay-size minerals play important roles in terrestrial biogeochemistry and atmospheric physics, but their data have been only partially compiled at global scale. We present a global dataset of clay-size minerals in the topsoil and subsoil at different spatial resolutions. The data of soil clay and its mineralogical composition were gathered through a literature survey and aggregated by soil orders of the Soil Taxonomy for each of the ten groups: gibbsite, kaolinite, illite/mica, smectite, vermiculite, chlorite, iron oxide, quartz, non-crystalline, and others. Using a global soil map, a global dataset of soil clay-size mineral distribution was developed at resolutions of 2' to 2° grid cells. The data uncertainty associated with data variability and assumption was evaluated using a Monte Carlo method, and validity of the clay-size mineral distribution obtained in this study was examined by comparing with other datasets. The global soil clay data offer spatially explicit studies on terrestrial biogeochemical cycles, dust emission to the atmosphere, and other interdisciplinary earth sciences.

Soil organic carbon (SOC) is the largest carbon pool in terrestrial ecosystems and may play a key role in biospheric feedbacks with elevated atmospheric carbon dioxide (CO2) in a warmer future world.We examined the simulation results of seven terrestrial biome models when forced with climate projections from four representative-concentrationpathways (RCPs)-based atmospheric concentration scenarios. The goal was to specify calculated uncertainty in global SOC stock projections from global and regional perspectives and give insight to the improvement of SOC-relevant processes in biome models. SOC stocks among the biome models varied from 1090 to 2650 PgC even in historical periods (ca. 2000). In a higher forcing scenario (i.e., RCP8.5), inconsistent estimates of impact on the total SOC (2099–2000) were obtained from different biome model simulations, ranging from a net sink of 347 Pg C to a net source of 122 Pg C. In all models, the increasing atmospheric CO2 concentration in the RCP8.5 scenario considerably contributed to carbon accumulation in SOC. However, magnitudes varied from 93 to 264 PgC by the end of the 21st century across biome models. Using the time-series data of total global SOC simulated by each biome model, we analyzed the sensitivity of the global SOC stock to global mean temperature and global precipitation anomalies (1T and 1P respectively) in each biome model using a state-space model. This analysis suggests that 1T explained global SOC stock changes in most models with a resolution of 1–2 C, and the magnitude of global SOC decomposition from a 2 C rise ranged from almost 0 to 3.53 Pg Cyr−1 among the biome models. However, 1P had a negligible impact on change in the global SOC changes. Spatial heterogeneity was evident and inconsistent among the biome models, especially in boreal to arctic regions. Our study reveals considerable climate uncertainty in SOC decomposition responses to climate and CO2 change among biome models. Further research is required to improve our ability to estimate biospheric feedbacks through both SOC-relevant and vegetation-relevant processes

Future changes in runoff can have important implications for water resources and flooding. In this study, runoff projections from ISI-MIP (Inter-sectoral Impact Model Intercomparison Project) simulations forced with HadGEM2-ES bias-corrected climate data under the Representative Concentration Pathway 8.5 have been analysed for differences between impact models. Projections of change from a baseline period (1981–2010) to the future (2070–2099) from 12 impacts models which contributed to the hydrological and biomes sectors of ISI-MIP were studied. The biome models differed from the hydrological models by the inclusion of CO2 impacts and most also included a dynamic vegetation distribution. The biome and hydrological models agreed on the sign of runoff change for most regions of the world. However, in West Africa, the hydrological models projected drying, and the biome models a moistening. The biome models tended to produce larger increases and smaller decreases in regionally averaged runoff than the hydrological models, although there is large inter-model spread. The timing of runoff change was similar, but there were differences in magnitude, particularly at peak runoff. The impact of vegetation distribution change was much smaller than the projected change over time, while elevated CO2 had an effect as large as the magnitude of change over time projected by some models in some regions. The effect of CO2 on runoff was not consistent across the models, with two models showing increases and two decreases. There was also more spread in projections from the runs with elevated CO2 than with constant CO2. The biome models which gave increased runoff from elevated CO2 were also those which differed most from the hydrological models. Spatially, regions with most difference between model types tended to be projected to have most effect from elevated CO2, and seasonal differences were also similar, so elevated CO2 can partly explain the differences between hydrological and biome model runoff change projections. Therefore, this shows that a range of impact models should be considered to give the full range of uncertainty in impacts studies