Late Quaternary palaeoenvironmental reconstruction from Lakes Ohrid and Prespa (Macedonia/Albania border) using stable isotopes

Abstract. Here we present stable isotope data from three sediment records from lakes that lie along the Macedonian- Albanian border (Lake Prespa: 1 core, and Lake Ohrid: 2 cores). The records only overlap for the last 40 kyr, although the longest record contains the MIS 5/6 transition (Lake Ohrid). The sedimentary characteristics of both lakes differ significantly between the glacial and interglacial phases. At the end of MIS 6 Lake Ohrid’s water level was low (high �18Ocalcite) and, although productivity was increasing (high calcite content), the carbon supply was mainly from inorganic catchment rock sources (high �13Ccarb). During the last interglacial, calcite and TOC production and preservation increased, progressively lower �18Ocalcite suggest increase in humidity and lake levels until around 115 ka. During ca. 80 ka to 11 ka the lake records suggest cold conditions as indicated by negligible calcite precipitation and low organic matter content. In Lake Ohrid, �13Corg are complacent; in contrast, Lake Prespa shows consistently higher �13Corg suggesting a low oxidation of 13C-depleted organic matter in agreement with a general deterioration of climate conditions during the glacial. From 15 ka to the onset of the Holocene, calcite and TOC begin to increase, suggesting lake levels were probably low (high �18Ocalcite). In the Holocene (11 ka to present) enhanced productivity is manifested by high calcite and organic matter content. All three cores show an early Holocene characterised by low �18Ocalcite, apart from the very early Holocene phase in Prespa where the lowest �18Ocalcite occurs at ca. 7.5 ka, suggesting a phase of higher lake level only in (the more sensitive) Lake Prespa. From 6 ka, �18Ocalcite suggest progressive aridification, in agreement with many other records in the Mediterranean, although the uppermost sediments in one core records low �18Ocalcite which we interpret as a result of human activity. Overall, the isotope data present here confirm that these two big lakes have captured the large scale, low frequency palaeoclimate variation that is seen in Mediterranean lakes, although in detail there is much palaeoclimate information that could be gained, especially small scale, high frequency differences between this region and the Mediterranean.

Predictive Margins and Marginal Effects in Stata

Tables of estimated regression coefficients, usually accompanied by additional information such as standard errors, t-statistics, p-values, confidence intervals or significance stars, have long been the preferred way of communicating results from statistical models. In recent years, however, the limits of this form of exposition have been increasingly recognized. For example, interpretation of regression tables can be very challenging in the presence of complications such as interaction effects, categorical variables, or nonlinear functional forms. Furthermore, while these issues might still be manageable in the case of linear regression, interpretational difficulties can be overwhelming in nonlinear models such as logistic regression. To facilitate sensible interpretation of such models it is often necessary to compute additional results such as marginal effects, predictive margins, or contrasts. Moreover, smart graphical displays of results can be very valuable in making complex relations accessible. A number of helpful commands geared at supporting these tasks have been recently introduced in Stata, making elaborate interpretation and communication of regression results possible without much extra effort. Examples of such commands are -margins-, -contrasts-, and -marginsplot-. In my talk, I will discuss the capabilities of these commands and present a range of examples illustrating their use.

Survey estimation and marginal effects in Stata

Stata is a general purpose software package that has become popular among various disciplines such as epidemiology, economics, or social sciences. Users like Stata for its scientific approach, its robustness and reliability, and the ease with which its functionality can be extended by user written programs. In this talk I will first give a brief overview of the functionality of Stata and then discuss two specific features: survey estimation and predictive margins/marginal effects. Most surveys are based on complex samples that contain multiple sampling stages, are stratified or clustered, and feature unequal selection probabilities. Standard estimators can produce misleading results in such samples unless the peculiarities of the sampling plan are taken into account. Stata offers survey statistics for complex samples for a wide variety of estimators and supports several variance estimation procedures such as linearization, jackknife, and balanced repeated replication (see Kreuter and Valliant, 2007, Stata Journal 7: 1-21). In the talk I will illustrate these features using applied examples and I will also show how user written commands can be adapted to support complex samples. Complex can also be the models we fit to our data, making it difficult to interpret them, especially in case of nonlinear or non-additive models (Mood, 2010, European Sociological Review 26: 67-82). Stata provides a number of highly useful commands to make results of such models accessible by computing and displaying predictive margins and marginal effects. In my talk I will discuss these commands provide various examples demonstrating their use.

