Search

We present new high-resolution pollen records combined with palaeoceanographic proxies from the same samples in deep-sea cores SHAK06-5K and MD01-2444 on the southwestern Iberian Margin, documenting regional vegetation responses to orbital and millennial-scale climate changes over the last 28 ka. The chronology of these records is based on high-resolution radiocarbon dates of monospecific samples of the planktonic foraminifera Globigerina bulloides, measured from SHAK06-5K and MD01-2444 and aligned using an automated stratigraphical alignment method. Changes in temperate and steppe vegetation during Marine Isotope Stage 2 are closely coupled with sea surface temperature (SST) and global ice-volume changes. The peak expansion of thermophilous woodland between ~10.1 and 8.4 cal ka bp lags behind the boreal summer insolation maximum by ~2 ka, possibly arising from residual high-latitude ice-sheets into the Holocene. Rapid changes in pollen percentages are coeval with abrupt transitions in SSTs, precipitation and winter temperature at the onset and end of Heinrich Stadial 2, the ice-rafted debris event and end of Heinrich Stadial 1, and the onset of the Younger Dryas, suggesting extrinsically forced southwestern Iberian ecosystem changes by abrupt North Atlantic climate events. In contrast, the abrupt decline in thermophilous elements at ~7.8 cal ka bp indicates an intrinsically mediated abrupt vegetation response to the gradually declining boreal insolation, potentially resulting from the crossing of a seasonality of precipitation threshold.

