Interglacial Sea Level Research Articles

Contrasting the Penultimate Glacial Maximum and the Last Glacial Maximum (140 and 21 ka) using coupled climate–ice sheet modelling

Abstract. The configuration of the Northern Hemisphere ice sheets during the Penultimate Glacial Maximum differed to the Last Glacial Maximum. However, the reasons for this are not yet fully understood. These differences likely contributed to the varied deglaciation pathways experienced following the glacial maxima and may have had consequences for the interglacial sea level rise. To understand the differences between the North American Ice Sheet at the Last and Penultimate glacial maxima (21 and 140 ka), we perform two perturbed-physics ensembles of 62 simulations using a coupled atmosphere–ice sheet model, FAMOUS-ice, with prescribed surface ocean conditions, in which the North American and Greenland ice sheets are dynamically simulated with the Glimmer ice sheet model. We apply an implausibility metric to find ensemble members that match reconstructed ice extent and volumes at the Last and Penultimate glacial maxima. We use a resulting set of “plausible” parameters to perform sensitivity experiments to decompose the role of climate forcings (orbit, greenhouse gases) and initial conditions on the final ice sheet configurations. This confirms that the initial ice sheet conditions used in the model are extremely important in determining the difference in final ice volumes between both periods due to the large effect of the ice–albedo feedback. In contrast to evidence of a smaller Penultimate North American Ice Sheet, our results show that the climate boundary conditions at these glacial maxima, if considered in isolation, imply a larger Penultimate Glacial Maximum North American Ice Sheet than at the Last Glacial Maximum by around 6 m sea level equivalent. This supports the notion that the growth of the ice sheet prior to the glacial maxima is key in explaining the differences in North American ice volume.

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.

Quantifying the uncertainty in the Eurasian ice-sheet geometry at the Penultimate Glacial Maximum (Marine Isotope Stage 6)

Abstract. The North Sea Last Interglacial sea level is sensitive to the fingerprint of mass loss from polar ice sheets. However, the signal is complicated by the influence of glacial isostatic adjustment driven by Penultimate Glacial Period ice-sheet changes, and yet these ice-sheet geometries remain significantly uncertain. Here, we produce new reconstructions of the Eurasian ice sheet during the Penultimate Glacial Maximum (PGM) by employing large ensemble experiments from a simple ice-sheet model that depends solely on basal shear stress, ice extent, and topography. To explore the range of uncertainty in possible ice geometries, we use a parameterised shear-stress map as input that has been developed to incorporate bedrock characteristics and the influence of ice-sheet basal processes. We perform Bayesian uncertainty quantification, utilising Gaussian process emulation, to calibrate against global ice-sheet reconstructions of the Last Deglaciation and rule out combinations of input parameters that produce unrealistic ice sheets. The refined parameter space is then applied to the PGM to create an ensemble of constrained 3D Eurasian ice-sheet geometries. Our reconstructed PGM Eurasian ice-sheet volume is 48±8 m sea-level equivalent (SLE). We find that the Barents–Kara Sea region displays both the largest mean volume and volume uncertainty of 24±8 m SLE while the British–Irish sector volume of 1.7±0.2 m SLE is the smallest. Our new workflow may be applied to other locations and periods where ice-sheet histories have limited empirical data.

Refining patterns of melt with forward stratigraphic models of stable Pleistocene coastlines

Abstract. 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 efforts have been made to better constrain ice sheet contributions to the 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, which also provide 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 of unraveling 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 stratigraphic model (from Beicip-Franlab) 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 modeling may provide further quantitative control when interpreting Last Interglacial reef sequences.

Last interglacial global mean sea level from high-precision U-series ages of Bahamian fossil coral reefs

