Sea-level Reconstructions Research Articles

Holocene history of Fiskesø, Prinsesse Ingeborg Halvø, eastern North Greenland

Sediment cores up to 92 cm in length were recovered with gravity and Russian peat corers from Fiskesø, Prinsesse Ingeborg Halvø, eastern North Greenland, during the summer 2014. The correlated sediment succession consists of clastic sediments that are interspersed in the upper part with layers of organic material and likely record the environmental history of the lake since deglaciation. The paucity of macrofossil remains hampered radiocarbon dating of the sediments from Fiskesø. According to published data, the deglaciation of the region took place c. 10 cal. ka BP. Relative sea-level reconstructions from the region suggest that the Fiskesø basin, which is today located at 33 m above sea level, was characterised by marine conditions until c. 8.1 cal. ka BP. Marine fossils in the lower part of the sediment succession support the prevalence of marine conditions. A reliable radiocarbon age from 5 cm above the isolation horizon in Fiskesø sediments indicates an age of 6.1 cal. ka BP and supports the isolation of the basin prior to this time. Cooling is indicated in the upper part of the sediment succession and is reported to have taken place in the region stepwise shortly after 6.1 cal. ka BP and at c. 4.5–4.0 cal. ka BP. Despite the poor chronology the data from Fiskesø support existing terrestrial and marine reconstructions from the region.

Comprehensive comparative analysis of reconstructed sea level datasets in the China Seas: insights from tide gauge and satellite altimetry

At present, there are many reconstructed datasets at the global scale. To test the applicability of these datasets in the China seas, the study comprehensively analyzes the reliability and accuracy of reconstructed sea level datasets in capturing nuanced temporal patterns of sea level changes in the China Seas. This study applied analysis methods or indicators such as time series, Taylor plots, correlation coefficients, growth rates, and standard deviations. Ocean Data Assimilations (ODAs) outperform Tide Gauge Reconstructions (TGRs) in terms of correlation with measured data in the nearshore, while TGRs exhibit superior capability in capturing oceanic sea level variability. Although the ODAs and TGRs both suffer from the underestimation of sea level variability in China as well as in neighboring seas, the TGRs perform better than the former. ODAs show inconsistency in reflecting the rate of sea level rise, but they, particularly the China Ocean Reanalysis (CORA), demonstrate a better correlation with satellite altimetry datasets. Meanwhile, both of them can reflect the Pacific Decadal Oscillation (PDO) well. TGRs, relying on oceanic tide gauge stations, suffer from poor correlation with tide gauge stations due to limited coverage. Reconstruction discrepancies are attributed to methodological differences and data assimilation techniques. Future studies should explore alternative variables like sea surface temperature and so on to enhance sea-level reconstruction, especially in regions with sparse tide gauge coverage.

Ice Sheet‐Albedo Feedback Estimated From Most Recent Deglaciation

AbstractIce sheet feedbacks are underrepresented in model assessments of climate sensitivity and their magnitudes are still poorly constrained. We combine a recently published record of Earth's Energy Imbalance (EEI) with existing reconstructions of temperature, atmospheric composition, and sea level to estimate both the magnitude and timescale of the ice sheet‐albedo feedback since the Last Glacial Maximum. This facilitates the first opportunity to quantify this feedback over the most recent deglaciation using a proxy data‐driven approach. We find the ice sheet‐albedo feedback to be amplifying, increasing the total climate feedback parameter by 42% and reaching an equilibrium magnitude of 0.55 Wm−2K−1, with a 66% confidence interval of 0.45–0.63 Wm−2K−1. The timescale to equilibrium is estimated as 3.6 ka (66% confidence: 1.9–5.5 ka). These results provide new evidence for the timescale and magnitude of the amplifying ice sheet‐albedo feedback that will drive anthropogenic warming for millennia to come.

Probabilistic reconstruction of sea-level changes and their causes since 1900

