New Research In
Physical Sciences
Social Sciences
Featured Portals
Articles by Topic
Biological Sciences
Featured Portals
Articles by Topic
- Agricultural Sciences
- Anthropology
- Applied Biological Sciences
- Biochemistry
- Biophysics and Computational Biology
- Cell Biology
- Developmental Biology
- Ecology
- Environmental Sciences
- Evolution
- Genetics
- Immunology and Inflammation
- Medical Sciences
- Microbiology
- Neuroscience
- Pharmacology
- Physiology
- Plant Biology
- Population Biology
- Psychological and Cognitive Sciences
- Sustainability Science
- Systems Biology
The multimillennial sea-level commitment of global warming
Edited by John C. Moore, College of Global Change and Earth System Science, Beijing, China, and accepted by the Editorial Board June 13, 2013 (received for review November 7, 2012)

Abstract
Global mean sea level has been steadily rising over the last century, is projected to increase by the end of this century, and will continue to rise beyond the year 2100 unless the current global mean temperature trend is reversed. Inertia in the climate and global carbon system, however, causes the global mean temperature to decline slowly even after greenhouse gas emissions have ceased, raising the question of how much sea-level commitment is expected for different levels of global mean temperature increase above preindustrial levels. Although sea-level rise over the last century has been dominated by ocean warming and loss of glaciers, the sensitivity suggested from records of past sea levels indicates important contributions should also be expected from the Greenland and Antarctic Ice Sheets. Uncertainties in the paleo-reconstructions, however, necessitate additional strategies to better constrain the sea-level commitment. Here we combine paleo-evidence with simulations from physical models to estimate the future sea-level commitment on a multimillennial time scale and compute associated regional sea-level patterns. Oceanic thermal expansion and the Antarctic Ice Sheet contribute quasi-linearly, with 0.4 m °C−1 and 1.2 m °C−1 of warming, respectively. The saturation of the contribution from glaciers is overcompensated by the nonlinear response of the Greenland Ice Sheet. As a consequence we are committed to a sea-level rise of approximately 2.3 m °C−1 within the next 2,000 y. Considering the lifetime of anthropogenic greenhouse gases, this imposes the need for fundamental adaptation strategies on multicentennial time scales.
Sea-level projections show a robust, albeit highly uncertain, increase by the end of this century (1, 2), and there is strong evidence that sea level will continue to rise beyond the year 2100 unless the current global mean temperature trend is reversed (3⇓⇓–6). At the same time, inertia in the climate and global carbon system causes the global mean temperature to decline slowly even after greenhouse gas emissions have ceased (6), raising the question of how much sea-level rise we are committed to on a multimillennial time scale for different levels of global mean temperature increase. During the 20th century, sea level rose by approximately 0.2 m (7, 8), and it is estimated to rise by significantly less than 2 m by 2100, even for the strongest scenarios considered (9). At the same time, past climate records suggest a sea-level sensitivity of as much as several meters per degree of warming during previous intervals of Earth history when global temperatures were similar to or warmer than present (10, 11). Although sea-level rise over the last century has been dominated by ocean warming and loss of glaciers (7), the sensitivity suggested from records of past sea level indicates important contributions from the Greenland and Antarctic Ice Sheets. Because of the uncertainties in the paleo-reconstructions, however, additional strategies are required to better constrain the sea-level commitment. Here we describe the models used and the resulting estimates of long-term sea-level rise from each component of the Earth system. We combine simulations from process-based physical models for the four main components that contribute to sea-level changes to give a robust estimate of the sea-level commitment on multimillennial time scales up to a global mean temperature increase of 4 °C. Our results are then compared with paleo-evidence, with the good agreement providing an independent validation of our modeling results.
Modeled Components of Sea Level
Thermal Expansion.
The thermal expansion of the ocean has been investigated by a spectrum of climate models of different complexity, ranging from zero-dimensional diffusion models (12, 13) via Earth System Models of Intermediate Complexity (EMIC) (6, 14) to comprehensive general circulation models (15, 16). Although uncertainty remains, especially owing to uncertainty in the ocean circulation and thereby the distribution of heat within the ocean, the physical processes are relatively well understood even if not fully represented in all models. On multimillennial time scales as applied here, the application of comprehensive climate models is not feasible because of the required computational effort. Because the general processes responsible for oceanic expansion are, however, also integrated into lower resolution ocean models as used in EMICs, the range of long-term thermal expansion is likely to be covered by these models.
We take the thermal expansion of the ocean on multimillennial time scales from 10,000-y integrations with six coupled climate models. The results, which were used in the Fourth Assessment Report of the Intergovernmental Panel on Climate Change (figure 10.34 in ref. 1), yield a rate of sea-level change in the range of 0.20–0.63 m °C−1 (Fig. 1A). For reference, a homogenous increase of ocean temperature by 1 °C would yield a global mean sea-level rise of 0.38 m when added to reanalysis data (17). Uncertainty arises owing to the different spatial distribution of the warming in models and the dependence of the expansion on local temperature and salinity.
Sea-level commitment per degree of warming as obtained from physical model simulations of (A) ocean warming, (B) mountain glaciers and ice caps, and (C) the Greenland and (D) the Antarctic Ice Sheets. (E) The corresponding total sea-level commitment, which is consistent with paleo-estimates from past warm periods (PI, preindustrial, Plio, mid-Pliocene; see text for discussion). Temperatures are relative to preindustrial. Dashed lines and large dots provide linear approximations: (A) sea-level rise for a spatially homogeneous increase in ocean temperature; (A, D, E) constant slopes of 0.42, 1.2, and 1.8 and 2.3 m/°C. Shading as well as boxes represent the uncertainty range as discussed in the text. (A–C) Thin lines provide the individual simulation results from different models (A and B) or different parameter combinations (C). The small black dots in D represent 1,000-y averages of the 5-million-year simulation of Antarctica following ref. 36.
Glaciers.
A number of different approaches have been used to estimate the contribution from glaciers to global sea level for the 21st century, including extrapolations of either current rates of mass loss or current acceleration of rates of mass loss (18), estimates based on either constant or constantly declining accumulation area ratios in the future (19), and estimates based on explicit modeling of future glacier mass balances (20, 21). Although these approaches differ to a certain extent in the outcome for the 21st century, the total sea-level contribution from glaciers is limited by their total ice volume, which is small compared with that of the ice sheets on Greenland and Antarctica. Glaciers thus have a minor role on multimillennial time scales compared with the large contribution from oceanic expansion and the ice sheets. To our knowledge, there are no published global estimates of glacier contribution to sea-level change on time scales longer than a few centuries, or of the equilibrium response of glaciers to climate change.
The total possible contribution of glaciers (i.e., all of the land ice excluding the ice sheets) is limited to ∼0.6 m (22). We use two models (20, 21) to compute the long-term contribution for different levels of global mean temperature. Both models couple surface mass balance with simplified ice-dynamics models and are forced by temperature and precipitation scenarios for each glacier in the world. Radić and Hock’s model (20) is forced by modified monthly temperature and precipitation series for 2001–2300 from four general circulation models (GCMs) of the coupled model intercomparison project CMIP-3 (U.K.MO-HadCM3, ECHAM5/MPI-OM, GFDL-CM2.0, and CSIRO-Mk3.0), using the A1B emission scenario, whereas the model by Marzeion et al. (21) is forced by temperature and precipitation anomalies from 15 GCMs of CMIP-5, using the representative concentration pathway (RCP)-8.5.
To obtain an estimate of the sea-level contribution on long time scales, temperature and precipitation patterns were kept constant at different levels of global warming. The uncertainty range is obtained as the model spread across the 19 different climate forcing and glacier model combinations. The resulting sensitivity of sea-level commitment decreases for increasing temperatures from 0.21 m °C−1 at preindustrial temperature levels to 0.04 m °C−1 at 4 °C of warming (Fig. 1B). The decline in sensitivity with higher temperatures is to a large extent explained by loss of low-lying glacier surface area and to a lesser extent by increasing precipitation adding mass to high-elevation glaciers. Although glaciers and thermal expansion have contributed approximately equally to the sea-level increase of the last 40 y (7), the sea-level commitment from glaciers is relatively small compared with thermal expansion.
The Greenland Ice Sheet.
Although it remains a challenge to simulate rapid ice discharge from the Greenland Ice Sheet in response to oceanic forcing (23), these fast ice fluxes are not crucial for a multimillennial estimate as attempted here. On a time scale of tens of thousands of years, the Greenland Ice Sheet shows threshold behavior with respect to the surrounding atmospheric temperature (24⇓⇓–27). Because summer temperatures around the ice sheet’s margins are warm enough to produce melt over a large area of the ice sheet, perturbations in the climate strongly affect the surface mass balance of the ice sheet. Although the associated changes are most likely not abrupt in a temporal sense, they self-amplify owing to positive feedbacks, particularly that between surface elevation and temperature (28).
For the multimillennial contribution of the Greenland Ice Sheet, we apply the recent results from an ensemble of simulations from a regional energy–moisture balance climate model coupled to an ice-sheet model that accounts for the positive feedback between temperature and surface elevation (25). The model’s parameters were constrained by comparison with surface mass balance estimates and topographical data for the present day and with estimated summit-elevation changes from ice-core records for the Last Interglacial period (LIG) (29), to ensure that the coupled model ensemble has a realistic sensitivity to climatic changes. In the transient response to global warming, simulated ice-sheet melting is comparable in timing and distribution to that of a GCM coupled to an ice sheet (24, 30).
The contribution to sea-level commitment from the Greenland Ice Sheet is relatively weak (on average 0.18 m °C−1 up to 1 °C and 0.34 m °C−1 between 2 and 4 °C) apart from the abrupt threshold of ice loss between 0.8 and 2.2 °C above preindustrial (90% credible interval) (Fig. 1C). This corresponds to a transition from a fully ice-covered Greenland to an essentially ice-free state (i.e., a reduction in ice volume of approximately 10% of the present-day volume, corresponding to a sea-level contribution of more than 6 m). Compared with previous studies (24, 30) the model applied here shows a threshold at lower temperatures than would correspond to a negative surface mass balance of the entire ice sheet because of dynamic ice motion. This explains the lower threshold temperatures compared with earlier studies. However, the uncertainty range plotted here shows a significant overlap with these earlier estimates.
The Antarctic Ice Sheet.
The total volume of the Antarctic Ice Sheet is equivalent to 55–60 m of global mean sea level (31). Most of the ice loss on Antarctica is due to solid-ice discharge, and thus the representation of the ice flow and sliding is particularly important. For continental-scale models, as required for estimates of the total sea-level contribution, a realistic representation of the grounding line position is a challenge at horizontal resolutions that allow for a simulation of whole Antarctica (32⇓–34). The model used here applies a parameterization of the ice flux at the grounding line (35) that captures grounding line motion realistically on multimillennial time scales as suggested by comparison of the model’s simulation of the last 5 million years (36) with local sediment records (37). Here we extract the sensitivity of the ice sheet from this model simulation by correlating the ice volume with the global mean temperature that forces the simulation. To this end we average the yearly temperature- and sea-level contribution data from the 5-million-year simulation over periods of 1,000 y (small dots in Fig. 1D). The sea-level contribution is then binned into temperature intervals of 0.2 °C. The median of the distribution within each temperature bin is provided as an estimate of the temperature sensitivity of the Antarctic ice sheet (thick red line in Fig. 1D). The uncertainty range is estimate as the 66% percentile of the distribution around the median (red shading). This uncertainty arises from uncertainty in the forcing data, the ice physics representation but also from the time-dependent nature of the simulation. For example, the existence of hysteresis behavior on the subcontinental scale can lead to different contributions for the same temperature increase.
The Antarctic Ice Sheet shows a relatively constant commitment of 1.2 m °C−1 (dashed line in Fig. 1D). This mostly involves retreat of marine portions (grounded below sea level) of the ice sheet due to increased sub-ice shelf oceanic melt, predominantly in West Antarctica; it does not include future expansion of East Antarctic volume due to anthropogenic increases in snowfall, which is expected to be a minor effect on millennial scales (see below). Although threshold behavior can occur during marine retreat, the response shown is quite linear owing to temporal smoothing and the rapidly varying paleo-climatic forcing in the transient simulation. Future anthropogenic changes will occur on even faster time scales, so this approach is arguably more relevant than fully equilibrated sensitivity tests.
Similar to Greenland, the much larger East Antarctic Ice Sheet, which is predominantly grounded above sea level, has the potential of threshold behavior in response to increased surface melting (38). Current models (4, 5) suggest that the required warming is much larger than anticipated in the near future (equivalent to ∼4 times preindustrial level of CO2 or higher). Here we only discuss the commitment up to 4 °C of warming, neglecting the possible deglaciation of East Antarctica. However, there is an apparent conflict with geologic evidence of up to 20 m sea-level rise during the Miocene (39), which implies substantial contributions from East Antarctica at times when estimated CO2 levels seem to have remained far below the model-derived thresholds. More work is needed to confidently rule out the possibility of East Antarctic threshold loss in the next few millennia.
Paleo-Evidence
To compare the model results with past sea-level anomalies for the temperature range up to 4 °C, we focus on three previous periods for which the geological record provides reasonable constraints on warmer climates and higher sea levels than preindustrial: the middle Pliocene, marine isotope stage 11, and the LIG (Fig. 1E).
During the obliquity-paced warm intervals of the middle Pliocene (∼3.3 to ∼3 Ma), reconstructions of sea-surface temperatures (40) and climate model simulations (11, 41) suggest that peak global surface air temperatures were 1.8–3.5 °C warmer than preindustrial. Estimates of peak mid-Pliocene sea levels based on a variety of geological records are consistent in suggesting higher-than-present sea levels, but they range widely (5–40 m) and are each subject to large uncertainties. For example, coastal records (shorelines, continental margin sequences) are influenced by glacio-hydro-isostatic (42) or global mantle dynamic processes (43). Both signals are large (5–30 m) and uncertain, and significant differences in published predictions of dynamic topography suggest that the latter is particularly poorly constrained.
Benthic δ18O records are better dated than many coastal records and provide a continuous time series, but the δ18O signal reflects some unknown combination of ice volume and temperature, with some contribution from regional hydrographic variability also possible. During the mid-Pliocene warm intervals, the LR04 benthic δ18O stack (44) identifies low values that exceeded present values by 0.1–0.25 per mil. If the signal only records ice volume, then on the basis of the relationship of the δ18O of seawater to ice volume derived from pore water chemistry (∼0.08 per mil per 10 m of sea-level equivalent) (45), and assuming this relationship still held during the Pliocene, the lowest δ18O values in the LR04 stack suggest peak sea levels ∼12–31 m higher than present. Dowsett et al. (40), however, reconstructed Pliocene deepwater (>2,000 m) temperature anomalies at 20 sites that ranged from −0.9 °C to 4.2 °C, with an average of 1 ± 1.2 °C (1 SD) warmer than present. A δ18O-temperature sensitivity of 0.28 per mil per °C suggests that, given these values, the lowest mid-Pliocene benthic δ18O values may be entirely explained by warmer deepwater temperature. Attempts to constrain the temperature component in benthic δ18O records indicate higher-than-present sea level during mid-Pliocene warm periods (46, 47), but these have large uncertainties (±15–25 m) (48).
Records of ice-rafted debris off Prydz Bay, East Antarctica (49) and the sedimentary record from the Ross Embayment of the West Antarctic Ice Sheet (37) provide the most direct evidence that mid-Pliocene sea levels must have been higher than present. This is consistent with mid-Pliocene Ross Sea surface temperatures reconstructed from biological and geochemical proxy methods, which have a range from 2 to 8 °C warmer than present (50), with values of >5 °C being, according to one ice-sheet model, above the stability threshold for ice shelves and marine portions of the West Antarctic Ice Sheet (36). Finally, ice-sheet models forced by warm-interval mid-Pliocene temperature reconstructions and orbital forcing suggest near-complete deglaciation of Greenland (7 m sea-level equivalent) (51, 52).
We conclude that the balance of evidence indicates higher sea levels during the mid-Pliocene warm intervals, with the most robust lines of evidence coming from sedimentary records that suggest episodic deglaciation of the West Antarctic Ice Sheet and from ice-sheet models that suggest near-complete deglaciation of the Greenland Ice Sheet under warm Pliocene boundary conditions. We thus conclude that sea level was at least 7 m above present during mid-Pliocene warm periods, while, allowing for uncertainties in less direct sea-level proxies, did not exceed 20 m above present (39).
During marine isotope stage 11 (MIS 11; ∼411–401 ka), Antarctic ice-core and tropical Pacific paleo-temperature estimates suggest that the global temperature was 1.5–2.0 °C warmer than preindustrial (53). Published studies of the magnitude of emergent shorelines attributed to MIS 11 have generated highly divergent estimates, from ∼22 m above present (54) to levels near present (55). After correcting sites for glacial isostatic adjustment (GIA), Raymo and Mitrovica (56) estimated that peak sea level during MIS 11 was 6 m to 13 m higher than present. Muhs et al. (57) derived similar estimates for Curaçao of 8.4 m to 10 m higher than present, a location for which the GIA effects are small (56, 58). We thus assign a possible sea-level range of 6–13 m higher than present for the MIS 11. A substantial fraction of this sea-level rise can be directly attributed to a substantial reduction in the Greenland Ice Sheet, as suggested by pollen records, which indicate development of a spruce forest over Greenland at this time (59).
During the LIG (∼130 ka to 116 ka), greenhouse gas concentrations were similar to preindustrial (280 ppm), but it experienced stronger orbital forcing than the present interglaciation, resulting in large positive insolation anomalies during boreal summer in the Northern Hemisphere and austral spring in the Southern Hemisphere. Temperature estimates suggest that global mean surface temperature was ∼1–2 °C warmer than preindustrial (60), with sea surface temperature warming being 0.7 ± 0.6 °C greater than present (50).
Kopp et al. (10) made probabilistic estimates of the global mean sea level from a large database of LIG sea-level estimates and concluded that global mean sea level was at least 6.4 m higher than present (95% probability) and that it is unlikely to have exceeded 8.8 m (33% probability). Dutton and Lambeck (58) restricted their analysis to well-dated sites but similarly concluded that global mean sea level was 5.5 m to 9 m higher than present. We thus assign a possible sea-level range of 5–9 m higher than present for the LIG.
Sea-Level Commitment for the Next 2,000 Years
The ability of the physical models to reproduce paleo-climatic records on a multimillennial time scale provides confidence in applying them to a time frame that is relevant for the societal discussion on future warming scenarios, as well as long enough for the model to be a valid representation of the dynamics involved. On a 2,000-y time scale, the sea-level contribution will be largely independent of the exact warming path during the first century. At the same time, 2,000 y is a relevant time scale, for example, for society’s cultural heritage. As can be seen from figure 10.34 in ref. 1, the oceanic heat content will be largely equilibrated after 2,000 y (Fig. 2A); the same is true for the glacier component (Fig. 2B). The situation for Antarctica is slightly more complicated, but as can be inferred from ref. 36 much of the West Antarctic retreat will have already occurred by 2,000 y, especially if the warming occurs on a decadal to centennial time scale. The opposite and smaller trend in East Antarctic ice volume due to increased snowfall in a warmer environment will also have largely equilibrated (61). The most significant difference arises from the contribution of Greenland. Consistent with previous estimates (3⇓⇓–6, 62), the rate of the sea-level contribution from Greenland increases with temperature. The transient simulations for an instantaneous temperature increase show a superlinear dependence of the sea-level contribution after 2,000 y (Fig. 2C). The results are quantitatively consistent with previous estimates on a millennial time scale (3).
Sea-level commitment per degree of warming as in Fig. 1 except for sea-level contributions after 2,000 y. Results were obtained from physical model simulations of (A) ocean warming, (B) mountain glaciers and ice caps, and (C) the Greenland and (D) the Antarctic Ice Sheets. (E) The corresponding total sea-level commitment. Temperatures are relative to preindustrial. Dashed lines and large black dots provide linear approximations: (A) sea-level rise for a homogeneous increase in temperature; (D and E) constant slopes of 1.2 and 2.3 m/°C. Shading represents the uncertainty range as discussed in the text. The large green dot in C corresponds to the sea-level contribution of Greenland (following ref. 3) after 1,000 y of integration with 560 ppm atmospheric CO2 concentration and a standard polar amplification of 1.5, which translates Greenland temperatures to the global mean temperature.
Mass loss from land ice results in a spatially variable sea-level change due to the resulting isostatic deformation and changes in gravity (63, 64). Using a model that simulates these processes (65⇓–67), we computed global patterns of sea-level change associated with mass loss from the Greenland and Antarctic Ice Sheets based on the volume contributions in Table 1 (Fig. 3). Information about the spatial distribution of mass loss from within each ice sheet is not known, and so the model runs assumed a uniform thinning of ice over a 2,000-y period with no margin retreat. Although this simplification will affect the accuracy of the results (68, 69), the general characteristics of the pattern will not be affected. Over millennial time scales, the dominant component of solid Earth deformation is nonelastic, and so we used a viscoelastic Maxwell rheology (70). The patterns will be sensitive to the adopted Earth viscosity profile (Fig. 3 legend) and the time span of the simulation (68).
Sea-level commitment from the different cryospheric components and ocean warming for different levels of global warming above preindustrial
Regional patterns of sea-level change computed using an isostatic surface loading model (see text for details) for scenarios of 1, 2, 3, and 4 °C of warming. These results are based on the assumption that ice thinning was uniform over the ice sheets and progressed linearly for 2,000 y. The contribution of glaciers and ice caps and ocean thermal expansion to the spatial pattern is not included; however, their global mean sea-level contribution is included (Table 1). These results are based on a spherically symmetric Earth model with 1D, depth-dependent viscosity structure. This structure comprises an elastic outer layer of thickness 96 km (to simulate the lithosphere) and two sublithosphere layers with viscosity values 5 × 1020 Pas (96 km to 660 km depth) and 1022 Pas (660 km to 2,900 km). These parameters are compatible with those inferred in a number of studies (e.g. refs 72 and 73). The white contour in each frame denotes the global mean sea-level value (Table 1).
We note that the spatial contributions from oceanic warming and glaciers and ice caps are not considered in Fig. 3 (only the global mean of these contributions is shown). Compared with the global mean values (white contours in Fig. 3), the key characteristic of the patterns is a fall near the ice margins, a reduced sea-level rise (relative to the global mean) at intermediate to high latitudes, and an enhanced rise at intermediate to low latitudes. Note that the 1 °C scenario includes only a small amount of melting from the Greenland ice sheet (Table 1) and so sea-level rise in the Northern Hemisphere is relatively uniform.
Discussion and Conclusion
Although the ice-sheet and climate models used here are not necessarily applicable to the fast time scales required to project sea-level rise of this century, their ability to simulate the long-term signal has been validated against paleo-evidence on longer time scales and is consistent with results from other models on these time scales. In particular, model simulations providing the separate contributions from ocean warming, glaciers, and the Greenland and Antarctic Ice Sheets are consistent with paleo-estimates of total sea-level rise in warmer periods. Together with their underlying physical theories, this provides confidence in our understanding of long-term sea-level commitments. It should be noted that the transient response within the next century, especially of the Antarctic Ice Sheet, might be different even qualitatively, owing to increased snowfall in a warming environment (61), which may induce dynamic responses of the ice sheet (71). Consistent with paleo-records, however, Antarctica will contribute to global sea-level on a multimillennial time scale.
Although the Antarctic Ice Sheet and the thermal expansion of ocean water contribute quasi-linearly, with approximately 1.2 m °C−1 and 0.4 m °C−1, respectively, Greenland shows a threshold behavior with a stepwise increase of approximately 6 m on a time scale on the order of several ten thousand years. After 2,000 y the model shows a superlinear response to the temperature increase. The contribution of glaciers declines at higher temperatures owing to the limited and decreasing ice stored there. As a consequence, the total sea-level commitment after 2,000 y is quasi-linear, with a sensitivity of 2.3 m °C−1.
Acknowledgments
This study was supported by the German Federal Ministry of Education and Research. P.U.C. and D.P. receive support from National Science Foundation Grants 1043517 and 1043018, respectively. B.M. was supported by Austrian Science Fund (FWF): P25362-N26. G.A.M. acknowledges support from the Canada Research Chairs program. A.R. was supported by the Gobierno de España, Ministerio de Economía y Competitividad under project CGL2011-29672-C02-01.
Footnotes
- ↵1To whom correspondence should be addressed. E-mail: Anders.levermann{at}pik-potsdam.de.
Author contributions: A.L. designed research; A.L., P.U.C., B.M., G.A.M., D.P., V.R., and A.R. performed research; B.M., G.A.M., D.P., V.R., and A.R. contributed new reagents/analytic tools; and A.L. and P.U.C. wrote the paper.
The authors declare no conflict of interest.
This article is a PNAS Direct Submission. J.C.M. is a guest editor invited by the Editorial Board.
See Commentary on page 13699.
Freely available online through the PNAS open access option.
References
- ↵
- Solomon S,
- et al.
- ↵
- Rahmstorf S
- ↵
- ↵
- Goelzer H,
- et al.
- ↵
- ↵
- Solomon S,
- Plattner GK,
- Knutti R,
- Friedlingstein P
- ↵
- ↵
- ↵
- Pfeffer WT,
- Harper JT,
- O’Neel S
- ↵
- ↵
- ↵
- Williams R,
- Goodwin P,
- Ridgwell A,
- Woodworth P
- ↵
- ↵
- ↵
- Pardaens AK,
- Lowe JA,
- Brown S,
- Nicholls RJ,
- de Gusmão D
- ↵
- Yin J
- ↵
- ↵
- Meier MF,
- et al.
- ↵
- ↵
- ↵
- ↵
- ↵
- Price SF,
- Payne AJ,
- Howat IM,
- Smith BE
- ↵
- Ridley J,
- Gregory JM,
- Huybrechts P,
- Lowe J
- ↵
- ↵
- ↵
- Letreguilly A,
- Huybrechts P,
- Reeh N
- ↵
- ↵
- Robinson A,
- Calov R,
- Ganopolski A
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- Foster GL,
- Rohling EJ
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- Schrag D
- ↵
- ↵
- Sosdian S,
- Rosenthal Y
- ↵
- Miller KG,
- et al.
- ↵
- ↵
- McKay NP,
- Overpeck JT,
- Otto-Bliesner BL
- ↵
- ↵
- Hill D,
- Dolan A,
- Haywood A,
- Hunter S,
- Stoll D
- ↵
- ↵
- Hearty P,
- Kindler P,
- Cheng H,
- Edwards R
- ↵
- Bowen D
- ↵
- ↵
- ↵
- Dutton A,
- Lambeck K
- ↵
- de Vernal A,
- Hillaire-Marcel C
- ↵
- ↵
- Uotila P,
- Lynch AH,
- Cassano JJ,
- Cullather RI
- ↵
- ↵
- ↵
- ↵
- Mitrovica JX,
- Milne GA
- ↵
- Kendall RA,
- Mitrovica JX,
- Milne GA
- ↵
- Mitrovica JX,
- Wahr J,
- Matsuyama I,
- Paulson A
- ↵
- Gomez N,
- Mitrovica JX,
- Tamisiea ME,
- Clark PU
- ↵
- Mitrovica JX,
- et al.
- ↵
- ↵
- ↵
- ↵
Citation Manager Formats
Sign up for Article Alerts
Article Classifications
- Physical Sciences
- Environmental Sciences


















See related content: