Search

Enhanced summer insolation over North Africa induced a monsoon precipitation increase during the mid-Holocene, about 6000 years ago, and led to a widespread expansion of lakes and wetlands in the present-day Sahara. This expansion of lakes and wetlands is documented in paleoenvironmental sediment records, but the spatially sparse and often discontinuous sediment records provide only a fragmentary picture. Previous simulation studies prescribed either a small lake and wetland extent from reconstructions or focused on documented mega-lakes only to investigate their effect on the mid-Holocene climate. In contrast to these studies, we investigate the possible range of mid-Holocene precipitation changes in response to a small-lake extent and a potential maximum lake and wetland extent. Our study shows that during the summer monsoon season, the African rain belt is shifted about 2 to 7∘ farther north in simulations with a maximum lake or wetland extent than in simulations with a small lake extent. This northward extent is caused by a stronger and prolonged monsoon rainfall season over North Africa which is associated with an increased monsoon precipitation over the southern Sahara and an increased precipitation from tropical plumes over the northwestern Sahara. Replacing lakes with vegetated wetlands causes an enhanced precipitation increase, which is likely due to the high surface roughness of the wetlands. A moisture budget analysis reveals that both lakes and wetlands cause a local precipitation increase not only by enhanced evaporation but also by a stronger inland moisture transport and local moisture recycling to the south of Lake Chad and the west Saharan lakes. Analysis of the dynamic response shows that lakes and wetlands cause a circulation response inverse to the one associated with the Saharan heat low. Depending on the latitudinal position of the lakes and wetlands, they predominantly cause a northward shift or a decay of the African Easterly Jet. These results indicate that the latitudinal position of the lakes and wetlands strongly affects the northward extension of the African summer monsoon.

The present study investigates the response of the high-latitude carbon cycle to changes in atmospheric greenhouse gas (GHG) concentrations in idealized climate change scenarios. To this end we use an adapted version of JSBACH – the land surface component of the Max Planck Institute for Meteorology Earth System Model (MPI-ESM) – that accounts for the organic matter stored in the permafrost-affected soils of the high northern latitudes. The model is run under different climate scenarios that assume an increase in GHG concentrations, based on the Shared Socioeconomic Pathway 5 and the Representative Concentration Pathway 8.5, which peaks in the years 2025, 2050, 2075 or 2100, respectively. The peaks are followed by a decrease in atmospheric GHGs that returns the concentrations to the levels at the beginning of the 21st century, reversing the imposed climate change. We show that the soil CO2 emissions exhibit an almost linear dependence on the global mean surface temperatures that are simulated for the different climate scenarios. Here, each degree of warming increases the fluxes by, very roughly, 50 % of their initial value, while each degree of cooling decreases them correspondingly. However, the linear dependence does not mean that the processes governing the soil CO2 emissions are fully reversible on short timescales but rather that two strongly hysteretic factors offset each other – namely the net primary productivity and the availability of formerly frozen soil organic matter. In contrast, the soil methane emissions show a less pronounced increase with rising temperatures, and they are consistently lower after the peak in the GHG concentrations than prior to it. Here, the net fluxes could even become negative, and we find that methane emissions will play only a minor role in the northern high-latitude contribution to global warming, even when considering the high global warming potential of the gas. Finally, we find that at a global mean temperature of roughly 1.75 K (±0.5 K) above pre-industrial levels the high-latitude ecosystem turns from a CO2 sink into a source of atmospheric carbon, with the net fluxes into the atmosphere increasing substantially with rising atmospheric GHG concentrations. This is very different from scenario simulations with the standard version of the MPI-ESM, in which the region continues to take up atmospheric CO2 throughout the entire 21st century, confirming that the omission of permafrost-related processes and the organic matter stored in the frozen soils leads to a fundamental misrepresentation of the carbon dynamics in the Arctic.

Northern peatlands have been a persistent natural carbon sink since the Last Glacial Maximum. The continued growth and expansion of these carbon-rich ecosystems could offset a large portion of anthropogenic carbon emissions before the end of the present interglacial period. Here we used an impeded drainage model and gridded data on the depth to bedrock and the fraction of histosol-type soils to evaluate the limits to the growth of northern peatland carbon stocks. Our results show that the potential carbon stock in northern peatlands could reach a total of 875±125 Pg C before the end of the present interglacial, which could, as a result, remove 330±200 Pg C of carbon from the atmosphere. We argue that northern peatlands, together with the oceans, will potentially play an important role in reducing the atmospheric carbon dioxide concentration over the next 5000 years.

During the last million years, northern Africa has alternated between arid and humid conditions, as recorded by different kinds of climate archives, including fossil pollen, lake sediments, marine sediments and archaeological remains. Variations occur at millennial scale, with dry phases being similar to the current desert state in the region, and with wet phases, known as African Humid Periods (AHPs), characterised by a strong summer monsoon which can carry enough moisture inland to support rivers, lakes and lush vegetation further north than seen today. Recent sediment records from the Mediterranean Sea revealed that the previous five AHPs had different intensities, in relation to rainfall and vegetation extent. Motivated by these findings, our work focuses on explaining what caused such differences in intensity. To this end, we use the CLIMBER-2 climate model to study the AHP response to changes in three drivers of atmospheric dynamics: Earth's orbit variations, atmospheric concentration of CO2 and inland ice extent. Global transient simulations of the last 190,000 years are used in new factorisation analyses, which allow us to separate the individual contributions of the forcings to the AHP intensity, as well as those of their synergies. We confirm the predominant role of the orbital forcing in the strength of the last five AHPs, and our simulations agree with previous estimates of a threshold in orbital forcing above which an AHP develops. Moreover, we show that atmospheric CO2 and the extent of ice sheets can also add up to be as important as the orbital parameters. High values of CO2, past a 205 ppm threshold, and low values of ice sheets extent, below an 8 % of global land surface threshold, yield the AHPs with the most precipitation and vegetation. Additionally, our results show that AHPs differ not only in amplitude, but also in their speed of change, and we find that the non-linear vegetation response of AHPs does not correlate with a single forcing and that the vegetation growth response is faster than its subsequent decline. In regards to future change, an extension of the simulations until the next 50,000 years, shows CO2 to be the main driver of AHPs, with orbital forcing only setting the pace and their intensities being scenario-dependent.

