The Somme valley terrace system (northern France); a model of river response to Quaternary climatic variations since 800,000 bp

Recent research on the Somme Valley terrace deposits, that began in 1987, has resulted in an integrated sequence of geological and geomorphological evidence for this area that is of particular significance for prehistoric investigations. This has been achieved using previous evidence and new results from archaeological and geological excavations. These include longitudinal profiles, cartography and detailed synthesis of sections of alluvial deposits, stratigraphy of river deposits and overlying loess sequences. The work, based on a study of the geometry of the terrace system, has resulted in the identification of ninestepped alluvial formations that occur along 70 km of the valley. Detailed study permits the integration of the various alluvial sequences. This suggests that each alluvial sequence results from sedimentation through a single glacial-interglacial climatic cycle. This hypothesis is confirmed by investigation of soil and loess sequences of the deposits overlying the fluvial sediments. The synthesis presented provides a basis for understanding the deposition of the terrace system and an initial chronostratigraphic interpretation for the deposits.     

Exponential time integrators for stochastic partial differential equations in 3D reservoir simulation

The transport of chemically reactive solutes (e.g. surfactants, CO₂ or dissolved minerals) is of fundamental importance to a wide range of applications in oil and gas reservoirs such as enhanced oil recovery and mineral scale formation. In this work, we investigate exponential time integrators, in conjunction with an upwind weighted finite volume discretisation in space, for the efficient and accurate simulation of advection–dispersion processes including non-linear chemical reactions in highly heterogeneous 3D oil reservoirs. We model sub-grid fluctuations in transport velocities and uncertainty in the reaction term by writing the advection–dispersion–reaction equation as a stochastic partial differential equation with multiplicative noise. The exponential integrators are based on the variation of constants solution and solve the linear system exactly. While this is at the expense of computing the exponential of the stiff matrix representing the finite volume discretisation, the use of real Léja point or the Krylov subspace technique to approximate the exponential makes these methods competitive compared to standard finite difference-based time integrators. For the deterministic system, we investigate two exponential time integrators, the second-order accurate exponential Euler midpoint (EEM) scheme and exponential time differencing of order one (ETD1). All our numerical examples demonstrate that our methods can compete in terms of efficiency and accuracy compared with standard first-order semi-implicit time integrators when solving (stochastic) partial differential equations that model mixing and chemical reactions in 3D heterogeneous porous media. Our results suggest that exponential time integrators such as the ETD1 and EEM schemes could be applied to typical 3D reservoir models comprising tens to hundreds of thousands unknowns.     

Rapid and cyclic aeolian deposition during the Last Glacial in European loess: a high-resolution record from Nussloch, Germany

This paper reports the results of an investigation of the Weichselian Upper Pleniglacial loess sequences of Nussloch (Rhine Valley, Germany) based on stratigraphy, palaeopedology, sedimentology, palynology, malacology and geochemistry (δ¹³C), supported by radiocarbon, TL and OSL dating. Grain-size and magnetic susceptibility records are taken at 5 cm intervals from the Upper Pleniglacial (UPG) loess. The data indicate cyclic variations in loess deposition between ca 34 and 17 ka, when the sedimentation rate is especially high (1.0–1.2 m per ka for more than 10 m). The grain-size index (GSI: ratio of coarse silt versus fine silt and clay) shows variations, which are assumed to be an indirect measurement of wind intensity. The sedimentation rate, interpreted from the profiles, indicates high values in loess (Loess events LE-1 to LE-7) and low or negligible values in tundra gley horizons G1 to G8. OSL ages from the loess and ¹⁴C dates from organic matter in the loess show that loess deposition was rapid but was interrupted by shorter periods of reduced aeolian sedimentation. Comparison between the data from Nussloch and other European sequences demonstrates a progressive coarsening of the loess deposits between ca 30 and 22 ka. This coarsening trend ends with a short but major grain-size decrease and is followed by an increase to a new maximum at 20 ± 2 ka (“W” shape). Correlation between the loess GSI and the Greenland ice-core dust records, suggests a global connection between North Atlantic and Western European global atmospheric circulation and wind regimes. In addition, the typical Upper Pleniglacial loess deposition begins at ca 30–31 ka, close to Heinrich event (HE) 3, and the main period of loess sedimentation at about 25 ± 2 ka is coeval to HE 2. Correlation of magnetic susceptibility and grain-size records shows that the periods, characterised by high GSI, coincide with an increase in the amount of ferromagnetic minerals reworked from the Rhine alluvial plain. They suggest enhancement in the frequency of the storms from N–NW. These results are integrated within a palaeogeographical model of dust transport and deposition in Western Europe for the Weichselian Upper Pleniglacial (or Late Pleniglacial).     

