Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Talc is a hydrous magnesium rich layered silicate that is widely disseminated in the Earth from the seafloor to over 100 km depth, in ultra-high pressure metamorphism of oceanic crust. In this paper we determine the single crystal elastic constants at pressures from 0 to 12 GPa of talc triclinic (C 1¯) and monoclinic (C2/c) polytypes using ab initio methods. We find that talc has an extraordinarily high elastic anisotropy at zero pressure that reduces with increasing pressure. The exceptional anisotropy is complemented by a negative Poisson's ratio for many directions in crystal space. Calculations show that talc is not only one of very few common minerals to exhibit auxetic behaviour, but the magnitude of this effect may be the largest reported so far for a mineral. The compression (Vp) and shear (Vs) wave velocity anisotropy is 80% and 85% for the triclinic polytype. At pressures where talc is known be stable in the Earth (up to 5 GPa) the Vp and Vs anisotropy is reduced to about 40% for both velocities, which is still a very high value. Vp is slow parallel to the c-axis and fast perpendicular to it. This remains unchanged with increasing pressure and is observed in both polytypes. The shear wave splitting (difference between fast and slow S-wave velocities) at low pressure has high values in the plane normal to the c-axis, with a maximum near the alow asterisk-axis in the triclinic and the b-axis in the monoclinic polytype. The c-axis is the direction of minimum splitting. The pattern of shear wave splitting does not change significantly with pressure. The volume fraction of talc varies between 11 and 41% for hydrated mantle rocks, but the lack of data on the crystallographic preferred orientation (CPO) precludes a detailed analysis of the impact of talc on seismic anisotropy in subduction zones. However, it is highly likely that CPO can easily develop in zones of deformation due to the platy habit of talc crystals. For random aggregates of talc, the isotropic Vp, Vs and Vp/Vs ratio have significantly lower values than those of antigorite and may explain low-velocity regions in the mantle wedge. Vp/Vs ratios are more complex in anisotropic media because there are fast and slow S-waves, resulting in Vp/Vs1 and Vp/Vs2 ratios for every propagation direction, making interpretation difficult in deformed polycrystalline talc with a CPO. Talc on the subduction plate boundary can strongly influence guided wave velocity as CPO would develop in this region of intense shearing. The very low coefficient of friction (< 0.1) of talc above 100 °C could also explain silent earthquakes at shallow depths (ca 30 km) along the subduction plate boundaries, frequently responsible for tsunami.
ISSN: 0012-821X

hal-00411481
https://hal.archives-ouvertes.fr/hal-00411481
DOI : 10.1016/j.epsl.2008.07.047

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Laboratory models have been conducted to improve our understanding of the role that the resistance of the slab to bending and its coupling to the ambient mantle play in subduction dynamics over geological time scales. Our models are set up with a viscous plate of silicone (lithosphere) subducting under negative buoyancy in a viscous layer of glucose syrup (mantle). For our study, the lithosphere/upper mantle viscosity contrast has been systematically varied, from similar to 10 to similar to 10(5) in order to explore the parameter space between weak and strong slab dynamics. We found that subduction is characterized by a retreating mode for viscosity ratios > 10(4), by the coexistence of a retreating mode and an advancing mode for viscosity ratios between similar to 10(4) and similar to 10(2), and quasi-stationary, Rayleigh-Taylor like behaviour for ratios < 10(2). By combining our experimental results and kinematic data from current subduction zones in four reference frames which differ in the amount of net rotation, we infer that a lithosphere/upper mantle viscosity contrast of 150-500 is necessary to obtain realistic trench/subducting plate velocity ratios as well as the variability of subduction styles observed in nature.
ISSN: 0012-821X

hal-00412008
https://hal.archives-ouvertes.fr/hal-00412008
DOI : 10.1016/j.epsl.2008.04.006

