Relative sea-level history of Marguerite Bay, Antarctic Peninsula derived from optically stimulated luminescence-dated beach cobbles

Calmette Bay within Marguerite Bay along the western side of the Antarctic Peninsula contains one of the most continuous flights of raised beaches described to date in Antarctica. Raised beaches extend to 40.8 m above sea level (masl) and are thought to reflect glacial isostatic adjustment due to the retreat of the Antarctic Peninsula Ice Sheet. Using optically stimulated luminescence (OSL), we dated quartz extracts from cobble surfaces buried in raised beaches at Calmette Bay. The beaches are separated into upper and lower beaches based on OSL ages, geomorphology, and sedimentary fabric. The two sets of beaches are separated by a prominent scarp. One of our OSL ages from the upper beaches dates to 9.3 thousand years ago (ka; as of 1950) consistent with previous extrapolation of sea-level data and the time of ice retreat from inner Marguerite Bay. However, four of the seven ages from the upper beaches date to the timing of glaciation. We interpret these ages to represent reworking of beaches deposited prior to the Last Glacial Maximum (LGM) by advancing and retreating LGM ice. Ages from the lower beaches record relative sea-level fall due to Holocene glacial-isostatic adjustment. We suggest a Holocene marine limit of 21.7 masl with an age of 5.5-7.3 ka based on OSL ages from Calmette Bay and other sea-level constraints in the area. A marine limit at 21.7 masl implies half as much relative sea-level change in Marguerite Bay during the Holocene as suggested by previous sea-level reconstructions. No evidence for a relative sea-level signature of neoglacial events, such as a decrease followed by an increase in RSL fall due to ice advance and retreat associated with the Little Ice Age, is found within Marguerite Bay indicating either: (1) no significant neoglacial advances occurred within Marguerite Bay; (2) rheological heterogeneity allows part of the Antarctic Peninsula (i.e. the South Shetland Islands) to respond to rapid ice mass changes while other regions are incapable of responding to short-lived ice advances; or (3) the magnitude of neoglacial events within Marguerite Bay is too small to resolve through relative sea-level reconstructions. Although the application of reconstructing sea-level histories using OSL-dated raised beach deposits provides a better understanding of the timing and nature of relative sea-level change in Marguerite Bay, we highlight possible problems associated with using raised beaches as sea-level indices due to post-depositional reworking by storm waves.

Probability inundation map for a scenario combining storm surge and sea-level rise

Coastal managers require reliable spatial data on the extent and timing of potential coastal inundation, particularly in a changing climate. Most sea level rise (SLR) vulnerability assessments are undertaken using the easily implemented bathtub approach, where areas adjacent to the sea and below a given elevation are mapped using a deterministic line dividing potentially inundated from dry areas. This method only requires elevation data usually in the form of a digital elevation model (DEM). However, inherent errors in the DEM and spatial analysis of the bathtub model propagate into the inundation mapping. The aim of this study was to assess the impacts of spatially variable and spatially correlated elevation errors in high-spatial resolution DEMs for mapping coastal inundation. Elevation errors were best modelled using regression-kriging. This geostatistical model takes the spatial correlation in elevation errors into account, which has a significant impact on analyses that include spatial interactions, such as inundation modelling. The spatial variability of elevation errors was partially explained by land cover and terrain variables. Elevation errors were simulated using sequential Gaussian simulation, a Monte Carlo probabilistic approach. 1,000 error simulations were added to the original DEM and reclassified using a hydrologically correct bathtub method. The probability of inundation to a scenario combining a 1 in 100 year storm event over a 1 m SLR was calculated by counting the proportion of times from the 1,000 simulations that a location was inundated. This probabilistic approach can be used in a risk-aversive decision making process by planning for scenarios with different probabilities of occurrence. For example, results showed that when considering a 1% probability exceedance, the inundated area was approximately 11% larger than mapped using the deterministic bathtub approach. The probabilistic approach provides visually intuitive maps that convey uncertainties inherent to spatial data and analysis.

Benthic foraminifera and sea level analyses from IODP Hole 310-M0005D

