Search

The atmospheric CO2 concentration increased by about 20 ppm from 6000 BCE to the pre-industrial period (1850 CE). Several hypotheses have been proposed to explain mechanisms of this CO2 growth based on either ocean or land carbon sources. Here, we apply the Earth system model MPI-ESM-LR for two transient simulations of climate and carbon cycle dynamics during this period. In the first simulation, atmospheric CO2 is prescribed following ice-core CO2 data. In response to the growing atmospheric CO2 concentration, land carbon storage increases until 2000 BCE, stagnates afterwards, and decreases from 1 CE, while the ocean continuously takes CO2 out of the atmosphere after 4000 BCE. This leads to a missing source of 166 Pg of carbon in the ocean–land–atmosphere system by the end of the simulation. In the second experiment, we applied a CO2 nudging technique using surface alkalinity forcing to follow the reconstructed CO2 concentration while keeping the carbon cycle interactive. In that case the ocean is a source of CO2 from 6000 to 2000 BCE due to a decrease in the surface ocean alkalinity. In the prescribed CO2 simulation, surface alkalinity declines as well. However, it is not sufficient to turn the ocean into a CO2 source. The carbonate ion concentration in the deep Atlantic decreases in both the prescribed and the interactive CO2 simulations, while the magnitude of the decrease in the prescribed CO2 experiment is underestimated in comparison with available proxies. As the land serves as a carbon sink until 2000 BCE due to natural carbon cycle processes in both experiments, the missing source of carbon for land and atmosphere can only be attributed to the ocean. Within our model framework, an additional mechanism, such as surface alkalinity decrease, for example due to unaccounted for carbonate accumulation processes on shelves, is required for consistency with ice-core CO2 data. Consequently, our simulations support the hypothesis that the ocean was a source of CO2 until the late Holocene when anthropogenic CO2 sources started to affect atmospheric CO2 .

Reconstructions of the global mean annual temperature evolution during the Holocene yield conflicting results. One temperature reconstruction shows global cooling during the late Holocene. The other reconstruction reveals global warming. Here we show that both a global warming mode and a cooling mode emerge when performing a spatio-temporal analysis of annual temperature variability during the Holocene using data from a transient climate model simulation. The warming mode is most pronounced in the tropics. The simulated cooling mode is determined by changes in the seasonal cycle of Arctic sea-ice that are forced by orbital variations and volcanic eruptions. The warming mode dominates in the mid-Holocene, whereas the cooling mode takes over in the late Holocene. The weighted sum of the two modes yields the simulated global temperature trend evolution. Our findings have strong implications for the interpretation of proxy data and the selection of proxy locations to compute global mean temperatures., Proxy reconstructions show a decreasing trend from the Middle to Late Holocene, which conflicts with model results showing an increasing trend. Statistical analysis of model output shows that these conflicting results originate from two distinct modes of variability, which dominate at different regions and times.

ICON‐A is the new icosahedral nonhydrostatic (ICON) atmospheric general circulation model in a configuration using the Max Planck Institute physics package, which originates from the ECHAM6 general circulation model, and has been adapted to account for the changed dynamical core framework. The coupling scheme between dynamics and physics employs a sequential updating by dynamics and physics, and a fixed sequence of the physical processes similar to ECHAM6. To allow a meaningful initial comparison between ICON‐A and the established ECHAM6‐LR model, a setup with similar, low resolution in terms of number of grid points and levels is chosen. The ICON‐A model is tuned on the base of the Atmospheric Model Intercomparison Project (AMIP) experiment aiming primarily at a well balanced top‐of atmosphere energy budget to make the model suitable for coupled climate and Earth system modeling. The tuning addresses first the moisture and cloud distribution to achieve the top‐of‐atmosphere energy balance, followed by the tuning of the parameterized dynamic drag aiming at reduced wind errors in the troposphere. The resulting version of ICON‐A has overall biases, which are comparable to those of ECHAM6. Problematic specific biases remain in the vertical distribution of clouds and in the stratospheric circulation, where the winter vortices are too weak. Biases in precipitable water and tropospheric temperature are, however, reduced compared to the ECHAM6. ICON‐A will serve as the basis of further development and as the atmosphere component to the coupled model, ICON‐Earth system model (ESM).

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

