Stratigraphic analysis of lake level fluctuations in Lake Ohrid: an integration of high resolution hydro-acoustic data and sediment cores

Abstract. Ancient Lake Ohrid is a steep-sided, oligotrophic, karst lake that was tectonically formed most likely within the Pliocene and often referred to as a hotspot of endemic biodiversity. This study aims on tracing significant lake level fluctuations at Lake Ohrid using high-resolution acoustic data in combination with lithological, geochemical, and chronological information from two sediment cores recovered from sub-aquatic terrace levels at ca. 32 and 60m water depth. According to our data, significant lake level fluctuations with prominent lowstands of ca. 60 and 35m below the present water level occurred during Marine Isotope Stage (MIS) 6 and MIS 5, respectively. The effect of these lowstands on biodiversity in most coastal parts of the lake is negligible, due to only small changes in lake surface area, coastline, and habitat. In contrast, biodiversity in shallower areas was more severely affected due to disconnection of today sublacustrine springs from the main water body. Multichannel seismic data from deeper parts of the lake clearly image several clinoform structures stacked on top of each other. These stacked clinoforms indicate significantly lower lake levels prior to MIS 6 and a stepwise rise of water level with intermittent stillstands since its existence as water-filled body, which might have caused enhanced expansion of endemic species within Lake Ohrid.

High-resolution sequence stratigraphic correlation in the Upper Jurassic (Kimmeridgian)–Upper Cretaceous (Cenomanian) peritidal carbonate deposits (Western Taurides, Turkey)

Upper Jurassic (Kimmeridgian)±Upper Cretaceous (Cenomanian) inner platform carbonates in the Western Taurides are composed of metre-scale upward-shallowing cyclic deposits (parasequences) and important karstic surfaces capping some of the cycles. Peritidal cycles (shallow subtidal facies capped by tidal-¯at laminites or fenestrate limestones) are regressive- and transgressive-prone (upward-deepening followed by upward-shallowing facies trends). Subtidal cycles are of two types and indicate incomplete shallowing. Submerged subtidal cycles are composed of deeper subtidal facies overlain by shallow subtidal facies. Exposed subtidal cycles consist of deeper subtidal facies overlain by shallow subtidal facies that are capped by features indicative of prolonged subaerial exposure. Subtidal facies occur characteristically in the Jurassic, while peritidal cycles are typical for the Lower Cretaceous of the region. Within the foraminiferal and dasyclad algal biostratigraphic framework, four karst breccia levels are recognized as the boundaries of major second-order cycles, introduced for the ®rst time in this study. These levels correspond to the Kimmeridgian±Portlandian boundary, mid-Early Valanginian, mid-Early Aptian and mid-Cenomanian and represent important sea level falls which affected the distribution of foraminiferal fauna and dasyclad ¯ora of the Taurus carbonate platform. Within the Kimmeridgian±Cenomanian interval 26 third-order sequences (types 1 and 2) are recognized. These sequences are the records of eustatic sea level ¯uctuations rather than the records of local tectonic events because the boundaries of the sequences representing 1±4 Ma intervals are correlative with global sea level falls. Third-order sequences and metre-scale cyclic deposits are the major units used for long-distance, high-resolution sequence stratigraphic correlation in the Western Taurides. Metre-scale cyclic deposits (parasequences) in the Cretaceous show genetical stacking patterns within third-order sequences and correspond to fourth-order sequences representing 100±200 ka. These cycles are possibly the E2 signal (126 ka) of the orbital eccentricity cycles of the Milankovitch band. The slight deviation of values, calculated for parasequences, from the mean value of eccentricity cycles can be explained by the currently imprecise geochronology established in the Cretaceous and missed sea level oscillations when the platform lay above fluctuating sea level.

