A palaeolimnological meta-database for assessing the ecological status of lakes

We describe a palaeolimnological meta-database (“LakeCores”) that includes information for lakes in Europe that have been cored for studies of recent environmental change. Only cores with intact mud-water interfaces and good chronologies for the last ~100–200 years are currently included. The information provided in the database includes the general characteristics of the lakes cored, the methods used for dating and analysis and the bibliographic source of the information selected. At present the database contains information on 1847 studies from 975 sites. Here we use the database to identify cores for which diatom-based pH and total phosphorus (TP) reconstructions have been derived and, from an inspection of the relevant papers, we compile summary data for the timing and extent of acidification and eutrophication of lakes across Europe. We show that there were few cases of acidification and eutrophication in Europe prior to about 1850, and therefore that conditions prior to 1850 can be used for reference with respect to current attempts to restore surface waters to “good ecological status”. We also show the extent to which pH and TP levels have changed across the continent both for all lakes, and, in the case of TP, for lakes of different types covering a range of alkalinities, altitudes, water depths and surface areas.     

Reconstructing 4000 years of mass movement and tsunami history in a deep peri‐Alpine lake (Lake Geneva,France‐Switzerland)

The study of mass movements in lake sediments provides insights into past natural hazards at historic and prehistoric timescales. Sediments from the deep basin of Lake Geneva reveal a succession of six large‐scale (volumes of 22 × 10⁶ to 250 × 10⁶ m³) mass‐transport deposits, associated with five mass‐movement events within 2600 years (4000 cal bp to 563 ad ). The mass‐transport deposits result from: (i) lateral slope failures (mass‐transport deposit B at 3895 ± 225 cal bp and mass‐transport deposits A and C at 3683 ± 128 cal bp ); and (ii) Rhône delta collapses (mass‐transport deposits D to G dated at 2650 ± 150 cal bp , 2185 ± 85 cal bp , 1920 ± 120 cal bp and 563 ad , respectively). Mass‐transport deposits A and C were most probably triggered by an earthquake, whereas the Rhône delta collapses were likely to be due to sediment overload with a rockfall as the external trigger (mass‐transport deposit G, the Tauredunum event in 563 ad known from historical records), an earthquake (mass‐transport deposit E) or unknown external triggers (mass‐transport deposits D and F). Independent of their origin and trigger mechanisms, numerical simulations show that all of these recorded mass‐transport deposits are large enough to have generated at least metre‐scale tsunamis during mass movement initiation. Since the Tauredunum event in 563 ad , two small‐scale (volumes of 1 to 2 × 10⁶ m³) mass‐transport deposits (H and I) are present in the seismic record, both of which are associated with small lateral slope failures. Mass‐transport deposits H and I might be related to earthquakes in Lausanne/Geneva (possibly) 1322 ad and Aigle 1584 ad , respectively. The sedimentary record of the deep basin of Lake Geneva, in combination with the historical record, show that during the past 3695 years, at least six tsunamis were generated by mass movements, indicating that the tsunami hazard in the Lake Geneva region should not be neglected, although such events are not frequent with a recurrence time of 0·0016 yr^?1.     

Assessing risk to groundwater quality using an integrated risk framework

Source-water protection strategies are ideally focused where the greatest amount of harm reduction can occur. This process of risk management requires an assessment of the spatial variability of risk to water. The assessment methodology presented herein combines aquifer susceptibility with a hazard threat inventory and an analysis of the consequence of contamination to assess the risk to water quality. Aquifer susceptibility combines the intrinsic susceptibility of the physical system with anthropogenic features that locally increase susceptibility. Hazard threats are assessed based on the properties of the chemicals (toxicity and environmental fate), the potential magnitude (extent and quantity of release) and the likelihood of release. The consequence is herein considered as the financial costs of the loss of the resource, including the replacement of a water source and the economic loss where water intensive businesses are lost. A second scenario is included that analyses health issues related to pathogen sources as well as the financial impact to the community where people fall ill and present a financial burden to the public health care system. The risk assessment methodology is applied to the Township of Langley, in southwestern British Columbia, Canada. The results outline the most vulnerable areas as those where susceptible aquifers coexist with potential chemical and biological threats. The risk is greatest where these vulnerable areas coincide with those with the greatest potential for financial loss: within the capture zones of major municipal production wells and where private wells serve agricultural operations with high annual farm sales.     

