Timing of the most recent Neoglacial advance and retreat in Potter Cove, Antarctic Peninsula

The timing of the most recent Neoglacial advance in the Antarctic Peninsula is important for establishing global climate teleconnections and providing important post-glacial rebound corrections to gravity-based satellite measurements of ice loss. However, obtaining accurate ages from terrestrial geomorphic and sedimentary indicators of the most recent Neoglacial advance in Antarctica has been hampered by the lack of historical records and the difficulty of dating materials in Antarctica. Here we use a new approach to dating flights of raised beaches in the South Shetland Islands of the northern Antarctic Peninsula to bracket the age of a Neoglacial advance that occurred between 1500 and 1700 AD, broadly synchronous with compilations for the timing of the Little Ice Age in the northern hemisphere. Our approach is based on optically stimulated luminescence of the underside of buried cobbles to obtain the age of beaches previously shown to have been deposited immediately inside and outside the moraines of the most recent Neoglacial advance. In addition, these beaches mark the timing of an apparent change in the rate of isostatic rebound thought to be in response to the same glacial advance within the South Shetland Islands. We use a Maxwell viscoelastic model of glacial-isostatic adjustment (GIA) to determine whether the rates of uplift calculated from the raised beaches are realistic given the limited constraints on the ice advance during this most recent Neoglacial advance. Our rebound model suggests that the subsequent melting of an additional 16-22% increase in the volume of ice within the South Shetland Islands would result in a subsequent uplift rate of 12.5 mm/yr that lasted until 1840 AD resulting in a cumulative uplift of 2.5 m. This uplift rate and magnitude are in close agreement with observed rates and magnitudes calculated from the raised beaches since the most recent Neoglacial advance along the South Shetland Islands and falls within the range of uplift rates from similar settings such as Alaska.

Validated sample ages from 823 EPD sites

The number of well-dated pollen diagrams in Europe has increased considerably over the last 30 years and many of them have been submitted to the European Pollen Database (EPD). This allows for the construction of increasingly precise maps of Holocene vegetation change across the continent. Chronological information in the EPD has been expressed in uncalibrated radiocarbon years, and most chronologies to date are based on this time scale. Here we present new chronologies for most of the datasets stored in the EPD based on calibrated radiocarbon years. Age information associated with pollen diagrams is often derived from the pollen stratigraphy itself or from other sedimentological information. We reviewed these chronological tie points and assigned uncertainties to them. The steps taken to generate the new chronologies are described and the rationale for a new classification system for age uncertainties is introduced. The resulting chronologies are fit for most continental-scale questions. They may not provide the best age model for particular sites, but may be viewed as general purpose chronologies. Taxonomic particularities of the data stored in the EPD are explained. An example is given of how the database can be queried to select samples with appropriate age control as well as the suitable taxonomic level to answer a specific research question.

Age dating and isotopic analyses of sediment core RC11-120

Using the concept of 'orbital tuning', a continuous, high-resolution deep-sea chronostratigraphy has been developed spanning the last 300,000 yr. The chronology is developed using a stacked oxygen-isotope stratigraphy and four different orbital tuning approaches, each of which is based upon a different assumption concerning the response of the orbital signal recorded in the data. Each approach yields a separate chronology. The error measured by the standard deviation about the average of these four results (which represents the 'best' chronology) has an average magnitude of only 2500 yr. This small value indicates that the chronology produced is insensitive to the specific orbital tuning technique used. Excellent convergence between chronologies developed using each of five different paleoclimatological indicators (from a single core) is also obtained. The resultant chronology is also insensitive to the specific indicator used. The error associated with each tuning approach is estimated independently and propagated through to the average result. The resulting error estimate is independent of that associated with the degree of convergence and has an average magnitude of 3500 yr, in excellent agreement with the 2500-yr estimate. Transfer of the final chronology to the stacked record leads to an estimated error of +/-1500 yr. Thus the final chronology has an average error of +/-5000 yr.

Physical properties and sedimentology of 4 profiles fom the South Atlantic

