Search

Investigation of sea-level positions during the highly-dynamic Marine Isotope Stage 3 (MIS 3: 29–61 kyrs BP) proves difficult because: (i) in stable and subsiding areas, coeval coastal sediments are currently submerged at depths of few to several tens of meters below the present sea level, (ii) in uplifting areas, the preservation of geomorphic features and sedimentary records is limited due to the erosion occurred during the Last Glacial Maximum (LGM) with sea level at a depth of −130 m, followed by marine transgression that determined the development of ravinement surfaces. This study discusses previous research in the Mediterranean and elsewhere, and describes new fossiliferous marine deposits overlaying the metamorphic bedrock at Cannitello (Calabria, Italy). Radiocarbon ages of marine shells (about 43 kyrs cal BP) indicate that these deposits, presently between 28 and 30 m above sea level, formed during MIS 3.1. Elevation correction of the Cannitello outcrops (considered in an intermediate-to-far-field position with respect to the ice sheet) with the local vertical tectonic rate and Glacial Isostatic Adjustment (GIA) rate allows the proposal of a revision of the eustatic depth for this highstand. Our results are consistent with recently proposed estimates based on a novel ice sheet modelling technique.

This study presents a world review as well as new additional data in form of submerged speleothems that are used for paleo sea level reconstructions. Speleothems significantly contributed to the understanding of the global and regional sea level variations during the Middle and Late Quaternary. The studied speleothems cover the last 1.4 Myr and are focused mainly on Marine Isotope Stages (MIS) 1, 2, 3, 5.1, 5.3, 5.5, 7.1, 7.2, 7.3 and 7.5. Results reveal that submerged speleothems represent extraordinary archives providing detailed information on former sea level changes. We present also new results from stalactites collected in central Mediterranean sea, at Favignana and Ustica islands (Sicily, Italy), both characterized by continental, phreatic or marine layers. The study and analysis of the latter speleothems provide results of great interest for relative sea level changes over the last 1000 years.

The investigation of submerged speleothems for sea level studies has made significant contributions to the understanding of the global and regional sea level variations during the Middle and Late Quateranry. This has been especially the case for the Mediterranean Sea, where more than 300 submerged speleothems sampled in 32 caves have been analysed so far. Here, we present a comprehensive review of the results obtained from the study of submerged speleothems since 1978. The studied speleothems cover the last 1.4 Ma and are focused mainly on Marine Isotope Stages (MIS) 1, 2, 3, 5.1, 5.3, 5.5, 7.1, 7.2, 7.3 and 7.5. Results reveal that submerged speleothems represent extraordinary archives providing accurate information on former sea level changes, also considering that the Mediterranean Sea is devoid of any tropical corals since the Miocene. New results from a stalagmite collected at Palinuro (Campania, Italy) characterized by marine overgrowth are also reported. The measured elevations of speleothems are contaminated by the local response to glacial- and hydro-isostatic adjustment (GIA), and thus might significantly deviate from the global eustatic signal. Age and altitude comparation between Mediterranean speleothems, flowstone from Bahamas with local GIA provide a new scenario for MIS 5 and 7 sea level reconstrutions.

The investigation of submerged speleothems for sea level studies has made significant contributions to the understanding of the global and regional sea level variations during the Middle and Late Quaternary. This has especially been the case for the Mediterranean Sea, where more than 300 submerged speleothems sampled in 32 caves have been analysed so far. Here, we present a comprehensive review of the results obtained from the study of submerged speleothems since 1978. The studied speleothems cover the last 1.4 Myr and are mainly focused on Marine Isotope Stages (MIS) 1, 2, 3, 5.1, 5.3, 5.5, 7.1, 7.2, 7.3, and 7.5. The results reveal that submerged speleothems represent extraordinary archives providing accurate information on former sea level changes. New results from a stalagmite collected at Palinuro (Campania, Italy) and characterized by marine overgrowth are also reported. The measured elevations of speleothems are affected by the local response to glacial and hydro-isostatic adjustment (GIA), and thus might significantly deviate from the global eustatic signal. A comparison of the ages and altitude values of the Mediterranean speleothems and flowstone from the Bahamas with local GIA provides a new scenario for MIS 5 and 7 sea level reconstructions.

The evolution of past global ice sheets is highly uncertain. One example is the missing ice problem during the Last Glacial Maximum (LGM, 26 000-19 000 years before present) – an apparent 8-28 m discrepancy between far-field sea level indicators and modelled sea level from ice sheet reconstructions. In the absence of ice sheet reconstructions, researchers often use marine δ18O proxy records to infer ice volume prior to the LGM. We present a global ice sheet reconstruction for the past 80 000 years, called PaleoMIST 1.0, constructed independently of far-field sea level and δ18O proxy records. Our reconstruction is compatible with LGM far-field sea-level records without requiring extra ice volume, thus solving the missing ice problem. However, for Marine Isotope Stage 3 (57 000-29 000 years before present) - a pre-LGM period - our reconstruction does not match proxy-based sea level reconstructions, indicating the relationship between marine δ18O and sea level may be more complex than assumed., The configuration of past ice sheets, and therefore sea level, is highly uncertain. Here, the authors provide a global reconstruction of ice sheets for the past 80,000 years that allows to test proxy based sea level reconstructions and helps to reconcile disagreements with sea level changes inferred from models.