A major link between climate and humans in tropical northern Africa, and the Sahel in particular, is land use and associated land cover change, mainly where subsistence farming prevails. Here we assess possible feedbacks between the type of land use and harvest intensity and climate by analysing a series of idealized GCM experiments using the Max Planck Institute Earth System Model (MPI-ESM). The baseline for these experiments is a simulation forced by the RCP8.5 (radiation concentration pathway) scenario, which includes strong greenhouse gas emissions and anthropogenic land cover changes. The anthropogenic land cover changes in the RCP8.5 scenario include a mixture of pasture and agriculture. In subsequent simulations, we replace the entire area affected by anthropogenic land cover change in the region between the Sahara in the north and the Guinean Coast in the south (4 to 20\\degree N) with either pasture or agriculture. In a second set-up we vary the amount of harvest in the case of agriculture. The RCP8.5 baseline simulation reveals strong changes in the area mean agriculture and monsoon rainfall. In comparison with these changes, any variation of the type of land use in the study area leads to very small, mostly insignificantly small, additional differences in mean temperature and annual precipitation change in this region. These findings are only based on the specific set-up of our experiments, which only focuses on variations in the kind of land use, and not the increase in land use, over the 21st century, nor whether land use is considered at all. Within the uncertainty of the representation of land use in current ESMs, our study suggests marginal feedback between land use changes and climate changes triggered by strong greenhouse gas emissions. Hence as a good approximation, climate can be considered as an external forcing: models investigating land-use–conflict dynamics can run offline by prescribing seasonal or mean values of climate as a boundary condition for climate.

Recent advances in the representation of soil carbon decomposition (Goll et al., 2015) and carbon-nitrogen interactions (Parida, 2011; Goll et al., 2012) implemented previously into separate versions of the land surface scheme JSBACH are here combined in a single version which is set to be used in the upcoming 6th phase of coupled model intercomparison project (CMIP6) (Eyring et al., 2016). Here we demonstrate that the new version of JSBACH is able to reproduce the spatial variability in the reactive nitrogen loss pathways as derived from a compilation of δ15N data (r=.63, RMSE=.26, Taylor score=.81). The inclusion of carbon-nitrogen interactions leads to a moderate reduction (−10 %) of the carbon-concentration feedback (βL) and has a negligible effect on the sensitivity of the land carbon cycle to warming (γL) compared to the same version of the model without carbon-nitrogen interactions in idealized simulations (1 % increase in atmospheric carbon dioxide per yr). In line with evidence from elevated carbon dioxide manipulation experiments (Shi et al., 2015; Liang et al., 2016), pronounced nitrogen scarcity is alleviated by (1) the accumulation of nitrogen due to enhanced nitrogen inputs by biological nitrogen fixation and reduced losses by leaching and volatilization as well as the (2) enhanced turnover of organic nitrogen. The strengths of the land carbon feedbacks of the recent version of JSBACH, with βL=0.61 Pg ppm−1 and γL=−27.5 Pg °C−1, are 34 % and 53 % less than the averages of CMIP5 models (Arora et al., 2013), although the CMIP5 version of JSBACH simulated βL and γL which are 59 % and 42 % higher than multi-model average. These changes are primarily due to the new decomposition model, stressing the importance of getting the basics right (here: the decomposition of soil carbon) before increasing the complexity of the model (here: carbon-nitrogen interactions).