Late Holocene air temperature variability reconstructed from the sediments of Laguna Escondida, Patagonia, Chile (45 degrees 30 ' S)

Climate and environmental reconstructions from natural archives are important for the interpretation of current climatic change. Few quantitative high-resolution reconstructions exist for South America which is the only land mass extending from the tropics to the southern high latitudes at 56°S. We analyzed sediment cores from two adjacent lakes in Northern Chilean Patagonia, Lago Castor (45°36′S, 71°47′W) and Laguna Escondida (45°31′S, 71°49′W). Radiometric dating (210Pb, 137Cs, 14C-AMS) suggests that the cores reach back to c. 900 BC (Laguna Escondida) and c. 1900 BC (Lago Castor). Both lakes show similarities and reproducibility in sedimentation rate changes and tephra layer deposition. We found eight macroscopic tephras (0.2–5.5 cm thick) dated at 1950 BC, 1700 BC, at 300 BC, 50 BC, 90 AD, 160 AD, 400 AD and at 900 AD. These can be used as regional time-synchronous stratigraphic markers. The two thickest tephras represent known well-dated explosive eruptions of Hudson volcano around 1950 and 300 BC. Biogenic silica flux revealed in both lakes a climate signal and correlation with annual temperature reanalysis data (calibration 1900–2006 AD; Lago Castor r = 0.37; Laguna Escondida r = 0.42, seven years filtered data). We used a linear inverse regression plus scaling model for calibration and leave-one-out cross-validation (RMSEv = 0.56 °C) to reconstruct sub decadal-scale temperature variability for Laguna Escondida back to AD 400. The lower part of the core from Laguna Escondida prior to AD 400 and the core of Lago Castor are strongly influenced by primary and secondary tephras and, therefore, not used for the temperature reconstruction. The temperature reconstruction from Laguna Escondida shows cold conditions in the 5th century (relative to the 20th century mean), warmer temperatures from AD 600 to AD 1150 and colder temperatures from AD 1200 to AD 1450. From AD 1450 to AD 1700 our reconstruction shows a period with stronger variability and on average higher values than the 20th century mean. Until AD 1900 the temperature values decrease but stay slightly above the 20th century mean. Most of the centennial-scale features are reproduced in the few other natural climate archives in the region. The early onset of cool conditions from c. AD 1200 onward seems to be confirmed for this region.

Martian sub-surface ionising radiation: biosignature and geology

The surface of Mars, unshielded by thick atmosphere or global magnetic field, is exposed to high levels of cosmic radiation. This ionising radiation field is deleterious to the survival of dormant cells or spores and the persistence of molecular biomarkers in the subsurface, and so its characterisation is of prime astrobiological interest. Here, we present modelling results of the absorbed radiation dose as a function of depth through the Martian subsurface, suitable for calculation of biomarker persistence. A second major implementation of this dose accumulation rate data is in application of the optically stimulated luminescence technique for dating Martian sediments. We present calculations of the dose-depth profile in the Martian subsurface for various scenarios: variations of surface composition (dry regolith, ice, layered permafrost), solar minimum and maximum conditions, locations of different elevation (Olympus Mons, Hellas basin, datum altitude), and increasing atmospheric thickness over geological history. We also model the changing composition of the subsurface radiation field with depth compared between Martian locations with different shielding material, determine the relative dose contributions from primaries of different energies, and discuss particle deflection by the crustal magnetic fields.

Surface uplift and climate change - the geomorphic evolution of the western escarpment of the Andes of northern Chile between the Miocene and present

The Western Escarpment of the Andes at 18.30°S (Arica area, northern Chile) is a classical example for a transient state in landscape evolution. This part of the Andes is characterized by the presence of >10,000 km2 plains that formed between the Miocene and the present, and >1500 m deeply incised valleys. Although processes in these valleys scale the rates of landscape evolution, determinations of ages of incision, and more importantly, interpretations of possible controls on valley formation have been controversial. This paper uses morphometric data and observations, stratigraphic information, and estimates of sediment yields for the time interval between ca. 7.5 Ma and present to illustrate that the formation of these valleys was driven by two probably unrelated components. The first component is a phase of base-level lowering with magnitudes of∼300–500 m in the Coastal Cordillera. This period of base-level change in the Arica area, that started at ca. 7.5 Ma according to stratigraphic data, caused the trunk streams to dissect headward into the plains. The headward erosion interpretation is based on the presence of well-defined knickzones in stream profiles and the decrease in valley widths from the coast toward these knickzones. The second component is a change in paleoclimate. This interpretation is based on (1) the increase in the size of the largest alluvial boulders (from dm to m scale) with distal sources during the last 7.5 m.y., and (2) the calculated increase in minimum fluvial incision rates of ∼0.2 mm/yr between ca. 7.5 Ma and 3 Ma to ∼0.3 mm/yr subsequently. These trends suggest an increase in effective water discharge for systems sourced in the Western Cordillera (distal source). During the same time, however, valleys with headwaters in the coastal region (local source) lack any evidence of fluvial incision. This implies that the Coastal Cordillera became hyperarid sometime after 7.5 Ma. Furthermore, between 7.5 Ma and present, the sediment yields have been consistently higher in the catchments with distal sources (∼15 m/m.y.) than in the headwaters of rivers with local sources (<7 m/m.y.). The positive correlation between sediment yields and the altitude of the headwaters (distal versus local sources) seems to reflect the effect of orographic precipitation on surface erosion. It appears that base-level change in the coastal region, in combination with an increase in the orographic effect of precipitation, has controlled the topographic evolution of the northern Chilean Andes.

Age modeling of young non-varved lake sediments: methods and limits. Examples from two lakes in Central Chile

High-resolution and highly precise age models for recent lake sediments (last 100–150 years) are essential for quantitative paleoclimate research. These are particularly important for sedimentological and geochemical proxies, where transfer functions cannot be established and calibration must be based upon the relation of sedimentary records to instrumental data. High-precision dating for the calibration period is most critical as it determines directly the quality of the calibration statistics. Here, as an example, we compare radionuclide age models obtained on two high-elevation glacial lakes in the Central Chilean Andes (Laguna Negra: 33°38′S/70°08′W, 2,680 m a.s.l. and Laguna El Ocho: 34°02′S/70°19′W, 3,250 m a.s.l.). We show the different numerical models that produce accurate age-depth chronologies based on 210Pb profiles, and we explain how to obtain reduced age-error bars at the bottom part of the profiles, i.e., typically around the end of the 19th century. In order to constrain the age models, we propose a method with five steps: (i) sampling at irregularly-spaced intervals for 226Ra, 210Pb and 137Cs depending on the stratigraphy and microfacies, (ii) a systematic comparison of numerical models for the calculation of 210Pb-based age models: constant flux constant sedimentation (CFCS), constant initial concentration (CIC), constant rate of supply (CRS) and sediment isotope tomography (SIT), (iii) numerical constraining of the CRS and SIT models with the 137Cs chronomarker of AD 1964 and, (iv) step-wise cross-validation with independent diagnostic environmental stratigraphic markers of known age (e.g., volcanic ash layer, historical flood and earthquakes). In both examples, we also use airborne pollutants such as spheroidal carbonaceous particles (reflecting the history of fossil fuel emissions), excess atmospheric Cu deposition (reflecting the production history of a large local Cu mine), and turbidites related to historical earthquakes. Our results show that the SIT model constrained with the 137Cs AD 1964 peak performs best over the entire chronological profile (last 100–150 years) and yields the smallest standard deviations for the sediment ages. Such precision is critical for the calibration statistics, and ultimately, for the quality of the quantitative paleoclimate reconstruction. The systematic comparison of CRS and SIT models also helps to validate the robustness of the chronologies in different sections of the profile. Although surprisingly poorly known and under-explored in paleolimnological research, the SIT model has a great potential in paleoclimatological reconstructions based on lake sediments

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.