Search

The warmest peak of the Last Interglacial (ca. 128–116 ka) is considered a process analogue, and is often studied to better understand the effects of a future warmer climate on the Earth's system. In particular, significant effort has been made to better constrain ice sheet contributions to peak Last Interglacial sea level through field observation of paleo relative sea level indicators. Along tropical coastal margins, these observations are predominantly based on fossil shallow coral reef sequences, also thanks to the possibility of gathering reliable U-series chronological constraints. However, the preservation of many Pleistocene reef sequences is often limited to a series of discrete relative sea-level positions within the interglacial, where corals suitable for dating were preserved. This in turn, limits our ability to understand the continuous evolution of paleo relative sea-level through an entire interglacial, also affecting the possibility to unravel the existence and pattern of sub-stadial sea level oscillations. While the interpretation of lithostratigraphic and geomorphologic properties is often used to overcome this hurdle, geological interpretation may present issues related to subjectivity when dealing with missing facies or incomplete sequences. In this study, we try to step back from a conventional approach generating a spectrum of synthetic Quaternary subtropical fringing reefs for a site in southwestern Madagascar (Indian Ocean). We use the DIONISOS forward (Beicip Franlab) stratigraphic model to build a fossil reef at this location. In each model run, we use distinct Greenland and Antarctica Ice Sheet melt scenarios produced by a coupled ANICE-SELEN glacial isostatic adjustment model. The resulting synthetic reef sequences are then used test these melt scenarios against the stratigraphic record. We propose that this sort of stratigraphic modelling may provide further quantitative control when interpreting Last Interglacial reef sequences.

Understanding sea level during the warmest peak of the Last Interglacial (125,000 yrs ago; Marine Isotope Stage 5e) is important for assessing future ice-sheet dynamics in response to climate change, and relies on the measurement and interpretation of paleo sea-level indicators, corrected for post-depositional vertical land motions. The coasts and continental shelves of northeastern Australia (Queensland) preserve an extensive Last Interglacial record in the facies of coastal strandplains onland and fossil reefs offshore. However, there is a discrepancy (amounting to tens of meters) in the elevation of sea-level indicators between offshore and onshore sites. Here, we assess the influence of geophysical processes that may have changed the elevation of these sea-level indicators since the Last Interglacial. We modeled sea-level change due to: i) dynamic topography; ii) glacial isostatic adjustment, and iii) isostatic adjustment due to coral reef loading, which we term "reef isostasy". These processes caused relative sea-level changes on the order of, respectively, 10 m, 5 m, and 0.3 m since the Last Interglacial. Of these geophysical processes, the dynamic topography predictions most closely match the tilting observed between onshore and offshore sea-level markers. However, we found that these combined geophysical processes cannot explain the full amplitude of the observed discrepancy between onshore and offshore sea-level indicators.

The Last Interglacial (LIG; ca. 125 ka) is a period of interest for climate research as it is the most recent period of the Earth's history when the boreal climate was warmer than at present. Previous research, based on models and geological evidence, suggests that the LIG may have featured enhanced patterns of ocean storminess, but this remains hotly debated. Here, we apply state-of-the-art climate and hydrodynamic modeling to simulate changes in sea level extremes caused by storm surges, under LIG and pre-industrial climate forcings. Significantly higher seasonal LIG sea level extremes emerge for coastlines along northern Australia, the Indonesian archipelago, much of northern and eastern Africa, the Mediterranean Sea, the Gulf of Saint Lawrence, the Arabian Sea, the east coast of North America, and islands of the Pacific Ocean and of the Caribbean. Lower seasonal LIG sea level extremes emerge for coastlines along the North Sea, the Bay of Bengal, China, Vietnam, and parts of Central America. Most of these anomalies are associated with anomalies in seasonal sea level pressure minima and in eddy kinetic energy calculated from near-surface wind fields, and therefore seem to originate from anomalies in the meridional position and intensity of the predominant wind bands. In a qualitative comparison, LIG sea level extremes seem generally higher than those projected for future warmer climates. These results help to constrain the interpretation of coastal archives of LIG sea level indicators.

