ALBERT

Description: We present a compilation and analysis of 1099 Holocene relative shore-level (RSL) indicators including 867 relative sea-level data points and 232 data points from the Ancylus Lake and the following transitional phase from 10.7 to 8.5 ka BP located around the Baltic Sea. The spatial distribution covers the Baltic Sea and near-coastal areas fairly well, but some gaps remain mainly in Sweden. RSL data follow the standardized HOLSEA format and, thus, are ready for spatially comprehensive applications in, e.g., glacial isostatic adjustment (GIA) modelling. Sampling method The data set is a compilation of rather different samples from geological, geomorphological and archaeological studies. Most of the data was already published in different formats. In this compilation we homogenized the meta information of the available information according to the HOLSEA database format, https://www.holsea.org/archive-your-data, which is a modification of the recommendations given in Hijma et al. (2015). In addition to the reformatting, the majority of samples with radiocarbon dating were recalibrated with oxcal-software using the calib13 and marine13 curves. Furthermore, all sample descriptions were critically checked for consistency in positioning, levelling and indicative meaning by experts of the respective geographic region see Supplement 2. Analytical method In principle, it is a compilation, recalibration and revision of already published data. Data Processing Data of individual compilations were revised and imported into a relational database system. Therein, the data was transferred into the HOLSEA format by specified rules. By this procedure, a homogeneous categorisation was achieved without losing the original data. Also this is stored in the relational database system allowing for later updates of the transfer procedure or a recalibration of the data. Description of data table HOLSEA-baltic-yymmdd.xlsx The workbook in excel format contains 5 sheets, see https://www.holsea.org/archive-your-data: · Long-form, containing the complete information available for each sample · Short-form, a subset of attributes of the Long-form sheet · Radiocarbon, containing the radiocarbon dating information of the respective samples · U-series, a corresponding table containing the respective information of Uranium dating · References, a complete reference list of the primary publications in which the individual data sampling is described. All online sources for the compilation are included in the metadata. A full list of source references is provided in the data description file.

Description: We provide the model results of the manuscript "Glacial-isostatic adjustment models using geodynamically constrained 3D Earth structures" (Bagge et al. 2020, Paper) including the (1) predicted relative sea-level and (2) applied sea-level data. The predicted relative-sea level is calculated with the VIscoelastic Lithosphere and MAntle model VILMA (Klemann et al. 2008, 2015, Martinec et al. 2018, Hagedoorn et al. 2007, Martinec & Hagedoorn 2005, Kendall et al. 2005). The glacial-isostatic adjustment models uses different Earth structures (3D, 1D global mean and 1D regionally adapted; Bagge et al. 2020, Paper; Bagge et al. 2020, https://doi.org/10.5880/GFZ.1.3.2020.004) and ice histories (ICE-5G, Peltier 2004; ICE-6G, Peltier et al. 2015, Argus et al. 2014; NAICE, Gowan et al. 2016) resulting in 44 3D models, 54 1D global mean models and 162 1D regionally adapted models. For more information on model description and input data see Bagge et al. (2020, Paper) and Bagge at al. (2020, https://doi.org/10.5880/GFZ.1.3.2020.004). The provided output data include (1a) the global distribution of predicted relative-sea level at 14 kilo years before present as ensemble range of the 3D GIA models for three ice histories as netCDF files, (1b) the predicted relative-sea level at eight locations at 14 kilo years before present for all models as ASCII file and (1c) the predicted relative sea-level for the deglaciation period for all models as ASCII files. Eight locations include Churchill, Angermanland, Ross Sea (Antarctica), San Jorge Gulf (Patagonia), Central Oregon Coast, Rao-Gandon Area (Senegal), Singapore and Pioneer Bay (Queensland, Australia). (2) The about 520 applied sea-level data provide information on time, relative sea-level and type of sea-level data. They are extracted for the eight locations from the GFZ database using SLIVisu (Unger et al. 2012, 2018) and provided as ACSII files.

Description: We provide 18 3D Earth structures on a global grid. This supplementary material of the manuscript includes (1) 18 netcdf files of the 3D Earth structures and (2) 72 figures that visualize the lithospheric thickness, lateral average viscosity of the asthenosphere, transition zone and upper mantle for all 18 3D Earth structures. The Earth structures were derived from seismic tomography models (Schaeffer & Lebedev 2013, 2010 update of Grand 2002) and, under consideration of geodynamic constrains, transferred to viscosity (Steinberger 2016, Steinberger & Calderwood 2006). The 18 Earth structures vary in conversion from seismic velocity to viscosity. Detailed description of the procedure can be found in the corresponding manuscript (Bagge et al. 2020a), where the Earth structure data were applied to the glacial-isostatic adjustment model VILMA (Klemann et al. 2008, 2015, Martinec et al. 2018) to predict the relative sea-level during the last deglaciation. The netCDF files are provided on a Gaussian grid of 256x512 grid points. Each Earth structure consist of 167 layers, while lateral variations in Earth structure are considered for 114 layers between surface and 870 km depth and radially symmetric layers are considered for 50 layers from 870 km to the Earth’s core. The Earth structure is given as logarithmic viscosity in log10[Viscosity(Pa s)]. To visualize the global 3D structures, we calculated the lithospheric thickness and average viscosity of the asthenosphere, upper mantle and transition zone. The lithospheric thickness is defined as minimum depth with a viscosity 〈 10^23.5 Pa s, the asthenosphere is defined between the base of the lithosphere and 225 km depth, the upper mantle between 225 km and 410 km and the transition zone between 410 km and 670 km depth.

Description: We provide present-day glacial isostatic adjustment (GIA) gravity changes simulated with the numerical model VILMA. The effects of Earth and ocean pole tide due to rotational deformation (considered in VILMA) were removed. The dataset contains the solutions for 56 GIA model ensemble members including 54 3D models and 2 1D models. The results are provided as Stokes coefficients with a resolution of degree/order 170.