Search

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.

The first records of Greenland Vikings date to 985 CE. Archaeological evidence yields insight into how Vikings lived, yet drivers of their disappearance in the 15th century remain enigmatic. Research suggests a combination of environmental and socioeconomic factors, and the climatic shift from the Medieval Warm Period (~900 to 1250 CE) to the Little Ice Age (~1250 to 1900 CE) may have forced them to abandon Greenland. Glacial geomorphology and paleoclimate research suggest that the Southern Greenland Ice Sheet readvanced during Viking occupation, peaking in the Little Ice Age. Counterintuitively, the readvance caused sea-level rise near the ice margin due to increased gravitational attraction toward the ice sheet and crustal subsidence. We estimate ice growth in Southwestern Greenland using geomorphological indicators and lake core data from previous literature. We calculate the effect of ice growth on regional sea level by applying our ice history to a geophysical model of sea level with a resolution of ~1 km across Southwestern Greenland and compare the results to archaeological evidence. The results indicate that sea level rose up to ~3.3 m outside the glaciation zone during Viking settlement, producing shoreline retreat of hundreds of meters. Sea-level rise was progressive and encompassed the entire Eastern Settlement. Moreover, pervasive flooding would have forced abandonment of many coastal sites. These processes likely contributed to the suite of vulnerabilities that led to Viking abandonment of Greenland. Sea-level change thus represents an integral, missing element of the Viking story.

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., Competing Interests: Competing interestsThe authors declare no competing interests., (© The Author(s) 2023.)

Rapid melting of ice sheets and glaciers drives a unique geometry, or fingerprint, of sea level change. However, the detection of individual fingerprints has been challenging because of sparse observations at high latitudes and the difficulty of disentangling ocean dynamic variability from the signal. We predict the fingerprint of Greenland Ice Sheet (GrIS) melt using recent ice mass loss estimates from radar altimetry data and model reconstructions of nearby glaciers and compare this prediction to an independent, altimetry-derived sea surface height trend corrected for ocean dynamic variability in the region adjacent to the ice sheet. A statistically significant correlation between the two fields ( P < 0.001) provides an unambiguous observational detection of the near-field sea level fingerprint of recent GrIS melting in our warming world.

An accurate record of preindustrial (pre-1900 CE) sea level is necessary to contextualize modern global mean sea level (GMSL) rise with respect to natural variability. Precisely dated phreatic overgrowths on speleothems (POS) provide detailed rates of Late Holocene sea-level rise in Mallorca. Statistical analysis indicates that sea level rose locally by 0.12 to 0.31 m (95% confidence) from 3.26 to 2.84 thousand years (ka) ago (2σ) and remained within 0.08 m (95% confidence) of preindustrial levels from 2.84 ka to 1900 CE. This sea-level history is consistent with glacial isostatic adjustment models adopting relatively weak upper mantle viscosities of ~10 20 Pa s. There is virtual certainty (>0.999 probability) that the average GMSL rise since 1900 CE has exceeded even the high average rate of sea-level rise between 3.26 and 2.84 ka inferred from the POS record. We conclude that modern GMSL rise is anomalous relative to any natural variability in ice volumes over the past 4000 years.

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., (© 2021. The Authors.)

Geodetic, seismic, and geological evidence indicates that West Antarctica is underlain by low-viscosity shallow mantle. Thus, as marine-based sectors of the West Antarctic Ice Sheet (WAIS) retreated during past interglacials, or will retreat in the future, exposed bedrock will rebound rapidly and flux meltwater out into the open ocean. Previous studies have suggested that this contribution to global mean sea level (GMSL) rise is small and occurs slowly. We challenge this notion using sea level predictions that incorporate both the outflux mechanism and complex three-dimensional viscoelastic mantle structure. In the case of the last interglacial, where the GMSL contribution from WAIS collapse is often cited as ~3 to 4 meters, the outflux mechanism contributes ~1 meter of additional GMSL change within ~1 thousand years of the collapse. Using a projection of future WAIS collapse, we also demonstrate that the outflux can substantially amplify GMSL rise estimates over the next century., (Copyright © 2021 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).)

