Establishing a chronology for lacustrine sediments using a multiple dating approach—A case study from the Frickenhauser See, central Germany

In order to establish a reliable chronology for lacustrine sediments of the Frickenhauser See (central Germany) different dating methods have been applied. A total of 17 AMS ¹⁴C dates, all within the last 2000 years, were supplemented with ¹³⁷Cs/²¹⁰Pb dating and varve counting of the uppermost sediments (131 years). The age–depth model for the Frickenhauser See has to cope with highly variable sedimentation rates and overlapping probability distributions of calibrated ¹⁴C dates. The uncertainty of calibrated ¹⁴C dates could be considerably reduced by including the stratigraphic relationship of the dated samples, the age information derived from short-lived isotopes and varve counting as well as an upper and lower limit of realistic sedimentation rates as ‘a priori’ information in the calibration procedure. Sets of possible age combinations obtained by repeated sampling from the modified probability distributions were used to calculate continuous age–depth relationships based on monotonic smoothing splines. The obtained age–depth model for the sediment record of the Frickenhauser See represents the average of over 16,000 such model runs and suggests a drastic increase in sedimentation rates from around 1–2 mm a⁻¹ (200–1000 AD) to over 25 mm a⁻¹ for the period between 1100 and 1300 AD. From then on, sedimentation rates exhibit relatively stable values around 3–9 mm a⁻¹. ‘Conventional’ age–depth models such as general polynomial regression or cubic splines either do not include the obtained age-information in a satisfying manner (the model being too “stiff”) or exhibit “swings” causing age-reversals in the model. Although the age–depth relationships obtained for monotonic smoothing splines and mixed-effect regression are generally very similar, they differ in their respective sedimentation rates as well as in their uncertainties. Mixed-effect regression resulted in much higher sedimentation rates of more than 37 mm a⁻¹. These results suggest that monotonic smoothing splines give better control of the age–depth model characteristics and are well suited in situations, where the integrity of ¹⁴C dates is high, i.e. the dated material represents the age of the respective layer.     

A continuum representation of fracture networks. Part II: application to the Äspö Hard Rock laboratory

A method to represent fracture networks in a continuum model is evaluated by an application to the Äspö Hard Rock Laboratory (HRL). The basic approach of the method is to first generate a fracture network, based on data from field measurements and then represent the network as grid cell conductivities, using a staggered grid arrangement. It is argued that the resulting conductivity fields capture the correlation and anisotropy structures, as well as the mean conductivity, of the rock in a realistic way. The main objective is to evaluate if the suggested method can be used in an application to a real site. The objective is formulated in three questions: ‘Can the method utilize available field data?’, ‘Can the generated fields be calibrated and conditioned?’ and ‘Are realistic distributions of pressure, flow and density obtained?’. From the application to the Äspö HRL, it is concluded that the hydraulic conductivity fields, generated by the suggested method, can utilize field data and be calibrated. It is also demonstrated that realistic distributions of pressure, flow and density can be simulated.     

Pelagic tar, dissolved/dispersed petroleum hydrocarbons and plastic distribution in the Cretan Sea, Greece

During the first cruise of R/V ‘Philia’ in July 1997 within the framework of the TALOS programme supported by the Greek Ministry of Physical Planning and Public Works, the sampling of floating tar, litter and sea water for dissolved/dispersed petroleum hydrocarbons (DDPH) was carried out in the Cretan Sea. Analysis of these data has allowed a first assessment of the status of floating marine pollution in the region. DDPH measurements showed a mean concentration of 0.145 μg/l of chrysene equivalents (n = 24). Tar and plastics concentrations were in the range of 1–4280 and 0–1160 μg/m², respectively. Mean pelagic tar concentration was 318 μg/m², more than two times higher than what was reported for the area in previous studies. Based on in situ hydrographic observations there is strong evidence that most of the floating tar enters the Cretan Sea through the Ionian Sea.     

The Socio-Economic Costs and Benefits of Coastal Habitat Rehabilitation and Creation

