Meltwater Pulse Research Articles

10Be exposure age dating of Late Quaternary relative sea level changes and deglaciation of W Jura and NE Islay, Scottish Inner Hebrides

ABSTRACTNew 10Be exposure age dating and geomorphological mapping of emerged shoreline features in W Jura and NE Islay throw new light on the regional pattern of ice sheet deglaciation and late-glacial relative sea level change. We conclude that the oldest and highest emerged shorelines in this area were produced ~15.7–16.3 ka, shortly after ice sheet deglaciation ~16.5 ka. It is envisaged that the first incursion of marine waters into coastal areas took place close to a former ice sheet margin that oscillated in position across this part of the Scottish Inner Hebrides. The first evidence of late-glacial marine sedimentation following deglaciation consists of emerged marine terrace fragments and unvegetated gravel beach ridges, the former represented by a prominent glacio-isostatically tilted shoreline that declines in altitude NE to SW, from ~40 m above ordnance datum (OD) in NW Jura to ~19 m OD in central Islay. In W Jura, north of Loch Tarbert, spectacular staircases of up to 55 unvegetated gravel beach ridges were formed shortly after regional deglaciation, possibly within 1 ka. A preliminary estimate of the average rate of relative sea level lowering across W Jura between deglaciation and the Younger Dryas is in the order of ~7 mmyr−1. Geomorphological evidence from Shian Bay, W Jura, indicates a truncation of the late-glacial beach ridge staircases by a large 480-m-long beach ridge (the Colonsay Ridge) at ~14.9 ka, when former relative sea level was at ~18 m OD. This ridge may represent the product of either a stillstand in the progressive lowering of relative sea level during the late-glacial or a reversal. This raises the intriguing possibility of an association between ridge formation and the timing of the well-established global meltwater pulse 1A between ~14.65 and ~14.8 ka.

Influence of Lateglacial meltwater influx from alpine glaciers on the aquatic ecosystem of Kanas Lake in the Altai Mountains, China

Glacier retreat in high altitude areas triggered by anthropogenic warming exerts a profound impact on the ecology of glacier-fed lakes. In this paper, we report on the results of a multi-proxy analysis of a Lateglacial (13.3–11.3 cal kyr B.P.) core obtained from the glacier-fed Kanas Lake in the Altai Mountains, China, to understand the responses of aquatic ecosystems to meltwater influxes. Core intervals with high sand fractions and enriched silicon dioxide and zirconium (~11.5, 11.7, 11.8, 12.2 and 13.0 cal kyrs B.P.) reflect intense meltwater pulses from alpine glaciers. In the core, these pulses generally coincide with low abundance of littoral cladocerans, high abundance of planktonic Pediastrum, and high phosphorus concentrations. Our results suggest that intense meltwater influxes not only brought nutrients (nitrogen and phosphorus) into the lake and promoted phytoplankton growth, they also generated cold, turbid water masses and restricted the growth of benthic algae and invertebrates. Therefore, we infer that future anthropogenic warming may lead to significant changes in the structure and function of aquatic ecosystems in glacier-fed lakes worldwide.

Rapid northern hemisphere ice sheet melting during the penultimate deglaciation

The rate and consequences of future high latitude ice sheet retreat remain a major concern given ongoing anthropogenic warming. Here, new precisely dated stalagmite data from NW Iberia provide the first direct, high-resolution records of periods of rapid melting of Northern Hemisphere ice sheets during the penultimate deglaciation. These records reveal the penultimate deglaciation initiated with rapid century-scale meltwater pulses which subsequently trigger abrupt coolings of air temperature in NW Iberia consistent with freshwater-induced AMOC slowdowns. The first of these AMOC slowdowns, 600-year duration, was shorter than Heinrich 1 of the last deglaciation. Although similar insolation forcing initiated the last two deglaciations, the more rapid and sustained rate of freshening in the eastern North Atlantic penultimate deglaciation likely reflects a larger volume of ice stored in the marine-based Eurasian Ice sheet during the penultimate glacial in contrast to the land-based ice sheet on North America as during the last glacial.

Relative sea level and abrupt mass unloading in Barbados during the Holocene

