Éditeur(s) : HAL CCSD
Résumé : International audience
Bulletin de la Société Géologique de France special issue "Caribbean Geosciences"

hal-01313406
https://hal.archives-ouvertes.fr/hal-01313406

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
The association of experimental data showing that the plastic deformation of olivine, the main constituent of the upper mantle, is highly anisotropic and the ubiquitous seismic anisotropy in the upper mantle, which indicates that olivine crystals show coherent orientations over scales of tens to hundreds of kilometers, implies that the long-term deformation in the upper mantle is anisotropic. We propose a multiscale approach, based on a combination of finite element and homogenization techniques, to model the deformation of a lithospheric plate while fully considering the mechanical anisotropy stemming from a strain-induced orientation of olivine crystals in the mantle. This multiscale model explicitly takes into account the evolution of crystal preferred orientations (CPO) of olivine and of the mechanical anisotropy during the deformation. We performed a series of numerical experiments simulating the uniaxial extension of a homogeneous (100% olivine) but anisotropic plate to test the role of the olivine CPO on the plate mechanical behavior and the link between CPO and mechanical anisotropy evolution. Even for this simple solicitation, different orientations and intensity of the initial olivine CPO result in variable plate strengths and deformation regimes. A plate with an initial CPO where the olivine [100] and [010] axes are concentrated at 45 degrees to the extension direction has high resolved shear stresses on the easy (010)[100] and (001)[100] slip systems of olivine. This results in low strength and in deformation by transtension. Plates with an initial CPO where the maximum of [100] axes is parallel or normal to the extension direction show a high initial strength. Isotropic plates have an intermediate behavior. The progressive rotation of olivine [100] axes toward the imposed stretching direction results in hardening in all models, except in those characterized by an initial concentration of olivine [100] axes normal to the imposed extension, in which softening is followed by hardening.
ISSN: 1525-2027

hal-00420063
https://hal.archives-ouvertes.fr/hal-00420063
DOI : 10.1029/2009GC002423

Éditeur(s) : HAL CCSD
2012 AGU Fall Meeting

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

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : Results are reported for an acidic water flow-through experiment conducted in a fractured argillaceous limestone sample (73% carbonates). The change in fracture geometry and related parameters is reported for six data sets obtained from synchrotron X-ray microtomography experiments. High-resolution three-dimensional images of the sample allowed quantification of the changes in fracture morphology at a spatial resolution of 6 ;μm. Mineral mass loss and permeability changes in the sample were also determined. Several physico-chemical phenomena were identified during the experiment. Initial smooth fracture surfaces evolved rapidly toward rough surfaces with uneven clay coverage due to the preferential dissolution of carbonate minerals compared to clay minerals whose dissolution rate is about 106 slower. A microporous clay coating progressively developed at the fluid-rock interface during heterogeneous dissolution of the fracture, while the global dissolution rate of the fracture walls exponentially decreased. The increase in surface roughness and the presumed reorganization of clays caused a progressive reduction in permeability. During the last flow-through stage, a large decrease in sample permeability was attributed to the large removal of clay particles; this process was responsible for a dramatic collapse of the fracture walls near the sample inlet and led to the development of preferential flow pathways. The development of the clay coating also acted as a barrier to flow and mass transfer between calcite grains and bulk solution and affected transport processes within the fracture
ISSN: 0043-1397

hal-00405333
https://hal.archives-ouvertes.fr/hal-00405333
DOI : 10.1029/2006WR005379