Éditeur(s) : HAL CCSD Springer Verlag (Germany) Springer Verlag
Résumé : Direct numerical simulations of flow through two millimeter-scale rock samples of limestone and sandstone are performed using three diverse fluid dynamic simulators. The resulting steady-state velocity fields are compared in terms of the associated empirical probability density functions (PDFs) and key statistics of the velocity fields. The pore space geometry of each sample is imaged at 5.06−μm voxel size resolution using X-ray microtomography. The samples offer contrasting characteristics in terms of total connected porosity (about 0.31 for the limestone and 0.07 for the sandstone) and are typical of several applications in hydrogeology and petroleum engineering. The three-dimensional fluid velocity fields within the explicit pore spaces are simulated using ANSYS® FLUENT® ANSYS Inc. (2009), EULAG Prusa et al. (Comput. Fluids 37, 1193–1207 2008), and SSTOKES Sarkar et al. (2002). These computational approaches are highly disperse in terms of algorithmic complexity, differ in terms of their governing equations, the adopted numerical methodologies, the enforcement of internal no-slip boundary conditions at the fluid-solid interface, and the computational mesh structure. As metrics of comparison to probe in a statistical sense the internal similarities/differences across sample populations of velocities obtained through the computational systems, we consider (i) integral quantities, such as the Darcy flux and (ii) main statistical moments of local velocity distributions including local correlations between velocity fields. Comparison of simulation results indicates that mutually consistent estimates of the state of flow are obtained in the analyzed samples of natural pore spaces despite the considerable differences associated with the three computational approaches. We note that in the higher porosity limestone sample, the structures of the velocity fields obtained using ANSYS FLUENT and EULAG are more alike than either compared against the results obtained using SSTOKES. In the low-porosity sample, the structures of the velocity fields obtained by EULAG and SSTOKES are more similar than either is to the fields obtained using ANSYS FLUENT. With respect to macroscopic quantities, ANSYS FLUENT and SSTOKES provide similar results in terms of the average vertical velocity for both of the complex microscale geometries considered, while EULAG tends to render the largest velocity values. The influence of the pore space structure on fluid velocity field characteristics is also discussed.
ISSN: 1420-0597

hal-01174152
https://hal.archives-ouvertes.fr/hal-01174152
DOI : 10.1007/s10596-015-9486-7

Éditeur(s) : HAL CCSD
Résumé : International audience
We present a new 3-D GPS velocity solution for 182 sites for the region encompassing the Western Alps, Pyrenees, and southern France. The velocity field is based on a Precise Point Positioning (PPP) solution, to which we apply a common-mode filter, defined by the 26 longest time series, in order to correct for network-wide biases (reference frame, unmodeled large-scale processes, etc.). We show that processing parameters, such as troposphere delay modeling, can lead to systematic velocity variations of 0.1–0.5 mm yr−1 affecting both accuracy and precision, especially for short (< 5 years) time series. A velocity convergence analysis shows that minimum time-series lengths of  ∼  3 and  ∼  5.5 years are required to reach a velocity stability of 0.5 mm yr−1 in the horizontal and vertical components, respectively. On average, horizontal residual velocities show a stability of  ∼  0.2 mm yr−1 in the Western Alps, Pyrenees, and southern France. The only significant horizontal strain rate signal is in the western Pyrenees with up to 4  ×  10−9 yr−1 NNE–SSW extension, whereas no significant strain rates are detected in the Western Alps (< 1  ×  10−9 yr−1). In contrast, we identify significant uplift rates up to 2 mm yr−1 in the Western Alps but not in the Pyrenees (0.1 ± 0.2 mm yr−1). A correlation between site elevations and fast uplift rates in the northern part of the Western Alps, in the region of the Würmian ice cap, suggests that part of this uplift is induced by postglacial rebound. The very slow uplift rates in the southern Western Alps and in the Pyrenees could be accounted for by erosion-induced rebound.
ISSN: 1869-9529

hal-01395116
https://hal.archives-ouvertes.fr/hal-01395116
https://hal.archives-ouvertes.fr/hal-01395116/document
https://hal.archives-ouvertes.fr/hal-01395116/file/nguyen2016.pdf
DOI : 10.5194/se-7-1349-2016

