Description: The terrestrial carbon (C) cycle has been commonly represented by a series of C balance equations to track C influxes into and effluxes out of individual pools in earth system models (ESMs). This representation matches our understanding of C cycle processes well but makes it difficult to track model behaviors. It is also computationally expensive, limiting the ability to conduct comprehensive parametric sensitivity analyses. To overcome these challenges, we have developed a matrix approach, which reorganizes the C balance equations in the original ESM into one matrix equation without changing any modeled C cycle processes and mechanisms. We applied the matrix approach to the Community Land Model (CLM4.5) with vertically resolved biogeochemistry. The matrix equation exactly reproduces litter and soil organic carbon (SOC) dynamics of the standard CLM4.5 across different spatial-temporal scales. The matrix approach enables effective diagnosis of system properties such as C residence time and attribution of global change impacts to relevant processes. We illustrated, for example, the impacts of CO2 fertilization on litter and SOC dynamics can be easily decomposed into the relative contributions from C input, allocation of external C into different C pools, nitrogen regulation, altered soil environmental conditions, and vertical mixing along the soil profile. In addition, the matrix tool can accelerate model spin-up, permit thorough parametric sensitivity tests, enable pool-based data assimilation, and facilitate tracking and benchmarking of model behaviors. Overall, the matrix approach can make a broad range of future modeling activities more efficient and effective.

Description: Ponds and lakes are abundant in Arctic permafrost lowlands. They play an important role in Arctic wetland ecosystems by regulating carbon, water, and energy fluxes and providing freshwater habitats. However, ponds, i.e., waterbodies with surface areas smaller than 1.0 × 10**4 m**2, have not been inventoried on global and regional scales. The Permafrost Region Pond and Lake (PeRL) database presents the results of a circum-Arctic effort to map ponds and lakes from modern (2002-2013) high-resolution aerial and satellite imagery with a resolution of 5 m or better. The database also includes historical imagery from 1948 to 1965 with a resolution of 6 m or better. PeRL includes 69 maps covering a wide range of environmental conditions from tundra to boreal regions and from continuous to discontinuous permafrost zones. Waterbody maps are linked to regional permafrost landscape maps which provide information on permafrost extent, ground ice volume, geology, and lithology. This paper describes waterbody classification and accuracy, and presents statistics of waterbody distribution for each site. Maps of permafrost landscapes in Alaska, Canada, and Russia are used to extrapolate waterbody statistics from the site level to regional landscape units. PeRL presents pond and lake estimates for a total area of 1.4 × 10**6 km**2 across the Arctic, about 17 % of the Arctic lowland ( 〈 300 m a.s.l.) land surface area. PeRL waterbodies with sizes of 1.0 × 10**6 m**2 down to 1.0 × 10**2 m**2 contributed up to 21 % to the total water fraction. Waterbody density ranged from 1.0 × 10 to 9.4 × 10**1/km². Ponds are the dominant waterbody type by number in all landscapes representing 45-99 % of the total waterbody number. The implementation of PeRL size distributions in land surface models will greatly improve the investigation and projection of surface inundation and carbon fluxes in permafrost lowlands.

Description: The Zero Emissions Commitment (ZEC) is the change in global mean temperature expected to occur following the cessation of net CO2 emissions and as such is a critical parameter for calculating the remaining carbon budget. The Zero Emissions Commitment Model Intercomparison Project (ZECMIP) was established to gain a better understanding of the potential magnitude and sign of ZEC, in addition to the processes that underlie this metric. A total of 18 Earth system models of both full and intermediate complexity participated in ZECMIP. All models conducted an experiment where atmospheric CO2 concentration increases exponentially until 1000 PgC has been emitted. Thereafter emissions are set to zero and models are configured to allow free evolution of atmospheric CO2 concentration. Many models conducted additional second-priority simulations with different cumulative emission totals and an alternative idealized emissions pathway with a gradual transition to zero emissions. The inter-model range of ZEC 50 years after emissions cease for the 1000 PgC experiment is −0.36 to 0.29 ∘C, with a model ensemble mean of −0.07 ∘C, median of −0.05 ∘C, and standard deviation of 0.19 ∘C. Models exhibit a wide variety of behaviours after emissions cease, with some models continuing to warm for decades to millennia and others cooling substantially. Analysis shows that both the carbon uptake by the ocean and the terrestrial biosphere are important for counteracting the warming effect from the reduction in ocean heat uptake in the decades after emissions cease. This warming effect is difficult to constrain due to high uncertainty in the efficacy of ocean heat uptake. Overall, the most likely value of ZEC on multi-decadal timescales is close to zero, consistent with previous model experiments and simple theory.

Description: The remaining carbon budget represents the total amount of CO2 that can still be emitted in the future while limiting global warming to a given temperature target. Remaining carbon budget estimates range widely, however, and this uncertainty can be used to either trivialize the most ambitious mitigation targets by characterizing them as impossible, or to argue that there is ample time to allow for a gradual transition to a low-carbon economy. Neither of these extremes is consistent with our best understanding of the policy implications of remaining carbon budgets. Understanding the scientific and socio-economic uncertainties affecting the size of the remaining carbon budgets, as well as the methodological choices and assumptions that underlie their calculation, is essential before applying them as a policy tool. Here we provide recommendations on how to calculate remaining carbon budgets in a traceable and transparent way, and discuss their uncertainties and implications for both international and national climate policies.

