Dating of Late Pleistocene terrace deposits of the River Rhine using Uranium series and luminescence methods: Potential and limitations

Uranium-series dating of pedogenic carbonate crusts from fluvial gravels is tested using Optically Stimulated Luminescence (OSL) ages as references. OSL dating yielded ages of 30–15 ka and 13–11 ka, which correlate with the cold periods of the Last Glacial Maximum and the Younger Dryas. These ages are internally coherent and consistent with the geological background and are thus regarded as reliable. Most of the U/Th results scatter widely in the ²³⁰Th/²³²Th vs. ²³⁴U/²³²Th isochron diagram, making regression unrealistic. Semiquantitative age estimates from the data were found to be mostly older than the OSL ages and the geological context. It is suggested that a heterogeneous initial ²³⁰Th input, not relatable to a detrital component, is responsible for the observed discrepancies. This input may be due to bacterial activities and Th transport on organic colloids. It appears necessary to avoid samples where bacteria could have contributed to carbonate precipitation. Further, the relative importance of this problem decreases with sample age, so that U/Th dating of sinters is expected to be more reliable in the >100 ka age range.     

Fine particulate organic matter (FPOM) transport and processing in littoral interstices – use of fluorescent markers

A fluorescent labelling method is presented as a new tool for the investigation of organic particle transport and biogenic carbon cycling processes in sandy littoral interstices at Lake Tegel, Berlin, Germany. Passive particle transport through the pore system was studied by in situ exposition of 2.4 μm monodisperse polymeric resin microparticles stained with 7-amino-4-methylcoumarin (AMC). Uptake of fluorescein-5-isothiocyanate (FITC)-labelled Chlorella vulgaris and fine particulate organic matter (FPOM) by the interstitial fauna was investigated in laboratory and field experiments. The major portion (>85%) of the FITC-labelled particles added to sediment cores was recovered from the topmost centimetre of sediment during the study period of 14 days. Uptake of FITC-labelled FPOM was observed in several benthic groups, e.g. chironomids, microcrustaceans, oligochaetes and tardigrads, whereas C. vulgaris was ingested by oligochaetes only. There is evidence to suggest that both are suitable materials for investigating the fragmentation and ingestion of organic material by herbivorous and detritivorous fauna. Field experiments with inert microparticles and FITC-labelled FPOM (<1 mm) prepared from dried alder leaves were carried out in plexiglass tubes as in situ whole core technique. Within the investigation period of two weeks, the transport of FPOM was restricted to the topmost 2–3 cm of sediment in conjunction with a distinct fragmentation to finer size classes with respect to increasing sediment depth. Vertical FPOM transport was hindered by a high interstitial concentration of natural POM and an intensive settlement of the interstices by algae (mainly epispammic algae, 65–96% of algae cell number) and extra-cellular polymeric substances (EPS), which formed a dense three-dimensional structure.     

Meltwater discharge through the subglacial bed and its land‐forming consequences from numerical experiments in the Polish lowland during the last glaciation

Numerical experiments suggest that the last glaciation severely affected the upper lithosphere groundwater system in NW Poland: primarily its flow pattern, velocities and fluxes. We have simulated subglacial groundwater flow in two and three spatial dimensions using finite difference codes for steady‐state and transient conditions. The results show how profoundly the ice sheet modifies groundwater pressure heads beneath and some distance beyond the ice margin. All model runs show water discharge at the ice forefield driven by ice‐sheet‐thickness‐modulated, down‐ice‐decreasing hydraulic heads. In relation to non‐glacial times, the transient 3D model shows significant changes in the groundwater flow directions in a regionally extensive aquifer ca. 90 m below the ice–bed interface and up to 40 km in front of the glacier. Comparison with empirical data suggests that, depending on the model run, only between 5 and 24% of the meltwater formed at the ice sole drained through the bed as groundwater. This is consistent with field observations documenting abundant occurrence of tunnel valleys, indicating that the remaining portion of basal meltwater was evacuated through a channelized subglacial drainage system. Groundwater flow simulation suggests that in areas of very low hydraulic conductivity and adverse subglacial slopes water ponding at the ice sole was likely. In these areas the relief shows distinct palaeo‐ice lobes, indicating fast ice flow, possibly triggered by the undrained water at the ice–bed interface. Owing to the abundance of low‐permeability strata in the bed, the simulated groundwater flow depth is less than ca. 200 m. Copyright © 2009 John Wiley & Sons, Ltd.     

Age and nature of eclogites in the Huwan shear zone,and the multi-stage evolution of the Qinling-Dabie-Sulu orogen,central China