Éditeur(s) : HAL CCSD Springer Verlag (Germany)
Résumé : International audience
Surface displacements due to temporal changes in environmental mass redistributions are observable in the coordinate time series of many Global Navigation Satellite System (GNSS) sites. In this study, we investigated the effect of loading on estimates of tectonic velocity computed from campaign-style GNSS observations. The study region is in the Pyrenees mountain range between France and Spain (ResPyr campaigns). In this area, seismic activity is continuous and moderate and the expected amplitude of the horizontal tectonic velocity is less than 0.5 mm/yr. In order to determine the velocity, 4 sparse GNSS campaigns were carried out from 1995 to 2010. Considering this small rate of deformation, loading phenomena can contribute a non-negligible artifact to the velocity computation that could affect our geodynamical interpretation. In this investigation, we specifically considered the atmospheric, hydrological, and non-tidal ocean loading phenomena. The computed loading deformations for this region show the horizontal displacements are dominated by the non-tidal ocean loading (maximum 4 mm for the North and 3.1 mm for the East components); the main vertical contributions are due to the atmospheric and continental water storage loading (maximum 14.3 for the atmosphere and 8.1 mm for the hydrology, respectively). We have found that the dominant loading effect on the horizontal velocity is the non-tidal ocean loading (mean of 0.11 mm/yr), whereas the vertical component is dominated by the hydrological loading (mean of 0.21 mm/yr). Since the study area is in a mountainous region, we also analyzed the difference between the atmospheric and the topography dependent atmospheric loading models at our GNSS campaign sites. We did not find any significant difference between the two atmospheric loading models in terms of horizontal velocity. Finally, we performed simulations to identify the optimum timing and frequency of future GNSS campaigns in this area that would minimize the loading effects on tectonic velocity estimates.
ISSN: 0039-3169

hal-00963629
https://hal.archives-ouvertes.fr/hal-00963629
DOI : 10.1007/s11200-012-0458-2

Éditeur(s) : American Geophysical Union
Résumé : As part of the French CLIPPER project, an eddy permitting model of the Atlantic circulation has been run for 22 years. The domain has open boundaries at Drake passage and at 30 degreesE, from Africa to Antarctica. The simulated mean circulation, as well as the eddy activity, is satisfactory for a 1/3 degrees model resolution, and the meridional heat transport at 30 degreesS is within the range estimated from observations. We use the "mixed" open boundary algorithm of Barnier et al. [1998], which has both a radiation condition and a relaxation to climatology. The climatological boundary forcing strongly constrains the solution in the whole domain. The model heat balance adjusts through the surface (heat flux retroaction term) more than the open boundaries. The radiation phase velocities calculated within the algorithm are analyzed. This shows, quite surprisingly, that both the eastern and western boundaries have a similar behavior, regardless of the preferred directions for advection (mainly eastward) and wave propagation (mainly westward). Our results confirm that open boundary algorithms behave differently according to the dynamics of the region considered. The passive boundary condition that Penduff et al. [2000] applied successfully in the north eastern Atlantic does not work in the present South Atlantic model. We emphasize the need for a careful prescription of the climatology at the open boundary, for which a new approach based on synoptic sections is implemented.
Journal of Geophysical Union - Research C - Oceans (0148-0227) (American Geophysical Union), 2001-10 , Vol. 106 , N. C10 , P. 22115-22129

Droits : 2001 American Geophysical Union
http://archimer.ifremer.fr/doc/2001/publication-454.pdf
DOI:10.1029/2000JC000376
http://archimer.ifremer.fr/doc/00000/454/

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
1] Recent statistical analysis by Lallemand et al. (2008) of subduction zone parameters revealed that the back-arc deformation mode depends on the combination between the subducting (nu(sub)) and upper (nu(up)) plate velocities. No significant strain is recorded in the arc area if plate kinematics verifies nu(up) = 0.5 vsub - 2.3 (cm/a) in the HS3 reference frame. Arc spreading ( shortening) occurs if nu(up) is greater ( lower) than the preceding relationship. We test this statistical law with numerical models of subduction, by applying constant plate velocities far away from the subduction zone. The subducting lithosphere is free to deform at all depths. We quantify the force applied on the two converging plates to sustain constant surface velocities. The simulated rheology combined viscous (non-Newtonian) and brittle behaviors, and depends on water content. The influence of subduction rate vs is first studied for a fixed upper plate. After 950 km of convergence ( steady state slab pull), the transition from extensional to compressive stresses in the upper plate occurs for vs similar to 1.4 cm/a. The effect of upper plate velocity is then tested at constant subduction rate. Upper plate retreat ( advance) with respect to the trench increases extension ( compression) in the arc lithosphere and increases ( decreases) the subducting plate dip. Our modeling confirms the statistical kinematic relationship between vsub and nu(up) that describes the transition from extensional to compressive stresses in the arc lithosphere, even if the modeled law is shifted toward higher rates of upper plate retreat, using our set of physical parameters ( e. g., 100 km thick subducting oceanic plate) and short- term simulations. Our results make valid the choice of the HS3 reference frame for assessing plate velocity influence on arc tectonic regime. The subduction model suggests that friction along the interplate contact and the mantle Stokes reaction could be the two main forces competing against slab pull for upper mantle subductions. Besides, our simulations show that the arc deformation mode is strongly time dependent.
ISSN: 1525-2027

