Carling C. Hay

Temperature-driven global sea-level variability in the Common Era

Significance We present the first, to our knowledge, estimate of global sea-level (GSL) change over the last ∼3,000 years that is based upon statistical synthesis of a global database of regional sea-level reconstructions. GSL varied by ∼±8 cm over the pre-Industrial Common Era, with a notable decline over 1000–1400 CE coinciding with ∼0.2 °C of global cooling. The 20th century rise was extremely likely faster than during any of the 27 previous centuries. Semiempirical modeling indicates that, without global warming, GSL in the 20th century very likely would have risen by between −3 cm and +7 cm, rather than the ∼14 cm observed. Semiempirical 21st century projections largely reconcile differences between Intergovernmental Panel on Climate Change projections and semiempirical models. We assess the relationship between temperature and global sea-level (GSL) variability over the Common Era through a statistical metaanalysis of proxy relative sea-level reconstructions and tide-gauge data. GSL rose at 0.1 ± 0.1 mm/y (2σ) over 0–700 CE. A GSL fall of 0.2 ± 0.2 mm/y over 1000–1400 CE is associated with ∼0.2 °C global mean cooling. A significant GSL acceleration began in the 19th century and yielded a 20th century rise that is extremely likely (probability P≥0.95) faster than during any of the previous 27 centuries. A semiempirical model calibrated against the GSL reconstruction indicates that, in the absence of anthropogenic climate change, it is extremely likely (P=0.95) that 20th century GSL would have risen by less than 51% of the observed 13.8±1.5 cm. The new semiempirical model largely reconciles previous differences between semiempirical 21st century GSL projections and the process model-based projections summarized in the Intergovernmental Panel on Climate Change’s Fifth Assessment Report.

Geographic Variability of Sea-Level Change

Local sea-level changes differ significantly from global-mean sea-level change as a result of (1) non-climatic, geological background processes; (2) atmosphere/ocean dynamics; and (3) the gravitational, elastic, and rotational “fingerprint” effects of ice and ocean mass redistribution. Though the research communities working on these different effects each have a long history, the integration of all these different processes into interpretations of past changes and projections of future change is an active area of research. Fully characterizing the past contributions of these processes requires information from sources covering a range of timescales, including geological proxies, tide-gauge observations from the last ~3 centuries, and satellite-altimetry data from the last ~2 decades. Local sea-level rise projections must account for the different spatial patterns of different processes, as well as potential correlations between different drivers.

The sea-level fingerprints of ice-sheet collapse during interglacial periods

Studies of sea level during previous interglacials provide insight into the stability of polar ice sheets in the face of global climate change. Commonly, these studies correct ancient sea-level highstands for the contaminating effect of isostatic adjustment associated with past ice age cycles, and interpret the residuals as being equivalent to the peak eustatic sea level associated with excess melting, relative to present day, of ancient polar ice sheets. However, the collapse of polar ice sheets produces a distinct geometry, or fingerprint, of sea-level change, which must be accounted for to accurately infer peak eustatic sea level from site-specific residual highstands. To explore this issue, we compute fingerprints associated with the collapse of the Greenland Ice Sheet, West Antarctic Ice Sheet, and marine sectors of the East Antarctic Ice Sheet in order to isolate regions that would have been subject to greater-thaneustatic sea-level change for all three cases. These fingerprints are more robust than those associated with modern melting events, when applied to infer eustatic sea level, because: (1) a significant collapse of polar ice sheets reduces the sensitivity of the computed fingerprints to uncertainties in the geometry of the melt regions; and (2) the sea-level signal associated with the collapse will dominate the signal from steric effects. We evaluate these fingerprints at a suite of sites where sea-level records from interglacial marine isotopes stages (MIS) 5e and 11 have been obtained. Using these results, we demonstrate that previously discrepant estimates of peak eustatic sea level during MIS5e based on sealevel markers in Australia and the Seychelles are brought into closer accord.

Estimating the sources of global sea level rise with data assimilation techniques