Abstract. Coastal communities around the world are increasingly exposed to extreme events that have been exacerbated by rising sea levels. Sustainable adaptation strategies to cope with the associated threats require a comprehensive understanding of past and possible future changes. Yet, many coastlines lack accurate long-term sea-level observations. Here, we introduce a novel probabilistic near-global reconstruction of relative sea-level changes and their causes over the period from 1900 to 2021. The reconstruction is based on tide gauge records and incorporates prior knowledge about physical processes from ancillary observations and geophysical model outputs, allowing us, for the first time, to resolve individual processes and their uncertainties. We demonstrate good agreement between the reconstruction and satellite altimetry and tide gauges (if local vertical land motion is considered). Validation against steric height estimates based on independent temperature and salinity observations over their overlapping periods shows moderate to good agreement in terms of variability, though with larger reconstructed trends in three out of six regions. The linear long-term trend in the resulting global-mean sea-level (GMSL) record is 1.5 ± 0.19 mm yr−1 since 1900, a value consistent with central estimates from the 6th Assessment Report of the Intergovernmental Panel on Climate Change. Multidecadal trends in GMSL have varied; for instance, there were enhanced rates in the 1930s and near-zero rates in the 1960s, although a persistent acceleration (0.08 ± 0.04 mm yr−2) has occurred since then. As a result, most recent rates have exceeded 4 mm yr−1 since 2019. The largest regional rates (>10 mm yr−1) over the same period have been detected in coastal areas near western boundary currents and the larger tropical Indo-Pacific region. Barystatic mass changes due to ice-melt and terrestrial-water-storage variations have dominated the sea-level acceleration at global scales, but sterodynamic processes are the most crucial factor locally, particularly at low latitudes and away from major melt sources. These results demonstrate that the new reconstruction provides valuable insights into historical sea-level change and its contributing causes, complementing observational records in areas where they are sparse or absent. The Kalman smoother sea-level reconstruction dataset can be accessed at https://doi.org/10.5281/zenodo.10621070 (Dangendorf, 2024).

Pliocene shorelines and the epeirogenic motion of continental margins: a target dataset for dynamic topography models

Abstract. Global mean sea level during the mid-Pliocene epoch (∼3 Ma), when CO2 and temperatures were above present levels, was notably higher than today due to reduced global ice sheet coverage. Nevertheless, the extent to which ice sheets responded to Pliocene warmth remains in question owing to high levels of uncertainty in proxy-based sea level reconstructions as well as solid Earth dynamic models that have been used to evaluate a limited number of data constraints. Here, we present a global dataset of 10 wave-cut scarps that formed by successive Pliocene sea level oscillations and which are observed today at elevations ranging from ∼6 to 109 m above sea level. The present-day elevations of these features have been identified using a combination of high-resolution digital elevation models and field mapping. Using the MATLAB interface TerraceM, we extrapolate the cliff and platform surfaces to determine the elevation of the scarp toe, which in most settings is buried under meters of talus. We correct the scarp-toe elevations for glacial isostatic adjustment and find that this process alone cannot explain observed differences in Pliocene paleo-shoreline elevations around the globe. We next determine the signal associated with mantle dynamic topography by back-advecting the present-day three-dimensional buoyancy structure of the mantle and calculating the difference in radial surface stresses over the last 3 Myr using the convection code ASPECT. We include a wide range of present-day mantle structures (buoyancy and viscosity) constrained by seismic tomography models, geodynamic observations, and rock mechanics laboratory experiments. Finally, we identify preferred dynamic topography change predictions based on their agreement with scarp elevations and use our most confident result to estimate a Pliocene global mean sea level based on one scarp from De Hoop, South Africa. This inference (11.6 ± 5.2 m) is a downward revision and may imply that ice sheets were relatively resistant to warm Pliocene climate conditions. We also conclude, however, that more targeted model development is needed to more reliably infer mid-Pliocene global mean sea level based on all scarps mapped in this study.

Re-evaluating Marine Isotope Stage 5a paleo-sea-level trends from across the Florida Keys reef tract