Challenging conservation agriculture on marginal slopes in Sehoul, Morocco

In Sehoul, Morocco, the use of marginal land for agriculture became a necessity for the local population due to increased poverty and the occupation of the best land by new owners. Desertification poses an additional threat to agricultural production on marginal slopes, which are often stony and degraded. In a participatory process embedded in the EU DESIRE research project, potential sustainable land management measures were selected to address land degradation and desertification. Promising experiences with no-tillage practices elsewhere in Morocco had motivated the Moroccan government to promote conservation agriculture throughout the country. This combination of crop rotation, minimal soil disturbance and soil cover maintenance, however, had not yet been tested on sloping degraded land. Field trials of grazing enclosure combined with no or minimum tillage were conducted on the plots of two farmers, and trial results were analyzed based on stakeholders’ criteria. Results suggest that increased soil cover with barley residues improved rainwater use efficiency and yields only slightly, although soil water was generally enhanced. Soil moisture measurements revealed that no-tillage was favorable mainly at soil depths of 5 cm and in connection with low-rainfall events (<20 mm); under these circumstances, moisture content was generally higher under no-tillage than under conventional tillage. Moreover, stakeholder discussion confirmed that farmers in Sehoul remain primarily interested in animal husbandry and are reluctant to change the current grazing system. Implementation of conservation agriculture is thus challenged both by the degraded, sloping and stony nature of the land, and by the socio-economic circumstances in Sehoul.

The origin of oriented lakes: evidence from the Bolivian Amazon

The presence of hundreds of rectangular and oriented lakes is one of the most striking characteristics of the Llanos de Moxos (LM) landscape in the Bolivian Amazon. Oriented lakes also occur in the Arctic coastal plains of Russia, Alaska and Canada and along the Atlantic Coastal Plain from northeast Florida to southeast New Jersey and along the coast of northeast Brazil. Many different mechanisms have been proposed for their formation. In the LM, Plafker's (1964) tectonic model, in which subsidence results from the propagation of bedrock faults through the foreland sediments, is the most accepted. However, this model has not been verified. Here, we present new results from stratigraphic transects across the borders of three rectangular and oriented lakes in the LM. A paleosol buried under mid-Holocene sediments is used as a stratigraphic marker to assess the vertical displacement of sediments on both sides of the alleged faults. Our results show that there is no vertical displacement and, therefore, that Plafker's model can be ruled out. We suggest that, among all the proposed mechanisms behind lake formation, the combined action of wind and waves is the most likely. The evidence from the LM provides new hints for the formation of oriented lakes worldwide.

Increasing Silicon Concentrations in Bohemian Forest Lakes

Long-term trends of dissolved silicon (Si) concentrations in five glacial lakes in the Bohemian Forest, Czech Republic, recovering from acidification show higher mobility of Si from the soil to surface waters despite lower atmospheric deposition of acids. Si increased by 0.95 to 1.95 mu mol yr(-1) (36 to 51%) from 1986-2004 and with increasing pH. A change in soil solution conditions because of a sharp decrease in acidic deposition has led to marked decline in Al mobility and to considerable decreases in dissolved Al, especially Al3+. The increase in Si may be related to: (1) unblocking of the inhibitory effect of dissolved Al on weathering of alummosilicates, (2) biogenic opal (phytoliths) dissolving faster, and/or (3) lower Si precipitation as secondary alummosilicates in soil. The change in Al speciation on the dissolution rate of biogenic silica is critical. A lack of change in Si at sites outside central Europe may be explained by small or no decline in mobility of dissolved Al. The effect of a long-term increase in temperature was probably minor.

Understanding the performance of the FLake model over two African Great Lakes