In situ LA-ICPMS U-Pb, trace element, and Hf isotope data in zircon demonstrate a Carboniferous age for eclogite-facies metamorphism in Siluro-Devonian protoliths in the Huwan shear zone, Dabie Mountains, Central China. This age contrasts with the more prevailing Triassic age for high- to ultrahigh pressure (HP to UHP) metamorphism in the Qinling-Dabie-Sulu orogen. Metamorphic zircon in two eclogite samples from Sujiahe is characterized by low Th/U ratios, small negative Eu anomalies, flat HREE patterns, and low ¹⁷⁶Lu/¹⁷⁷Hf ratios. These geochemical signatures suggest that the zircon crystallized in the presence of garnet and in the absence of plagioclase feldspar. Furthermore, temperatures of ~ 655 and ~ 638 °C, calculated using the Ti content of zircon, are consistent with their formation during eclogite-facies metamorphism. The weighted mean ²⁰⁶Pb/²³⁸U age of 309 ± 4 Ma (2δ) for this zircon improves previous age estimates for eclogite-facies metamorphism in the Huwan shear zone, ranging from 420 to 220 Ma. Metamorphic zircon from one eclogite sample from Hujiawan, most likely formed during prograde metamorphism, yields an equivalent age estimate of 312 ± 11 Ma. Magmatic zircon cores in the three samples yield ages for the magmatic protoliths of the eclogites ranging from 420 ± 7 to 406 ± 5 Ma, and post-dating the middle Paleozoic collision of the North China and the Qinling terrain. The zircon crystals in the three eclogite samples display a large variation of ε_Hf (t) values of ? 4.9 to 21.3. The metamorphic zircon overgrowths show the same range of ε_Hf (t) values as those of the inherited magmatic crystal interiors. This suggests that the metamorphic zircon overgrowths may have formed by dissolution-reprecipitation of pre-existing magmatic zircon thereby preserving their original Hf isotopic composition. The high ε_Hf (t) values suggest that the protoliths were derived from depleted mantle sources, most likely Paleotethyan oceanic crust; while the low ε_Hf (t) values are attributed to crustal contamination. Some eclogites in the Huwan shear zone have a distinctive signature of continental crust most probably derived from the Yangtze Craton. The coexistence of Paleozoic oceanic crust and Neoproterozoic continental crust with similar metamorphic ages in the Huwan shear zone implies that Paleozoic Paleotethyan oceanic crust was produced within a marginal basin of the northern Yangtze Craton. The opening of the Paleo-Tethyan ocean was slightly younger than the collision of the North China Craton and the Qinling terrain during the Late Paleozoic in the Qinling-Dabie-Sulu orogen. Subduction of the Paleo-Tethyan oceanic crust and associated continental basement resulted in the 309 ± 2 Ma (2σ) eclogite-facies metamorphism in the Huwan shear zone. The subsequent Triassic continent-continent collision led to the final coalescence of the Yangtze and Sino-Korean cratons. Amalgamation of the Yangtze and North China cratons was, therefore, a multistage process extending over at least 200 Ma.     

Real-time envelope cross-correlation detector: application to induced seismicity in the Insheim and Landau deep geothermal reservoirs

We develop and test a real-time envelope cross-correlation detector for use in seismic response plans to mitigate hazard of induced seismicity. The incoming seismological data are cross-correlated in real-time with a set of previously recorded master events. For robustness against small changes in the earthquake source locations or in the focal mechanisms we cross-correlate the envelopes of the seismograms rather than the seismograms themselves. Two sequenced detection conditions are implemented: After passing a single trace cross-correlation condition, a network cross-correlation is calculated taking amplitude ratios between stations into account. Besides detecting the earthquake and assigning it to the respective reservoir, real-time magnitudes are important for seismic response plans. We estimate the magnitudes of induced microseismicity using the relative amplitudes between master event and detected event. The real-time detector is implemented as a SeisComP3 module. We carry out offline and online performance tests using seismic monitoring data of the Insheim and Landau geothermal power plants (Upper Rhine Graben, Germany), also including blasts from a nearby quarry. The comparison of the automatic real-time catalogue with a manually processed catalogue shows, that with the implemented parameters events are always correctly assigned to the respective reservoir (4 km distance between reservoirs) or the quarry (8 km and 10 km distance, respectively, from the reservoirs). The real-time catalogue achieves a magnitude of completeness around 0.0. Four per cent of the events assigned to the Insheim reservoir and zero per cent of the Landau events are misdetections. All wrong detections are local tectonic events, whereas none are caused by seismic noise.     

Extracting drainage networks and their connectivity using LiDAR data

Policies, measures, and models geared towards flood prevention and managing surface waters benefit from high quality data on the presence and characteristics of drainage ditches. As a cost and labour effective alternative for acquiring such data through field surveys, we propose a method (a) to extract vector data representing ditch drainage networks based on local morphologic features derived from high resolution digital elevation models (DEM) and (b) to identify possible connections in the ditch network by calculating a probability of the connectivity using a logistic regression where the predictor variables are characteristics of the ditch centre lines or derived from the DEM. Using Light Detection and Ranging (LiDAR) derived DEMs with a 1 m resolution, the method was developed and tested for a mixed agricultural residential area in north‐eastern Belgium. The derived ditch segments had an error of omission of 8% and an error of commission of 5%. The original positional accuracy of the centre lines of the extracted ditches was 0.6 m and could be improved to 0.4 m by shifting each vertex to the position of the lowest LiDAR point located within a radius equal to the spatial resolution of the used DEM. About 69% of the false disconnections in the network were identified and corrected leading to a reduction of the unconnected parts of the ditch network by 71%. The extracted and connected network approximated the reference ditch network fairly well.     