The course of sea-level fluctuations during Termination II (TII; the penultimate deglaciation), which is critical for understanding ice-sheet dynamics and suborbital climate variability, has yet to be established. This is partly because most shallow-water sequences encompassing TII were eroded during sea-level lowstands of the last glacial period or were deposited below the present sea level. Here we report a new sequence recording sea-level changes during TII in the Pleistocene sequence at Hole M0005D (water depth: 59.63 m below sea level [mbsl]) off Tahiti, French Polynesia, which was drilled during Integrated Ocean Drilling Program Expedition 310. Lithofacies variations and stratigraphic changes in the taxonomic composition, preservation states, and intraspecific test morphology of large benthic foraminifers indicate a deepening-upward sequence in the interval from Core 310-M0005D-26R (core depth: 134 mbsl) through -16R (core depth: 106 mbsl). Reconstruction of relative sea levels, based on paleodepth estimations using large benthic foraminifers, indicated a rise in sea level of about 90 m during this interval, suggesting its correlation with one of the terminations. Assuming that this rise in sea level corresponds to that during TII, after correcting for subsidence since the time of deposition, a highstand sea-level position would be 2 ± 15 m above present sea level (masl), which is generally consistent with highstand sea-level positions in MIS 5e (4 ± 2 masl). If this rise in sea level corresponds to that during older terminations, the subsidence-corrected highstand sea-level positions (30 ± 15 masl for Termination III and 54 ± 15 masl for Termination IV) are not consistent with reported ranges of interglacial sea-level highstands (-18 to 15 masl). Therefore, the studied interval likely records the rise in sea level and associated environmental changes during TII. In particular, the intervening cored materials between the two episodes of sea-level rise found in the studied interval might record the sea-level reversal event during TII. This conclusion is consistent with U/Th ages of around 133 ka, which were obtained from slightly diagenetically altered (i.e., < 1% calcite) in situ corals in the studied interval (Core 310-M0005D-20R [core depth: 118 mbsl]). This study also suggests that our inverse approach to correlate a stratigraphic interval with an approximate time frame could be useful as an independent check on the accuracy of uranium-series dating, which has been applied extensively to fossil corals in late Quaternary sea-level studies.

Solid phase analysis of sediment core GeoB9309-1 from the Zambezi (or Mozambique) deep-sea fan

The Zambezi deep-sea fan, the largest of its kind along the east African continental margin, is poorly studied to date, despite its potential to record marine and terrestrial climate signals in the southwest Indian Ocean. Therefore, gravity core GeoB 9309-1, retrieved from 1219 m water depth, was investigated for various geophysical (magnetic susceptibility, porosity, colour reflectance) and geochemical (pore water and sediment geochemistry, Fe and P speciation) properties. Onboard and onshore data documented a sulphate/methane transition (SMT) zone at ~ 450-530 cm sediment depth, where the simultaneous consumption of pore water sulphate and methane liberates hydrogen sulphide and bi-carbonate into the pore space. This leads to characteristic changes in the sediment and pore water chemistry, as the reduction of primary Fe (oxyhydr)oxides, the precipitation of Fe sulphides, and the mobilization of Fe (oxyhydr)oxide-bound P. These chemical processes also lead to a marked decrease in magnetic susceptibility. Below the SMT, we find a reduction of porosity, possibly due to pore space cementation by authigenic minerals. Formation of the observed geochemical, magnetic and mineralogical patterns requires a fixation of the SMT at this distinct sediment depth for a considerable time-which we calculated to be ~ 10 000 years assuming steady-state conditions-following a period of rapid upward migration towards this interval. We postulate that the worldwide sea-level rise at the last glacial/interglacial transition (~ 10 000 years B.P.) most probably caused the fixation of the SMT at its present position, through drastically reduced sediment delivery to the deep-sea fan. In addition, we report an internal redistribution of P occurring around the SMT, closely linked to the (de)coupling of sedimentary Fe and P, and leaving a characteristic pattern in the solid P record. By phosphate re-adsorption onto Fe (oxyhydr)oxides above, and formation of authigenic P minerals (e.g. vivianite) below the SMT, deep-sea fan deposits may potentially act as long-term sinks for P.

Radiocarbon dating, oxygen isotopes, and sea-level reconstruction of sediment core GeoB5844-2 from the northern end of the Red Sea