Relative sea level records climatic change as well as vertical land movement. In Barbados, uplift variation is necessary to interpret one of the most complete coral reef records. Here we show that an abrupt mass unloading of 30 ± 10 km3 caused an uplift variation of ~0.45 ± 0.15 mm/yr using a modelling approach. Simulations have been conducted for different volumes and elastic thicknesses. Isostatic adjustment in relation with an abrupt mass unloading explains the observed uplift rate increased from 0.34 mm/yr to 0.8 mm/yr that occurred 11.2 kyr ago. The reconstructed sea-level curve highlights a sea-level jump of 4.8 m, with a delay of 150 yr from the termination of Younger Dryas cold event and 300 yr before the abrupt mass unloading. This sea-level jump corresponds to meltwater pulse MWP-1B and is not an artefact. A stagnation of 500 yr occurred from 12 to 11.5 kyr BP. Relative sea level records are useful to detect past landslides and erosion. Accurate analysis and reconstruction of sea-level permits to determine sea-level abrupt rise caused by climate warming during the last thousand years.

Sedimentological and morphological evidences of Meltwater Pulse 1B in the Southwestern Atlantic Margin

During the post-LGM deglaciation, global temperatures and rates of sea-level rise were influenced by centennial events driven by climate and ice volume variations, leading to cooling/warming climate stages, sea-level slowstands and meltwater pulses (MWP). Sea level fluctuations have shaped the shelf morphology by drowning coastal features and imposing shifts in coastal/shelf sedimentary environments worldwide. The influence of meltwater pulses in the South Atlantic is still not well established. Thus, here, we investigate morphological and sedimentological records of MWP-1B along two distinct areas of the Eastern Brazilian shelf (Espirito Santo and Abrolhos shelf), far-field sites in the western South Atlantic. Sedimentary core analysis, C14 dating, and multibeam data are combined to provide stratigraphic and morphological evidences of sea-level rise acceleration at the end of Younger Dryas, in the South Atlantic. Results of core analysis indicate a coastal environment around 12 kyrs BP and a major shift in sedimentary facies after 11.5 kyrs BP. In both cores, Unit 1 comprises terrigenous (low carbonate content) and mixed facies, and are interpreted as terrigenous/coastal/estuarine facies. Unit 2 is characterized by carbonate/mixed facies, interpreted as open marine deposits. Terrigenous and mixed estuarine/coastal facies (Unit 1) are replaced by open marine facies (unit 2) at 57.3 m in the Abrolhos Depression core (ABDR) and at 47.6 m in the Espirito Santo paleovalley shelf (PVS core). Morphological evidence shows a paleoshoreline at 60 m water depth associated with five drowned and partially filled incised shelf valleys. A depth x age analysis indicates that the sedimentary environment shift and the depth of paleoshorelines are associated with sea-level slowstand during Younger Dryas, followed by a sea-level rise acceleration, representing MWP-1B. These evidences in the western South Atlantic are combined with other morphological and stratigraphic results presented for South Hemisphere far field sites, indicating that the 60 m isobath is an average water depth marking sea-level stabilization during Younger Dryas, and that MWP-1B is potentially a global event.

Rapid sea-level rise during the first phase of the Littorina transgression in the western Baltic Sea

Glacio-isostatic adjustment (GIA) corrected limiting non-marine and marine sea-level data from the western Baltic Sea indicate a rapid relative sea-level (RSL) rise of 18 m for the first phase of the Littorina transgression using a linear age-elevation model. This accelerated RSL rise occurred in the time period 8.57 to 8.0 ka before present (BP) with an RSL rise rate of 31.5 mm/year. RSL rose from 28 m to 10 m below the mean relative sea level (bmrsl) during the time period of less than 600 years. The first phase of Littorina transgression documents the last big meltwater pulse (MWP), which has been reported also from various other regions in the world as a consequence of the final decay of the Laurentide ice sheet and additional meltwater sources from Antarctica.

When did Goliath collapse? Dating the recent events of a major submarine slide complex, Southeastern Mediterranean margin

