Search Results
On the effect of orbital forcing on mid-Pliocene climate, vegetation and ice sheets
We present results from modelling of the mid-Pliocene warm period (3.3–3 million years ago) using the Earth system model of intermediate complexity CLIMBER-2 analysing the effect of changes in boundary conditions as well as of orbital forcing on climate. First we performed equilibrium experiments following the PlioMIP (Pliocene Model Intercomparison Project) protocol with a CO2 concentration of 405 ppm, reconstructed mid-Pliocene orography and vegetation and a present-day orbital configuration. Simulated global Pliocene warming is about 2.5 °C, fully consistent with results of atmosphere–ocean general circulation model simulations performed for the same modelling setup. A factor separation analysis attributes 1.5 °C warming to CO2, 0.3 °C to orography, 0.2 °C to ice sheets and 0.4 °C to vegetation. Transient simulations for the entire mid-Pliocene warm period with time-dependent orbital forcing as well as interactive ice sheets and vegetation give a global warming varying within the range 1.9–2.8 °C. Ice sheet and vegetation feedbacks in synergy act as amplifiers of the orbital forcing, transforming seasonal insolation variations into an annual mean temperature signal. The effect of orbital forcing is more significant at high latitudes, especially during boreal summer, when the warming over land varies in the wide range from 0 to 10 °C. The modelled ice-sheet extent and vegetation distribution also show significant temporal variations. Modelled and reconstructed data for Northern Hemisphere sea-surface temperatures and vegetation distribution show the best agreement if the reconstructions are assumed to be representative for the warmest periods during the orbital cycles. This suggests that low-resolution Pliocene palaeoclimate reconstructions can reflect not only the impact of increased CO2 concentrations and topography changes but also the effect of orbital forcing. Therefore, the climate (Earth system) sensitivity estimates from Pliocene reconstructions which do not account for the effect of orbital forcing can be biased toward high values.
Heinrich event 1: An example of dynamical ice-sheet reaction to oceanic changes
Heinrich events, identified as enhanced ice-rafted detritus (IRD) in North Atlantic deep sea sediments (Heinrich, 1988; Hemming, 2004) have classically been attributed to Laurentide ice-sheet (LIS) instabilities (MacAyeal, 1993; Calov et al., 2002; Hulbe et al., 2004) and assumed to lead to important disruptions of the Atlantic meridional overturning circulation (AMOC) and North Atlantic deep water (NADW) formation. However, recent paleoclimate data have revealed that most of these events probably occurred after the AMOC had already slowed down or/and NADW largely collapsed, within about a thousand years (Hall et al., 2006; Hemming, 2004; Jonkers et al., 2010; Roche et al., 2004), implying that the initial AMOC reduction could not have been caused by the Heinrich events themselves. Here we propose an alternative driving mechanism, specifically for Heinrich event 1 (H1; 18 to 15 ka BP), by which North Atlantic ocean circulation changes are found to have strong impacts on LIS dynamics. By combining simulations with a coupled climate model and a three-dimensional ice sheet model, our study illustrates how reduced NADW and AMOC weakening lead to a subsurface warming in the Nordic and Labrador Seas resulting in rapid melting of the Hudson Strait and Labrador ice shelves. Lack of buttressing by the ice shelves implies a substantial ice-stream acceleration, enhanced ice-discharge and sea level rise, with peak values 500–1500 yr after the initial AMOC reduction. Our scenario modifies the previous paradigm of H1 by solving the paradox of its occurrence during a cold surface period, and highlights the importance of taking into account the effects of oceanic circulation on ice-sheets dynamics in order to elucidate the triggering mechanism of Heinrich events.
The role of orbital forcing, carbon dioxide and regolith in 100 kyr glacial cycles
The origin of the 100 kyr cyclicity, which dominates ice volume variations and other climate records over the past million years, remains debatable. Here, using a comprehensive Earth system model of intermediate complexity, we demonstrate that both strong 100 kyr periodicity in the ice volume variations and the timing of glacial terminations during past 800 kyr can be successfully simulated as direct, strongly nonlinear responses of the climate-cryosphere system to orbital forcing alone, if the atmospheric CO2 concentration stays below its typical interglacial value. The existence of long glacial cycles is primarily attributed to the North American ice sheet and requires the presence of a large continental area with exposed rocks. We show that the sharp, 100 kyr peak in the power spectrum of ice volume results from the long glacial cycles being synchronized with the Earth's orbital eccentricity. Although 100 kyr cyclicity can be simulated with a constant CO2 concentration, temporal variability in the CO2 concentration plays an important role in the amplification of the 100 kyr cycles.
Coupled Northern Hemisphere permafrost-ice-sheet evolution over the last glacial cycle
Permafrost influences a number of processes which are relevant for local and global climate. For example, it is well known that permafrost plays an important role in global carbon and methane cycles. Less is known about the interaction between permafrost and ice sheets. In this study a permafrost module is included in the Earth system model CLIMBER-2, and the coupled Northern Hemisphere (NH) permafrost–ice-sheet evolution over the last glacial cycle is explored. The model performs generally well at reproducing present-day permafrost extent and thickness. Modeled permafrost thickness is sensitive to the values of ground porosity, thermal conductivity and geothermal heat flux. Permafrost extent at the Last Glacial Maximum (LGM) agrees well with reconstructions and previous modeling estimates. Present-day permafrost thickness is far from equilibrium over deep permafrost regions. Over central Siberia and the Arctic Archipelago permafrost is presently up to 200–500 m thicker than it would be at equilibrium. In these areas, present-day permafrost depth strongly depends on the past climate history and simulations indicate that deep permafrost has a memory of surface temperature variations going back to at least 800 ka. Over the last glacial cycle permafrost has a relatively modest impact on simulated NH ice sheet volume except at LGM, when including permafrost increases ice volume by about 15 m sea level equivalent in our model. This is explained by a delayed melting of the ice base from below by the geothermal heat flux when the ice sheet sits on a porous sediment layer and permafrost has to be melted first. Permafrost affects ice sheet dynamics only when ice extends over areas covered by thick sediments, which is the case at LGM.
The modern and glacial overturning circulation in the Atlantic ocean in PMIP coupled model simulations
This study analyses the response of the Atlantic meridional overturning circulation (AMOC) to LGM forcings and boundary conditions in nine PMIP coupled model simulations, including both GCMs and Earth system Models of Intermediate Complexity. Model results differ widely. The AMOC slows down considerably (by 20–40%) during the LGM as compared to the modern climate in four models, there is a slight reduction in one model and four models show a substantial increase in AMOC strength (by 10–40%). It is found that a major controlling factor for the AMOC response is the density contrast between Antarctic Bottom Water (AABW) and North Atlantic Deep Water (NADW) at their source regions. Changes in the density contrast are determined by the opposing effects of changes in temperature and salinity, with more saline AABW as compared to NADW consistently found in all models and less cooling of AABW in all models but one. In only two models is the AMOC response during the LGM directly related to the response in net evaporation over the Atlantic basin. Most models show large changes in the ocean freshwater transports into the basin, but this does not seem to affect the AMOC response. Finally, there is some dependence on the accuracy of the control state.
Simulation of the last glacial cycle with a coupled climate ice-sheet model of intermediate complexity
A new version of the Earth system model of intermediate complexity, CLIMBER-2, which includes the three-dimensional polythermal ice-sheet model SICOPOLIS, is used to simulate the last glacial cycle forced by variations of the Earth's orbital parameters and atmospheric concentration of major greenhouse gases. The climate and ice-sheet components of the model are coupled bi-directionally through a physically-based surface energy and mass balance interface. The model accounts for the time-dependent effect of aeolian dust on planetary and snow albedo. The model successfully simulates the temporal and spatial dynamics of the major Northern Hemisphere (NH) ice sheets, including rapid glacial inception and strong asymmetry between the ice-sheet growth phase and glacial termination. Spatial extent and elevation of the ice sheets during the last glacial maximum agree reasonably well with palaeoclimate reconstructions. A suite of sensitivity experiments demonstrates that simulated ice-sheet evolution during the last glacial cycle is very sensitive to some parameters of the surface energy and mass-balance interface and dust module. The possibility of a considerable acceleration of the climate ice-sheet model is discussed.
MIS-11 duration key to disappearance of the Greenland ice sheet
Palaeo data suggest that Greenland must have been largely ice free during Marine Isotope Stage 11 (MIS-11). However, regional summer insolation anomalies were modest during this time compared to MIS-5e, when the Greenland ice sheet likely lost less volume. Thus it remains unclear how such conditions led to an almost complete disappearance of the ice sheet. Here we use transient climate-ice sheet simulations to simultaneously constrain estimates of regional temperature anomalies and Greenland's contribution to the MIS-11 sea-level highstand. We find that Greenland contributed 6.1 m (3.9-7.0 m, 95% credible interval) to sea level, ∼7 kyr after the peak in regional summer temperature anomalies of 2.8 °C (2.1-3.4 °C). The moderate warming produced a mean rate of mass loss in sea-level equivalent of only around 0.4 m per kyr, which means the long duration of MIS-11 interglacial conditions around Greenland was a necessary condition for the ice sheet to disappear almost completely.
Modeling sensitivity study of the possible impact of snow and glaciers developing over Tibetan Plateau on Holocene African-Asian summer monsoon climate
The impacts of various scenarios of a gradual snow and glaciers developing over the Tibetan Plateau on climate change in Afro-Asian monsoon region and other regions during the Holocene (9 kyr BP–0 kyr BP) are studied by using the Earth system model of intermediate complexity, CLIMBER-2. The simulations show that the imposed snow and glaciers over the Tibetan Plateau in the mid-Holocene induce global summer temperature decreases over most of Eurasia but in the Southern Asia temperature response is opposite. With the imposed snow and glaciers, summer precipitation decreases strongly in North Africa and South Asia as well as northeastern China, while it increases in Southeast Asia and the Mediterranean. For the whole period of Holocene (9 kyr BP–0 kyr BP), the response of vegetation cover to the imposed snow and glaciers cover over the Tibetan Plateau is not synchronous in South Asia and in North Africa, showing an earlier and a more rapid decrease in vegetation cover in North Africa from 9 kyr BP to 6 kyr BP while it has only minor influence on that in South Asia until 5 kyr BP. The precipitation decreases rapidly in North Africa and South Asia while it decreases slowly or unchanged during 6 kyr BP to 0 kyr BP with imposed snow and glacier cover over the Tibetan Plateau. The different scenarios of snow and glacier developing over the Tibetan Plateau would result in differences in variation of temperature, precipitation and vegetation cover in North Africa, South Asia and Southeast Asia. The model results suggest that the development of snow and ice cover over Tibetan Plateau represents an additional important climate feedback, which amplify orbital forcing and produces a significant synergy with the positive vegetation feedback.
Glacial CO 2 cycle as a succession of key physical and biogeochemical processes
During glacial-interglacial cycles, atmospheric CO2 concentration varied by about 100 ppmv in amplitude. While testing mechanisms that have led to the low glacial CO2 level could be done in equilibrium model experiments, an ultimate goal is to explain CO2 changes in transient simulations through the complete glacial-interglacial cycle. The computationally efficient Earth System model of intermediate complexity CLIMBER-2 is used to simulate global biogeochemistry over the last glacial cycle (126 kyr). The physical core of the model (atmosphere, ocean, land and ice sheets) is driven by orbital changes and reconstructed radiative forcing from greenhouses gases, ice, and aeolian dust. The carbon cycle model is able to reproduce the main features of the CO2 changes: a 50 ppmv CO2 drop during glacial inception, a minimum concentration at the last glacial maximum 80 ppmv lower than the Holocene value, and an abrupt 60 ppmv CO2 rise during the deglaciation. The model deep ocean δ13C also resembles reconstructions from deep-sea cores. The main drivers of atmospheric CO2 evolve in time: changes in sea surface temperatures and in the volume of bottom water of southern origin control atmospheric CO2 during the glacial inception and deglaciation; changes in carbonate chemistry and marine biology are dominant during the first and second parts of the glacial cycle, respectively. These feedback mechanisms could also significantly impact the ultimate climate response to the anthropogenic perturbation.