[1] The new Max-Planck-Institute Earth System Model (MPI-ESM) is used in the Coupled Model Intercomparison Project phase 5 (CMIP5) in a series of climate change experiments for either idealized CO2-only forcing or forcings based on observations and the Representative Concentration Pathway (RCP) scenarios. The paper gives an overview of the model configurations, experiments related forcings, and initialization procedures and presents results for the simulated changes in climate and carbon cycle. It is found that the climate feedback depends on the global warming and possibly the forcing history. The global warming from climatological 1850 conditions to 2080–2100 ranges from 1.5°C under the RCP2.6 scenario to 4.4°C under the RCP8.5 scenario. Over this range, the patterns of temperature and precipitation change are nearly independent of the global warming. The model shows a tendency to reduce the ocean heat uptake efficiency toward a warmer climate, and hence acceleration in warming in the later years. The precipitation sensitivity can be as high as 2.5% K−1 if the CO2 concentration is constant, or as small as 1.6% K−1, if the CO2 concentration is increasing. The oceanic uptake of anthropogenic carbon increases over time in all scenarios, being smallest in the experiment forced by RCP2.6 and largest in that for RCP8.5. The land also serves as a net carbon sink in all scenarios, predominantly in boreal regions. The strong tropical carbon sources found in the RCP2.6 and RCP8.5 experiments are almost absent in the RCP4.5 experiment, which can be explained by reforestation in the RCP4.5 scenario.

The scientific understanding of the Earth's climate system, including the central question of how the climate system is likely to respond to human-induced perturbations, is comprehensively captured in GCMs and Earth System Models(ESM). Diagnosing the simulated climate response, and comparing responses across different models, is crucially dependent on transparent assumptions of how the GCM/ESM has been driven – especially because the implementation can involve subjective decisions and may differ between modelling groups performing the same experiment. This paper outlines the climate forcings and setup of the Met Office Hadley Centre ESM, HadGEM2-ES for the CMIP5 set of centennial experiments. We document the prescribed greenhouse gas concentrations, aerosol precursors, stratospheric and tropospheric ozone assumptions, as well as implementation of land-use change and natural forcings for the HadGEM2-ES historical and future experiments following the Representative Concentration Pathways. In addition, we provide details of how HadGEM2-ES ensemble members were initialised from the control run and how the palaeoclimate and AMIP experiments, as well as the "emission-driven" RCP experiments were performed.

Coordinated experimental design and implementation has become a cornerstone of global climate modelling. So-called Model Intercomparison Projects (MIPs) enable systematic and robust analysis of results across many models to identify common signals and understand model similarities and differences without being hindered by ad-hoc differences in model set-up or experimental boundary conditions. The activity known as the Coupled Model Intercomparison Project (CMIP) has thus grown significantly in scope and as it enters its 6th phase, CMIP6, the design and documentation of individual simulations has been devolved to individual climate science communities. The Coupled Climate-Carbon Cycle Model Intercomparison Project (C4MIP) takes responsibility for design, documentation and analysis of carbon cycle feedbacks and interactions in climate simulations. These feedbacks are potentially large and play a leading order contribution in determining the atmospheric composition in response to human emissions of CO2 and in the setting of emissions targets to stabilise climate or avoid dangerous climate change. For over a decade C4MIP has coordinated coupled climate-carbon cycle simulations and in this paper we describe the C4MIP simulations that will be formally part of CMIP6. While the climate-carbon cycle community has formed this experimental design the simulations also fit into the wider CMIP activity and conform to some common standards such as documentation and diagnostic requests and are designed to complement the CMIP core experiments known as the DECK. C4MIP has 3 key strands of scientific motivation and the requested simulations are designed to satisfy their needs: (1) pre-industrial and historical simulations (formally part of the common set of CMIP6 experiments) to enable model evaluation; (2) idealised coupled and partially-coupled simulations with 1 % per year increases in CO2 to enable diagnosis of feedback strength and its components; (3) future scenario simulations to project how the Earth System will respond over the 21st century and beyond to anthropogenic activity. This paper documents in detail these simulations, explains their rationale and planned analysis, and describes how to set-up and run the simulations. Particular attention is paid to boundary conditions and input data required, and also the output diagnostics requested. It is important that modelling groups participating in C4MIP adhere as closely as possible to this experimental design.