Unraveling the failure history of submarine slides and delineating its environmental controls are crucial for practical hazard assessments, but are commonly challenging and costly. Here we address this challenge by investigating the Goliath slide complex, one of the largest documented slide scars on the ~200 km long continental slope offshore Israel. Geophysical evidence suggest that the Goliath complex is a relatively young feature. However, despite its proximity to submarine infrastructure and to one of the most populated coastal areas in the Mediterranean, its recent temporal development was never determined.Based on a detailed analysis of sub-meter resolution multi-channel seismic profiles, four sediment cores were retrieved from the northern head scarp and toe domains of the Goliath slide complex. Selective sedimentological and chronological analyses were carried out following a detailed computed tomography analysis of the cores. The results reveal two generations of deposits on the northern head scarp of Goliath, separated by a detachment surface. Two 14C ages of 7.46 ± 0.13 cal ka BP (core PHS-1) and 7.38 ± 0.16 cal ka BP (core PHS-5) were obtained for the sediments immediately overlying this surface, implying that the Goliath northern head scarp was formed, at least in part, during a substantial ~7.4 cal ka BP event. A ~3.7 kyrs hiatus within the late Holocene sedimentary sequence of core PHS-5 may suggest an additional, more recent and presumably smaller, slide event that occurred at ~2.4 ka BP. The sedimentary succession from Goliath's toe domain (core PTL-3), located ~32 km downslope from the head scarp, show an upper 1.2 m long continuous undisturbed sequence dated to the last ~14 ka BP and including sapropel S1 deposits. This sequence overlies three disturbed units, interpreted as mass transport deposits (MTDs) containing various deformation features, large-scale burrows, sharp contacts and fold structures. A 14C age dates the base of the hemipelagic sequence to 13.91 ± 0.23 cal ka BP, representing the minimum age of the slide event that deposited this MTD. Ages within the MTD units are substantially older (18.21 ± 0.27 to 28.66 ± 0.33 cal ka BP) than in the overlying sequence and an age-inversion within the MTDs supports their interpretation as transported sediments.The timing of the prominent slide events in the Goliath complex corresponds with the late stages of meltwater pulses 1A (~14.6–13.7 ka BP) and 1C (~8.2–7.5 ka BP), periods of accelerated sea level rise and probably increased sediments supply. Our results suggest that climate-driven factors likely contributed significantly to slope instability in this area, while earthquakes on the Dead Sea Transform may have been the triggering factor for the Goliath Slide events.

Double-plunge structure of the East Asian summer monsoon during Heinrich stadial 1 recorded in Xianyun Cave, southeastern China

The interhemispheric anti-phase relationship of monsoon precipitation on the millennial timescale has been widely verified and is generally attributed to the migration of the mean position of the Intertropical Convergence Zone (ITCZ). However, previous high-resolution records showed that the interhemispheric monsoon climate could have been decoupled within Heinrich stadial 1 (HS1), creating a two-fold subdivision in the Asian summer monsoon and a double-plunge structure in the South American summer monsoon. In this study, we present a high-resolution and precisely dated composite stalagmite oxygen-isotope record exhibiting a clear double-plunge structure during HS1 from southeastern China, the frontal zone predominantly influenced by the East Asian summer monsoon (EASM). Our new record characterizes a weak EASM period corresponding to HS1, divided into two weak periods (HS1c, 17.89–17.01 ka B.P.; HS1a, 15.94–14.76 ka B.P.) by a relatively strong phase (HS1b, 17.01–15.94 ka B.P.). Comparisons with other records from global monsoon systems suggest the double-plunge structure of HS1 is widespread in the low latitudes of monsoon domains, and the interhemispheric anti-phase relationship of the monsoonal hydroclimate still holds on the sub-millennial timescale. This structure could be closely linked with two meltwater pulses in the North Atlantic, separated by a weak period and a brief rerouting of the main freshwater input produced by the Laurentide ice sheet from subpolar North Atlantic to the Gulf of Mexico. The change in the interhemispheric thermal gradient during HS1b likely pulled/pushed the ITCZ northward briefly and induced the synchronous and reverse change in the interhemispheric monsoon climate.

A statistical framework for integrating nonparametric proxy distributions into geological reconstructions of relative sea level