Unraveling how Global Mean Sea Level (GMSL) fluctuated during past warm periods can improve our understanding of linkages between sea-level fluctuations, orbital forcing, and ice-sheet dynamics. Current estimates of GMSL for Marine Isotope Stages (MIS) 5a and 5c — two warm intervals following the relatively well-documented MIS 5e — contain meters of uncertainty and fewer data due to several challenges. These challenges include concealment of datable in-situ coral facies by MIS 1 deposits and inaccessibility due to submergence by modern sea level. We present a comprehensive dataset based on U–Th dating and stratigraphic correlation of 23 cores totaling over 170 m of recovered coral-reef deposits across the tectonically stable Florida Keys Reef Tract (FKRT). Following detailed facies descriptions, 34 in-situ, minimally altered aragonitic coral samples (≤2.7% calcite) below the Holocene-Pleistocene boundary were targeted for U–Th geochronology. Fourteen closed-system coral U–Th ages from MIS 5a include the commonly used sea-level indicator Acropora palmata, but also the massive coral taxa Pseudodiploria strigosa, Siderastrea siderea, Orbicella spp., and Porites astreoides. Dating yielded ages in the range of 88–81 ka (average 2σ uncertainty of less than 200 years). These ages suggest MIS 5a reef initiation at ∼88 ka BP, a peak near 83 ka with minimum elevations between −6.0 ± 0.5 and −5.6 ± 0.5 m MSL (2σ uncertainty and subsidence-corrected), and reef termination and sea-level fall by ∼81 ka BP. Notably, the range of peak MIS 5a relative sea-level estimates of −6.5 to −5.1 m MSL are more than 2 m shallower (higher) than previous estimates of −11 to −9 m. Our higher resolution regional sea-level reconstruction across four subregions of the Florida Keys reef tract aligns with changes in July insolation at 65° N: a trend that most other records, such as deep-sea sediments, do not have the accuracy and precision to resolve. Three massive coral samples from MIS 5c, consisting of Pseudodiploria clivosa, and Orbicella spp., yielded ages in the range of 104 to 99 ka (average 2σ uncertainty less than 200 years); however, because only one sample met the closed-system criteria, our ability to estimate MIS 5c sea level is relatively limited. More empirical estimates of sea-level from the MIS 5a and MIS 5c intervals based on numerical dating of reliable local sea-level constraints are critical for GMSL calculations and relating changes in sea-level amplitude and timing to global ice volume modeling and glacio-isostatic effects, all of which can improve predictions of future sea-level changes in coastal regions.

Effects of using the consistent boundary flux method on dynamic topography estimates

SUMMARY Dynamic topography is defined as the deflection of Earth's surface due to the convecting mantle. ASPECT (Advanced Solver for Planetary Evolution, Convection, and Tectonics) is a continually evolving, finite element code that uses modern numerical methods to investigate problems in mantle convection. With ASPECT version 2.0.0 a consistent boundary flux (CBF) algorithm, used to calculate radial stresses at the model boundaries, was implemented into the released version of ASPECT. It has been shown that the CBF algorithm improves the accuracy of dynamic topography calculations by approximately one order of magnitude. We aim to evaluate the influence of the CBF algorithm and explore the geophysical implications of these improved estimates of dynamic topography changes along the East Coast of the United States. We constrain our initial temperature conditions using the tomography models SAVANI, S40RTS and TX2008, and combine them with a corresponding radial viscosity profile (2 for TX2008) and two different boundary conditions for a total of eight experiments. We perform simulations with and without the CBF method, which takes place during post-processing and does not affect the velocity solution. Our dynamic topography calculations are spatially consistent in both approaches, but generally indicate an increase in magnitude using the CBF method (on average ∼15 and ∼76 per cent absolute change in present-day instantaneous and rate of change of dynamic topography, respectively). This enhanced accuracy in dynamic topography calculations can be used to better evaluate the effects of mantle convection on surface processes including vertical land motions, sea level changes, and sedimentation and erosion. We explore results along the US East Coast, where a Pliocene shoreline has been deformed by dynamic topography change. An increased accuracy in estimates of dynamic topography can improve Pleistocene and Pliocene sea level reconstructions, which allow for a better understanding of past sea level changes and ice sheet stability.