This paper provides a comprehensive overview of the merits and limitations of using an economics based approach to assess and implement initiatives for coastal habitat rehabilitation and creation. A review of the literature indicates that habitat rehabilitation/creation costs vary widely between and within ecosystems. For coral reefs, costs range from US$ 10,000 to 6.5 million/hectare (ha);¹ for mangroves US$ 3000–510,000/ha; for seagrasses US$ 9000–680,000/ha and for saltmarshes US$ 2000–160,000/ha. A review of the economic benefits derived from various coastal habitats based on a ‘Total Economic Value’ approach (i.e. accounting for direct and indirect uses, and ‘non-uses’) reveals that many thousands of US$ per hectare could ultimately accrue from their rehabilitation/creation. The paper concludes that despite its limitations, the ‘benefit-cost analysis’ framework can play an important role both in assessing the justification of coastal habitat rehabilitation/creation initiatives, and by helping to improve the overall effectiveness of such initiatives.     

An acidification model (MAGIC) with organic acids evaluated using whole-catchment manipulations in Norway

The results from two whole-catchment manipulation experiments in Norway are used to evaluate MAGIC (Model of Acidification of Groundwater In Catchments), a model of ecosystem biogeochemical response to acid deposition. MAGIC is an aggregated catchment-scale model of acidification that has been widely used in assessment activities in Europe and North America. The experiments involved artificial decrease (Reversing Acidification In Norway, RAIN) and increases (Humid Lake Acidification Experiment, HUMEX) in acidic deposition. Runoff from both sites is influenced by moderate levels of organic acids. At each site the model was calibrated to the control catchment and then applied to the manipulated catchments with only minor adjustments. The major responses in runoff chemistry to the manipulations were closely simulated by the model. Differences between simulated and observed volume-weighted annual average concentrations of all major ions were less than ±6 μequiv 1⁻¹ for the entire 4–8 year period of prediction. Trends in response to the manipulations were correctly simulated. Most of the residual error resulted from an inability of the model to reproduce the year-to-year variability (noise) around the trends. Although such model evaluations cannot ‘prove’ the correctness of the model structure, good fits to experimental data increase confidence in model applications for assessment and management purposes. Evaluations of this sort can also identify aspects of the model that need further development. For MAGIC, these are primarily a need for improvement in the calibration of aluminium solubility and a better process basis for nitrogen dynamics.     

Efficient methods for solving water flow in variably saturated soils under prescribed flux infiltration

The non-linear solvers in numerical solutions of water flow in variably saturated soils are prone to convergence difficulties. Many aspects can give rise to such difficulties and in this paper we address the gravity term and the prescribed-flux boundary in the Picard iteration. The problem of the gravity term in the Picard iteration is iteration-to-iteration oscillation as the gravity term is treated, by analogy with the time-step advance technique, ‘explicitly’ in the iteration. The proposed method for the gravity term is an improvement of the ‘implicit’ approach of Zhang and Ewen [Water Resour. Res. 36 (2000) 2777] by extending it to heterogeneous soil and approximating the inter-nodal hydraulic conductivity in the diffusive term and the gravity term with the same scheme. The prescribed-flux boundary in traditional methods also gives rise to iteration-to-iteration oscillation because there is no feedback to the flux in the solution at the new iteration. To reduce such oscillation, a new method is proposed to provide such a feedback to the flux. Comparison with traditional Picard and Newton iteration methods for a wide range of problems show that a combination of these two proposed methods greatly improves the stability and consequently the computational efficiency, making the use of small time step and/or under-relaxation solely for convergence unnecessary.     

New estimates of future changes in extreme rainfall across the UK using regional climate model integrations. 2. Future estimates and use in impact studies

