Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
Although magmatism may occur during the earliest stages of continental rifting, its role in strain accommodation remains weakly constrained by largely 2-D studies. We analyze seismicity data from a 13 month, 39-station broadband seismic array to determine the role of magma intrusion on state-of-stress and strain localization, and their along-strike variations. Precise earthquake locations using cluster analyses and a new 3-D velocity model reveal lower crustal earthquakes beneath the central basins and along projections of steep border faults that degas CO2. Seismicity forms several disks interpreted as sills at 6-10 km below a monogenetic cone field. The sills overlie a lower crustal magma chamber that may feed eruptions at Oldoinyo Lengai volcano. After determining a new ML scaling relation, we determine a b-value of 0.87±0.03. Focal mechanisms for 65 earthquakes, and 13 from a catalogue prior to our array reveal an along-axis stress rotation of ∼60° in the magmatically active zone. New and prior mechanisms show predominantly normal slip along steep nodal planes, with extension directions ∼N90°E north and south of an active volcanic chain consistent with geodetic data, and ∼N150°E in the volcanic chain. The stress rotation facilitates strain transfer from border fault systems, the locus of early-stage deformation, to the zone of magma intrusion in the central rift. Our seismic, structural, and geochemistry results indicate that frequent lower crustal earthquakes are promoted by elevated pore pressures from volatile degassing along border faults, and hydraulic fracture around the margins of magma bodies. Results indicate that earthquakes are largely driven by stress state around inflating magma bodies.
ISSN: 1525-2027

hal-01667198
https://hal.archives-ouvertes.fr/hal-01667198
DOI : 10.1002/2017GC007027

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
Rifting in a cratonic lithosphere is strongly controlled by several interacting processes including crust/mantle rheology, magmatism, inherited structure and stress regime. In order to better understand how these physical parameters interact, a 2 yr long seismological experiment has been carried out in the North Tanzanian Divergence (NTD), at the southern tip of the eastern magmatic branch of the East African rift, where the southward-propagating continental rift is at its earliest stage. We analyse teleseismic data from 38 broad-band stations ca. 25 km spaced and present here results from their receiver function (RF) analysis. The crustal thickness and Vp/Vs ratio are retrieved over a ca. 200 × 200 km2 area encompassing the South Kenya magmatic rift, the NTD and the Ngorongoro-Kilimanjaro transverse volcanic chain. Cratonic nature of the lithosphere is clearly evinced through thick (up to ca. 40 km) homogeneous crust beneath the rift shoulders. Where rifting is present, Moho rises up to 27 km depth and the crust is strongly layered with clear velocity contrasts in the RF signal. The Vp/Vs ratio reaches its highest values (ca. 1.9) beneath volcanic edifices location and thinner crust, advocating for melting within the crust. We also clearly identify two major low-velocity zones (LVZs) within the NTD, one in the lower crust and the second in the upper part of the mantle. The first one starts at 15–18 km depth and correlates well with recent tomographic models. This LVZ does not always coexist with high Vp/Vs ratio, pleading for a supplementary source of velocity decrease, such as temperature or composition. At a greater depth of ca. 60 km, a mid-lithospheric discontinuity roughly mimics the step-like and symmetrically outward-dipping geometry of the Moho but with a more slanting direction (NE–SW) compared to the NS rift. By comparison with synthetic RF, we estimate the associated velocity reduction to be 8–9 per cent. We relate this interface to melt ponding, possibly favouring here deformation process such as grain-boundary sliding (EAGBS) due to lithospheric strain. Its geometry might have been controlled by inherited lithospheric fabrics and heterogeneous upper mantle structure. We evidence that crustal and mantle magmatic processes represent first order mechanisms to ease and locate the deformation during the first stage of a cratonic lithospheric breakup.
ISSN: 0956-540X

hal-01622078
https://hal.archives-ouvertes.fr/hal-01622078
https://hal.archives-ouvertes.fr/hal-01622078/document
https://hal.archives-ouvertes.fr/hal-01622078/file/plasmanGJI2017.pdf
DOI : 10.1093/gji/ggx177