Accurate characterization of Last Interglacial (MIS 5e; ∼129–116 ka) sea level is important for understanding ice sheet sensitivity to climate change, with implications for predicting future sea-level rise. Here we present a record of MIS 5e sea level based on high-precision U-series ages of 23 corals with precise elevation measurements from reefs around Crooked Island, Long Cay, Long Island, and Eleuthera, The Bahamas. Rigorous screening criteria identified the most pristine samples, and nearly all samples show a narrow δ234Uinitial range between 143.8 and 151.3‰. We infer global mean sea level (GMSL) from these local observations by correcting them for glacial isostatic adjustment (GIA) and long-term subsidence. For GIA, we consider a range of ice histories and Earth viscosity structures. We identify, via Bayesian inference, the range of isostatic and GMSL histories that are consistent with MIS 5e observations across The Bahamas. When applying an open-system correction to our ages, we find that MIS 5e GMSL likely peaked higher than 1 m, but very unlikely exceeded 2.7 m. Our posterior GMSL is lower than previous estimates, but consistent with recent results of modeling and observations. Additionally, sea level observations at other locations (Seychelles, Western Australia, Yucatan) are only slightly above/within the 95% range of predicted local sea level, i.e., GIA plus GMSL, for our open-system/closed-system results. Our relatively constant MIS 5e GMSL indicates that Greenland and Antarctica melted beyond their present extents and, given the insolation forcing, that their contributions to GMSL were likely out-of-phase. These results indicate that the ice sheets may be very sensitive to regional temperature, which has important implications for their combined impact on global sea levels at a time when greenhouse gases increases are causing simultaneous warming at both poles.

Record of the last interglacial sea level highstand based on new coastal deposits in the Cantabrian margin (Northern Iberian Peninsula)

Quaternary sea level fluctuations have been responsible for significant remodelling of coastal landscapes, generating distinctive sea level indicators at different positions and altitudes that are still visible today. Raised coastal deposits are among the characteristic features of significant value in reconstructing the tectonic and sea-level history of the areas in which they are preserved. To date, few remnants of Quaternary coastal deposits have been described and directly dated along the coast of the Cantabrian margin (on the north of the Iberian Peninsula). In this work, two new deposits located in the towns of Castro-Urdiales and Mendexa were analysed and dated using U/Th geochronology and Optically Stimulated Luminescence. Considering the ages and altitudes of these deposits, it seems most likely that they correlate to the Marine Isotope Stage (MIS) 5e. Furthermore, the data acquired suggest that most of the Cantabrian margin underwent a similar tectonic evolution, with no significant tectonic uplift, since at least MIS 5e.

Constraining the contribution of the Antarctic Ice Sheet to Last Interglacial sea level

Polar temperatures during the Last Interglacial [LIG; ~129 to 116 thousand years (ka)] were warmer than today, making this time period an important testing ground to better understand how ice sheets respond to warming. However, it remains debated how much and when the Antarctic and Greenland ice sheets changed during this period. Here, we present a combination of new and existing absolutely dated LIG sea-level observations from Britain, France, and Denmark. Because of glacial isostatic adjustment (GIA), the LIG Greenland ice melt contribution to sea-level change in this region is small, which allows us to constrain Antarctic ice change. We find that the Antarctic contribution to LIG global mean sea level peaked early in the interglacial (before 126 ka), with a maximum contribution of 5.7 m (50th percentile, 3.6 to 8.7 m central 68% probability) before declining. Our results support an asynchronous melt history over the LIG, with an early Antarctic contribution followed by later Greenland Ice Sheet mass loss.

New datings and elevations of a fossil reef in Lembetabe, southwest Madagascar: eustatic and tectonic implications

The study of geological sea-level proxies formed during previous interglacials is a common approach to assess how global sea level will evolve under warmer climate conditions. Over the last decades, technical advancements in both survey and geochronology have allowed improving our knowledge of past sea-level highstands. This is of prime importance to refine our understanding of processes contributing to sea-level changes, and ultimately to improve both local and global sea-level projections. Last Interglacial sea-level proxies in the Western Indian Ocean (and more specifically the island nation of Madagascar), have been less investigated than in other intertropical oceans over the last decades. As a result, paleo sea-level data in this region are less abundant and less precise than elsewhere. Here, we report the results of two field campaigns aimed at studying the site of Lembetabe, southwest Madagascar, where a fossil reef was first described by the researcher René Battistini more than 50 years ago. We estimate paleo relative sea level history in space and time from 15 new U-series ages from a fossil reef platform mapped with differential GNSS and drone photogrammetry. Our data suggest that, between 129 ka and 115 ka, paleo relative sea level at this location was about 3.4 ± 1.4 m above modern. Once corrected for glacial isostatic adjustment, we find that paleo global mean sea level did not exceed 3 m above modern. Only slight crustal subsidence would reconcile the peak Last Interglacial sea level measured at Lembetabe with the 5–10 m range reported in the literature.

Robbins Island: The index site for regional Last Interglacial sea level, wave climate and the subtropical ridge around Bass Strait, Australia