Under enhanced greenhouse conditions, climate models suggest an increase in rainfall intensities in the northern Hemisphere. Major flood events in the UK during autumn 2000 and central Europe in August 2002, have focussed attention on the dramatic impacts these changes may have on many sectors of society. In the companion paper [Fowler et al., J. Hydrol. (2004) this issue], we suggested that the HadRM3H model may be used with some confidence to estimate extreme rainfall distributions, showing good predictive skill in estimating statistical properties of extreme rainfall during the baseline period, 1961–1990. In this study, we use results from the future integration of HadRM3H (following the IPCC SRES scenario A2 for 2070–2100) to assess possible changes in extreme rainfall across the UK using two methods: regional frequency analysis and individual grid box analysis. Results indicate that for short duration events (1–2 days), event magnitude at a given return period will increase by 10% across the UK. For longer duration events (5–10 days), event magnitudes at given return periods show large increases in Scotland (up to +30%), with greater relative change at higher return periods (25–50 years). In the rest of the UK, there are small increases in the magnitude of more frequent events (up to +10%) but reductions at higher return periods (up to −20%). These results provide information to alter design storm depths to examine climate change impacts on various structures. The uncertainty bounds of the estimated changes and a ‘scaling’ methodology are additionally detailed. This allows the estimation of changes for the 2020s, 2050s and 2080s, and gives some confidence in the use of these estimates in impact studies.     

The errors in surface runoff prediction by neglecting the relationship between infiltration rate and overland flow depth

Many simplifications are used in modeling surface runoff over a uniform slope. A very common simplification is to determine the infiltration rate independent of the overland flow depth and to combine it afterward with the kinematic-wave equation to determine the overland flow depth. Another simplication is to replace the spatially variable infiltration rates along the slope i(x, t) due to the water depth variations h(x,t) with an infiltration rate that is determined at a certain location along the slope. The aim of this study is to evaluate the errors induced by these simplications on predicted infiltration rates, overland flow depths, and total runoff volume. The error analysis is accomplished by comparing a simplified model with a model where the interaction between the overland flow depth and infiltration rate is counted. In this model, the infiltration rate is assumed to vary along the slope with the overland flow depth, even for homogeneous soil profiles. The kinematic-wave equation with interactive infiltration rate, calculated along the slopy by Richard's equation, are then solved by a finite difference scheme for a 100-m-long uniform slope. In the first error analysis, we study the effect of combining an ‘exact’ and ‘approximate’ one-dimensional infiltration rate with the kinematic-wave equation for three different soil surface roughness coefficients. The terms ‘exact’ and ‘approximate’ stand for the solution of Richard's equation with and without using the overland flow depth in the boundary condition, respectively. The simulations showed that higher infiltration rates and lower overland flow depths are obtained during the rising stage of the hydrograph when overland flow depth is used in the upper boundary condition of the one-dimensional Richard's equation. During the recession period, the simplified model predicts lower infiltration rates and higher overland flow depths. The absolute relative errors between the ‘exact’ and ‘approximate’ solutions are positively correlated to the overland flow depths which increase with the soil surface roughness coefficient. For this error analysis, the relative errors in surface runoff volume per unit slope width throughout the storm are much smaller than the relative errors in momentary overland flow depths and discharges due to the alternate signs of the deviations along the rising and falling stages. In the second error analysis, when the spatially variable infiltration rate along the slope i(x, t) is replaced in the kinematic-wave equation by i(t), calculated at the slope outlet, the overland flow depth is underestimated during the rising stage of the hydrograph and overestimated during the falling stage. The deviations during the rising stage are much smaller than the deviations during the falling stage, but they are of a longer duration. This occurs because the solution with i(x, t) recognizes that part of the slope becomes dry after rainfall stops, while overland flow still exists with i(t) determined at the slope outlet. As obtained for the first error analysis, the relative errors in surface runoff volume per unit slope width are also much smaller than the relative errors in momentary overland flow depths and discharges. The relation between the errors in overland flow depth and discharge to different mathematical simplifications enables to evaluate whether certain simplifications are justified or more computational efforts should be used.     

Chemical change in Antarctica—Significance? A perspective