Sea-level rise due to ice loss in the Northern Hemisphere in response to insolation and greenhouse gas forcing is thought to have caused grounding-line retreat of marine-based sectors of the Antarctic Ice Sheet (AIS) 1-3 . Such interhemispheric sea-level forcing may explain the synchronous evolution of global ice sheets over ice-age cycles. Recent studies that indicate that the AIS experienced substantial millennial-scale variability during and after the last deglaciation 4-7 (roughly 20,000 to 9,000 years ago) provide further evidence of this sea-level forcing. However, global sea-level change as a result of mass loss from ice sheets is strongly nonuniform, owing to gravitational, deformational and Earth rotational effects 8 , suggesting that the response of AIS grounding lines to Northern Hemisphere sea-level forcing is more complicated than previously modelled 1,2,6 . Here, using an ice-sheet model coupled to a global sea-level model, we show that AIS dynamics are amplified by Northern Hemisphere sea-level forcing. As a result of this interhemispheric interaction, a large or rapid Northern Hemisphere sea-level forcing enhances grounding-line advance and associated mass gain of the AIS during glaciation, and grounding-line retreat and mass loss during deglaciation. Relative to models without these interactions, the inclusion of Northern Hemisphere sea-level forcing in our model increases the volume of the AIS during the Last Glacial Maximum (about 26,000 to 20,000 years ago), triggers an earlier retreat of the grounding line and leads to millennial-scale variability throughout the last deglaciation. These findings are consistent with geologic reconstructions of the extent of the AIS during the Last Glacial Maximum and subsequent ice-sheet retreat, and with relative sea-level change in Antarctica 3-7,9,10 .

Global sea level provides an important indicator of the state of the warming climate, but changes in regional sea level are most relevant for coastal communities around the world. With improvements to the sea-level observing system, the knowledge of regional sea-level change has advanced dramatically in recent years. Satellite measurements coupled with in situ observations have allowed for comprehensive study and improved understanding of the diverse set of drivers that lead to variations in sea level in space and time. Despite the advances, gaps in the understanding of contemporary sea-level change remain and inhibit the ability to predict how the relevant processes may lead to future change. These gaps arise in part due to the complexity of the linkages between the drivers of sea-level change. Here we review the individual processes which lead to sea-level change and then describe how they combine and vary regionally. The intent of the paper is to provide an overview of the current state of understanding of the processes that cause regional sea-level change and to identify and discuss limitations and uncertainty in our understanding of these processes. Areas where the lack of understanding or gaps in knowledge inhibit the ability to provide the needed information for comprehensive planning efforts are of particular focus. Finally, a goal of this paper is to highlight the role of the expanded sea-level observation network-particularly as related to satellite observations-in the improved scientific understanding of the contributors to regional sea-level change., Competing Interests: The authors declare no competing interests., (©2020. American Geophysical Union. All Rights Reserved.)

The artificial impoundment of water behind dams causes global mean sea level (GMSL) to fall as reservoirs fill but also generates a local rise in sea level due to the increased mass in the reservoir and the crustal deformation this mass induces. To estimate spatiotemporal fluctuations in sea level due to water impoundment, we use a historical data set that includes 6,329 reservoirs completed between 1900 and 2011, as well as projections of 3,565 reservoirs that are expected to be completed by 2040. The GMSL change associated with the historical data (-0.2 mm yr -1 from 1900-2011) is consistent with previous studies, but the temporal and spatial resolution allows for local studies that were not previously possible, revealing that some locations experience a sea level rise of as much as 40 mm over less than a decade. Future construction of reservoirs through ~2040 is projected to cause a GMSL fall whose rate is comparable to that of the last century (-0.3 mm yr -1 ) but with a geographic distribution that will be distinct from the last century, including a rise in sea level in more coastal areas. The analysis of expected construction shows that significant impoundment near coastal communities in the coming decades could enhance the flooding risk already heightened by global sea level rise., (©2020. The Authors.)

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)., (Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution License 4.0 (CC BY).)

Sea-level histories during the two most recent deglacial-interglacial intervals show substantial differences 1-3 despite both periods undergoing similar changes in global mean temperature 4,5 and forcing from greenhouse gases 6 . Although the last interglaciation (LIG) experienced stronger boreal summer insolation forcing than the present interglaciation 7 , understanding why LIG global mean sea level may have been six to nine metres higher than today has proven particularly challenging 2 . Extensive areas of polar ice sheets were grounded below sea level during both glacial and interglacial periods, with grounding lines and fringing ice shelves extending onto continental shelves 8 . This suggests that oceanic forcing by subsurface warming may also have contributed to ice-sheet loss 9-12 analogous to ongoing changes in the Antarctic 13,14 and Greenland 15 ice sheets. Such forcing would have been especially effective during glacial periods, when the Atlantic Meridional Overturning Circulation (AMOC) experienced large variations on millennial timescales 16 , with a reduction of the AMOC causing subsurface warming throughout much of the Atlantic basin 9,12,17 . Here we show that greater subsurface warming induced by the longer period of reduced AMOC during the penultimate deglaciation can explain the more-rapid sea-level rise compared with the last deglaciation. This greater forcing also contributed to excess loss from the Greenland and Antarctic ice sheets during the LIG, causing global mean sea level to rise at least four metres above modern levels. When accounting for the combined influences of penultimate and LIG deglaciation on glacial isostatic adjustment, this excess loss of polar ice during the LIG can explain much of the relative sea level recorded by fossil coral reefs and speleothems at intermediate- and far-field sites.