hal-00412328
https://hal.archives-ouvertes.fr/hal-00412328
https://hal.archives-ouvertes.fr/hal-00412328/document
https://hal.archives-ouvertes.fr/hal-00412328/file/arcayG32008.pdf
DOI : 10.1029/2007GC001875

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
Two-dimensional simulations using a thermomechanical model based on a finite-difference method on a staggered grid and a marker in cell method are performed to study the plume-lithosphere interaction beneath moving plates. The plate and the convective mantle are modelled as a homogeneous peridotite with a Newtonian temperature- and pressure-dependent viscosity. A constant velocity, ranging from 5 to 12.5 cm yr−1, is imposed at the top of the plate. Plumes are generated by imposing a thermal anomaly of 150 to 350 K on a 50 km wide domain at the base of the model (700 km depth); the plate atop this thermal anomaly is 40 Myr old. We analyse (1) the kinematics of the plume as it impacts the moving plate, (2) the dynamics of time-dependent small-scale convection (SSC) instabilities developing in the low-viscosity layer formed by spreading of hot plume material at the base of the lithosphere and (3) the resulting thermal rejuvenation of the lithosphere. The spreading of the plume material at the base of the lithosphere, characterized by the ratio between the maximum down- and upstream horizontal (dimensionless) velocities in the plume-fed sublithospheric layer, Peup/Pedown depends on the ratio between the maximum plume upwelling velocity and the plate velocity, Peplume/Peplate. For fast plate velocities and sluggish plumes (low Peplume/Peplate), plate motion drags most plume material and downstream flow is dominant. As Peplume/Peplate increases, an increasing part of the plume material flows upstream. SSC systematically develops in the plume-fed sublithospheric layer, downstream from the plume. Onset time of SSC decreases with the Rayleigh number. For vigorous plumes, it does not depend on plate velocity. For more sluggish plumes, however, variations in the plume spreading behaviour at the base of the lithosphere result in a decrease in the onset time of SSCs with increasing plate velocity. In any case, SSC results in uplift of the isotherm 1573 K by up to 20 km relative to its initial equilibrium depth at the impact point.
ISSN: 0956-540X

hal-00856763
https://hal.archives-ouvertes.fr/hal-00856763
DOI : 10.1093/gji/ggt128

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
The three-dimensional S-wave velocity structure beneath the South Pacific Superswell is obtained from joint broadband seismic experiments on the ocean floor and islands. We collected only approximately 800 relative times of long-period teleseismic SH-waves by using a waveform cross-correlation from 76 events occurring from January 2003 to May 2005. We conducted relative time tomography to obtain a 3D structure to depths of 1600 km. In the resultant image, we find a characteristic distribution of low- velocity regions. The most prominent features are a large doughnut-shaped low-velocity region at 800 km depth, and an elongated large low-velocity region beneath the Society to Pitcairn hotspots at 1200 km depth. Our model suggests that a large low shear velocity province rooted in the D00 extends upwards and culminates near the top of the lower mantle beneath the central part of the South Pacific Superswell although its perfect continuity is not still confirmed.
ISSN: 0094-8276

hal-01389114
http://hal.univ-reunion.fr/hal-01389114
http://hal.univ-reunion.fr/hal-01389114/document
http://hal.univ-reunion.fr/hal-01389114/file/Tanaka_PLUME_lower_mantle_GRL2009.pdf
DOI : 10.1029/2009GL037568