Éditeur(s) : HAL CCSD Nature Publishing Group
Résumé : International audience
Lateral variations along the Himalayan arc are suggested by an increasing number of studies and carry important information about the orogen’s segmentation. Here we compile the hitherto most complete land gravity dataset in the region which enables the currently highest resolution plausible analysis. To study lateral variations in collisional structure we compute arc-parallel gravity anomalies (APaGA) by subtracting the average arc-perpendicular profile from our dataset; we compute likewise for topography (APaTA). We find no direct correlation between APaGA, APaTA and background seismicity, as suggested in oceanic subduction context. In the Himalayas APaTA mainly reflect relief and erosional effects, whereas APaGA reflect the deep structure of the orogen with clear lateral boundaries. Four segments are outlined and have disparate flexural geometry: NE India, Bhutan, Nepal & India until Dehradun, and NW India. The segment boundaries in the India plate are related to inherited structures, and the boundaries of the Shillong block are highlighted by seismic activity. We find that large earthquakes of the past millennium do not propagate across the segment boundaries defined by APaGA, therefore these seem to set limits for potential rupture of megathrust earthquakes.
ISSN: 2045-2322

hal-01420404
https://hal.archives-ouvertes.fr/hal-01420404
https://hal.archives-ouvertes.fr/hal-01420404/document
https://hal.archives-ouvertes.fr/hal-01420404/file/hetenyi2016.pdf
DOI : 10.1038/srep33866

Éditeur(s) : HAL CCSD Seismological Society of America
Résumé : International audience
We present an improvement of the maximum intersection (MAXI) method in which absolute earthquake location is defined by the maximum intersection number of hyperbolic equal differential time (EDT) volumes (one EDT volume being described as all grid nodes satisfying the arrival‐time differences between two stations, ± a tolerance value known as TERR). This 3D technique is well adapted to a strongly heterogeneous environment, avoids the depth versus origin‐time trade‐off (even using P arrivals alone), and objectively filters possible erroneous arrival times. Improvements consist of introducing an iterative multiscale approach on the TERR parameter and conducting a hierarchy of grid discretization to refine the preliminary solution (multigrid algorithm). The outlier cleaning process is strengthened by conducting statistical analyses based on the number of EDT volume intersections per station, rather than on travel‐time residues. Processing two series of synthetic arrival times computed in the 3D seismological configuration of the southernmost Ryukyu subduction zone, offshore Taiwan, allows us to (1) confirm MAXI robustness, even when arrival times are highly disturbed or azimuthal coverage largely reduced, and (2) show that MAXI, by limiting trade‐offs between origin time and depth and between epicentral position and depth, is less affected by a wrong velocity model than a technique that minimizes travel‐time residues. To conclude, we apply the MAXI method to locate earthquakes recorded by the combined Central Weather Bureau of Taiwan (CWB) and the Japanese Meteorological Agency (JMA) networks (1991-2008) within a georealistic a priori 3D P‐velocity model. Hypocenter determinations are selected based on MAXI confidence factors, and earthquake clustering along well‐known tectonic features witnesses hypocenter determination accuracy.
ISSN: 0037-1106

hal-00725227
https://hal.archives-ouvertes.fr/hal-00725227
DOI : 10.1785/0120100311

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The southernmost part of the Ryukyu subduction, where the Philippine Sea Plate is subducting under the Eurasian Plate, is known to be a very seismically active region of transition from a north-dipping subduction along the Ryukyu subduction to an ~ SE-NW collision along the Taiwanese orogenic wedge. In this paper, we will focus on the Ryukyu forearc area close to Taiwan where the deformation is paroxysmal. In order to decipher the nature of the seismic deformation in this region, a three month passive experiment, combining 22 Ocean Bottom Seismometers and 51 onland stations, has been led. Starting from an a-priori heterogeneous model, we have obtained 801 well-located earthquake hypocenters, a precise P-wave tomography model and 14 focal mechanisms. The seismicity along the Ryukyu forearc is mainly located not only in the vicinity of the Interplate Seismogenic Zone (ISZ) but also within both the subducting PSP and the overriding plate. Seismicity within the upper-plate is essentially localized east of Nanao basin where NW-SE extension occurs, and northwest of the Hoping basin where strike-slip dominates. As revealed by both the P-wave velocity structure and the newly derived seismicity, we argue that a sub-vertical step offsetting the subducting PSP around 10 km may support the presence of a trench-parallel tear. The PSP also undergoes extension in its upper part that is probably caused by buckling and slab pull. The P-wave velocity structure reveals three other major features: (1) a continuity between the Central Range and the Ryukyu Arc with a shallower Moho (~ 30 km depth) between ~ 122.3°N and ~ 122.5°N along the Ryukyu Arc, (2) high P-wave velocities along the eastern side of the Central Range and, (3) two bodies with similar high crustal velocities (6.5-7.0 km/s) at 12-18 km depths, embedded within the Ryukyu arc basement, just north of Hoping Basin and north of the Nanao Basin.
ISSN: 0040-1951