This study presents a global overview of the submerged speleothems used to reconstruct paleo sea levels and reports new results from two stalactites collected in the Mediterranean Sea. Coastal cave deposits significantly contributed to the understanding of global and regional sea-level variations during the Middle and Late Quaternary. The studied speleothems cover the last 1.4 Myr and focused mainly on Marine Isotope Stages (MIS) 1, 2, 3, 5.1, 5.3, 5.5, 7.1, 7.2, 7.3 and 7.5. The results indicate that submerged speleothems represent extraordinary archives that can provide detailed information on former sea-level changes. The two stalactites collected in the central Mediterranean Sea, at Favignana and Ustica islands (Sicily, Italy), are both characterized by continental, phreatic or marine layers. The U-Th and 14C ages of the new speleothems provide results of great interest for relative sea-level changes over the last 1000 years.

We use a standardized template for Pleistocene sea-level data to review last interglacial (Marine Isotope Stage 5 – MIS 5) sea-level indicators along the coasts of the western Atlantic and southwestern Caribbean, on a transect spanning from Brazil to Honduras and including the islands of Aruba, Bonaire, and Curaçao. We identified six main types of sea-level indicators (beach deposits, coral reef terraces, lagoonal deposits, marine terraces, Ophiomorpha burrows, and tidal notches) and produced 55 standardized data points, each constrained by one or more geochronological methods. Sea-level indicators are well preserved along the Brazilian coasts, providing an almost continuous north-to-south transect. However, this continuity disappears north of the Rio Grande do Norte Brazilian state. According to the sea-level index points (discrete past position of relative sea level in space and time) the paleo sea-level values range from ∼ 5.6 to 20 m above sea level (a.s.l.) in the continental sector and from ∼ 2 to 10 m a.s.l. in the Caribbean islands. In this paper, we address the uncertainties surrounding these values. From our review, we identify that the coasts of northern Brazil, French Guiana, Suriname, Guyana, and Venezuela would benefit from a renewed study of Pleistocene sea-level indicators, as it was not possible to identify sea-level index points for the last interglacial coastal outcrops of these countries. Future research must also be directed at improving the chronological control at several locations, and several sites would benefit from the re-measurement of sea-level index points using more accurate elevation measurement techniques. The database compiled in this study is available in spreadsheet format at the following link: https://doi.org/10.5281/zenodo.5516444 (Version 1.02; Rubio-Sandoval et al., 2021).

The timing, duration and evolution of sea level during the Marine Isotope Stage 5e (MIS 5e) highstand is a subject of intense debate. A major problem in resolving this debate is the difficulty of chronologically constraining the sea level fall that followed the peak of the highstand. This was mainly controlled by ice-sheet dynamics, the understanding of which is relevant for assessing future sea-level behavior due to global warming. Here we use stratigraphical and geochoronological (U/Th dating and tephra fingerprinting) evidence from the Infreschi archaeological cave (Marina di Camerota, Southern Italy) to constrain relative sea level (RSL) evolution during the MIS 5e highstand and younger stages. Uranium-thorium dating of speleothem deposition phases places the maximum highstand RSL at 8.90 ± 0.6m a.s.l., as indicated by the near-horizontal upper limit of Lithophaga boreholes measured for along a ∼3.5 km coastal cliff section. Geochronological data show that RSL fell more than 6 m before ∼120 ka, suggesting a duration of the Last Interglacial highstand significantly shorter than proposed in some previous studies. Modelling shows that the RSL trend predicted by the ICE-5G and ICE-6G ice-sheet simulations is consistent with our data, but requires an additional significant reduction of both Greenland and Antarctic ice sheets to match the height of the local maximum highstand if no correction for tectonics is applied. Reconciling field data and models requires an earlier and likely shorter duration of the MIS 5e highstand. This suggests that our new data can constrain global ice-volume variations during the penultimate deglaciation, as well as glacial inception at the end of the Last Interglacial. According to our chronology, there is no local evidence of higher-than-present-day sea levels after 120 ka.