Éditeur(s) : HAL CCSD American Geophysical Union (AGU) AGU and the Geochemical Society
Résumé : International audience
To study how an impacting plume modifies the mantle lithosphere, we analyzed the microstructures and crystal preferred orientations (CPO) of 29 peridotites and 37 pyroxenites that sample the mantle root of the Ontong Java Plateau (OJP) from 60 to 120 km depth. The peridotites show a strong compositional variability, but homogeneous coarse granular to tabular microstructures, except for those equilibrated at the shallowest and deepest depths, which are porphyroclastic. All peridotites have clear olivine CPO, with dominant fiber-[010] patterns. Low intragranular misorientations and straight grain boundaries in olivine suggest that, above 100 km depth, annealing often followed deformation. Calculated density and P wave velocities of the peridotites decrease weakly with depth. S wave velocities decrease faster, resulting in increasing Vp/Vs ratio with depth. Calculated densities and seismic velocity profiles are consistent with those estimated for normal mantle compositions under a cold oceanic geotherm. Enrichment in pyroxenites may further increase seismic velocities. The calculated seismic properties cannot therefore explain the low S waves velocities predicted by Rayleigh wave tomography and ScS data in the mantle beneath the OJP. Calculated P and S waves anisotropy is variable (2–12%). It is higher on average in the deeper section of the lithosphere. Because olivine has dominantly [010]-fiber CPO patterns, if foliations are horizontal, vertically propagating S waves and Rayleigh waves will sample very weak anisotropy in the OJP mantle lithosphere. Moreover, if the orientation of the lineation changes with depth, the anisotropy-induced contrast in seismic properties might produce an intralithospheric reflector marking the stratification of the OJP mantle root.
ISSN: 1525-2027

hal-01116428
https://hal.archives-ouvertes.fr/hal-01116428
DOI : 10.1002/2014GC005452

Éditeur(s) : North-Holland Publishing Company, Amsterdam
Résumé : The Barents Sea is underlain by a thick (up to 5 km) sedimentary basin. Seismic refraction has outlined four main velocity discontinuities which have been correlated with geological units on the basis of the geological history of the region. The basin is underlain by a crystalline basement, the nature of which cannot be determined on the basis of seismic velocity alone. Metamorphosed Paleozoic units (velocities around 5-4 km/sec) lie over this basement. Their thickness is not well established but appears to reach 1 or 2 km in some cases. A very distinct and thick (up to 2.5 km) layer (4.1 km/sec) is found almost everywhere and is thought to correspond to a major discontinuity at the end ot the Paleozoic. This discontinuity is followed by a variable velocity layer (2.4-4 km/sec, up to 1.2 km thick) which includes Mesozoic and Cenozoic sediments and is limited to the south of the Barents Sea. The absence of thick Tertiary deposits support the idea that the Barents shelf was emergent at that time. Recent low-velocity sediments are found close to the shelf edge. [NOT CONTROLLED OCR]
Earth and Planetary Science Letters (North-Holland Publishing Company, Amsterdam), 1974 , Vol. 24 , P. 33-47

Droits : North-Holland Publishing Company, Amsterdam
http://archimer.ifremer.fr/doc/1974/publication-5447.pdf
http://archimer.ifremer.fr/doc/00000/5447/

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : A differential effective medium (DEM) model is used to predict elastic properties for a set of porous and anisotropic aggregates, comprised of mixtures of calcite and muscovite. The DEM takes into consideration an anisotropic background medium with triclinic or higher symmetry, in which inclusions of idealized ellipsoidal shape are added incrementally. In general, the calculated elastic properties of a solid that contains inclusions representing "dry" pores/cracks are strongly dependent on the orientation and aspect ratio of the inclusions. Aspect ratios of inclusions in the synthetic aggregates, which consist of air-filled pores, are estimated from anisotropy of magnetic susceptibility (AMS) of samples whose pore space has been impregnated with a colloidal ferrofluid. The AMS derived pore shape geometry is used as an input value for inclusions in the DEM. Modeling results are compared with laboratory determined elastic properties, measured with ultrasonic waves. Calculated shear wave velocities agree in general well with laboratory measured S wave velocities, whereas calculated P wave velocities are typically 0.5-1.1 km/s higher than measured values. Differences between calculated and measured P wave velocities are attributed mainly to incomplete and biased ferrofluid saturation of pores. Spherical pores are preferably filled during imbibition, in comparison to thin cracks, which leads to overprediction of the calculated P wave velocities. The amount of ferrofluid that fills the pore space is dependent on the ratio of calcite to muscovite and the load used for compaction during sample synthesis. The permeability decreases with increasing muscovite content and increasing compaction load. Incomplete saturation of samples with high-muscovite content is confirmed by X-ray microtomography density contrast imaging of dry and ferrofluid saturated specimens.
ISSN: 0148-0227