Permafrost ponds are a steady source of methane. However, it is difficult to assess the sensitivity of pond methane emissions to ongoing warming and climate-change-induced drainage, because pond methane emissions show large temporal and spatial variability already on local scale.We study this sensitivity on the landscape level with a new process-based model for Methane Emissions from Ponds (MeEP model), which simulates the three main pathways of methane emissions (diffusion, plant-mediated transport and ebullition) as well as the temperature profile of the water column and the surrounding soils. The model was set up for the polygonal tundra in the Lena River Delta. Due to a temporal resolution of one hour, it is capable of capturing the diurnal, day-to-day and seasonal variability in methane fluxes. MeEP also considers one of the main drivers of spatial variability - ground heterogeneity. Depending on where ponds form in the polygonal tundra, they can be classified as ice-wedge, polygonal-centre or merged-polygonal ponds. In MeEP, each of these pond types is simulated separately and the representation of these ponds was informed by dedicated measurements.The model performance is validated against eddy-covariance measurements of methane fluxes and against in-situ measurements of the aqueous methane concentration, both obtained on Samoylov Island. We will present results regarding the sensitivity of modeled methane emissions from ponds to warming and drainage on the landscape scale.

Between the last glacial maximum (LGM) and preindustrial times (PI), the atmospheric concentration of CH4, as shown by reconstructions from ice cores, roughly doubled. It then doubled again from PI to the present. Ice cores, however, cannot tell us how that development will continue in the future, and ice cores also cannot shed light on the causes of the rise in methane, as well as the rapid fluctuations during periods such as the Bolling-Allerod and Younger Dryas.We use a methane-enabled version of MPI-ESM, the Max Planck Institute for Meteorology Earth System Model, to investigate changes in methane cycling in a transient ESM experiment from the LGM to the present, continuing onwards into the future for the next millennium. The model is driven by prescribed orbit, greenhouse gases and ice sheets, with all other changes to the climate system determined internally. Methane cycling is modelled by modules representing the atmospheric transport and sink of methane, as well as terrestrial sources and sinks from soils, termites, and fires. Thus, the full natural methane cycle – with the exception of geological and animal emissions – is represented in the model. For historical and future climate, anthropogenic emissions of methane are considered, too.We show that the methane increase since the LGM is largely driven by source changes, with LGM emissions substantially reduced in comparison to the early Holocene and preindustrial states due to lower temperature, CO2, and soil carbon. Depending on the future climate scenario, these dependencies then lead to further increases in CH4, with a further doubling of atmospheric CH4 easily possible if one of the higher radiative forcing scenarios is followed. Furthermore, the future increases in CH4 will persist for a long time, as CH4 only decreases when the climate system cools again.

Peatlands are impacted by climate and land-use changes, with feedback to warming by acting as either sources or sinks of carbon. Expert elicitation combined with literature review reveals key drivers of change that alter peatland carbon dynamics, with implications for improving models. The carbon balance of peatlands is predicted to shift from a sink to a source this century. However, peatland ecosystems are still omitted from the main Earth system models that are used for future climate change projections, and they are not considered in integrated assessment models that are used in impact and mitigation studies. By using evidence synthesized from the literature and an expert elicitation, we define and quantify the leading drivers of change that have impacted peatland carbon stocks during the Holocene and predict their effect during this century and in the far future. We also identify uncertainties and knowledge gaps in the scientific community and provide insight towards better integration of peatlands into modelling frameworks. Given the importance of the contribution by peatlands to the global carbon cycle, this study shows that peatland science is a critical research area and that we still have a long way to go to fully understand the peatland-carbon-climate nexus.

We investigate the changes in terrestrial natural methane emissions between the Last Glacial Maximum (LGM) and preindustrial (PI) periods by performing time-slice experiments with a methane-enabled version of MPI-ESM, the Max Planck Institute Earth System Model. We consider all natural sources of methane except for emissions from wild animals and geological sources, i.e. emissions from wetlands, fires, and termites. Changes are dominated by changes in tropical wetland emissions, with mid-to-high-latitude wetlands playing a secondary role, and all other natural sources being of minor importance. The emissions are determined by the interplay of vegetation productivity, a function of CO2 and temperature; source area size, affected by sea level and ice sheet extent; and the state of the West African monsoon, with increased emissions from northern Africa during strong monsoon phases. We show that it is possible to explain the difference in atmospheric methane between LGM and PI purely by changes in emissions. As emissions more than double between LGM and PI, changes in the atmospheric lifetime of CH4, as proposed in other studies, are not required.