This paper presents the current state-of-knowledge of the New Zealand (Aotearoa) last interglacial (MIS 5 sensu lato) sea-level record compiled within the framework of the World Atlas of Last Interglacial Shorelines (WALIS) database. Seventy-seven total relative sea-level (RSL) indicators (direct, marine-, and terrestrial-limiting points), commonly in association with marine terraces, were identified from over 120 studies reviewed. Extensive coastal deformation around New Zealand has resulted in a significant range of elevation measurements on both the North Island (276.8 to −94.2 msl) and South Island (173.1 to −70.0 msl) and prompted the use of RSL indicators to estimate rates of vertical land movement; however, indicators lack adequate description and age constraint. Identified RSL indicators are correlated with MIS 5, MIS 5e, MIS 5c, and MIS 5a and indicate the potential for the New Zealand sea-level record to inform sea-level fluctuation and climatic change within MIS 5 (sensu lato). The Northland (North Island) and Otago (South Island) regions, historically considered stable, have the potential to provide a regional sea-level curve in a remote location of the South Pacific across broad degrees of latitude. Future work requires modern analogue information, heights above a defined sea-level datum, better stratigraphic descriptions, and use of improved geochronological methods. The database presented in this study is available open-access at this link: http://doi.org/10.5281/zenodo.4056376 (Ryan et al., 2020a).