Impact of coseismic uplifting on relative sea level change in the Southern Adriatic during the past 4.500 years – New evidence from Dubrovnik epicentral area based on analysis of algal rims and tidal notches

The Southern Adriatic is recognized as one of the most seismically active regions of the Central Mediterranean. It hosted the strongest historically known and instrumentally recorded earthquakes in Croatia and Montenegro. We conducted a detailed study of algal rims and tidal notches along the Dubrovnik–Konavle coastal area, aimed to contribute in reconstruction of relative sea level (RSL) change, palaeoearthquakes and neotectonics during the last 4500 years. The RSL reconstruction based on high-resolution geochronology showed that the spatial variability of RSL change is largely controlled by local tectonics, leading generally to an alternation of periods of RSL rise and rapid falls caused by coseismic uplift events. During the interseismic periods RSL rose at rates from 1.03 to 0.7 mm/yr before AD, followed by a slowdown to ~0.43 mm/yr, and reaching ~1.3 mm/yr after 1850s. The only exception relates to a period of slow RSL drop probably related to the 3.2 ka cold event. Rapid RSL falls caused by coseismic uplift indicate the occurrence of several large-scale events during the 4th–6th centuries AD, the 800–1360 cal AD period, and in 1395, with hints of earlier events between 1395 cal BC and 565 cal AD. These instances correspond to either unknown or insufficiently documented earthquakes. Conversely, the earthquakes in 1520 and 1667 have been clearly identified. The most compelling evidence of coseismic uplift is associated with the 1667 Dubrovnik earthquake, with an estimated uplift of 40–60 ± 15 cm documented along 40 km of coastline. Analyses of multiple sites with the same high-resolution marker show that Holocene tectonic deformations vary along different fault segments with minimum uplift rates ranging from ~0.3 to ~0.9 mm/yr over ~1.5 ka. The main seismogenic source responsible for coseismic uplifts is attributed to the Dalmatian unit basal thrust and its NE-dipping splays at the Adria-Dinarides boundary.

Influence of Foraminifera Count Size and Rare Species on Transfer Function Results Used in Sea-Level Reconstructions

Influence of Foraminifera Count Size and Rare Species on Transfer Function Results Used in Sea-Level Reconstructions

Sea-level and monsoonal control on the Maldives carbonate platform (Indian Ocean) over the last 1.3 million years

Abstract. The Maldives Archipelago (Indian Ocean), composed of two rows of atolls that enclose an inner sea, offers an excellent study site to explore the forcings of carbonate production at platforms. Glacial–interglacial sea-level changes have been claimed to be the main factor controlling the carbonate platform factories; however, climatic factors may also have an impact. In this work we used geochemical compositional records, obtained by X-ray fluorescence (XRF) core-scanning from the International Ocean Discovery Program (IODP) Site U1467 in the Maldives Inner Sea, to analyze the orbitally driven fluctuations on the carbonate production and export from the neritic environment into the Maldives Inner Sea over the last 1.3 million years. High Sr aragonite-rich carbonates (HSAC) from neritic settings were deposited in the Maldives Inner Sea during sea-level highstand intervals, increasing the Sr/Ca values. In contrast, low Sr/Ca values are observed coincident with sea-level lowstand periods, suggesting that large areas of the atolls were exposed or unable to grow, and therefore, there was a demise in the carbonate production and sediment export to the Maldives Inner Sea. However, comparison of the Sr/Ca values and the sea-level reconstructions for different interglacial periods before and after the mid-Brunhes event (MBE, ∼ 430 ka) indicates that sea level is not the only factor controlling the production of HSAC during sea-level highstands. The study of monsoon and primary productivity proxies (Fe-normalized, Fe/K, and Br-normalized records) from the same site suggests that the intensity of the summer monsoon and the Indian Ocean dipole probably modulated the carbonate production at the atolls. Moreover, Marine Isotope Stage 11 stands out as a period with high sea level and extraordinary carbonate production in the Maldives platform. This outstanding carbonate production in the Maldives atolls (and in other low-latitude carbonate platforms) probably contributed to the mid-Brunhes dissolution event through a strong shelf-to-basin fractionation of carbonate deposition.

