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.

Bias associated with detrital zircon geochronology and thermochronometry

Detrital studies that utilize zircon U–Pb geochronology and fission-track (FT) thermochronometry are subject to a range of potential sources of bias that should be properly evaluated and minimized. Some of them are common to any single-grain mineral analysis (e.g., variable bedrock mineral fertility, hydraulic sorting during transport, selective grain loss during sample processing), whereas others are intrinsic to zircon, and are related to radiation damage and age discordance. In this article, we quantify the impact of intrinsic bias on detrital studies thanks to the analysis of modern detritus shed from the European Alps, and illustrate the general implications on geological interpretations. We show that detrital zircon U–Pb age distributions based on statistically robust datasets are highly reproducible and representative of the parent bedrock ages in the catchment. Arbitrary or selective removal of discordant grain ages can be minimized by using the Kolmogorov–Smirnov test to identify an appropriate cutoff level. Loss of metamict (α-damaged) zircon has a minor impact on data representativeness, and is mainly controlled by regional metamorphism rather than by mechanical abrasion during river transport. Zircon FT grain-age distributions were found to have poor reproducibility, although age spectra are consistent with bedrock data. However, unlike the U–Pb datasets, U-rich zircon grains (> 1000 ppm) are systematically missed, and undatable grains may exceed 50%. We identify two major sources of distribution bias specific to zircon FT datasets: (i) sediment sources dominated by U-rich zircon grains are markedly underrepresented in the detrital record, because such grains often have uncountable high densities of fission tracks (“U concentration bias”); (ii) sediment sources that shed zircon grains with high levels of α-damage are underrepresented, because these grains are lost during chemical etching for FT revelation (“etching bias”). In the case of multimethod dating on the same grains (e.g., FT and U–Pb double dating), bias affecting detrital zircon FT dating propagates to the entire dataset. These effects may not impact on exhumation-rate studies that utilize the youngest grain ages (i.e., lag-time approach). However, they represent a limiting factor for conventional provenance studies, and generally preclude application of zircon FT dating to sediment budget calculations.

Vegetation and environmental changes in Northern Anatolia between 134 and 119 ka recorded in Black Sea sediments

This multiproxy study on SE Black Sea sediments provides the first detailed reconstruction of vegetation and environmental history of Northern Anatolia between 134 and 119 ka. Here, the glacial–interglacial transition is characterized by several short-lived alternating cold and warm events preceding a meltwater pulse (~ 130.4–131.7 ka). The latter is reconstructed as a cold arid period correlated to Heinrich event 11. The initial warming is evidenced at ~ 130.4 ka by increased primary productivity in the Black Sea, disappearance of ice-rafted detritus, and spreading of oaks in Anatolia. A Younger Dryas-type event is not identifiable. The Eemian vegetation succession corresponds to the main climatic phases in Europe: i) the Quercus–Juniperus phase (128.7–126.4 ka) indicates a dry continental climate; ii) the Ostrya–Corylus–Quercus–Carpinus phase (126.4–122.9 ka) suggests warm summers, mild winters, and high year-round precipitation; iii) the Fagus–Carpinus phase (122.9–119.5 ka) indicates cooling and high precipitation; and iv) increasing Pinus at ~ 121 ka marks the onset of cooler/drier conditions. Generally, pollen reconstructions suggest altitudinal/latitudinal migrations of vegetation belts in Northern Anatolia during the Eemian caused by increased transport of moisture. The evidence for the wide distribution of Fagus around the Black Sea contrasts with the European records and is likely related to climatic and genetic factors.

The Value of Adding Optics to Ecosystem Models: A Case Study

