scholarly journals Simulation of climate, ice sheets and CO<sub>2</sub> evolution during the last four glacial cycles with an Earth system model of intermediate complexity

2017 ◽  
Author(s):  
Andrei Ganopolski ◽  
Victor Brovkin

Abstract. In spite of significant progress in paleoclimate reconstructions and modeling of different aspects of the past glacial cycles, the mechanisms which transform regional and seasonal variations in solar insolation into long-term and global-scale glacial-interglacial cycles are still not fully understood, in particular, for CO2 variability. Here using the Earth system model of intermediate complexity CLIMBER-2 we performed simulations of co-evolution of climate, ice sheets and carbon cycle over the last 400,000 years using the orbital forcing as the only external forcing. The model simulates temporal dynamics of CO2, global ice volume and other climate system characteristics in good agreement with paleoclimate reconstructions. Using simulations performed with the model in different configurations, we also analyze the role of individual processes and sensitivity to the choice of model parameters. While many features of simulated glacial cycles are rather robust, some details of CO2 evolution, especially during glacial terminations, are rather sensitive to the choice of model parameters. Specifically, we found two major regimes of CO2 changes during terminations: in the first one, when the recovery of the Atlantic meridional overturning circulation (AMOC) occurs only at the end of the termination, a pronounced overshoot in CO2 concentration occurs at the beginning of the interglacial and CO2 remains almost constant during interglacial or even decline towards the end, resembling Eemian CO2 dynamics. However, if the recovery of the AMOC occurs in the middle of the glacial termination, CO2 concentration continues to rise during interglacial, similar to Holocene. We also discuss potential contribution of the brine rejection mechanism for the CO2 and carbon isotopes in the atmosphere and the ocean during the past glacial termination.

2017 ◽  
Vol 13 (12) ◽  
pp. 1695-1716 ◽  
Author(s):  
Andrey Ganopolski ◽  
Victor Brovkin

Abstract. In spite of significant progress in paleoclimate reconstructions and modelling of different aspects of the past glacial cycles, the mechanisms which transform regional and seasonal variations in solar insolation into long-term and global-scale glacial–interglacial cycles are still not fully understood – in particular, in relation to CO2 variability. Here using the Earth system model of intermediate complexity CLIMBER-2 we performed simulations of the co-evolution of climate, ice sheets, and carbon cycle over the last 400 000 years using the orbital forcing as the only external forcing. The model simulates temporal dynamics of CO2, global ice volume, and other climate system characteristics in good agreement with paleoclimate reconstructions. These results provide strong support for the idea that long and strongly asymmetric glacial cycles of the late Quaternary represent a direct but strongly nonlinear response of the Northern Hemisphere ice sheets to orbital forcing. This response is strongly amplified and globalised by the carbon cycle feedbacks. Using simulations performed with the model in different configurations, we also analyse the role of individual processes and sensitivity to the choice of model parameters. While many features of simulated glacial cycles are rather robust, some details of CO2 evolution, especially during glacial terminations, are sensitive to the choice of model parameters. Specifically, we found two major regimes of CO2 changes during terminations: in the first one, when the recovery of the Atlantic meridional overturning circulation (AMOC) occurs only at the end of the termination, a pronounced overshoot in CO2 concentration occurs at the beginning of the interglacial and CO2 remains almost constant during the interglacial or even declines towards the end, resembling Eemian CO2 dynamics. However, if the recovery of the AMOC occurs in the middle of the glacial termination, CO2 concentration continues to rise during the interglacial, similar to the Holocene. We also discuss the potential contribution of the brine rejection mechanism for the CO2 and carbon isotopes in the atmosphere and the ocean during the past glacial termination.


2012 ◽  
Vol 9 (7) ◽  
pp. 9425-9451 ◽  
Author(s):  
P. B. Holden ◽  
N. R. Edwards ◽  
D. Gerten ◽  
S. Schaphoff