The ability of the one-dimensional lake model FLake to represent the mixolimnion temperatures for tropical conditions was tested for three locations in East Africa: Lake Kivu and Lake Tanganyika's northern and southern basins. Meteorological observations from surrounding automatic weather stations were corrected and used to drive FLake, whereas a comprehensive set of water temperature profiles served to evaluate the model at each site. Careful forcing data correction and model configuration made it possible to reproduce the observed mixed layer seasonality at Lake Kivu and Lake Tanganyika (northern and southern basins), with correct representation of both the mixed layer depth and water temperatures. At Lake Kivu, mixolimnion temperatures predicted by FLake were found to be sensitive both to minimal variations in the external parameters and to small changes in the meteorological driving data, in particular wind velocity. In each case, small modifications may lead to a regime switch, from the correctly represented seasonal mixed layer deepening to either completely mixed or permanently stratified conditions from similar to 10 m downwards. In contrast, model temperatures were found to be robust close to the surface, with acceptable predictions of near-surface water temperatures even when the seasonal mixing regime is not reproduced. FLake can thus be a suitable tool to parameterise tropical lake water surface temperatures within atmospheric prediction models. Finally, FLake was used to attribute the seasonal mixing cycle at Lake Kivu to variations in the near-surface meteorological conditions. It was found that the annual mixing down to 60m during the main dry season is primarily due to enhanced lake evaporation and secondarily to the decreased incoming long wave radiation, both causing a significant heat loss from the lake surface and associated mixolimnion cooling.

Resin composites: Modulus of elasticity and marginal quality

OBJECTIVE To investigate how the modulus of elasticity of resin composites influences marginal quality in restorations submitted to thermocyclic and mechanical loading. METHODS Charisma, Filtek Supreme XTE and Grandio were selected as they were found to possess different moduli of elasticity but quite similar polymerization contraction. MOD cavities (n=30) were prepared in extracted premolars, restored and then subjected to thermocyclic and mechanical loading. Marginal quality of the restorations before and after loading was analyzed on epoxy replicas under a scanning electron microscope. The percentage of gap-free margins and occurrence of paramarginal fractures were registered. Modulus of elasticity and polymerization contraction were analyzed with parametric and margins with nonparametric ANOVA and post hoc Tukey HSD or Wilcoxon rank-sum tests, respectively. The number of paramarginal fractures was analyzed with exact Fisher tests (α=0.05). RESULTS Grandio demonstrated significantly more gap-free enamel margins than Charisma and Filtek Supreme XTE, before and after loading (p<0.01), whereas there was no difference between Charisma and Filtek Supreme XTE (p>0.05). No significant effect of resin composite (p=0.81) on the quality of dentine margins was observed, before or after loading. Deterioration of all margins was evident after loading (p<0.0001). More paramarginal enamel fractures were observed after loading in teeth restored with Grandio when compared to Charisma (p=0.008). CONCLUSIONS The resin composite with the highest modulus of elasticity resulted in the highest number of gap-free enamel margins but with an increased incidence of paramarginal enamel fractures. CLINICAL SIGNIFICANCE The results from this study suggest that the marginal quality of restorations can be improved by the selection of a resin composite with modulus of elasticity close to that of dentine, although an increase in paramarginal enamel fractures can result as a consequence.

Quantitative morphometric analysis of lakes using GIS: rectangularity R , ellipticity E , orientation O , and the rectangularity vs. ellipticity index, REi

Quantitative measures of polygon shapes and orientation are important elements of geospatial analysis. These kinds of measures are particularly valuable in the case of lakes, where shape and orientation patterns can help identifying the geomorphological agents behind lake formation and evolution. However, the lack of built-in tools in commercial geographic information system (GIS) software packages designed for this kind of analysis has meant that many researchers often must rely on tools and workarounds that are not always accurate. Here, an easy-to-use method to measure rectangularity R, ellipticity E, and orientation O is developed. In addition, a new rectangularity vs. ellipticity index, REi, is defined. Following a step-by-step process, it is shown how these measures and index can be easily calculated using a combination of GIS built-in functions. The identification of shapes and estimation of orientations performed by this method is applied to the case study of the geometric and oriented lakes of the Llanos de Moxos, in the Bolivian Amazon, where shape and orientation have been the two most important elements studied to infer possible formation mechanisms. It is shown that, thanks to these new indexes, shape and orientation patterns are unveiled, which would have been hard to identify otherwise.

Lake surface water temperatures of European Alpine lakes (1989–2013) based on the Advanced Very High Resolution Radiometer (AVHRR) 1 km data set