Many ecosystem models have been developed to study the ocean's biogeochemical properties, but most of these models use simple formulations to describe light penetration and spectral quality. Here, an optical model is coupled with a previously published ecosystem model that explicitly represents two phytoplankton (picoplankton and diatoms) and two zooplankton functional groups, as well as multiple nutrients and detritus. Surface ocean color fields and subsurface light fields are calculated by coupling the ecosystem model with an optical model that relates biogeochemical standing stocks with inherent optical properties (absorption, scattering); this provides input to a commercially available radiative transfer model (Ecolight). We apply this bio-optical model to the equatorial Pacific upwelling region, and find the model to be capable of reproducing many measured optical properties and key biogeochemical processes in this region. Our model results suggest that non-algal particles largely contribute to the total scattering or attenuation (> 50% at 660 nm) but have a much smaller contribution to particulate absorption (< 20% at 440 nm), while picoplankton dominate the total phytoplankton absorption (> 95% at 440 nm). These results are consistent with the field observations. In order to achieve such good agreement between data and model results, however, key model parameters, for which no field data are available, have to be constrained. Sensitivity analysis of the model results to optical parameters reveals a significant role played by colored dissolved organic matter through its influence on the quantity and quality of the ambient light. Coupling explicit optics to an ecosystem model provides advantages in generating: (1) a more accurate subsurface light-field, which is important for light sensitive biogeochemical processes such as photosynthesis and photo-oxidation, (2) additional constraints on model parameters that help to reduce uncertainties in ecosystem model simulations, and (3) model output which is comparable to basic remotely-sensed properties. In addition, the coupling of biogeochemical models and optics paves the road for future assimilation of ocean color and in-situ measured optical properties into the models.

Mycobacterium ulcerans persistence at a village water source of Buruli ulcer patients

Buruli ulcer (BU), a neglected tropical disease of the skin and subcutaneous tissue, is caused by Mycobacterium ulcerans and is the third most common mycobacterial disease after tuberculosis and leprosy. While there is a strong association of the occurrence of the disease with stagnant or slow flowing water bodies, the exact mode of transmission of BU is not clear. M. ulcerans has emerged from the environmental fish pathogen M. marinum by acquisition of a virulence plasmid encoding the enzymes required for the production of the cytotoxic macrolide toxin mycolactone, which is a key factor in the pathogenesis of BU. Comparative genomic studies have further shown extensive pseudogene formation and downsizing of the M. ulcerans genome, indicative for an adaptation to a more stable ecological niche. This has raised the question whether this pathogen is still present in water-associated environmental reservoirs. Here we show persistence of M. ulcerans specific DNA sequences over a period of more than two years at a water contact location of BU patients in an endemic village of Cameroon. At defined positions in a shallow water hole used by the villagers for washing and bathing, detritus remained consistently positive for M. ulcerans DNA. The observed mean real-time PCR Ct difference of 1.45 between the insertion sequences IS2606 and IS2404 indicated that lineage 3 M. ulcerans, which cause human disease, persisted in this environment after successful treatment of all local patients. Underwater decaying organic matter may therefore represent a reservoir of M. ulcerans for direct infection of skin lesions or vector-associated transmission.

Decadal Variability in the Northeast Pacific in a Physical-Ecosystem Model: Role of Mixed Layer Depth and Trophic Interactions

A basin-wide interdecadal change in both the physical state and the ecology of the North Pacific occurred near the end of 1976. Here we use a physical-ecosystem model to examine whether changes in the physical environment associated with the 1976-1977 transition influenced the lower trophic levels of the food web and if so by what means. The physical component is an ocean general circulation model, while the biological component contains 10 compartments: two phytoplankton, two zooplankton, two detritus pools, nitrate, ammonium, silicate, and carbon dioxide. The model is forced with observed atmospheric fields during 1960-1999. During spring, there is a similar to 40% reduction in plankton biomass in all four plankton groups during 1977-1988 relative to 1970-1976 in the central Gulf of Alaska (GOA). The epoch difference in plankton appears to be controlled by the mixed layer depth. Enhanced Ekman pumping after 1976 caused the halocline to shoal, and thus the mixed layer depth, which extends to the top of the halocline in late winter, did not penetrate as deep in the central GOA. As a result, more phytoplankton remained in the euphotic zone, and phytoplankton biomass began to increase earlier in the year after the 1976 transition. Zooplankton biomass also increased, but then grazing pressure led to a strong decrease in phytoplankton by April followed by a drop in zooplankton by May: Essentially, the mean seasonal cycle of plankton biomass was shifted earlier in the year. As the seasonal cycle progressed, the difference in plankton concentrations between epochs reversed sign again, leading to slightly greater zooplankton biomass during summer in the later epoch.