Abstract. We derive a constraint on the strength of CO2 fertilisation of the terrestrial biosphere through a "top-down" approach, calibrating Earth System Model parameters constrained only by the post-industrial increase of atmospheric CO2 concentration. We derive a probabilistic prediction for the globally averaged strength of CO2 fertilisation in nature, implicitly net of other limiting factors such as nutrient availability. The approach yields an estimate that is independent of CO2 enrichment experiments and so provides a new constraint that can in principal be combined with data-driven priors. To achieve this, an essential requirement was the incorporation of a Land Use Change (LUC) scheme into the GENIE earth system model, which we describe in full. Using output from a 671-member ensemble of transient GENIE simulations we build an emulator of the change in atmospheric CO2 concentration change over the preindustrial period (1850 to 2000). We use this emulator to sample the 28-dimensional input parameter space. A Bayesian calibration of the emulator output suggests that the increase in Gross Primary Productivity in response of a doubling of CO2 from preindustrial values is likely to lie in the range 11 to 53%, with a most likely value of 28%. The present-day land-atmosphere flux (1990–2000) is estimated at −0.6 GTC yr−1 (likely in the range 0.9 to −2.0 GTC yr−1). The present-day land-ocean flux (1990–2000) is estimated at −2.2 GTC yr−1 (likely in the range −1.6 to −2.8 GTC yr−1). We estimate cumulative net land emissions over the post-industrial period (land use change emissions net of the CO2 fertilisation sink) to be 37 GTC, likely to lie in the range 130 to −20 GTC.


2010 ◽  
Vol 3 (2) ◽  
pp. 603-633 ◽  
Author(s):  
H. Goosse ◽  
V. Brovkin ◽  
T. Fichefet ◽  
R. Haarsma ◽  
P. Huybrechts ◽  
...  

Abstract. The main characteristics of the new version 1.2 of the three-dimensional Earth system model of intermediate complexity LOVECLIM are briefly described. LOVECLIM 1.2 includes representations of the atmosphere, the ocean and sea ice, the land surface (including vegetation), the ice sheets, the icebergs and the carbon cycle. The atmospheric component is ECBilt2, a T21, 3-level quasi-geostrophic model. The ocean component is CLIO3, which consists of an ocean general circulation model coupled to a comprehensive thermodynamic-dynamic sea-ice model. Its horizontal resolution is of 3° by 3°, and there are 20 levels in the ocean. ECBilt-CLIO is coupled to VECODE, a vegetation model that simulates the dynamics of two main terrestrial plant functional types, trees and grasses, as well as desert. VECODE also simulates the evolution of the carbon cycle over land while the ocean carbon cycle is represented by LOCH, a comprehensive model that takes into account both the solubility and biological pumps. The ice sheet component AGISM is made up of a three-dimensional thermomechanical model of the ice sheet flow, a visco-elastic bedrock model and a model of the mass balance at the ice-atmosphere and ice-ocean interfaces. For both the Greenland and Antarctic ice sheets, calculations are made on a 10 km by 10 km resolution grid with 31 sigma levels. LOVECLIM1.2 reproduces well the major characteristics of the observed climate both for present-day conditions and for key past periods such as the last millennium, the mid-Holocene and the Last Glacial Maximum. However, despite some improvements compared to earlier versions, some biases are still present in the model. The most serious ones are mainly located at low latitudes with an overestimation of the temperature there, a too symmetric distribution of precipitation between the two hemispheres, and an overestimation of precipitation and vegetation cover in the subtropics. In addition, the atmospheric circulation is too weak. The model also tends to underestimate the surface temperature changes (mainly at low latitudes) and to overestimate the ocean heat uptake observed over the last decades.


2013 ◽  
Vol 10 (1) ◽  
pp. 339-355 ◽  
Author(s):  
P. B. Holden ◽  
N. R. Edwards ◽  
D. Gerten ◽  
S. Schaphoff