‘Diatoms and pH reconstruction’ (1990) revisited

The 167 sample lake-water pH-diatom calibration data-set created as part of the Palaeolimnology Programme within the Surface Water Acidification Project (SWAP) is re-analysed numerically using nine different numerical methods, six based on simple two-way weighted-averaging (WA), and the other three involving Gaussian logit regression (GLR) and maximum-likelihood (ML) calibration, the modern analogue technique, or weighted-averaging partial least-squares regression and calibration. Root mean squared error of prediction and maximum bias were estimated for all nine methods based on 10,000 internal and 10,000 external cross-validations involving a training-set, an optimisation-set, and a test-set. The results show that WA with a monotonic deshrinking spline equals or slightly outperforms WA with linear inverse deshrinking, especially in external cross-validation. Methods that employ tolerance downweighting generally have an inferior performance except when combined with monotonic deshrinking. It appears that simple two-way WA extensively used in SWAP cannot be significantly bettered. Thanks to increased computing power, better software, and more rigorous cross-validations, GLR shows good performance, especially in external cross-validation.     

Stratigraphy,distribution and geochemistry of widespread felsic volcanic units in the Mesoproterozoic Gawler Range Volcanics,South Australia

Three widespread felsic volcanic units, the Eucarro Rhyolite, Pondanna Dacite Member and Moonaree Dacite Member, have been distinguished in the Mesoproterozoic Gawler Range Volcanics. These three units are the largest in the Gawler Range Volcanics, each in excess of 500 km³. Each unit is ～300 m thick and includes a black, formerly glassy base, a granophyric columnar‐jointed interior, and an amygdaloidal outer margin. The units are very gently dipping and locally separated by thin (<20 m) lenses of either ignimbrite (Mt Double Ignimbrite), tuffaceous sandstone or faults. The youngest unit, the Moonaree Dacite Member, covers a central area with a diameter greater than 80 km. The southern two units have east‐west extents in the order of 180 km, but are much less extensive from south to north (5–60 km). All three units are dominated by euhedral phenocrysts and are relatively crystal rich. Both the Eucarro Rhyolite and Moonaree Dacite Member contain clasts of basement granitoid and other lithologies and are locally heterogeneous in texture and composition. Some granitoid clasts have disintegrated, liberating feldspar and quartz crystals into the surrounding host. These liberated crystals cause textural variations, but can be identified on the basis of shape (amoeboid or skeletal) and/or size (megacrysts). Textural and lithofacies characteristics are consistent with the interpretation that these units are lavas; the strongly elongate distribution and wide extent of the Eucarro Rhyolite and Pondanna Dacite Member could indicate that vents were aligned along an extensive east‐west‐trending fissure system. Stratigraphic nomenclature has been revised to better reflect the presence of the three emplacement units. The oldest unit, the Eucarro Rhyolite, is dominated by plagioclase‐phyric rhyolite that locally includes granitoid clasts and megacrysts. Along the northern margin, the rhyolite is amygdaloidal and has mingled with a quartz‐rich rhyolite (Paney Rhyolite Member). The Eucarro Rhyolite and Paney Rhyolite Member replace the formerly defined ‘Eucarro Dacite’, ‘Nonning Rhyodacite’, ‘Yannabie Rhyodacite’ and ‘Paney Rhyolite’. The two younger units, Pondanna Dacite Member and Moonaree Dacite Member, are compositionally and spatially distinct, newly defined members of the Yardea Dacite.     

Considerations for modeling burn probability across landscapes with steep environmental gradients: an example from the Columbia Mountains, Canada