Processing and interpretation of the gravity field of the East African Rift: implication for crustal extension

Compilation of new and existing gravity data were undertaken to assess the nature of the crust beneath the East African Rift System. Using 3D gravity modeling code crustal model of gravity profiles across two sectors of the rift were computed. The results are discussed in light of the structure of the rift system.The results of the 3D modeling of gravity profiles across the two rift zones revealed northward thinning of the crust. The maximum crustal attenuation occurs beneath the Afar depression, indicating the Afar rift undergoes an intense fragmentation of the crust resulting from faulting and magmatic activity. However, our computed crustal thickness below the Afar depression falls within an upper bound compared to elsewhere below tectonically active rift zones. This can be explained in terms of crustal accretion resulting from an impact of the Afar mantle plume since 30 Ma ago.The residual gravity obtained using high-cut filtering techniques reveals significant density contrast between the northern and southern sectors of the rift. The northern part of the rift is characterized by regular patterns of positive gravity anomalies, which can be interpreted in terms of a zone of crustal thinning through which relatively dense materials have intruded the overlying crust. In contrast, south of the Main Ethiopian Rift, the anomalies are characterized by random patterns and low amplitudes. The along-rift-axis variation in gravity anomalies implies that the style of crustal deformation changed progressively, beginning with regionally distributed crustal deformation, such as the one we observe within the more juvenile and wider southern segment of the rift, to localized deformation within the active and narrow rift zones of the northern sector of the Ethiopian Rift. We suggest that the key parameters controlling along-rift-axis variation in gravity anomalies are the rate of crustal extension, faulting and magmatic activities.     

Precambrian crystalline basement in southern Mongolia as revealed by SHRIMP zircon dating

Single zircon ages determined by ion microprobe (SHRIMP II) for granitoid gneisses from the southern slope of the Baga Bogd massif (Gobi-Altai, southern Mongolia) reveal several episodes of zircon growth, ranging from late Palaeoproterozoic to late Cambrian. The oldest events are documented by a zircon crystallization age for a gneiss protolith at 1519 ± 11 Ma and by a xenocrystic zircon from a dark grey augen-gneiss yielding an age of c. 1701 Ma. Discrete igneous events are recorded in granite-gneisses with protolith emplacement ages of 983 ± 6, 956 ± 3 and 954 ± 8 Ma. These ages provide the first record of early Neoproterozoic magmatic activity in this region. A much younger and discrete magmatic event is recorded by several dioritic to granitic orthogneisses which are tectonically interlayered with the older gneisses and have protolith emplacement ages between 502 and 498 Ma. These late Cambrian granitoids of calc-alkaline affinity are likely to have been emplaced along an active continental margin and suggest that the Baga Bogd Precambrian crustal fragment was either docked against the southward (present-day coordinates) growing margin of the CAOB or was a large enough crustal entity to develop an arc along its margin. We speculate that the Precambrian gneisses of this massif may be part of a crustal fragment rifted off the Tarim Craton.     

Atténuation de l'intensité macrosismique pour la France métropolitaine : importance de l'intensité épicentrale

The main result of this work is to show that macroseismic intensity decay with distance strongly depends on the epicentral intensity. An attenuation law that takes this parameter into account is proposed for Metropolitan France, from the analysis of SISFRANCE macroseismic database. Such a model significantly reduces the difference between observed and theoretical intensities. A map of the attenuation variations is also set up for Metropolitan France. No major site effects are observed, but at a broad scale, young Alpine regions display a stronger attenuation than old Hercynian regions. To cite this article: P. Arroucau et al., C. R. Geoscience 338 (2006).     