hal-00617717
https://hal.archives-ouvertes.fr/hal-00617717
DOI : 10.1029/2010JB007750

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
A set of synthetic aggregates, containing mixtures of calcite and muscovite, were compacted with varying uniaxial loads ranging from 20 MPa up to 400 MPa. Their elastic properties have been measured using compressional and shear waves in confining hydrostatic pressures up to 475 MPa. Measured seismic velocities are shown to depend on 1) the ratio of calcite to muscovite, 2) the uniaxial load used during sample manufacturing, and 3) the porosity amount prior to velocity measurements. The matrix framework may also affect the seismic velocities, but this effect is not easily quantifiable. In general, measured seismic velocities decrease with increasing muscovite content and porosity. Elastic properties have been calculated based on texture measurements, which were obtained by neutron diffraction. The calculated velocities are based on the calcite and muscovite single crystal elastic tensors, their orientation distribution functions, and their modal fractions. A large discrepancy is apparent between calculated and measured velocities, where the former always overestimate the actual velocities. Shear waves display less of a difference between calculated and measured values, than do compressional waves, indicating that waves propagating by shearing are less affected by sample porosity and matrix framework. A better agreement between measured and calculated seismic velocities arises when the initial porosity is taken into account, using a differential effective medium model. Seismic anisotropy is evident from both measured and calculated velocities, and is low for samples containing pure calcite, but becomes prominent as the muscovite concentration increases. The intensity of anisotropy further depends on the uniaxial load used during sample compaction and the initial porosity.
ISSN: 0148-0227

hal-00512666
https://hal.archives-ouvertes.fr/hal-00512666
DOI : 10.1029/2009JB006523

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
We analyse P-wave receiver functions across the western Gulf of Aden and southern Red Sea continental margins in Western Yemen to constrain crustal thickness, internal crustal structure and the bulk seismic velocity characteristics in order to address the role of magmatism, faulting and mechanical crustal thinning during continental breakup. We analyse teleseismic data from 21 stations forming the temporary Young Conjugate Margins Laboratory (YOCMAL) network together with GFZ and Yemeni permanent stations. Analysis of computed receiver functions shows that (1) the thickness of unextended crust on the Yemen plateau is ∼35km; (2) this thins to ∼22km in coastal areas and reaches less than 14km on the Red Sea coast, where presence of a high-velocity lower crust is evident. The average Vp/Vs ratio for the western Yemen Plateau is 1.79, increasing to ∼1.92 near the Red Sea coast and decreasing to 1.68 for those stations located on or near the granitic rocks. Thinning of the crust, and by inference extension, occurs over a ∼130-km-wide transition zone from the Red Sea and Gulf of Aden coasts to the edges of the Yemen plateau. Thinning of continental crust is particularly localized in a <30-km-wide zone near the coastline, spatially co-incident with addition of magmatic underplate to the lower crust, above which on the surface we observe the presence of seaward dipping reflectors (SDRs) and thickened Oligo-Miocene syn-rift basaltic flows. Our results strongly suggest the presence of high-velocity mafic intrusions in the lower crust, which are likely either synrift magmatic intrusion into continental lower crust or alternatively depleted upper mantle underplated to the base of the crust during the eruption of the SDRs. Our results also point towards a regional breakup history in which the onset of rifting was synchronous along the western Gulf of Aden and southern Red Sea volcanic margins followed by a second phase of extension along the Red Sea margin.
ISSN: 0956-540X

hal-00825194
https://hal.archives-ouvertes.fr/hal-00825194
https://hal.archives-ouvertes.fr/hal-00825194/document
https://hal.archives-ouvertes.fr/hal-00825194/file/RF_Hakim2012_GJI.pdf
DOI : 10.1093/gji/ggt072