The Australian Community Climate and Earth System Simulator (ACCESS) has contributed to the World Climate Research Programme’s Coupled Model Intercomparison Project Phase 6 (CMIP6) using two fully coupled model versions (ACCESS-CM2 and ACCESS-ESM1.5) and two ocean–sea-ice model versions (1° and 0.25° resolution versions of ACCESS-OM2). The fully coupled models differ primarily in the configuration and version of their atmosphere components (including the aerosol scheme), with smaller differences in their sea-ice and land model versions. Additionally, ACCESS-ESM1.5 includes biogeochemistry in the land and ocean components and can be run with an interactive carbon cycle. CMIP6 comprises core experiments and associated thematic Model Intercomparison Projects (MIPs). This paper provides an overview of the CMIP6 submission, including the methods used for the preparation of input forcing datasets and the post-processing of model output, along with a comprehensive list of experiments performed, detailing their initialisation, duration, ensemble number and computational cost. A small selection of model output is presented, focusing on idealised experiments and their variants at global scale. Differences in the climate simulation of the two coupled models are highlighted. ACCESS-CM2 produces a larger equilibrium climate sensitivity (4.7°C) than ACCESS-ESM1.5 (3.9°C), likely a result of updated atmospheric parameterisation in recent versions of the atmospheric component of ACCESS-CM2. The idealised experiments run with ACCESS-ESM1.5 show that land and ocean carbon fluxes respond to both changing atmospheric CO2 and to changing temperature. ACCESS data submitted to CMIP6 are available from the Earth System Grid Federation (https://doi.org/10.22033/ESGF/CMIP6.2281 and https://doi.org/10.22033/ESGF/CMIP6.2288). The information provided in this paper should facilitate easier use of these significant datasets by the broader climate community.

The Australian Community Climate and Earth System Simulator (ACCESS) has contributed to the World Climate Research Programme’s Coupled Model Intercomparison Project Phase 6 (CMIP6) using two fully coupled model versions (ACCESS-CM2 and ACCESS-ESM1.5) and two ocean–sea-ice model versions (1° and 0.25° resolution versions of ACCESS-OM2). The fully coupled models differ primarily in the configuration and version of their atmosphere components (including the aerosol scheme), with smaller differences in their sea-ice and land model versions. Additionally, ACCESS-ESM1.5 includes biogeochemistry in the land and ocean components and can be run with an interactive carbon cycle. CMIP6 comprises core experiments and associated thematic Model Intercomparison Projects (MIPs). This paper provides an overview of the CMIP6 submission, including the methods used for the preparation of input forcing datasets and the post-processing of model output, along with a comprehensive list of experiments performed, detailing their initialisation, duration, ensemble number and computational cost. A small selection of model output is presented, focusing on idealised experiments and their variants at global scale. Differences in the climate simulation of the two coupled models are highlighted. ACCESS-CM2 produces a larger equilibrium climate sensitivity (4.7°C) than ACCESS-ESM1.5 (3.9°C), likely a result of updated atmospheric parameterisation in recent versions of the atmospheric component of ACCESS-CM2. The idealised experiments run with ACCESS-ESM1.5 show that land and ocean carbon fluxes respond to both changing atmospheric CO2 and to changing temperature. ACCESS data submitted to CMIP6 are available from the Earth System Grid Federation (https://doi.org/10.22033/ESGF/CMIP6.2281 and https://doi.org/10.22033/ESGF/CMIP6.2288). The information provided in this paper should facilitate easier use of these significant datasets by the broader climate community.

Model intercomparison studies of coupled carbon-climate simulations have the potential to improve our understanding of the processes explaining the (Formula presented.) drawdown at the Last Glacial Maximum (LGM) and to identify related model biases. Models participating in the Paleoclimate Modeling Intercomparison Project (PMIP) now frequently include the carbon cycle. The ongoing PMIP-carbon project provides the first opportunity to conduct multimodel comparisons of simulated carbon content for the LGM time window. However, such a study remains challenging due to differing implementation of ocean boundary conditions (e.g., bathymetry and coastlines reflecting the low sea level) and to various associated adjustments of biogeochemical variables (i.e., alkalinity, nutrients, dissolved inorganic carbon). After assessing the ocean volume of PMIP models at the pre-industrial and LGM, we investigate the impact of these modeling choices on the simulated carbon at the global scale, using both PMIP-carbon model outputs and sensitivity tests with the iLOVECLIM model. We show that the carbon distribution in reservoirs is significantly affected by the choice of ocean boundary conditions in iLOVECLIM. In particular, our simulations demonstrate a (Formula presented.) GtC effect of an alkalinity adjustment on carbon sequestration in the ocean. Finally, we observe that PMIP-carbon models with a freely evolving (Formula presented.) and no additional glacial mechanisms do not simulate the (Formula presented.) drawdown at the LGM (with concentrations as high as 313, 331, and 315 ppm), especially if they use a low ocean volume. Our findings suggest that great care should be taken on accounting for large bathymetry changes in models including the carbon cycle.

The Antarctic Ice Sheet (AIS) is out of equilibrium with the current anthropogenic‐enhanced climate forcing. Paleoenvironmental records and ice sheet models reveal that the AIS has been tightly coupled to the climate system during the past and indicate the potential for accelerated and sustained Antarctic ice mass loss into the future. Modern observations by contrast suggest that the AIS has only just started to respond to climate change in recent decades. The maximum projected sea level contribution from Antarctica to 2100 has increased significantly since the Intergovernmental Panel on Climate Change (IPCC) 5th Assessment Report, although estimates continue to evolve with new observational and theoretical advances. This review brings together recent literature highlighting the progress made on the known processes and feedbacks that influence the stability of the AIS. Reducing the uncertainty in the magnitude and timing of the future sea level response to AIS change requires a multidisciplinary approach that integrates knowledge of the interactions between the ice sheet, solid Earth, atmosphere, and ocean systems and across time scales of days to millennia. We start by reviewing the processes affecting AIS mass change, from atmospheric and oceanic processes acting on short time scales (days to decades), through to ice processes acting on intermediate time scales (decades to centuries) and the response to solid Earth interactions over longer time scales (decades to millennia). We then review the evidence of AIS changes from the Pliocene to the present and consider the projections of global sea level rise and their consequences. We highlight priority research areas required to improve our understanding of the processes and feedbacks governing AIS change.

The Zero Emissions Commitment (ZEC) is the change in global mean temperature expected to occur following the cessation of net CO2 emissions and as such is a critical parameter for calculating the remaining carbon budget. The Zero Emissions Commitment Model Intercomparison Project (ZECMIP) was established to gain a better understanding of the potential magnitude and sign of ZEC, in addition to the processes that underlie this metric. A total of 18 Earth system models of both full and intermediate complexity participated in ZECMIP. All models conducted an experiment where atmospheric CO2 concentration increases exponentially until 1000 PgC has been emitted. Thereafter emissions are set to zero and models are configured to allow free evolution of atmospheric CO2 concentration. Many models conducted additional second-priority simulations with different cumulative emission totals and an alternative idealized emissions pathway with a gradual transition to zero emissions. The inter-model range of ZEC 50 years after emissions cease for the 1000 PgC experiment is −0.36 to 0.29 ∘C, with a model ensemble mean of −0.07 ∘C, median of −0.05 ∘C, and standard deviation of 0.19 ∘C. Models exhibit a wide variety of behaviours after emissions cease, with some models continuing to warm for decades to millennia and others cooling substantially. Analysis shows that both the carbon uptake by the ocean and the terrestrial biosphere are important for counteracting the warming effect from the reduction in ocean heat uptake in the decades after emissions cease. This warming effect is difficult to constrain due to high uncertainty in the efficacy of ocean heat uptake. Overall, the most likely value of ZEC on multi-decadal timescales is close to zero, consistent with previous model experiments and simple theory.

Irreversible shifts of large-scale components of the Earth system (so-called ‘tipping elements’) on policy-relevant timescales are a major source of uncertainty for projecting the impacts of future climate change. The high latitudes are particularly vulnerable to positive feedbacks that amplify change through atmosphere-ocean-ice interactions. Unfortunately, the short instrumental record does not capture the full range of past or projected climate scenarios (a situation particularly acute in the high latitudes). Natural archives from past periods warmer than present day, however, can be used to explore drivers and responses to forcing, and provide data against which to test models, thereby offering insights into the future. The Last Interglacial (129–116,000 years before present) — the warmest interglacial of the last 800,000 years — was the most recent period during which global temperatures were comparable with low-end 21st Century projections (up to 2 °C warmer, with temperature increase amplified over polar latitudes), providing a potentially useful analogue for future change. Substantial environmental changes happened during this time. Here we synthesise the nature and timing of potential high-latitude tipping elements during the Last Interglacial, including sea ice, extent of the boreal forest, permafrost, ocean circulation, and ice sheets/sea level. We also review the thresholds and feedbacks that likely operated through this period. Notably, substantial ice mass loss from Greenland, the West Antarctic, and possibly sectors of the East Antarctic drove a 6–9 m rise in global sea level. This was accompanied by reduced summer sea-ice extent, poleward-extended boreal forest, and reduced areas of permafrost. Despite current chronological uncertainties, we find that tipping elements in the high latitudes all experienced rapid and abrupt change (within 1–2 millennia of each other) across both hemispheres, while recovery to prior conditions took place over multi-mill

The Southern Ocean occupies 14% of the Earth’s surface and plays a fundamental role in the global carbon cycle and climate. It provides a direct connection to the deep ocean carbon reservoir through biogeochemical processes that include surface primary productivity, remineralization at depth and the upwelling of carbon-rich water masses. However, the role of these different processes in modulating past and future air–sea carbon flux remains poorly understood. A key period in this regard is the Antarctic Cold Reversal (ACR, 14.6–12.7 kyr bp), when mid- to high-latitude Southern Hemisphere cooling coincided with a sustained plateau in the global deglacial increase in atmospheric CO2. Here we reconstruct high-latitude Southern Ocean surface productivity from marine-derived aerosols captured in a highly resolved horizontal ice core. Our multiproxy reconstruction reveals a sustained signal of enhanced marine productivity across the ACR. Transient climate modelling indicates this period coincided with maximum seasonal variability in sea-ice extent, implying that sea-ice biological feedbacks enhanced CO2 sequestration and created a substantial regional marine carbon sink, which contributed to the plateau in CO2 during the ACR. Our results highlight the role Antarctic sea ice plays in controlling global CO2, and demonstrate the need to incorporate such feedbacks into climate–carbon models.

One of the major phases of sea level rise during the last deglaciation (∼19–11 thousand years before present [ka BP]) is Meltwater Pulse-1A (MWP-1A; ∼14.5 ka BP), when sea levels rose by 8.6 to 18 m in less than 400 years. Whether the meltwater originated from the partial disintegration of northern hemispheric ice sheets, from Antarctica, or both, remains controversial. Here we perform a series of idealized transient simulations of the last deglaciation, focusing on MWP-1A, with a three-dimensional oxygen-isotope enabled Earth System Climate Model. Three meltwater scenarios are considered during MWP-1A: a sole northern hemispheric source discharging into the North Atlantic, a sole Antarctic source, and a combined northern hemispheric-Antarctic source. A comparison of simulated changes in the oxygen-isotope composition (δ18O) of seawater and calcite with published marine sediment records points to a significant contribution from Antarctica. The best model-data fit is obtained with a contribution from both hemispheres. While the simulated changes over the 350 years of MWP-1A are overestimated in our simulations, the millennial-scale changes (∼14.6–13 ka BP) are underestimated, potentially alluding to a longer and sustained meltwater input over the whole period. Meltwater was not applied in the Arctic, the Gulf of Mexico, or the North Pacific in our simulations, and therefore, scenarios with meltwater originating from these regions cannot be excluded.

Benthic δ13C is often used to infer past changes in ocean circulation, though the interpretation of this proxy is difficult due to data scarcity and uncertainties. We present two methods for reconstructing the δ13C signal of North Atlantic Deep Water (NADW) and Antarctic Bottom Water and calculating the average oceanic δ13C values for the Atlantic Ocean based on δ13C from benthic foraminifera. The two simple statistical models are described and tested for the Holocene and the Last Glacial Maximum. The first statistical model consists of regressions of the δ13C data, which vary quadratically with depth and linearly with latitude. It differentiates between two regions, one for NADW and another for Antarctic Bottom Water. The second method consists of a hyperbolic tangent regression, which is bound asymptotically by the water mass source region averages (end-members). To test the robustness of the statistical models, two isotope-enabled climate models, the UVic ESCM and LOVECLIM, are sampled randomly, generating “pseudoproxies.” These are then used for testing the accuracy of the statistical models against the complete climate model δ13C outputs. We quantitatively compare the average δ13C and NADW depth against the original climate model outputs. We find that both statistical approaches are robust, regardless of the spatial distribution of the pseudoproxies, with the quadratic approach better able to capture the shape of NADW δ13C signal. Hence, this method can potentially be applied to different δ13C data sets to evaluate past changes in NADW.

The last glacial period was marked by pronounced millennial-scale variability in ocean circulation and global climate. Shifts in sea ice cover within the Nordic Seas are believed to have amplified the glacial climate variability in northern high latitudes and contributed to abrupt, high-amplitude temperature changes over Greenland. We present unprecedented empirical evidence that resolves the nature, timing, and role of sea ice fluctuations for abrupt ocean and climate change 32 to 40 thousand years ago, using biomarker sea ice reconstructions from the southern Norwegian Sea. Our results document that initial sea ice reductions at the core site preceded the major reinvigoration of convective deep-water formation in the Nordic Seas and abrupt Greenland warming; sea ice expansions preceded the buildup of a deep oceanic heat reservoir. Our findings suggest that the sea ice variability shaped regime shifts between surface stratification and deep convection in the Nordic Seas during abrupt climate changes.

During the Last Glacial Maximum (LGM; ~20,000 years ago), the global ocean sequestered a large amount of carbon lost from the atmosphere and terrestrial biosphere. Suppressed CO2 outgassing from the Southern Ocean is the prevailing explanation for this carbon sequestration. By contrast, the North Atlantic Ocean—a major conduit for atmospheric CO2 transport to the ocean interior via the overturning circulation—has received much less attention. Here we demonstrate that North Atlantic carbon pump efficiency during the LGM was almost doubled relative to the Holocene. This is based on a novel proxy approach to estimate air–sea CO2 exchange signals using combined carbonate ion and nutrient reconstructions for multiple sediment cores from the North Atlantic. Our data indicate that in tandem with Southern Ocean processes, enhanced North Atlantic CO2 absorption contributed to lowering ice-age atmospheric CO2.

Considerable ambiguity remains over the extent and nature of millennial/centennial-scale climate instability during the Last Interglacial (LIG). Here we analyse marine and terrestrial proxies from a deep-sea sediment sequence on the Portuguese Margin and combine results with an intensively dated Italian speleothem record and climate-model experiments. The strongest expression of climate variability occurred during the transitions into and out of the LIG. Our records also document a series of multi-centennial intra-interglacial arid events in southern Europe, coherent with cold water-mass expansions in the North Atlantic. The spatial and temporal fingerprints of these changes indicate a reorganization of ocean surface circulation, consistent with low-intensity disruptions of the Atlantic meridional overturning circulation (AMOC). The amplitude of this LIG variability is greater than that observed in Holocene records. Episodic Greenland ice melt and runoff as a result of excess warmth may have contributed to AMOC weakening and increased climate instability throughout the LIG.

Understanding feedbacks between the Greenland Ice Sheet (GrIS) and the Atlantic Meridional Overturning Circulation (AMOC) is crucial for reducing uncertainties over future sea level and ocean circulation change. Reconstructing past GrIS dynamics can extend the observational record and elucidate mechanisms that operate on multi-decadal timescales. We report a highly-constrained last glacial vertical profile of cosmogenic isotope exposure ages from Sermilik Fjord, a marine-terminating ice stream in the southeast sector of the GrIS. Our reconstruction reveals substantial ice-mass loss throughout the Younger Dryas (12.9-11.7 ka), a period of marked atmospheric and sea-surface cooling. Earth-system modelling reveals that southern GrIS marginal melt was likely driven by strengthening of the Irminger Current at depth due to a weakening of the AMOC during the Younger Dryas. This change in North Atlantic circulation appears to have drawn warm subsurface waters to southeast Greenland despite markedly cooler sea surface temperatures, enhancing thermal erosion at the grounding lines of palaeo ice-streams, supporting interpretation of regional marine-sediment cores. Given current rates of GrIS meltwater input into the North Atlantic and the vulnerability of major ice streams to water temperature changes at the grounding line, this mechanism has important implications for future AMOC changes and northern hemisphere heat transport.

The early part of the last deglaciation is characterised by a ~40 ppm atmospheric CO2 rise occurring in two abrupt phases. The underlying mechanisms driving these increases remain a subject of intense debate. Here, we successfully reproduce changes in CO2, δ 13C and Δ14C as recorded by paleo-records during Heinrich stadial 1 (HS1). We show that HS1 CO2 increase can be explained by enhanced Southern Ocean upwelling of carbon-rich Pacific deep and intermediate waters, resulting from intensified Southern Ocean convection and Southern Hemisphere (SH) westerlies. While enhanced Antarctic Bottom Water formation leads to a millennial CO2 outgassing, intensified SH westerlies induce a multi-decadal atmospheric CO2 rise. A strengthening of SH westerlies in a global eddy-permitting ocean model further supports a multi-decadal CO2 outgassing from the Southern Ocean. Our results highlight the crucial role of SH westerlies in the global climate and carbon cycle system with important implications for future climate projections.

Reconstructing the dynamic response of the Antarctic ice sheets to warming during the Last Glacial Termination (LGT; 18,000–11,650 yrs ago) allows us to disentangle ice-climate feedbacks that are key to improving future projections. Whilst the sequence of events during this period is reasonably well-known, relatively poor chronological control has precluded precise alignment of ice, atmospheric and marine records, making it difficult to assess relationships between Antarctic ice-sheet (AIS) dynamics, climate change and sea level. Here we present results from a highly-resolved ‘horizontal ice core’ from the Weddell Sea Embayment, which records millennial-scale AIS dynamics across this extensive region. Counterintuitively, we find AIS mass-loss across the full duration of the Antarctic Cold Reversal (ACR; 14,600–12,700 yrs ago), with stabilisation during the subsequent millennia of atmospheric warming. Earth-system and ice-sheet modelling suggests these contrasting trends were likely Antarctic-wide, sustained by feedbacks amplified by the delivery of Circumpolar Deep Water onto the continental shelf. Given the anti-phase relationship between inter-hemispheric climate trends across the LGT our findings demonstrate that Southern Ocean-AIS feedbacks were controlled by global atmospheric teleconnections. With increasing stratification of the Southern Ocean and intensification of mid-latitude westerly winds today, such teleconnections could amplify AIS mass loss and accelerate global sea-level rise.

© 2017 The Author(s). Solar variability has been hypothesized to be a major driver of North Atlantic millennial-scale climate variations through the Holocene along with orbitally induced insolation change. However, another important climate driver, volcanic forcing has generally been underestimated prior to the past 2,500 years partly owing to the lack of proper proxy temperature records. Here, we reconstruct seasonally unbiased and physically constrained Greenland Summit temperatures over the Holocene using argon and nitrogen isotopes within trapped air in a Greenland ice core (GISP2). We show that a series of volcanic eruptions through the Holocene played an important role in driving centennial to millennial-scale temperature changes in Greenland. The reconstructed Greenland temperature exhibits significant millennial correlations with K+ and Na+ ions in the GISP2 ice core (proxies for atmospheric circulation patterns), and δ18O of Oman and Chinese Dongge cave stalagmites (proxies for monsoon activity), indicating that the reconstructed temperature contains hemispheric signals. Climate model simulations forced with the volcanic forcing further suggest that a series of large volcanic eruptions induced hemispheric-wide centennial to millennial-scale variability through ocean/sea-ice feedbacks. Therefore, we conclude that volcanic activity played a critical role in driving centennial to millennial-scale Holocene temperature variability in Greenland and likely beyond.