Climate spectrum estimation in the presence of timescale errors

We introduce an algorithm (called REDFITmc2) for spectrum estimation in the presence of timescale errors. It is based on the Lomb-Scargle periodogram for unevenly spaced time series, in combination with the Welch's Overlapped Segment Averaging procedure, bootstrap bias correction and persistence estimation. The timescale errors are modelled parametrically and included in the simulations for determining (1) the upper levels of the spectrum of the red-noise AR(1) alternative and (2) the uncertainty of the frequency of a spectral peak. Application of REDFITmc2 to ice core and stalagmite records of palaeoclimate allowed a more realistic evaluation of spectral peaks than when ignoring this source of uncertainty. The results support qualitatively the intuition that stronger effects on the spectrum estimate (decreased detectability and increased frequency uncertainty) occur for higher frequencies. The surplus information brought by algorithm REDFITmc2 is that those effects are quantified. Regarding timescale construction, not only the fixpoints, dating errors and the functional form of the age-depth model play a role. Also the joint distribution of all time points (serial correlation, stratigraphic order) determines spectrum estimation.

Geomorphological and sedimentary evidence for late Pleistocene to Holocene hydrological change along the Río Mamoré, Bolivian Amazon

Aerial photography and satellite imagery reveal manifold geomorphological evidence of a dynamic evolution of past and present rivers in the Bolivian Amazon. Comparison of oxbow lake and meander scar dimensions along an inactive meander belt of the Río Mamoré (Llanos de Moxos, north-eastern Bolivia) and its modern counterpart suggests significant regional paleohydrological variability. We used these features as geomorphological and sedimentary archives to enhance our understanding of longer-term variations of the fluvial system. Late Pleistocene to Holocene hydrological changes of the Río Mamoré are inferred from: (i) the analysis of satellite imagery, (ii) discharge estimates from meander morphology, (iii) stratigraphic, and (iv) chronological information based on luminescence and radiocarbon dating. The combined data from three oxbows indicate that the now abandoned meander belt – the paleo-Mamoré – continued to be active at least until ∼5 ka, and likely even postdating 3 ka. An up to threefold increase in discharge is estimated for the modern Río Mamoré versus the paleo-Mamoré. The altered runoff regime may have triggered an avulsive shift towards the currently active Río Mamoré. The preceding increase in discharge in turn, was possibly related to a shift in climatic conditions, which changed markedly between the mid- and late Holocene in tropical South America. In addition, it may have been the indirect result of capturing the avulsive Río Grande system to the east of the Río Mamoré. Alternative explanations for the differences in dimensions of the paleo versus the modern Río Mamoré, i.e. contemporaneous activity of both rivers or alteration of site factors such as the channel/floodplain relationship, are considered to be unlikely.

Are great Cascadia earthquakes recorded in the sedimentary records from small forearc lakes?

Here we investigate sedimentary records from four small inland lakes located in the southern Cascadia forearc region for evidence of earthquakes. Three of these lakes are in the Klamath Mountains near the Oregon–California border, and one is in the central Oregon Coast range. The sedimentary sequences recovered from these lakes are composed of normal lake sediment interbedded with disturbance event layers. The thickest of these layers are graded, and appear to be turbidites or linked debrites (turbidites with a basal debris-flow deposit), suggesting rapid deposition. Variations in particle size and organic content of these layers are reflected in the density and magnetic susceptibility data. The frequency and timing of these events, based on radiocarbon ages from detrital organics, is similar to the offshore seismogenic turbidite record from trench and slope basin cores along the Cascadia margin. Stratigraphic correlation of these anomalous deposits based on radiocarbon ages, down-core density, and magnetic susceptibility data between lake and offshore records suggest synchronous triggering. The areal extent and multiple depositional environments over which these events appear to correlate suggest that these deposits were most likely caused by shaking during great Cascadia earthquakes.

Chronology of Lake El'gygytgyn sediments – a combined magnetostratigraphic, palaeoclimatic and orbital tuning study based on multi-parameter analyses

