Constraints on global mean sea level during Pliocene warmth.

Reconstructing the evolution of sea level during past warmer epochs such as the Pliocene provides insight into the response of sea level and ice sheets to prolonged warming1. Although estimates of the global mean sea level (GMSL) during this time do exist, they vary by several tens of metres2-4, hindering the assessment of past and future ice-sheet stability. Here we show that during the mid-Piacenzian Warm Period, which was on average two to three degrees Celsius warmer than the pre-industrial period5, the GMSL was about 16.2 metres higher than today owing to global ice-volume changes, and around 17.4 metres when thermal expansion of the oceans is included. During the even warmer Pliocene Climatic Optimum (about four degrees Celsius warmer than pre-industrial levels)6, our results show that the GMSL was 23.5 metres above the present level, with an additional 1.6 metres from thermal expansion. We provide six GMSL data points, ranging from 4.39 to 3.27 million years ago, that are based on phreatic overgrowths on speleothems from the western Mediterranean (Mallorca, Spain). This record is unique owing to its clear relationship to sea level, its reliable U-Pb ages and its long timespan, which allows us to quantify uncertainties on potential uplift. Our data indicate that ice sheets are very sensitive to warming and provide important calibration targets for future ice-sheet models7.

GIA imaging of 3-D mantle viscosity based on palaeo sea level observations - Part I: Sensitivity kernels for an Earth with laterally varying viscosity.

A key initial step in geophysical imaging is to devise an effective means of mapping the sensitivity of an observation to the model parameters, that is to compute its Fréchet derivatives or sensitivity kernel. In the absence of any simplifying assumptions and when faced with a large number of free parameters, the adjoint method can be an effective and efficient approach to calculating Fréchet derivatives and requires just two numerical simulations. In the Glacial Isostatic Adjustment problem, these consist of a forward simulation driven by changes in ice mass and an adjoint simulation driven by fictitious loads that are applied at the observation sites. The theoretical basis for this approach has seen considerable development over the last decade. Here, we present the final elements needed to image 3-D mantle viscosity using a dataset of palaeo sea-level observations. Developments include the calculation of viscosity Fréchet derivatives (i.e. sensitivity kernels) for relative sea-level observations, a modification to the numerical implementation of the forward and adjoint problem that permits application to 3-D viscosity structure, and a recalibration of initial sea level that ensures the forward simulation honours present-day topography. In the process of addressing these items, we build intuition concerning how absolute sea-level and relative sea-level observations sense Earth's viscosity structure and the physical processes involved. We discuss examples for potential observations located in the near field (Andenes, Norway), far field (Seychelles), and edge of the forebulge of the Laurentide ice sheet (Barbados). Examination of these kernels: (1) reveals why 1-D estimates of mantle viscosity from far-field relative sea-level observations can be biased; (2) hints at why an appropriate differential relative sea-level observation can provide a better constraint on local mantle viscosity and (3) demonstrates that sea-level observations have non-negligible 3-D sensitivity to deep mantle viscosity structure, which is counter to the intuition gained from 1-D radial viscosity Fréchet derivatives. Finally, we explore the influence of lateral variations in viscosity on relative sea-level observations in the Amundsen Sea Embayment and at Barbados. These predictions are based on a new global 3-D viscosity inference derived from the shear-wave speeds of GLAD-M25 and an inverse calibration scheme that ensures compatibility with certain fundamental geophysical observations. Use of the 3-D viscosity inference leads to: (1) generally greater complexity within the kernel; (2) an increase in sensitivity and presence of shorter length-scale features within lower viscosity regions; (3) a zeroing out of the sensitivity kernel within high-viscosity regions where elastic deformation dominates and (4) shifting of sensitivity at a given depth towards distal regions of weaker viscosity. The tools and intuition built here provide the necessary framework to explore inversions for 3-D mantle viscosity based on palaeo sea-level data.

Sea-level trends across The Bahamas constrain peak last interglacial ice melt.