Abstract. Robust, proxy-based reconstructions of relative sea-level (RSL) change are critical to distinguishing the processes that drive spatial and temporal sea-level variability. The relationships between individual proxies and RSL can be complex and are often poorly represented by traditional methods that assume Gaussian likelihood distributions. We develop a new statistical framework to estimate past RSL change based on nonparametric, empirical modern distributions of proxies in relation to RSL, applying the framework to corals and mangroves as an illustrative example. We validate our model by comparing its skill in reconstructing RSL and rates of change to two previous RSL models using synthetic time-series datasets based on Holocene sea-level data from South Florida. The new framework results in lower bias, better model fit, and greater accuracy and precision than the two previous RSL models. We also perform sensitivity tests using sea-level scenarios based on two periods of interest – meltwater pulses (MWPs) and the Holocene – to analyze the sensitivity of the statistical reconstructions to the quantity and precision of proxy data; we define high-precision indicators, such as mangroves and the reef-crest coral Acropora palmata, with 2σ vertical uncertainties within ± 3 m and lower-precision indicators, such as Orbicella spp., with 2σ vertical uncertainties within ± 10 m. For reconstructing rapid rates of change in RSL of up to ∼ 40 m kyr−1, such as those that may have characterized MWPs during deglacial periods, we find that employing the nonparametric model with 5 to 10 high-precision data points per kiloyear enables us to constrain rates to within ± 3 m kyr−1 (1σ). For reconstructing RSL with rates of up to ∼ 15 m kyr−1, as observed during the Holocene, we conclude that employing the model with 5 to 10 high-precision (or a combination of high- and low-precision) data points per kiloyear enables precise estimates of RSL within ±∼ 2 m (2σ) and accurate RSL reconstructions with errors ≲ 0.7 m. Employing the nonparametric model with only lower-precision indicators also produces fairly accurate estimates of RSL with errors ≲1.50 m, although with less precision, only constraining RSL to ±∼ 3–4 m (2σ). Although the model performs better than previous models in terms of bias, model fit, accuracy, and precision, it is computationally expensive to run because it requires inverting large matrices for every sample. The new model also provides minimal gains over similar models when a large quantity of high-precision data are available. Therefore, we recommend incorporating the nonparametric likelihood distributions when no other information (e.g., reef facies or epibionts indicative of shallow-water environments to refine coral elevational uncertainties) or no high-precision data are available at a location or during a given time period of interest.

Local oceanic CO<sub>2</sub> outgassing triggered by terrestrial carbon fluxes during deglacial flooding

Abstract. Exchange of carbon between the ocean and the atmosphere is a key process that influences past climates via glacial–interglacial variations of the CO2 concentration. The melting of ice sheets during deglaciations induces a sea level rise which leads to the flooding of coastal land areas, resulting in the transfer of terrestrial organic matter to the ocean. However, the consequences of such fluxes on the ocean biogeochemical cycle and on the uptake and release of CO2 are poorly constrained. Moreover, this potentially important exchange of carbon at the land–sea interface is not represented in most Earth system models. We present here the implementation of terrestrial organic matter fluxes into the ocean at the transiently changing land–sea interface in the Max Planck Institute for Meteorology Earth System Model (MPI-ESM) and investigate their effect on the biogeochemistry during the last deglaciation. Our results show that during the deglaciation, most of the terrestrial organic matter inputs to the ocean occurs during Meltwater Pulse 1a (between 15–14 ka) which leads to the transfer of 21.2 Gt C of terrestrial carbon (mostly originating from wood and humus) to the ocean. Although this additional organic matter input is relatively small in comparison to the global ocean inventory (0.06 %) and thus does not have an impact on the global CO2 flux, the terrestrial organic matter fluxes initiate oceanic outgassing in regional hotspots like in Indonesia for a few hundred years. Finally, sensitivity experiments highlight that terrestrial organic matter fluxes are the drivers of oceanic outgassing in flooded coastal regions during Meltwater Pulse 1a. Furthermore, the magnitude of outgassing is rather insensitive to higher carbon-to-nutrient ratios of the terrestrial organic matter. Our results provide a first estimate of the importance of terrestrial organic matter fluxes in a transient deglaciation simulation. Moreover, our model development is an important step towards a fully coupled carbon cycle in an Earth system model applicable to simulations at glacial–interglacial cycles.

Slow and soft passage through tipping point of the Atlantic Meridional Overturning Circulation in a changing climate

Paleo-proxy records suggest that the Atlantic Meridional Overturning Circulation (AMOC) exhibits a threshold for an abrupt change, a so-called tipping point. A classical bifurcation theory, a basis of the tipping dynamics of AMOC implicitly assumes that the tipping point is fixed. However, when a system is subjected to time-varying forcing (e.g., AMOC exposed to ice meltwater) an actual tipping point can be overshot due to delayed tipping, referred to as the slow passage effect. Here, using an Earth system model of intermediate complexity and a low-order model with freshwater forcing, we show that the tipping point of AMOC is largely delayed by the slow passage effect. It causes a large tipping lag of up to 1300 years, and strongly relaxes the abruptness of tipping as well. We further demonstrate that the tipping modulation can actively occur in past, present, and future climates by quantifying the effect during Dansgaard-Oeschger events, meltwater pulse 1A (MWP-1A), and current Greenland ice sheet melting. The suggested slow passage effect may explain the observed lagged AMOC collapse to MWP-1A of about 1000 years and provides implications tipping risk in the future.

Morphology and Late Pleistocene–Holocene sedimentation of the Strait of Istanbul (Bosphorus): a review

