Eocene-Pliocene time scale and stratigraphy of the Upper Rhine Graben (URG) and the Swiss Molasse Basin (SMB)

We present a general stratigraphic synthesis for the Upper Rhine Graben (URG) and the Swiss Molasse Basin (SMB) from Eocene to Pliocene times. The stratigraphic data were compiled both from literature and from research carried out by the authors during the past 6 years ; an index of the stratigraphically most important localitites is provided. We distinguish 14 geographical areas from the Helvetic domain in the South to the Hanau Basin in the North. For each geographical area, we give a synthesis of the biostratigraphy, lithofacies, and chronostratigraphic ranges. The relationships between this stratigraphic record and the global sea-level changes are generally disturbed by the geodynamic (e.g., subsidence) evolution of the basins. However, global sea-level changes probably affected the dynamic of transgression–regression in the URG (e.g., Middle Pechelbronn Beds and Serie Grise corresponding with sea-level rise between Ru1/Ru2 and Ru2/Ru3 sequences, respectively) as well as in the Molasse basin (regression of the UMM corresponding with the sea-level drop at the Ch1 sequence). The URGENT-project (Upper Rhine Graben evolution and neotectonics) provided an unique opportunity to carry out and present this synthesis. Discussions with scientists addressing sedimentology, tectonics, geophysics and geochemistry permitted the comparison of the sedimentary history and stratigraphy of the basin with processes controlling its geodynamic evolution. Data presented here back up the palaeogeographic reconstructions presented in a companion paper by the same authors (see Berger et al. in Int J Earth Sci 2005).     

Diurnal changes in the bio-optical properties of the phytoplankton in the Alborán Sea (Mediterranean Sea)

The changes in the phytoplankton absorption properties during a diurnal cycle were investigated at one station located in the north-western area of the Alborán Sea. The experiment was performed in spring when the water column was strongly stratified. This hydrological situation permitted the establishment of a deep chlorophyll a (chl a) fluorescence maximum (DFM) which was located on average close to the lower limit of the mixed layer and the nutricline. The relative abundance of pico-phytoplankton (estimated as its contribution to the total chl a) was higher in the surface, however, micro-phytoplankton dominated the community at the DFM level. Chl a specific absorption coefficient (a*(λ)) also varied with optical depth, with a* (the spectrally average specific absorption coefficient) decreasing by 30% at the DFM depth with respect to the surface. A significant negative correlation between the contribution of the micro-phytoplankton to the total chl a and a* was obtained indicating that a* reduction was due to changes in the packaging effect. Below the euphotic layer, a* increased three-fold with respect to the DFM, which agrees with the expected accumulation of accessory pigments relative to chl a as an acclimation response to the low available irradiance. The most conspicuous change during the diurnal cycle was produced in the euphotic layer where the chl a concentration decreased significantly in the afternoon (from a mean concentration of 1.1 μg L⁻¹ to 0.7 μg L⁻¹) and increased at dusk when it averaged 1.4 μg L⁻¹. In addition, a* and the blue-to-red absorption band ratio increased in the afternoon. These results suggest that a*(λ) diurnal variability was due to increase in photo-protective and accessory pigments relative to chl a. The variation ranges of a*(λ) at 675 and 440 nm (the absorption peaks in the red and blue spectral bands, respectively) in the euphotic layer were 0.01–0.04 and 0.02–0.10 m² mg⁻¹ chl a, respectively. Approximately 30% out of this variability can be attributed to the diurnal cycle. This factor should therefore be taken into account in refining primary production models based on phytoplankton light absorption.     

The effect of mountainous topography on moisture exchange between the “surface” and the free atmosphere