International audience; The quantification of sources and sinks of carbon from land use and land cover changes (LULCC) is uncertain. We investigated how the parametrization of LULCC and of organic matter decomposition, as well as initial land cover, affects the historical and future carbon fluxes in an Earth System Model (ESM). Using the land component of the Max Planck Institute ESM, we found that the historical (1750-2010) LULCC flux varied up to 25% depending on the fraction of biomass which enters the atmosphere directly due to burning or is used in short-lived products. The uncertainty in the decadal LULCC fluxes of the recent past due to the parametrization of decomposition and direct emissions was 0.6 Pg C yr −1 , which is 3 times larger than the uncertainty previously attributed to model and method in general. Preindustrial natural land cover had a larger effect on decadal LULCC fluxes than the aforementioned parameter sensitivity (1.0 Pg C yr −1). Regional differences between reconstructed and dynamically computed land covers, in particular, at low latitudes, led to differences in historical LULCC emissions of 84-114 Pg C, globally. This effect is larger than the effects of forest regrowth, shifting cultivation, or climate feedbacks and comparable to the effect of differences among studies in the terminology of LULCC. In general, we find that the practice of calibrating the net land carbon balance to provide realistic boundary conditions for the climate component of an ESM hampers the applicability of the land component outside its primary field of application.

Climatic Change, 109 (3-4), ISSN:0165-0009, ISSN:1573-1480

Thum, T., et al. (2011), 'Soil carbon model alternatives for ECHAM5/JSBACH climate model: Evaluation and impacts on global carbon cycle estimates', Journal of Geophysical Research, Vol. 116, G02028, published 29 June 2011. The version of record is available at doi:10.1029/2010JG001612 © 2011 American Geophysical Union.

We investigate the relative importance of climate change (CC) and anthropogenic land cover change (ALCC) for the dust emissions and burden changes between the late nineteenth century and today. For this purpose, the climate-aerosol model ECHAM6-HAM2 is complemented by a new scheme to derive potential dust sources at runtime using the vegetation cover provided by the land component JSBACH of ECHAM6. Dust emissions are computed online using information from the ECHAM6 atmospheric component. This allows us to account for changes in land cover and climate interactively and to distinguish between emissions from natural and agricultural dust sources. For today's climate we find that nearly 10% of dust particles are emitted from agricultural areas. According to our simulations, global annual dust emissions have increased by 25% between the late nineteenth century and today (e. g., from 729 Tg/a to 912 Tg/a). Globally, CC and ALCC (e. g., agricultural expansion) have both contributed to this change (56% and 40%, respectively). There are however large regional differences. For example, change in dust emissions in Africa are clearly dominated by CC. Global dust burden have increased by 24.5% since the late nineteenth century, which results in a clear-sky radiative forcing at top of the atmosphere of similar to 0.14W/m(2). Based on these findings, we recommend that both climate changes and anthropogenic land cover changes should be considered when investigating long-term changes in dust emissions.

Eleven coupled climate–carbon cycle models used a common protocol to study the coupling between climate change and the carbon cycle. The models were forced by historical emissions and the Intergovernmental Panel on Climate Change (IPCC) Special Report on Emissions Scenarios (SRES) A2 anthropogenic emissions of CO2 for the 1850–2100 time period. For each model, two simulations were performed in order to isolate the impact of climate change on the land and ocean carbon cycle, and therefore the climate feedback on the atmospheric CO2 concentration growth rate. There was unanimous agreement among the models that future climate change will reduce the efficiency of the earth system to absorb the anthropogenic carbon perturbation. A larger fraction of anthropogenic CO2 will stay airborne if climate change is accounted for. By the end of the twenty-first century, this additional CO2 varied between 20 and 200 ppm for the two extreme models, the majority of the models lying between 50 and 100 ppm. The higher CO2 levels led to an additional climate warming ranging between 0.1° and 1.5°C. All models simulated a negative sensitivity for both the land and the ocean carbon cycle to future climate. However, there was still a large uncertainty on the magnitude of these sensitivities. Eight models attributed most of the changes to the land, while three attributed it to the ocean. Also, a majority of the models located the reduction of land carbon uptake in the Tropics. However, the attribution of the land sensitivity to changes in net primary productivity versus changes in respiration is still subject to debate; no consensus emerged among the models.