Abstract The Bosphorus (Istanbul) Strait is a natural strait that connects the Black Sea with the Aegean Sea via the Sea of Marmara and the Dardanelles Strait. It is a 31 km-long and 3.5 km-wide winding channel, with an irregular bottom morphology. It has depressions up to 110 m deep, and two sills with depths of 35 and 58 m in the south and north, respectively. Presently, a two-layer water exchange exists through the strait, with the Mediterranean and Black Sea waters forming the lower and upper layers, respectively. The Bosphorus channel extends as shelf valleys on the Black Sea and Sea of Marmara shelves. However, it operated as a river valley or an estuary during the stadial lowstand periods. The infill sedimentary succession of the Bosphorus channel is up to 100 m thick above the Paleozoic–Cretaceous basement with an irregular topography. The oldest sediments are sandy to muddy fluvial–lacustrine facies of late Pleistocene age, which are preserved only in up to 160 m-deep scoured depressions of the basement. They are overlain by mid–late Holocene estuarine–marine shelly sandy to muddy sediments with patches of bioherms and shelly lag deposits. The Bosphorus outlet areas of the Black Sea and Sea of Marmara are characterized by a submarine fan and a shelf valley, respectively. The fan system in the Black Sea started depositing c. 900 years after the initial vigorous marine water incursion at c. 8.4 14 C ka BP. On the Marmara shelf, extension of the Bosphorus channel is a sinuous shelf valley with a channel–levee complex that was deposited by the Black Sea outflow during 11–10 14 C ka BP. Catastrophic floodings of the Sea of Marmara by torrential Black Sea outflows during the Greenland Interstadial melt-water pulses, as well as the strong Mediterranean current towards the Black Sea during the interglacial periods, were responsible for carving out the Bosphorus channel and the shelf valleys, as well as removing the sediments belonging to the earlier periods.

Glacial isostatic adjustment near the center of the former Patagonian Ice Sheet (48°S) during the last 16.5 kyr

Our understanding of glacial isostatic rebound across Patagonia is highly limited, despite its importance to constrain past ice volume estimates and better comprehend relative sea-level variations. With this in mind, our research objective is to reconstruct the magnitude and rate of Late Glacial and Holocene glacial isostatic adjustment near the center of the former Patagonian Ice Sheet. We focus on Larenas Bay (48°S; 1.26 km2), which is connected to Baker Channel through a shallow (7.4 m) and narrow (ca. 150 m across) inlet, and hence has the potential to record periods of basin isolation and marine ingression. The paleoenvironmental evolution of the bay was investigated through a sedimentological analysis of a 9.2 m long radiocarbon-dated sediment core covering the last 16.8 kyr. Salinity indicators, including diatom paleoecology, alkenone concentrations and CaCO3 content, were used to reconstruct the bay's connectivity to the fjord. Results indicate that Larenas Bay was a marine environment before 16.5 cal kyr BP and after 9.1 cal kyr BP, but that it was disconnected from Baker Channel in-between. We infer that the postglacial rebound started before 16.5 cal kyr BP and that it outpaced global sea-level rise until slightly before 9.1 cal kyr BP. During the Late Glacial and early Holocene, the center of the former Patagonian Ice Sheet experienced an absolute uplift of ca. 96 m, at an average rate of 1.3 cm/yr. During the remainder of the Holocene, glacial isostatic adjustment continued (ca. 20 m), but at a slower average pace of 0.2 cm/yr. Comparisons between multi-millennial variations in the salinity indicators and existing records of global sea-level rise suggest that the glacial isostatic adjustment rate also fluctuated within these time intervals, likely in response to glacier dynamics. More specifically, most of the glacial isostatic adjustment registered between 16.5 – 9.1 cal kyr BP seems to have occurred before meltwater pulse 1A (14.5 – 14.0 kyr BP). Likewise, it appears that the highest glacial isostatic rebound rates of the last 9.1 kyr occurred during the late Holocene, most likely in response to glacier recession from their Neoglacial maxima. This implies a relatively rapid response of the local solid earth to ice unloading, which agrees with independent modelling studies investigating contemporary uplift. We conclude that the center of the former Patagonian Ice Sheet experienced a glacial isostatic adjustment of ca. 116 m over the last 16.5 kyr, and that >80% occurred during the Late Glacial and early Holocene.

Deglacial Land-Ocean Linkages at the Alaskan Continental Margin in the Bering Sea