The imprint of glacial isostatic adjustment has long been recognized in shoreline elevations of oceans and proglacial lakes, but to date, its signature has not been identified in river long profiles. Here, we reveal that the buried bedrock valley floor of the upper Mississippi River exhibits a 110-m-deep, 300-km-long overdeepening that we interpret to be a partial cast of the Laurentide Ice Sheet forebulge, the ring of flexurally raised lithosphere surrounding the ice sheet. Incision through this forebulge occurred during a single glacial cycle at some time between 2.5 and 0.8 million years before present, when ice-sheet advance forced former St. Lawrence River tributaries in Minnesota and Wisconsin to flow southward. This integrated for the first time the modern Mississippi River, permanently changing continental-scale hydrology and carving a bedrock valley through the migrating forebulge with sediment-poor water. The shape of the inferred forebulge is consistent with an ice sheet ~1 km thick near its margins, similar to the Laurentide Ice Sheet at the Last Glacial Maximum, and provides evidence of the impact of geodynamic processes on geomorphology even in the midst of a stable craton.

Identifying the causes of historical trends in relative sea level-the height of the sea surface relative to Earth's crust-is a prerequisite for predicting future changes. Rates of change along the eastern coast of the USA (the US East Coast) during the past century were spatially variable, and relative sea level rose faster along the Mid-Atlantic Bight than along the South Atlantic Bight and the Gulf of Maine. Past studies suggest that Earth's ongoing response to the last deglaciation 1-5 , surface redistribution of ice and water 5-9 and changes in ocean circulation 9-13 contributed considerably to this large-scale spatial pattern. Here we analyse instrumental data 14,15 and proxy reconstructions 4,12 using probabilistic methods 16-18 to show that vertical motions of Earth's crust exerted the dominant control on regional spatial differences in relative sea-level trends along the US East Coast during 1900-2017, explaining most of the large-scale spatial variance. Rates of coastal subsidence caused by ongoing relaxation of the peripheral forebulge associated with the last deglaciation are strongest near North Carolina, Maryland and Virginia. Such structure indicates that Earth's elastic lithosphere is thicker than has been assumed in other models 19-22 . We also find a substantial coastal gradient in relative sea-level trends over this period that is unrelated to deglaciation and suggests contributions from twentieth-century redistribution of ice and water. Our results indicate that the majority of large-scale spatial variation in long-term rates of relative sea-level rise on the US East Coast is due to geological processes that will persist at similar rates for centuries.

Earth's body tide-also known as the solid Earth tide, the displacement of the solid Earth's surface caused by gravitational forces from the Moon and the Sun-is sensitive to the density of the two Large Low Shear Velocity Provinces (LLSVPs) beneath Africa and the Pacific. These massive regions extend approximately 1,000 kilometres upward from the base of the mantle and their buoyancy remains actively debated within the geophysical community. Here we use tidal tomography to constrain Earth's deep-mantle buoyancy derived from Global Positioning System (GPS)-based measurements of semi-diurnal body tide deformation. Using a probabilistic approach, we show that across the bottom two-thirds of the two LLSVPs the mean density is about 0.5 per cent higher than the average mantle density across this depth range (that is, its mean buoyancy is minus 0.5 per cent), although this anomaly may be concentrated towards the very base of the mantle. We conclude that the buoyancy of these structures is dominated by the enrichment of high-density chemical components, probably related to subducted oceanic plates or primordial material associated with Earth's formation. Because the dynamics of the mantle is driven by density variations, our result has important dynamical implications for the stability of the LLSVPs and the long-term evolution of the Earth system.

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.

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.

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 [Formula: see text]) 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 ([Formula: see text]) that 20th century GSL would have risen by less than 51% of the observed [Formula: see text] 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., Competing Interests: The authors declare no conflict of interest.

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.