Based on a recent reconstruction of anthropogenic land cover change (ALCC), we derive the associated CO 2 emissions since 800 AD by two independent methods: a bookkeeping approach and a process model. The results are compared with the pre-industrial development of atmospheric CO 2 known from antarctic ice cores. Our results show that pre-industrial CO 2 emissions from ALCC have been relevant for the pre-industrial carbon cycle, although before 1750 AD their trace in atmospheric CO 2 is obscured by other processes of similar magnitude. After 1750 AD, the situation is different: the steep increase in atmospheric CO 2 until 1850 AD—this is before fossil fuel emissions rose to significant values—is to a substantial part explained by growing emissions from ALCC. DOI: 10.1111/j.1600-0889.2010.00479.x To access the supplementary material to this article please see Supplementary files in the column to the right (under Article Tools).

The sensitivity of the climate–biogeochemistry system to volcanic eruptions is investigated using the comprehensive Earth System Model developed at the Max Planck Institute for Meteorology. The model includes an interactive carbon cycle with modules for terrestrial biosphere as well as ocean biogeochemistry. The volcanic forcing is based on a recent reconstruction for the last 1200 yr. An ensemble of five simulations is performed and the averaged response of the system is analysed in particular for the largest eruption of the last millennium in the year 1258. After this eruption, the global annual mean temperature drops by 1 K and recovers slowly during 10 yr. Atmospheric CO 2 concentration declines during 4 yr after the eruption by ca. 2 ppmv to its minimum value and then starts to increase towards the pre-eruption level. This CO 2 decrease is explained mainly by reduced heterotrophic respiration on land in response to the surface cooling, which leads to increased carbon storage in soils, mostly in tropical and subtropical regions. The ocean acts as a weak carbon sink, which is primarily due to temperature-induced solubility. This sink saturates 2 yr after the eruption, earlier than the land uptake. DOI: 10.1111/j.1600-0889.2010.00471.x

The Max Planck Institute Grand Ensemble (MPI‐GE) is the largest ensemble of a single comprehensive climate model currently available, with 100 members for the historical simulations (1850–2005) and four forcing scenarios. It is currently the only large ensemble available that includes scenario representative concentration pathway (RCP) 2.6 and a 1% CO2 scenario. These advantages make MPI‐GE a powerful tool. We present an overview of MPI‐GE, its components, and detail the experiments completed. We demonstrate how to separate the forced response from internal variability in a large ensemble. This separation allows the quantification of both the forced signal under climate change and the internal variability to unprecedented precision. We then demonstrate multiple ways to evaluate MPI‐GE and put observations in the context of a large ensemble, including a novel approach for comparing model internal variability with estimated observed variability. Finally, we present four novel analyses, which can only be completed using a large ensemble. First, we address whether temperature and precipitation have a pathway dependence using the forcing scenarios. Second, the forced signal of the highly noisy atmospheric circulation is computed, and different drivers are identified to be important for the North Pacific and North Atlantic regions. Third, we use the ensemble dimension to investigate the time dependency of Atlantic Meridional Overturning Circulation variability changes under global warming. Last, sea level pressure is used as an example to demonstrate how MPI‐GE can be utilized to estimate the ensemble size needed for a given scientific problem and provide insights for future ensemble projects.