Based on a radiocarbon and paleomagnetically dated sediment record from the northern Red Sea and the exceptional sensitivity of the regional changes in the oxygen isotope composition of sea water to the sea-level-dependent water exchange with the Indian Ocean, we provide a new global sea-level reconstruction spanning the last glacial period. The sea-level record has been extracted from the temperature-corrected benthic stable oxygen isotopes using coral-based sea-level data as constraints for the sea-level/oxygen isotope relationship. Although, the general features of this millennial-scale sea-level records have strong similarities to the rather symmetric and gradual Southern Hemisphere climate patterns, we observe, in constrast to previous findings, pronounced sea level rises of up to 25 m to generally correspond with Northern Hemisphere warmings as recorded in Greenland ice-core interstadial intervals whereas sea-level lowstands mostly occur during cold phases. Corroborated by CLIMBER-2 model results, the close connection of millennial-scale sea-level changes to Northern Hemisphere temperature variations indicates a primary climatic control on the mass balance of the major Northern Hemisphere ice sheets and does not require a considerable Antarctic contribution.

Stable isotope and geochemical record of sediments from the Bering Sea

here is controversy over the role of marine methane hydrates in atmospheric methane concentrations and climate change during the last glacial period. In this study of two sediment cores from the southeast Bering Sea (700 m and 1467 m water depth), we identify multiple episodes during the last glacial period of intense methane flux reaching the seafloor. Within the uncertainty of the radiocarbon age model, the episodes are contemporaneous in the two cores and have similar timing and duration as Dansgaard-Oeschger events. The episodes are marked by horizons of sediment containing 13C-depleted authigenic carbonate minerals; 13C-depleted archaeal and bacterial lipids, which resemble those found in ANME-1 type anaerobic methane oxidizing microbial consortia; and changes in the abundance and species distribution of benthic foraminifera. The similar timing and isotopic composition of the authigenic carbonates in the two cores is consistent with a region-wide increase in the upward flux of methane bearing fluids. This study is the first observation outside Santa Barbara Basin of pervasive, repeated methane flux in glacial sediments. However, contrary to the "Clathrate Gun Hypothesis" (Kennett et al., 2003), these coring sites are too deep for methane hydrate destabilization to be the cause, implying that a much larger part of the ocean's sedimentary methane may participate in climate or carbon cycle feedback at millennial timescales. We speculate that pulses of methane in these opal-rich sediments could be caused by the sudden release of overpressure in pore fluids that builds up gradually with silica diagenesis. The release could be triggered by seismic shaking on the Aleutian subduction zone caused by hydrostatic pressure increase associated with sea level rise at the start of interstadials.

Orbital-resolution Pliocene-Pleistocene simulations of CO2, sea level, global temperature and benthic d18O, and d11B-based proxy-CO2 data

During the past five million yrs, benthic d18O records indicate a large range of climates, from warmer than today during the Pliocene Warm Period to considerably colder during glacials. Antarctic ice cores have revealed Pleistocene glacial-interglacial CO2 variability of 60-100 ppm, while sea level fluctuations of typically 125 m are documented by proxy data. However, in the pre-ice core period, CO2 and sea level proxy data are scarce and there is disagreement between different proxies and different records of the same proxy. This hampers comprehensive understanding of the long-term relations between CO2, sea level and climate. Here, we drive a coupled climate-ice sheet model over the past five million years, inversely forced by a stacked benthic d18O record. We obtain continuous simulations of benthic d18O, sea level and CO2 that are mutually consistent. Our model shows CO2 concentrations of 300 to 470 ppm during the Early Pliocene. Furthermore, we simulate strong CO2 variability during the Pliocene and Early Pleistocene. These features are broadly supported by existing and new d11B-based proxy CO2 data, but less by alkenone-based records. The simulated concentrations and variations therein are larger than expected from global mean temperature changes. Our findings thus suggest a smaller Earth System Sensitivity than previously thought. This is explained by a more restricted role of land ice variability in the Pliocene. The largest uncertainty in our simulation arises from the mass balance formulation of East Antarctica, which governs the variability in sea level, but only modestly affects the modeled CO2 concentrations.

Gridded climate data from 5 Global Climate Models (GCM) of the Last Glacial Maximum (LGM) downscaled to 30 arc seconds for Europe, with links to NetCDF files