The bitter confrontations of 20–25 years ago between a small environmental science community and several industries manufacturing chlorinated hydrocarbons are now only an historical footnote in the rapid evolution of an environmental consciousness. Regulations developed in industrial countries are being applied to protect the most pristine of the world's environments in Antarctica. Many of the proposed regulations, however, would serve no useful purpose; some of the reports of environmental damage are grossly exaggerated. A certain skepticism might be appropriate. Nevertheless, although monitoring is not a ‘scientific’ activity, the case is made for a monitoring programme in Antarctica that would measure temporal changes in the environmental levels of contaminants from both local and distant sources.     

Calibration and uncertainty analysis of the SWAT model using Genetic Algorithms and Bayesian Model Averaging

In this paper, the Genetic Algorithms (GA) and Bayesian Model Averaging (BMA) were used to simultaneously conduct calibration and uncertainty analysis for the Soil and Water Assessment Tool (SWAT). In this combined method, several SWAT models with different structures are first selected; next GA is used to calibrate each model using observed streamflow data; finally, BMA is applied to combine the ensemble predictions and provide uncertainty interval estimation. This method was tested in two contrasting basins, the Little River Experimental Basin in Georgia, USA, and the Yellow River Headwater Basin in China. The results obtained in the two case studies show that this combined method can provide deterministic predictions better than or comparable to the best calibrated model using GA. The 66.7% and 90% uncertainty intervals estimated by this method were analyzed. The differences between the percentage of coverage of observations and the corresponding expected coverage percentage are within 10% for both calibration and validation periods in these two test basins. This combined methodology provides a practical and flexible tool to attain reliable deterministic simulation and uncertainty analysis of SWAT.     

Cadmium and manganese distributions in the Hudson River estuary: interannual and seasonal variability

Surface waters collected along the salinity gradient of the Hudson River estuary in four cruises between 1995 and 1997 were size-fractionated into particulate (>0.45 μm), ‘dissolved' (<0.45 μm), colloidal (10 kDa, 0.45 μm) and low molecular weight (<10 kDa) phases. Dissolved Cd concentrations (range: 0.11–1.19 nM) in surface waters of the estuary appear to have decreased fourfold (from an average of 2.36 to 0.61 nM) over a 23-year period, since the initial analysis of samples collected in the 1970s by Klinkhammer and Bender [Estuar. Coastal Shelf Sci. 12 (1981) 629–643]. This interannual decline reflects improvement in sewage treatment and the elimination of industrial Cd sources to the Hudson River estuary. In contrast, dissolved Mn levels (range: 0.033–1.46 μM) have remained relatively constant over the same period of time, suggesting that anthropogenic sources have very limited impact on Mn concentrations in the estuary. The concentrations of both Cd and Mn appeared to strongly depend on the season and/or river discharge. The highest concentrations were detected under low freshwater discharge, implying that limited hydraulic flushing allows a build-up of metals in the water column. Although the decline in Cd levels within the estuary reflects a reduction in the magnitude of anthropogenic inputs, mass balance estimates indicated that current sources of Cd to the estuary include sewage discharges (in the lower estuary around Manhattan) and diagenetic remobilization from industrial Cd deposited in sediments nearly 2 decades ago (in the upper estuary near Foundry Cove). Moreover, under low river discharge, the sources considered in our model (sewage, riverine input, atmospheric deposition, and benthic fluxes) could account for no more than 60% of the Cd exported from the lower estuary to the ocean. This suggests that undefined sources such as ground water and inputs from other watersheds (e.g., Long Island Sound and Newark Bay) may potentially influence the water quality of the New York Harbor. The size-fractionated metal concentrations indicated that most of the traditionally defined ‘dissolved' Cd and Mn consisted of <10 kDa molecular weight species. High molecular weight colloidal species of Mn accounted for about 50% of the dissolved fraction at the riverine end-member and <5% at intermediate and high salinities. Colloidal Cd accounted for <6% of the dissolved phase throughout the estuary. Unlike the non-conservative excess (relative to ideal dilution of river water and seawater) of dissolved Mn observed along the estuary, high molecular weight colloidal Mn appeared to be removed at the head of the estuary. The small contribution of colloidal Cd and Mn to the ‘dissolved' phase suggests that remobilization from suspended particulate phases and/or from sediments occurs through the formation of small molecular weight species.     