Typical numerical weather and climate prediction models apply parameterizations to describe the subgrid-scale exchange of moisture, heat and momentum between the surface and the free atmosphere. To a large degree, the underlying assumptions are based on empirical knowledge obtained from measurements in the atmospheric boundary layer over flat and homogeneous topography. It is, however, still unclear what happens if the topography is complex and steep. Not only is the applicability of classical turbulence schemes questionable in principle over such terrain, but mountains additionally induce vertical fluxes on the meso-γ scale. Examples are thermally or mechanically driven valley winds, which are neither resolved nor parameterized by climate models but nevertheless contribute to vertical exchange. Attempts to quantify these processes and to evaluate their impact on climate simulations have so far been scarce. Here, results from a case study in the Riviera Valley in southern Switzerland are presented. In previous work, measurements from the MAP-Riviera field campaign have been used to evaluate and configure a high-resolution large-eddy simulation code (ARPS). This model is here applied with a horizontal grid spacing of 350 m to detect and quantify the relevant exchange processes between the valley atmosphere (i.e. the ground “surface” in a coarse model) and the free atmosphere aloft. As an example, vertical export of moisture is evaluated for three fair-weather summer days. The simulations show that moisture exchange with the free atmosphere is indeed no longer governed by turbulent motions alone. Other mechanisms become important, such as mass export due to topographic narrowing or the interaction of thermally driven cross-valley circulations. Under certain atmospheric conditions, these topographical-related mechanisms exceed the “classical” turbulent contributions a coarse model would see by several times. The study shows that conventional subgrid-scale parameterizations can indeed be far off from reality if applied over complex topography, and that large-eddy simulations could provide a helpful tool for their improvement.     

Investigating parabolic and nebkha dune formation using a cellular automaton modelling approach

Vegetation plays an important role in shaping the morphology of aeolian dune landscapes in coastal and semi‐arid environments, where ecogeomorphic interactions are complex and not well quantified. We present a Discrete ECogeomorphic Aeolian Landscape model (DECAL) capable of simulating realistic looking vegetated dune forms, permitting exploration of relationships between ecological and morphological processes at different temporal and spatial scales. The cellular automaton algorithm applies three simple rules that lead to self‐organization of complex dune environments, including nebkhas with distinctive deposition tails that form in association with mesquite‐type shrubs, and hairpin (long‐walled) parabolic dunes with trailing ridges that evolve from blowouts in association with vegetation succession. Changing the conditions of simulations produces differing landscapes that conform qualitatively to observations of real‐world dunes. The model mimics the response of the morphology to changes in sediment supply, vegetation distribution, density and growth characteristics, as well as initial disturbances. The introduction of vegetation into the model links spatial and temporal scales, previously dimensionless in bare‐sand cellular automata. Grid resolutions coarser than the representative size of the modelled vegetation elements yield similar morphologies, but when cell size is reduced to much smaller dimensions, the resultant landscape evolution is dramatically different. The model furthermore demonstrates that the relative response characteristics of the multiple vegetation types and their mutual feedback with geomorphological processes impart a significant influence on landscape equilibria, suggesting that vegetation induces a characteristic length scale in aeolian environments. This simple vegetated dune model illustrates the power and versatility of a cellular automaton approach for exploring the effects of interactions between ecology and geomorphology in complex earth surface systems. Copyright © 2007 John Wiley & Sons, Ltd.     

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.     

The constancy of galactic cosmic rays as recorded by cosmogenic nuclides in iron meteorites

We measured the He, Ne, and Ar isotopic concentrations and the ¹⁰Be, ²⁶Al, ³⁶Cl, and ⁴¹Ca concentrations in 56 iron meteorites of groups IIIAB, IIAB, IVA, IC, IIA, IIB, and one ungrouped. From ⁴¹Ca and ³⁶Cl data, we calculated terrestrial ages indistinguishable from zero for six samples, indicating recent falls, up to 562 ± 86 ka. Three of the studied meteorites are falls. The data for the other 47 irons confirm that terrestrial ages for iron meteorites can be as long as a few hundred thousand years even in relatively humid conditions. The ³⁶Cl‐³⁶Ar cosmic ray exposure (CRE) ages range from 4.3 ± 0.4 Ma to 652 ± 99 Ma. By including literature data, we established a consistent and reliable CRE age database for 67 iron meteorites. The high quality of the CRE ages enables us to study structures in the CRE age histogram more reliably. At first sight, the CRE age histogram shows peaks at about 400 and 630 Ma. After correction for pairing, the updated CRE age histogram comprises 41 individual samples and shows no indications of temporal periodicity, especially not if one considers each iron meteorite group separately. Our study contradicts the hypothesis of periodic GCR intensity variations (Shaviv 2002, 2003), confirming other studies indicating that there are no periodic structures in the CRE age histogram (e.g., Rahmstorf et al. 2004; Jahnke 2005). The data contradict the hypothesis that periodic GCR intensity variations might have triggered periodic Earth climate changes. The ³⁶Cl‐³⁶Ar CRE ages are on average 40% lower than the ⁴¹K‐K CRE ages (e.g., Voshage 1967). This offset can either be due to an offset in the ⁴¹K‐K dating system or due to a significantly lower GCR intensity in the time interval 195–656 Ma compared to the recent past. A 40% lower GCR intensity, however, would have increased the Earth temperature by up to 2 °C, which seems unrealistic and leaves an ill‐defined ⁴¹K‐K CRE age system the most likely explanation. Finally, we present new ²⁶Al/²¹Ne and ¹⁰Be/²¹Ne production rate ratios of 0.32 ± 0.01 and 0.44 ± 0.03, respectively.     