During the last interglacial (LIG) period, global mean sea level (GMSL) was higher than at present, likely driven by greater high-latitude insolation. Past sea-level estimates require elevation measurements and age determination of marine sediments that formed at or near sea level, and those elevations must be corrected for glacial isostatic adjustment (GIA). However, this GIA correction is subject to uncertainties in the GIA model inputs, namely, Earth's rheology and past ice history, which reduces precision and accuracy in estimates of past GMSL. To better constrain the GIA process, we compare our data and existing LIG sea-level data across the Bahamian archipelago with a suite of 576 GIA model predictions. We calculated weights for each GIA model based on how well the model fits spatial trends in the regional sea-level data and then used the weighted GIA corrections to revise estimates of GMSL during the LIG. During the LIG, we find a 95% probability that global sea level peaked at least 1.2 m higher than today, and it is very unlikely (5% probability) to have exceeded 5.3 m. Estimates increase by up to 30% (decrease by up to 20%) for portions of melt that originate from the Greenland ice sheet (West Antarctic ice sheet). Altogether, this work suggests that LIG GMSL may be lower than previously assumed.

Glacial isostatic adjustment reduces past and future Arctic subsea permafrost.

Sea-level rise submerges terrestrial permafrost in the Arctic, turning it into subsea permafrost. Subsea permafrost underlies ~ 1.8 million km2 of Arctic continental shelf, with thicknesses in places exceeding 700 m. Sea-level variations over glacial-interglacial cycles control subsea permafrost distribution and thickness, yet no permafrost model has accounted for glacial isostatic adjustment (GIA), which deviates local sea level from the global mean due to changes in ice and ocean loading. Here we incorporate GIA into a pan-Arctic model of subsea permafrost over the last 400,000 years. Including GIA significantly reduces present-day subsea permafrost thickness, chiefly because of hydro-isostatic effects as well as deformation related to Northern Hemisphere ice sheets. Additionally, we extend the simulation 1000 years into the future for emissions scenarios outlined in the Intergovernmental Panel on Climate Change's sixth assessment report. We find that subsea permafrost is preserved under a low emissions scenario but mostly disappears under a high emissions scenario.

Geodynamically corrected Pliocene shoreline elevations in Australia consistent with midrange projections of Antarctic ice loss.

The Mid-Pliocene represents the most recent interval in Earth history with climatic conditions similar to those expected in the coming decades. Mid-Pliocene sea level estimates therefore provide important constraints on projections of future ice sheet behavior and sea level change but differ by tens of meters due to local distortion of paleoshorelines caused by mantle dynamics. We combine an Australian sea level marker compilation with geodynamic simulations and probabilistic inversions to quantify and remove these post-Pliocene vertical motions at continental scale. Dynamic topography accounts for most of the observed sea level marker deflection, and correcting for this effect and glacial isostatic adjustment yields a Mid-Pliocene global mean sea level of +16.0 (+10.4 to +21.5) m (50th/16th to 84th percentiles). Recalibration of recent high-end sea level projections using this revised estimate implies a more stable Antarctic Ice Sheet under future warming scenarios, consistent with midrange forecasts of sea level rise that do not incorporate a marine ice cliff instability.

Constraining the contribution of the Antarctic Ice Sheet to Last Interglacial sea level.

Polar temperatures during the Last Interglacial [LIG; ~129 to 116 thousand years (ka)] were warmer than today, making this time period an important testing ground to better understand how ice sheets respond to warming. However, it remains debated how much and when the Antarctic and Greenland ice sheets changed during this period. Here, we present a combination of new and existing absolutely dated LIG sea-level observations from Britain, France, and Denmark. Because of glacial isostatic adjustment (GIA), the LIG Greenland ice melt contribution to sea-level change in this region is small, which allows us to constrain Antarctic ice change. We find that the Antarctic contribution to LIG global mean sea level peaked early in the interglacial (before 126 ka), with a maximum contribution of 5.7 m (50th percentile, 3.6 to 8.7 m central 68% probability) before declining. Our results support an asynchronous melt history over the LIG, with an early Antarctic contribution followed by later Greenland Ice Sheet mass loss.

Influence of reef isostasy, dynamic topography, and glacial isostatic adjustment on sea-level records in Northeastern Australia.

Understanding sea level during the peak of the Last Interglacial (125,000 yrs ago) is important for assessing future ice-sheet dynamics in response to climate change. The coasts and continental shelves of northeastern Australia (Queensland) preserve an extensive Last Interglacial record in the facies of coastal strandplains onland and fossil reefs offshore. However, there is a discrepancy, amounting to tens of meters, in the elevation of sea-level indicators between offshore and onshore sites. Here, we assess the influence of geophysical processes that may have changed the elevation of these sea-level indicators. We modeled sea-level change due to dynamic topography, glacial isostatic adjustment, and isostatic adjustment due to coral reef loading. We find that these processes caused relative sea-level changes on the order of, respectively, 10 m, 5 m, and 0.3 m. Of these geophysical processes, the dynamic topography predictions most closely match the tilting observed between onshore and offshore sea-level markers.