Ultrasonic P wavc transmission seismograms recorded on sediment cores have been analyzed to study the acoustic and estimate the clastic properties of marine sediments from different provinces dominated by terrigenous, calcareous, amI diatomaceous sedimentation. Instantaneous frequencies computed from the transmission seismograms are displayed as gray-shaded images to give an acoustic overview of the lithology of each core. Ccntirneter-scale variations in the ultrasonic waveforms associated with lithological changes are illustrated by wiggle traces in detail. Cross-correlation, multiple-filter, and spectral ratio techniques are applied to derive P wave velocities and attenuation coefficients. S wave velocities and attenuation coefficients, elastic moduli, and permeabilities are calculated by an inversion scheme based on the Biot-Stoll viscoelastic model. Together wilh porosity measurements, P and S wave scatter diagrams are constructed to characterize different sediment types by their velocity- and attenuation-porosity relationships. They demonstrate that terrigenous, calcareous, and diatomaceous sediments cover different velocity- and attenuation-porosity ranges. In terrigcnous sediments, P wave vclocities and attenuation coefficients decrease rapidly with increasing porosity, whereas S wave velocities and shear moduli are very low. Calcareous sediments behave similarly at relatively higher porosities. Foraminifera skeletons in compositions of terrigenous mud and calcareous ooze cause a stiffening of the frame accompanied by higher shear moduli, P wave velocities, and attenuation coefficients. In diatomaceous ooze the contribution of the shear modulus becomes increasingly important and is controlled by the opal content, whereas attenuation is very low. This leads to the opportunity to predict the opal content from nondestructive P wave velocity measurements at centimeter-scale resolution.

Delta 13C measured on alkenones of surface sediment samples GeoB1008-6 to GeoB3603-1 from the South Atlantic

We have analyzed the stable carbon isotopic composition of the diunsaturated C37 alkenone in 29 surface sediments from the equatorial and South Atlantic Ocean. Our study area covers different oceanographic settings, including sediments from the major upwelling regions off South Africa, the equatorial upwelling, and the oligotrophic western South Atlantic. In order to examine the environmental influences on the sedimentary record the alkenone-based carbon isotopic fractionation (Ep) values were correlated with the overlying surface water concentrations of aqueous CO2 ([CO2(aq)]), phosphate, and nitrate. We found Ep positively correlated with 1/[CO2(aq)] and negatively correlated with [PO43-] and [NO3-]. However, the relationship between Ep and 1/[CO2(aq)] is opposite of what is expected from a [CO2(aq)] controlled, diffusive uptake model. Instead, our findings support the theory of Bidigare et al. (1997, doi:10.1029/96GB03939) that the isotopic fractionation in haptophytes is related to nutrient-limited growth rates. The relatively high variability of the Ep-[PO4] relationship in regions with low surface water nutrient concentrations indicates that here other environmental factors also affect the isotopic signal. These factors might be variations in other growth-limiting resources such as light intensity or micronutrient concentrations.

Revised composite depth scale, Fe content and orbital tuning of Middle Eocene Climate Optimum samples from ODP Site 207-1260

A high-resolution stratigraphy is essential toward deciphering climate variability in detail and understanding causality arguments of events in earth history. Because the highly dynamic middle to late Eocene provides a suitable testing ground for carbon cycle models for a waning warm world, an accurate time scale is needed to decode climate-driving mechanisms. Here we present new results from ODP Site 1260 (Leg 207) which covers a unique expanded middle Eocene section (magnetochrons C18r to C20r, late Lutetian to early Bartonian) of the tropical western Atlantic including the chron C19r transient hyperthermal event and the Middle Eocene Climate Optimum (MECO). To establish a detailed cyclostratigraphy we acquired a distinctive iron intensity records by XRF scanning Site 1260 cores. We revise the shipboard composite section, establish a cyclostratigraphy and use the exceptional eccentricity modulated precession cycles for orbital tuning. The new astrochronology revises the age of magnetic polarity chrons C19n to C20n, validates the position of very long eccentricity minima at 40.2 and 43.0 Ma in the orbital solutions, and extends the Astronomically Tuned Geological Time Scale back to 44 Ma. For the first time the new data provide clear evidence for an orbital pacing of the chron C19r event and a likely involvement of the very long eccentricity cycle contributing to the evolution of the MECO.

Secchi depth maps for the Baltic Sea, 2003-2012, with links to GeoTIFs

Secchi depth is a measure of water transparency. In the Baltic Sea region, Secchi depth maps are used to assess eutrophication and as input for habitat models. Due to their spatial and temporal coverage, satellite data would be the most suitable data source for such maps. But the Baltic Sea's optical properties are so different from the open ocean that globally calibrated standard models suffer from large errors. Regional predictive models that take the Baltic Sea's special optical properties into account are thus needed. This paper tests how accurately generalized linear models (GLMs) and generalized additive models (GAMs) with MODIS/Aqua and auxiliary data as inputs can predict Secchi depth at a regional scale. It uses cross-validation to test the prediction accuracy of hundreds of GAMs and GLMs with up to 5 input variables. A GAM with 3 input variables (chlorophyll a, remote sensing reflectance at 678 nm, and long-term mean salinity) made the most accurate predictions. Tested against field observations not used for model selection and calibration, the best model's mean absolute error (MAE) for daily predictions was 1.07 m (22%), more than 50% lower than for other publicly available Baltic Sea Secchi depth maps. The MAE for predicting monthly averages was 0.86 m (15%). Thus, the proposed model selection process was able to find a regional model with good prediction accuracy. It could be useful to find predictive models for environmental variables other than Secchi depth, using data from other satellite sensors, and for other regions where non-standard remote sensing models are needed for prediction and mapping. Annual and monthly mean Secchi depth maps for 2003-2012 come with this paper as Supplementary materials.