A 318-metre-long sedimentary profile drilled by the International Continental Scientific Drilling Program (ICDP) at Site 5011-1 in Lake El’gygytgyn, Far East Russian Arctic, has been analysed for its sedimentologic response to global climate modes by chronostratigraphic methods. The 12 km wide lake is sited off-centre in an 18 km large crater that was created by the impact of a meteorite 3.58 Ma ago. Since then sediments have been continuously deposited. For establishing their chronology, major reversals of the earth’s magnetic field provided initial tie points for the age model, confirming that the impact occurred in the earliest geomagnetic Gauss chron. Various stratigraphic parameters, reflecting redox conditions at the lake floor and climatic conditions in the catchment were tuned synchronously to Northern Hemisphere insolation variations and the marine oxygen isotope stack, respectively. Thus, a robust age model comprising more than 600 tie points could be defined. It could be shown that deposition of sediments in Lake El’gygytgyn occurred in concert with global climatic cycles. The upper �160m of sediments represent the past 3.3 Ma, equivalent to sedimentation rates of 4 to 5 cm ka−1, whereas the lower 160m represent just the first 0.3 Ma after the impact, equivalent to sedimentation rates in the order of 45 cm ka−1. This study also provides orbitally tuned ages for a total of 8 tephras deposited in Lake El’gygytgyn.

A paleoclimate record with tephrochronological age control for the last glacial-interglacial cycle from Lake Ohrid, Albania and Macedonia

Abstract Lake Ohrid is probably of Pliocene age, and the oldest extant lake in Europe. In this study climatic and environmental changes during the last glacial-interglacial cycle are reconstructed using lithological, sedimentological, geochemical and physical proxy analysis of a 15-m-long sediment succession from Lake Ohrid. A chronological framework is derived from tephrochronology and radiocarbon dating, which yields a basal age of ca. 136 ka. The succession is not continuous, however, with a hiatus between ca. 97.6 and 81.7 ka. Sediment accumulation in course of the last climatic cycle is controlled by the complex interaction of a variety of climate-controlled parameters and their impact on catchment dynamics, limnology, and hydrology of the lake. Warm interglacial and cold glacial climate conditions can be clearly distinguished from organic matter, calcite, clastic detritus and lithostratigraphic data. During interglacial periods, short-term fluctuations are recorded by abrupt variations in organic matter and calcite content, indicating climatically-induced changes in lake productivity and hydrology. During glacial periods, high variability in the contents of coarse silt to fine sand sized clastic matter is probably a function of climatically-induced changes in catchment dynamics and wind activity. In some instances tephra layers provide potential stratigraphic markers for short-lived climate perturbations. Given their widespread distribution in sites across the region, tephra analysis has the potential to provide insight into variation in the impact of climate and environmental change across the Mediterranean.

Sediment core fossils in ancient Lake Ohrid: testing for faunal change since the Last Interglacial

Ancient Lake Ohrid is probably of early Pleistocene or Pliocene origin and amongst the few lakes in the world harbouring an outstanding degree of endemic biodiversity. Although there is a long history of evolutionary research in Lake Ohrid, particularly on molluscs, a mollusc fossil record has been missing up to date. For the first time, gastropod and bivalve fossils are reported from the basal, calcareous part of a 2.6 m long sediment succession (core Co1200) from the north-eastern part of Lake Ohrid. Electron spin resonance (ESR) dating of mollusc shells from the same stratigraphic level yielded an age of 130 ± 28 ka. Lithofacies III sediments, i.e. a stratigraphic subdivision comprising the basal succession of core Co1200 between 181.5–263 cm, appeared solid, greyish-white, and consisted almost entirely of silt-sized endogenic calcite (CaCO3>70%) and intact and broken mollusc shells. Here we compare the faunal composition of the thanatocoenosis with recent mollusc associations in Lake Ohrid. A total of 13 mollusc species (9 gastropod and 4 bivalve species) could be identified within Lithofacies III sediments. The value of sediment core fossils for reconstructing palaeoenvironmental settings was evaluated and the agreement between sediment and palaeontological proxies was tested. The study also aims at investigating major faunal changes since the Last Interglacial and searching for signs of extinction events. The combined findings of the ecological study and the sediment characteristics suggest deposition in a shallow water environment during the Last Interglacial. The fossil fauna exclusively included species also found in the present fauna, i.e. no extinction events are evident for this site since the Last Interglacial. The thanatocoenosis showed the highest similarity with recent Intermediate Layer (5–25 m water depth) mollusc assemblages. The demonstrated existence of a mollusc fossil record in Lake Ohrid sediment cores also has great significance for future deep drilling projects. It can be hoped that a more far reaching mollusc fossil record will then be obtained, enabling insight into the early evolutionary history of Lake Ohrid.

