Transfer and early diagenesis of biogenic silica oxygen isotope signals during settling and sedimentation of diatoms in a temperate freshwater lake (Lake Holzmaar, Germany)

We have investigated the transfer of oxygen isotope signals of diatomaceous silica (δ¹⁸O_diatom) from the epilimnion (0-7 m) through the hypolimnion to the lake bottom (∼20 m) in freshwater Lake Holzmaar, Germany. Sediment-traps were deployed in 2001 at depths of 7 and 16 m to harvest fresh diatoms every 28 days. The 7 m trap collected diatoms from the epilimnion being the main zone of primary production, while the 16 m trap collected material already settled through the hypolimnion. Also a bottom sediment sample was taken containing diatom frustules from approximately the last 25 years. The δ¹⁸O_diatom values of the 7 m trap varied from 29.4‰ in spring/autumn to 26.2‰ in summer according to the temperature dependence of oxygen isotope fractionation and represent the initial isotope signal in this study. Remarkably, despite the short settling distance δ¹⁸O_diatom values of the 7 and the 16 m trap were identical only during spring and autumn seasons while from April to September δ¹⁸O_diatom values of the 16 m trap were roughly ∼1.5‰ enriched in ¹⁸O compared to those of the 7 m trap. Isotopic exchange with the isotopically lighter water of the hypolimnion would shift the δ¹⁸O_diatom value to lower values during settling from 7 to 16 m excluding this process as a cause for the deviation. Dissolution of opal during settling with intact organic coatings of the diatom cells and near neutral pH of the water should only cause a minor enrichment of the 16 m values. Nevertheless, opal from the bottom sediment was found to be 2.5‰ enriched in ¹⁸O compared to the weighted average of the opal from the 7 m trap. Thus, resuspension of bottom material must have contributed to the intermediate δ¹⁸O_diatom signal of the 16 m trap during summer. Dissolution experiments allowed further investigation of the cause for the remarkably enriched δ¹⁸O_diatom value of the bottom sediment. Experiments with different fresh diatomaceous materials show an increase of opaline ¹⁸O at high pH values which is remarkably reduced when organic coatings of the cells still exist or at near neutral pH. In contrast, high pH conditions do not affect the δ¹⁸O_diatom values of sub-fossil and even fossil opal. IR analyses show that the ¹⁸O enrichment of the sedimentary silica is associated with a decrease in Si-OH groups and the formation of Si-O-Si linkages. This indicates a silica dehydroxylation process as cause for the isotopic enrichment of the bottom sediment. Silica dissolution and dehydroxylation clearly induce a maturation process of the diatom oxygen isotope signal presumably following an exponential behaviour with a rapid initial phase of signal alteration. The dynamics of this process is of particular importance for the quantitative interpretation of sedimentary δ¹⁸O_diatom values in terms of palaeothermometry.     

Ductile deformation of garnet in mylonitic gneisses from the Münchberg Massif (Germany)

Mylonitic gneisses from the Münchberg Massif contain single grains (type I) and polycrystalline aggregates (type II) of garnet displaying a distinct elongation parallel to a macroscopic lineation which is interpreted as the result of ductile deformation. Lattice-preferred orientations of quartz (textures) symmetrical to the macroscopic foliation and lineation and the lack of rotational microfabrics indicate that the bulk deformation was pure shear at least during the latest strain increments. Garnet textures measured by EBSD together with microprobe analyses demonstrate that these two structural types of garnet can be related to two different processes of ductile deformation: (1) For the single grains stretching can be attributed to diffusion creep along grain boundary zones (Coble creep). The related mass transfer is indicated by the fact that primary growth zones are cut off at the long faces of the grains while the related strain shadow domains do not show comparable chemical zoning. Pressure solution and precipitation suitable to produce similar structures can be largely ruled out because retrogressive reactions pointing to the presence of free hydrous fluids are missing. (2) For the polycrystalline garnet aggregates consisting of cores grading into fine-grained mantles, dislocation creep and associated rotation recrystallization can be assumed. Continuous lattice rotation from the core to the outer polycrystalline rim allow a determination of the related dominant slip systems which are {100}<010> and equivalent systems according to the cubic lattice symmetry. The same holds for garnets which appear to be completely recrystallized. For this type of fine-grained aggregates an alternative nucleation model is discussed. Due to penetrative dislocation glide in connection with short range diffusion and the resulting lattice rotation, primary growth zones are strongly disturbed.Since for the considered rock unit of the Münchberg Massif peak metamorphic temperatures between 630 and 670 °C can be assumed, this study clearly demonstrates that the inferred processes of ductile garnet deformation can occur not only in HT regimes as often suggested in the literature even if embedded within a matrix of “low-strength” minerals like quartz, feldspars and micas.     