For about three decades, the Global Navigation Satellite System (GNSS) has been used for high-precision positioning in scientific and engineering applications, such as deformation monitoring for seismicity and volcano eruption. Such high-precision positioning applications require millimeter-level positioning accuracy. There are many man-made and natural reflective surfaces near the GNSS receiving antennas. GNSS signals can be reflected and then arrive at the GNSS antenna. The multipath effect occurs when the direct signal is mixed with the reflected signal at the GNSS receiver. Choke-ring antennas are designed to mitigate the multipath effect of reflected signals from below the horizontal plane of the GNSS receiving antenna. Moreover, GNSS receiving antennas at network/permanent stations are usually installed on tall pillars or monuments to prevent multipath from “ground” reflected signals. However, part of the reflected signals can still arrive at the GNSS antenna center and cause multipath errors in GNSS measurements. How much can the multipath effect be on the real-time GNSS-measured displacements in studies on seismicity and volcano eruption? This work investigates the below-the-horizon multipath effect of choke-ring antennas on GNSS carrier-phase measurements. Here we show the differenced carrier-phase multipath errors of two commonly used GNSS antennas at the International GNSS Service (IGS) tracking stations can reach 8 mm, the maximum, with the mean and SD in a few millimeters at the 95% confidence level. The findings of this work should be applicable to other choke-ring antennas with similar architecture. [ABSTRACT FROM AUTHOR]

Ocean tide loading (OTL) displacements observed by the Global Navigation Satellite System (GNSS) are generally less accurate for solar (S2 and P1) and lunisolar (K2 and K1) tidal constituents than for lunar tidal constituents (M2, N2, O1 and Q1). We propose to construct the tidal admittance function (TAF) in the complex domain for tidal inference, aiming to refine noisier GNSS-observed solar and lunisolar constituents. The vertical OTL displacements are estimated from the single GNSS (GPS and Galileo) and multi-GNSS (GPS + Galileo + GLONASS) solutions, using kinematic precise point positioning under undifferenced ambiguity resolution to process 2.5 years of continuous GNSS observations from 49 global GNSS stations. Results show that TAF inferred P1 and K2 agree better with FES2014b model predictions than GNSS estimates for most GNSS stations. We find that the P1 and K2 inferences are only perturbed by less than 0.3 times the observed K1 and S2 errors, respectively. In contrast, the K1 and S2 inferences suffer from over three times the observed P1 and K2 errors. As a result, the tidal inference can only provide limited refinement of GNSS K1 and S2 estimates for a few stations. By comparing the two different TAF approaches for the tidal inference, we find that direct fitting of phase lags is susceptible to phase addition and subtraction of 360°, without disrupting our approach. In addition, the OTL displacements predicted at most GNSS stations are found to satisfy the smoothness of the tidal admittance, which is fundamental to the robustness of the tidal inference. Our approach is also applicable to improve the horizontal solar and lunisolar constituents estimated by GNSS, but should be used with caution in cases where smooth admittances may be distorted by the tidal resonance. [ABSTRACT FROM AUTHOR]

Global mean sea-level (GMSL) change can shed light on how the Earth system responds to warming. Glaciological evidence indicates that Earth's ice sheets retreated inland of early industrial (1850 CE) extents during the Holocene (11.7-0 ka), yet previous work suggests that Holocene GMSL never surpassed early industrial levels. We merge sea-level data with a glacial isostatic adjustment model ensemble and reconstructions of postglacial thermosteric sea-level and mountain glacier evolution to estimate Holocene GMSL and ice volume. We show it is likely (probability P = 0.75) GMSL exceeded early industrial levels after 7.5ka, reaching 0.24 m (−3.3 to 1.0 m, 90% credible interval) above present by 3.2ka; Antarctica was likely (P = 0.78) smaller than present after 7ka; GMSL rise by 2150 will very likely (P = 0.9) be the fastest in the last 5000 years; and by 2060, GMSL will as likely than not (P = 0.5) be the highest in 115,000 years. Sea level data and models of solid Earth effects and ocean thermal expansion are used to show that global mean sea level and Antarctic ice volume were likely higher and smaller, respectively, in the mid-late Holocene than at the preindustrial. [ABSTRACT FROM AUTHOR]

Ammonia volatilized from penguin excreta is a significant nitrogen source in Antarctic ecosystems, influencing climate through new particle formation (NPF). Freeze-thaw events can trigger ammonia emissions, but their impact on penguin-derived ammonia is understudied and overlooked in models. Here we investigate the contribution of penguins to ammonia and their climatic impacts using cruise observations and GEOS-Chem-APM simulations. High ammonia concentrations, with a maximum exceeding 7000 ng/m3, were observed over the Southern Ocean and Prydz Bay, driven by air masses from penguin colonies. Simulations showed that incorporating freeze-thaw impact improves model performance, with penguin-derived ammonia emissions enhanced by up to 20-fold and reaching a total of 49 Gg across Antarctica in November. Elevated ammonia increased simulated secondary particle number concentrations by 30−300% through NPF, enhancing simulated cloud droplet number concentrations by 10−20% and altering cloud properties. This study underscores the importance of incorporating penguin emissions into models, particularly during freeze-thaw events. [ABSTRACT FROM AUTHOR]