How do we know when we have done enough to protect the environment?

The changing concepts and challenges that the process industry has to face in relation to environmental protection are discussed. The original requirements on industry to reduce emissions of contaminants through the installation of filters and waste water treatment facilities or through process modification and introduction of ‘clean production’ strategies were mainly based on economic and technical considerations rather than on an effort to avoid environmental impacts. The introduction of concepts such as ‘assimilative capacity of the environment’ and ‘critical load’ provided useful instruments for setting effect-related emission standards, resulting in acceptable environmental protection. However, the actual requirement for ‘sustainability’ has forced industry to focus not only on controlling emissions of contaminants and pollutants, but to take a global environmental approach, including the choice of raw materials and energy sources, recycling and re-utilization of wastes and to take responsibility for the fate of their products, during their whole lifecycle.     

Temporary and definitive fixation of atmospheric lead in deep-sea sediments of the Western Mediterranean Sea

Most lead brought to the Mediterranean Sea has an anthropogenic origin and is mainly transported through the atmosphere. Atmospheric Pb was continuously collected at Cap Ferrat in 1986 and 1987. From this study, the estimation of the anthropogenic Pb flux on the whole Western Mediterranean was, averaged on 1986 and 1987 data, 4080 t. Assuming that the atmospheric anthropogenic Pb input varied in this course of time similarly to the consumption of Pb added to gasolines in France, the mean annual flux could be calculated: 3.95 kg km⁻² yr⁻¹, that is an annual input of 3360 t yr⁻¹. Reaching the sea, this metal seems to become rapidly bound to phytoplankton. Grazing by zooplankton leads to the production of faecal pellets which frequently contain rather high metal concentrations. The sinking rate of pellets of various zooplankton species is high; within a few days pellets may reach deep-sea sediments. After deposition, Pb is released from this organic-rich material during early diagenesis. In most cases, it, therefore, returns to the overlaying water body by ascending diffusion. But, in a deep-sea area of approximately 80 000 km² where Mn oxide precipitation occurs in surficial sediments, Pb seems to remain stored by coprecipitation processes. By considering the lead stored in ‘excess’ in the surficial sediment of the deep-sea area, we estimate that a mean annual anthropogenic Pb amount ranging from 800 up to 1080 t was stored every year from 1950. On the same area, taking into account the Pb loss at the straits, the ‘direct’ atmospheric input to the sea bottom is, on average, 184 t yr⁻¹. The remaining part, that is (800–1080)−184=(616–896) t yr⁻¹, corresponds to an additional ‘indirect’ Pb flux in water due to Pb released from sediments of the surrounding areas where it does not remain stored.     

Isotope geochemistry of Pantelleria volcanic fluids, Sicily Channel rift: a mantle volatile end-member for volcanism in southern Europe