Probability of earthquake occurrence and magnitude estimation in the post shut-in phase of geothermal projects

Induced seismicity in geothermal projects is observed to continue after shut-in of the fluid injection. Recent experiments show that the largest events tend to occur after the termination of injection. We use a probabilistic approach based on Omori??s law and the Gutenberg?CRichter magnitude?Cfrequency distribution to demonstrate that the probability of exceeding a certain maximum magnitude still increases after shut-in. This increase is governed by the exponent of Omori??s law q and the Gutenberg?CRichter b value. For a reduced b value in the post-injection phase, the probability of occurrence directly after shut-in can be even higher than the corresponding probability for an ongoing injection. For the reference case of q?=?2 and a 10% probability at shut-in time t _S to exceed a given maximum magnitude, we obtain an increase to 14.6% for t?=?2t _S at a constant Gutenberg?CRichter b value after shut-in. A reduction of the b value by one quarter leads to a probability of 20.5%. If we consider a constant probability level of occurrence for an event larger than a given magnitude at shut-in time, this maximum magnitude increases by 0.12?units for t?=?2t _S (0.26?units for a reduced b value). For the Soultz-sous-Forêts (France) injection experiment in 2000, recent studies reveal q?=?9.5 and a b value reduction by 14%. A magnitude 2.3 event 9?h after shut-in falls in the phase with a probability higher than for the continued injection. The probability of exceeding the magnitude of this post-injection event is determined to 97.1%.     

The Younger Dryas in the English Lake District: reconciling geomorphological evidence with numerical model outputs

The geomorphology of the south‐western and central Lake District, England is used to reconstruct the mountain palaeoglaciology pertaining to the Lateglacial and Younger Dryas. Limitations to previous ice‐mass reconstructions and consequent palaeoclimatic inferences include: (i) the use of static (steady‐state) glacier reconstructions, (ii) the assumption of a single‐stage Younger Dryas advance, (iii) greatly varying ice‐volume estimates, (iv) inexplicable spatial variations in ELA (Equilibrium Line Altitude), and (v) a lack of robust extent chronology. Here we present geomorphological mapping based on aerial photography and the NextMap Britain Digital Elevation Model, checked by ground survey. Former glacier extents were inferred and ELAs were calculated using the Balance Ratio method of Osmaston. Independently, a time‐dependant 2‐D ice‐flow model was forced by a regional ELA history that was scaled to the GRIP record. This provided a dynamic reconstruction of a mountain ice field that allowed for non‐steady‐state glacier evolution. Fluctuations in climate during the Younger Dryas resulted in multiple glacial advance positions that show agreement with the location of mapped moraines, and may further explain some of the ELA variations found in previous local and static reconstructions. Modelling based on the GRIP record predicts three phases: an initial maximum extent, a middle minor advance or stillstand, and a pronounced but less extensive final advance. The comparisons find that the reconstructions derived from geomorphological evidence are effective representations of steady‐state glacier geometries, but we do propose different extents for some glaciers and, in particular, a large former glacier in Upper Eskdale.