Ongoing work shows that species richness patterns on volcanic oceanic islands are shaped by surface area changes driven by longer time scale (>1 ka) geological processes and natural sea level fluctuations. A key question is: what are the rates and magnitudes of the forces driving spatial changes on volcanic oceanic islands which in turn affect evolutionary and biogeographic processes? We quantified the rates of surface-area changes of a whole island resulting from both volcanogenic flows and sea level change over the last glacial-interglacial (GI) cycle (120 ka) for the volcanically active island of Terceira, (Azores, Macaronesia, Portugal). Volcanogenic activity led to incidental but long-lasting surface area expansions by the formation of a new volcanic cone and lava-deltas, whereas sea level changes led to both contractions and expansions of area. The total surface area of Terceira decreased by as much as 24% per time step due to changing sea levels and increased by 37% per time step due to volcanism per time step of 10 ka. However, while sea levels nearly continuously changed the total surface area, volcanic activity only impacted total surface area during two time steps over the past 120 ka. The surface area of the coastal and lowland region (here defined as area

We surveyed the elevation and age (14C) of paleo sea-level indicators in five islands of the Spermonde Archipelago. We describe 24 new sea-level index points from fossil microatolls, and we compare our dataset with both previously published proxies and sea-level predictions from a set of 54 Glacial Isostatic Adjustment (GIA) models, using different assumptions on both ice melting histories and mantle structure and viscosity. We then investigate the implications of our data and models in terms of vertical land movements in the study area, with two main results.First, data from the heavily populated island of Barrang Lompo are significantly lower (ca. 80 cm) than those at all the other islands. In absence of instrumental data (e.g., GPS or tide gauges) in any of the islands, we advance the hypothesis that this difference may be due to groundwater extraction and loading of buildings on Barrang Lompo, that would cause this island to subside at rates in the order of ~3-11 mm/a.Second, Common Era data (0-400 a BP) seem to indicate that the islands in the archipelago may be affected by tectonically-driven vertical land motions in the order of -0.88±0.61 mm/a (1-sigma), albeit slight uplift cannot be excluded. Different assumptions on vertical land motions affect, in turn, the assessment of which GIA model shows the best match with Late Holocene (ca. 4-5 ka) sea level data. Tectonic stability or slight uplift would favor iterations of ANICESELEN (De Boer et al., 2014), while subsidence would cause the sea level data to fit better with iterations of ICE-6G (Peltier et al., 2015).ReferencesDe Boer, Bas, Paolo Stocchi, and Roderik Van De Wal. A fully coupled 3-D ice-sheet-sea-level model: algorithm and applications." Geoscientific Model Development 7.5 (2014): 2141-2156.Peltier, W. R., D. F. Argus, and R. Drummond. Space geodesy constrains ice age terminal deglaciation: The global ICE‐6G_C (VM5a) model. Journal of Geophysical Research: Solid Earth 120.1 (2015): 450-487.AcknowledgmentsThis project is funded by SEASCHANGE (RO-5245/1-1) and HAnsea (MA-6967/2-1) from the Deutsche Forschungsgemeinschaft (DFG), part of the Special Priority Program (SPP)-1889 "Regional Sea Level Change and Society". Parts of this study are under review in Climate of the Past (https://www.clim-past-discuss.net/cp-2019-63/)

Indonesia is a country composed of several thousand islands, many of them small, low-lying and densely inhabited. These are, in particular, subject to high risk of inundation due to future relative sea level changes. The Spermonde Archipelago, off the coast of Southwest Sulawesi, consists of more than 100 small islands. This study presents a dataset of 24 sea-level index points from fossil microatolls, surveyed on five islands in the Spermonde Archipelago and compares these new results with published data from the same region and with relative sea level predictions from different Glacial Isostatic Adjustment (GIA) models. The newly surveyed fossil microatolls are located around the islands of Tambakulu, Suranti (both ~ 60 km offshore of Makassar city), Bone Batang and Kodingareng Keke (both located in the center of the Archipelago) and Sanrobengi (located ~ 20 km south-southwest of Makassar). Results from the near- and mid-shelf islands indicate that relative sea level between 4 to 6 ka BP was less than one meter above present sea level. The only exception to this pattern is the heavily populated island of Barrang Lompo, where we record a significant subsidence when compared to the other islands. These new results support the conclusions from a previous dataset and are relevant to constrain late Holocene ice melting scenarios. Samples from the two outer islands (Tambakulu and Suranti) yielded ages spanning the Common Era that represent, to our knowledge, the first reported for the entire Southeast Asian region.

Sea level changes can be driven by either variations in the masses or volume of the oceans, or by changes of the land with respect to the sea surface. In the first case, a sea level change is defined ‘eustatic’; otherwise, it is defined ‘relative’. Several techniques can be used to observe changes in sea level, from satellite data to tide gauges to geological or archeological proxies. Regardless of the technique used, ‘eustasy’ cannot be measured directly, but only calculated after perturbing factors of different origins are taken into account. In this paper, we review the meaning and main processes that contribute to eustatic and relative sea level changes, and we give an overview of the different techniques used to observe them.

We assembled a database of Holocene relative sea-level index points (n = 213) and marine (n = 211) and terrestrial (n = 122) limiting points for the broader South and Southeast Asian region including the Maldives, India and Sri Lanka. The standardized review of published age-elevation information from corals, deltaic, estuarine and mangrove deposits, beachrocks and tidal notches, yielded a new suite of relative sea-level index and limiting points produced according to a standardized protocol. Expected spatial variability in Holocene relative sea-level change due to glacial isostatic adjustment was accounted for, by first subdividing the study area into ten geographic sub-regions from the Central Indian Ocean to the Western Tropical Pacific, and second by comparing sub-regional relative sea-level data to model predictions of glacial isostatic adjustment. Results show that some of the regionally constrained relative sea-level data are characterized by significant inconsistencies that cannot be explained by glacial isostatic adjustment. Such inconsistencies of standardized relative sea-level data become particularly obvious in areas around the Red River Delta in Vietnam, the Gulf of Thailand, the northwest coast of Malaysia and the Spermonde Archipelago in Indonesia. Based on a critical evaluation of the reviewed relative sea-level indicators, we discuss possible sources of local divergence and identify regions where data are currently insufficient to constrain glacial isostatic adjustment predictions. The remaining quality-controlled and consistent relative sea-level data show that glacial isostatic adjustment and syn-/post-formational influences such as tectonic uplift, subsidence and compaction were the dominant local drivers of Holocene relative sea-level change. Collectively, the results of this review suggest that Holocene sea levels in South and Southeast Asia and surrounding regions have been controlled by a variety of global and local drivers and imply that additional index points from the Java Sea in Indonesia would be valuable to better assess the spatial variability, and to calibrate geophysical models of glacial isostatic adjustment. Published version

Earth is heading towards a climate that last existed more than three million years ago (Ma) during the ‘mid-Pliocene warm period’1, when atmospheric carbon dioxide concentrations were about 400 parts per million, global sea level oscillated in response to orbital forcing2,3 and peak global-mean sea level (GMSL) may have reached about 20 metres above the present-day value4,5. For sea-level rise of this magnitude, extensive retreat or collapse of the Greenland, West Antarctic and marine-based sectors of the East Antarctic ice sheets is required. Yet the relative amplitude of sea-level variations within glacial–interglacial cycles remains poorly constrained. To address this, we calibrate a theoretical relationship between modern sediment transport by waves and water depth, and then apply the technique to grain size in a continuous 800-metre-thick Pliocene sequence of shallow-marine sediments from Whanganui Basin, New Zealand. Water-depth variations obtained in this way, after corrections for tectonic subsidence, yield cyclic relative sea-level (RSL) variations. Here we show that sea level varied on average by 13 ± 5 metres over glacial–interglacial cycles during the middle-to-late Pliocene (about 3.3–2.5 Ma). The resulting record is independent of the global ice volume proxy3 (as derived from the deep-ocean oxygen isotope record) and sea-level cycles are in phase with 20-thousand-year (kyr) periodic changes in insolation over Antarctica, paced by eccentricity-modulated orbital precession6 between 3.3 and 2.7 Ma. Thereafter, sea-level fluctuations are paced by the 41-kyr period of cycles in Earth’s axial tilt as ice sheets stabilize on Antarctica and intensify in the Northern Hemisphere3,6. Strictly, we provide the amplitude of RSL change, rather than absolute GMSL change. However, simulations of RSL change based on glacio-isostatic adjustment show that our record approximates eustatic sea level, defined here as GMSL unregistered to the centre of the Earth. Nonetheless, under conservative assumptions, our estimates limit maximum Pliocene sea-level rise to less than 25 metres and provide new constraints on polar ice-volume variability under the climate conditions predicted for this century. Sea level varied by 13 ± 5 metres on average, but up to 25 metres, over glacial–interglacial cycles during the Pliocene epoch, due to partial collapses of Antarctic Ice Sheets.

As global climate warms and sea level rises, coastal areas will be subject to more frequent extreme flooding and hurricanes. Geologic evidence for extreme coastal storms during past warm periods has the potential to provide fundamental insights into their future intensity. Recent studies argue that during the Last Interglacial (MIS 5e, ∼128–116 ka) tropical and extratropical North Atlantic cyclones may have been more intense than at present, and may have produced waves larger than those observed historically. Such strong swells are inferred to have created a number of geologic features that can be observed today along the coastlines of Bermuda and the Bahamas. In this paper, we investigate the most iconic among these features: massive boulders atop a cliff in North Eleuthera, Bahamas. We combine geologic field surveys, wave models, and boulder transport equations to test the hypothesis that such boulders must have been emplaced by storms of greater-than-historical intensity. By contrast, our results suggest that with the higher relative sea level (RSL) estimated for the Bahamas during MIS 5e, boulders of this size could have been transported by waves generated by storms of historical intensity. Thus, while the megaboulders of Eleuthera cannot be used as geologic proof for past “superstorms,” they do show that with rising sea levels, cliffs and coastal barriers will be subject to significantly greater erosional energy, even without changes in storm intensity.

Paleo relative sea-level (RSL) indicators formed during the Marine Isotope Stage (MIS) 5e have been reported by a large number of studies worldwide. Despite this, three main aspects are seldom reported: (1) use of high-precision survey techniques applied to MIS 5e RSL indicators; (2) application of modern analogs to understand the indicative meaning of MIS 5e RSL indicators; (3) estimates of the effects of glacial isostatic adjustment (GIA) on the MIS 5e records. In this study, we show how the three points above have been addressed in a focused study on Last Interglacial outcrops on the island of Mallorca. We measured the elevation of several RSL indicators with high-accuracy differential GPS (vertical accuracies down to 0.1 m) and we established the relationship between each RSL indicator and the paleo sea level through calculation of the indicative meaning for each RSL indicator. In particular, we present a novel technique to calculate the indicative meaning of fossil beach deposits with a phase-averaged morphodynamic model (CSHORE). We show how this approach helps overcoming difficulties with the survey of the modern analogs for these indicators. Our results show that two paleo RSLs are imprinted in Mallorca at + 2.9 ± 0.8 m and + 11.3 ± 1.0 m. We then compare our field-based results with modelled paleo RSL, calculated from the predictions of the ice-earth coupled ANICE-SELEN model, using few different ice-sheet melting scenarios during MIS 5e. We conclude that indicative ranges can be derived from relatively simple morphodynamic models and that the comparison of field-derived and modelled RSL values is a good method to validate possible scenarios of MIS 5e sea-level variability, especially in absence of precise dating.

This paper provides new relative sea level data inferred from coastal archaeological sites located along the Turkish coasts of the Gulf of Fethye (8 sites), and Israel, between Akziv and Caesarea (5 sites). The structures selected are those that, for effective functioning, can be accurately related to sea level at the time of their construction. Thus their positions with respect to present sea level provide a measure of the relative sea level change since their time of construction. Useful information was obtained from the investigated sites spanning an age range of ˜2.3–˜1.6 ka BP. The inferred changes in relative sea level for the two areas are distinctly different, from a rise of 2.41 to 4.50 m in Turkey and from 0 to 0.18 m in Israel. Sea level change is the combination of several processes, including vertical tectonics, glacio–hydro-isostatic signals associated with the last glacial cycle, and changes in ocean volume. For the Israel section, the present elevations of the MIS-5.5 Tyrrhenian terraces occur at a few meters above present sea level and vertical tectonic displacements are small. Data from GPS and tide gauge measurements also indicate that any recent vertical movements are small. The MIS-5.5 shorelines are absent from the investigated section of the Turkish coast, consistent with crustal subsidence associated with the Hellenic Arc. The isostatic signals for the Israel section of the coast are also small (ranging from −0.11 mm/yr to 0.14 mm/yr, depending on site and earth model) and the observed (eustatic) average sea level change, corrected for this contribution, is a rise of 13.5 ± 2.6 cm during the past ˜2 ka. This is attributed to the time-integrated contribution to sea level from a combination of thermal expansion and other increases in ocean volume. The observed sea levels from the Turkish sites, in contrast, indicate a much greater rise of up to 2.2 mm/yr since 2.3 ka BP occurring in a wide area between Knidos and Kekova. The isostatic signal here is also one of a rising sea level (of up to ˜1 mm/yr and site and earth-model dependent) and the corrected tectonic rate of land subsidence is ˜1.48 mm/yr. This is the primary cause of dramatic relative sea level rise for this part of the coast.

Vertical movements in the Po plain (northern Italy) are controlled by natural and anthropogenic effects. Since Italy is located in the far–field of the former late Pleistocene ice sheets, isostatic deformations are primarily driven by melt water loading and represent a major component of long–term natural movements across the entire Mediterranean. In addition to far–field sources, here we consider the isostatic effects of melting of the nearby Wurm Alpine ice–sheet, suggesting that it is possible to put bounds on its maximum thickness, extent and chronology by Holocene relative sea level observations from the northern Adriatic. Using various plausible ice models, and adopting a viscosity profile that matches Holocene relative sea level observations in the Mediterranean, we find that melting of the Alpine ice sheet is always responsible for upward movements in the Po plain, currently at rates of ~ 0.5 mm/yr. When both far– and near–field sources are considered, the rate of sea level change in the Venetian Lagoon for the most reasonable mantle rheology and melting chronology is negative, i.e., opposite to that attributed to human activity and recent climatic variations. However, its amplitude (fractions of mm/yr) is small compared to the secular signal observed by tide gauges (~2 mm/yr), which makes glacial isostasy a second–order mechanism of sea level variation in this region.

Knowledge on the timing and magnitude of past sea-level changes is essential to understand modern and future sea-level variability. Holocene sea-level data from literature on the west coast of Sulawesi, central Indonesia, suggest that this region experienced two relative sea-level highstands over the last 6000 years, with magnitudes exceeding two meters. However, recent datasets from the Indo-Pacific region do not support high-magnitude sea-level oscillations during the Holocene in tectonically stable far-field locations. Here we present a new, high-precision, mid-Holocene sea-level dataset from the Spermonde Shelf off southwest Sulawesi. We surveyed 21 fossil microatolls on the reef flats of two coral islands (Pulau Panambungan and Pulau Barrang Lompo) and referred their elevations to local mean sea level and to the height of living coral. Radiometrically calibrated ages from emergent fossil microatolls on Pulau Panambungan indicate a relative sea-level highstand not exceeding 0.5 m above present at ca. 5600 cal. yr BP. The highstand is followed by a relatively rapid sea-level fall towards present sea level that was reached at around 4000 cal. yr BP. Fossil microatolls from nearby Pulau Barrang Lompo show the same trend, however with a coherent negative vertical offset of about 0.8 m compared to their equivalents on Pulau Panambungan. The largely consistent gradients of both trends (~− 0.14 mm yr− 1), the consistent elevation of living microatolls in the Spermonde, and a number of instructive geomorphic features indicate a localized, post-formational and probably recent drop of the fossil microatolls on the densely populated island Pulau Barrang Lompo. The relative sea-level trend inferred from Pulau Panambungan is well within the range of geophysical predictions based on ANICE-SELEN ice sheet model, which predict a highstand that is significantly lower than those predicted by other GIA models for this area. Although a complete interpretation of the Holocene sea-level history will require additional high-resolution datasets from this and surrounding territories in SE Asia, our results suggest that there was merely a single Holocene highstand in central Indonesia, the magnitude of which was substantially lower than hitherto assumed.

International audience; The Last Interglacial (MIS 5e, 128-116 ka) is among the most studied past periods in Earth's history. The climate at that time was warmer than today, primarily due to different orbital conditions, with smaller ice sheets and higher sea-level. Field evidence for MIS 5e sea-level was reported from thousands of sites, but often paleo shorelines were measured with low-accuracy techniques and, in some cases, there are contrasting interpretations about paleo sea-level reconstructions. For this reason, large uncertainties still surround both the maximum sea-level attained as well as the pattern of sea-level change throughout MIS 5e. Such uncertainties are exacerbated by the lack of a uniform approach to measuring and interpreting the geological evidence of paleo sea-levels. In this review, we discuss the characteristics of MIS 5e field observations, and we set the basis for a standardized approach to MIS 5e paleo sea-level reconstructions, that is already successfully applied in Holocene sea-level research. Application of the standard definitions and methodologies described in this paper will enhance our ability to compare data from different research groups and different areas, in order to gain deeper insights into MIS 5e sea-level changes. Improving estimates of Last Interglacial sea-level is, in turn, a key to understanding the behavior of ice sheets in a warmer world. (C) 2016 Elsevier B.V. All rights reserved.

The present-day sea level variations and geodetically observed ground deformations in the Mediterranean area are normally ascribed to the combined effect of tectonic or human-driven subsidence and postglacial uplift as a result of the melting of the major Pleistocene ice sheets. However, another potential cause of deformation, only marginally considered to date, is the melting of the glacier that covered the Alps during the last glacial maximum (LGM). The aim of this paper is to predict the long-term sea level variations induced by the melting of both the late-Pleistocene and Alpine ice sheets and compare our results with the relative sea level (RSL) observations available in the Mediterranean region. This task is accomplished solving the sea level equation (SLE) for a spherically symmetric viscoelastic Earth. Our analysis shows that the melting of the Alpine glacier has marginally affected the Holocene sea level variations in the near-field sites in southern France (Marseilles and Roussillon) and the central Tyrrhenian sea (Civitavecchia), and that the RSL predictions are significantly sensitive to the chronology of the remote ice aggregates. The computations, which are performed using a specific mantle viscosity profile consistent with global observations of RSL rise, show that the uplift rate driven by the Alpine isostatic readjustment may account for up to 1/3 or the rates observed at GPS stations in the western portion of the chain. Our results suggest that a thorough modelization of both near- and far-field ice sheets is necessary to gain a better insight into the present-day deformations and sea level variations in the Mediterranean region.

Precisely quantifying the current climate-related sea level change requires accurate knowledge of long-term geological processes known as Glacial Isostatic Adjustments (GIA). Although the major postglacial melting phase is likely to have ended ∼6–4 ka BP (before present), GIA is still significantly affecting the present-day vertical position of the mean sea surface and the sea bottom. Here we present empirical rsl (relative sea level) data based on U/Th dated fossil corals from reef platforms of the Society Islands, French Polynesia, together with the corresponding GIA-modeling. Fossil coral data constrain the timing and amplitude of rsl variations after the Holocene sea level maximum (HSLM). Upon correction for isostatic island subsidence, we find that local rsl was at least ∼1.5 ± 0.4 m higher than present at ∼5.4 ka. Later, minor amplitude variations occurred until ∼2 ka, when the rsl started dropping to its present position with a rate of ∼0.4 mm/yr. The data match with predicted rsl curves based on global ice-sheet chronologies confirming the role of GIA-induced ocean siphoning effect throughout the mid to late Holocene. A long lasting Late Holocene highstand superimposed with second-order amplitudinal fluctuations as seen from our data suggest that the theoretical predicted timing of rsl change can still be refined pending future calibration.

Models for glacial isostatic adjustment (GIA) can provide constraints on rheology of the mantle if past ice thickness variations are assumed to be known. The Pleistocene ice loading histories that are used to obtain such constraints are based on an a priori 1-D mantle viscosity profile that assumes a single deformation mechanism for mantle rocks. Such a simplified viscosity profile makes it hard to compare the inferred mantle rheology to inferences from seismology and laboratory experiments. It is unknown what constraints GIA observations can provide on more realistic mantle rheology with an ice history that is not based on an a priori mantle viscosity profile. This paper investigates a model for GIA with a new ice history for Fennoscandia that is constrained by palaeoclimate proxies and glacial sediments. Diffusion and dislocation creep flow law data are taken from a compilation of laboratory measurements on olivine. Upper-mantle temperature data sets down to 400 km depth are derived from surface heatflow measurements, a petrochemical model for Fennoscandia and seismic velocity anomalies. Creep parameters below 400 km are taken from an earlier study and are only varying with depth. The olivine grain size and water content (a wet state, or a dry state) are used as free parameters. The solid Earth response is computed with a global spherical 3-D finite-element model for an incompressible, self-gravitating Earth. We compare predictions to sea level data and GPS uplift rates in Fennoscandia. The objective is to see if the mantle rheology and the ice model is consistent with GIA observations. We also test if the inclusion of dislocation creep gives any improvements over predictions with diffusion creep only, and whether the laterally varying temperatures result in an improved fit compared to a widely used 1-D viscosity profile (VM2). We find that sea level data can be explained with our ice model and with information on mantle rheology from laboratory experiments, heatflow and seismology and a pure olivine rheology above 400 km. Moreover, laterally heterogeneous models provide a significantly better fit to relative sea level data than the VM2 viscosity, for our ice model as well as for the ICE-5G model that is based on the VM2 profile. The new ice model gives different constraints on mantle rheology than the ICE-5G model, indicating a possible bias towards mantle viscosity in the latter or shortcomings in our ice model. Present-day uplift rates for a dry rheology are close to GPS observed uplift rate for certain combinations of grain size and temperature fields. Sea level data show a preference for a wet olivine rheology, but in that case uplift rates are too low for all grain sizes and temperature fields. The difficulty to fit sea level data and uplift rate data simultaneously can not be resolved by varying creep parameters below 400 km. Uncertainties in the flow law and the neglect of other materials in the upper mantle, as well as the neglect of flow in the crust could affect our conclusions.

Two-dimensional reconstructions of the Mediterranean sea level corrected for the atmospheric effects are proposed at monthly interval over the period 1970-2006 using 14 tide gage records and 33-year long (1970-2002) sea level grids from the NEMOMED8 regional ocean circulation model (NM8) and the PROTHEUS System Atmosphere-Ocean coupled model (PS). They are compared with a similar reconstruction using decade-long sea level grids from altimetry (Topex/Poseidon and Jason1) and a published reconstruction by Calafat and Gomis (2009). Tests with extra tide gages, not used in the computation, show that interannual variability is better captured when using long (33-year) spatial grids. In particular the NM8-based reconstruction reproduces better the sea level variability at all independent tide gages. An empirical Orthogonal Function decomposition of this reconstruction over 1970-2006 shows that the temporal curve of the two first modes are highly correlated with the North Atlantic Oscillation. We note in particular different behaviors over the 1970-1994 and 1994-2006 time spans. Results suggest that the North Atlantic Oscillation forcing modified the spatial patterns of the Mediterranean sea level around the year 1993 close to the date of occurrence of the Eastern Mediterranean Transient (a major change in the deep water formation of the Levantine and Ionian basin that occurred in the early 1990s). © 2011 Elsevier B.V., We would like to acknowledge the support and financial contribution of the Topo-Europe programme and CECILE programme.

Published and new data exist for relative sea-level change for 105 locations (127 samples) during the late Holocene, along the Italian (and Istrian) coasts. These data, compared with predictions (derived from two different models associated with the last glacial cycle) allowed the calculation of the tectonic vertical movements. They are based on precise measures of geomorphological and archaeological markers between 0.4 and 12.6 ka cal. BP, sampled at elevations between +7 and −51 m. In order to decipher the broad pattern of Holocene tectonic vertical movements along the Italian coastline, these data were compared with predicted sea-level curves using the most recent models published for the Mediterranean sea. Tectonic rates varied from −4.85 mm/a, in a core at Sybaris, to 5 mm/a, in the volcanic areas of Pozzuoli and Pantelleria. New MIS 5.5 (125 ka) data, mostly from the Venetian plain, are reported. In particular the depth of the base of MIS 5.5 paralic deposits found in four cores near Venezia provides a mean subsidence of 0.62 mm/a. New, precise mass spectrometer U-Th analyses on Cladocora layers from the bottom of a long core (named ENEA), indicate older ages (195.7 ± 1.6 and 161.2 ± 1.2 ka, respectively), relative to the published MIS 5.5 ages, which were based on alpha-counting U-Th data. Instrumental data obtained from tide gauges and repeated levelling measurements from the NE Adriatic and Sicily are also considered. These methods have one great advantage with respect to continuous GPS measurements and the satellite altimetric observations, in that a much greater time span is available. Although the altimetric measurements are available for 16 years, and the GPS for less than a decade, repeated levelling lines cover up to 50 years and tide gauge observations in some cases to 100 years or more. The greater time span allows for more stable differential rate estimates. The repeated levelling shows that the plain east of Mestre is subsiding (to −4 mm/a). The Messina tidal gauge demonstrates a total coseismic and post-seismic subsidence of 77 cm associated with the event of 1908, the post-seismic phase lasting for at least 13 years. The Reggio Calabria tidal station points to an uplift of this station relative to Palermo in the order of 1–2 mm/a.

Solving the sea–level equation for a Maxwell Earth, we analyze the sensitivity of Holocene sea–level records in SE Tunisia to the time–history of remote ice sheets. Assuming that mantle viscosity increases moderately with depth, we find that in this region the sea–level variations driven by the Northern Hemisphere ice sheets cancel, so that the late–Holocene sea–level high–stand suggested by the geological record merely reflects the melting history of Antarctica. New insight into the history of this ice sheet is obtained analyzing the information contained in a revised set of relative sea–level observations for sites across the Mediterranean covering the last 8 kyrs. From a trial–and–error misfit analysis, it holds true that in this region the match between model predictions and observations improves when the volume of water released from Antarctica is well below the value imposed by the ICE3G chronology and when a sudden meltwater pulse is allowed between 8 and 7 kyrs before present, corresponding to the epoch of the catastrophic rise event known as CRE3.

This paper presents the current state of knowledge of the Aotearoa New Zealand last interglacial (marine isotope stage 5, MIS 5, sensu lato) sea-level record compiled within the framework of the World Atlas of Last Interglacial Shorelines (WALIS) database. A total of 77 relative sea-level (RSL) indicators (direct, marine-limiting, and terrestrial-limiting points), commonly in association with marine terraces, were identified from over 120 studies reviewed. Extensive coastal deformation around New Zealand has prompted research focused on active tectonics, the scale of which overprints the sea-level record in most regions. The ranges of last interglacial palaeo-shoreline elevations are significant on both the North Island (276.8 ± 10.0 to −94.2 ± 10.6 ma.m.s.l., above mean sea level) and South Island (165.8 ± 2.0 to −70.0 ± 10.3 ma.m.s.l.) and have been used to estimate rates of vertical land movement; however, in many instances there is a lack of adequate description and age constraint for high-quality RSL indicators. Identified RSL indicators are correlated with MIS 5, MIS 5e, MIS 5c, and MIS 5a and indicate the potential for the New Zealand sea-level record to inform sea-level fluctuation and climatic change within MIS 5. The Northland Region of the North Island and southeastern South Island, historically considered stable, have the potential to provide a regional sea-level curve, minimally impacted by glacio- and hydro-isostatic adjustment (GIA) and reflecting near-eustatic fluctuations in a remote location of the South Pacific, across broad degrees of latitude; however, additional records from these regions are needed. Future work requires modern analogue information, heights above a defined sea-level datum, better stratigraphic descriptions, and use of improved geochronological methods. The database presented in this study is available open access at this link: https://doi.org/10.5281/zenodo.4590188 (Ryan et al., 2020a).