A Review of the Contribution of Satellite Altimetry and Tide Gauge Data to Evaluate Sea Level Trends in the Adriatic Sea within a Mediterranean and Global Context

The relatively new sea level satellite altimetry and secular coastal tide gauge data made the reconstruction of sea levels on regional and global scales possible about one century back. Due to better estimations of the Earth’s crustal, glacial, tectonic, and other possible motion biases in tide gauge data, some additional improvements can be expected in sea level reconstructions, analysis, and predictions. A more detailed review of published sea level-related results was conducted for the Eastern Adriatic coast, including the operation of the tide gauge network and data processing, crustal movement estimations, and the establishment of a new reference height system in Croatia, based on five tide gauge sea level data. It was shown that sea level variation and trend-related indicators are spatially homogeneous, especially on a sub-Adriatic scale. The regional Adriatic Sea mean sea level rise rate of +2.6 mm/year for the satellite altimetry era (1993–2019) is less than the global mean sea level (GMSL) rise rate of +3.3 mm/year for the period of 1993–2022. Several empirical methods for GMSL projections and expected IPCC (Intergovernmental Panel on Climate Change) assessments until the end of the 21st century are considered.

Glacial seismic geomorphology and Plio-Pleistocene ice sheet history offshore NW Europe

Abstract Plio-Pleistocene records of ice-rafted detritus suggest NW European ice sheets regularly reached coastlines. However, these records provide limited insight on the frequency, extent and dynamics of ice sheets delivering the detritus. Three-dimensional reflection seismic data of the NW European glaciated margin have previously documented buried landforms that inform us on these uncertainties. This paper reviews and combines these existing records with new seismic geomorphological observations to catalogue landform occurrence along the European glaciated margin and considers how they relate to ice sheet history. The compilation shows Early Pleistocene ice sheets regularly advanced across the continental shelves. Early Pleistocene sea-level reconstructions demonstrate lower magnitude fluctuations compared to the Middle–Late Pleistocene, and more extensive/frequent Early Pleistocene glaciation provides a possible mismatch with sea-level reconstructions. This evidence is discussed with global records of glaciation to consider possible impacts on our wider understanding of Plio-Pleistocene climate changes, in particular how well Early Pleistocene sea-level records capture ice sheet volume changes. Resolving such issues relies on how well landforms are dated, whether they can be correlated with other proxy datasets, and how accurately these proxies reconstruct the magnitudes of past climatic changes. Many questions about Pleistocene glaciation in Europe and elsewhere remain.

Imprint of relative sea level histories on Last Interglacial coral preservation

SUMMARY Fossil corals are commonly used to reconstruct Last Interglacial (∼125 ka, LIG) sea level. Sea level reconstructions assume the water depth at which the coral lived, called the ‘relative water depth’. However, relative water depth varies in time and space due to coral reef growth in response to relative sea level (RSL) changes. RSL changes can also erode coral reefs, exposing older reef surfaces with different relative water depths. We use a simplified numerical model of coral evolution to investigate how sea level history systematically influences the preservation of corals in the Bahamas and western Australia, regions which house >100 LIG coral fossils. We construct global ice histories spanning the uncertainty of LIG global mean sea level (GMSL) and predict RSL with a glacial isostatic adjustment model. We then simulate coral evolution since 132 ka. We show that preserved elevations and relative water depths of modelled LIG corals are sensitive to the magnitude, timing and number of GMSL highstand(s). In our simulations, the influence of coral growth and erosion (i.e. the ‘growth effect’) can have an impact on RSL reconstructions that is comparable to glacial isostatic adjustment. Thus, without explicitly accounting for the growth effect, additional uncertainty is introduced into sea level reconstructions. Our results suggest the growth effect is most pronounced in western Australia due to Holocene erosion, but also plays a role in the Bahamas, where LIG RSL rose rapidly due to the collapsing peripheral bulge associated with Laurentide Ice Sheet retreat. Despite the coral model's simplicity, our study highlights the utility of process-based RSL reconstructions.