OSL and TL dating of the Middle Stone Age sequence at Diepkloof Rock Shelter (South Africa): a clarification

Diepkloof Rock Shelter offers an exceptional opportunity to study the onset and evolution of both Still Bay (SB) and Howiesons Poort (HP) techno-complexes. However, previous age estimates based on luminescence dating of burnt quartzites (Tribolo et al., 2009) and of sediments (Jacobs et al., 2008) were not in agreement. Here, we present new luminescence ages for 17 rock samples (equivalent dose estimated with a SAR-ITL protocol instead of classical MAAD-TL) as well as for 5 sediment samples (equivalent dose estimated with SAR-single grain OSL protocol) and an update of the 22 previous age estimates for burnt lithics (modified calibration and beta dose estimates). While a good agreement between the rock and sediment ages is obtained, these estimates are still significantly older than those reported by Jacobs et al. (2008). After our own analyses of the sediment from Diepkloof, it is suspected that these authors did not correctly chose the parameters for the equivalent dose determination, leading to an underestimate of the equivalent doses, and thus of the ages. From bottom to top, the mean ages are 100 ± 10 ka for stratigraphic unit (SU) Noël and 107 ± 11 ka for SU Mark (uncharacterized Lower MSA), 100 ± 10 ka for SU Lynn-Leo (Pre-SB type Lynn), 109 ± 10 ka for SUs Kim-Larry (SB), 105 ± 10 ka for SUs Kerry-Kate and 109 ± 10 ka for SU Jess (Early HP), 89 ± 8 ka for SU Jude (MSA type Jack), 77 ± 8 ka for SU John, 85 ± 9 ka for SU Fox, 83 ± 8 ka for SU Fred and 65 ± 8 ka for SU OB5 (Intermediate HP), 52 ± 5 ka for SUs OB2-4 (Late HP). This chronology, together with the technological analyses, greatly modifies the current chrono-cultural model regarding the SB and the HP and has important archaeological implications. Indeed, SB and HP no longer appear as short-lived techno-complexes with synchronous appearances for each and restricted to Oxygen Isotopic Stage (OIS) 4 across South Africa, as suggested by Jacobs et al. (2008, 2012). Rather, the sequence of Diepkloof supports a long chronology model with an early appearance of both SB and HP in the first half of OIS 5 and a long duration of the HP into OIS 3. These new dates imply that different technological traditions coexisted during OIS 5 and 4 in southern Africa and that SB and HP can no longer be considered as horizon markers.

Lake sediments as archives of recurrence rates and intensities of past flood events

Palaeoflood hydrology is an expanding field as the damage potential of flood and flood-related processes are increasing with the population density and the value of the infrastructure. Assessing the risk of these hazards in mountainous terrain requires knowledge about the frequency and severness of such events in the past. A wide range of methods is employed using diverse biologic, geomorphic or geologic evidences to track past flood events. Impact of floods are studied and dated on alluvial fans and cones using for example the growth disturbance of trees (Stoffel and Bollschweiler 2008; Schneuwly-Bollschweiler and Stoffel 2012: this volume) or stratigraphic layers deposited by debris flows, allowing to reconstruct past flood frequencies (Bardou et~al. 2003). Further downstream, the classical approach of palaeoflood hydrology (Kochel and Baker 1982) utilizes geomorphic indicators such as overbank sediments, silt lines and erosion features of floods along a river (e.g. Benito and Thorndycraft 2005). Fine-grained sediment settles out of the river suspension in eddies or backwater areas, where the flow velocity of the river is reduced. Records of these deposits at different elevations across a river’s profile can be used to assess the discharge of the past floods. This approach of palaeoflood hydrology studies was successfully applied in several river catchments (e.g. Ely et al. 1993; Macklin and Lewin 2003; O’Connor et al. 1994; Sheffer et al. 2003; Thorndycraft et al. 2005; Thorndycraft and Benito 2006). All these different reconstruction methods have their own advantages and disadvantages, but often these studies have a limited time coverage and the records are potentially incomplete due to lateral limits of depositional areas and due to the erosional power of fluvial processes that remove previously deposited flood witnesses. Here, we present a method that follows the sediment particle transported by a flood event to its final sink: the lacustrine basin.

