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

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