Assessment of the ecological quality status of soft-bottoms in Reunion Island (tropical Southwest Indian Ocean) using AZTI marine biotic indices

The ability of the two synthetic marine biotic indices, AMBI and M-AMBI, to account for changes in the ecological quality of coastal soft bottoms of Reunion Island according to disturbances was assessed from macrobenthic samples collected in five sectors between 1994 and 2004. Samples were collected under non-perturbed conditions and at two sites subjected to heavy organic enrichment. Both indices are based on a classification of macrofauna into ecological groups (EG), and their transfer to tropical waters required some adaptations. These indices proved efficient in detecting a degradation of habitat quality. Their use resulted in the classification of all sites sampled between 1996 and 1998 as "good" or "high". M-AMBI nevertheless tended to result in the attribution of a slightly worse ecological quality status than AMBI. Together with an update of the EG species list for the Indian Ocean area, our results support the extension of both indices for the assessment of tropical soft bottoms.     

Etude des Déformations et Inclinaisons du Sous-sol: Développements Instrumentaux et Applications dans la Région Sismique d’Arette,France

We present a modified version of a tilt and strain meter conceived to have great sensibility and good stability. These instruments were deployed in a cave at the seismic Pyrenean region of Arette (France). We set in evidence several types of phenomena, of which deformations brought about by water recharge filling up fissures on the massif are of importance. Owing to the redundance in the data, the study of the co-seismic deformations allows one to separate instrumental displacements and local readjustments from regional motions.     

Colonisation of epilithic diatom population of the River Taff, South Wales, U.K.

1. A light and scanning electron microscopic study of the structure of epilithic diatom communities was carried out at two sites on the River Taff, South Wales, an upstream unpolluted site (Site 1) and a downstream polluted site (Site 2). 2. Community structure at the micro-scale revealed that, although communities at both sites were similar in terms of species present (similarity indices ranging between 0.44 and 0.96), they were very different with respect to their physical structure. 3. The traditional approach in studying epilithic diatom communities ignores information relating to community structure. The combination of light and scanning electron microscopy proved to be most useful for studying the epilithic diatom communities. 4. The examination of community succession on stone substrata showed that communities were highly dynamic on a weekly time scale, which suggests that the “age” of a substratum is a very important determinant of the attaching communities. A change from a prostate, two-dimensional community structure to an erect, three-dimensional community structure, and a shift from the dominance of epilithic diatoms to the dominance of epiphytic diatoms were noted. Such a change was more noticed at the polluted downstream site, possibly due to a higher rate of surface pre-conditioning and biofilm production at this site than at the unpolluted upstream site.     

Analysis and interpretation of the October 21, 1766 earthquake in the Southeastern Caribbean

The October 21, 1766 earthquake is the most widely felt event in the seismic history of Trinidad and Venezuela. Previous works diverged on the interpretation of the historical data available for this event. They associated the earthquake either with the Lesser Antilles subduction zone, with strike-slip motion along El Pilar fault, or with intraplate deformation at the edge of Guyana shield. Isoseismal areas are proposed after a new search and analysis of primary and secondary sources of historical information. Two of the largest earthquakes of the twentieth century which occurred in the region, the 1968 (M _S 6.4, h = 103 km), and the 1997 (M _W 6.9, h = 25 km) events, for which both intensity data and instrumentally determined source parameters are available, are used to calibrate the isoseismal areas and to interpret them in terms of source depth and magnitude. It is concluded that the large extent of intensity values higher than V is diagnostic of the depth (85 ± 20 km) of the 1766 source, and of local amplifications of ground motion due to soft soil conditions and to strong contrasts of impedance at the edge of Guyana shield. It is proposed that the event occurred either in slab, or close to the bottom lithospheric interface between the Caribbean and South American plates (∼11°N; ∼62.5°W). The value of the magnitude is estimated at 6.5 < M _S < 7.5 depending on the source depth and on the decay of ground motion as a function of distance. Deep and intermediate depth earthquakes can induce important casualties in Trinidad, Venezuela, and Guyana, possibly more damaging than those induced by shallower earthquakes along the strike of El Pilar Fault.