A ∼200-year relative sea-level reconstruction from the Wellington region (New Zealand) reveals insights into vertical land movement trends

The sea-level rise threat to New Zealand's coastal cities is regionally exacerbated due to spatially varying vertical land movement (VLM). At Wellington, the capital city, situated adjacent to a major active plate boundary, strong regional spatial and temporal variability of VLM is indicated by the relatively short (∼25 year-long) continuous Global Navigation Satellite Systems (GNSS) network, but until now longer records of VLM have been lacking. Here, a ∼ 200-year-long relative sea-level reconstruction is presented from Pāuatahanui salt marsh in the northern Wellington region. The foraminifera-based relative sea-level reconstruction indicates that ∼1 ± 0.45 m of sudden uplift occurred during the 1855CE Mw 8.2 Wairarapa earthquake. Following this, Pāuatahanui has experienced a mean rate of relative sea-level rise (1855CE to present) of 1.5 ± 0.6 mm/yr, or 2.4 ± 0.8 mm/yr since the start of the twentieth century, consistent with ongoing subsidence in concert with climate-driven sea-level rise. Further acceleration to >3 mm/yr since the 1990s (with 4 mm/yr also possible if the full 95% confidence range is taken into consideration) is consistent with the globally documented acceleration in sea-level rise, although low model precision hampers confidence in this interpretation. This record is the first of its kind from a tectonically complex setting in New Zealand, shedding light on the effects of the historically significant 1855 earthquake, and fills a gap between millennial-scale and contemporary records of VLM with important implications for future sea-level projections in the region.

Sudden disappearance of yew (Taxus baccata) woodlands from eastern England coincides with a possible climate event around 4.2 ka ago

Tree-ring chronologies form the backbone of high-resolution palaeoclimatology. However, their number declines drastically prior to medieval times, and only a few such records worldwide extend back to the mid-Holocene. Here, we present a collection of more than 400 subfossil yew (Taxus baccata L.) trees excavated from near sea-level peat-rich sediments in the Fenland region of eastern England. The well-preserved yew trunks are between two and eight metres long, often exhibit adventitious root layers, and contain up to 400 rings of highly irregular growth. Combined dendrochronological and radiocarbon dating resulted in two tree-ring width chronologies that comprise 36 and 32 trees, span 413 and 418 years, exhibit mean inter-series correlations (Rbar) of 0.50 and 0.51, and were provisionally dated to 5225−4813 (±4) and 4612−4195 (±6) years cal BP. Together with a total of 63 radiocarbon dates (14C), our subfossil tree-ring evidence suggests that yew establishment (or onset of preservation) began ∼5250 years cal BP and resulted in extensive climax forests between 5200 and 4200 years cal BP. A first stage of yew decline ∼4800−4600 years cal BP was possibly caused by oxygen deprivation from soil wetting, whereas yew disappearance ∼4200 years cal BP is suggested to have been triggered by marine inundation as a consequence of rapid sea-level rise. Both phases of yew decline in eastern England coincide with marked reductions in subfossil oak and pine from peatbogs in Ireland, Germany and the Netherlands. Our results emphasise the potential to develop a dendrochronological network in coastal England for better understanding of larger-scale climate and environmental changes during the mid-Holocene, including the still debated 4.2 ka climate event. Moreover, we expect our subfossil yew chronologies to facilitate the dating of local archaeological remains, refine sea-level reconstructions around the British Isles, and contribute to the international radiocarbon calibration curve IntCal.

Late Pleistocene Evolution of Tides and Tidal Dissipation