Identification of erosional mechanisms during past glaciations based on a bedrock surface model of the central European Alps

The bedrock topography beneath the Quaternary cover provides an important archive for the identification of erosional processes during past glaciations. Here, we combined stratigraphic investigations of more than 40,000 boreholes with published data to generate a bedrock topography model for the entire plateau north of the Swiss Alps including the valleys within the mountain belt. We compared the bedrock map with data about the pattern of the erosional resistance of Alpine rocks to identify the controls of the lithologic architecture on the location of overdeepenings. We additionally used the bedrock topography map as a basis to calculate the erosional potential of the Alpine glaciers, which was related to the thickness of the LGM ice. We used these calculations to interpret how glaciers, with support by subglacial meltwater under pressure, might have shaped the bedrock topography of the Alps. We found that the erosional resistance of the bedrock lithology mainly explains where overdeepenings in the Alpine valleys and the plateau occur. In particular, in the Alpine valleys, the locations of overdeepenings largely overlap with areas where the underlying bedrock has a low erosional resistance, or where it was shattered by faults. We also found that the assignment of two end-member scenarios of erosion, related to glacial abrasion/plucking in the Alpine valleys, and dissection by subglacial meltwater in the plateau, may be adequate to explain the pattern of overdeepenings in the Alpine realm. This most likely points to the topographic controls on glacial scouring. In the Alps, the flow of LGM and previous glaciers were constrained by valley flanks, while ice flow was mostly divergent on the plateau where valley borders are absent. We suggest that these differences in landscape conditioning might have contributed to the contrasts in the formation of overdeepenings in the Alpine valleys and the plateau.

On the stratigraphic integrity of leaf-wax biomarkers in loess paleosols

Paleoenvironmental and paleoclimate reconstructions based on molecular proxies, such as those derived from leaf-wax biomarkers, in loess-paleosol sequences represent a promising line of investigation in Quaternary research. The main premise of such reconstructions is the synsedimentary deposition of biomarkers and dust, which has become a debated subject in recent years. This study uses two independent approaches to test the stratigraphic integrity of leaf-wax biomarkers: (i) long-chain n-alkanes and fatty acids are quantified in two sediment-depth profiles in glacial till on the Swiss Plateau, consisting of a Holocene topsoil and the underlying B and C horizons. Since glacial sediments are initially very poor in organic matter, significant amounts of leaf-wax biomarkers in the B and C horizons of those profiles would reflect postsedimentary root-derived or microbial contributions. (ii) Compound-specific radiocarbon measurements are conducted on n-alkanes and n-alkanoic (fatty) acids from several depth intervals in the loess section "Crvenka", Serbia, and the results are compared to independent estimates of sediment age. We find extremely low concentrations of plant-wax n-alkanes and fatty acids in the B and C horizons below the topsoils in the sediment profiles. Moreover, compound-specific radiocarbon analysis yields plant-wax 14C ages that agree well with published luminescence ages and stratigraphy of the Serbian loess deposit. Both approaches confirm that postsedimentary, root-derived or microbial contributions are negligible in the two investigated systems. The good agreement between the ages of odd and even homologues also indicates that reworking and incorporation of fossil leaf waxes is not particularly relevant either.

Geochemical variability in distal and proximal glass from the Youngest Toba Tuff eruption