Éditeur(s) : John Wiley & Sons
Résumé : Based on all available seismic refraction data, the seismic structure of the Western Approaches to the English Channel and the South Armorican continental shelf has been established. The correlation of seismic velocities showed a three-fold structure, related to geological layers as follows: 1. High-velocity, deep sequence (4.4-6.6km/sec): metamorphic basement and Lower to Middle Palaeozoic. 2. Intermediate sequence (3.6-4.3km/sec): Permo-Triassic (and Jurassic?) and possibly also Lower Cretaceous on the margin of the Western Approaches. 3. Upper sequence (1.9-3.6km/sec): Cretaceous and Cenozoic sediments. An acoustic unconformity appears to be present between the Upper sequence and the High velocity and Intermediate sequences. Different geomorphology characterizes these layers: the Palaeozoic is tectonized and generally follows the Hercynian structural trends; mainly Permo-Triassic sediments fill the depressions of the Palaeozoic floor, but only in the Western Approaches; eastwards-transgressive Cretaceous and Cenozoic sediments cover discordantly the preceding layers. [NOT CONTROLLED OCR]
Petroleum and the continental shelf of North-west Europe (John Wiley & Sons), 1975 , P. 29-43

Droits : John Wiley & Sons
http://archimer.ifremer.fr/doc/1975/publication-5465.pdf
http://archimer.ifremer.fr/doc/00000/5465/

Éditeur(s) : HAL CCSD Oxford University Press (OUP): Policy P - Oxford Open Option A Oxford University Press (OUP)
Résumé : International audience
The Pyrenean mountain range is a slowly deforming belt with continuous and moderate seismic activity. To quantify its deformation field, we present the velocity field estimated from a GPS survey of the Pyrenees spanning 18 yr. The PotSis and ResPyr networks, including a total of 85 GPS sites, were installed and first measured in 1992 and 1995–1997, respectively, and remeasured in 2008 and 2010. We obtain a deformation field with velocities less than 1 mm yr−1 across the range. The estimated velocities for individual stations do not differ significantly from zero with 95 per cent confidence. Even so, we estimate a maximum extensional horizontal strain rate of 2.0 ± 1.7 nanostrain per year in a N–S direction in the western part of the range. We do not interpret the vertical displacements due to their large uncertainties. In order to compare the horizontal strain rates with the seismic activity, we analyse a set of 194 focal mechanisms using three methods: (i) the ‘r’ factor relating their P and T axes, (ii) the stress tensors obtained by fault slip inversion and (iii) the strain-rate tensors. Stress and strain-rate tensors are estimated for: (i) the whole data set, (ii) the eastern and western parts of the range separately, and (iii) eight zones, which are defined based on the seismicity and the tectonic patterns of the Pyrenees. Each of these analyses reveals a lateral variation of the deformation style from compression and extension in the east to extension and strike-slip in the west of the range. Although the horizontal components of the strain-rate tensors estimated from the seismic data are slightly smaller in magnitude than those computed from the GPS velocity field, they are consistent within the 2σ uncertainties. Furthermore, the orientations of their principal axes agree with the mapped active faults.
ISSN: 0956-540X

hal-01171752
https://hal.archives-ouvertes.fr/hal-01171752
DOI : 10.1093/gji/ggv052

É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 Elsevier
Résumé : International audience
Melt generation in the continental crust is thought to significantly influence seismic bulk velocities and anisotropy, although existing laboratory data provide limited constraints on such seismic attributes. In this study we measured ultrasonic compressional wave speeds (Vp) and anisotropy during sample compaction, mica-breakdown and melt-generation in synthetic, foliated quartz–muscovite aggregates. Measurements were performed at peak conditions of 300 MPa hydrostatic confining pressure and 750 °C, over a six hour time period, for three separate samples where wave propagation directions were aligned at 0°, 45° and 90° with respect to the foliation plane. The experiments are initially marked by sample compaction and rapid reduction in porosity (from 25% to <2 vol%), with a corresponding increase in Vp. During initial stages of the experiment Vp ranges from 4.65 to 5.45 km/s, depending on the direction of the P wave propagation, resulting in up to ∼16% anisotropy, which originates mainly from the preferred orientation of muscovite and possible incipient melt lenses oriented parallel to foliation. As the experiment progresses the bulk Vp decreases gradually and Vp anisotropy drops from ∼16% to ∼2% after five hours into the experiment. The drop in anisotropy is caused by muscovite-breakdown and increasing amount of partial melting of the aggregate. The bulk velocities are considerably lower than predicted for a mixture of muscovite and β-quartz (stable at 300 MPa and 750 °C). Thermodynamic modeling indicates that mica breakdown and production of melt with ∼2 wt% dissolved water in the aggregate can explain the low Vp, even though β-quartz is stable; the modeled Vp is ∼4.75 km/s at 300 MPa and 750 °C. These ultrasonic attributes signify movement of the melt to fill vacant pores and coating grain boundaries, effectively connecting melt films and reducing seismic anisotropy. The results may apply to regions where crustal melting is on-going, in particular regions of elevated temperature that are only weakly deforming, with coupled anomalously low bulk velocities and seismic anisotropy.
ISSN: 0012-821X