The time between the last glacial maximum (LGM) and the present is highly interesting with regard to atmospheric methane. Between the LGM and 10 ka BP atmospheric CH4, as reconstructed from ice cores, nearly doubled, with very rapid concentration changes of about 200 ppb occurring during the Bølling Allerød (BA) and Younger Dryas (YD) transitions. During the Holocene, atmospheric CH4 is very similar for 10 ka BP and PI, but CH4 is about 15% lower in between at 5 ka BP.We use a methane-enabled version of MPI-ESM, the Max Planck Institute Earth System Model, to investigate changes in methane cycling in a transient ESM experiment from the LGM to the present. The model is driven by prescribed orbit, greenhouse gases and ice sheets, with all other changes to the climate system determined internally. Methane cycling is modelled by modules representing the atmospheric transport and sink of methane, as well as terrestrial sources and sinks from soils, termites, and fires. Thus, the full natural methane cycle – with the exception of geological and animal emissions – is represented in the model.Model results are compared to methane concentrations from ice cores, and key periods in climate/methane evolution are highlighted by detailed analyses. Methane concentrations can mainly be explained by emission changes, with LGM emissions substantially reduced in comparison to the early Holocene and preindustrial states due to lower temperature, CO2, and soil carbon. For the large transitions during the deglaciation, such as the transitions from Older Dryas to BA, BA to YD, and YD to Holocene, ocean circulation changes are required to obtain atmospheric methane changes of sufficient magnitude and rapidity.

Northern peatlands store large amounts of carbon (C) and have played an important role in the global carbon cycle since the Last Glacial Maximum. Most northern peatlands have established since the end of the deglaciation and accumulated C over the Holocene, leading to a total present-day stock of 500 ± 100 GtC. This is a consolidated estimate, emerging from a diversity of methods using observational data. Recently, Nichols and Peteet (2019 Nature Geoscience 12: 917-921) presented an estimate of the northern peat C stock of 1055 GtC—exceeding previous estimates by a factor of two. Here, we will review various approaches and estimates of northern peatlands C storage in the literature and consider peat C storage in the context of the Holocene global C budget. We argue that the estimate by Nichols and Peteet is an overestimate, caused by systematic bias introduced by their inclusion of data that are representative for the major peatland regions and of records that lack direct measurements of C density. In particular, some “peatland” sites and data that were included in their synthesis were likely from lacustrine sediments prior to the onset of peat deposits. Furthermore, we argue that their estimate cannot be reconciled within the constraints offered by ice-core and marine records of stable C isotopes and estimated contributions from other processes that affected the terrestrial C storage during the Holocene.

Methane (CH4) is the second most important naturally occurring greenhouse gas (GHG) after carbon dioxide (Myhre G et al 2013 Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (Cambridge: Cambridge University Press) pp 659–740). For both GHGs, the present-day budget is dominated by anthropogenic emissions (Friedlingstein P et al 2019 Earth Syst. Sci. Data 11 1783–838; Saunois M et al 2020 Earth Syst. Sci. Data 12 1561–623). For CO2 it is well established that the projected future rise in atmospheric concentration is near exclusively determined by anthropogenic emissions (Ciais P et al 2013 Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Inter-governmental Panel on Climate Change (Cambridge: Cambridge University Press) pp 465–570). For methane, this appears to be the common assumption, too, but whether this assumption is true has never been shown conclusively. Here, we investigate the evolution of atmospheric methane until 3000 CE under five Shared Socioeconomic Pathway (SSP) scenarios, for the first time using a methane-enabled state-of-the-art Earth System Model (ESM). We find that natural methane emissions, i.e. methane emissions from the biosphere, rise strongly as a reaction to climate warming, thus leading to atmospheric methane concentrations substantially higher than assumed in the scenarios used for CMIP6. We also find that the natural emissions become larger than the anthropogenic ones in most scenarios, showing that natural emissions cannot be neglected.

We investigate the changes in terrestrial natural methane emissions between the Last Glacial Maximum (LGM) and preindustrial (PI) by performing time-slice experiments with a methane-enabled version of MPI-ESM, the Max Planck Institute for Meteorology Earth System Model. We consider all natural sources of methane except for emissions from wild animals and geological sources, i.e. emissions from wetlands, fires, and termites. Changes are dominated by changes in tropical wetland emissions, with mid-to-high latitude wetlands playing a secondary role, and all other natural sources being of minor importance. The emissions are determined by the interplay of vegetation productivity, a function of CO2 and temperature, source area size, affected by sea level and ice sheet extent, and the state of the West African Monsoon, with increased emissions from north Africa during strong monsoon phases. We show that it is possible to explain the difference in atmospheric methane between LGM and PI purely by changes in emissions. As emissions more than double between LGM and PI, changes in the atmospheric lifetime of CH4, as proposed in other studies, are not required.

Northern peatlands have been a persistent natural carbon sink since the last glacial maximum. If there were no limits to their growth, carbon accumulation in these ecosystems could offset a large portion of anthropogenic carbon emissions until the end of the present interglacial period. Evaluation of the limits to northern peatland carbon stocks shows that northern peatlands will potentially play an important role, second only to the oceans, in reducing the atmospheric carbon dioxide concentration to the level that is typical of interglacial periods if cumulative anthropogenic carbon emissions will be kept below 1000 Pg of carbon.

In the version of this Analysis originally published, the following affiliation for A. Lohila was missing: ‘Finnish Meteorological Institute, Climate System Research, Helsinki, Finland’. This affiliation has now been added, and subsequent affiliations renumbered accordingly, in the online versions of the Analysis.