Studies on the impact of historical, current and future global change require very high-resolution climate data (less or equal 1km) as a basis for modelled responses, meaning that data from digital climate models generally require substantial rescaling. Another shortcoming of available datasets on past climate is that the effects of sea level rise and fall are not considered. Without such information, the study of glacial refugia or early Holocene plant and animal migration are incomplete if not impossible. Sea level at the last glacial maximum (LGM) was approximately 125m lower, creating substantial additional terrestrial area for which no current baseline data exist. Here, we introduce the development of a novel, gridded climate dataset for LGM that is both very high resolution (1km) and extends to the LGM sea and land mask. We developed two methods to extend current terrestrial precipitation and temperature data to areas between the current and LGM coastlines. The absolute interpolation error is less than 1°C and 0.5 °C for 98.9% and 87.8% of all pixels for the first two 1 arc degree distance zones. We use the change factor method with these newly assembled baseline data to downscale five global circulation models of LGM climate to a resolution of 1km for Europe. As additional variables we calculate 19 'bioclimatic' variables, which are often used in climate change impact studies on biological diversity. The new LGM climate maps are well suited for analysing refugia and migration during Holocene warming following the LGM.

Sedimentological, mineralogical, biogeochemical and micropaleontological parameters on sediment core PS69/849-2, Burton Basin, MacRobertson Shelf, East Antarctica

A multi-proxy study including sedimentological, mineralogical, biogeochemical and micropaleontological methods was conducted on sediment core PS69/849-2 retrieved from Burton Basin, MacRobertson Shelf, East Antarctica. The goal of this study was to depict the deglacial and Holocene environmental history of the MacRobertson Land-Prydz Bay region. A special focus was put on the timing of ice-sheet retreat and the variability of bottom-water formation due to sea ice formation through the Holocene. Results from site PS69/849-2 provide the first paleo-environmental record of Holocene variations in bottom-water production probably associated to the Cape Darnley polynya, which is the second largest polynya in the Antarctic. Methods included end-member modeling of laser-derived high-resolution grain size data to reconstruct the depositional regimes and bottom-water activity. The provenance of current-derived and ice-transported material was reconstructed using clay-mineral and heavy-mineral analysis. Conclusions on biogenic production were drawn by determination of biogenic opal and total organic carbon. It was found that the ice shelf front started to retreat from the site around 12.8 ka BP. This coincides with results from other records in Prydz Bay and suggests warming during the early Holocene optimum next to global sea level rise as the main trigger. Ice-rafted debris was then supplied to the site until 5.5 cal. ka BP, when Holocene global sea level rise stabilized and glacial isostatic rebound on MacRobertson Land commenced. Throughout the Holocene, three episodes of enhanced bottom-water activity probably due to elevated brine rejection in Cape Darnley polynya occured between 11.5 and 9 cal. ka BP, 5.6 and 4.5 cal. ka BP and since 1.5 cal. ka BP. These periods are related to shifts from warmer to cooler conditions at the end of Holocene warm periods, in particular the early Holocene optimum, the mid-Holocene warm period and at the beginning of the neoglacial. In contrast, between 7.7 and 6.7 cal. ka BP, brine rejection shut down, maybe owed to warm conditions and pronounced open-water intervals.

Location and ages of terrestrial peatlands located on continental shelves and in coastal sediments

This synthesis dataset contains records of freshwater peat and lake sediments from continental shelves and coastal areas. Information included is site location (when available), thickness and description of terrestrial sediments as well as underlying and overlying sediments, dates (when available), and references.

Sea-surface temperature and salinity reconstruction for sediment core GeoB5844-2

We present high-resolution paleoceanographic records of surface and deep water conditions within the northern Red Sea covering the last glacial maximum and termination I using alkenone paleothermometry, stable oxygen isotopes, and sediment compositional data. Paleoceanographic records in the restricted desert-surrounded northern Red Sea are strongly affected by the stepwise sea level rise and appear to record and amplify well-known millennial-scale climate events from the North Atlantic realm. During the last glacial maximum (LGM), sea surface temperatures were about 4°C cooler than the late Holocene. Pronounced coolings associated with Heinrich event 1 (~2°C below the LGM level) and the Younger Dryas imply strong atmospheric teleconnections to the North Atlantic. Owing to the restricted exchange with the Indian Ocean, Red Sea salinity is particularly sensitive to changes in global sea level. Paleosalinities exceeded 50 psu during the LGM. A pronounced freshening of the surface waters is associated with the meltwater peaks MWP1a and MWP1b owing to an increased surface-near inflow of "normal" saline water from the Indian Ocean. Vertical delta18O gradients are also increased during these phases, indicating stronger surface water stratification. The combined effect of deglacial changes in sea surface temperature and salinity on water column stratification initiated the formation of two sapropel layers, which were deposited under almost anoxic condition in a stagnant water body.

Grain size distribution and sea level changes in Gåsfjärden fjord