A unique index-record of Last Interglacial (Marine Isotope Stage 5e MIS5e) relative sea level (RSL) and wave climate history in South-east Australia is presented from Robbins Island, in western Bass Strait. This is applied to interpret the wider MIS5e coastal evidence around Bass Strait. At Robbins Island, the combination of low wave and wind energy, a tide-modified regime and a sand supply resulted in the shoreline progradation throughout MIS5e. This preserved a time-series of paleo-sea level across a 7 km wide strandplain (Remarkable Banks). After a highstand, MIS5e RSL attained a stillstand of +5.75 ± 0.5 m above modern mean sea level during 126 to ∼119 ka BP. The MIS5e RSL interpretation is underpinned by modern analogues and hydrodynamic modelling of waves, tides and currents. A high resolution LiDAR Digital Elevation Model (DEM) supported by morpho-sedimentary studies, ground-penetration radar (GPR) surveys and a geochronology based upon Optically Stimulated Luminescence (OSL) methods were used to define the proxy RSL record. The observed RSL history was compared to modelled RSL history that accounted for the theoretical fall in RSL (regression) throughout MIS5e, due to the Glacio-Isostatic Adjustment (GIA) forcing.

Relict inland mangrove ecosystem reveals Last Interglacial sea levels

Climatic oscillations during the Pleistocene played a major role in shaping the spatial distribution and demographic dynamics of Earth's biota, including our own species. The Last Interglacial (LIG) or Eemian Period (ca. 130 to 115 thousand years B.P.) was particularly influential because this period of peak warmth led to the retreat of all ice sheets with concomitant changes in global sea level. The impact of these strong environmental changes on the spatial distribution of marine and terrestrial ecosystems was severe as revealed by fossil data and paleogeographic modeling. Here, we report the occurrence of an extant, inland mangrove ecosystem and demonstrate that it is a relict of the LIG. This ecosystem is currently confined to the banks of the freshwater San Pedro Mártir River in the interior of the Mexico-Guatemala El Petén rainforests, 170 km away from the nearest ocean coast but showing the plant composition and physiognomy typical of a coastal lagoon ecosystem. Integrating genomic, geologic, and floristic data with sea level modeling, we present evidence that this inland ecosystem reached its current location during the LIG and has persisted there in isolation ever since the oceans receded during the Wisconsin glaciation. Our study provides a snapshot of the Pleistocene peak warmth and reveals biotic evidence that sea levels substantially influenced landscapes and species ranges in the tropics during this period.

The effect of lateral variations in Earth structure on Last Interglacial sea level

SUMMARYIt is generally agreed that the Last Interglacial (LIG; ∼130–115 ka) was a time when global average temperatures and global mean sea level were higher than they are today. However, the exact timing, magnitude and spatial pattern of ice melt is much debated. One difficulty in extracting past global mean sea level from local observations is that their elevations need to be corrected for glacial isostatic adjustment (GIA), which requires knowledge of Earth’s internal viscoelastic structure. While this structure is generally assumed to be radially symmetric, evidence from seismology, geodynamics and mineral physics indicates that large lateral variations in viscosity exist within the mantle. In this study, we construct a new model of Earth’s internal structure by converting shear wave speed into viscosity using parametrizations from mineral physics experiments and geodynamic constraints on Earth’s thermal structure. We use this 3-D Earth structure, which includes both variations in lithospheric thickness and lateral variations in viscosity, to calculate the first 3-D GIA prediction for LIG sea level. We find that the difference between predictions with and without lateral Earth structure can be metres to 10s of metres in the near field of former ice sheets, and up to a few metres in their far field. We demonstrate how forebulge dynamics and continental levering are affected by laterally varying Earth structure, with a particular focus on those sites with prominent LIG sea level records. Results from four 3-D GIA calculations show that accounting for lateral structure can act to increase local sea level by up to ∼1.5 m at the Seychelles and minimally decrease it in Western Australia. We acknowledge that this result is only based on a few simulations, but if robust, this shift brings estimates of global mean sea level from these two sites into closer agreement with each other. We further demonstrate that simulations with a suitable radial viscosity profile can be used to locally approximate the 3-D GIA result, but that these radial profiles cannot be found by simply averaging viscosity below the sea level indicator site.

The marine terraces of Santa Cruz Island, California: Implications for glacial isostatic adjustment models of last-interglacial sea-level history