AbstractStudies of the Last Glacial Maximum (LGM; 26.5–19 ka) tides showed strong enhancements in open ocean tidal amplitudes and dissipation rates; however, changes prior to the LGM remain largely unexplored. Using two different ice sheet and sea level reconstructions, we explicitly simulate the evolution of the leading semi‐diurnal and diurnal tidal constituents (M2, S2, K1, and O1) over the last glacial cycle with a global tide model. Both sets of simulations show that global changes, dominated by the Atlantic, take place for the semi‐diurnal constituents, while changes for the diurnal constituents are mainly regional. Irrespective of the reconstruction, open ocean dissipation peaks during the sea level lowstands of MIS 2 (∼20 ka) and MIS 4 (∼60 ka), although dissipation values prior to MIS 2 are sensitive to differences in reconstructed ice sheet extent. Using the statistically significant relationship between global mean sea level and dissipation, we apply regression analysis to infer open ocean and shelf dissipation, respectively, over the last four glacial cycles back to 430 ka. Our analysis shows that open ocean tidal energy was probably increased for most of this period, peaking during glacial maxima, and returning to near‐present‐day values during interglacials. Due to tidal resonance during glacial phases, small changes in bathymetry could have caused large changes in tidal amplitudes and dissipation, emphasizing the need for accurate ice margin reconstructions. During glacial phases, once global mean sea level decreased by more than ∼100 m, the amount of open ocean tidal energy available for ocean mixing approximately doubled.

The Coupled Ice Sheet–Earth System Model Bern3D v3.0

Abstract Understanding climate variability from millennial to glacial–interglacial time scales remains challenging due to the complex and nonlinear feedbacks between ice, ocean, sediments, biosphere, and atmosphere. Complex climate models generally struggle to dynamically and comprehensively simulate such long time periods as a result of the large computational costs. Here, we therefore coupled a dynamical ice sheet model to the Bern3D Earth system model of intermediate complexity, which allows for simulating multiple glacial–interglacial cycles. The performance of the model is first validated against modern observations and its response to abrupt perturbations, such as atmospheric CO2 changes and North Atlantic freshwater hosing, is investigated. To further test the fully coupled model, the climate evolution over the entire last glacial cycle is explored in a transient simulation forced by variations in the orbital configuration and greenhouse gases and aerosols. The model simulates global mean surface temperature in fair agreement with reconstructions, exhibiting a gradual cooling trend since the last interglacial that is interrupted by two more rapid cooling events during the early Marine Isotope Stage (MIS) 4 and Last Glacial Maximum (LGM). Simulated Northern Hemispheric ice sheets show pronounced variability on orbital time scales, and ice volume more than doubles from MIS3 to the LGM in good agreement with recent sea level reconstructions. At the LGM, the Atlantic overturning has a strength of about 14 Sv (1 Sv ≡ 106 m3 s−1), which is a reduction by about one-quarter compared to the preindustrial. We thus demonstrate that the new coupled model is able to simulate large-scale aspects of glacial–interglacial cycles.

Early–mid Holocene relative sea-level rise in the Yangtze River Delta, China

Regional-scale Holocene sea-level reconstruction is the key to understanding natural climatic variability. The tidal flat, salt marsh, and tidal floodplain in the Yangtze River Delta were very sensitive to morphological changes and sea-level variation during the early Holocene. In this study, the lithology, radiocarbon ages, sediment grain size, benthic foraminifera, and ostracods of three new cores were analyzed. Twenty-four sea-level index points were extracted from incised-valley fills beneath the westernmost part of the Yangtze River Delta and used to construct a detailed relative sea-level curve for 11.03–7.25 ka. The relative sea level were − 38.90 ± 3.48 and 1.59 ± 3.28 m at 11.03 and 7.25 ka, respectively, and the average rate of sea-level rise was 10.71 mm/yr. The relative sea level gradually increased at 11.03–10.10 and 9.29–8.33 ka from 4.80 to 13.30 and 3.19 to 19.59 mm/yr, respectively. The rate of relative sea-level rise gradually decreased from 19.59 mm/yr at 8.33 ka to 13.40 mm/yr at 7.82 ka and further decreased to 5.53 mm/yr at 7.25 ka. The Yangtze River Delta recorded accelerations in the rate of sea-level rise from 8.72 to 8.18 ka of 14.2–19.59 mm/yr, corresponding to the pre-8.2 ka event. A comparison of the sea-level histories of the Yangtze River Delta and other regions indicates tectonics, glacial isostatic adjustment, and coastal levering effects caused by the marine inundation of the continental shelves. These new sea-level data contribute to the understanding of the difference in the sea-level rise rate during the early–mid Holocene.