Description: The amount of additional future temperature change following a complete cessation of CO2 emissions is a measure of the unrealized warming to which we are committed due to CO2 already emitted to the atmosphere. This "zero emissions commitment" (ZEC) is also an important quantity when estimating the remaining carbon budget - a limit on the total amount of CO2 emissions consistent with limiting global mean temperature at a particular level. In the recent IPCC Special Report on Global Warming of 1.5 degrees C, the carbon budget framework used to calculate the remaining carbon budget for 1.5 degrees C included the assumption that the ZEC due to CO2 emissions is negligible and close to zero. Previous research has shown significant uncertainty even in the sign of the ZEC. To close this knowledge gap, we propose the Zero Emissions Commitment Model Intercomparison Project (ZECMIP), which will quantify the amount of unrealized temperature change that occurs after CO2 emissions cease and investigate the geophysical drivers behind this climate response. Quantitative information on ZEC is a key gap in our knowledge, and one that will not be addressed by currently planned CMIP6 simulations, yet it is crucial for verifying whether carbon budgets need to be adjusted to account for any unrealized temperature change resulting from past CO2 emissions. We request only one top-priority simulation from comprehensive general circulation Earth system models (ESMs) and Earth system models of intermediate complexity (EMICs) - a branch from the 1% CO2 run with CO2 emissions set to zero at the point of 1000 PgC of total CO2 emissions in the simulation - with the possibility for additional simulations, if resources allow. ZECMIP is part of CMIP6, under joint sponsorship by C4MIP and CDR-MIP, with associated experiment names to enable data submissions to the Earth System Grid Federation. All data will be published and made freely available.

Description: Increasing atmospheric methane (CH4) concentrations have contributed to approximately 20% of anthropogenic climate change. Despite the importance of CH4 as a greenhouse gas, its atmospheric growth rate and dynamics over the past two decades, which include a stabilization period (1999–2006), followed by renewed growth starting in 2007, remain poorly understood. We provide an updated estimate of CH4 emissions from wetlands, the largest natural global CH4 source, for 2000–2012 using an ensemble of biogeochemical models constrained with remote sensing surface inundation and inventory-based wetland area data. Between 2000–2012, boreal wetland CH4 emissions increased by 1.2 Tg yr−1 (−0.2–3.5 Tg yr−1), tropical emissions decreased by 0.9 Tg yr−1 (−3.2−1.1 Tg yr−1), yet globally, emissions remained unchanged at 184 ± 22 Tg yr−1. Changing air temperature was responsible for increasing high-latitude emissions whereas declines in low-latitude wetland area decreased tropical emissions; both dynamics are consistent with features of predicted centennial-scale climate change impacts on wetland CH4 emissions. Despite uncertainties in wetland area mapping, our study shows that global wetland CH4 emissions have not contributed significantly to the period of renewed atmospheric CH4 growth, and is consistent with findings from studies that indicate some combination of increasing fossil fuel and agriculture-related CH4 emissions, and a decrease in the atmospheric oxidative sink.

Description: Soils contain a large and dynamic fraction of global terrestrial carbon stocks. The distribution of soil carbon (SC) with depth varies among ecosystems and land uses and is an important factor in calculating SC stocks and their vulnerabilities. Systematic analysis of SC depth distributions across databases of SC profiles has been challenging due to the heterogeneity of soil profile measurements, which vary in depth sampling. Here, we fit over 40,000 SC depth profiles to an exponential decline relationship with depth to determine SC concentration at the top of the mineral soil, minimum SC concentration at depth, and the characteristic “length” of SC concentration decline with depth. Fitting these parameters allowed profile characteristics to be analyzed across a large and heterogeneous dataset. We then assessed the differences in these depth parameters across soil orders and land cover types and between soil profiles with or without a history of tillage, as represented by the presence of an Ap horizon. We found that historically tilled soils had more gradual decreases of SC with depth (greater e-folding depth or Z∗), deeper SC profiles, lower SC concentrations at the top of the mineral soil, and lower total SC stocks integrated to 30 cm. The large database of profiles allowed these results to be confirmed across different land cover types and spatial areas within the Continental United States, providing robust evidence for systematic impacts of historical tillage on SC stocks and depth distributions.

Description: The permafrost zone is expected to be a substantial carbon source to the atmosphere, yet large-scale models currently only simulate gradual changes in seasonally thawed soil. Abrupt thaw will probably occur in 〈20% of the permafrost zone but could affect half of permafrost carbon through collapsing ground, rapid erosion and landslides. Here, we synthesize the best available information and develop inventory models to simulate abrupt thaw impacts on permafrost carbon balance. Emissions across 2.5 million km2 of abrupt thaw could provide a similar climate feedback as gradual thaw emissions from the entire 18 million km2 permafrost region under the warming projection of Representative Concentration Pathway 8.5. While models forecast that gradual thaw may lead to net ecosystem carbon uptake under projections of Representative Concentration Pathway 4.5, abrupt thaw emissions are likely to offset this potential carbon sink. Active hillslope erosional features will occupy 3% of abrupt thaw terrain by 2300 but emit one-third of abrupt thaw carbon losses. Thaw lakes and wetlands are methane hot spots but their carbon release is partially offset by slowly regrowing vegetation. After considering abrupt thaw stabilization, lake drainage and soil carbon uptake by vegetation regrowth, we conclude that models considering only gradual permafrost thaw are substantially underestimating carbon emissions from thawing permafrost.