A 70 km laterally extending tidal notch was carved during the Last Interglacial (127-116 ka BP) along the East-Sardinian coast, Orosei Gulf, with a maximum height of 9.1 m. We present the results of the detailed mapping of the Orosei tidal notch using the Structure from Motion Multi-View Stereo (SfM-MVS) reconstruction method. Although its geometry is laterally constant, the notch depth differs due to local factors. Results of SfM-MVS reconstructions, together with local parameters, were used as input in a geometric model that simulates the notch shape based on randomly sampled relative sea-level curves. The modeled profiles best matching the real shape of the notch are those characterized by a bimodal geometry. Extracting information from this rocky imprint is considered crucial to better understand past changes in sea level, which are in turn significant to better gauge sea-level variations in warmer climate conditions., https://app.dimensions.ai/details/grant/grant.7819993

Several physical, chemical, and biological processes shape coastal environments close to sea level. Acting through time, these processes create a variety of coastal landforms. When found outside their environment of formation, these landforms can be used by geoscientists as geomorphological indicators of former relative sea levels. In this chapter, we outline the main processes acting on coastal areas, and link them to distinct types of sea-level indicators, defining general elevation/depth formation ranges for each, which are essential to calculate the paleo relative sea-level position starting from the measured elevation or depth of the landform.

The last interglacial (LIG), ca. 128-116 ka, is widely considered a process analogue in understanding Earth’s systems in a future warmer climate. In particular, significant effort has been made to better constrain ice sheet contributions to sea level rise through direct field observation of relative sea level (RSL) indicators. In order to extract the RSL, a series of corrections for formational parameters and post-depositional processes need to be applied. Along tropical coastal margins, LIG RSL observations are predominately based on exposed shallow coral reef sequences due to their relatively narrow indicative range and reliable U-series chronological constraints. However, the often-limited sub-stadial temporal preservation of many Pleistocene reef sequences on stable coastlines restrict many reported RSLs to a series of distinct points in within the LIG. This in turn, limits ability to elucidate different commonly reported meter-scale sub-stadial sea level peak patterns and their associated uncertainties. In order to address this shortcoming, lithostratigraphic and geomorphologic traces are often used to place RSLs into a broader context. Unfortunately, this is often subjective, with significant reliance on field observations where missing facies and incomplete sequences can distort interpretations. Stepping back from a conventional approach, in this study we generate a spectrum of synthetic Quaternary subtropical fringing reefs in southwestern Madagascar within the DIONISOS forward stratigraphic model environment. Each reef sequence has been subjected to distinct Greenland and Antarctica melt scenarios produced by a coupled ANICE-SELEN global isostatic adjustment model, matching previously hypothesized LIG sea level curves in the Indo-Pacific Basin. The resulting suite of synthetic reef sequences provides the ability to probabilistically test any number of melt scenarios against the sensitivity of the stratigraphic record. We propose this accessible additional quantitative quality control during the final interpretation phase of establishing emergent reef sequence based LIG RSL indicators can assist in narrowing down the wide uncertainty surrounding inter-stadial ice sheet behaviors.