Chemical and isotopic ratio (He, C, H and O) analysis of hydrothermal manifestations on Pantelleria island, the southernmost active volcano in Italy, provides us with the first data upon mantle degassing through the Sicily Channel rift zone, south of the African–European collision plate boundary. We find that Pantelleria fluids contain a CO₂–He-rich gas component of mantle magmatic derivation which, at shallow depth, variably interacts with a main thermal (100°C) aquifer of mixed marine–meteoric water. The measured ³He/⁴He ratios and δ¹³C of both the free gases (4.5–7.3 R_a and −5.8 to −4.2‰, respectively) and dissolved helium and carbon in waters (1.0–6.3 R_a and −7.1 to −0.9‰), together with their covariation with the He/CO₂ ratio, constrain a ³He/⁴He ratio of 7.3±0.1 R_a and a δ¹³C of ca. −4‰ for the magmatic end-member. These latter are best preserved in fluids emanating inside the active caldera of Pantelleria, in agreement with a higher heat flow across this structure and other indications of an underlying crustal magma reservoir. Outside the caldera, the magmatic component is more affected by air dilution and, at a few sites, by mixing with either organic carbon and/or radiogenic ⁴He leached from the U–Th-rich trachytic host rocks of the aquifer. Pantelleria magmatic end-member is richer in ³He and has a lower (closer to MORB) δ¹³C than all fluids yet analyzed in volcanic regions of Italy and southern Europe, including Mt. Etna in Sicily (6.9±0.2 R_a, δ¹³C=−3±1‰). This observation is consistent with a south to north increasing imprint of subducted crustal material in the products of Italian volcanoes, whose He and C (but also O and Sr) isotopic ratios gradually evolve towards crustal values northward of the African–Eurasian plate collision boundary. Our results for Pantelleria extend this regional isotopic pattern further south and suggest the presence of a slightly most pristine or ‘less contaminated’, ³He-richer mantle source beneath the Sicily Channel rift zone. The lower than MORB ³He/⁴He ratio but higher than MORB CO₂/³He ratio of Pantelleria volatile end-member are compatible with petro-geochemical evidence that this mantle source includes an upwelling HIMU–EM1-type asthenospheric plume component whose origin, according to recent seismic data, may be in the lower mantle.     

Simulating sensitivity of soil moisture and evapotranspiration under heterogeneous soils and land uses

Soil moisture (SM) plays an important role in land surface and atmospheric interactions. It modifies energy balance at the surface and the rate of water cycling between the land and atmosphere. In this paper we provide a sensitivity assessment of SM and ET for heterogeneous soil physical properties and for three land uses including irrigated maize, rainfed maize, and grass at a climatological time-scale by using a water balance model. Not surprisingly, the study finds increased soil water content in the root zone throughout the year under irrigated farming. Soil water depletes to its lowest level under rainfed maize cultivation. We find a ‘land use’ effect as high as 36 percent of annual total evapotranspiration, under irrigated maize compared to rainfed maize and grass, respectively. Sensitivity analyses consisting of comparative simulations using the model show that soil characteristics, like water holding capacity, influence SM in the root zone and affect seasonal total ET estimates at the climatological time-scale. This ‘soils’ effect is smaller than the ‘land use’ effect associated with irrigation but, it is a source of consistent bias for both SM and ET estimates. The ‘climate’ effect basically masks the ‘soils’ effect under wet conditions. These results lead us to conclude that appropriate representation of land use, soils, and climate are necessary to accurately represent the water and energy balance in real landscapes.     

Formation of Paleoproterozoic eclogitic mantle, Slave Province (Canada): Insights from in-situ Hf and U–Pb isotopic analyses of mantle zircons

In-situ Hf isotope analyses and U–Pb dates were obtained by laser ablation-MC-ICP-MS for a zircon-bearing mantle eclogite xenolith from the diamondiferous Jericho kimberlite located within the Archean Slave Province (Nunavut), Canada. The U–Pb zircon results yield a wide range of ages (2.0 to 0.8 Ga) indicating a complex geological history. Of importance, one zircon yields a U–Pb upper intercept date of 1989 ± 67 Ma, providing a new minimum age constraint for zircon crystallization and eclogite formation. In contrast, Hf isotope systematics for the same zircons display an intriguing uniformity, and corresponding Hf depleted mantle model ages range between 2.1 ± 0.1 and 2.3 ± 0.1 Ga; the youngest Hf model age is within error to the oldest U–Pb date.