Branchinecta gaini Daday, 1910 is the southernmost extant Anostracan species occurring in various localities across Maritime Antarctica. Since this region was almost completely glaciated during the Last Glacial Maximum, the identification of potential glacial refugia and the colonization sequence of freshwater habitats within Maritime Antarctica remains primary biogeographic objective for understanding the evolution of the Antarctic environment. To evaluate the history of B. gaini in Maritime Antarctica, we performed phylogeographic analyses of B. gaini and its closest Patagonian congener Branchinecta granulosa Daday, 1902. We sampled 47 populations from all three Antarctic bioregions where B. gaini occurs – South Orkney Islands, North-west and North-east Antarctic Peninsula. B. granulosa was represented by 8 populations from Patagonia. Molecular analyses of two mitochondrial (16S, COI) and one nuclear (ITS2) marker showed low overall mitochondrial variability and the ambiguous nature of ITS2 variability. The species assessment revealed insufficient genetic differentiation between B. gaini and B. granulosa to consider them two separate species. Therefore, the widely accepted idea that there is a species of Branchinecta endemic to Antarctica and sub-Antarctic islands should be omitted, and B. gaini should be synonymized with B. granulosa. Molecular clock analysis of COI for two substitution rates dated the separation of B. gaini from B. granulosa to the Pleistocene (170 – 502 ka BP and 17 – 50 ka BP, respectively). Therefore, Antarctic populations of B. gaini separated from Patagonian B. granulosa before or during the Last Glacial Maximum, indicating an older association with the southern continent than previously anticipated. The geographic distribution of haplotypes further suggests that B. gaini most likely survived at least the last glacial period in an as-yet unknown refugium in Antarctica. [ABSTRACT FROM AUTHOR]

Estimating and accounting for twentieth-century global mean sea level (GMSL) rise is critical to characterizing current and future human-induced sea-level change. Several previous analyses of tide gauge records--employing different methods to accommodate the spatial sparsity and temporal incompleteness of the data and to constrain the geometry of long-term sea-level change--have concluded that GMSL rose over the twentieth century at a mean rate of 1.6 to 1.9 millimetres per year. Efforts to account for this rate by summing estimates of individual contributions from glacier and ice-sheet mass loss, ocean thermal expansion, and changes in land water storage fall significantly short in the period before 1990. The failure to close the budget of GMSL during this period has led to suggestions that several contributions may have been systematically underestimated. However, the extent to which the limitations of tide gauge analyses have affected estimates of the GMSL rate of change is unclear. Here we revisit estimates of twentieth-century GMSL rise using probabilistic techniques and find a rate of GMSL rise from 1901 to 1990 of 1.2 ± 0.2 millimetres per year (90% confidence interval). Based on individual contributions tabulated in the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, this estimate closes the twentieth-century sea-level budget. Our analysis, which combines tide gauge records with physics-based and model-derived geometries of the various contributing signals, also indicates that GMSL rose at a rate of 3.0 ± 0.7 millimetres per year between 1993 and 2010, consistent with prior estimates from tide gauge records.The increase in rate relative to the 1901-90 trend is accordingly larger than previously thought; this revision may affect some projections of future sea-level rise.

At the Last Glacial Maximum (LGM), about 21,000 years before present, land-based ice sheets held enough water to reduce global mean sea level by 130 metres. Yet after decades of study, major uncertainties remain as to the distribution of that ice. Here we test four reconstructions of North American deglacial ice-sheet history by quantitatively connecting them to high-resolution oxygen isotope (δ(18)O) records from the Gulf of Mexico using a water mixing model. For each reconstruction, we route meltwater and seasonal runoff through the time-evolving Mississippi drainage basin, which co-evolves with ice geometry and changing topography as ice loads deform the solid Earth and produce spatially variable sea level in a process known as glacial isostatic adjustment. The δ(18)O records show that the Mississippi-drained southern Laurentide ice sheet contributed only 5.4 ± 2.1 metres to global sea level rise, of which 0.66 ± 0.07 metres were released during the meltwater pulse 1A event 14,650-14,310 years before present, far less water than previously thought. In contrast, the three reconstructions based on glacial isostatic adjustment overpredict the δ(18)O-based post-LGM meltwater volume by a factor of 1.6 to 3.6. The fourth reconstruction, which is based on ice physics, has a low enough Mississippi-routed meltwater discharge to be consistent with δ(18)O constraints, but also contains the largest LGM North American ice volume. This suggests that modelling based on ice physics may be the best way of matching isotopic records while also sequestering enough water in the North American ice sheets to match the observed LGM sea level fall.