Toward a Cenozoic history of atmospheric CO2.

The geological record encodes the relationship between climate and atmospheric carbon dioxide (CO2) over long and short timescales, as well as potential drivers of evolutionary transitions. However, reconstructing CO2 beyond direct measurements requires the use of paleoproxies and herein lies the challenge, as proxies differ in their assumptions, degree of understanding, and even reconstructed values. In this study, we critically evaluated, categorized, and integrated available proxies to create a high-fidelity and transparently constructed atmospheric CO2 record spanning the past 66 million years. This newly constructed record provides clearer evidence for higher Earth system sensitivity in the past and for the role of CO2 thresholds in biological and cryosphere evolution.

Total isostatic response to the complete unloading of the Greenland and Antarctic Ice Sheets.

The land surface beneath the Greenland and Antarctic Ice Sheets is isostatically suppressed by the mass of the overlying ice. Accurate computation of the land elevation in the absence of ice is important when considering, for example, regional geodynamics, geomorphology, and ice sheet behaviour. Here, we use contemporary compilations of ice thickness and lithospheric effective elastic thickness to calculate the fully re-equilibrated isostatic response of the solid Earth to the complete removal of the Greenland and Antarctic Ice Sheets. We use an elastic plate flexure model to compute the isostatic response to the unloading of the modern ice sheet loads, and a self-gravitating viscoelastic Earth model to make an adjustment for the remaining isostatic disequilibrium driven by ice mass loss since the Last Glacial Maximum. Feedbacks arising from water loading in areas situated below sea level after ice sheet removal are also taken into account. In addition, we quantify the uncertainties in the total isostatic response associated with a range of elastic and viscoelastic Earth properties. We find that the maximum change in bed elevation following full re-equilibration occurs over the centre of the landmasses and is +783 m in Greenland and +936 m in Antarctica. By contrast, areas around the ice margins experience up to 123 m of lowering due to a combination of sea level rise, peripheral bulge collapse, and water loading. The computed isostatic response fields are openly accessible and have a number of applications for studying regional geodynamics, landscape evolution, cryosphere dynamics, and relative sea level change.

The role of gravitational body forces in the development of metamorphic core complexes.

Within extreme continental extension areas, ductile middle crust is exhumed at the surface as metamorphic core complexes. Sophisticated quantitative models of extreme extension predicted upward transport of ductile middle-lower crust through time. Here we develop a general model for metamorphic core complexes formation and demonstrate that they result from the collapse of a mountain belt supported by a thickened crustal root. We show that gravitational body forces generated by topography and crustal root cause an upward flow pattern of the ductile lower-middle crust, facilitated by a detachment surface evolving into low-angle normal fault. This detachment surface acquires large amounts of finite strain, consistent with thick mylonite zones found in metamorphic core complexes. Isostatic rebound exposes the detachment in a domed upwarp, while the final Moho discontinuity across the extended region relaxes to a flat geometry. This work suggests that belts of metamorphic core complexes are a fossil signature of collapsed highlands.

Coupled influence of tectonics, climate, and surface processes on landscape evolution in southwestern North America.

The Cenozoic landscape evolution in southwestern North America is ascribed to crustal isostasy, dynamic topography, or lithosphere tectonics, but their relative contributions remain controversial. Here we reconstruct landscape history since the late Eocene by investigating the interplay between mantle convection, lithosphere dynamics, climate, and surface processes using fully coupled four-dimensional numerical models. Our quantified depth-dependent strain rate and stress history within the lithosphere, under the influence of gravitational collapse and sub-lithospheric mantle flow, show that high gravitational potential energy of a mountain chain relative to a lower Colorado Plateau can explain extension directions and stress magnitudes in the belt of metamorphic core complexes during topographic collapse. Profound lithospheric weakening through heating and partial melting, following slab rollback, promoted this extensional collapse. Landscape evolution guided northeast drainage onto the Colorado Plateau during the late Eocene-late Oligocene, south-southwest drainage reversal during the late Oligocene-middle Miocene, and southwest drainage following the late Miocene.