Lability of Proteinaceous Material in Estuarine Seston and Subcellular Fractions of Phytoplankton

The rate of proteolysis of amino acids was used to assess the nutritional lability of various materials making up estuarine seston in 3 Maine, USA, estuaries. Physical separations of subcellular fractions of phytoplankton cells led to higher proteolysis rate constants for the cytoplasmic fraction (>1.2 h(-1)) than for the membrane fraction (0.2 to 1 h(-1)). Whole cells, copepod fecal pellets, bottom sediments, and estuarine seston had overlapping ranges of rate constants of 0.17 to 1.3 h(-1), which were indistinguishable from one another. Protein pools in the seston of these estuaries throughout the seasons were dominated by phytoplankton production and its fresh detrital products. Inverse relationships between proteolysis rate constants for estuarine seston and the ratios of pheopigments to chlorophyll indicates that the average lability of seston decreases with the disappearance of cytoplasmic material in suspension. This kinetic approach to the quality of food resources implies the existence of different pools of digestible protein for estuarine heterotrophs with different gut residence times. Preferential enrichment of membrane components in sestonic detritus may result from the differential lability of proteins in cytoplasm versus membrane components of cells.

Provenance versus weathering control on the composition of tropical river mud (southern Africa)

This study presents an integrated mineralogical-geochemical data base on fine-grained sediments transported by all major rivers of southern Africa, including the Zambezi, Okavango, Limpopo, Olifants, Orange and Kunene. Clay mineralogy, bulk geochemistry, Sr and Nd isotopic signatures of river mud, considered as proxy of suspended load, are used to investigate the influence of source-rock lithology and weathering intensity on the composition of clay and silt produced in subequatorial to subtropical latitudes. Depletion in mobile alkali and alkaline-earth metals, minor in arid Namibia, is strong in the Okavango, Kwando and Upper Zambezi catchments, where recycling is also extensive. Element removal is most significant for Na, and to a lesser extent for Sr. Depletion in K, Ca and other elements, negligible in Namibia, is moderate elsewhere. The most widespread clay minerals are smectite, dominant in muds derived from Karoo or Etendeka flood basalts, or illite and chlorite, dominant in muds derived from metasedimentary rocks of the Damara Orogen or Zimbabwe Craton. Kaolinite represents 30-40% of clay minerals only in Okavango and Upper Zambezi sediments sourced in humid subequatorial Angola and Zambia. After subtracting the effects of recycling and of local accumulation of authigenic carbonates in soils, the regional distribution of clay minerals and chemical indices consistently reflect weathering intensity primarily controlled by climate. Bulk geochemistry identifies most clearly volcaniclastic sediments and mafic sources in general, but cannot discriminate the other sources of detritus in detail. Instead, Sr and Nd isotopic fingerprints are insensitive to weathering, and thus mirror faithfully the tectonic structure of the southern African continent. Isotopic tools thus represent a much firmer basis than bulk geochemistry or clay mineralogy in the provenance study of mudrocks.

Slab rollback orogeny in the Alps and evolution of the Swiss Molasse basin

The stratigraphies of foreland basins have been related to orogeny, where continent–continent collision causes the construction of topography and the downwarping of the foreland plate. These mechanisms have been inferred for the Molasse basin, stretching along the northern margin of the European Alps. Continuous flexural bending of the subducting European lithosphere as a consequence of topographic loads alone would imply that the Alpine topography would have increased at least between 30 Ma and ca. 5–10 Ma when the basin accumulated the erosional detritus. This, however, is neither consistent with observations nor with isostatic mass balancing models because paleoaltimetry estimates suggest that the topography has not increased since 20 Ma. Here we show that a rollback mechanism for the European plate is capable of explaining the construction of thick sedimentary successions in the Molasse foreland basin where the extra slab load has maintained the Alpine surface at low, but constant, elevations.

Sand geochemistry on sediment core CRP-3 from the Ross Sea, Antarctica