Sedimentary rocks from Virginia through Florida record marine flooding during the mid-Pliocene. Several wave-cut scarps that at the time of deposition would have been horizontal are now draped over a warped surface with a maximum variation of 60 meters. We modeled dynamic topography by using mantle convection simulations that predict the amplitude and broad spatial distribution of this distortion. The results imply that dynamic topography and, to a lesser extent, glacial isostatic adjustment account for the current architecture of the coastal plain and proximal shelf. This confounds attempts to use regional stratigraphic relations as references for longer-term sea-level determinations. Inferences of Pliocene global sea-level heights or stability of Antarctic ice sheets therefore cannot be deciphered in the absence of an appropriate mantle dynamic reference frame.

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.

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.

Deciphering the evolution of global climate from the end of the Last Glacial Maximum approximately 19 ka to the early Holocene 11 ka presents an outstanding opportunity for understanding the transient response of Earth's climate system to external and internal forcings. During this interval of global warming, the decay of ice sheets caused global mean sea level to rise by approximately 80 m; terrestrial and marine ecosystems experienced large disturbances and range shifts; perturbations to the carbon cycle resulted in a net release of the greenhouse gases CO(2) and CH(4) to the atmosphere; and changes in atmosphere and ocean circulation affected the global distribution and fluxes of water and heat. Here we summarize a major effort by the paleoclimate research community to characterize these changes through the development of well-dated, high-resolution records of the deep and intermediate ocean as well as surface climate. Our synthesis indicates that the superposition of two modes explains much of the variability in regional and global climate during the last deglaciation, with a strong association between the first mode and variations in greenhouse gases, and between the second mode and variations in the Atlantic meridional overturning circulation.

Contentious observations of Pleistocene shoreline features on the tectonically stable islands of Bermuda and the Bahamas have suggested that sea level about 400,000 years ago was more than 20 metres higher than it is today. Geochronologic and geomorphic evidence indicates that these features formed during interglacial marine isotope stage (MIS) 11, an unusually long interval of warmth during the ice age. Previous work has advanced two divergent hypotheses for these shoreline features: first, significant melting of the East Antarctic Ice Sheet, in addition to the collapse of the West Antarctic Ice Sheet and the Greenland Ice Sheet; or second, emplacement by a mega-tsunami during MIS 11 (ref. 4, 5). Here we show that the elevations of these features are corrected downwards by ∼10 metres when we account for post-glacial crustal subsidence of these sites over the course of the anomalously long interglacial. On the basis of this correction, we estimate that eustatic sea level rose to ∼6-13 m above the present-day value in the second half of MIS 11. This suggests that both the Greenland Ice Sheet and the West Antarctic Ice Sheet collapsed during the protracted warm period while changes in the volume of the East Antarctic Ice Sheet were relatively minor, thereby resolving the long-standing controversy over the stability of the East Antarctic Ice Sheet during MIS 11.

The timing of the last maximum extent of the Antarctic ice sheets relative to those in the Northern Hemisphere remains poorly understood. We develop a chronology for the Weddell Sea sector of the East Antarctic Ice Sheet that, combined with ages from other Antarctic ice-sheet sectors, indicates that the advance to and retreat from their maximum extent was within dating uncertainties synchronous with most sectors of Northern Hemisphere ice sheets. Surface climate forcing of Antarctic mass balance would probably cause an opposite response, whereby a warming climate would increase accumulation but not surface melting. Our new data support teleconnections involving sea-level forcing from Northern Hemisphere ice sheets and changes in North Atlantic deep-water formation and attendant heat flux to Antarctic grounding lines to synchronize the hemispheric ice sheets.

With polar temperatures approximately 3-5 degrees C warmer than today, the last interglacial stage (approximately 125 kyr ago) serves as a partial analogue for 1-2 degrees C global warming scenarios. Geological records from several sites indicate that local sea levels during the last interglacial were higher than today, but because local sea levels differ from global sea level, accurately reconstructing past global sea level requires an integrated analysis of globally distributed data sets. Here we present an extensive compilation of local sea level indicators and a statistical approach for estimating global sea level, local sea levels, ice sheet volumes and their associated uncertainties. We find a 95% probability that global sea level peaked at least 6.6 m higher than today during the last interglacial; it is likely (67% probability) to have exceeded 8.0 m but is unlikely (33% probability) to have exceeded 9.4 m. When global sea level was close to its current level (>or=-10 m), the millennial average rate of global sea level rise is very likely to have exceeded 5.6 m kyr(-1) but is unlikely to have exceeded 9.2 m kyr(-1). Our analysis extends previous last interglacial sea level studies by integrating literature observations within a probabilistic framework that accounts for the physics of sea level change. The results highlight the long-term vulnerability of ice sheets to even relatively low levels of sustained global warming.