Abstract. We derive a constraint on the strength of CO2 fertilisation of the terrestrial biosphere through a "top-down" approach, calibrating Earth system model parameters constrained by the post-industrial increase of atmospheric CO2 concentration. We derive a probabilistic prediction for the globally averaged strength of CO2 fertilisation in nature, for the period 1850 to 2000 AD, implicitly net of other limiting factors such as nutrient availability. The approach yields an estimate that is independent of CO2 enrichment experiments. To achieve this, an essential requirement was the incorporation of a land use change (LUC) scheme into the GENIE Earth system model. Using output from a 671-member ensemble of transient GENIE simulations, we build an emulator of the change in atmospheric CO2 concentration change since the preindustrial period. We use this emulator to sample the 28-dimensional input parameter space. A Bayesian calibration of the emulator output suggests that the increase in gross primary productivity (GPP) in response to a doubling of CO2 from preindustrial values is very likely (90% confidence) to exceed 20%, with a most likely value of 40–60%. It is important to note that we do not represent all of the possible contributing mechanisms to the terrestrial sink. The missing processes are subsumed into our calibration of CO2 fertilisation, which therefore represents the combined effect of CO2 fertilisation and additional missing processes. If the missing processes are a net sink then our estimate represents an upper bound. We derive calibrated estimates of carbon fluxes that are consistent with existing estimates. The present-day land–atmosphere flux (1990–2000) is estimated at −0.7 GTC yr−1 (likely, 66% confidence, in the range 0.4 to −1.7 GTC yr−1). The present-day ocean–atmosphere flux (1990–2000) is estimated to be −2.3 GTC yr−1 (likely in the range −1.8 to −2.7 GTC yr−1). We estimate cumulative net land emissions over the post-industrial period (land use change emissions net of the CO2 fertilisation and climate sinks) to be 66 GTC, likely to lie in the range 0 to 128 GTC.


2010 ◽  
Vol 3 (1) ◽  
pp. 309-390 ◽  
Author(s):  
H. Goosse ◽  
V. Brovkin ◽  
T. Fichefet ◽  
R. Haarsma ◽  
P. Huybrechts ◽  
...  

Abstract. The main characteristics of the new version 1.2 of the three-dimensional Earth system model of intermediate complexity LOVECLIM are briefly described. LOVECLIM 1.2 includes representations of the atmosphere, the ocean and sea ice, the land surface (including vegetation), the ice sheets, the icebergs and the carbon cycle. The atmospheric component is ECBilt2, a T21, 3-level quasi-geostrophic model. The oceanic component is CLIO3, which is made up of an ocean general circulation model coupled to a comprehensive thermodynamic-dynamic sea-ice model. Its horizontal resolution is 3° by 3°, and there are 20 levels in the ocean. ECBilt-CLIO is coupled to VECODE, a vegetation model that simulates the dynamics of two main terrestrial plant functional types, trees and grasses, as well as desert. VECODE also simulates the evolution of the carbon cycle over land while the oceanic carbon cycle is represented in LOCH, a comprehensive model that takes into account both the solubility and biological pumps. The ice sheet component AGISM is made up of a three-dimensional thermomechanical model of the ice sheet flow, a visco-elastic bedrock model and a model of the mass balance at the ice-atmosphere and ice ocean interfaces. For both the Greenland and Antarctic ice sheets, calculations are made on a 10 km by 10 km resolution grid with 31 sigma levels. LOVECLIM 1.2 reproduces well the major characteristics of the observed climate both for present-day conditions and for key past periods such as the last millennium, the mid-Holocene and the Last Glacial Maximum. However, despite some improvements compared to earlier versions, some biases are still present in the model. The most serious ones are mainly located at low latitudes with an overestimation of the temperature there, a too symmetric distribution of precipitation between the two hemispheres, an overestimation of precipitation and vegetation cover in the subtropics. In addition, the atmospheric circulation is too weak. The model also tends to underestimate the surface temperature changes (mainly at low latitudes) and to overestimate the ocean heat uptake observed over the last decades.