A rapidly melting ice sheet produces a distinctive geometry, or fingerprint, of sea level (SL) change. Thus, a network of SL observations may, in principle, be used to infer sources of meltwater flux. We outline a formalism, based on a modified Kalman smoother, for using tide gauge observations to estimate the individual sources of global SL change. We also report on a series of detection experiments based on synthetic SL data that explore the feasibility of extracting source information from SL records. The Kalman smoother technique iteratively calculates the maximum-likelihood estimate of Greenland ice sheet (GIS) and West Antarctic ice sheet (WAIS) melt at each time step, and it accommodates data gaps while also permitting the estimation of nonlinear trends. Our synthetic tests indicate that when all tide gauge records are used in the analysis, it should be possible to estimate GIS and WAIS melt rates greater than ∼0.3 and ∼0.4 mm of equivalent eustatic sea level rise per year, respectively. We have also implemented a multimodel Kalman filter that allows us to account rigorously for additional contributions to SL changes and their associated uncertainty. The multimodel filter uses 72 glacial isostatic adjustment models and 3 ocean dynamic models to estimate the most likely models for these processes given the synthetic observations. We conclude that our modified Kalman smoother procedure provides a powerful method for inferring melt rates in a warming world.

Reconciling past changes in Earth’s rotation with 20th century global sea-level rise: Resolving Munk’s enigma

Reanalysis of Earth rotation observations reconciles them with globally averaged sea-level change in the 20th century. In 2002, Munk defined an important enigma of 20th century global mean sea-level (GMSL) rise that has yet to be resolved. First, he listed three canonical observations related to Earth’s rotation [(i) the slowing of Earth’s rotation rate over the last three millennia inferred from ancient eclipse observations, and changes in the (ii) amplitude and (iii) orientation of Earth’s rotation vector over the last century estimated from geodetic and astronomic measurements] and argued that they could all be fit by a model of ongoing glacial isostatic adjustment (GIA) associated with the last ice age. Second, he demonstrated that prevailing estimates of the 20th century GMSL rise (~1.5 to 2.0 mm/year), after correction for the maximum signal from ocean thermal expansion, implied mass flux from ice sheets and glaciers at a level that would grossly misfit the residual GIA-corrected observations of Earth’s rotation. We demonstrate that the combination of lower estimates of the 20th century GMSL rise (up to 1990) improved modeling of the GIA process and that the correction of the eclipse record for a signal due to angular momentum exchange between the fluid outer core and the mantle reconciles all three Earth rotation observations. This resolution adds confidence to recent estimates of individual contributions to 20th century sea-level change and to projections of GMSL rise to the end of the 21st century based on them.

Evolving Understanding of Antarctic Ice-Sheet Physics and Ambiguity in Probabilistic Sea-Level Projections

Mechanisms such as ice-shelf hydrofracturing and ice-cliff collapse may rapidly increase discharge from marine-based ice sheets. Here, we link a probabilistic framework for sea-level projections to a small ensemble of Antarctic ice-sheet (AIS) simulations incorporating these physical processes to explore their influence on global-mean sea-level (GMSL) and relative sea-level (RSL). We compare the new projections to past results using expert assessment and structured expert elicitation about AIS changes. Under high greenhouse gas emissions (Representative Concentration Pathway [RCP] 8.5), median projected 21st century GMSL rise increases from 79 to 146 cm. Without protective measures, revised median RSL projections would by 2100 submerge land currently home to 153 million people, an increase of 44 million. The use of a physical model, rather than simple parameterizations assuming constant acceleration of ice loss, increases forcing sensitivity: overlap between the central 90% of simulations for 2100 for RCP 8.5 (93–243 cm) and RCP 2.6 (26–98 cm) is minimal. By 2300, the gap between median GMSL estimates for RCP 8.5 and RCP 2.6 reaches >10 m, with median RSL projections for RCP 8.5 jeopardizing land now occupied by 950 million people (versus 167 million for RCP 2.6). The minimal correlation between the contribution of AIS to GMSL by 2050 and that in 2100 and beyond implies current sea-level observations cannot exclude future extreme outcomes. The sensitivity of post-2050 projections to deeply uncertain physics highlights the need for robust decision and adaptive management frameworks.

Sea Level Fingerprints in a Region of Complex Earth Structure: The Case of WAIS