hal-00795521
https://hal.archives-ouvertes.fr/hal-00795521
DOI : 10.1016/j.tecto.2012.04.011

Éditeur(s) : HAL CCSD Hindawi Publishing Corporation
Résumé : During the 1982-1984 bradyseismic crises in the Campi Flegrei area (Italy), the University of Wisconsin deployed a network of seismological stations to record local earthquakes. In order to analyse the potential of the recorded data in terms of tomographic imaging, a blind test was recently set up and carried out in the framework of a research project. A model representing a hypothetical 3D structure of the area containing the Campi Flegrei caldera was also set up, and a synthetic dataset of time arrivals was in turn computed. The synthetic dataset consists of several thousand P- and S-time arrivals, computed at about fourteen stations. The tomographic inversion was performed by four independent teams using different methods. The teams had no knowledge of either the input velocity model or the earthquake hypocenters used to create the synthetic dataset. The results obtained by the different groups were compared and analysed in light of the true model. This work provides a thorough analysis of the earthquake tomography potential of the dataset recording the seismic activity at Campi Flegrei in the 1982-1984 period. It shows that all the tested earthquake tomography methods provide reliable low-resolution images of the background velocity field of the Campi Flegrei area, but with some differences. However, none of them succeeds in detecting the hypothetical structure details (i.e. with a size smaller than about 1.5-2 km), such as a magmatic chamber 4 km deep and especially the smaller, isolated bodies, which represent possible magmatic chimneys and intrusions.
ISSN: 1687-885X

hal-00795584
https://hal.archives-ouvertes.fr/hal-00795584
DOI : 10.1155/2012/505286