We analysed long-term variations in grain-size distribution in sediments from Gåsfjärden, a fjord-like inlet on the south-west Baltic Sea, and explored potential drivers of the recorded changes in sediment grain-size data. Over the last 5.4 thousand years (ka), the relative sea level decreased 17 m in the study region, caused by isostatic land uplift. As a consequence, Gåsfjärden has been transformed from an open coastal setting into a semi-closed inlet surrounded on the east by numerous small islands. To quantitatively estimate the morphological changes in Gåsfjärden over the last 5.4 ka and to further link the changes to our grain-size data, a digital elevation model (DEM)-based openness index was calculated. In the period between 5.4 and 4.4 ka BP, the inlet was characterised by the largest openness index. During this interval, the highest sand contents (~0.4 %) and silt/clay ratios (~0. 3) in the sediment sequence were recorded, indicating relatively high bottom water energy. After 4.4 ka BP, the average sand content was halved to ~0.2 % and the silt/clay ratios showed a significant decreasing trend over the last 4 ka. These changes are found to be associated with the gradual embayment of Gåsfjärden as represented in the openness index. The silt/clay ratios exhibited a delayed and slower change compared with the sand contents, which further suggest that finer particles are less sensitive to changes in hydrodynamic energy. Our DEM-based coastal openness index has proved to be a useful tool for interpreting the sedimentary grain-size record.

Facies distribution in CRP sediment cores from the Ross Sea, Antarctica

Clasts from the Cape Roberts Project cores CRP-2/2A and CRP-3 provide indications of glacially influenced depositional environments in Oligocene and Miocene strata in the western Victoria Land Basin, Antarctica. CRP-2/2A is interpreted to represent strongly glacially influenced, unconformity bound depositional sequences produced by repeated advance and retreat of floating and grounded ice across the shelf. A similar interpretation is extended to the upper 330 meters of the CRP-3 core, but the lower part of the core records shallow marine deposition with significantly less glacial influence. Clast shape analysis from selected coarse-grained facies throughout the cored interval indicates that most clasts are glacially sourced, with little distinction between diamictite and conglomeratic facies. Three dimensional clast fabric analysis from units immediately above sequence boundaries generally display weak or random fabrics and do not suggest that grounded ice actually reached the drillsite at these intervals. Striated and outsized clasts present in fine-grained lithofacies throughout the cores provide further evidence of sub-glacially transported sediment and iceberg rafting. The distribution of these striated and out-sized clasts indicate that a significant glacial influence persisted through most of the time represented by the cores with glaciers actively calving at sea-level introducing ice-berg rafted glacial debris even in the earliest Oligocene.

Age determinations from Bykovsky Peninsula, Siberia

A wide variety of environmental records is necessary for analysing and understanding the complex Late Quaternary dynamics of permafrost-dominated Arctic landscapes. A NE Siberian periglacial key region was studied in detail using sediment records, remote sensing data, and terrain modelling, all incorporated in a geographical information system (GIS). The study area consists of the Bykovsky Peninsula and the adjacent Khorogor Valley in the Kharaulakh Ridge situated a few kilometres southeast of the Lena Delta. In this study a comprehensive cryolithological database containing information from 176 sites was compiled. The information from these sites is based on the review of previously published borehole data, outcrop profiles, surface samples, and our own field data. These archives cover depositional records of three periods: from Pliocene to Early Pleistocene, the Late Pleistocene and the Holocene. The main sediment sequences on the Bykovsky Peninsula consist of up to 50 m thick ice-rich permafrost deposits (Ice Complex) that were accumulated during the Late Pleistocene. They were formed as a result of nival processes around extensive snowfields in the Kharaulakh Ridge, slope processes in these mountains (such as in the Khorogor Valley), and alluvial/proluvial sedimentation in a flat accumulation plain dominated by polygonal tundra in the mountain foreland (Bykovsky Peninsula). During the early to middle Holocene warming, a general landscape transformation occurred from an extensive Late Pleistocene accumulation plain to a strongly thermokarst-dominated relief dissected by numerous depressions. Thermokarst subsidence had an enormous influence on the periglacial hydrological patterns, the sediment deposition, and on the composition and distribution of habitats. Climate deterioration, lake drainage, and talik refreezing occurred during the middle to late Holocene. The investigated region was reached by the post-glacial sea level rise during the middle Holocene, triggering thermo-abrasion of ice-rich coasts and the marine inundation of thermokarst depressions.