hal-01243569
https://hal.archives-ouvertes.fr/hal-01243569
DOI : 10.1016/j.epsl.2015.05.039

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
While vertical motion induced by long-term geological loads is often used to estimate the flexural rigidity of the lithosphere, we intend to evaluate the shear rigidity of the lithosphere using horizontal motion. Our approach considers that the rigidity of the lithosphere may be defined as its resistance to horizontal tectonic lateral forces. In this case, a spatial distribution of an effective shear rigidity can be estimated from the analysis of the interseismic velocity fields. We consider the Western United States zone where weakly strained areas (e.g., the Sierra Nevada) are connected with areas of large strain rate (e.g. San Andreas Fault system). By inverting interseismic strain distribution measured by geodetic methods, we infer the effective shear rigidity of the lithosphere. The forward problem is defined using the equations of linear elasticity. The inversion relies on the minimization of the sum of a quadratic measure of the differences between measured and modelled velocity fields. The functional also includes regularization terms for the parameters of the model. The gradient of the functional with respect to the minimization parameters is computed using an adjoint formulation. This permits the treatment of large dimensional minimization problems. Finally, a measure of the uncertainty of our inversion is illustrated through the covariance matrix of the parameters at the optimum. The optimization chart is validated on two synthetic velocity distributions. Then, the effective shear rigidity variations of the Western United States are estimated using the CMM3 interseismic velocities. The inversion displays low effective rigidities along the San Andreas Fault system, the Mojave Desert and in the Eastern California Shear Zone, while rigid areas are found in the Sierra Nevada and in the South Basin and Range. Finally, we discuss the differences between our strain rate and rigidity maps with previously published results for the Western United States.
ISSN: 0040-1951

Droits : http://hal.archives-ouvertes.fr/licences/copyright/
hal-01505024
https://hal.archives-ouvertes.fr/hal-01505024
https://hal.archives-ouvertes.fr/hal-01505024/document
https://hal.archives-ouvertes.fr/hal-01505024/file/Paper_tectophysics.pdf
DOI : 10.1016/j.tecto.2017.03.015

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
We image the lithospheric and upper asthenospheric structure of western continental Yemen with 24 broadband stations to evaluate the role of the Afar plume on the evolution of the continental margin and its extent eastwards along the Gulf of Aden. We use teleseismic tomography to compute relative P-wave velocity variations in southwestern Yemen down to 300 km depth. Published receiver function analysis suggest a dramatic and localized thinning of the crust in the vicinity of the Red Sea and the Gulf of Aden, consistent with the velocity structure that we retrieve in our model. The mantle part of the model is dominated by the presence of a low velocity anomaly in which we infer partial melting just below thick Oligocene flood basalts and recent off-axis volcanic events (from 15 Ma to present). This low velocity anomaly could correspond to an abnormally hot mantle and could be responsible for dynamic topography and recent magmatism in western Yemen. Our new P-wave velocity model beneath western Yemen suggests the young rift flank volcanoes beneath margins and on the flanks of the Red Sea rift are caused by focused small-scale diapiric upwelling from a broad region of hot mantle beneath the area. Our work shows that relatively hot mantle, along with partial melting of the mantle, can persist beneath rifted margins after breakup has occurred.
ISSN: 1525-2027

hal-01015496
https://hal.archives-ouvertes.fr/hal-01015496
https://hal.archives-ouvertes.fr/hal-01015496/document
https://hal.archives-ouvertes.fr/hal-01015496/file/Korostelev_etal_GGG_2014.pdf
DOI : 10.1002/2014GC005316