Glacial isostatic adjustment (GIA) models hypothesize that along coastal California, last interglacial (LIG, broadly from ~130 to ~115 ka) sea level could have been as high as +11 m to +13 m, relative to present, substantially higher than the commonly estimated elevation of +6 m. Areas with low uplift rates can test whether such models are valid. Marine terraces on Santa Cruz Island have previously been reported to occur at low (<10 m) elevations, but ages of many such localities are not known. Using lidar imagery as a base, marine terraces on Santa Cruz Island were newly mapped, elevations were measured, fossils were collected for U-series dating (corals), strontium isotope compositions and amino acid geochronology (mollusks), and paleozoogeography (all taxa). Sr isotope compositions of mollusks from the highest of three marine terraces give ages of ~2.5 Ma to 1.9 Ma, along with Pliocene ages, from shells interpreted to be reworked. U-series ages of corals from the western part of the island indicate that low-elevation terraces north of the Santa Cruz Island fault correlate to the LIG. Where corals are lacking, amino acid ratios and faunal aspects support terrace correlation to the LIG high stand of sea. Elevations of most terrace localities north of the east-west trending Santa Cruz Island fault, in both the western and eastern parts of the island, range from 5.75 m to 8 m above sea level, well below the modeled paleo-sea-level range. Subsidence is ruled out as a mechanism for explaining the lower-than-modeled elevations, because higher-elevation terraces are present along much of the Santa Cruz Island coast north of the fault, indicating long-term tectonic uplift. The low elevations of the LIG terrace fragments are, however, consistent with a low rate of uplift derived from the higher, ~2.5–1.9 Ma terrace. A number of other localities on the Pacific Coast, also dated to the LIG, have marine terrace elevations below the modeled level. GIA models may have overestimated last interglacial sea level by a substantial amount and need to be revised if used for forecasts for future sea-level rise.

Last interglacial sea levels within the Gulf of Mexico and northwestern Caribbean Sea

Abstract. During the last interglacial (LIG) the volume of additional water in the world's oceans was large enough to raise global sea levels about 6–9 m higher than present levels. However, LIG sea levels vary regionally and those regional differences hold clues about the past distribution of ice sheets and local rates of subsidence and tectonic uplift. In this study, I used a standardized database template to review and summarize the existing constraints on LIG sea levels across the northern Gulf of Mexico and Caribbean shoreline of the Yucatán Peninsula. In total, I extracted 32 sea-level indicators including the insertion of 16 U-series ages on corals, 1 electron spin resonance age, 2 amino acid racemization ages, and 26 luminescence ages. Most dated sea-level indicators for the northern Gulf of Mexico are based on optically stimulated luminescence (OSL) ages of beach deposits of a mappable LIG shoreline. This shoreline extends from the Florida Panhandle through south Texas but is buried or removed by the Mississippi River across most of Louisiana. A similar feature is observed in satellite images south of the Rio Grande within the Mexican portions of the Gulf of Mexico but has yet to be dated. Elevations measured on portions of this feature close to the modern coast point to sea levels less than 1 m to ∼5 m higher than present for much of the northern Gulf of Mexico. However, a few, albeit undated, portions of the same shoreline located at more inland locations point to sea levels up to +7.2 m, attesting to up to 7 m of differential subsidence between the inland and coastal sites. Across the Yucatán Peninsula, U-series dating of corals has provided the main index points for LIG sea levels. Other carbonate coastal features such as beach ridges and eolianites have also been described but rely on corals for their dating. The maximum elevation of the LIG coral-based relative sea-level (RSL) estimates decrease from around +6 m across the Caribbean shoreline of the Yucatán Peninsula near Cancún, Mexico, to as low as −6 m to the south beneath the southern atolls of Belize, although discussion continues as to the validity of the ages for these southern corals. If these lower-elevation corals are LIG in age, their below-present elevations may be a result of vertical motion along faults dipping into the Cayman Trough. South of Belize only one purported LIG coral has been dated on the Isla de Roatán off the coast of Honduras at a likely tectonically uplifted elevation of 37.2 m. Thus the elevation of LIG sea levels within the inland siliciclastic shorelines of Guatemala and Honduras as well as the southwestern Gulf of Mexico remains poorly constrained and a potential venue for future research. The database described in this paper is available open access in spreadsheet format as Simms (2020), at this link: https://doi.org/10.5281/zenodo.4556163.

Late Quaternary dynamics of the Lambert Glacier-Amery Ice Shelf system, East Antarctica