International audience; the accumulation of magma within the Monti Sabatini Volcanic District (MSVD), italy, coupled with the extensional tectonics of the region, pose both volcanic and tectonic hazards to the city of Rome, located 20 km to the southeast. We combine 40 Ar/ 39 Ar geochronology of volcanic deposits and a geomorphologic/stratigraphic/paleomagnetic study of fluvial terraces to determine the recurrence interval and the time elapsed since the last eruption of the MSVD. Moreover, we provide a date for the youngest known eruption of the MSVD and assess the timing of the most recent volcanic phase. Results of this study show: (i) The most recent eruptive phase occurred between 100 ka and 70 ka; (ii) the anomalously high elevation of the MIS 5 terrace indicates that it was concurrent with 50 m of uplift in the volcanic area; (iii) the time since the last eruption (70 ka) exceeds the average recurrence interval (39 ky) in the last 300 ky, as well as the longest previous dormancy (50 ky) in that time span. (iv) the current duration of dormancy is similar to the timespan separating the major explosive phase that occurred 590-450 ka. The magnitude and patterns of deformation of a volcanic field provide constraints on the magmatic processes operating beneath a volcano. A recent geomorphological study 1 reconstructed a series of paleo-surfaces along a 40-km-long stretch of the Tiber River Valley north of Rome, between Magliano Sabina and Monterotondo, located at the eastern margin of the Monti Sabatini Volcanic District (MSVD) (Fig. 1a). These paleo-surfaces are interpreted as fluvial terraces formed through the interplay between regional uplift and glacio-eustasy (e.g. 2) during a regressive phase that formed the hydrographic network of the Paleo-Tiber River since the end of the Santernian (1.78-1.5 Ma; lower Calabrian) 3. The reconstructed rates of uplift during the last 1.8 Ma recognized two major pulses: 0.86 through 0.5 Ma, and 0.25 Ma through the present time 1. The coincidence of the uplift pulses with the ages of the main volcanic phases 1 are interpreted as mainly driven by uprising magma bodies from a metasomatized mantle source of the Roman Magmatic Province (e.g. 4-6). This "magmatic" uplift overlaps a smaller isostatic component on the Tyrrhenian Sea Margin of central Italy. In the present study, we further refine the chronostratigraphy of the paleo-surfaces in the area previously investigated, providing five new 40 Ar/ 39 Ar age determinations on volcanic products intercalated within the

We have developed a new module to calculate soil organic carbon (SOC) accumulation in perennially frozen ground in the land surface model JSBACH. Running this offline version of MPI-ESM we have modelled permafrost carbon accumulation and release from the Last Glacial Maximum (LGM) to the Pre-industrial (PI). Our simulated near-surface PI permafrost extent of 16.9 Mio km2 is close to observational evidence. Glacial boundary conditions, especially ice sheet coverage, result in profoundly different spatial patterns of glacial permafrost extent. Deglacial warming leads to large-scale changes in soil temperatures, manifested in permafrost disappearance in southerly regions, and permafrost aggregation in formerly glaciated grid cells. In contrast to the large spatial shift in simulated permafrost occurrence, we infer an only moderate increase of total LGM permafrost area (18.3 Mio km2) – together with pronounced changes in the depth of seasonal thaw. Reconstructions suggest a larger spread of glacial permafrost towards more southerly regions, but with a highly uncertain extent of non-continuous permafrost. Compared to a control simulation without describing the transport of SOC into perennially frozen ground, the implementation of our newly developed module for simulating permafrost SOC accumulation leads to a doubling of simulated LGM permafrost SOC storage (amounting to a total of ~ 150 PgC). Despite LGM temperatures favouring a larger permafrost extent, simulated cold glacial temperatures – together with low precipitation and low CO2 levels – limit vegetation productivity and therefore prevent a larger glacial SOC build-up in our model. Changes in physical and biogeochemical boundary conditions during deglacial warming lead to an increase in mineral SOC storage towards the Holocene (168 PgC at PI), which is below observational estimates (575 PgC in continuous and discontinuous permafrost). Additional model experiments clarified the sensitivity of simulated SOC storage to model parameters, affecting long-term soil carbon respiration rates and simulated active layer depths. Rather than a steady increase in carbon release from the LGM to PI as a consequence of deglacial permafrost degradation, our results suggest alternating phases of soil carbon accumulation and loss as an effect of dynamic changes in permafrost extent, active layer depths, soil litter input, and heterotrophic respiration.

The climate of Marine Isotope Stage (MIS) 11, the interglacial roughly 400,000 years ago, is investigated for four time slices, 416, 410, 400, and 396 ka. We compare results from two climate models, the earth system model of intermediate complexity CLIMBER2-LPJ and the general circulation model CCSM3, to reconstructions of MIS 11 temperature, precipitation and vegetation, mainly from terrestrial records. The overall picture is that MIS 11 was a relatively warm interglacial in comparison to preindustrial, with Northern Hemisphere (NH) summer temperatures early in MIS 11 (416–410 ka) warmer than preindustrial, though winters were cooler. Later in MIS 11, especially around 400 ka, conditions were cooler in the NH summer, mainly in the high latitudes. Climate changes simulated by the models were mainly driven by insolation changes, with the exception of two local feedbacks that amplify climate changes. Here, the NH high latitudes, where reductions in sea ice cover lead to a winter warming early in MIS 11, as well as the tropics, where monsoon changes lead to stronger climate variations than one would expect on the basis of latitudinal mean insolation change alone, are especially prominent. Both models used in this study support a northward expansion of trees at the expense of grasses in the high northern latitudes early during MIS 11, especially in northern Asia and North America, in line with the available pollen-based reconstructions. With regard to temperature and precipitation changes, there is general agreement between models and reconstructions, but reconstructed precipitation changes are often larger than those simulated by the models. The very limited number of records of sufficiently high resolution and dating quality hinders detailed comparisons between models and reconstructions.