The Youngest Toba Tuff (YTT, erupted ca. 74 ka ago) is a distinctive and widespread tephra marker across south and southeast Asia. The climatic, human and environmental consequences of the YTT eruption are widely debated. Although a considerable body of geochemical data is available for this unit, there has not been a systematic study of the variability of the ash geochemistry. Intrinsic (magmatic) and extrinsic (post-depositional) chemical variations bring fundamental information regarding the petrogenesis of the magma, the distribution of the tephra and the interaction between the ash and the receiving environment. Considering the importance of the geochemistry of the YTT for stratigraphic correlations and eruptive models, it is central to the YTT debate to quantify and interpret such variations. Here we collate all published geochemical data on the YTT glass, including analyses from 68 sites described in the literature and three new samples. Two principal sources of chemical variation are investigated: (i) compositional zonation of the magma reservoir, and (ii) post-depositional alteration. Post-depositional leaching is responsible for up to ca. 11% differences in Na2O/K2O and ca. 1% differences in SiO2/Al2O3 ratios in YTT glass from marine sites. Continental tephra are 2% higher in Na2O/K2O and 3% higher in SiO2/Al2O3 respect to the marine tephra. We interpret such post-depositional glass alteration as related to seawater induced alkali migration in marine environments, or to site-specific water pH. Crystal fractionation and consequential magmatic differentiation, which produced order-of-magnitude variations in trace element concentrations reported in the literature, also produced major element differences in the YTT glass. FeO/Al2O3 ratios vary by about 50 %, which is analytically significant. These variations represent magmatic fractionation involving Fe-bearing phases. We also compared major element concentrations in YTT and Oldest Toba Tuff (OTT) ash samples, to identify potential compositional differences that could constrain the stratigraphic identity of the Morgaon ash (Western India); no differences between the OTT and YTT samples were observed.

A stratigraphic framework for abrupt climatic changes during the Last Glacial period based on three synchronized Greenland ice-core records: refining and extending the INTIMATE event stratigraphy

Due to their outstanding resolution and well-constrained chronologies, Greenland ice-core records provide a master record of past climatic changes throughout the Last Interglacial–Glacial cycle in the North Atlantic region. As part of the INTIMATE (INTegration of Ice-core, MArine and TErrestrial records) project, protocols have been proposed to ensure consistent and robust correlation between different records of past climate. A key element of these protocols has been the formal definition and ordinal numbering of the sequence of Greenland Stadials (GS) and Greenland Interstadials (GI) within the most recent glacial period. The GS and GI periods are the Greenland expressions of the characteristic Dansgaard–Oeschger events that represent cold and warm phases of the North Atlantic region, respectively. We present here a more detailed and extended GS/GI template for the whole of the Last Glacial period. It is based on a synchronization of the NGRIP, GRIP, and GISP2 ice-core records that allows the parallel analysis of all three records on a common time scale. The boundaries of the GS and GI periods are defined based on a combination of stable-oxygen isotope ratios of the ice (δ18O, reflecting mainly local temperature) and calcium ion concentrations (reflecting mainly atmospheric dust loading) measured in the ice. The data not only resolve the well-known sequence of Dansgaard–Oeschger events that were first defined and numbered in the ice-core records more than two decades ago, but also better resolve a number of short-lived climatic oscillations, some defined here for the first time. Using this revised scheme, we propose a consistent approach for discriminating and naming all the significant abrupt climatic events of the Last Glacial period that are represented in the Greenland ice records. The final product constitutes an extended and better resolved Greenland stratotype sequence, against which other proxy records can be compared and correlated. It also provides a more secure basis for investigating the dynamics and fundamental causes of these climatic perturbations.

“Unconformity-Bounded” Stratigraphic Units in the South Polar Layered Deposits (Promethei Lingula, Mars)

We conducted a stratigraphic analysis of the South Polar Layered Deposits (SPLDs) in Promethei Lingula (PL, Mars) based on the identification of regional unconformities at visible and radar wavelengths. According to the terrestrial classification, this approach constrains the stratigraphy of the region and remedies the ambiguous interpretation of stratigraphy through marker layers, bypassing the problem related to the morphologic and radiometric appearance of the layers. Thus, the approach does not exclude diverse classifications, but complements them, whereas other discriminant elements are doubtful or difficult/impossible to be defined. Using this approach, we defined two stratigraphic units (or synthems: PL1 and PL2) in PL, which are morphologically different and divided by a regional unconformity (AuR1). This stratigraphic architecture implies that the geological history of PL has been conditioned by periodic changes in climate, which in turn are related to orbital variations of Mars.