The Lambert Glacier-Amery Ice Shelf system (LGAISS) is the largest outlet glacier system in East Antarctica but its response to past climate variability is poorly constrained. In this study, we explore its dynamics over the last ∼520 thousand years using new high-resolution sedimentary records retrieved off Prydz Bay. Episodic occurrences of iceberg rafted debris indicate a dynamic ice sheet throughout this interval, while changes in clay mineral compositions provide detailed evidence on the waxing and waning of the LGAISS. Our data indicate that advance and retreat of the LGAISS was sensitive to oceanic forcing, but also responded to local summer insolation when insolation peaks combined with interglacial sea level high stands. Subglacial bed topography may have played an additional role in modulating the LGAISS behaviour, which could explain sectoral differences in the response of the East Antarctic Ice Sheet to climate change. Overall, our records indicate that the LGAISS advanced more extensively during previous late Quaternary glacial periods than during Marine Isotope Stages 2 and 4. Furthermore, the LGAISS retreated more significantly than present during Marine Isotope Stage 13, and only moderately during Marine Isotope Stage 7, which suggests that the duration of warm climatic states could be a key factor affecting the dynamics of the East Antarctic Ice Sheet.

Early to Middle Holocene Estuarine Shellfish Collecting on the Islands and Mainland Coast of the Santa Barbara Channel, California, USA

Terminal Pleistocene to Middle Holocene sea level rise resulted in a number of changes to coastal ecosystems around the world, providing new challenges and opportunities for coastal peoples. In California, glacial to interglacial sea level rise resulted in some reductions in rocky shore kelp forests, but it also resulted in the formation of estuaries. Estuaries were important for terminal Pleistocene peoples in the Santa Barbara Channel region (SBC), a pattern that persisted through the Early to Middle Holocene, and sometimes later. While paleoestuaries appear to have been relatively common along the SBC mainland coast, they were rare to absent on the Channel Islands. The Abalone Rocks Paleoestuary on Santa Rosa Island is the only well documented island estuary. However, questions remain about the size and productivity of this estuary and its importance for human subsistence and settlement relative to the more extensive mainland estuaries. Faunal data from two previously unreported site components and synthesis of shellfish data from other Abalone Rocks sites and similarly aged sites near mainland estuaries illustrate the importance of SBC mainland versus island estuaries. Estuarine shellfish were considerably more abundant at most Early and Middle Holocene mainland sites, with the Abalone Rocks Paleoestuary largely supplementary to rocky shore habitats. At island estuary sites, taxonomic richness was fairly consistent during the Early to Middle Holocene, although diversity and evenness decline slightly through time, with estuarine shellfish largely disappearing from island assemblages prior to 5000 years ago. These data demonstrate the power of archaeological research to evaluate the relationships between past environmental change and human behavior.

Towards assessing the influence of sediment loading on Last Interglacial sea level

SUMMARY Locally, the elevation of last interglacial (LIG; ∼122 ka) sea level markers is modulated by processes of vertical displacement, such as tectonic uplift or glacial isostatic adjustment, and these processes must be accounted for in deriving estimates of global ice volumes from geological sea level records. The impact of sediment loading on LIG sea level markers is generally not accounted for in these corrections, as it is assumed that the impact is negligible except in extremely high depositional settings, such as the world's largest river deltas. Here we perform a generalized test to assess the extent to which sediment loading may impact global variability in the present-day elevation of LIG sea level markers. We numerically simulate river sediment deposition using a diffusive model that incorporates a migrating shoreline to construct a global history of sedimentation over the last glacial cycle. We then calculate sea level changes due to this sediment loading using a gravitationally self-consistent model of glacial isostatic adjustment, and compare these predictions to a global compilation of LIG sea level data. We perform a statistical analysis, which accounts for spatial autocorrelation, across a global compilation of 1287 LIG sea level markers. Though limited by uncertainties in the LIG sea level database and the precise history of river deposition, this analysis suggests there is not a statistically significant global signal of sediment loading in LIG sea level markers. Nevertheless, at sites where LIG sea level markers have been measured, local sea level predicted using our simulated sediment loading history is perturbed up to 16 m. More generally, these predictions establish the relative sensitivity of different regions to sediment loading. Finally, we consider the implications of our results for estimates of tectonic uplift rates derived from LIG marine terraces; we predict that sediment loading causes 5–10 m of subsidence over the last glacial cycle at specific locations along active margin regions such as California and Barbados, where deriving long-term tectonic uplift rates from LIG shorelines is a common practice.

Nonlinear response of the Antarctic Ice Sheet to late Quaternary sea level and climate forcing