The end of the African Humid Period about 6,000 years ago was associated with vegetation change and decreased precipitation. Conceptual modelling suggests that the nature of the feedback between climate and vegetation is dependent on vegetation type and diversity.

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.

The responses of carbon dioxide (CO2) and other climate variables to an emission pulse of CO2 into the atmosphere are often used to compute the Global Warming Potential (GWP) and Global Temperature change Potential (GTP), to characterize the response timescales of Earth System models, and to build reduced-form models. In this carbon cycle-climate model intercomparison project, which spans the full model hierarchy, we quantify responses to emission pulses of different magnitudes injected under different conditions. The CO2 response shows the known rapid decline in the first few decades followed by a millennium-scale tail. For a 100 Gt-C emission pulse added to a constant CO2 concentration of 389 ppm, 25 ± 9% is still found in the atmosphere after 1000 yr; the ocean has absorbed 59 ± 12% and the land the remainder (16 ± 14%). The response in global mean surface air temperature is an increase by 0.20 ± 0.12 °C within the first twenty years; thereafter and until year 1000, temperature decreases only slightly, whereas ocean heat content and sea level continue to rise. Our best estimate for the Absolute Global Warming Potential, given by the time-integrated response in CO2 at year 100 multiplied by its radiative efficiency, is 92.5 × 10−15 yr W m−2 per kg-CO2. This value very likely (5 to 95% confidence) lies within the range of (68 to 117) × 10−15 yr W m−2 per kg-CO2. Estimates for time-integrated response in CO2 published in the IPCC First, Second, and Fourth Assessment and our multi-model best estimate all agree within 15% during the first 100 yr. The integrated CO2 response, normalized by the pulse size, is lower for pre-industrial conditions, compared to present day, and lower for smaller pulses than larger pulses. In contrast, the response in temperature, sea level and ocean heat content is less sensitive to these choices. Although, choices in pulse size, background concentration, and model lead to uncertainties, the most important and subjective choice to determine AGWP of CO2 and GWP is the time horizon.

