Coleopteran evidence for a mosaic of environments at high altitude in the eastern Pyrénées,France, during the climatic transition between the Allerød and Younger Dryas

Late‐glacial environmental and climatic implications are inferred from an insect fauna from organic sediments infilling a palaeochannel on the banks of the River Têt, eastern Pyrénées, France. A pine cone in association with the insect fauna has been radiocarbon dated to 10 920 ± 60 yr BP, namely close to the Allerød – Younger Dryas boundary. Two distinct insect associations appear to be recognisable here. One is an assemblage typical of the high altitude forest and a second is characteristic of an alpine grassland. The close coexistence of these two assemblages is attributed to the climatic cooling towards the start of the Younger Dryas Stadial, when the forest cover broke up into remnant patches interspersed by alpine grassland. It is suggested that in a region of such high relief a mosaic of habitats may have been caused by patchy differences in insolation aspect, especially during a period of climatic deterioration. Copyright © 1999 John Wiley & Sons, Ltd.     

Hydrothermal alterations in the Echassières granitic cupola (Massif central,france)

Detailed petrographic and mineralogic investigations of an albite-lepidolite granite at Echassières (Massif Central, France; scientific deep drill program) shows the existence of hydrothermal stages which are closely related to the magmatic and structural history. According to fluid inclusion data, K-Ar datations and ¹⁸O/¹⁶O-D/H compositions of secondary minerals, two successive hydrothermal periods have been recognized. The early one (273–268 million years) produced a series of aluminous phyllosilicates: muscovite, pyrophyllite, donbassite, tosudite, kaolinite which are observed as vein deposits (<10 mm wide) and alteration products of primary minerals in wall-rocks. The vein system was sealed by monomineralic assemblages during a cooling period (400–150°C). This early hydrothermal alteration stage was controlled by interactions of rock with low salinity (1–10 wt% NaCl equivalent) fluids expelled from the granitic body during the cooling processes. The chemical properties of these fluids were the following: low pH, very low Mg and Fe and high Li, Na and K contents. Thermodynamic calculations show that the sequence pyrophyllite, Li-bearing donbassite, tosudite is mostly temperature dependent. From the chemical composition of secondary minerals and isotopic data it can be deduced that these fluids, which have a meteoric origin, have been expelled from the granite body during its cooling period and after interaction with it at high temperature. The late hydrothermal stage corresponds to deposits of fluorite and Fe-Mg rich illite (151 million years) in subvertical fractures. Temperature conditions did not exceed 250° C and fluids came through the surrounding metamorphic rocks into the granitic body. IIlite/smectite mixed-layer minerals have been identified in subvertical fractures which were opened during Tertiary periods. In the host micaschists, successive hydrothermal alterations took place during the cooling of the Beauvoir granite. Early magmatic fluids interacted with these micaschists. Locally, the metamorphic assemblage is replaced by a metasomatic one. Secondary topaz and (F, Li)-rich mica crystals were formed over a range of 450 of 150°C. Later hydrothermal fluids reacted with the country rocks to form phengite-biotite, chlorite-illite and kaolinite over a range of 300 to 150°C. Illite/smectite mixed-layer minerals crystallized in the roof micaschists and within the Beauvoir granite during the Tertiary alteration period. Meteoric water invaded open fractures producing supergene alteration mineral assemblages.     

Identification of Groundwater Parameters Using an Adaptative Multiscale Method

The identification of groundwater parameters in heterogeneous systems is a major challenge in groundwater modeling. Flexible parameterization methods are needed to assess the complexity of the spatial distributions of these parameters in real aquifers. In this article, we introduce an adaptative parameterization to identify the distribution of hydraulic conductivity within the large‐scale (4400 km²) Upper Rhine aquifer. The method is based on adaptative multiscale triangulation (AMT) coupled with an inverse problem procedure that identifies the parameters' distributions by reducing the error between measured and simulated heads. The AMT method has the advantage of combining both zonation and interpolation approaches. The AMT method uses area‐based interpolation rather than an interpolation based on stochastic features. The method is applied to a standard 2D groundwater model that takes into account the interactions between the aquifer and surface water bodies, groundwater recharge, and pumping wells. The simulation period covers 204 months, from January 1986 to December 2002. Recordings at 109 piezometers are used for model calibration. The simulated heads are globally quite accurate and reproduce the main dynamics of the system. The local hydraulic conductivities resulting from the AMT method agree qualitatively with existing local experimental observations across the Rhine aquifer.     

