Search Results
Role of CO2, climate and land use in regulating the seasonal amplitude increase of carbon fluxes in terrestrial ecosystems: A multimodel analysis
We examined the net terrestrial carbon flux to the atmosphere (FTA) simulated by nine models from the TRENDY dynamic global vegetation model project for its seasonal cycle and amplitude trend during 1961–2012. While some models exhibit similar phase and amplitude compared to atmospheric inversions, with spring drawdown and autumn rebound, others tend to rebound early in summer. The model ensemble mean underestimates the magnitude of the seasonal cycle by 40 % compared to atmospheric inversions. Global FTA amplitude increase (19 ± 8 %) and its decadal variability from the model ensemble are generally consistent with constraints from surface atmosphere observations. However, models disagree on attribution of this long-term amplitude increase, with factorial experiments attributing 83 ± 56 %, −3 ± 74 and 20 ± 30 % to rising CO2, climate change and land use/cover change, respectively. Seven out of the nine models suggest that CO2 fertilization is the strongest control – with the notable exception of VEGAS, which attributes approximately equally to the three factors. Generally, all models display an enhanced seasonality over the boreal region in response to high-latitude warming, but a negative climate contribution from part of the Northern Hemisphere temperate region, and the net result is a divergence over climate change effect. Six of the nine models show that land use/cover change amplifies the seasonal cycle of global FTA: some are due to forest regrowth, while others are caused by crop expansion or agricultural intensification, as revealed by their divergent spatial patterns. We also discovered a moderate cross-model correlation between FTA amplitude increase and increase in land carbon sink (R2 =  0.61). Our results suggest that models can show similar results in some benchmarks with different underlying mechanisms; therefore, the spatial traits of CO2 fertilization, climate change and land use/cover changes are crucial in determining the right mechanisms in seasonal carbon cycle change as well as mean sink change.
The importance of snow albedo for ice sheet evolution over the last glacial cycle
The surface energy and mass balance of ice sheets strongly depends on the amount of solar radiation absorbed at the surface, which is mainly controlled by the albedo of snow and ice. Here, using an Earth system model of intermediate complexity, we explore the role played by surface albedo for the simulation of glacial cycles. We show that the evolution of the Northern Hemisphere ice sheets over the last glacial cycle is very sensitive to the representation of snow albedo in the model. It is well known that the albedo of snow depends strongly on snow grain size and the content of light-absorbing impurities. Excluding either the snow aging effect or the dust darkening effect on snow albedo leads to an excessive ice build-up during glacial times and consequently to a failure in simulating deglaciation. While the effect of snow grain growth on snow albedo is well constrained, the albedo reduction due to the presence of dust in snow is much more uncertain because the light-absorbing properties of dust vary widely as a function of dust mineral composition. We also show that assuming slightly different optical properties of dust leads to very different ice sheet and climate evolutions in the model. Conversely, ice sheet evolution is less sensitive to the choice of ice albedo in the model. We conclude that a proper representation of snow albedo is a fundamental prerequisite for a successful simulation of glacial cycles.
Simulation of climate, ice sheets and CO2 evolution during the last four glacial cycles with an Earth system model of intermediate complexity
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.
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.
Prominent role of volcanism in Common Era climate variability and human history
Climate reconstructions for the Common Era are compromised by the paucity of annually-resolved and absolutely-dated proxy records prior to medieval times. Where reconstructions are based on combinations of different climate archive types (of varying spatiotemporal resolution, dating uncertainty, record length and predictive skill), it is challenging to estimate past amplitude ranges, disentangle the relative roles of natural and anthropogenic forcing, or probe deeper interrelationships between climate variability and human history. Here, we compile and analyse updated versions of all the existing summer temperature sensitive tree-ring width chronologies from the Northern Hemisphere that span the entire Common Era. We apply a novel ensemble approach to reconstruct extra-tropical summer temperatures from 1 to 2010 CE, and calculate uncertainties at continental to hemispheric scales. Peak warming in the 280s, 990s and 1020s, when volcanic forcing was low, was comparable to modern conditions until 2010 CE. The lowest June–August temperature anomaly in 536 not only marks the beginning of the coldest decade, but also defines the onset of the Late Antique Little Ice Age (LALIA). While prolonged warmth during Roman and medieval times roughly coincides with the tendency towards societal prosperity across much of the North Atlantic/European sector and East Asia, major episodes of volcanically-forced summer cooling often presaged widespread famines, plague outbreaks and political upheavals. Our study reveals a larger amplitude of spatially synchronized summer temperature variation during the first millennium of the Common Era than previously recognised.
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.
Mechanisms and time scales of glacial inception simulated with an Earth system model of intermediate complexity
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.
Geoengineering climate by stratospheric sulfur injections: Earth system vulnerability to technological failure
We use a coupled climate-carbon cycle model of intermediate complexity to investigate scenarios of stratospheric sulfur injections as a measure to compensate for CO2-induced global warming. The baseline scenario includes the burning of 5,000 GtC of fossil fuels. A full compensation of CO2-induced warming requires a load of about 13 MtS in the stratosphere at the peak of atmospheric CO2 concentration. Keeping global warming below 2°C reduces this load to 9 MtS. Compensation of CO 2 forcing by stratospheric aerosols leads to a global reduction in precipitation, warmer winters in the high northern latitudes and cooler summers over northern hemisphere landmasses. The average surface ocean pH decreases by 0.7, reducing the calcifying ability of marine organisms. Because of the millennial persistence of the fossil fuel CO2 in the atmosphere, high levels of stratospheric aerosol loading would have to continue for thousands of years until CO2 was removed from the atmosphere. A termination of stratospheric aerosol loading results in abrupt global warming of up to 5°C within several decades, a vulnerability of the Earth system to technological failure. © 2008 The Author(s).
Historical greenhouse gas concentrations for climate modelling (CMIP6)
Atmospheric greenhouse gas (GHG) concentrations are at unprecedented, record-high levels compared to the last 800000 years. Those elevated GHG concentrations warm the planet and – partially offset by net cooling effects by aerosols – are largely responsible for the observed warming over the past 150 years. An accurate representation of GHG concentrations is hence important to understand and model recent climate change. So far, community efforts to create composite datasets of GHG concentrations with seasonal and latitudinal information have focused on marine boundary layer conditions and recent trends since the 1980s. Here, we provide consolidated datasets of historical atmospheric concentrations (mole fractions) of 43 GHGs to be used in the Climate Model Intercomparison Project – Phase 6 (CMIP6) experiments. The presented datasets are based on AGAGE and NOAA networks, firn and ice core data, and archived air data, and a large set of published studies. In contrast to previous intercomparisons, the new datasets are latitudinally resolved and include seasonality. We focus on the period 1850–2014 for historical CMIP6 runs, but data are also provided for the last 2000 years. We provide consolidated datasets in various spatiotemporal resolutions for carbon dioxide (CO2), methane (CH4) and nitrous oxide (N2O), as well as 40 other GHGs, namely 17 ozone-depleting substances, 11 hydrofluorocarbons (HFCs), 9 perfluorocarbons (PFCs), sulfur hexafluoride (SF6), nitrogen trifluoride (NF3) and sulfuryl fluoride (SO2F2). In addition, we provide three equivalence species that aggregate concentrations of GHGs other than CO2, CH4 and N2O, weighted by their radiative forcing efficiencies. For the year 1850, which is used for pre-industrial control runs, we estimate annual global-mean surface concentrations of CO2 at 284.3ppm, CH4 at 808.2ppb and N2O at 273.0ppb. The data are available at https://esgf-node.llnl.gov/search/input4mips/ and http://www.climatecollege.unimelb.edu.au/cmip6. While the minimum CMIP6 recommendation is to use the global- and annual-mean time series, modelling groups can also choose our monthly and latitudinally resolved concentrations, which imply a stronger radiative forcing in the Northern Hemisphere winter (due to the latitudinal gradient and seasonality).
Significance of uncertain phasing between the onsets of stadial–interstadial transitions in different Greenland ice core proxies
Different paleoclimate proxy records evidence repeated abrupt climate transitions during previous glacial intervals. These transitions are thought to comprise abrupt warming and increase in local precipitation over Greenland, sudden reorganization of the Northern Hemisphere atmospheric circulation, and retreat of sea ice in the North Atlantic. The physical mechanism underlying these so-called Dansgaard–Oeschger (DO) events remains debated. A recent analysis of Greenland ice core proxy records found that transitions in Na+ concentrations and δ18O values are delayed by about 1 decade with respect to corresponding transitions in Ca2+ concentrations and in the annual layer thickness during DO events. These delays are interpreted as a temporal lag of sea-ice retreat and Greenland warming with respect to a synoptic- and hemispheric-scale atmospheric reorganization at the onset of DO events and may thereby help constrain possible triggering mechanisms for the DO events. However, the explanatory power of these results is limited by the uncertainty of the transition onset detection in noisy proxy records. Here, we extend previous work by testing the significance of the reported lags with respect to the null hypothesis that the proposed transition order is in fact not systematically favored. If the detection uncertainties are averaged out, the temporal delays in the δ18O and Na+ transitions with respect to their counterparts in Ca2+ and the annual layer thickness are indeed pairwise statistically significant. In contrast, under rigorous propagation of uncertainty, three statistical tests cannot provide evidence against the null hypothesis. We thus confirm the previously reported tendency of delayed transitions in the δ18O and Na+ concentration records. Yet, given the uncertainties in the determination of the transition onsets, it cannot be decided whether these tendencies are truly the imprint of a prescribed transition order or whether they are due to chance. The analyzed set of DO transitions can therefore not serve as evidence for systematic lead–lag relationships between the transitions in the different proxies, which in turn limits the power of the observed tendencies to constrain possible physical causes of the DO events.