Gravity reduction with three-dimensional atmospheric pressure data for precise ground gravity measurements

The redistribution of air masses induces gravity variations (atmospheric pressure effect) up to about 20 μgal. These variations are disturbing signals in gravity records and they must be removed very carefully for detecting weak gravity signals. In the past, different methods have been developed for modelling of the atmospheric pressure effect. These methods use local or two-dimensional (2D) surface atmospheric pressure data and a standard height-dependent air density distribution. The atmospheric pressure effect is consisting of the elastic deformation and attraction term. The deformation term can be well modelled with 2D surface atmospheric pressure data, for instance with the Green's function method. For modelling of the attraction term, three-dimensional (3D) data are required. Results with 2D data are insufficient.From European Centre for Medium-Range Weather Forecasts (ECMWF) 3D atmospheric pressure data are now available. The ECMWF data used here are characterised by a spacing of Δ and Δλ = 0.5°, 60 pressure levels up to a height of 60 km and an interval of 6 h. These data are used for modelling of the atmospheric attraction term. Two attraction models have been developed based on the point mass attraction of air segments and the gravity potential of the air masses. The modelling shows a surface pressure-independent part of gravity variations induced by mass redistribution of the atmosphere in the order of some μgal. This part can only be determined by using 3D atmospheric pressure data. It has been calculated for the Vienna Superconducting Gravimeter site.From this follows that the gravity reduction can be improved by applying the 3D atmospheric attraction model for analysing long-periodic tidal waves including the polar tide. The same improvement is expected for reduction of long-term absolute gravity measurements or comparison of gravity measurements at different seasonal times. By using 3D atmospheric pressure data, the gravity correction can be improved up to some μgal.     

Combining long- and short-term probabilistic volcanic hazard assessment with cost-benefit analysis to support decision making in a volcanic crisis from the Auckland Volcanic Field,New Zealand

By using BET_VH, we propose a quantitative probabilistic hazard assessment for base surge impact in Auckland, New Zealand. Base surges resulting from phreatomagmatic eruptions are among the most dangerous phenomena likely to be associated with the initial phase of a future eruption in the Auckland Volcanic Field. The assessment is done both in the long-term and in a specific short-term case study, i.e. the simulated pre-eruptive unrest episode during Exercise Ruaumoko, a national civil defence exercise. The most important factors to account for are the uncertainties in the vent location (expected for a volcanic field) and in the run-out distance of base surges. Here, we propose a statistical model of base surge run-out distance based on deposits from past eruptions in Auckland and in analogous volcanoes. We then combine our hazard assessment with an analysis of the costs and benefits of evacuating people (on a 1 × 1-km cell grid). In addition to stressing the practical importance of a cost-benefit analysis in creating a bridge between volcanologists and decision makers, our study highlights some important points. First, in the Exercise Ruaumoko application, the evacuation call seems to be required as soon as the unrest phase is clear; additionally, the evacuation area is much larger than what is recommended in the current contingency plan. Secondly, the evacuation area changes in size with time, due to a reduction in the uncertainty in the vent location and increase in the probability of eruption. It is the tradeoff between these two factors that dictates which cells must be evacuated, and when, thus determining the ultimate size and shape of the area to be evacuated.     

An ecohydrological approach to predicting regional woody species distribution patterns in dryland ecosystems

This paper presents a quantitative ecohydrological framework for predicting regional distribution patterns of woody species in dryland ecosystems. The framework is based on an existing stochastic model for the daily mass balance of water that represents the interactions between soils, climate, and vegetation. Individual species selection is based on an optimality trade-off hypothesis, which states that dryland vegetation patterns are constrained by maximization of water use and simultaneous minimization of water stress. The relative importance of water use and stress avoidance to the overall fitness of three Acacia species is determined from the heterogeneous basin, the Upper Ewaso Ng’iro river basin, of the central Kenya highlands. The model results indicate that overall fitness is more strongly influenced by water use than stress avoidance but that consideration of both stress avoidance and water use is critical to predicting basin-scale patterns of species distribution. We identify a linear trend in the frequency and intensity of storms with the same annual total using a basin-wide gauge precipitation dataset. After calibration, we apply the basin average linear trends in time for average rain per storm and storm arrival rates. The model results indicate the upslope migration of two species, Acacia tortilis and Acacia xanthophloea to areas with higher total rainfall. Lastly, we explore the modeled changes of species cover in the basin influenced by changes in rainfall total holding growing season rainfall variability constant and changes in growing season rainfall variability holding total rainfall constant. We find that changes in dryland species distribution patterns and relative abundance may be as sensitive to growing season rainfall variability as they are to changes in total rainfall amounts.