The Jericho eclogites have previously been interpreted as representing remnants of metamorphosed oceanic crust, and their formation related to Paleoproterozoic subduction regimes along the western margin of the Archean Slave craton during the Wopmay orogeny. Hf isotope compositions and U–Pb results for the Jericho zircons reported here are in good agreement with a Paleoproterozoic subduction model, suggesting that generation of oceanic crust and eclogite formation occurred between 2.0 and 2.1 Ga. The slightly older Hf depleted mantle model ages (2.1 to 2.3 Ga) may be reconciled with this model by invoking mixing between ‘crustal’-derived Hf from sediments and more radiogenic Hf associated with the oceanic crust during the 2 Ga subduction event. This results in intermediate Hf isotope compositions for the Jericho zircons that yield ‘fictitiously’ older Hf model ages.     

Evolution of Earth's redox state during upwelling of carbon-bearing mantle

The oxygen fugacity (f(O₂)) values recorded by diamondiferous peridotite and eclogite xenoliths from Siberia indicate that the redox state of the ancient lithosphere is heterogeneous on a scale of at least four log units, mainly in the range between the wüstite-magnetite (WM) and iron-wüstite (IW) oxygen buffers. Highly reduced peridotites can be interpreted as relict from earlier lower f(O₂). The f(O₂) values recorded by ‘fertile’ and less modified spinel peridotites from Mongolia, Baikal and Tien-Shan show that the redox state of the lithosphere beneath central Asia and Tien-Shan is heterogeneous on a scale of 2–3 log units, mainly in the range between the WM and IW + 1 oxygen buffers. These data provide evidence for the presence of a lower-f(O₂) regime of carbon-bearing mantle beneath the Baikal rift zone and Tien-Shan, and the oxidation of diapirs ascending from the asthenosphere. The ‘dry’ xenoliths from Mongolia primarily reflect closed system behavior in the upper mantle, the f(O₂) of which is buffered by ferric-ferrous redox equilibrium. The observed evolution of f(O₂) values is closely linked to the distribution of volatile species in the mantle. H₂O and CO₂ are the dominant volatiles for the more depleted and oxidized part of peridotites, and CH₄ for the more reduced and less modified part. It is proposed that the upper mantle was originally more reduced and has become progressively more oxidized, resulting perhaps largely from the preferential loss of hydrogen and carbon during melting. The oxygen budget of the upper mantle results from the opposing contributions of crustal recycling and transfer of carbon-bearing material from the deep mantle.     

Assessing morphological degradation in running waters using Blackfly communities (Diptera, Simuliidae): Can habitat quality be predicted from land use?

Based on benthic macroinvertebrate samples from 57 sites of streams in the European Central Highlands (Ecoregion 9; Illies (1978). Limnofauna Europaea, Stuttgart), the composition and the abundance of Simuliid species were analysed in relation to hydromorphological and land use parameters. Sampling sites were located at two stream types differing in catchment geology, stream morphology and channel width. Land use data were taken from the official German information system for cartography and topography (ATKIS) and the German River Habitat Survey. Land use categories suitable to describe the sampling sites were ‘agricultural land’, ‘urban’ and near-natural areas. Hydromorphological parameters of the sites were recorded using the AQEM site protocol (AQEM Consortium (2002). Manual for the application of the AQEM system. A comprehensive method to assess European streams using benthic macroinvertebrates, developed for the purpose of the Water Framework Directive. EVK1-CT1999-00027, Version 1.0. Available via the Internet from www.aqem.de). Relevant parameters controlling Simuliid distribution in streams were the mineral substrates and the biotic microhabitats. Blackfly larvae and pupae were sampled at each site for 15 min. Statistical analysis was performed by CANOCO^® (Ter Braak & Smilauer (1997). CANOCO Version 4.5. Biometrics Plant Research International, Wageningen, The Netherlands) using redundancy analysis (RDA).

Our results show a stream type-specific composition of the Blackfly fauna. The analyses reveal sensitivity of the Simuliid species to morphological degradation, which is indicated by the shift in the longitudinal zonation of the Simuliid communities. Especially, Prosimulium hirtipes (Fries, 1824) and Simulium argyreatum Meigen, 1838 are typical representatives of headwater streams. While they seem to indicate undisturbed conditions of this stream type and react sensitively to the degradation of stream habitats, Simulium ornatum Meigen, 1818 and Simulium equinum (Linnaeus, 1758) are more tolerant to stream degradation.