É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 Springer Verlag
Résumé : International audience
This paper presents a synthesis of a multidisciplinary study carried out along the Nîmes Roman aqueduct, located in the southeast part of France. The study was motivated by archaeologists attempting to explain the partial destructions of only one aerial bridge of the aqueduct (Pont de la Lône). Given its close proximity to the Nîmes seismically active fault, a possible seismic origin for the destructions was invoked. Seismologists and structural engineers thus carried out a variety of field and numerical investigations to test the seismic hypothesis. Supporting field evidence was found first along the aerial bridge section of the aqueduct: broken stalactites, arch warping, cracks and destruction of the bridge just above a breast wall shortly after its construction. Secondly, the underground part of the canal was analysed: irregularity in the thickness of calcite deposits of the canal walls, presence of numerous cracks, horizontal shift of the otherwise linear structure of the canal and presence of calcite twins in the deposits, found only where the canal crosses the Nîmes fault system. Numerical modelling and experimental results show that (1) only the Pont de la Lône would have suffered serious damage under seismic solicitation (assuming an M6 earthquake at 10 km distance). The second aerial bridge, the three levels arches Pont du Gard, an historical monument still standing today, would have suffered less damage, due to its very different fundamental frequency of around 0.4 Hz, far from the amplified seismic near-field spectrum; (2) the numerical models also show that a fall of the canal wall would require a higher magnitude event (M > 6); (3) the presence of calcite twinning requires a differential static stress of 4 MPa or greater, which can only be achieved very close to the fault rupture of a M > 6 earthquake; (4) finally, local surface rupturing of such a fault would also corroborate the hypothesis that the observed offset of the canal may be partly seismically induced. Although other possible origins for each individual evidence may not be excluded, the observed spatiotemporal concentration of architectural/geological anomalies together with the numerical results allow us to support a possible co-seismic origin for these disorders, indirectly attesting to the potential seismic activity (M > 6) of the nearby Nîmes fault. Furthermore, following the conclusions of archaeological studies, the disorders occurred between 250 and 350 year AD, thus providing timing for this possible seismic event, an essential parameter in seismic hazard analysis. These results validate the archaeo-seismological approach as a tool that may help improve the knowledge of major infrequent earthquakes in areas of moderate seismic activity.
ISSN: 0921-030X

hal-00379064
https://hal.archives-ouvertes.fr/hal-00379064
DOI : 10.1007/s11069-008-9276-9

Éditeur(s) : HAL CCSD Nature Publishing Group
Résumé : International audience
Reactivation of structures inherited from previous collisional or rifting events, especially lithospheric-scale faults, is a major feature of plate tectonics. Its expression ranges from continental break-up along ancient collisional belts(1,2) to linear arrays of intraplate magmatism and seismicity(3,4). Here we use multiscale numerical models to show that this reactivation can result from an anisotropic mechanical behaviour of the lithospheric mantle due to an inherited preferred orientation of olivine crystals. We explicitly consider an evolving anisotropic viscosity controlled by the orientation of olivine crystals in the mantle. We find that strain is localized in domains where shear stresses on the inherited mantle fabric are high, and that this leads to shearing parallel to the inherited fabric. During rifting, structural reactivation induced by anisotropy results in oblique extension, followed by either normal extension or failure. Our results suggest that anisotropic viscosity in the lithospheric mantle controls the location and orientation of intraplate deformation zones that may evolve into new plate boundaries, and causes long-lived lithospheric-scale wrench faults, contributing to the toroidal component of plate motions on Earth.
ISSN: 1752-0894

hal-00413033
https://hal.archives-ouvertes.fr/hal-00413033
DOI : 10.1038/ngeo528

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
This study explores the dynamics of porosity and reactive surface area changes during porous limestone dissolution by CO2-rich water. The Sr and Ca concentrations in both the rock and the outlet solution are used to evaluate the reactive surface area changes of the two rock-forming calcites, i.e. micrite grains and sparite crystals, which have different trace element signatures. The geometric surface area measured with X-ray microtomography decreases slightly whereas the reactive surface area increases continuously with increasing porosity from 20.3 to 30.2%. Surprisingly, changes in reactive surface areas are very different between the two calcites. The reactive surface area changes in the micrite are parabolic while the reactive surface area of sparite increases greatly. The numerical model HYTEC is used to model the change in reactive surface areas during the experiment. Different geometrical models are tested. The model based on spherical-grain dissolution and spherical-pore growth fails to reproduce the experimental results, while the sugar-lump model provides reasonable agreement with the experiment.
ISSN: 0009-2541

hal-00445274
https://hal.archives-ouvertes.fr/hal-00445274
DOI : 10.1016/j.chemgeo.2009.01.032

Éditeur(s) : HAL CCSD
Résumé : International audience
Joint Meeting of the Statistical Society of Canada and the Société Française de Statistique
Ottawa, Canada

hal-00652079
https://hal.inria.fr/hal-00652079

Éditeur(s) : HAL CCSD
Résumé : National audience
Journées MAS de la SMAI
Rennes, France

hal-00652081
https://hal.inria.fr/hal-00652081

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
ISSN: 0040-1951

hal-00843380
https://hal-mines-paristech.archives-ouvertes.fr/hal-00843380
DOI : 10.1016/j.tecto.2013.05.043

Éditeur(s) : HAL CCSD IJCA
Résumé : International audience
International Journal of Computer Applications

hal-01313407
https://hal.archives-ouvertes.fr/hal-01313407

É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
2012 AGU Fall Meeting

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

Éditeur(s) : HAL CCSD
Résumé : International audience
RST 2008
Nancy, France

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