Éditeur(s) : HAL CCSD
Résumé : International audience
The Corinth Rift is the most seismically active zone in Europe. Since the year 2000, the Corinth Rift Laboratory (CRL, http://crlab.eu) consisting in a multidisciplinary natural observatory, aims at understanding the mechanics of faulting and earthquake nucleation in the Rift. Recent studies have improved our view about fault geometry and mechanics within CRL but there is still a critical need for a better knowledge of the structure at depth.In this project, we aim to analyze the complete seismological database (13 years of recordings) of CRL by using recently developed methodologies of structural imaging in order to determine, with high resolution, the local 3D structure and the earthquake locations. We perform an iterative joint determination of 3D velocity model and earthquake coordinates. P and S velocity models are estimated through a first arrival traveltime tomography by using the adjoint state method proposed by Taillandier et al. (2009). It consists in the minimization of the cost function between observed and theoretical arrival times,which is achieved by estimating the gradient of the cost function. Earthquake locations are also estimated with a linearized inversion approach by computing the gradient of the misfit function. Iterations are repeated until the cost function no longer decreases. We present preliminary results consisting in: (1) the adjustment of a 1D velocity model that is used as initial model of the 3D tomography and (2) a first attempt of the joint determination of 3D velocity model and earthquake locations.We also perform a preliminary receiver function analysis of teleseismic data recorded by the broadband stations of the CRL network. The multiscale approach (Gesret et. al., 2010) analysis of teleseismic P to S converted waves allows to analyze the signal toward the higher frequencies and thus to characterize the structure at depth with a high resolution. Our multiscale approach should give us reliable depth of interfaces beneath seismometers in the CRL area. In this first attempt, we adjust the 1D velocity model that produces a synthetic RF as similar as possible to the stack of observed RF for a subset of data. The identified interfaces could then be compared with heterogeneities imaged by the tomography.
AGU Fall Meeting
San Fransico, United States

hal-01417184
https://hal-mines-paristech.archives-ouvertes.fr/hal-01417184

Éditeur(s) : HAL CCSD Wiley
Résumé : International audience
In this paper we investigate finite-frequency effects in crustal tomography. We developed an inversion procedure based on an exact numerical computation of the sensitivity kernels. In this approach we compute the 3D travel-time sensitivity kernels by using (1) graph theory and an additional bending to estimate accurately both rays and travel-times between source/receiver and diffraction points and (2) paraxial ray theory to estimate the amplitude along theses rays. We invert both the velocity and the hypocentre parameters, using these so-called banana-doughnut kernels and the LSQR iterative solver. We compare the ray-theoretical and the finite-frequency tomography to image the intermediate structures beneath the Gulf of Corinth (Greece), which has long been recognized as the most active continental rifting zone in the Mediterranean region. Our dataset consists of 451 local events with 9233 P- firstarrival times recorded in the western part of the Gulf (Aigion area) in the framework of the 3F-Corinth European project. Previous tomographic images showed a complex velocity crustal model and a low-dip surface that may accommodate the deformation. Accurate velocity models will help to better constrain the rifting process, which is still a subject of debate. The main results of this study show that finite-frequency tomography improves crustal tomographic images by providing better resolved images of the 3D complicated velocity structure. Because the kernels spread the information over a volume, finite-frequency tomography results in a sharpening of layer boundaries as we observed for the shallower part of the crust (down to 5 km depth) beneath the Gulf of Corinth.
ISSN: 0016-8025

insu-00354702
https://hal-insu.archives-ouvertes.fr/insu-00354702
DOI : 10.1111/j.1365-2478.2007.00683.x

Éditeur(s) : HAL CCSD Elsevier Masson
Résumé : On the basis of an existing multiparameter borehole fluid sensor(p, T Cw, pH, Eh), a new downhole tool designed for 200 bars and 75 degrees C was developed to measure the spontaneous electrical potential in situ with great precision (better than a mu V). To this end, the new sensor is based on the use of unpolarizable Pb/PbCl2 electrodes either at the surface or downhole. In situ testing has demonstrated a capacity to identify several subsurface sources of natural electrical potential such as diffusion ones (membrane potential in the presence of clays, Fickian processes due to pore fluid salinity gradients), or else the electrokinetic mechanisms with gradients in pore fluid pressure. As a consequence, Multi-Sensors Electrical Tool (MuSET) might be used as an "electrical flow meter" sensitive to both horizontal and vertical fluid flow in a vertical borehole.
ISSN: 1631-0713

hal-00445251
https://hal.archives-ouvertes.fr/hal-00445251
DOI : 10.1016/j.crte.2009.07.009

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
An exploration 3D seismic data set from the Gjallar Ridge off mid-Norway images a giant fluid seep structure 3 × 5 km wide, which connects to late Palaeocene magmatic sills at depth. Two of the pipes that have developed as hydrothermal vents reach all the way to the modern seafloor implying that they either were active much longer than the original hydrothermal activity or have been reactivated. We combine detailed seismic analysis of the northern pipe and sandbox modeling to constrain pipe initiation and propagation. Although both the seismic data and the sandbox models suggest that fluids at depth are focused through a vertical conduit, sandbox models show that fluids ascend and reach a critical depth migration where focused migration abruptly transforms into distributed fluid flow through unconsolidated sediments. This indicates that at this level the sediments are intensely deformed during pipe propagation, creating a V-shaped structure, i.e. an inverted cone at depth and a positive relief anomaly, 5 to 10 m high, at the seafloor, which is clearly identified on 3D seismic data. Comparison of the geometries observed in sandbox modeling with the seismically observed geometries of the Giant Gjallar Vent suggests that the Giant Gjallar Vent may be a proto-fluid seep at an early stage of its development, preceding the future collapse of the structure forming a seafloor depression. Our results imply that the Gjallar Giant Vent can be used as a window into the geological processes active in the deep parts of the Vøring Basin.
ISSN: 0025-3227

hal-00772465
https://hal.archives-ouvertes.fr/hal-00772465
DOI : 10.1016/j.margeo.2012.08.010

Éditeur(s) : HAL CCSD
Résumé : International audience
The Corinth Rift is the most seismically active zone in Europe. The area is characterized by very localized NSextension at a rate of1.5cm/year, the occurrence of frequent and intensive microseismic crises and occasionalmoderate to large earthquakes like in 1995 (Mw=6.1).Since the year 2000, the Corinth Rift Laboratory (CRL, http://crlab.eu) consisting in a multidisciplinary naturalobservatory, aims at understanding the mechanics of faulting and earthquake nucleation in the Rift. Recent studieshave improved our view about fault geometry and mechanics within CRL, but there is still a critical need for a betterknowledge of the structure at depth both for the accuracy of earthquake locations and for mechanical interpretationof the seismicity.In this project, we aim to analyze the complete seismological database (13 years of recordings) of CRL by usingrecently developed methodologies of structural imaging, in order to determine at the same time and with highresolution, the local 3D structure and the earthquake locations.We perform an iterative joint determination of 3D velocity model and earthquake coordinates. In a first step, Pand S velocity models are determined using first arrival time tomography method proposed by Taillandier et al.(2009). It consists in the minimization of the cost function between observed and theoretical arrival times which isachieved by the steepest descent method (e.g. Tarantola 1987). This latter requires computing the gradient of thecost function by using the adjoint state method (Chavent 1974). In a second step, earthquakes are located in thenew velocity model with a non-linear inversion method based on a Bayesian formulation (Gesret et al. 2015). Step1 and 2 are repeated until the cost function no longer decreases.We present preliminary results consisting in: (1) the adjustement of a 1D velocity model that is used as initialmodel of the 3D tomography and (2) a first attempt of the joint determination of 3D velocity model and earthquakelocation.In addition to the tomographic imaging, we perform a preliminary receiver function analysis of teleseismic datarecorded by the broadband stations of the CRL network. The RF analysis should provide the interface depthsbeneath seismometers and increase the imaging resolution of the upper crustal structures provided by the 3Dtomography. In this first attempt, we adjust the 1D velocity model that produces a synthetic RF as similar aspossible to the observed RF for a subset of data. We compare the identified interfaces with structures imaged bythe tomography.
EGU general assembly 2016
Vienne, Austria

hal-01306224
https://hal-mines-paristech.archives-ouvertes.fr/hal-01306224

Éditeur(s) : HAL CCSD Wiley
Résumé : Clay shale landscapes are sensitive to landslide, erosion, and re-sedimentation processes. In this context, tectonized clay shales have peculiar hydrological and mechanical behaviours. In order to improve our understanding of the processes involved in such systems, an experimental site has been settled in the Draix Observatory (ORE DRAIX) on a black marl hill slope of the Southern Alps (France) in the framework of the ECOU-PREF project. The site is a natural laboratory to study the evolution of clay shale landscapes at several scales (regional to µm). In this study, we aim to characterize the internal structure and fissure patterns of a stable interfluve, in order to locate preferential water flows within the system, and anticipate its evolution. This work is based on a series of mutli-disciplinary approaches including petrophysical analyses, well pulse injection tests and downhole geophysical measurements. Borehole geophysical results were combined with laboratory measurements on core plugs (permeability, density, porosity, and acoustic velocity). Optical and acoustical images have been used as an original tool for internal discontinuities characterisation and potential active flow path detection. They allowed the identification of three main lithological units interpreted as several stages of shale alteration. In addition, several main structural plane discontinuities (open fractures, schistosity planes, and fractures infilled with calcite or clays) were identified from the images and their potential transmissivity discussed with regard to piezometric and tracer measurements. In all, the integration of measurements leads to propose a simple scenario of fluid circulation and chemical alteration of the interfluve.
ISSN: 0885-6087

hal-00757581
https://hal.archives-ouvertes.fr/hal-00757581
DOI : 10.1002/hyp.7997

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
The Eastern Rift System (ERS) of northern Tanzania and southern Kenya, where a cratonic lithosphere is in the early stages of rifting, offers an ideal venue for investigating the roles of magma and other fluids in such an environment. To illuminate these roles, we jointly invert arrival times of locally recorded P and S body waves, phase delays of ambient noise generated Rayleigh waves and Bouguer anomalies from gravity observations to generate a 3-D image of P and S wave speeds in the upper 25 km of the crust. While joint inversion of gravity and arrival times requires a relationship between density and wave speeds, the improvement in resolution obtained by the combination of these disparate data sets serves to further constrain models, and reduce uncertainties. The most significant features in the 3-D model are (1) P and S wave speeds that are 10–15 per cent lower beneath the rift zone than in the surrounding regions, (2) a relatively high wave speed tabular feature located along the western edge of the Natron and Manyara rifts, and (3) low (∼1.71) values of Vp/Vs throughout the upper crust, with the lowest ratios along the boundaries of the rift zones. The low P and S wave speeds at mid-crustal levels beneath the rift valley are an expected consequence of active volcanism, and the tabular, high-wave speed feature is interpreted to be an uplifted footwall at the western edge of the rift. Given the high levels of CO2 outgassing observed at the surface along border fault zones, and the sensitivity of Vp/Vs to pore-fluid compressibility, we infer that the low Vp/Vs values in and around the rift zone are caused by the volcanic plumbing in the upper crust being suffused by a gaseous CO2 froth on top of a deeper, crystalline mush. The repository for molten rock is likely located in the lower crust and upper mantle, where the Vp/Vs ratios are significantly higher.
ISSN: 0956-540X

hal-01622055
https://hal.archives-ouvertes.fr/hal-01622055
https://hal.archives-ouvertes.fr/hal-01622055/document
https://hal.archives-ouvertes.fr/hal-01622055/file/roeckerGJI2017.pdf
DOI : 10.1093/gji/ggx220