Early Jurassic rift structures associated with the Soapaga and Boyacá faults of the Eastern Cordillera, Colombia: Sedimentological inferences and regional implications

The NW-trending Bucaramanga fault links, at its southern termination, with the Soapaga and Boyacá faults, which by their NW trend define an ample horsetail structure. As a result of their Neogene reactivation as reverse faults, they bound fault-related anticlines that expose the sedimentary fill of two Early Jurassic rift basins. These sediments exhibit the wedge-like geometry of rift fills related to west-facing normal faults. Their structural setting was controlled further by segmentation of the bounding faults at approximately 10 km intervals, in which each segment is separated by a transverse basement high. Isopach contours and different facies associations suggest these transverse anticlines may have separated depocenters of their adjacent subbasins, which were shaped by a slightly different subsidence history and thereby decoupled. The basin fill of the relatively narrow basin associated with the Soapaga fault is dominated by fanglomeratic successions organized in two coarsening-upward cycles. In the larger basin linked to the Boyacá fault, the sedimentary fill consists of two coarsening-upward sequences that, when fully developed, vary from floodplain to alluvial fan deposits. These Early Jurassic rift fills temporally constrain the evolution of the Bucaramanga fault, which accommodated right-lateral displacement during the early Mesozoic rift event.     

Assessing the Vibrational Frequencies of the Cathedral of Cologne (Germany) by Means of Ambient Seismic Noise Analysis

Ambient seismic noise measurements were conducted inside the Cathedral of Cologne (Germany) for assessing its frequencies of vibration and for checking whether these occur in the range where soil amplification is expected. If this is the case, damages may increase in case of an earthquake due to an increased structural response of the building. Analysis of the ratio between the horizontal and vertical components of the spectra recorded at stations located inside the building as well as the ratio between the corresponding components of the spectra recorded simultaneously inside the building and at a reference station placed in the basement of the cathedral indicated several modes of vibration. Facilitated by these results an assessment of the seismic vulnerability was attempted for a 2D ground motion scenario using the finite element method.     

Tracing leachates at waste sites using geophysical and geochemical modelling

Tracing leachates at landfills is usually carried out using either geophysical methods or chemical analyses of groundwater. There are often problems with fingerprinting pollution sources or clarifying the spreading pattern due to a wide variety of possibilities giving similar anomalies. The aims of the project were to evaluate the advantages of combining results from multigeophysical modelling and statistical/chemical modelling in order to identify pollution sources and the spreading pattern and to test a new technique for chemical fingerprinting. The project was carried out at a landfill in central Sweden using geophysical measurements and modelling of CVES, GPR and VLF as well as chemical modelling using M3 (multivariate mixing and mass balance calculations). The results indicate that by combining geophysical modelling and chemical calculations, the possibilities of fingerprinting the origin of pollution as well as delineating the spreading pattern are significantly increased.     

Early Holocene water budget of the Nakuru-Elmenteita basin, Central Kenya Rift

The Nakuru-Elmenteita basin in the Central Kenya Rift, contains two shallow, alkaline lakes, Lake Nakuru (1770 m above sea level) and Lake Elmenteita (1786 m). Ancient shorelines and lake sediments at 1940 m suggest that these two lakes formed a single large and deep lake as a result of a wetter climate during the early Holocene. Here, we used a hydrological model to compare the precipitation–evaporation balance during the early Holocene to today. Assuming that the Nakuru-Elmenteita basin was hydrologically closed, as it is today, the most likely climate scenario includes a 45% increase in mean-annual precipitation, a 0.5°C decrease in air temperature, and an increase of 9% in cloud coverage from the modern values. Compared to the modeling results from other East African lake basins, this dramatic increase in precipitation seems to be unrealistic. Therefore, we propose a significant flow of water from the early Holocene Lake Naivasha in the south towards the Nakuru-Elmenteita basin to compensate the extremely negative hydrological budget of this basin. Since we did not find any field evidence for a surface connection, as often proposed during the last 70 years, the hydrological deficit of the Nakuru-Elmenteita basin could have also been compensated by a subsurface water exchange.     