Paleocene–Eocene boundary in the Gabal El-Gir, east Qena, Egypt: Paleoenvironment and sequence stratigraphy interpretations

Comprehensive field, litho-, biostratigraphy, and paleoenvironmental studies of the upper Paleocene–lower Eocene (P–E) boundary interval in the Gabal El-Gir area were conducted. Two selected sections El-Gir-1 and El-Gir-2 were investigated. The prevailed paleoenvironmental conditions were discussed in the frame of the nannofossils, planktonic, and benthic foraminiferal assemblages. Benthic foraminifera assemblages indicate a remarkable turnover across the PETM. The extinction of Angulogavelinella avnimelechi at the base of the PETM represents a significant marker for the base of the benthic foraminifera extinction event at this site and other neritic sites in the Tethyan region. The dominance of calcareous planktonic excursion taxa at the base of the PETM indicates adaptation to extremely warm sea surface temperatures and low nutrient levels. Also, our data suggests that the biotic responses to environmental change in the early Eocene resembled those observed in the early Eocene hyperthermals at the nearbay Dababiya and El-Ballas sections. This could indicate that similar processes characterized the Egyptian sedimentary basin during the early Eocene. Stratigraphically, despite the limited distance between the studied sections, the Dababiya Quarry Member (DQM) in El-Gir-1 resemble those of the Global Stratotype Section and Point (GSSP) but reduced in thickness. While in El-Gir-2, the lower part of the DQM was missing. These findings could indicate that the DQM was deposited in a submarine channel. Three third-order depositional sequences were recognized covering the P–E interval. The reconstructed sea level curve shows the presence of eustatic similarities with the global sea level.

Unravelling the morphogenesis of coastal terraces at Cape Laundi (Sumba Island, Indonesia): Insights from numerical models

AbstractThe morphology of coastal sequences provides fundamental observations to unravel past sea level (SL) variations. For that purpose, converting morphometric observations into a SL datum requires understanding their morphogenesis. The long‐lasting sequence of coral reef terraces (CRTs) at Cape Laundi (Sumba Island, Indonesia) could serve as a benchmark. Yet, it epitomizes a pitfall that challenges the ultimate goal: the overall chronology of its development remains poorly constrained. The polycyclic nature of the terraces, involving marine erosion and reoccupation of old coral colonies by more recent ones hinders any clear assignment of Marine Isotope Stages (MIS) to specific terraces, in particular the reference datum corresponding to the last Interglacial maximum (i.e., MIS 5e). Thus, to overcome these obstacles, we numerically model the genesis of the sequence, testing a range of eustatic SL (ESL) reconstructions and uplift rates, as well as exploring the parameter space to address reef growth, erosion and sedimentation. A total of 625 model runs allowed us to improve the morpho‐chronological constraints of the coastal sequence and, more particularly, to explain the morphogenesis of the several CRTs associated with MIS 5e. Our results suggest that the lowermost main terrace was first constructed during the marine transgression of MIS 5e and was later reshaped during the marine regression of MIS 5e, as well as during the MIS 5c and MIS 5a highstands. Finally, we discuss the general morphology of the sequence and the implications it may have on SL reconstructions. At Cape Laundi, as elsewhere, we emphasize the necessity of addressing the development of CRT sequences with a dynamic approach, that is, considering that a CRT is a landform built continuously throughout the history of SL oscillations, and not simply during a singular SL maximum.