Lake water temperature (LWT) is an important driver of lake ecosystems and it has been identified as an indicator of climate change. Consequently, the Global Climate Observing System (GCOS) lists LWT as an essential climate variable. Although for some European lakes long in situ time series of LWT do exist, many lakes are not observed or only on a non-regular basis making these observations insufficient for climate monitoring. Satellite data can provide the information needed. However, only few satellite sensors offer the possibility to analyse time series which cover 25 years or more. The Advanced Very High Resolution Radiometer (AVHRR) is among these and has been flown as a heritage instrument for almost 35 years. It will be carried on for at least ten more years, offering a unique opportunity for satellite-based climate studies. Herein we present a satellite-based lake surface water temperature (LSWT) data set for European water bodies in or near the Alps based on the extensive AVHRR 1 km data record (1989–2013) of the Remote Sensing Research Group at the University of Bern. It has been compiled out of AVHRR/2 (NOAA-07, -09, -11, -14) and AVHRR/3 (NOAA-16, -17, -18, -19 and MetOp-A) data. The high accuracy needed for climate related studies requires careful pre-processing and consideration of the atmospheric state. The LSWT retrieval is based on a simulation-based scheme making use of the Radiative Transfer for TOVS (RTTOV) Version 10 together with ERA-interim reanalysis data from the European Centre for Medium-range Weather Forecasts. The resulting LSWTs were extensively compared with in situ measurements from lakes with various sizes between 14 and 580 km2 and the resulting biases and RMSEs were found to be within the range of −0.5 to 0.6 K and 1.0 to 1.6 K, respectively. The upper limits of the reported errors could be rather attributed to uncertainties in the data comparison between in situ and satellite observations than inaccuracies of the satellite retrieval. An inter-comparison with the standard Moderate-resolution Imaging Spectroradiometer (MODIS) Land Surface Temperature product exhibits RMSEs and biases in the range of 0.6 to 0.9 and −0.5 to 0.2 K, respectively. The cross-platform consistency of the retrieval was found to be within ~ 0.3 K. For one lake, the satellite-derived trend was compared with the trend of in situ measurements and both were found to be similar. Thus, orbital drift is not causing artificial temperature trends in the data set. A comparison with LSWT derived through global sea surface temperature (SST) algorithms shows lower RMSEs and biases for the simulation-based approach. A running project will apply the developed method to retrieve LSWT for all of Europe to derive the climate signal of the last 30 years. The data are available at doi:10.1594/PANGAEA.831007.

The stable carbon isotopic composition of Daphnia ephippia in small, temperate lakes reflects in-lake methane availability

Daphnia can ingest methane-oxidizing bacteria and incorporate methanogenic carbon into their biomass, leading to low stable carbon isotope ratios (expressed as δ13C values) of their tissue. Therefore, δ13C analysis of Daphnia resting eggs (ephippia) in lake sediment records can potentially be used to reconstruct past in-lake availability of methane (CH4). However, detailed multilake studies demonstrating that δ13C values of recently deposited Daphnia ephippia (δ13Cephippia) are systematically related to in-lake CH4 concentrations (CH4aq) are still missing. We measured δ13Cephippia from surface sediments of 15 small lakes in Europe, and compared these values with late-summer CH4aq. δ13Cephippia ranged from −51.6‰ to −25.9‰, and was strongly correlated with CH4aq in the surface water and above the sediment (r −0.73 and −0.77, respectively), whereas a negative rather than the expected positive correlation was found with δ13C values of carbon dioxide (CO2) (r −0.54), and no correlation was observed with CO2aq. At eight sites, offsets between δ13 CCO2 and δ13Cephippia exceeded offsets between δ13 CCO2 and δ13Calgae reported in literature. δ13Cephippia was positively correlated with δ13C values of sedimentary organic matter (r 0.54), but up to 20.7‰ lower in all except one of the lakes (average −6.1‰). We conclude that incorporation of methanogenic carbon prior to ephippia formation must have been widespread by Daphnia in our study lakes, especially those with high CH4aq. Our results suggest a systematic relationship between δ13Cephippia values and CH4aq in small temperate lakes, and that δ13Cephippia analysis on sediment records may provide insights into past changes in in-lake CH4aq.