阿尔金断裂西段首期GPS观测资料处理及初步结果

本文介绍了中国科学院青藏高原研究所在阿尔金断裂西段进行GPS观测的剖面布设和首期观测情况.利用GAMIT/GLOBK软件对首期观测数据进行了处理,并将结果归算到区域参考框架下.GAMIT计算得到的"准观测值"均满足所须精度,其中单天解标准化均方差小于0.3,基线的平均相对精度为1.8×10-9.GLOBK在区域参考框架...     

Controlling factors on the morphology and spatial distribution of methane-related tubular concretions – Case study of an Early Eocene seep system

Up to 10 m in length and >1 m in diameter tubular, calcite-cemented sandstone concretions are hosted by the faulted Dikilitash unconsolidated sands and sandstones. These structures document shallow subsurface pathways of Early Eocene methane seepage in the Balkan Mountains foreland (NE Bulgaria). Their exceptional exposure allowed a unique study of the factors governing the morphology and spatial distribution of such fossilized fluid conduits. The large dimensions and subvertical, cylindrical shape of the most common tube type primarily reflects the buoyancy-driven, vertical path of an ascending gas-bearing fluid through permeable, mainly unconsolidated sandy host sediments. Tube morphology was also influenced by local stratigraphic anisotropies and might as well document differences in former seepage conditions. Mapping of >800 tubular concretions showed the NNW–NNE elongation and alignment of tube clusters and massive cemented sandstone structures. This suggests that Paleogene fault systems played a major role in directing the movement of fluids. However, within a single tube cluster, tubes are preferentially aligned, over distances up to 50 m along directions at an angle between 10° and 36° with respect to the inferred NNW–NNE, cluster parallel fault traces. In addition, cylindrical tubes of analogue dimensions are aligned over distances >100 m along N15° to N25°-oriented directions. It is hypothesized that this spatial geometry of tubular concretions reflects the complex geometry of deformations structures in fault damage zones along which fluids were preferentially channelled.     

Evolution of the Axial Geometry of the Southwest Indian Ocean Ridge between the Melville Fracture Zone and the Indian Ocean Triple Junction – Implications for Segmentation on Very Slow-Spreading Ridges

Geophysical data from 900 km of the Southwest Indian Ridge are used todescribe the pattern of evolution of the plate boundary between 61° Eand 70° E over the past 20 million years. The SWIR is anobliquely-opening, ultra slow-spreading axis, and east of61° E comprises a series of ridge sections, each about 100–120 kmin length. The orientation of these sections varies fromsub-orthogonal to oblique to the approximately N–S spreadingdirection. In general, the suborthogonal sections are shallower, commonlysubdivided into an array of discrete axial segments, and carry recognisablecentral magnetic anomalies. The majority of the oblique sections are single,continuous rifts without continuous axial magnetic signatures.Morphotectonics of the Southwest Indian Ridge crust have not previously beenwell constrained off-axis, and we here present sidescan sonar andswath bathymetric data up to 100 km from the ridge to demonstrate the complexities of its spatial and temporal evolution.A model is proposed that the segmentation style correlates with analong-axis variation between: (a) relatively thick crustal sections which overlie mantle sections with higher magmatic supply created in orthogonally-spreading segments and (b) those oblique sections associated with cooler, magmatically-starved mantle and thinner crust. These latter sections are formed at broad offset zones in theplate boundary, more precisely defined on faster-spreading ridges asnontransform discontinuities. The nonsystematic pattern of crustalconstruction, extensional basin formation and the absence of extension-parallel traces of discontinuities off-axis suggest that the oblique spreading sections are not fixed in space or time.     

Propagation of the Southwest Indian Ridge at the Rodrigues Triple Junction