A marine sediment record from the central Bering Sea, spanning the last 20 thousand years (ka), was studied to unravel the depositional history with regard to terrigenous sediment supply and biogenic sedimentation. Methodic approaches comprised the inference of accumulation rates of siliciclastic and biogenic components, grain-size analysis, and (clay) mineralogy, as well as paleoclimatic modelling. Changes in the depositional history provides insight into land-ocean linkages of paleoenvironmental changes. During the finale of the Last Glacial Maximum, the depositional environment was characterized by hemipelagic background sedimentation. A marked change in the terrigenous sediment provenance during the late Heinrich 1 Stadial (15.7–14.5 ka), indicated by increases in kaolinite and a high glaciofluvial influx of clay, gives evidence of the deglaciation of the Brooks Range in the hinterland of Alaska. This meltwater pulse also stimulated the postglacial onset of biological productivity. Glacial melt implies regional climate warming during a time of widespread cooling on the northern hemisphere. Our simulation experiment with a coupled climate model suggests atmospheric teleconnections to the North Atlantic, with impacts on the dynamics of the Aleutian Low system that gave rise to warmer winters and an early onset of spring during that time. The late deglacial period between 14.5 and 11.0 ka was characterized by enhanced fluvial runoff and biological productivity in the course of climate amelioration, sea-level rise, seasonal sea-ice retreat, and permafrost thaw in the hinterland. The latter processes temporarily stalled during the Younger Dryas stadial (12.9-11.7 ka) and commenced again during the Preboreal (earliest Holocene), after 11.7 ka. High river runoff might have fertilized the Bering Sea and contributed to enhanced upper ocean stratification. Since 11.0 ka, advanced transgression has shifted the coast line and fluvial influence of the Yukon River away from the study site. The opening of the Bering Strait strengthened contour currents along the continental slope, leaving behind winnowed sand-rich sediments through the early to mid-Holocene, with non-deposition occurring since about 6.0 ka.

Penultimate deglaciation Asian monsoon response to North Atlantic circulation collapse

During glacial terminations, massive iceberg discharges and meltwater pulses in the North Atlantic triggered a shutdown of the Atlantic Meridional Overturning Circulation (AMOC). Speleothem calcium carbonate oxygen isotope records (δ18OCc) indicate that the collapse of the AMOC caused dramatic changes in the distribution and variability of the East Asian and Indian monsoon rainfall. However, the mechanisms linking changes in the intensity of the AMOC and Asian monsoon δ18OCc are not fully understood. Part of the challenge arises from the fact that speleothem δ18OCc depends on not only the δ18O of precipitation but also temperature and kinetic isotope effects. Here we quantitatively deconvolve these parameters affecting δ18OCc by applying three geochemical techniques in speleothems covering the penultimate glacial termination. Our data suggest that the weakening of the AMOC during meltwater pulse 2A caused substantial cooling in East Asia and a shortening of the summer monsoon season, whereas the collapse of the AMOC during meltwater pulse 2B (133,000 years ago) also caused a dramatic decrease in the intensity of the Indian summer monsoon. These results reveal that the different modes of the AMOC produced distinct impacts on the monsoon system.

Decadal-scale onset and termination of Antarctic ice-mass loss during the last deglaciation

Emerging ice-sheet modeling suggests once initiated, retreat of the Antarctic Ice Sheet (AIS) can continue for centuries. Unfortunately, the short observational record cannot resolve the tipping points, rate of change, and timescale of responses. Iceberg-rafted debris data from Iceberg Alley identify eight retreat phases after the Last Glacial Maximum that each destabilized the AIS within a decade, contributing to global sea-level rise for centuries to a millennium, which subsequently re-stabilized equally rapidly. This dynamic response of the AIS is supported by (i) a West Antarctic blue ice record of ice-elevation drawdown >600 m during three such retreat events related to globally recognized deglacial meltwater pulses, (ii) step-wise retreat up to 400 km across the Ross Sea shelf, (iii) independent ice sheet modeling, and (iv) tipping point analysis. Our findings are consistent with a growing body of evidence suggesting the recent acceleration of AIS mass loss may mark the beginning of a prolonged period of ice sheet retreat and substantial global sea level rise.

Hydroclimate response in Texas and Gulf of Mexico to rapid warming during the last deglacial: High-resolution speleothem proxy and monitoring evidence