On catchment scale, ‘% natural forest’ and ‘% agricultural land use’ illustrate the degree of degradation of the two selected stream types. ‘Average stream depth’ and ‘relation riffles/pools’ account for hydromorphological degradation reflected by Simuliid species on the smaller reach scale. The analysis of habitat quality revealed that ‘% woody debris’ represents an important parameter of morphological degradation reflected by the Blackfly community.     

Seasonal variability of benthic indices: an approach to test the applicability of different indices for ecosystem quality assessment

In recent studies, several benthic biological indices were developed or used to assess the ecological quality status of marine environments. In the present study the seasonal variability of several univariate and multimetric indices was studied on a monthly scale (September 2000 until May 2002) in different areas of the North Sea such as the German Bight, the Oyster Ground and the Dogger Bank. The stations were chosen to reflect a gradient in the hydrographic regime, temperature and organic matter supply. The seasonal variability was highest for the univariate indices such as the Shannon–Wiener Index and the Hurlbert Index. Thus, due to sensitivity to recruitment the corresponding ecological status ranged from ‘good’ to ‘poor’ depending on the season. For the multimetric indices such as the AMBI or the BQI the seasonal variability and the corresponding ecological status were low. The results are discussed concerning possible consequences for ecological quality assessment especially related to the Water Framework Directive (WFD).     

Noble gas systematics of the Réunion mantle plume source and the origin of primordial noble gases in Earth’s mantle

New noble gas data of ultramafic xenoliths from Réunion Island, Indian Ocean, further constrain the characteristics of primordial and radiogenic noble gases in Earth’s mantle plume reservoirs. The mantle source excess of nucleogenic ²¹Ne is significantly higher than for the Hawaiian and Icelandic plume reservoirs, similar to excess of radiogenic ⁴He. ⁴⁰Ar/³⁶Ar of the Réunion mantle source can be constrained to range between 8000 and 12 000, significant ¹²⁹Xe and fission Xe excess are present. Regarding the relative contribution of primordial and radiogenic rare gas nuclides, the Réunion mantle source is intermediate between Loihi- and MORB-type reservoirs. This confirms the compositional diversity of plume sources recognized in other radioisotope systematics. Another major result of this study is the identification of the same basic primordial component previously found for the Hawaiian and Icelandic mantle plumes and the MORB reservoir. It is a hybrid of solar-type He and Ne, and ‘atmosphere-like’ or ‘planetary’ Ar, Kr, Xe (Science 288 (2000) 1036). ²⁰Ne/²²Ne ratios extend to maximum values close to 12.5 (Ne-B), which is the typical signature of solar neon implanted as solar corpuscular radiation. This suggests that Earth’s solar-type noble gas inventory was acquired by small (less than km-sized) precursor planetesimals that were irradiated by an active early sun in the accretion disk after nebular gas dissipation, or, alternatively, that planetesimals incorporated constituents irradiated in transparent regions of the solar nebula. Previously, such an early irradiation scenario was suggested for carbonaceous chondrites which follow common volatile depletion trends in the sequence CI–CM–CV–Earth. In turn, CV chondrites closely match Earth’s mantle composition in ²⁰Ne/²²Ne, ³⁶Ar/²²Ne and ³⁶Ar/³⁸Ar. This indicates that mantle Ar could well be a planetary component inherited from precursor planetesimals. However, a corresponding conclusion for mantle Kr and Xe is less convincing yet, but this may be just due to the lack of appropriate ‘meteoritic’ building blocks matching terrestrial composition. Alternatively, heavy noble gases in Earth’s mantle could be due to admixing of severely fractionated air, but this effect must have affected all mantle sources to a very similar extent, e.g. by global subduction before the last homogenization of the mantle reservoirs.