Sea level fingerprinting of the Bering Strait flooding history detects the source of the Younger Dryas climate event.

During the Last Glacial Maximum, expansive continental ice sheets lowered globally averaged sea level ~130 m, exposing a land bridge at the Bering Strait. During the subsequent deglaciation, sea level rose rapidly and ultimately flooded the Bering Strait, linking the Arctic and Pacific Oceans. Observational records of the Bering Strait flooding have suggested two apparently contradictory scenarios for the timing of the reconnection. We reconcile these enigmatic datasets using gravitationally self-consistent sea-level simulations that vary the timing and geometry of ice retreat between the Laurentide and Cordilleran Ice Sheets to the southwest of the Bering Strait to fit observations of a two-phased flooding history. Assuming the datasets are robust, we demonstrate that their reconciliation requires a substantial melting of the Cordilleran and western Laurentide Ice Sheet from 13,000 to 11,500 years ago. This timing provides a freshwater source for the widely debated Younger Dryas cold episode (12,900 to 11,700 years ago).

Sea-level records from the U.S. mid-Atlantic constrain Laurentide Ice Sheet extent during Marine Isotope Stage 3.

The U.S. mid-Atlantic sea-level record is sensitive to the history of the Laurentide Ice Sheet as the coastline lies along the ice sheet's peripheral bulge. However, paleo sea-level markers on the present-day shoreline of Virginia and North Carolina dated to Marine Isotope Stage (MIS) 3, from 50 to 35 ka, are surprisingly high for this glacial interval, and remain unexplained by previous models of ice age adjustment or other local (for example, tectonic) effects. Here, we reconcile this sea-level record using a revised model of glacial isostatic adjustment characterized by a peak global mean sea level during MIS 3 of approximately -40 m, and far less ice volume within the eastern sector of the Laurentide Ice Sheet than traditional reconstructions for this interval. We conclude that the Laurentide Ice Sheet experienced a phase of very rapid growth in the 15 kyr leading into the Last Glacial Maximum, thus highlighting the potential of mid-field sea-level records to constrain areal extent of ice cover during glacial intervals with sparse geological observables.

Mechanisms for oscillatory true polar wander.

Palaeomagnetic studies of Palaeoproterozoic to Cretaceous rocks propose a suite of large and relatively rapid (tens of degrees over 10 to 100 million years) excursions of the rotation pole relative to the surface geography, or true polar wander (TPW). These excursions may be linked in an oscillatory, approximately coaxial succession about the centre of the contemporaneous supercontinent. Within the framework of a standard rotational theory, in which a delayed viscous adjustment of the rotational bulge acts to stabilize the rotation axis, geodynamic models for oscillatory TPW generally appeal to consecutive, opposite loading phases of comparable magnitude. Here we extend a nonlinear rotational stability theory to incorporate the stabilizing effect of TPW-induced elastic stresses in the lithosphere. We demonstrate that convectively driven inertia perturbations acting on a nearly prolate, non-hydrostatic Earth with an effective elastic lithospheric thickness of about 10 kilometres yield oscillatory TPW paths consistent with palaeomagnetic inferences. This estimate of elastic thickness can be reduced, even to zero, if the rotation axis is stabilized by long-term excess ellipticity in the plane of the TPW. We speculate that these sources of stabilization, acting on TPW driven by a time-varying mantle flow field, provide a mechanism for linking the distinct, oscillatory TPW events of the past few billion years.

GRACE gravity data constrain ancient ice geometries and continental dynamics over Laurentia.

The free-air gravity trend over Canada, derived from the Gravity Recovery and Climate Experiment (GRACE) satellite mission, robustly isolates the gravity signal associated with glacial isostatic adjustment (GIA) from the longer-time scale mantle convection process. This trend proves that the ancient Laurentian ice complex was composed of two large domes to the west and east of Hudson Bay, in accord with one of two classes of earlier reconstructions. Moreover, GIA models that reconcile the peak rates contribute approximately 25 to approximately 45% to the observed static gravity field, which represents an important boundary condition on the buoyancy of the continental tectosphere.

Sea-level fingerprinting as a direct test for the source of global meltwater pulse IA.

The ice reservoir that served as the source for the meltwater pulse IA remains enigmatic and controversial. We show that each of the melting scenarios that have been proposed for the event produces a distinct variation, or fingerprint, in the global distribution of meltwater. We compare sea-level fingerprints associated with various melting scenarios to existing sea-level records from Barbados and the Sunda Shelf and conclude that the southern Laurentide Ice Sheet could not have been the sole source of the meltwater pulse, whereas a substantial contribution from the Antarctic Ice Sheet is consistent with these records.

Deep-mantle high-viscosity flow and thermochemical structure inferred from seismic and geodynamic data.

Surface geophysical data that are related to the process of thermal convection in the Earth's mantle provide constraints on the rheological properties and density structure of the mantle. We show that these convection-related data imply the existence of a region of very high effective viscosity near 2,000 km depth. This inference is obtained using a viscous-flow model based on recent high-resolution seismic models of three-dimensional structure in the mantle. The high-viscosity layer near 2,000 km depth results in a re-organization of flow from short to long horizontal length scales, which agrees with seismic tomographic observations of very long wavelength structures in the deep mantle. The high-viscosity region also strongly suppresses flow-induced deformation and convective mixing in the deep mantle. Here we predict compositional and thermal heterogeneity in this region, using viscous-flow calculations based on the new viscosity profile, together with independent mineral physics data. These maps are consistent with the anti-correlation of anomalies in seismic shear and bulk sound velocity in the deep mantle. The maps also show that mega-plumes in the lower mantle below the central Pacific and Africa are, despite the presence of compositional heterogeneity, buoyant and actively upwelling structures.

Space-geodetic constraints on glacial isostatic adjustment in Fennoscandia.

Analysis of Global Positioning System (GPS) data demonstrates that ongoing three-dimensional crustal deformation in Fennoscandia is dominated by glacial isostatic adjustment. Our comparison of these GPS observations with numerical predictions yields an Earth model that satisfies independent geologic constraints and bounds both the average viscosity in the upper mantle (5 x 10(20) to 1 x 10(21) pascal seconds) and the elastic thickness of the lithosphere (90 to 170 kilometers). We combined GPS-derived radial motions with Fennoscandian tide gauge records to estimate a regional sea surface rise of 2.1 +/- 0.3 mm/year. Furthermore, ongoing horizontal tectonic motions greater than approximately 1 mm/year are ruled out on the basis of the GPS-derived three-dimensional crustal velocity field.

Recent mass balance of polar ice sheets inferred from patterns of global sea-level change.

Global sea level is an indicator of climate change, as it is sensitive to both thermal expansion of the oceans and a reduction of land-based glaciers. Global sea-level rise has been estimated by correcting observations from tide gauges for glacial isostatic adjustment--the continuing sea-level response due to melting of Late Pleistocene ice--and by computing the global mean of these residual trends. In such analyses, spatial patterns of sea-level rise are assumed to be signals that will average out over geographically distributed tide-gauge data. But a long history of modelling studies has demonstrated that non-uniform--that is, non-eustatic--sea-level redistributions can be produced by variations in the volume of the polar ice sheets. Here we present numerical predictions of gravitationally consistent patterns of sea-level change following variations in either the Antarctic or Greenland ice sheets or the melting of a suite of small mountain glaciers. These predictions are characterized by geometrically distinct patterns that reconcile spatial variations in previously published sea-level records. Under the--albeit coarse--assumption of a globally uniform thermal expansion of the oceans, our approach suggests melting of the Greenland ice complex over the last century equivalent to -0.6 mm yr(-1) of sea-level rise.