Cave calcite (i.e., speleothem) proxy records that span abrupt warming events, such as those of the last deglacial, may prologue regional responses of hydroclimate to 21st century warming. Proxy reliability is advanced by an understanding of the modern hydrologic system and controls on calcite growth. Here we integrate monitoring of Cave CWN in central Texas with high-resolution imaging and δ18O analysis of Stalagmite CWN-4 to understand how regional hydroclimate changed during a rapid deglacial warming event - the transition from Heinrich Stadial 1 (HS1) to the Bølling-Allerød (BA) warm period. This event is recorded in CWN-4 as a rapid negative δ18O excursion of 1.7‰ over ∼210 years (HS1-BA δ18O excursion = ‘HBO’), followed by a 6-fold growth rate increase.Monitoring demonstrates that modern calcite grows during cool seasons, growth rates follow drip rate, and drip water and calcite δ18O reflect the δ18O of recharge. Because drip waters have interannual epikarst residence times, their δ18O values are well-mixed averages of their Gulf of Mexico (GoM) moisture source δ18O values, which rapidly decreased during Laurentide Ice Sheet (LIS) melting events. Confocal laser fluorescence microscopy imaging indicates that growth lamina are seasonal within U-series age constraints. The timing of the HBO (14.75 ± 0.08 to 14.54 ± 0.08 ka), corresponds with 1) the onset of Bølling warming in Greenland (14.64 ka), 2) AMOC reinvigoration (∼14.61–14.25 ka), 3) Meltwater Pulse 1a (14.65–14.31 ka), and 4) Meltwater Flood 3 (MWF-3) - the most intense episode of LIS meltwater discharge into the GoM (14.97–14.46 ka). The rapid increase in growth rates following the HBO indicates an abrupt transition to a wetter hydroclimate in central Texas following these regional and global events.The HBO δ18O excursion occurred monotonically, consistent with meltwater discharge being continuous. This excursion also corresponds with climate change recognized across the HS1-BA transition in marine and terrestrial records influenced by changes in North Atlantic climate, AMOC, and ITCZ migrations. We infer that central Texas deglacial hydroclimate largely followed changes in AMOC, linked via meltwater input to the GoM. Periods of fast (slow) speleothem growth coincide with warmer (cooler) North Atlantic episodes, with the exception of the HS1 to BA transition. For this transition, slow growth rates during MWF-3 may reflect protracted cooling of GoM SST and corresponding atmospheric moisture reduction. The common timing and duration of the HBO with shifts toward wetter conditions across the HS1-BA transition in many northern Hemisphere climate records supports common ocean-atmosphere teleconnections related to AMOC intensity.

Freshwater routing in eddy-permitting simulations of the last deglacial: the impact of realistic freshwater discharge

Abstract. Freshwater, in the form of glacial runoff, is hypothesized to play a critical role in centennial- to millennial-scale climate variability, such as the Younger Dryas and Dansgaard–Oeschger events, but this relationship is not straightforward. Large-scale glacial runoff events, such as Meltwater Pulse 1a (MWP1a), are not always temporally proximal to subsequent large-scale cooling. Moreover, the typical design of hosing experiments that support this relationship tends to artificially amplify the climate response. This study explores the impact that limitations in the representation of runoff in conventional “hosing” simulations has on our understanding of this relationship by examining where coastally released freshwater is transported when it reaches the ocean. We particularly focus on the impact of (1) the injection of freshwater directly over sites of deep-water formation (DWF) rather than at runoff locations (i.e. hosing), (2) excessive freshwater injection volumes (often by a factor of 5), and (3) the use of present-day (rather than palaeo) ocean gateways. We track the routing of glaciologically constrained freshwater volumes from four different inferred injection locations in a suite of eddy-permitting glacial ocean simulations using the Massachusetts Institute of Technology General Circulation Model (MITgcm) under both open and closed Bering Strait conditions. Restricting freshwater forcing values to realistic ranges results in less spreading of freshwater across the North Atlantic and indicates that the freshwater anomalies over DWF sites depend strongly on the geographical location of meltwater input. In particular, freshwater released into the Gulf of Mexico generates a very weak freshwater signal over DWF regions as a result of entrainment by the turbulent Gulf Stream. In contrast, freshwater released into the Arctic with an open Bering Strait or from the Eurasian ice sheet is found to generate the largest salinity anomalies over DWF regions in the North Atlantic and GIN (Greenland–Iceland–Norwegian) seas region respectively. Experiments show that when the Bering Strait is open, the Mackenzie River source exhibits more than twice as much freshening of the North Atlantic deep-water formation regions as when the Bering Strait is closed. Our results illustrate that applying freshwater hosing directly into the North Atlantic with even “realistic” freshwater amounts still overestimates the amount of terrestrial runoff reaching DWF regions. Given the simulated salinity anomaly distributions and the lack of reconstructed impact on deep-water formation during the Bølling–Allerød, our results support that the majority of the North American contribution to MWP1a was not routed through the Mackenzie River.