Using the framework of the World Atlas of Last Interglacial Shorelines (WALIS), we produced a standardized database of Last Interglacial (LIG) sea-level indicators in Southeast Asia after reviewing available studies on relative sea-level (RSL) proxies such as coral reef terraces and tidal notches in the Philippines; Sulawesi; and Sumba, Timor, and Alor regions. In total, we reviewed 43 unique RSL proxies in the region and highlighted sites for future studies. Following this work, we revisited a site in west Luzon, Philippines where LIG coral reef terraces were previously reported. In this paper, we present new geomorphic and stratigraphic data on the fossil coral reef terraces in Pangasinan, west Luzon which adds to the limited sea-level indicators in the region. The low-lying areas of western Pangasinan are underlain by sequences of calcareous sandstone-mudstone with minor pebbly conglomerate and tuffaceous sandstone units belonging to the Sta. Cruz Formation, with tentative age designation of Late Miocene to Early Pliocene. Unconformably overlying the tentatively assigned sandstone unit of Sta. Cruz Formation is the Plio-Pleistocene Bolinao Limestone, the youngest formational unit in the area. Based on previous literature, a sequence of coral reef terraces (possibly LIG) is cut onto the Bolinao Limestone. Rising to about 14 meters above mean sea level (m amsl) along the coast of western Pangasinan are previously dated Holocene coral reef terraces. While additional data is needed to shed more light on the RSL changes in the region, our work proves to be more challenging due to the difficulties of doing field surveys during a global pandemic. Nonetheless, we hope that data from this research will help us further understand the different drivers of past sea-level changes in SE Asia providing necessary geologic baseline data for projections of sea-level change in the future.

During the last decade there has been an increasing demand to improve models of present-day loading processes and glacial-isostatic adjustment (GIA). This is especially important when modelling the GIA process in tectonically active regions like the Pacific Northwest, Patagonia or West Antarctica. All these regions are underlain by zones of low-viscosity mantle. Although one-dimensional earth models may be sufficient to model local-scale uplift within these regions, modeling of the wider-scale deformation patterns requires consideration of three-dimensional viscosity structure that is consistent with other geophysical and laboratory findings. It is this wider-scale modeling that is necessary for earth-system model applications as well as for the validation or reduction of velocity fields determined by geodetic observation networks based on GNSS, for improving satellite gravimetry, and for present-day sea-level change as paleo sea-level reconstructions.There are a number of numerical GIA codes in the community, which can consider lateral variations in viscoelastic earth structure, but a proper benchmark focusing on lateral heterogeneity is missing to date. Accordingly, ambiguity remains when interpreting the modelling results. The numerical codes are based on rather different methods to solve the respective field equations applying, e.g., finite elements, finite volumes, finite differences or spectral elements. Aspects like gravity, compressibility and rheology are dealt with differently. In this regard, the set of experiments to be performed has to be agreed on carefully, and we have to accept that not all structural features can be considered in every code.We present a tentative catalogue of synthetic experiments. These are designed to isolate different aspects of lateral heterogeneity of the Earth's interior and investigate their impact on vertical and horizontal surface displacements, geocenter and polar motion, gravity, sea-level change and stress. The study serves as a follow up of the successful benchmarks of Spada et al. (2011) and Martinec et al. (2018) on 1D earth models and the sea-level equation. The study was initiated by the PALSEA-SERCE Workshop in 2021 (Austermann and Simms, 2022) and benefits from discussions inside different SCAR-INSTANT subcommittees, the IAG Joint Study Group 3.1 “Geodetic, Seismic and Geodynamic Constraints on Glacial Isostatic Adjustment", the IAG Subcommission 3.4 “Cryospheric Deformation" and PALSEA.References:Austermann, J., Simms, A., 2022 (in press). Unraveling the complex relationship between solid Earth deformation and ice sheet change. PAGES Mag., 30(1). doi:10.22498/pages.30.1.14Martinec, Z., Klemann, V., van der Wal, W., Riva, R. E. M., Spada, G., Sun, Y., Melini, D., Kachuck, S. B., Barletta, V., Simon, K., A, G., James, T. S., 2018. A benchmark study of numerical implementations of the sea level equation in GIA modelling. Geophys. J. Int., 215:389-414. doi:10.1093/gji/ggy280Spada, G., Barletta, V. R., Klemann, V., Riva, R. E. M., Martinec, Z., Gasperini, P., Lund, B., Wolf, D., Vermeersen, L. L. A., King, M. A. (2011). A benchmark study for glacial isostatic adjustment codes. Geophys. J. Int., 185:106-132. doi:10.1111/j.1365-246X.2011.04952.x