Abstract. Antarctic ice volume has varied substantially during the late Quaternary, with reconstructions suggesting a glacial ice sheet extending to the continental shelf break and interglacial sea level highstands of several meters. Throughout this period, changes in the Antarctic Ice Sheet were driven by changes in atmospheric and oceanic conditions and global sea level; yet, so far modeling studies have not addressed which of these environmental forcings dominate and how they interact in the dynamical ice sheet response. Here, we force an Antarctic Ice Sheet model with global sea level reconstructions and transient, spatially explicit boundary conditions from a 408 ka climate model simulation, not only in concert with each other but, for the first time, also separately. We find that together these forcings drive glacial–interglacial ice volume changes of 12–14 ms.l.e., in line with reconstructions and previous modeling studies. None of the individual drivers – atmospheric temperature and precipitation, ocean temperatures, or sea level – single-handedly explains the full ice sheet response. In fact, the sum of the individual ice volume changes amounts to less than half of the full ice volume response, indicating the existence of strong nonlinearities and forcing synergy. Both sea level and atmospheric forcing are necessary to create full glacial ice sheet growth, whereas the contribution of ocean melt changes is found to be more a function of ice sheet geometry than climatic change. Our results highlight the importance of accurately representing the relative timing of forcings of past ice sheet simulations and underscore the need for developing coupled climate–ice sheet modeling frameworks that properly capture key feedbacks.

What do we know about Last Interglacial sea level?

What do we know about Last Interglacial sea level?

The archaeological ensemble from Campoverde (Agro Pontino, central Italy): new constraints on the Last Interglacial sea level markers

We present a combined geomorphological and biochronological study aimed at providing age constraints to the deposits forming a wide paleo-surface in the coastal area of the Tyrrhenian Sea, south of Anzio promontory (central Italy). We review the faunal assemblage recovered in Campoverde, evidencing the occurrence of the modern fallow deer subspecies Dama dama dama, which in peninsular Italy is not present before MIS 5e, providing a post-quem terminus of 125 ka for the deposit hosting the fossil remains. The geomorphological reconstruction shows that Campoverde is located within the highest of three paleosurfaces progressively declining towards the present coast, at average elevations of 36, 26 and 15 m a.s.l. The two lowest paleosurfaces match the elevation of the previously recognized marine terraces in this area; we define a new, upper marine terrace corresponding to the 36 m paleosurface, which we name Campoverde complex. Based on the provided evidence of an age as young as MIS 5e for this terrace, we discuss the possibility that previous identification of a tectonically stable MIS 5e coastline ranging 10–8 m a.s.l. in this area should be revised, with significant implications on assessment of the amplitude of sea-level oscillations during the Last Interglacial in the Mediterranean Sea.

Beyond the Last Glacial Maximum: Island endemism is best explained by long‐lasting archipelago configurations

AbstractAimTo quantify the influence of past archipelago configuration on present‐day insular biodiversity patterns, and to compare the role of long‐lasting archipelago configurations over the Pleistocene to configurations of short duration such as at the Last Glacial Maximum (LGM) and the present‐day.Location53 volcanic oceanic islands from 12 archipelagos worldwide—Azores, Canary Islands, Cook Islands, Galápagos, Gulf of Guinea, Hawaii, Madeira, Mascarenes, Pitcairn, Revillagigedo, Samoan Islands and Tristan da Cunha.Time periodThe last 800 kyr, representing the nine most recent glacial–interglacial cycles.Major taxa studiedLand snails and angiosperms.MethodsSpecies richness data for land snails and angiosperms were compiled from existing literature and species checklists. We reconstructed archipelago configurations at the following sea levels: the present‐day high interglacial sea level, the intermediate sea levels that are representative of the Pleistocene and the low sea levels of the LGM. We fitted two alternative linear mixed models for each archipelago configuration using the number of single‐island endemic, multiple‐island endemic and (non‐endemic) native species as a response. Model performance was assessed based on the goodness‐of‐fit of the full model, the variance explained by archipelago configuration and model parsimony.ResultsSingle‐island endemic richness in both taxonomic groups was best explained by intermediate palaeo‐configuration (positively by area change, and negatively by palaeo‐connectedness), whereas non‐endemic native species richness was poorly explained by palaeo‐configuration. Single‐island endemic richness was better explained by intermediate archipelago configurations than by the archipelago configurations of the LGM or present‐day.Main conclusionsArchipelago configurations at intermediate sea levels—which are representative of the Pleistocene—have left a stronger imprint on single‐island endemic richness patterns on volcanic oceanic islands than extreme archipelago configurations that persisted for only a few thousand years (such as the LGM). In understanding ecological and evolutionary dynamics of insular biota it is essential to consider longer‐lasting environmental conditions, rather than extreme situations alone.