Sulfur isotope composition of pore fluids at ODP Sites 128-798 and 128-799

Micro-crystalline barites recovered by deep-sea drilling from Site 684 on the Peru margin and Site 799 in the Japan Sea are highly enriched in the heavy sulfur isotope relative to seawater ( d34S up to +84?). This isotopic composition is consistent with remobilization of biogenic barite triggered by sulfate reduction, and subsequent reprecipitation as a diagenetic barite front. The high levels of barium sulfate in these deposits (10-50%) cannot be explained by a diffusive transport model in sediments experiencing a constant rate of sedimentation. When sedimentation rates change radically, the barite front will remain at a given depth interval leading to large accumulations of barium sulfate. Such conditions may have generated the barite deposits at Site 799. At Site 684, on the other hand, there is evidence that the barite deposits are a result of the tectonically-driven advection of sulfate-bearing fluids through the sediment column.

E-RUN version 1.1: Observational gridded runoff estimates for Europe, link to data in NetCDF format (69 MB)

River runoff is an essential climate variable as it is directly linked to the terrestrial water balance and controls a wide range of climatological and ecological processes. Despite its scientific and societal importance, there are to date no pan-European observation-based runoff estimates available. Here we employ a recently developed methodology to estimate monthly runoff rates on regular spatial grid in Europe. For this we first assemble an unprecedented collection of river flow observations, combining information from three distinct data bases. Observed monthly runoff rates are first tested for homogeneity and then related to gridded atmospheric variables (E-OBS version 12) using machine learning. The resulting statistical model is then used to estimate monthly runoff rates (December 1950 - December 2015) on a 0.5° x 0.5° grid. The performance of the newly derived runoff estimates is assessed in terms of cross validation. The paper closes with example applications, illustrating the potential of the new runoff estimates for climatological assessments and drought monitoring.

E-RUN version 1.0: Observational gridded runoff estimates for Europe, link to data in NetCDF format (68 MB)

River runoff is an essential climate variable as it is directly linked to the terrestrial water balance and controls a wide range of climatological and ecological processes. Despite its scientific and societal importance, there are to date no pan-European observation-based runoff estimates available. Here we employ a recently developed methodology to estimate monthly runoff rates on regular spatial grid in Europe. For this we first collect an unprecedented collection of river flow observations, combining information from three distinct data bases. Observed monthly runoff rates are first tested for homogeneity and then related to gridded atmospheric variables (E-OBS version 11) using machine learning. The resulting statistical model is then used to estimate monthly runoff rates (December 1950-December 2014) on a 0.5° × 0.5° grid. The performance of the newly derived runoff estimates is assessed in terms of cross validation. The paper closes with example applications, illustrating the potential of the new runoff estimates for climatological assessments and drought monitoring.

Global LAI uncertainty and relative uncertainty datasets derived from a Triple Collocation Error Model (TCEM)

Uncertainty information for global leaf area index (LAI) products is important for global modeling studies but usually difficult to systematically obtain at a global scale. Here, we present a new method that cross-validates existing global LAI products and produces consistent uncertainty information. The method is based on a triple collocation error model (TCEM) that assumes errors among LAI products are not correlated. Global monthly absolute and relative uncertainties, in 0.05° spatial resolutions, were generated for MODIS, CYCLOPES, and GLOBCARBON LAI products, with reasonable agreement in terms of spatial patterns and biome types. CYCLOPES shows the lowest absolute and relative uncertainties, followed by GLOBCARBON and MODIS. Grasses, crops, shrubs, and savannas usually have lower uncertainties than forests in association with the relatively larger forest LAI. With their densely vegetated canopies, tropical regions exhibit the highest absolute uncertainties but the lowest relative uncertainties, the latter of which tend to increase with higher latitudes. The estimated uncertainties of CYCLOPES generally meet the quality requirements (± 0.5) proposed by the Global Climate Observing System (GCOS), whereas for MODIS and GLOBCARBON only non-forest biome types have met the requirement. Nevertheless, none of the products seems to be within a relative uncertainty requirements of 20%. Further independent validation and comparative studies are expected to provide a fair assessment of uncertainties derived from TCEM. Overall, the proposed TCEM is straightforward and could be automated for the systematic processing of real time remote sensing observations to provide theoretical uncertainty information for a wider range of land products.