We present Version 1.0 of the World Atlas of Last Interglacial Shorelines (WALIS), a global database containing samples and sea-level proxies dated to Marine Isotope Stage 5 (~70 to 130 ka). The database was built through manuscripts and associated datasets compiled in a Special Issue of the journal Earth System Science data (https://essd.copernicus.org/articles/special_issue1055.html). We collated the single contributions (archived in Zenodo at this link: https://zenodo.org/communities/walis_database/) into an open-access standalone database. Database documentation is available at this link: https://doi.org/10.5281/zenodo.3961544. Version 1.0 of the database contains 4005 sea-level index points and 4390 dated samples connected with several tables containing relevant metadata (e.g., elevation measurement techniques, sea-level datums, and literature references).

Since the Pleistocene, the Mediterranean is a hot spot for climate change and human migrations, thus offering the opportunity to investigate how human populations have responded to environmental changes and sea level variations. This is the main topic of the SPHeritage Project (MUR grant: FIRS2019_00040, P.I.: M. Pappalardo) that proposes an interdisciplinary approach to investigate the human-environment interaction (in particular sea level variations) over the last 400,000 years using a combination of micro-invasive methods. The Project is re-investigating the well-known archaeological area of ​​the Balzi Rossi (Ventimiglia, at the border between Italy and France), which represents a unique assemblage of archaeological sites dating to the Palaeolithic, distributed in a geomorphological setting rich of markers of past sea level changes. As most of the local archaeological sequences have been extensively investigated at the beginning of the last century and large part of the deposits removed, we will combine the analyses of materials preserved in museums (including strips of sediments) and the remnants still preserved inside many rockshelters and caves of the archaeological complex. Moreover, our geomorphological survey identified new sedimentary sequences preserving information on relative sea level changes. This approach will permit to obtain innovative data submitting small samples to state-of-the-art methods for dating and palaeoenvironmental reconstruction, thus offering the opportunity to better constrain the time and steps of climate change, sea level oscillations, and human settlements. Moreover, data will converge into geoheritage analyses aimed at finding the best practices for promoting and protecting the site. Here, we present an overview of the project and preliminary results from some of the major archaeological sites.

International audience; In the context of industrial-era global change, Mediterranean coastal areas are threatened by relative sea level (RSL) rise. Shifts in the drivers of coastal dynamics are forecasted to trigger changes in the frequency of flooding of low-lying areas, with significant effects on marine-coastal environments, societies, economy and urban systems. Here, we probe coastal floods in the eastern part of the Gulf of Venice (coastal Croatia) to understand the drivers of saltwater intrusions. We reconstructed RSL rise in the north-eastern Adriatic during the Holocene based on 43 RSL index points and analyzed the evolution of coastal submersions on the Istrian Peninsula for the last 5200 years based on inputs of marine components and increases in supratidal scrubs. We produced pollen-based climate reconstructions to analyze the potential effects of air temperature and precipitation changes on submersions. We investigated the response of precipitation and temperature to mid-late Holocene summer/winter insolation forcing and insolation-induced changes in sea surface temperatures (SSTs). We found that during periods of warmer SST, coastal flooding increased markedly. This process seems to have been initiated by warmer atmospheric temperatures that led to increases in summer evaporation, counterbalancing the effects of heavy precipitation during winter. As a result, freshwater flows into coastal areas were reduced and resulted in recurrent inputs of saltwater inland. Our study suggests that present-day climate drivers (increases in SSTs and air temperatures, and decreases in precipitation) will probably favour frequent coastal flooding, a process that will be amplified by RSL rise.

Determining the rates, mechanisms, and geographic variability of relative sea-level (RSL) change following the Last Glacial Maximum (LGM) provides insight into the sensitivity of ice sheets to climate change, the response of the solid Earth and gravity field to ice-mass redistribution, and constrains statistical and physical models used to project future sea-level rise. To do so in a scientifically robust way requires standardized datasets that enable broad spatial comparisons that minimize bias. As part of a larger goal to develop a unified, spatially-comprehensive post-LGM global RSL database, in this special issue we provide a standardized global synthesis of regional RSL data that resulted from the first ‘Geographic variability of HOLocene relative SEA level (HOLSEA)’ meetings in Mt Hood, Oregon (2016) and St Lucia, South Africa (2017). The HOLSEA meetings brought together sea-level researchers to agree upon a consistent protocol to standardize, interpret, and incorporate realistic uncertainties of RSL data. This special issue provides RSL data from ten geographical regions including new databases from Atlantic Europe and the Russian Arctic and revised/expanded databases from Atlantic Canada, the British Isles, the Netherlands, the western Mediterranean, the Adriatic, Israel, Peninsular Malaysia, Southeast Asia, and the Indian Ocean. In total, the database derived from this special issue includes 5634 (5290 validated) index (n = 3202) and limiting points (n = 2088) that span from ∼20,000 years ago to present. Progress in improving the standardization of sea-level databases has also been accompanied by advancements in statistical and analytical methods used to infer spatial patterns and rates of RSL change from geological data that have a spatially and temporally sparse distribution and geochronological and elevational uncertainties. This special issue marks the inception of a unified, spatially-comprehensive post-LGM global RSL database. Ministry of Education (MOE) National Research Foundation (NRF) Published version NSK and BPH were funded by the Singapore Ministry of Edu- cation Academic Research Fund MOE2018-T2-1-030, the National Research Foundation Singapore, and the Singapore Ministry of Education, under the Research Centres of Excellence initiative. This work is Earth Observatory of Singapore contribution 252. AR ac- knowledges the Institutional Strategy of the University of Bremen, funded by the German Excellence Initiative (ABPZuK-03/2014) and the SEASCHANGE (RO-5245/1-1) from the Deutsche For- schungsgemeinschaft (DFG) as part of the Special Priority Program (SPP)-1889 00Regional Sea Level Change and Society”. TET ac- knowledges support from the Paleo Perspectives on Climate Change program of the US National Science Foundation (OCE- 1502588). We thank Maren Bender for her help in preparing the comparison of data from Southeast Asia. This article is a contribu- tion to PALSEA (Palaeo-Constraints on Sea-Level Rise) and Inter- national Geoscience Program (IGCP) Project 639, “Sea Level Change from Minutes to Millennia”.

Quaternary interglacials provide key observations of the Earth system's responses to orbital and greenhouse gas forcing. They also inform on the capabilities of Earth system models, used for projecting the polar ice-sheet and sea-level responses to a regional warmth comparable to that expected by 2100 C.E. However, a number of uncertainties remain regarding the processes and feedbacks linking climate, ice-sheet and sea-level changes during past warm intervals. Here, we delineate the major research questions that need to be resolved and future research directions that should be taken by the paleoclimate, sea-level and ice-sheet research communities in order to increase confidence in the use of past interglacial climate, ice-sheet and sea-level reconstructions to constrain future predictions. These questions were formulated during a joint workshop held by the PAGES-INQUA PALSEA (PALeo constraints on SEA level rise) and the PAGES-PMIP QUIGS (QUaternary InterGlacialS) Working Groups in September 2018.

Sea-level rise is predicted to cause major damage to tropical coastlines. While coral reefs can act as natural barriers for ocean waves, their protection hinges on the ability of scleractinian corals to produce enough calcium carbonate (CaCO3 ) to keep up with rising sea levels. As a consequence of intensifying disturbances, coral communities are changing rapidly, potentially reducing community-level CaCO3 production. By combining colony-level physiology and long-term monitoring data, we show that reefs recovering from major disturbances can produce 40% more CaCO3 than currently estimated due to the disproportionate contribution of juvenile corals. However, the buffering effect of highly productive juvenile corals is compromised by recruitment failures, which have been more frequently observed after large-scale, repeated bleaching events. While the size structure of corals can bolster a critical ecological function on reefs, climate change impacts on recruitment may undermine this buffering effect, thus further compromising the persistence of reefs and their provision of important ecosystem services.

In this paper, we describe a sea-level database compiled using published last interglacial, Marine Isotopic Stage 5 (MIS 5), geological sea-level proxies within East Africa and the Western Indian Ocean (EAWIO). Encompassing vast tropical coastlines and coralline islands, this region has many occurrences of well-preserved last interglacial stratigraphies. Most notably, islands almost entirely composed of Pleistocene reefs (such as Aldabra, the Seychelles) have provided reliable paleo relative sea-level indicators and well-preserved samples for U-series chronology. Other sea-level proxies include uplifted marine terraces in the north of Somalia and Pleistocene eolian deposits notched by the MIS 5 sea level in Mozambique to tidal notches in luminescence-limited eolian deposits in Mozambique. Our database has been compiled using the World Atlas of Last Interglacial Shorelines (WALIS) interface and contains 58 sea-level indicators and 2 terrestrial-limiting data points. The open-access database is available at https://doi.org/10.5281/zenodo.4302244 (Version 1.03; Boyden et al., 2020).

An open access database containing raw data of Last Interglacial sea-level proxies for the Western Mediterranean has been compiled by reviewing hundreds of original published papers in accordance with the WALIS template (https://warmcoasts.eu/world-atlas.html). WALIS allows collecting both the relative sea-level (RSL) indicators and ages data in a standardized format. Ca. 360 sea-level index points for the coasts of Spain, France, Italy, Albania, Algeria and Morocco have been included in the database. The sea-level index points of the database are related to ca. 350 samples dated by a wide range of dating techniques, e.g., U-series, Amino Acid Racemization, Luminescence (Tl/OSL) and Electron Spin Resonance methods or chronostratigraphically correlated to marine deposits bearing P. latus and “Senegalese fauna”. In fact, for some areas of the Mediterranean Sea, the “Senegalese fauna” is indicative of the MIS 5e.Among the eleven types of sea-level indicators of our database, the majority of them are represented by marine terraces, beach deposits (or beachrocks), and tidal notches. Whenever the relationship between the RSL indicators and the former sea-level could not be quantified, such indicators have been considered as marine or terrestrial limiting points. An indirect age of the tidal notches has been provided by correlation with the nearby dated deposit. In the case no precise elevation information has been reported by the Authors, the elevation error of RSL datapoints has been reassessed in the 20 % of the elevation value, more a 5% if the sea-level datum was lacking in the scientific papers we have reviewed. Overall, the quality of each RSL datapoints and the associated age have been ranked in a 0 to 5 scale score according to Rovere et al., (2020). ReferencesRovere, A., Ryan, D., Murray-Wallace, C., Simms, A., Vacchi, M., Dutton, A., Gowan, E., 2020. Descriptions of database fields for the World Atlas of Last Interglacial Shorelines (WALIS) (Version 1,0). Zenodo. https://doi.org/http://doi.org/10.5281/zenodo.3961544

The standardization of geological data, and their compilation into geodatabases, is essential to allow more coherent regional and global analyses. In sea-level studies, the compilation of databases containing details on geological paleo sea-level proxies has been the subject of decades of work. This was largely spearheaded by the community working on Holocene timescales. While several attempts were also made to compile data from older interglacials, a truly comprehensive approach was missing. Here, we present the ongoing efforts directed to create the World Atlas of Last Interglacial Shorelines (WALIS), a project spearheaded by the PALSEA (PAGES/INQUA) community and funded by the European Research Council (ERC StG 802414). The project aims at building a sea-level database centered on the Last Interglacial (Marine Isotope Stage 5e, 125 ka), a period of time considered as an "imperfect analog" for a future warmer climate. The database is composed of 17 tables embedded into a mySQL framework with a total of more than 500 single fields to describe several properties related to paleo sea-level proxies, dated samples and metadata. In this presentation, we will show the first results of the global compilation, which includes nearly 2000 data points and will discuss its relevance in answering some of the most pressing questions related to sea-level changes in past warmer worlds.

South-eastern Sicily is among the most seismically active areas of the central Mediterranean. As such, it is marked by a high level of crustal seismicity producing major earthquakes (up to Mw ∼7), and consequent several earthquake-generated tsunami, which have affected the Ionian coast of South-eastern Sicily in historical times. These tsunami events left geomorphic imprints such as large boulders or high-energy deposits along the Sicily coasts. In Ognina, a small town located 20 km south of Siracusa, high-energy deposits were correlated with three tsunami events that struck this coast on 21 July 365 Common Era (CE), 4 February 1169 CE, and 11 January 1693 CE. The deposits are detected in the inner part of a narrow channel, that is thought to have funnelled the tsunami flow energy. In this work, numerical models have been performed to simulate the tsunami impacts, considering the most probable tsunamogenic sources described in literature and integrating them with the past sea-level positions. To this end, we used Delft Dashboard, Delft 3d-FLOW and XBeach. A reconstruction of the past topography of Ognina coast was performed through geological and historical information, in order to model the tsunami wave propagation in the ancient landscape. Geological evidence with model results, under different scenarios, allow us to benchmark fault location and displacement scenarios. Modelling results indicate that the 1693 tsunami event was stronger than others impacting the Ognina area, determining significant inland flooding in the narrow channel. Moreover, simulations show that the most probable tsunamogenic sources of 1693 and 1169 tsunami events could be attributed to Western Fault dislocations occurred off-shore of Ognina area, rather than the other tsunamogenic sources described in literature, located off-shore of Catania and Siracusa. Modelling of 365 AD event shows a long period for the tsunami wave that determined the sedimentation on the lower units in the outcrop. For each of the three tsunami events, models of high-energy deposition match with position and thickness of high-energy layers detected in the field. The results of this study show how a combined approach between geological evidence and tsunami modelling could be a suitable tool for the attribution of tsunami deposits connected to specific tsunamogenic sources. Keyword: tsunami; earthquake; faults; flooding; sea-level

The Last Interglacial (LIG), as well as other warmer periods in the Earth’s geologic history, provides an analogue for predicted warming conditions in the near future. Analysis of sea-level indicators during this period is important in constraining regional drivers of relative sea-level change (RSL) and in modeling future trajectories of sea-level rise. In southeast Asia, several studies have been done to examine LIG sea-level indicators such as coral reef terraces and tidal notches. A synthesis of the state-of-the-art of the LIG RSL indicators in the region, meanwhile, has yet to be done. We reviewed over 50 published works on the LIG RSL indicators in southeast Asia and used the framework of the World Atlas of Last Interglacial Shorelines (WALIS) in building a standardized database of previously published LIG RSL indicators in the region. In total, we identified 38 unique RSL indicators and inserted almost 140 ages in the database. Available data from Indonesia, the Philippines, and East Timor points to variable elevation of sea-level indicators during the LIG highlighting the complex tectonic setting of this region. Variable uplift rates (from as low as 0.02 to as high as 1.1 m/ka) were reported in the study areas echoing various collision and subduction processes influencing these sites. Although several age constraints and elevation measurements have been provided by these studies, more data is still needed to shed more light on the RSL changes in the region. With this effort under the WALIS framework, we hope to identify gaps in the LIG RSL indicators literature in SE Asia and recognize potential areas that can be visited for future work. We also hope that this initiative will help us further understand the different drivers of past sea-level changes in SE Asia and will provide inputs for projections of sea-level change in the future.

Bivalve and gastropod shell beds deposited during the Early Pliocene (4.69-5.23 Ma) occur in uplifted outcrops (36 – 180 m above sea level) along the east coast of Patagonian Argentina. These rock units provide a record of sea level during a geologic period when atmospheric CO2 and temperatures were higher than today. As such, reconstructing the elevation of global mean sea level (GMSL) during this time allows us to better understand how sensitive ice sheets are to increased past and future warming. However, reconstructing GMSL from local sea level indicators is hindered by effects such as mantle dynamic topography and glacial isostatic adjustment (GIA) that cause local sea level to deviate from the global mean. Here we use geodynamic modeling to better understand this complex dynamic setting and quantify the amount of uplift along this coastline. Despite being located on a relatively stable passive margin, significant variations in the elevation of the paleo shoreline indicators imply that the underlying convecting mantle is deforming the coastline. In particular, the subduction of the Chile Rise beginning ~18 Ma beneath Patagonia has generated a slab window underneath this region through which hot asthenosphere ascends. However, the former slab is still present deeper in the mantle, which causes a complex interplay between the downwelling slab and the upwelling asthenosphere. To quantify the effects of dynamic topography change since the Pliocene, we run 3D mantle convection simulations using the code ASPECT. We initialize our global model with a composite temperature structure derived from recent tomographic studies and a calibrated parameterization of upper mantle anelasticity. Independent estimates of pressure and temperature from thermobarometric calculations of proximal Pali-Aike xenoliths agree with the thermal structure of the tomography-based Earth model. We back-advect temperature perturbations and extract the resulting change in dynamic topography. Pairing GIA models and a suite of convection simulations in which we vary the viscosity and buoyancy structure with the observed differential paleo shoreline elevations allows us to forward model the most likely scenario for uplift along this coast.

Reconstructing the topography of shallow underwater environments using Structure-from-Motion – Multi View Stereo (SfM-MVS) techniques applied to aerial imagery from Unmanned Aerial Vehicles (UAVs) is a challenging problem, as it involves non-linear distortions caused by water refraction. This study presents an experiment with aerial photographs collected with a consumer-grade UAV on the shallow-water reef of Fuvahmulah, the Maldives. Under conditions of rising tide, we surveyed the same portion of the reef in ten successive flights. For each flight, we used SfM-MVS to reconstruct the Digital Elevation Model (DEM) of the reef, and used the flight at low tide (where the reef is almost entirely dry) to compare the performance of DEM reconstruction under higher water levels. Our results show that differences with the reference DEM increase with increasing depth, but are substantially larger if no underwater ground control points are taken into account in the processing. Correcting our imagery with algorithms that account for refraction did not improve the overall accuracy of reconstruction. We conclude that reconstructing shallow-water reefs (less than 1 m depth) with consumer-grade SfM-MVSs and SfM-MVS is possible, but its precision is limited and strongly correlated with water depth. The best results are achieved when ground control points were placed both above and underwater and no refraction correction is used in our processing.

Sea-level rise represents a severe hazard for populations living within low-elevation coastal zones and is already largely affecting coastal communities worldwide. As sea level continues to rise following unabated greenhouse gas emissions, the exposure of coastal communities to inundation and erosion will increase exponentially. These impacts will be further magnified under extreme storm conditions. In this paper, we focus on one of the most valuable coastal real estate markets globally (Palm Beach, FL). We use XBeach, an open-source hydro and morphodynamic model, to assess the impact of a major tropical cyclone (Hurricane Matthew, 2016) under three different sea-level scenarios. The first scenario (modern sea level) serves as a baseline against which other model runs are evaluated. The other two runs use different 2100 sea-level projections, localized to the study site: (i) IPCC RCP 8.5 (0.83 m by 2100) and (ii) same as (i), but including enhanced Antarctic ice loss (1.62 m by 2100). Our results show that the effective doubling of future sea level under heightened Antarctic ice loss amplifies flow velocity and wave height, leading to a 46% increase in eroded beach volume and the overtopping of coastal protection structures. This further exacerbates the vulnerability of coastal properties on the island, leading to significant increases in parcel inundation.

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.

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).