Fire and land management in fire-prone areas can be greatly enhanced by estimating the likelihood of fire at every point on the landscape. In recent years, powerful fire simulation models, combined with an in-depth understanding of an area’s fire regime and fire environment, have allowed forest managers to estimate spatial burn probabilities. This study describes a methodology for selecting input data and model parameters when creating burn probability maps in difficult-to-model areas and reports the results of a case study for a large area of the Columbia Mountains, British Columbia, Canada. In addition to having particularly mountainous topography, the study area is covered by vegetation types that are poorly represented in fire behavior systems, even though these vegetation types have experienced considerable (if highly irregular) fire activity in premodern times (before 1920). Parameterization of the fire environment for simulation modeling was accomplished by combining various types of fire information (e.g., fire history studies, reconstructed fire climatologies), new technologies (high-resolution remotely sensed data, wind flow modeling), and—a must in data-limited areas—ample expert advice. In this study, we made extensive use of personal accounts from experienced fire behavior officers for the creation of model inputs. Despite difficulties in validating outputs of burn probability models, the multisource model-building approach described here provides a conservative, yet informative, means of estimating the likelihood of fire. Due to the data-intensive nature of the modeling and paucity of input data, an argument is made that modelers must focus on the inputs that are the most influential for their study area.     

A Bayesian palaeoenvironmental transfer function model for acidified lakes

A Bayesian approach to palaeoecological environmental reconstruction deriving from the unimodal responses generally exhibited by organisms to an environmental gradient is described. The approach uses Bayesian model selection to calculate a collection of probability-weighted, species-specific response curves (SRCs) for each taxon within a training set, with an explicit treatment for zero abundances. These SRCs are used to reconstruct the environmental variable from sub-fossilised assemblages. The approach enables a substantial increase in computational efficiency (several orders of magnitude) over existing Bayesian methodologies. The model is developed from the Surface Water Acidification Programme (SWAP) training set and is demonstrated to exhibit comparable predictive power to existing Weighted Averaging and Maximum Likelihood methodologies, though with improvements in bias; the additional explanatory power of the Bayesian approach lies in an explicit calculation of uncertainty for each individual reconstruction. The model is applied to reconstruct the Holocene acidification history of the Round Loch of Glenhead, including a reconstruction of recent recovery derived from sediment trap data. The Bayesian reconstructions display similar trends to conventional (Weighted Averaging Partial Least Squares) reconstructions but provide a better reconstruction of extreme pH and are more sensitive to small changes in diatom assemblages. The validity of the posteriors as an apparently meaningful representation of assemblage-specific uncertainty and the high computational efficiency of the approach open up the possibility of highly constrained multiproxy reconstructions.     

Quantification of braided river channel change using archival digital image analysis

Historical archives of grey‐scale river channel imagery are extensive. Here, we present and test a methodology to extract detailed quantitative topographic data from such imagery of sand‐bed rivers. Extracting elevation information from rivers is difficult as they are characterized by a low relative relief (<4 m); the area of interest may be spatially extensive (e.g. active channel widths >500 m in large braided rivers); the rate of change of surface elevation is generally low except in the vicinity of individual channel banks where the rate of change is very high; there is the complication that comes from inundation; and there may be an added complication caused by blockage of the field of view by vegetation. Here, we couple archival photogrammetric techniques with image processing methods and test these for quantification of sand‐bed braided river dynamics, illustrated for a 500 m wide, 3 km long reach of the South Saskatchewan River, Canada. Digital photogrammetry was used to quantify dry areas and water edge elevations. A methodology was then used to calibrate the spectral signature of inundated areas by combining established two media digital photogrammetric methods and image matching. This allowed determination of detailed depth maps for inundated areas and, when combined with dry area data, creation of complete digital elevation models. Error propagation methods were used to determine the erosion and deposition depths detectable from sequential digital elevation models. The result was a series of elevation models that demonstrate the potential for acquiring detailed and precise elevation data from any historical aerial imagery of rivers without needing associated calibration data, provided that imagery is of the necessary scale to capture the features of interest. We use these data to highlight several aspects of channel change on the South Saskatchewan River, including bar movement, bank erosion and channel infilling. Copyright © 2010 John Wiley & Sons, Ltd.