Strike-slip fault bridge fluid pumping mechanism: insights from field-based palaeostress analysis and numerical modelling

We present a finite-element study of stress perturbation in evolving compressive and extensional strike-slip fault bridges. The results are compared with a fracture study of a compressive bridge at St Donats, UK. Horizontally interbedded calcareous mudstone and bioclastic calcilutite at St Donats have a distinct vertical permeability anisotropy. This sedimentary sequence behaves as a set of horizontal aquifers. The fluid flow in these aquifers is sensitive to mean stress gradients. Paleostress analysis of field fracture data, verified by finite-element modelling, indicates a rotation of σ₁ towards parallelism with boundary faults inside the growing compressive bridge. Boundary faults and bridge faults recorded numerous fluid flow events. The modelled mean stress pattern shows a regional maximum within the bridge and local maxima/minima pairs at boundary fault tips.Finite-element modelling of an extensional bridge indicates that σ₃ rotates towards parallelism with boundary faults. The mean stress pattern is similar to the pattern in compressive bridge but with maxima and minima locations interchanged. The stress patterns are reestablished by each stress build-up preceding the rupturation of the boundary faults throughout the development stages of strike-slip fault bridges. Mean stress gradients developed pre-failure control the fluid flow in fractures of the strike-slip fault system at and after the end of each stress build-up and the fluid flow in boundary faults post-failure. Fracture reactivation and new fracture generation within an evolving bridge is a process consisting of multiple successive events that retain the storage capacity of the bridge. Rupture and sealing of the main bounding-faults is a step-wise process that opens and closes fluid conduits between areas with different pressures.     

Probability Density Functions as Botanical-Climatological Transfer Functions for Climate Reconstruction

We present a new procedure, the pdf method (pdf=probability density function), for reconstructing Quaternary climate utilizing botanical data. The procedure includes the advantages of the indicator species method by considering the fossil and modern presence and absence of taxa rather than their frequencies, thus avoiding the need for modern analog plant communities. Overcoming the problematic use of absolute limits to describe climate response ranges is the main progress of the pdf method in comparison to the indicator species method. This advantage results from estimating probability density functions (pdfs) for monthly mean January and July temperature conditional on the present day occurrence of single taxa. Gaussian distributions sufficiently approximate pdfs of many, although not all, studied taxa. On the assumption of statistical independence, the procedure calculates a joint pdf as the product of the pdfs of the individual taxa. This algorithm weights each taxon according to the extent of its climate response range expressed by its covariance structure. We interpret the maximum of the resulting pdf as the most likely climate and its confidence interval as the uncertainty range. To avoid an artificial reduction of uncertainty arising from the use of numerous similar pdfs, a preselection method is proposed based on the Mahalanobis distance between pdfs. The pdf method was applied to the Carpinus phase of a profile from Gröbern, Germany, that spans the last interglaciation (Eemian). The reconstructed most probable January and July temperatures of about 0.0°C and 18.4°C barely differ from the modern values of −0.5°C and 18.3°C.     

Inelastic dynamic analysis of RC bridges accounting for spatial variability of ground motion,site effects and soil–structure interaction phenomena. Part 1: Methodology and analytical tools

During strong ground motion it is expected that extended structures (such as bridges) are subjected to excitation that varies along their longitudinal axis in terms of arrival time, amplitude and frequency content, a fact primarily attributed to the wave passage effect, the loss of coherency and the role of local site conditions. Furthermore, the foundation interacts with the soil and the superstructure, thus significantly affecting the dynamic response of the bridge. A general methodology is therefore set up and implemented into a computer code for deriving sets of appropriately modified time histories and spring–dashpot coefficients at each support of a bridge with account for spatial variability, local site conditions and soil–foundation–superstructure interaction, for the purposes of inelastic dynamic analysis of RC bridges. In order to validate the methodology and code developed, each stage of the proposed procedure is verified using recorded data, finite‐element analyses, alternative computer programs, previous research studies, and closed‐form solutions wherever available. The results establish an adequate degree of confidence in the use of the proposed methodology and code in further parametric analyses and seismic design. Copyright © 2003 John Wiley & Sons, Ltd.