Using the climate model CLIMBER-X, we present an efficient method for assimilating the temporal evolution of surface temperatures for the last deglaciation covering the period 22000 to 6500 years ...before the present. The data assimilation methodology combines the data and the underlying dynamical principles governing the climate system to provide a state estimate of the system, which is better than that which could be obtained using just the data or the model alone. In applying an ensemble Kalman filter approach, we make use of the advances in the parallel data assimilation framework (PDAF), which provides parallel data assimilation functionality with a relatively small increase in computation time. We find that the data assimilation solution depends strongly on the background evolution of the decaying ice sheets rather than the assimilated temperatures. Two different ice sheet reconstructions result in a different deglacial meltwater history, affecting the large-scale ocean circulation and, consequently, the surface temperature. We find that the influence of data assimilation is more pronounced on regional scales than on the global mean. In particular, data assimilation has a stronger effect during millennial warming and cooling phases, such as the Bølling-Allerød and Younger Dryas, especially at high latitudes with heterogeneous temperature patterns. Our approach is a step toward a comprehensive paleo-reanalysis on multi-millennial time scales, including incorporating available paleoclimate data and accounting for their uncertainties in representing regional climates.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Abstract
With ongoing anthropogenic CO
2
emissions, the Greenland ice sheet approaches critical thresholds of inevitable, long-term mass loss. Future technologies might be able to efficiently remove ...CO
2
from the atmosphere and thereby cool down our planet. We explore whether and to what extent a realization of this concept could lead to a regrowth of the Greenland ice sheet once it has partly melted. Using the fully coupled Earth system model of intermediate complexity CLIMBER-X, emission pulses between 0 and 4000 GtC are released into the atmosphere, and 1 kyr, 2 kyr, and 5 kyr, the atmospheric CO
2
concentration is reduced back to its pre-industrial value. We find that independent of a specific trajectory, once the southern part of the Greenland ice sheet has partly melted with a total mass loss of more than 0.4 m sea level equivalent, regrowth is inhibited. Uncertainties preclude determination of precise thresholds, but model results indicate that cumulative industrial-era emissions approaching 1000 to 1500 GtC and beyond increasingly risk irreversible mass loss of the Greenland ice sheet. Once this threshold is passed, artificial atmospheric carbon removal would need to be utilised within the next centuries at massive scale. Beyond that, artificial atmospheric carbon removal has limited abilities to avoid long-term mass loss of the Greenland ice sheet. In conclusion, keeping cumulative anthropogenic emissions below 1000 to 1500 GtC is the only safe way to avoid irreversible mass loss of the Greenland ice sheet.
Understanding the future fate of the Greenland Ice Sheet (GIS) in the context of anthropogenic CO2 emissions is crucial to predict sea level rise. With the fully coupled Earth system model of ...intermediate complexity CLIMBER‐X, we study the stability of the GIS and its transient response to CO2 emissions over the next 10 Kyr. Bifurcation points exist at global temperature anomalies of 0.6 and 1.6 K relative to pre‐industrial. For system states in the vicinity of the equilibrium ice volumes corresponding to these temperature anomalies, mass loss rate and sensitivity of mass loss to cumulative CO2 emission peak. These critical ice volumes are crossed for cumulative emissions of 1,000 and 2,500 GtC, which would cause long‐term sea level rise by 1.8 and 6.9 m respectively. In summary, we find tipping of the GIS within the range of the temperature limits of the Paris agreement.
Plain Language Summary
With ongoing carbon dioxide emissions from the burning of fossil fuels, the atmosphere heats up, which has dramatic consequences for the ice sheets on Earth. In this study, we focus on the Greenland ice sheet (GIS), which holds so much ice that a complete melting would cause the global sea level to rise by 7 m. However, future mass loss of the GIS is challenging to predict because it is a non‐linear function of temperature and occurs over long timescales. For this reason, we use CLIMBER‐X, which is a coupled model of the whole Earth system. We find that the GIS features two critical volume thresholds, whose crossing would imply extensive further mass loss so that it would be difficult for the ice to grow back, even in thousands of years. Near these critical ice volumes, the mass loss rates are particularly high, and differences in the total carbon dioxide emission have a large impact. In summary, if cumulative emissions larger than 1,000 Gt carbon are released into the atmosphere, the GIS will shrink below a critical threshold and mass loss will inevitably continue until a substantial part of the ice sheet has melted.
Key Points
Bifurcation points exist at global mean temperature anomalies of 0.6 and 1.6 K relative to pre‐industrial
Mass loss rate and sensitivity to cumulative CO2 emission peak near the equilibrium ice volumes belonging to these temperature anomalies
Substantial long‐term mass loss of the Greenland ice sheet for cumulative emissions larger than 1,000 Gt carbon
Abrupt changes in climate, termed Dansgaard-Oeschger and Heinrich events, have punctuated the last glacial period (approximately 100-10 kyr ago) but not the Holocene (the past 10 kyr). Here we use an ...intermediate-complexity climate model to investigate the stability of glacial climate, and we find that only one mode of Atlantic Ocean circulation is stable: a cold mode with deep water formation in the Atlantic Ocean south of Iceland. However, a 'warm' circulation mode similar to the present-day Atlantic Ocean is only marginally unstable, and temporary transitions to this warm mode can easily be triggered. This leads to abrupt warm events in the model which share many characteristics of the observed Dansgaard-Oeschger events. For a large freshwater input (such as a large release of icebergs), the model's deep water formation is temporarily switched off, causing no strong cooling in Greenland but warming in Antarctica, as is observed for Heinrich events. Our stability analysis provides an explanation why glacial climate is much more variable than Holocene climate.
Celotno besedilo
Dostopno za:
DOBA, IJS, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
The Pliocene–Pleistocene Transition (PPT), from around 3.2 to 2.5 million years ago (Ma), represented a major shift in the climate system and was characterized by a gradual cooling trend and the ...appearance of large continental ice sheets over northern Eurasia and North America. Paleo evidence indicates that the PPT was accompanied and possibly caused by a decrease in atmospheric CO2, but the temporal resolution of CO2 reconstructions is low for this period of time and uncertainties remain large. Therefore, instead of applying existent CO2 reconstructions we solved an ‘inverse’ problem by finding a schematic CO2 concentration scenario that allows us to simulate the temporal evolution of key climate characteristics in agreement with paleoclimate records. To this end, we performed an ensemble of transient simulations with an Earth system model of intermediate complexity from which we derived a best guess transient CO2 scenario for the interval from 3.2 to 2.4 Ma that gives the best fit between the simulated and reconstructed benthic δ18O and global sea surface temperature evolution. Our data-constrained CO2 scenarios are consistent with recent CO2 reconstructions and suggest a gradual CO2 decline from 375–425 to 275–300 ppm, between 3.2 and 2.4 Ma. In addition to a gradual decline, the best fit to paleoclimate data requires the existence of pronounced CO2 variability coherent with the 41-kyr (1 kyr = 1000 years) obliquity cycle. In our simulations the long-term CO2 decline is accompanied by a relatively abrupt intensification of Northern Hemisphere glaciation at around 2.7 Ma. This is the result of a threshold behaviour of the ice sheets response to gradual CO2 decrease and orbital forcing. The simulated Northern Hemisphere ice sheets during the early Pleistocene glacial cycles reach a maximum volume equivalent to a sea level drop of about 40 m. Both ice volume and benthic δ18O are dominated by 41-kyr cyclicity. Our simulations suggest that before 2.7 Ma Greenland was ice free during summer insolation maxima and only partly ice covered during periods of minimum summer insolation. A fully glaciated Greenland comparable to its present-day ice volume is modelled only during glacial maxima after 2.7 Ma and more continuously after 2.5 Ma.
•Model-based approach to estimate CO2 over the Pliocene–Pleistocene transition.•We find a gradual CO2 decrease by about 100 ppm between 3.2 Ma and 2.4 Ma.•Model reproduces the intensification of Northern Hemisphere (NH) glaciation.•Threshold response of NH ice sheets to CO2 drop and orbital forcing.
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.
Glacial cycles of the late Quaternary are controlled by the asymmetrically varying mass balance of continental ice sheets in the Northern Hemisphere. Surface mass balance is governed by processes of ...ablation and accumulation. Here two ablation schemes, the positive-degree-day (PDD) method and the surface energy balance (SEB) approach, are compared in transient simulations of the last glacial cycle with the Earth system model of intermediate complexity CLIMBER-2. The standard version of the CLIMBER-2 model incorporates the SEB approach and simulates ice volume variations in reasonable agreement with paleoclimate reconstructions during the entire last glacial cycle. Using results from the standard CLIMBER-2 model version, we simulated ablation with the PDD method in offline mode by applying different combinations of three empirical parameters of the PDD scheme. We found that none of the parameter combinations allow us to simulate a surface mass balance of the American and European ice sheets that is similar to that obtained with the standard SEB method. The use of constant values for the empirical PDD parameters led either to too much ablation during the first phase of the last glacial cycle or too little ablation during the final phase. We then substituted the standard SEB scheme in CLIMBER-2 with the PDD scheme and performed a suite of fully interactive (online) simulations of the last glacial cycle with different combinations of PDD parameters. The results of these simulations confirmed the results of the offline simulations: no combination of PDD parameters realistically simulates the evolution of the ice sheets during the entire glacial cycle. The use of constant parameter values in the online simulations leads either to a buildup of too much ice volume at the end of glacial cycle or too little ice volume at the beginning. Even when the model correctly simulates global ice volume at the last glacial maximum (21 ka), it is unable to simulate complete deglaciation during the Holocene. According to our simulations, the SEB approach proves superior for simulations of glacial cycles.
Global mean surface temperatures are rising in response to anthropogenic greenhouse gas emissions. The magnitude of this warming at equilibrium for a given radiative forcing-referred to as specific ...equilibrium climate sensitivity (
)-is still subject to uncertainties. We estimate global mean temperature variations and
using a 784,000-year-long field reconstruction of sea surface temperatures and a transient paleoclimate model simulation. Our results reveal that
is strongly dependent on the climate background state, with significantly larger values attained during warm phases. Using the Representative Concentration Pathway 8.5 for future greenhouse radiative forcing, we find that the range of paleo-based estimates of Earth's future warming by 2100 CE overlaps with the upper range of climate simulations conducted as part of the Coupled Model Intercomparison Project Phase 5 (CMIP5). Furthermore, we find that within the 21st century, global mean temperatures will very likely exceed maximum levels reconstructed for the last 784,000 years. On the basis of temperature data from eight glacial cycles, our results provide an independent validation of the magnitude of current CMIP5 warming projections.
This paper presents a first assessment of the global effects of climate-induced sea-level rise on the erosion of sandy beaches, and its consequent impacts in the form of land loss and forced ...migration of people. We consider direct erosion on open sandy coasts and indirect erosion near selected tidal inlets and estuaries, using six global mean sea-level scenarios (in the range of 0.2–0.8m) and six SRES socio-economic development scenarios for the 21st century. Impacts are assessed both without and with adaptation in the form of shore and beach nourishment, based on cost-benefit analysis that includes the benefits of maintaining sandy beaches for tourism. Without nourishment, global land loss would amount to about 6000–17,000km2 during the 21st century, leading to 1.6–5.3million people being forced to migrate and migration costs of US$ 300–1000billion (not discounted). Optimal beach and shore nourishment would cost about US$ 65–220billion (not discounted) during the 21st century and would reduce land loss by 8–14%, forced migration by 56–68% and the cost of forced migration by 77–84% (not discounted). The global share of erodible coast that is nourished increases from about 4% in 2000 to 18–33% in 2100, with beach nourishment being 3–4 times more frequent than shore nourishment, reflecting the importance of tourism benefits. In absolute terms, with or without nourishment, large countries with long shorelines appear to have the largest costs, but in relative terms, small island states appear most impacted by erosion. Considerable uncertainty remains due to the limited availability of basic coastal geomorphological data and models on a global scale. Future work should also further explore the effects of beach tourism, including considering sub-national distributions of beach tourists.
•A first global analysis of erosion of sandy beaches due to sea-level rise•The analysis includes beach and shore nourishment as adaptation options.•6000–17,000km2 land may be lost during the 21st century.•1.6–5.3million people may be forced to migrate during the 21st century.•Nourishment based on benefit-cost analysis reduces forced migration by 56–68%.