2009 ◽  
Vol 34 (1) ◽  
pp. 151-151
Author(s):  
Marisa Montoya ◽  
Alexa Griesel ◽  
Anders Levermann ◽  
Juliette Mignot ◽  
Matthias Hofmann ◽  
...  

2009 ◽  
Vol 5 (2) ◽  
pp. 245-258 ◽  
Author(s):  
R. Calov ◽  
A. Ganopolski ◽  
C. Kubatzki ◽  
M. Claussen

Abstract. We investigate glacial inception and glacial thresholds in the climate-cryosphere system utilising the Earth system model of intermediate complexity CLIMBER-2, which includes modules for atmosphere, terrestrial vegetation, ocean and interactive ice sheets. The latter are described by the three-dimensional polythermal ice-sheet model SICOPOLIS. A bifurcation which represents glacial inception is analysed with two different model setups: one setup with dynamical ice-sheet model and another setup without it. The respective glacial thresholds differ in terms of maximum boreal summer insolation at 65° N (hereafter referred as Milankovitch forcing (MF)). The glacial threshold of the configuration without ice-sheet dynamics corresponds to a much lower value of MF compared to the full model. If MF attains values only slightly below the aforementioned threshold there is fast transient response. Depending on the value of MF relative to the glacial threshold, the transient response time of inland-ice volume in the model configuration with ice-sheet dynamics ranges from 10 000 to 100 000 years. Due to these long response times, a glacial threshold obtained in an equilibrium simulation is not directly applicable to the transient response of the climate-cryosphere system to time-dependent orbital forcing. It is demonstrated that in transient simulations just crossing of the glacial threshold does not imply large-scale glaciation of the Northern Hemisphere. We found that in transient simulations MF has to drop well below the glacial threshold determined in an equilibrium simulation to initiate glacial inception. Finally, we show that the asynchronous coupling between climate and inland-ice components allows one sufficient realistic simulation of glacial inception and, at the same time, a considerable reduction of computational costs.


2020 ◽  
Author(s):  
Christopher T. Reinhard ◽  
Stephanie Olson ◽  
Sandra Kirtland Turner ◽  
Cecily Pälike ◽  
Yoshiki Kanzaki ◽  
...  

Abstract. The methane (CH4) cycle is a key component of the Earth system that links planetary climate, biological metabolism, and the global biogeochemical cycles of carbon, oxygen, sulfur, and hydrogen. However, currently lacking is a numerical model capable of simulating a diversity of environments in the ocean where CH4 can be produced and destroyed, and with the flexibility to be able to explore not only relatively recent perturbations to Earth’s CH4 cycle but also to probe CH4 cycling and associated climate impacts under the very low-O2 conditions characteristic of most of Earth history and likely widespread on other Earth-like planets. Here, we present a refinement and expansion of the ocean-atmosphere CH4 cycle in the intermediate-complexity Earth system model cGENIE, including parameterized atmospheric O2-O3-CH4 photochemistry and schemes for microbial methanogenesis, aerobic methanotrophy, and anaerobic oxidation of methane (AOM). We describe the model framework, compare model parameterizations against modern observations, and illustrate the flexibility of the model through a series of example simulations. Though we make no attempt to rigorously tune default model parameters, we find that simulated atmospheric CH4 levels and marine dissolved CH4 distributions are generally in good agreement with empirical constraints for the modern and recent Earth. Finally, we illustrate the model’s utility in understanding the time-dependent behavior of the CH4 cycle resulting from transient carbon injection into the atmosphere, and present model ensembles that examine the effects of atmospheric pO2, oceanic dissolved SO42− and the thermodynamics of microbial metabolism on steady-state atmospheric CH4 abundance. Future model developments will address the sources and sinks of CH4 associated with the terrestrial biosphere and marine CH4 gas hydrates, both of which will be essential for comprehensive treatment of Earth’s CH4 cycle during geologically recent time periods.


Sign in / Sign up

Export Citation Format

Share Document