We used 5704 14C, 10Be, and 3He ages that span the interval from 10,000 to 50,000 years ago (10 to 50 ka) to constrain the timing of the Last Glacial Maximum (LGM) in terms of global ice-sheet and mountain-glacier extent. Growth of the ice sheets to their maximum positions occurred between 33.0 and 26.5 ka in response to climate forcing from decreases in northern summer insolation, tropical Pacific sea surface temperatures, and atmospheric CO2. Nearly all ice sheets were at their LGM positions from 26.5 ka to 19 to 20 ka, corresponding to minima in these forcings. The onset of Northern Hemisphere deglaciation 19 to 20 ka was induced by an increase in northern summer insolation, providing the source for an abrupt rise in sea level. The onset of deglaciation of the West Antarctic Ice Sheet occurred between 14 and 15 ka, consistent with evidence that this was the primary source for an abrupt rise in sea level approximately 14.5 ka.

Recent projections of sea-level rise after a future collapse of the West Antarctic Ice Sheet (for example, the Fourth Intergovernmental Panel on Climate Change Assessment Report) assume that meltwater will spread uniformly (that is, eustatically) across the oceans once marine-based sectors of the West Antarctic are filled. A largely neglected 1977 study predicted that peak values would be 20% higher than the eustatic in the North Pacific and 5 to 10% higher along the U.S. coastline. We show, with use of a state-of-the-art theory, that the sea-level rise in excess of the eustatic value will be two to three times higher than previously predicted for U.S. coastal sites.

A suite of observations suggests that the northern plains of Mars, which cover nearly one third of the planet's surface, may once have contained an ocean. Perhaps the most provocative evidence for an ancient ocean is a set of surface features that ring the plains for thousands of kilometres and that have been interpreted as a series of palaeoshorelines of different age. It has been shown, however, that topographic profiles along the putative shorelines contain long-wavelength trends with amplitudes of up to several kilometres, and these trends have been taken as an argument against the martian shoreline (and ocean) hypothesis. Here we show that the long-wavelength topography of the shorelines is consistent with deformation caused by true polar wander--a change in the orientation of a planet with respect to its rotation pole--and that the inferred pole path has the geometry expected for a true polar wander event that postdates the formation of the massive Tharsis volcanic rise.

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.

The planform is a defining feature of mantle convection. It can be gleaned from the stratigraphic record by mapping the continent-scale distribution of hiatus and no hiatus surfaces serving as a proxy for high and low dynamic topography. We carry this out for all continents apart from Antarctica for eight geological series since the Upper Jurassic, showing that: (i) the planform as indicated by our maps contains wavelengths of the order of 1000 km, smaller than the convective scales implied by the geoid. (ii) The planform changes on timescales of geological series (10–20 Myrs), smaller than the mantle transit time. (iii) Flood basalt eruptions are frequently preceded by hiatus surfaces. (iv) Some hiatus surfaces are not linked to any known plume, potentially reflecting the lateral transport of material in the asthenosphere. Our results reveal the importance of mantle viscosity stratification in shaping the convective planform and the resulting dynamic topography. Geodynamic Earth models should aim to reproduce the global characteristics of our maps, as well as specific regional events identified in this work. Finally, we separate the effects of sea-level variation from regional changes in base level induced by dynamic topography by contrasting the stratigraphic evolution of different regions. [ABSTRACT FROM AUTHOR]

Previous studies revealed extensive genetic introgression between Ovis species, which affects genetic adaptation and morphological traits. However, the exact evolutionary scenarios underlying the hybridization between sheep and allopatric wild relatives remain unknown. To address this problem, we here integrate the reference genomes of several ovine and caprine species: domestic sheep, argali, bighorn sheep, snow sheep, and domestic goats. Additionally, we use 856 whole genomes representing 169 domestic sheep populations and their six wild relatives: Asiatic mouflon, urial, argali, snow sheep, thinhorn sheep, and bighorn sheep. We implement a comprehensive set of analyses to test introgression among these species. We infer that the argali lineage originated ∼3.08 to 3.35 Mya and hybridized with the ancestor of Pachyceriforms (e.g. bighorn sheep and snow sheep) at ∼1.56 Mya. Previous studies showed apparent introgression from North American Pachyceriforms into the Bashibai sheep, a Chinese native sheep breed, despite of their wide geographic separation. We show here that, in fact, the apparent introgression from the Pachyceriforms into Bashibai can be explained by the old introgression from Pachyceriforms into argali and subsequent recent introgression from argali into Bashibai. Our results illustrate the challenges of estimating complex introgression histories and provide an example of how indirect and direct introgression can be distinguished. [ABSTRACT FROM AUTHOR]