Sedimentation and organic content in the mires and other sites of sediment accumulation in the Sydney region, eastern Australia, in the period after the Last Glacial Maximum

This research reports on a synthesis of radiocarbon (14C) dates and the organic content from a variety of depositional sites in a relatively small region in humid eastern Australia centred on Sydney. We focused on basal dates, changes in accumulation rates and the organic content of these sediment records in the period post-dating 21 thousand years ago to make inferences about past environmental conditions, or to infer the timing of past environmental change. We found that low rates of sediment and organic accumulation at the Last Glacial Maximum continued well into the late Pleistocene. The average rate of sediment deposition and organic content increased from about 14.5 calibrated kiloanni before present (cal ka BP), perhaps coincident with Meltwater Pulse 1A, but this was then checked by dry conditions during the Antarctic Climatic Reversal chronozone. An abrupt increase in the average rate of sediment accumulation coincided with the Pleistocene-Holocene boundary (at 11.7 cal ka BP), and this continued until about 9.2 cal ka BP, largely coinciding with the Antarctic thermal maximum, and probably reflecting reduced moisture availability in the Sydney region. Indication of depressed temperature and a more positive moisture availability, from about 9 to 7.6 cal ka BP, witnessed the fastest increase in the number of sites than at any other time. An overall drier early Holocene abruptly altered at about 7.5 cal ka BP and a mid-Holocene ‘optimum’ in moisture availability extended to 6.2 cal ka BP. The compilation of organic contents of the mid-late Holocene suggests environmental variability, but El Niño frequency or strength is not implicated as a driving force. Increased organic productivity, especially after 2.8 cal ka BP, corresponded to a peak in summer insolation, and associated seasonality of insolation and this is likely to have resulted in increased, or more consistent, summer rainfall and easterly air flow into the Sydney region. In contrast to previous generalisations about changes in south-eastern Australia during the Holocene, which have often relied on palaeoenvironmental records from further south, the southward migration of the intertropical convergence zone during the Holocene had vastly different synoptic outcomes in the Sydney region, with the early Holocene drier as continental air dominated, and the late Holocene wetter, as the synoptic conditions saw more tropical and easterly airflow. This work demonstrates the danger in making generalisations across regions and highlights the utility of the compilation of data at a homoclimatic, regional scale.

Morphological evidences of eustatic events in the last 14,000 years in a far-field site, East-Southeast Brazilian continental shelf

The post-LGM transgression was marked by eustatic variability events, comprising periods of increased sea-level rise, triggered by meltwater pulses (MWP) on polar regions, preceded by periods of stillstand. Several studies highlight the influence of such events on the geomorphology of far-field continental margins worldwide by mapping drowned relict features. However, there is still a lack of information in the South Atlantic and a major discussion about the global distribution of MWP-1B. The present work aims to investigate morphological evidences in the east-southeast Brazilian continental shelf (Abrolhos and Espírito Santo shelf), by examining shelf morphology using regional bathymetric data and high-resolution seismic sections. The analytical strategy involved comparing depth signatures of interpreted relict features with a set of sea-level curves. The results suggested the predominance of depositional process in shoreline sectors under closer influence of fluvial discharges during the Younger Dryas (YD). During YD stillstand or low sea-level rise rates (~60 m deep), it was possible to recognize the occurrence of a lagoonal system with shallow water/shoreface features (sandbanks, sand ribbons, or other bedforms) associated with tidal channels in the Abrolhos Shelf Depression and erosive features such as coastal cliffs. In the Espírito Santo outer shelf, the YD coastal features were characterized by barrier island systems associated with incised valley mouths and possible fringing reefs, in a context of a marginal shallow marine environment. All these features tend to be preserved in the modern shelf morphology due to a fast increase in the accommodation space that defines the end of the YD due to the onset of MWP-1B. In general, morphological interpretation revealed that MWP-1B caused the overstepping or drowning of those coastal features. Erosive features, defining the YD period, were also identified far from the influence of fluvial input. In these areas, coastal erosive processes prevailed and are morphologically expressed as abrasion terraces developed since the YD onset ~70 m and features like a paleo-sea cliff. During MWP-1B, the presence of shallow water reefs in the Abrolhos shelf is hypothesized, experiencing backstepping until their complete submersion. Despite the lack of additional data on composition, age and 3D geometry of the features, the morphological approach proposed here strongly supports the concept that eustatic events, triggered mainly in the Northern Hemisphere, influenced the shelf geomorphology and sedimentary processes evolution in the South Atlantic, during the last deglaciation.