The analysis of multibeam bathymetric data of the Southwest Indian Ridge(SWIR) domain between the triple junction traces from 68° E to theRodrigues Triple Junction (RTJ; 70° E) reveals the evolution of thisridge since magnetic anomaly 4 (8 Ma). Image processing has been used toshow that the horizontal component of strain due to a network of normal stepfaults increases dramatically between 69°30 E and the RTJ. Thisarea close to the RTJ is characterized by a deep graben at the foot of thetriple junction trace on the African plate and by a narrow fault-boundedridge that joins an offset of the trace on the Antarctic plate. In thatarea, spreading is primarily amagmatic and dominated by tectonic extensionprocesses. To the west of 69°30 E, some lobate bathymetricfeatures atop of a large topographic high suggest volcanic constructions.Between 68°10 E and 69°25 E the southern flank of theSWIR domain is wider than the northern one and is characterized by a series of 7 en echelon bathymetric highs similar in size,shape and orientation to the one centred at 69°30E near the present-day triple junction. Their en echelon organization along the triple junction trace on the Antarctic plate and the typical lack of conjugated parts on the northern flank show that these bathymetric highs have been shifted to the south by successive northward relocalisations of the SWIR rifting zone. This evolution results in the asymmetric spreading of the SWIR in the survey area. The off-axis bathymetric highs connect to the offsets of the triple junction trace on the Antarctic plate when the Southeast Indian Ridges lightly lengthenstoward the northwest and the triple junction is relocated to the north. We propose that the SWIR lengthens toward the northeast with two propagation modes: 1) a continuous and progressive propagation with distributed deformation in preexisting crust of the Central Indian Ridge, 2) a discontinuous propagation with focusing of the deformation in a rift zone when the triple junction migrates rapidly to the north. The modes of propagation of the SWIR are related to different localisation and distribution of strain which are in turn controlled by changes of the triple junction configurations due to propagation, recession or a symmetric spreading on the Central and Southeast Indian Ridges.     

Robust dual‐sensor data sea‐surface ghost attenuation for quality control purposes

Three‐dimensional receiver ghost attenuation (deghosting) of dual‐sensor towed‐streamer data is straightforward, in principle. In its simplest form, it requires applying a three‐dimensional frequency–wavenumber filter to the vertical component of the particle motion data to correct for the amplitude reduction on the vertical component of non‐normal incidence plane waves before combining with the pressure data. More elaborate techniques use three‐dimensional filters to both components before summation, for example, for ghost wavelet dephasing and mitigation of noise of different strengths on the individual components in optimum deghosting. The problem with all these techniques is, of course, that it is usually impossible to transform the data into the crossline wavenumber domain because of aliasing. Hence, usually, a two‐dimensional version of deghosting is applied to the data in the frequency–inline wavenumber domain. We investigate going down the “dimensionality ladder” one more step to a one‐dimensional weighted summation of the records of the collocated sensors to create an approximate deghosting procedure. We specifically consider amplitude‐balancing weights computed via a standard automatic gain control before summation, reminiscent of a diversity stack of the dual‐sensor recordings. This technique is independent of the actual streamer depth and insensitive to variations in the sea‐surface reflection coefficient. The automatic gain control weights serve two purposes: (i) to approximately correct for the geometric amplitude loss of the Z data and (ii) to mitigate noise strength variations on the two components. Here, Z denotes the vertical component of the velocity of particle motion scaled by the seismic impedance of the near‐sensor water volume. The weights are time‐varying and can also be made frequency‐band dependent, adapting better to frequency variations of the noise. The investigated process is a very robust, almost fully hands‐off, approximate three‐dimensional deghosting step for dual‐sensor data, requiring no spatial filtering and no explicit estimates of noise power. We argue that this technique performs well in terms of ghost attenuation (albeit, not exact ghost removal) and balancing the signal‐to‐noise ratio in the output data. For instances where full three‐dimensional receiver deghosting is the final product, the proposed technique is appropriate for efficient quality control of the data acquired and in aiding the parameterisation of the subsequent deghosting processing.     

Generation of 3D Spatially Variable Anisotropy for Groundwater Flow Simulations

Sedimentary units generally present anisotropy in their hydraulic properties, with higher hydraulic conductivity along bedding planes, rather than perpendicular to them. This common property leads to a modeling challenge if the sedimentary structure is folded. In this paper, we show that the gradient of the geological potential used by implicit geological modeling techniques can be used to compute full hydraulic conductivity tensors varying in space according to the geological orientation. For that purpose, the gradient of the potential, a vector normal to the bedding, is used to construct a rotation matrix that allows the estimation of the 3D hydraulic conductivity tensor in a single matrix operation. A synthetic 2D cross section example is used to illustrate the method and show that flow simulations performed in such a folded environment are highly influenced by this rotating anisotropy. When using the proposed method, the streamlines follow very closely the folded formation. This is not the case with an isotropic model.