Sea-level stands from the Western Mediterranean over the past 6.5 million years.

Sea-level reconstructions are important for understanding past ice sheet variability and its response to past and future warming. Here we present Neogene and Quaternary sea-level snapshots using phreatic overgrowths on speleothems (POS) from caves on Mallorca, Spain. POS are excellent sea level index points because of their clear relationship to sea level and precise U-Pb chronology. We find that local sea-level before and at the onset of the Messinian Salinity Crisis was at 33.3 ± 0.25 m (6.54 ± 0.37 Ma) and 31.8 ± 0.25 m (5.86 ± 0.60 Ma) above present levels, respectively. We further present global mean sea level (GMSL) estimates, i.e. local sea level corrected for glacial isostatic adjustment and long-term uplift, for three other POS. The results show that GMSL during the Pliocene-Pleistocene Transition was 6.4 m (- 2.0-8.8 m) at 2.63 ± 0.11 Ma and during the beginning and the end of the Mid-Pleistocene Transition was - 1.1 m (- 5.6-2.4 m) and 5 m (1.5-8.1 m), respectively. These estimates provide important constraints for the past evolution of sea level and show that local sea level prior to the MSC was similar to the highest stand during the Pliocene, with markedly lower position afterwards.

Long-Term Earth-Moon Evolution With High-Level Orbit and Ocean Tide Models.

Tides and Earth-Moon system evolution are coupled over geological time. Tidal energy dissipation on Earth slows E a r t h ' s rotation rate, increases obliquity, lunar orbit semi-major axis and eccentricity, and decreases lunar inclination. Tidal and core-mantle boundary dissipation within the Moon decrease inclination, eccentricity and semi-major axis. Here we integrate the Earth-Moon system backwards for 4.5 Ga with orbital dynamics and explicit ocean tide models that are "high-level" (i.e., not idealized). To account for uncertain plate tectonic histories, we employ Monte Carlo simulations, with tidal energy dissipation rates (normalized relative to astronomical forcing parameters) randomly selected from ocean tide simulations with modern ocean basin geometry and with 55, 116, and 252 Ma reconstructed basin paleogeometries. The normalized dissipation rates depend upon basin geometry and E a r t h ' s rotation rate. Faster Earth rotation generally yields lower normalized dissipation rates. The Monte Carlo results provide a spread of possible early values for the Earth-Moon system parameters. Of consequence for ocean circulation and climate, absolute (un-normalized) ocean tidal energy dissipation rates on the early Earth may have exceeded t o d a y ' s rate due to a closer Moon. Prior to ∼ 3 Ga , evolution of inclination and eccentricity is dominated by tidal and core-mantle boundary dissipation within the Moon, which yield high lunar orbit inclinations in the early Earth-Moon system. A drawback for our results is that the semi-major axis does not collapse to near-zero values at 4.5 Ga, as indicated by most lunar formation models. Additional processes, missing from our current efforts, are discussed as topics for future investigation.

Plasmon mediated confocal dark-field microscopy.

An efficient mode for scanning confocal dark-field microscopy through a thin gold film is established that takes advantage of the intermediate excitation of surface plasmons both in the excitation and in the emission process. This concept is verified by experimental investigation of the effective point-spread function, the intensity distribution of the scattered radiation and by comparison with a classical dark-field geometry. The wavelength-dependence of both the signal strength and the point-spread function are discussed.

Detection of a dynamic topography signal in last interglacial sea-level records.

Estimating minimum ice volume during the last interglacial based on local sea-level indicators requires that these indicators are corrected for processes that alter local sea level relative to the global average. Although glacial isostatic adjustment is generally accounted for, global scale dynamic changes in topography driven by convective mantle flow are generally not considered. We use numerical models of mantle flow to quantify vertical deflections caused by dynamic topography and compare predictions at passive margins to a globally distributed set of last interglacial sea-level markers. The deflections predicted as a result of dynamic topography are significantly correlated with marker elevations (>95% probability) and are consistent with construction and preservation attributes across marker types. We conclude that a dynamic topography signal is present in the elevation of last interglacial sea-level records and that the signal must be accounted for in any effort to determine peak global mean sea level during the last interglacial to within an accuracy of several meters.