A total of 167 samples distubuted throughout the CRP-3 drillhole from 5.77 to 787.68 mbsf and representing fine to coarse sandstones have been analysed by X-ray fluorescence spectrometry (XRF) Bulk sample geochemistry (major and trace elements) indicates a dominant provenance of detritus from the Ferrar Supergroup in the uppermost 200 mbsf of the core. A markedly increased contribution from the Beacon sandstones is recognized below 200 mbsf and down to 600 mbsf. In the lower part of CRP-3, down to 787.68 mbsf, geochemical evidence for influxes of Ferrar materials is again recorded. On the basis of preliminary magnetostratigraphic data reported for the lower 447 mbsf of the drillhole, we tentatively evaluated the main periodicities modulating the geochemical records. Our results identify a possible influence of the precession, obliquity and long-eccentricity astronomical components (21, 41, and 400 ky frequency bands) on the deposition mechanisms of the studied glaciomarine sediments.

Geochemical investigations of elements in sediment core CRP-1 from the Ross Sea, Antarctica

Geochemical data are presented for samples from strata, mainly of Miocene age, in the Cape Roberts-1 core (western McMurdo Sound, Antarctica) to assess the sediment provenance. Bulk (major and trace element) chemistry together with bulk mineralogy of fine-grained sandstones, siltstones, mudstones, and diamictites indicate that chemical alteration of source materials, fractionation due to sedimentary sorting, and diagenetic effects were not significant in the Cape Roberts sediment history. Relevant geochemical parameters are consistent with the Cape Roberts sediments being derived mainly from the crystalline basement and the Beacon Supergroup. On the basis of element distributions, an additional contribution from the Ferrar Dolerite and, mainly above about 60 m, influxes of detritus derived from basanitic to intermediate members of the McMurdo Volcanic Group are recognised.

Oceanography during TYDEMAN cruise DCM (DeepChlorMax) to the Equatorial Atlantic

The cruise with RV Tydeman was devoted to study permanently stratified plankton systems in the (sub)tropical ocean, which are characterised by a deep chlorophyll peak between 80 and 150 m. To minimise lateral effects by horizontal transport of nutrients and organic matter from river outflow and upwelling regions, stations were selected in the middle of the North Atlantic Ocean between the continents of America and Africa. (5 - 35° N and 50 - 15° W). Here the vertical distributions of light and nutrients control the abundance and growth of autotrophic algae in the thermically stratified water column. This phytoplankton is numerically dominated by the prokaryotic picoplankters Synechococcus spp. and Prochlorococcus spp., which are smaller than 2 ?m. The productivity of the 100 to 150 m deep euphotic zone can be high, because a high heterotrophic/autotrophic biomass ratio induces a rapid regeneration of nutrients and inorganic carbon. Primary grazers are mainly micro-organisms such as heterotrophic nannoflagellates and ciliates, which feed on the small algae and on bacteria. Heterotrophic bacteria can outnumber the autotrophic algae, because their number is related to the substrate pools of dissolved and particulate dead organic matter. These DOC and detritus pools reach equilibrium at a concentration, where the rate of their production (proportional to algal biomass) equals their mineralisation and sinking rate (proportional to the concentration and weight of POC and detritus). At a relatively low value of the weight-specific loss rates, the equilibrium concentration of these carbon pools and their load of bacteria can be high. The bacterial productivity is proportional to the mineralisation rate, which in a steady state can never be higher than the rate of primary production. Hence the ratio in turnover rate of bacteria and autotrophs tends to be reciprocally proportional to their biomass ratio.

Mineralogy and age determination of sediment core D11957P