Age and isotopes in bottom sediments from Core VER-99-1-St2, Baikal Lake

A simple, reliable, and efficient method has been elaborated for direct determination of isotopic composition of authigenic uranium in siliceous lacustrine sediments. The method is based on studying kinetics of selective extraction of authigenic uranium from sediments by weak solutions of ammonium hydrocarbonate followed by ICP-MS analysis of nuclides. To estimate contamination of authigenic uranium by terrigenous one contents of 232Th and some other clastogenic elements in the extracts were measured simultaneously. Selectivity of extraction of authigenic uranium from the sediments treated with 1% NH4HCO3 solution appeared to be no worse than 99%. The method was applied to analysis of isotopic composition of authigenic uranium at several key horizons of the earlier dated core from the Baikal Lake. Measurements directly show that 234U/238U values in Baikal water varied depending on climate, which contradicts existing hypotheses. Measured 234U/238U ratios in water of the paleo-Baikal match corresponding values reconstructed from isotopic data for total uranium in the sediments on supposition that U/Th ratio is constant in terrigenous fraction of the sediments. Direct experimental determination of total and authigenic nuclides in sediments enhances potentiality of the method for 234U-230Th dating of non-carbonate lacustrine sediments including those from the Baikal Lake within intervals corresponding to periods of glaciation, when sediments were rich in terrigenous components. Portions of terrigenous and authigenic uranium are well separated and we can study variability of sources of terrigenous matter and refine the earlier model for reconstructing climate humidity in the East Siberia.

Model - data comparison of western African rainfall evolutions during the Last Interglacial (130-115 ka)

This dataset characterizes the evolution of western African precipitation indicated by marine sediment geochemical records in comparison to transient simulations using CCSM3 global climate model throughout the Last Interglacial (130-115 ka). It contains (1) defined tie-points (age models), newly published stable isotopes of benthic foraminifera and Al/Si log-ratios of eight marine sediment cores from the western African margin and (2) annual and seasonal rainfall anomalies (relative to pre-industrial values) for six characteristic latitudinal bands in western Africa simulated by CCSM3 (two transient simulations: one non-accelerated and one accelerated experiment).

Mean dynamic topography and oceanographic parameters estimated from an inverse model and satellite geodesy, with link to model result in one single NetCDF file (392 MB), including the inverse data error covariance

Geostrophic surface velocities can be derived from the gradients of the mean dynamic topography-the difference between the mean sea surface and the geoid. Therefore, independently observed mean dynamic topography data are valuable input parameters and constraints for ocean circulation models. For a successful fit to observational dynamic topography data, not only the mean dynamic topography on the particular ocean model grid is required, but also information about its inverse covariance matrix. The calculation of the mean dynamic topography from satellite-based gravity field models and altimetric sea surface height measurements, however, is not straightforward. For this purpose, we previously developed an integrated approach to combining these two different observation groups in a consistent way without using the common filter approaches (Becker et al. in J Geodyn 59(60):99-110, 2012, doi:10.1016/j.jog.2011.07.0069; Becker in Konsistente Kombination von Schwerefeld, Altimetrie und hydrographischen Daten zur Modellierung der dynamischen Ozeantopographie, 2012, http://nbn-resolving.de/nbn:de:hbz:5n-29199). Within this combination method, the full spectral range of the observations is considered. Further, it allows the direct determination of the normal equations (i.e., the inverse of the error covariance matrix) of the mean dynamic topography on arbitrary grids, which is one of the requirements for ocean data assimilation. In this paper, we report progress through selection and improved processing of altimetric data sets. We focus on the preprocessing steps of along-track altimetry data from Jason-1 and Envisat to obtain a mean sea surface profile. During this procedure, a rigorous variance propagation is accomplished, so that, for the first time, the full covariance matrix of the mean sea surface is available. The combination of the mean profile and a combined GRACE/GOCE gravity field model yields a mean dynamic topography model for the North Atlantic Ocean that is characterized by a defined set of assumptions. We show that including the geodetically derived mean dynamic topography with the full error structure in a 3D stationary inverse ocean model improves modeled oceanographic features over previous estimates.