AbstractSea level fingerprints associated with rapid melting of the West Antarctic Ice Sheet (WAIS) have generally been computed under the assumption of a purely elastic response of the solid Earth. The authors investigate the impact of viscous effects on these fingerprints by computing gravitationally self-consistent sea level changes that adopt a 3D viscoelastic Earth model in the Antarctic region consistent with available geological and geophysical constraints. In West Antarctica, the model is characterized by a thin (~65 km) elastic lithosphere and sublithospheric viscosities that span three orders of magnitude, reaching values as low as approximately 4 × 1018 Pa s beneath WAIS. Calculations indicate that sea level predictions in the near field of WAIS will depart significantly from elastic fingerprints in as little as a few decades. For example, when viscous effects are included, the peak sea level fall predicted in the vicinity of WAIS during a melt event will increase by about 20% and about 50%, re...

On the Robustness of Bayesian Fingerprinting Estimates of Global Sea Level Change

AbstractGlobal mean sea level (GMSL) over the twentieth century has been estimated using techniques that include regional averaging of sparse tide gauge observations, combining satellite altimetry observations with tide gauge records in empirical orthogonal function (EOF) analyses, and most recently the Bayesian approaches of Kalman smoothing (KS) and Gaussian process regression (GPR). Estimated trends in GMSL over 1901–90 obtained using the Bayesian techniques are 1.1–1.2 mm yr−1, approximately 20% lower than previous estimates. It has been suggested that the adoption of a less restrictive subset of records biased the Bayesian-derived estimates. In this study, different subsets of records are used to demonstrate that GMSL estimates based on the Bayesian methodologies are robust to tide gauge selection. A method for determining the resolvability of individual sea level components estimated in a Bayesian framework is also presented and applied. It is found that the incomplete tide gauge observations result...

Minimal East Antarctic Ice Sheet retreat onto land during the past eight million years

The East Antarctic Ice Sheet (EAIS) is the largest potential contributor to sea-level rise. However, efforts to predict the future evolution of the EAIS are hindered by uncertainty in how it responded to past warm periods, for example, during the Pliocene epoch (5.3 to 2.6 million years ago), when atmospheric carbon dioxide concentrations were last higher than 400 parts per million. Geological evidence indicates that some marine-based portions of the EAIS and the West Antarctic Ice Sheet retreated during parts of the Pliocene1,2, but it remains unclear whether ice grounded above sea level also experienced retreat. This uncertainty persists because global sea-level estimates for the Pliocene have large uncertainties and cannot be used to rule out substantial terrestrial ice loss3, and also because direct geological evidence bearing on past ice retreat on land is lacking. Here we show that land-based sectors of the EAIS that drain into the Ross Sea have been stable throughout the past eight million years. We base this conclusion on the extremely low concentrations of cosmogenic 10Be and 26Al isotopes found in quartz sand extracted from a land-proximal marine sediment core. This sediment had been eroded from the continent, and its low levels of cosmogenic nuclides indicate that it experienced only minimal exposure to cosmic radiation, suggesting that the sediment source regions were covered in ice. These findings indicate that atmospheric warming during the past eight million years was insufficient to cause widespread or long-lasting meltback of the EAIS margin onto land. We suggest that variations in Antarctic ice volume in response to the range of global temperatures experienced over this period—up to 2–3 degrees Celsius above preindustrial temperatures4, corresponding to future scenarios involving carbon dioxide concentrations of between 400 and 500 parts per million—were instead driven mostly by the retreat of marine ice margins, in agreement with the latest models5,6.Analysis of cosmogenic isotopes from a marine sediment core shows that much of the land-based East Antarctic Ice Sheet has remained stable for the past eight million years, including during the warm Pliocene epoch.

Bias in Estimates of Global Mean Sea Level Change Inferred from Satellite Altimetry

abstractEstimates of regional and global average sea level change remain a focus of climate change research. One complication in obtaining coherent estimates is that geodetic datasets measure different aspects of the sea level field. Satellite altimetry constrains changes in the sea surface height (SSH; or absolute sea level), whereas tide gauge data provide a measure of changes in SSH relative to the crust (i.e., relative sea level). The latter is a direct measure of changes in ocean volume (and the combined impacts of ice sheet melt and steric effects), but the former is not since it does not account for crustal deformation. Nevertheless, the literature commonly conflates the two estimates by directly comparing them. We demonstrate that using satellite altimetry records to estimate global ocean volume changes can lead to biases that can exceed 15%. The level of bias will depend on the relative contributions to sea level changes from the Antarctic and Greenland Ice Sheets. The bias is also more sensitive...