Accelerating sea level rise (SLR) and changing storm patterns will increasingly expose barrier islands to coastal hazards, including flooding, erosion, and rising groundwater tables. We assess the exposure of Cape Lookout National Seashore, a barrier island system in North Carolina (USA), to projected SLR and storm hazards over the twenty-first century. We estimate that with 0.5 m of SLR, 47% of current subaerial barrier island area would be flooded daily, and the 1-year return period storm would flood 74%. For 20-year return period storms, over 85% is projected to be flooded for any SLR. The modelled groundwater table is already shallow (< 2 m deep), and while projected to shoal to the land surface with SLR, marine flooding is projected to overtake areas with emergent groundwater. Projected shoreline retreat reaches an average of 178 m with 1 m of SLR and no interventions, which is over 60% of the current island width at narrower locations. Compounding these hazards is subsidence, with one-third of the study area currently lowering at > 2 mm/yr. Our results demonstrate the difficulty of managing natural barrier systems such as those managed by federal park systems tasked with maintaining natural ecosystems and protecting cultural resources. [ABSTRACT FROM AUTHOR]

At least two global "Snowball Earth" glaciations occurred during the Neoproterozoic Era (1000-538.8 million years ago). Post-glacial surface environments during this time are recorded in cap carbonates: layers of limestone or dolostone that directly overlie glacial deposits. Postulated environmental conditions that created the cap carbonates lack consensus largely because single hypotheses fail to explain the cap carbonates' global mass, depositional timescales, and geochemistry of parent waters. Here, we present a global geologic carbon cycle model before, during, and after the second glaciation (i.e. the Marinoan) that explains cap carbonate characteristics. We find a three-stage process for cap carbonate formation: (1) low-temperature seafloor weathering during glaciation generates deep-sea alkalinity; (2) vigorous post-glacial continental weathering supplies alkalinity to a carbonate-saturated freshwater layer, rapidly precipitating cap carbonates; (3) mixing of post-glacial meltwater with deep-sea alkalinity prolongs cap carbonate deposition. We suggest how future geochemical data and modeling refinements could further assess our hypothesis. A new formation mechanism is proposed to explain the presence and characteristics of "cap carbonates": enigmatic rocks leftover from Snowball Earth events, during which Earth was globally ice-covered for millions of years in the distant past. [ABSTRACT FROM AUTHOR]

Ganymede has an ancient impact structure called a furrow system. The furrow system is the largest impact structure in the outer solar system, and the impact should have significantly affected Ganymede's early history; however, its effects are poorly understood. No attention has been given to the center of the furrow system coinciding with Ganymede's tidal axis, indicating that mass redistribution induced by the furrow-forming impact caused a reorientation (true polar wander) of Ganymede. We propose that the impact ejecta created a mass anomaly that reoriented the impact site toward the tidal axis. We found that an impactor with a radius of 150 km and an incidence angle between 60° and 90° most accurately reproduces the current location of the furrow system. We predict that future explorations would reveal remnant topographic profiles or gravity anomalies associated with the furrow-forming impact and reorientation. Additionally, various possible explanations for the reorientation of Ganymede, such as an impactor-origin mascon beneath the basin or a thickness variation in the lithosphere, should be studied. [ABSTRACT FROM AUTHOR]

We developed automated ablation stakes to measure colocated in situ changes in relative glacier-surface elevation and climatological drivers of ablation. The designs, refined over 10 years of development and deployments, implement open-source hardware and common building materials. The ablation stakes record distance to the snow/ice surface, air temperature and relative humidity every 1–15 min. Using these high-frequency data, we demonstrate that melt factors calculated using standard melt-rate vs temperature regressions converge over averaging windows of approximately 12 h or greater. Furthermore, we evaluate an integral approach to estimating temperature-index melt factors for ablation. In a test case on Glaciar Perito Moreno, Argentina, this integral approach reveals an overall positive-degree-day melt factor of 7.5 mm w.e. $^\circ$ C−1 d−1. We describe four deployments with iteratively improved designs and provide a list of materials required to construct an automated ablation stake. [ABSTRACT FROM AUTHOR]