The Tore Seamount is a circular, volcano-like feature 100 km in diameter with its summit at 2200 m water depth and a small, 5000 m deep basin in its interior. It is situated approximately 300 km west of Lisbon and is surrounded by deep abyssal plains. This site with a standard pelagic stratigraphy is the southernmost point where the so-called Heinrich events have so far been recorded. A succession of alternating interglacial/glacial periods reveals a stratigraphic record back to the beginning of isotopic stage 7 (225 kyr). Climatic changes are identifiable by coherent variations in colour, carbonate content and distribution of ice-rafted detritus in the carbonate-free fraction. Inputs of ice-rafted quartz are well defined. Characteristics in common with other sites showing Heinrich layers include a high terrigenous to biogenic ratio, a dramatic decrease in the accumulation rate of foraminifera shells, an increase in dolomite abundance and the occurrence of polar foraminiferal species indicating southwards penetration of cold waters which lead us to consider a wider southeastern extent of the North Atlantic ice-rafted detritus belt than hitherto. If the presently accepted position of the Polar Front is maintained, icebergs must have been swept southwards from the southern boundary of the pack ice in a current merging into the ancestral Canary Current, bringing ice-rafted material to the Tore Seamount. The coincidence of reddish-feldspar, probably derived from the northern Appalachian Triassic red facies, with the transparent quartz suggests at least a partial Labrador source for all the Heinrich layers here, including HL 3. In comparison to other sites in the entire North Atlantic, two exceptions stand out: the absence of HL 5 and the low detritus to biogenics ratio for HL 3. The simultaneous occurrence of these two types of ice-rafted minerals is a new piece in the puzzle of the origin of Heinrich layers.

(Table 1) Sediment descriptions and age at DSDP Legs 5 and 63 Holes

Sediments from immediately above basalt basement and from between sections of basalt recovered from Deep Sea Drilling Project Legs 5 and 63 were analyzed by atomic absorption spectroscopy for Mg, Al, Si, Ca, Mn, Fe, Co, Ni, Cu, Zn, and Ba. All of these sediments showed enrichment in Fe and Mn over values typical of detritus supplied to the northeastern Pacific Ocean. X-ray diffractometry and differential chemical leaching indicate that up to 50% of the sediment, by weight, is in amorphous phases and that these phases are rich in Mn, Co, Cu, Ni, and Zn. Multivariate statistical analysis and normative partitioning of the chemical data indicate that much of the excess Fe and other transition elements in the sediment originate from hydrothermal sources.

Carbonate content in lower Cretaceous sediments of ODP Hole 123-765C

Sediments recovered during Ocean Drilling Program (ODP) Leg 123 from the Argo Abyssal Plain (AAP) consist largely of turbidites derived from the adjacent Australian continental margin. The oldest abundant turbidites are Valanginian-Aptian in age and have a mixed (smarl) composition; they contain subequal amounts of calcareous and siliceous biogenic components, as well as clay and lesser quartz. Most are thin-bedded, fine sand- to mud-sized, and best described by Stow and Piper's model (1984) for fine-grained biogenic turbidites. Thicker (to 3 m), coarser-grained (medium-to-coarse sand-sized) turbidites fit Bouma's model (1962) for sandy turbidites; these generally are base-cut-out (BCDE, BDE) sequences, with B-division parallel lamination as the dominant structure. Parallel laminae most commonly concentrate quartz and/or calcispheres vs. lithic clasts or clay, but distinctive millimeter- to centimeter-thick, radiolarian-rich laminae occur in both fine- and coarse-grained Valanginian-Hauterivian turbidites. AAP turbidites were derived from relatively deep parts of the continental margin (outer shelf, slope, or rise) that lay below the photic zone, but above the calcite compensation depth (CCD). Biogenic components are largely pelagic (calcispheres, foraminifers, radiolarians, nannofossils); lesser benthic foraminifers are characteristic of deep-water (abyssal to bathyal) environments. Abundant nonbiogenic components are mostly clay and clay clasts; smectite is the dominant clay species, and indicates a volcanogenic provenance, most likely the Triassic-Jurassic volcanic suite exposed along the northern Exmouth Plateau. Lower Cretaceous smarl turbidites were generated during eustatic lowstands and may have reached the abyssal plain via Swan Canyon, a submarine canyon thought to have formed during the Late Jurassic. In contrast to younger AAP turbidites, however, Lower Cretaceous turbidites are relatively fine-grained and do not contain notably older reworked fossils. Early in its history, the northwest Australian margin provided mainly contemporaneous slope sediment to the AAP; marginal basins adjacent to the continent trapped most terrigenous detritus, and pronounced canyon incisement did not occur until Late Cretaceous and, especially, Cenozoic time.