The global methane (CH4) budget is becoming an increasingly important component for managing realistic pathways to mitigate climate change. This relevance, due to a shorter atmospheric lifetime and a stronger warming potential than carbon dioxide, is challenged by the still unexplained changes of atmospheric CH4 over the past decade. Emissions and concentrations of CH4 are continuing to increase, making CH4 the second most important human-induced greenhouse gas after carbon dioxide. Two major difficulties in reducing uncertainties come from the large variety of diffusive CH4 sources that overlap geographically, and from the destruction of CH4 by the very short-lived hydroxyl radical (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 research on the methane cycle, and producing regular (similar to biennial) updates of the global methane budget. This consortium includes atmospheric physicists and chemists, biogeochemists of surface and marine emissions, and socio-economists who study anthropogenic emissions. Following Kirschke et al. (2013), we propose here the first version of a living review paper that integrates results of top-down studies (exploiting atmospheric observations within an atmospheric inverse-modelling framework) and bottom-up models, inventories and data-driven approaches (including process-based models for estimating land surface emissions and atmospheric chemistry, and inventories for anthropogenic emissions, data-driven extrapolations). For the 2003-2012 decade, global methane emissions are estimated by top-down inversions at 558 TgCH(4) yr(-1), range 540-568. About 60% of global emissions are anthropogenic (range 50-65 %). Since 2010, the bottom-up global emission inventories have been closer to methane emissions in the most carbon-intensive Representative Concentrations Pathway (RCP8.5) and higher than all other RCP scenarios. Bottom-up approaches suggest larger global emissions (736 TgCH(4) yr(-1), range 596-884) mostly because of larger natural emissions from individual sources such as inland waters, natural wetlands and geological sources. Considering the atmospheric constraints on the top-down budget, it is likely that some of the individual emissions reported by the bottom-up approaches are overestimated, leading to too large global emissions. Latitudinal data from top-down emissions indicate a predominance of tropical emissions (similar to 64% of the global budget, amp;lt;30 degrees N) as compared to mid (similar to 32 %, 30-60 degrees N) and high northern latitudes (similar to 4 %, 60-90 degrees N). Top-down inversions consistently infer lower emissions in China (similar to 58 TgCH(4) yr(-1), range 51-72, -14 %) and higher emissions in Africa (86 TgCH(4) yr(-1), range 73-108, + 19 %) than bottom-up values used as prior estimates. Overall, uncertainties for anthropogenic emissions appear smaller than those from natural sources, and the uncertainties on source categories appear larger for top-down inversions than for bottom-up inventories and models. The most important source of uncertainty on the methane budget is attributable to emissions from wetland and other inland waters. We show that the wetland extent could contribute 30-40% on the estimated range for wetland emissions. Other priorities for improving the methane budget include the following: (i) the development of process-based models for inland-water emissions, (ii) the intensification of methane observations at local scale (flux measurements) to constrain bottom-up land surface models, and at regional scale (surface networks and satellites) to constrain top-down inversions, (iii) improvements in the estimation of atmospheric loss by OH, and (iv) improvements of the transport models integrated in top-down inversions. The data presented here can be downloaded from the Carbon Dioxide Information Analysis Center (http://doi.org/10.3334/CDIAC/GLOBAL_METHANE_BUDGET_2016_V1.1) and the Global Carbon Project. Funding Agencies|Swiss National Science Foundation; NASA [NNX14AF93G, NNX14AO73G]; National Environmental Science Program - Earth Systems and Climate Change Hub; European Commission [283576, 633080]; ESA Climate Change Initiative Greenhouse Gases Phase 2 project; US Department of Energy, BER [DE-AC02-05CH11231]; FAO member countries; Environment Research and Technology Development Fund of the Ministry of the Environment, Japan [2-1502]; ERC [322998]; NERC [NE/J00748X/1]; Swedish Research Council VR; Research Council of Norway [244074]; NSF [1243232, 1243220]; National Science and Engineering Research Council of Canada (NSERC); Chinas QianRen Program; CSIRO Australia; Australian Bureau of Meteorology; Australian Institute of Marine Science; Australian Antarctic Division; NOAA USA; Meteorological Service of Canada; National Aeronautic and Space Administration (NASA) [NAG5-12669, NNX07AE89G, NNX11AF17G, NNX07AE87G, NNX07AF09G, NNX11AF15G, NNX11AF16G]; Department of Energy and Climate Change (DECC, UK) [GA01081]; Commonwealth Scientific and Industrial Research Organization (CSIRO Australia); Bureau of Meteorology (Australia); Joint DECC/Defra Met Office Hadley Centre Climate Programme [GA01101]

A consistent, process-based methane module for a global land surface scheme has been developed which is general enough to be applied in permafrost regions as well as wetlands outside permafrost areas. Methane production, oxidation and transport by ebullition, diffusion and plants are represented. Oxygen has been explicitly incorporated in diffusion, transport by plants and two oxidation processes, of which one uses soil oxygen, while the other uses oxygen that is available via roots. Permafrost and wetland soils show special behaviour, such as variable soil pore space due to freezing and thawing or water table depths due to changing soil water content. This has been integrated directly into the methane-related processes. A detailed application at the polygonal tundra site Samoylov, Lena delta, Russia, is used for evaluation purposes. The application at Samoylov also shows differences in the importance of the several transport processes and in the methane dynamics under varying soil moisture, ice and temperature conditions during different seasons and on different microsites. These microsites are the elevated moist polygonal rim and the depressed wet polygonal center. The evaluation shows sufficiently good agreement with field observations despite the fact that the module has not been specifically calibrated to these data.

Boreal and subarctic peatlands are an important dynamical component of the earth system. They are sensitive to climate change and could either continue to serve as a carbon sink or become a carbon source. Climatic thresholds for switching peatlands from sink to source are not well defined and therefore, incorporating peatlands into Earth system models is a challenging task. Here we introduce a climatic index, warm precipitation excess, to delineate the potential geographic distribution of boreal peatlands for a given climate and landscape morphology. This allows us to explain the present-day distribution of peatlands in Western Siberia, their absence during the Last Glacial Maximum, their expansion during the mid-Holocene and to form a working hypothesis about the trend to peatland degradation in the southern taiga belt of Western Siberia under an RCP 8.5 scenario for the projected climate in year 2100.

International audience; In the version of this Article originally published, data given for total exceedance budgets of CO2 for 1.5 °C and 2 °C targets were incorrect in the main text, although the correct values were given in Supplementary Table 1. These errors also resulted in an incorrect estimation of the percentage difference between the authors' results and estimates by the IPCC. These errors have now been corrected in the online versions.

The role of a reduction in the Atlantic meridional overturning and that of a persistently negative North Atlantic Oscillation in explaining the coldness of the European Little Ice Age (LIA) has been assessed in two sets of numerical experiments. These experiments are performed using an intermediate complexity climate model and a full complexity GCM. The reduction in the Meridional Overturning Circulation (MOC) of ca. 25% is triggered by a conventional fresh-water hosing set-up. A persistently negative NAO winter circulation, at NAO-index value -0.5, is imposed using recently developed data-assimilation techniques applicable on paleoclimatic timescales. The hosing experiments lead to a reduction in oceanic meridional heat transport and cooler sea-surface temperatures. Next to a direct cooling effect on European climate, the change in ocean surface temperatures feedback on the atmospheric circulation modifying European climate significantly. The data-assimilation experiments showed a reduction of winter temperatures over parts of Europe, but there is little persistence into the summer season. The output of all model experiments are compared to reconstructions of winter and summer temperature based on the available temperature data for the LIA period. This demonstrates that the hypothesis of a persistently negative NAO as an explanation for the European LIA does not hold. The hosing experiments do not clearly support the hypothesis that a reduction in the MOC is the primary driver of LIA climate change. However, a reduction in the Atlantic overturning might have been a cause of the European LIA climate, depending on whether there is a strong enough feedback on the atmospheric circulation. © 2010 The Author(s).

An approach to considering changes in flooding probability in the integrated assessment of climate change is introduced. A reduced-form hydrological model for flood prediction and a downscaling approach suitable for integrated assessment modeling are presented. Based on these components, the fraction of world population living in river basins affected by changes in flooding probability in the course of climate change is determined. This is then used as a climate impact response function in order to derive emission corridors limiting the population affected. This approach illustrates the consideration of probabilistic impacts within the framework of the tolerable windows approach. Based on the change in global mean temperature, as calculated by the simple climate models used in integrated assessment, spatially resolved changes in climatic variables are determined using pattern scaling, while natural variability in these variables is considered using twentieth century deviations from the climatology. Driven by the spatially resolved climate change, the hydrological model then aggregates these changes to river basin scale. The hydrological model is subjected to a sensitivity analysis with regard to the water balance, and the uncertainty arising through the different projections of changes in mean climate by differing climate models is considered by presenting results based on different models. The results suggest that up to 20% of world population live in river basins that might inevitably be affected by increased flood events in the course of global warming, depending on the climate model used to estimate the regional distribution of changes in climate.

Anthropogenic land-cover change (ALCC) is one of the few climate forcings for which the net direction of the climate response over the last two centuries is still not known. The uncertainty is due to the often counteracting temperature responses to the many biogeophysical effects and to the biogeochemical versus biogeophysical effects. Palaeoecological studies show that the major transformation of the landscape by anthropogenic activities in the southern zone of the Baltic Sea basin occurred between 6000 and 3000/2500 cal year BP. The only modelling study of the biogeophysical effects of past ALCCs on regional climate in north-western Europe suggests that deforestation between 6000 and 200 cal year BP may have caused significant change in winter and summer temperature. There is no indication that deforestation in the Baltic Sea area since AD 1850 would have been a major cause of the recent climate warming in the region through a positive biogeochemical feedback. Several model studies suggest that boreal reforestation might not be an effective climate warming mitigation tool as it might lead to increased warming through biogeophysical processes.

Global wetlands are believed to be climate sensitive, and are the largest natural emitters of methane (CH4). Increased wetland CH4 emissions could act as a positive feedback to future warming. The Wetland and Wetland CH4 Inter-comparison of Models Project (WETCHIMP) investigated our present ability to simulate large-scale wetland characteristics and corresponding CH4 emissions. To ensure inter-comparability, we used a common experimental protocol driving all models with the same climate and carbon dioxide (CO2) forcing datasets. The WETCHIMP experiments were conducted for model equilibrium states as well as transient simulations covering the last century. Sensitivity experiments investigated model response to changes in selected forcing inputs (precipitation, temperature, and atmospheric CO2 concentration). Ten models participated, covering the spectrum from simple to relatively complex, including models tailored either for regional or global simulations. The models also varied in methods to calculate wetland size and location, with some models simulating wetland area prognostically, while other models relied on remotely sensed inundation datasets, or an approach intermediate between the two. Four major conclusions emerged from the project. First, the suite of models demonstrate extensive disagreement in their simulations of wetland areal extent and CH4 emissions, in both space and time. Simple metrics of wetland area, such as the latitudinal gradient, show large variability, principally between models that use inundation dataset information and those that independently determine wetland area. Agreement between the models improves for zonally summed CH4 emissions, but large variation between the models remains. For annual global CH4 emissions, the models vary by ±40% of the all-model mean (190 Tg CH4 yr−1). Second, all models show a strong positive response to increased atmospheric CO2 concentrations (857 ppm) in both CH4 emissions and wetland area. In response to increasing global temperatures (+3.4 °C globally spatially uniform), on average, the models decreased wetland area and CH4 fluxes, primarily in the tropics, but the magnitude and sign of the response varied greatly. Models were least sensitive to increased global precipitation (+3.9 % globally spatially uniform) with a consistent small positive response in CH4 fluxes and wetland area. Results from the 20th century transient simulation show that interactions between climate forcings could have strong non-linear effects. Third, we presently do not have sufficient wetland methane observation datasets adequate to evaluate model fluxes at a spatial scale comparable to model grid cells (commonly 0.5°). This limitation severely restricts our ability to model global wetland CH4 emissions with confidence. Our simulated wetland extents are also difficult to evaluate due to extensive disagreements between wetland mapping and remotely sensed inundation datasets. Fourth, the large range in predicted CH4 emission rates leads to the conclusion that there is both substantial parameter and structural uncertainty in large-scale CH4 emission models, even after uncertainties in wetland areas are accounted for.

International audience; The Wetland and Wetland CH4 Intercomparison of Models Project (WETCHIMP) was created to evaluate our present ability to simulate large-scale wetland characteristics and corresponding methane (CH4) emissions. A multi-model comparison is essential to evaluate the key uncertainties in the mechanisms and parameters leading to methane emissions. Ten modelling groups joined WETCHIMP to run eight global and two regional models with a common experimental protocol using the same climate and atmospheric carbon dioxide (CO2) forcing datasets. We reported the main conclusions from the intercomparison effort in a companion paper (Melton et al., 2012). Here we provide technical details for the six experiments, which included an equilibrium, a transient, and an optimized run plus three sensitivity experiments (temperature, precipitation, and atmospheric CO2 concentration). The diversity of approaches used by the models is summarized through a series of conceptual figures, and is used to evaluate the wide range of wetland extents and CH4 fluxes predicted by the models in the equilibrium run. We discuss relationships among the various approaches and patterns in consistencies of these model predictions. Within this group of models, there are three broad classes of methods used to estimate wetland extent: prescribed based on wetland distribution maps, prognostic relationships between hydrological states based on satellite observations, and explicit hydrological mass balances. A larger variety of approaches was used to estimate the net CH4 fluxes from wetland systems. Even though modelling of wetland extents and CH4 emissions has progressed significantly over recent decades, large uncertainties still exist when estimating CH4 emissions: there is little consensus on model structure or complexity due to knowledge gaps, different aims of the models, and the range of temporal and spatial resolutions of the models.

[1] A method is introduced to estimate the proximity of climate sub-systems to non-linear thresholds. We suggest to measure the smallest decay rate of the system under investigation and to consider its trend. We argue that this is the diagnostic variable most directly linked to the distance from a bifurcation threshold. With the climate model of intermediate complexity CLIMBER2 we demonstrate our method for the North Atlantic thermohaline circulation. It is shown that proper analysis of paleo information could significantly reduce the uncertainty which plagues current estimates of the distance from the shutdown of the thermohaline circulation.

Emission budgets are defined as the cumulative amount of anthropogenic CO2 emission compatible with a global temperature-change target. The simplicity of the concept has made it attractive to policy-makers, yet it relies on a linear approximation of the global carbon–climate system’s response to anthropogenic CO2 emissions. Here we investigate how emission budgets are impacted by the inclusion of CO2 and CH4 emissions caused by permafrost thaw, a non-linear and tipping process of the Earth system. We use the compact Earth system model OSCAR v2.2.1, in which parameterizations of permafrost thaw, soil organic matter decomposition and CO2 and CH4 emission were introduced based on four complex land surface models that specifically represent high-latitude processes. We found that permafrost carbon release makes emission budgets path dependent (that is, budgets also depend on the pathway followed to reach the target). The median remaining budget for the 2 °C target reduces by 8% (1–25%) if the target is avoided and net negative emissions prove feasible, by 13% (2–34%) if they do not prove feasible, by 16% (3–44%) if the target is overshot by 0.5 °C and by 25% (5–63%) if it is overshot by 1 °C. (Uncertainties are the minimum-to-maximum range across the permafrost models and scenarios.) For the 1.5 °C target, reductions in the median remaining budget range from ~10% to more than 100%. We conclude that the world is closer to exceeding the budget for the long-term target of the Paris Climate Agreement than previously thought. © 2018, The Author(s), under exclusive licence to Springer Nature Limited.

The ongoing development of the Global Carbon Project (GCP) global methane (CH4 ) budget shows a continuation of increasing CH4 emissions and CH4 accumulation in the atmosphere during 2000-2017. Here, we decompose the global budget into 19 regions (18land and 1 oceanic) and five key source sectors to spatially attribute the observed global trends. A comparison of top-down (TD) (atmospheric and transport model-based) and bottom-up (BU) (inventory- and process model-based) CH4 emission estimates demonstrates robust temporal trends with CH4 emissions increasing in 16 of the 19 regions. Five regions-China, Southeast Asia, USA, South Asia, and Brazil-account for >40% of the global total emissions (their anthropogenic and natural sources together totaling >270Tg CH4 yr-1 in 2008-2017). Two of these regions, China and South Asia, emit predominantly anthropogenic emissions (>75%) and together emit more than 25% of global anthropogenic emissions. China and the Middle East show the largest increases in total emission rates over the 2000 to 2017 period with regional emissions increasing by >20%. In contrast, Europe and Korea and Japan show a steady decline in CH4 emission rates, with total emissions decreasing by ~10% between 2000 and 2017. Coal mining, waste (predominantly solid waste disposal) and livestock (especially enteric fermentation) are dominant drivers of observed emissions increases while declines appear driven by a combination of waste and fossil emission reductions. As such, together these sectors present the greatest risks of further increasing the atmospheric CH4 burden and the greatest opportunities for greenhouse gas abatement.

Approximately 1700 Pg of soil carbon (C) are stored in the northern circumpolar permafrost zone, more than twice as much C than in the atmosphere. The overall amount, rate, and form of C released to the atmosphere in a warmer world will influence the strength of the permafrost C feedback to climate change. We used a survey to quantify variability in the perception of the vulnerability of permafrost C to climate change. Experts were asked to provide quantitative estimates of permafrost change in response to four scenarios of warming. For the highest warming scenario (RCP 8.5), experts hypothesized that C release from permafrost zone soils could be 19-45 Pg C by 2040, 162-288 Pg C by 2100, and 381-616 Pg C by 2300 in CO2 equivalent using 100-year CH4 global warming potential (GWP). These values become 50 % larger using 20-year CH4 GWP, with a third to a half of expected climate forcing coming from CH4 even though CH4 was only 2.3 % of the expected C release. Experts projected that two-thirds of this release could be avoided under the lowest warming scenario (RCP 2.6). These results highlight the potential risk from permafrost thaw and serve to frame a hypothesis about the magnitude of this feedback to climate change. However, the level of emissions proposed here are unlikely to overshadow the impact of fossil fuel burning, which will continue to be the main source of C emissions and climate forcing. © 2013 The Author(s).

Permafrost soils in the high northern latitudes contain a substantial amount of carbon which is not decomposed due to frozen conditions. Climate change will lead to a thawing of at least part of the permafrost, implying that the stored carbon will become accessible to decomposition and be released to the atmosphere. We use a land surface model to quantify the amount of carbon released up until 2300 and determine the net carbon balance of the northern hemisphere permafrost region under climate warming following the RCP scenarios 2.6, 4.5, and 8.5. Here we show for the first time that the net carbon balance of the permafrost region is not just strongly dependent on the overall warming, but also on the CO2 concentration pathway. As a result moderate warming scenarios may counter- intuitively lead to lower net carbon emissions from the permafrost region than low warming scenarios.