É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 American Geophysical Union
Résumé : International audience
The majority of the largest subduction megathrust earthquakes share the common characteristic of rupturing more than one asperity along strike of the margin. Understanding the factors that control coseismic failure of multiple asperities, and thus maximum magnitude, is central for seismic hazard assessment. To investigate the role of asperities size and spacing on maximum magnitude, seismicity rate, and percentage of synchronized ruptures, we use analog models simulating along-strike rupture behavior of megathrust earthquakes. We found negative correlations between the barrier-to-asperity length ratio Db/Da and maximum magnitude and seismicity rate. Db/Da also controls the process of asperities synchronization along the megathrust. A permanent barrier behavior is observed for Db/Da > 0.5. Comparing our experimental results to the Nankai Trough historical seismicity, we propose that the distribution of megathrust frictional heterogeneities likely explains the diversity of earthquakes which occurred there.
ISSN: 0094-8276

hal-01685546
https://hal.archives-ouvertes.fr/hal-01685546
https://hal.archives-ouvertes.fr/hal-01685546/document
https://hal.archives-ouvertes.fr/hal-01685546/file/Corbi_et_al-2017-Geophysical_Research_Letters.pdf
DOI : 10.1002/2017GL074182

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
We report measurements of long-period P wave polarization (P pol) in Australia and Tahiti made by combining modeling of the polarization deviation and harmonic analysis. The analysis of the deviation of the horizontal polarization of the P wave as a function of event back azimuth may be used to obtain information about (1) sensor misorientation, (2) dipping discontinuities, (3) seismic anisotropy, and (4) velocity heterogeneities beneath a seismic station. The results from harmonic analysis and a grid search using Snell's law suggest the presence of a dipping seismic discontinuity beneath stations CTAO and CAN in Australia. These results are consistent with published receiver function studies for these stations. The P pol fast axis orientation is close to the N–S absolute plate motion direction at station TAU (Tasmania), which may be due to plate-motion-driven alignment of olivine crystals in the asthenosphere. Interestingly, measurements of SKS splitting at Tahiti (French Polynesia) show an apparent isotropy, whereas an inversion of P pol observations at PPTL seismic station located in Tahiti suggests the presence of two anisotropic layers. The fast axis azimuth is oriented E–W in the upper layer, and it is close to the NW–SE orientation in the lower layer. Since P pol orientations are used for real-time earthquake locations, especially in poorly instrumented areas such as the South Pacific, we show that the bias from anisotropy and sensor misorientation determined here can be corrected to improve the location accuracy, which yields fundamental data for rapid location necessary for effective tsunami warning.
ISSN: 2169-9313

hal-01389113
http://hal.univ-reunion.fr/hal-01389113
http://hal.univ-reunion.fr/hal-01389113/document
http://hal.univ-reunion.fr/hal-01389113/file/Fontaine_Ppol_JGR2009.pdf
DOI : 10.1029/2008JB005709

Éditeur(s) : HAL CCSD Mineralogical Society of America
Résumé : International audience
In situ mineralization of CO2 in ultramafic rock-hosted aquifers is one of the promising solutions for decreasing CO2 concentrations in the atmosphere. Naturally altered ultramafic rocks suggest that carbonation processes are controlled by local heterogeneities in the structure of the rock and fluid transport at the water-rock interfaces. We studied the role of rock crystallographic anisotropy relative to the global fluid flow direction on the mineralization of CO2 by means of electron microscope analyses from the macro- to the micrometer scale (EBSD-FIB). The sample used for the measurements was a hot pressed olivine core percolated by water enriched in CO2 (pCO2 = 10 MPa) at 180 °C. During the percolation experiment, olivine was dissolved and two types of carbonates, dolomite, and magnesite, were precipitated on olivine surfaces. The results showed that the dissolution of olivine is controlled by its crystallographic properties as shown by the development of etch-pits only on the (010)ol planes and with elongated shapes parallel to the [010]ol axes. In contrast, the precipitation of carbonates is governed by hydrodynamic properties. Carbonates are heterogeneously distributed in the percolated rock. They are mainly located along the moderate (for dolomite) and the minor (for magnesite) flow paths, both oriented parallel to the principal fluid flow direction, which allow carbonates to be supplied with divalent cations (e.g., Ca2+, Mg2+, and Fe2+). In these flow paths, carbonate growth is systematically oriented normal to the flow that facilitates the development of chemical gradients with cationic supersaturation conditions for carbonate precipitation near the walls.In natural systems, the (010)ol planes are parallel to the Moho and the (100)ol planes are vertical; our study suggests that flow of CO2-rich fluids will induce precipitation of carbonates localized along, and preferentially clogging, vertical flow paths while favoring olivine dissolution along horizontal fluid pathways. This dual control of structure and fluid flow on carbonation mechanisms could be an important parameter allowing sustainable CO2 storage in peridotites, while limiting the risks of leakage toward the surface.
ISSN: 0003-004X

hal-01143393
https://hal.archives-ouvertes.fr/hal-01143393
DOI : 10.2138/am-2015-4913

Éditeur(s) : HAL CCSD Public Library of Science
Résumé : International audience
Local variation in the density of Anopheles mosquitoes and the risk of exposure to bites are essential to explain the spatial and temporal heterogeneities in the transmission of malaria. Vector distribution is driven by environmental factors. Based on variables derived from satellite imagery and meteorological observations, this study aimed to dynamically model and map the densities of Anopheles darlingi in the municipality of Saint-Georges de l’Oyapock (French Guiana). Longitudinal sampling sessions of An. darlingi densities were conducted between September 2012 and October 2014. Landscape and meteorological data were collected and processed to extract a panel of variables that were potentially related to An. darlingi ecology. Based on these data, a robust methodology was formed to estimate a statistical predictive model of the spatial-temporal variations in the densities of An. darlingi in Saint-Georges de l’Oyapock. The final cross-validated model integrated two landscape variables—dense forest surface and built surface—together with four meteorological variables related to rainfall, evapotranspiration, and the minimal and maximal temperatures. Extrapolation of the model allowed the generation of predictive weekly maps of An. darlingi densities at a resolution of 10-m. Our results supported the use of satellite imagery and meteorological data to predict malaria vector densities. Such fine-scale modeling approach might be a useful tool for health authorities to plan control strategies and social communication in a cost-effective, targeted, and timely manner.
ISSN: 1932-6203

hal-01383486
https://hal.archives-ouvertes.fr/hal-01383486
https://hal.archives-ouvertes.fr/hal-01383486/document
https://hal.archives-ouvertes.fr/hal-01383486/file/Adde_et_al_PlosONE_2016.pdf
DOI : 10.1371/journal.pone.0164685

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Past studies have shown that high coastal uplift rates are restricted to active areas, especially in a subduction context. The origin of coastal uplift in subduction zones, however, has not yet been globally investigated. Quaternary shorelines correlated to the last interglacial maximum (MIS 5e) were defined as a global tectonic benchmark (Pedoja et al. (2011)). In order to investigate the relationships between the vertical motion and the subduction dynamic parameters, we cross-linked this coastal uplift database with the "geodynamical" databases from Heuret (2005), Conrad and Husson (2009) and Müller et al. (2008). Our statistical study shows that: [1] the most intuitive parameters one can think responsible for coastal uplift (e.g., subduction obliquity, trench motion, oceanic crust age, interplate friction and force, convergence variation, dynamic topography, overriding and subducted plate velocity) are not related with the uplift (and its magnitude); [2] the only intuitive parameter is the distance to the trench which shows in specific areas a decrease from the trench up to a distance of ~300 km; [3] the slab dip (especially the deep slab dip), the position along the trench and the overriding plate tectonic regime are correlated with the coastal uplift, probably reflecting transient changes in subduction parameters. Finally we conclude that the first order parameter explaining coastal uplift is small-scale heterogeneities of the subducting plate, as for instance subducting aseismic ridges. The influence of large-scale geodynamic setting of subduction zones is secondary.
ISSN: 0264-3707

hal-00968538
https://hal.archives-ouvertes.fr/hal-00968538
https://hal.archives-ouvertes.fr/hal-00968538/document
https://hal.archives-ouvertes.fr/hal-00968538/file/henryterraces_v2.pdf
DOI : 10.1016/j.jog.2014.04.001

Éditeur(s) : HAL CCSD Frontiers Media
Résumé : International audience
We present high resolution observations of microstructure and texture evolution during dynamic recrystallization (DRX) of ice polycrystals deformed in the laboratory at high temperature (≈0.98Tm). Ice possesses a significant viscoplastic anisotropy that induces strong strain heterogeneities, which result in an early occurrence of DRX mechanisms. It is therefore a model material to explore these mechanisms. High resolution c-axis measurements at sample scale by optical techniques and full crystallographic orientation measurements by cryo-Electron Back Scattering Diffraction (EBSD) provide a solid database for analyzing the relative impact of the macroscopic imposed stress vs. the local and internal stress fields on DRX mechanisms. Analysis of misorientation gradients in the EBSD data highlights a heterogeneous dislocation distribution, which is quantified by the Nye tensor estimation. Joint analyses of the dislocation density maps and microstructural observations highlight spatial correlation between high dislocation density sites and the onset of nucleation taking place by grain-boundary bulging, subgrain rotation or by the formation of kink-bands.
ISSN: 2296-6463

hal-01515284
https://hal.univ-lorraine.fr/hal-01515284
DOI : 10.3389/feart.2015.00081

Éditeur(s) : HAL CCSD Springer Verlag
Résumé : Limestone dissolution by CO2 -rich brine induces critical changes of the pore network geometrical parameters such as the pore size distribution, the connectivity, and the tortuosity which govern the macroscopic transport properties (permeability and dispersivity) that are required to parameterize the models, simulating the injection and the fate of CO2 . A set of four reactive core-flood experiments reproducing underground conditions ( T=100∘C and P=12 MPa) has been conducted for different CO2 partial pressures (0.034
ISSN: 0169-3913

hal-00963597
https://hal.archives-ouvertes.fr/hal-00963597
DOI : 10.1007/s11242-013-0257-4

Éditeur(s) : HAL CCSD Elsevier
Résumé : The correct characterization of the effective reactive transport dynamics is an important issue for modeling reactive transport on the Darcy scale, specifically in situations in which reactions are localized, that is when different reactions occur in different portions of the porous medium. Under such conditions the conventional approach of homogenizing only the porous medium chemistry is not appropriate. We consider here reactive transport in a porous medium that is characterized by mass transfer between a mobile and a distribution of immobile regions. Chemical and physical heterogeneities are reflected by distributions of kinetic reaction rate constants and residence times in the immobile zones. We derive an effective reactive transport equation for the mobile solute that is characterized by non-local physical mass transfer and reaction terms. Specifically, chemical heterogeneity is upscaled in terms of a reactive memory function that integrates both chemical and physical heterogeneity. Mass transfer limitations due to physical heterogeneity yield effective kinetic rate coefficients that can be much smaller than the volumetric average of the local scale coefficients. These results help to explain and quantify the often reported discrepancy between observed field reaction rate constants and the ones obtained under well mixed laboratory conditions. Furthermore, these results indicate that transport under physical and chemical heterogeneity cannot be upscaled separately.
ISSN: 0169-7722

hal-00617674
https://hal.archives-ouvertes.fr/hal-00617674
DOI : 10.1016/j.jconhyd.2010.06.002

Éditeur(s) : HAL CCSD Springer Verlag
Résumé : International audience
Random walk (RW) methods are recurring Monte Carlo methods used to modelconvective and diffusive transport in complex heterogeneous media. Many applications canbe found, including fluid mechanic, hydrology and chemical reactors modeling. These methodsare easy to implement, very versatile and flexible enough to become appealing for manyapplications because they generally overlook or deeply simplify the building of explicit complexmeshes required by deterministic methods. RWprovides a good physical understandingof the interactions between the space scales of heterogeneities and the transport phenomenaunder consideration. In addition, they can result in efficient upscaling methods, especiallyin the context of flow and transport in fractured media. In the present study, we review theapplications of RW to several situations that cope with diverse spatial scales and differentinsights into upscaling problems. The advantages and downsides of RW are also discussed,thus providing a few avenues for further works and applications.
ISSN: 0169-3913

insu-01308060
https://hal-insu.archives-ouvertes.fr/insu-01308060
DOI : 10.1007/s11242-016-0693-z

Éditeur(s) : HAL CCSD Springer Verlag
Résumé : International audience
Chromitite pods in the Mayarí-Cristal ophiolitic massif (eastern Cuba) were formed in the Late Cretaceous when island arc tholeiites and MORB-like back-arc basin basalts reacted with residual mantle peridotites and generated chromite-rich bodies enclosed in dunite envelopes. Platinum-group minerals (PGM) in the podiform chromitites exhibit important Os-isotope heterogeneities at the kilometric, hand sample and thin section scales. 187Os/188Os calculated at the time of chromitite crystallization (~90 Ma) ranges between 0.1185 and 0.1295 (γOs = −7.1 to +1.6, relative to enstatite chondrite), and all but one PGM have subchondritic 187Os/188Os. Grains in a single hand sample have initial 187Os/188Os that spans from 0.1185 to 0.1274, and in one thin section it varies between 0.1185 and 0.1232 in two PGM included in chromite which are only several millimeters apart. As the Os budget of a single micrometric grain derives from a mantle region that was at least several m3 in size, the variable Os isotopic composition of PGM in the Mayarí-Cristal chromitites probably reflects the heterogeneity of their mantle sources on the 10-100 m scale. Our results show that this heterogeneity was not erased by pooling and mingling of individual melt batches during chromitite crystallization but was transferred to the ore deposits on mineral scale. The distribution of the Os model ages calculated for PGM shows four main peaks, at ~100, 500, 750 and 1,000 Ma. These variable Os model ages reflect the presence of different depleted domains in the oceanic (Pacific-related) upper mantle of the Greater Antilles paleo-subduction zone. The concordance between the age of crystallization of the Mayarí-Cristal chromitites and the most recent peak of the Os model age distribution in PGM supports that Os in several grains was derived from fertile domains of the upper mantle, whose bulk Os isotopic composition is best approximated by that of enstatite chondrites; on the other hand, most PGM are crystallized by melts that tapped highly refractory mantle sources.
ISSN: 0010-7999

hal-00822742
https://hal.archives-ouvertes.fr/hal-00822742
DOI : 10.1007/s00410-010-0575-2

Éditeur(s) : HAL CCSD American Geophysical Union (AGU)
Résumé : Tectonics, v. 27, p. TC5011, 2008. http://dx.doi.org/10.1029/2007TC002232
International audience
The seismic activity of the Ligurian Basin, the northeastern termination of the western Mediterranean basin, is larger than in surrounding regions, even though recent geodetic studies attest that this area is subject to very low levels of deformation. This basin is an example of a type of passive margins that cannot be considered solely as inert sites of sedimentation and of progressive subsidence and that are reactivated in a compressive pattern; other examples include the Kwanza basin (Angola) and the Brazil margin. We investigated, by means of 2-D thermomechanical modeling, the structural and rheological heterogeneities that can lead to concentration of strain in this marginal basin. We deduced that the deformation of the basin is due to its particular geometric features, narrow and with a thick surrounding continental crust, related to its position at the southern termination of the Alps. This sharp transition, in terms of both geometry and rheological contrast, is a main factor in explaining the weakness of the margin. We discuss the importance of buoyancy forces versus tectonic forces, as well as thermal effects, on the observed reactivation. Influence of contrast in rheology between an oceanic-type crust and continental crust is also studied. Geodynamical implications are proposed for the region. The good agreement between the predicted localized deformation and the observed seismicity distribution should help improve seismic hazard assessment in the region.
ISSN: 0278-7407

hal-00407859
https://hal.archives-ouvertes.fr/hal-00407859
DOI : 10.1029/2007TC002232

Éditeur(s) : HAL CCSD Springer
Résumé : The analysis of gamma ray measurements provides valuable information about the basement's lithospherical model and the radioactive elements contained in the subsurface. To be exploited, the raw data need pre-processing operations which consist of corrections of the Compton effect, altimetry and the aircraft background, and 2D interpolations. These transformations may modify the data's stochastic component which reflects heterogeneities of the investigated medium. In this study, a fractal analysis is carried out on airborne gamma ray data measured in the Hoggar's area (Algeria). First, it is shown that the analyzed data exhibit fractal properties. Due to the dependence of the latter on the spatial location of the measures, the data can be considered as paths of multifractional Brownian motions. Another result worth noting is that the fractal nature of the raw data is not modified by the pre-processing operations. The fractal analysis can be then performed on the raw spectrometric measurements without submitting them to any pre-processing; thus, the stochastic component of the measurements is kept intact.
ISSN: 1866-7411

hal-00817494
https://hal.archives-ouvertes.fr/hal-00817494
DOI : 10.1007/s12517-011-0407-3

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Geochemical studies of primitive basalts have documented the presence of crustal-derived garnet pyroxenite in their mantle sources. The processes whereby melts with the signature of garnet pyroxenite are produced in the mantle are, however, poorly understood and somewhat controversial. Here we investigate a natural example of the interaction between melts of garnet pyroxenite derived from recycled plagioclase-rich crust and surrounding mantle in the Ronda peridotite massif. Melting of garnet pyroxenite at ∼1.5 GPa generated spinel websterite residues with MREE/HREE fractionation and preserved the positive Eu anomaly of their garnet pyroxenite precursor in whole-rock and clinopyroxene. Reaction of melts from garnet pyroxenite with depleted surrounding peridotite generated secondary fertile spinel lherzolite. These secondary lherzolites differ from common spinel lherzolite from Ronda and elsewhere by their lower-Mg# in clinopyroxene, orthopyroxene and olivine, lower-Cr# in spinel and higher whole-rock Al2O3, CaO, Sm/Yb and FeO⁎ at a given SiO2. Remarkably, secondary spinel lherzolite shows the geochemical signature of ghost plagioclase in the form of positive Eu and Sr anomalies in whole-rock and clinopyroxene, reflecting the transfer of a low-pressure crustal imprint from recycled pyroxenite to hybridized peridotite. Garnet pyroxenite melting and melt-peridotite interaction, as shown in the Ronda massif, may explain how the signature of subducted or delaminated crust is transferred to the mantle and how a garnet pyroxenite component is introduced into the source region of basalts. The efficiency of these processes in conveying the geochemical imprint of crustal-derived garnet pyroxenite to extruded lavas depends on the reactivity of pyroxenite melt with peridotite and the mantle permeability, which may be controlled by prior refertilization reactions similar to those documented in the Ronda massif. Highly fertile heterogeneities produced by pyroxenite-peridotite interaction, such as secondary spinel lherzolite in Ronda, may nucleate magmatic channels that remain chemically isolated from the ambient mantle and act as preferential pathways for melts with the signature of recycled crust.
ISSN: 0012-821X

hal-00814010
https://hal.archives-ouvertes.fr/hal-00814010
DOI : 10.1016/j.epsl.2012.11.047

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Microstructures, U-Pb SHRIMP zircon ages and temperature estimates were used to constrain the fabric evolution and timing of mylonitisation in the Patos shear zone (NE Brazil). Melt-bearing HT-mylonitic orthogneiss displays solid-state fabrics with coarse quartz ribbons and sutured grain boundaries. K-feldspar is often fractured and shows peripheral myrmekite. The textures typically become magmatic at the contacts with the Espinho Branco anatexite. Zircon crystals from two samples of the leucosome show Paleoproterozoic (ca. 2.2 Ga) inherited cores that are enveloped by recrystallised Neoproterozoic rims. These zircon grains define a discordia with lower intercept ages of 558 and 562 Ma but with large analytical errors. A leucogranite with transitional contacts with host diatexites provides a well-constrained zircon mean age of 566 ± 6 Ma (n = 15, 95% conf.), which is considered the best estimate for the peak metamorphic conditions. The transition from HT-mylonites to fine-grained mylonites and ultramylonites is marked by a progressive grain size reduction assisted by recrystallisation. Quartz crystallographic fabrics display [0 0 0 1] concentrations between Z and Y. K-feldspar and plagioclase fabrics record the activity of the (0 1 0)[0 0 1] and (0 1 0)[1 0 0] slip systems and, in the lower-temperature tectonites, the (1 0 0)[0 1 0] slip system. These results suggest a diachronic evolution in which an initial HT northern domain was deformed by intracrystalline slip and melt-assisted granular flow, followed by medium- to low-temperature solid-state creep in the southern shear zone domain. The mylonites therefore record the rheological heterogeneities responsible for the formation of a late continental shear zone that nucleated under partially molten conditions and was later reactivated through melt-absent strain localization following exhumation.
ISSN: 0301-9268

hal-01053497
https://hal.archives-ouvertes.fr/hal-01053497
DOI : 10.1016/j.precamres.2013.12.020

Éditeur(s) : HAL CCSD Wiley
Résumé : International audience
Scientific interest in urban meteorology has increased because highly populated areas experience high vulnerability to pollution or heavy rain. However, compared to urban air quality or urban heat island (UHI) processes, the urban water vapour cycle is poorly understood because it has been investigated less due to the lack of upper-air measurements and the high sensitivity of surface measurements to local heterogeneities. In this paper, surface measurements of wind, temperature, pressure and humidity, as well as integrated water vapour (IWV) from GPS and MODIS and numerical simulations, have been used to investigate the urban cycle of water vapour in May and June 2004 during the VAPIC field experiment in the Paris area. The surface data show the typical characteristics of an urban area with the absence of water vapour sources and a UHI of about 6 °C at night. The urban IWV distribution differs completely, with an urban IWV excess on average between 1600 and 0600 UTC (with a maximum of about 1.5 kg m−2 at 0600 and 1700 UTC). No IWV difference between the urban and rural areas is found in the middle of the day. The numerical simulations reproduce accurately the urban IWV anomaly. Shallow surface wind convergence associated with the UHI during nighttime provides moisture from the rural areas. Urban areas are therefore under wind convergence for most of the time. The rural water vapour sources and the depth of the convergence control the amplitude of the urban IWV excess. At about 1200 UTC, entrainment at the top of the urban boundary layer is the key process that inhibits the urban IWV excess observed at night
ISSN: 0035-9009

hal-00406370
https://hal.archives-ouvertes.fr/hal-00406370
DOI : 10.1002/qj.446

Éditeur(s) : HAL CCSD
Résumé : In coastal carbonate aquifers submitted to a saline intrusion, water - rock interactions processes at the interface between freshwater and saltwater are a cause for heterogeneities development within the porous media. The dissolutions and precipitations induced by the mixing zone play a major role in the transport properties of the reservoir and so in the long term evolution of the saltwater intrusion. The characterization and the in situ observation of a current saline intrusion in the Llucmajor carbonate platform at the experimental site of Ses Sitjoles (South-Easth of Mallorca Island, Spain), allowed to study the temporal evolution of the saltwater intrusion and the biogeochemical processes active in the mixing zone. The influence of microstructure, characterized using X-Ray microtomography, on the petrophysical and flow properties was investigated in laboratory and flow-through experiments were also conducted to study the modifications of the microstructures when submitted to water - rock interactions processes. The mixing zone at the Ses Sitjoles site appears as quite thick and not very sensitive to climate and anthropic stress and is proved to be the locus of active microbiological activity. Laboratory studies highlight a strong control of the microporosity on transport processes and of pCO2 on water - rock interactions. The percolation experiments of slightly reactive waters reveal calcite dissolution together with a decrease in permeability due to a rearrangement of micro-grains within the microporous network.
Les processus d'interactions fluide-roche à l'interface eau douce-eau salée, dans les aquifères carbonatés côtiers soumis à une intrusion saline, sont une des causes du développement d'hétérogénéités dans la structure du milieu poreux. Les dissolutions et/ou précipitations engendrées par la présence de cette zone de mélange jouent un rôle déterminant pour les propriétés de transport du réservoir et donc l'évolution à long terme de l'intrusion d'eau salée. En prenant l'exemple de la plateforme carbonatée de Llucmajor au niveau du site expérimental de Ses Sitjoles (Sud-Est de Majorque, Espagne), un travail de caractérisation et d'observation in situ à moyen terme de l'intrusion saline actuelle de ce réservoir permet de décrire son évolution temporelle et les processus biogéochimiques actifs dans la zone de mélange. L'influence de la microstructure du milieu poreux, caractérisée par microtomographie RX, sur les propriétés pétrophysiques et hydrodynamiques est étudiée en laboratoire et des expériences de percolations d'eaux réactives (mélanges eau douce - eau salée et eau pure chargée en CO2) sont réalisées afin d'étudier les modifications de ces microstructures lorsqu'elles sont soumises à des interactions eau - roche. La zone de mélange eau douce - eau salée rencontrée au niveau du site de Ses Sitjoles est large et peu sensible aux contraintes extérieures (climatologiques et anthropiques) et constitue un lieu d'activité microbiologique particulièrement actif. En laboratoire, on observe un fort contrôle des phénomènes de transport par la microporosité et des processus d'interactions eau - roche par la pCO2. Dans le cas d'une percolation de fluides moyennement réactifs, les expériences mettent en évidence une dissolution de calcite accompagnée par une diminution de la perméabilité, liée à un réarrangement des micro-grains à l'intérieur du réseau microporeux.
https://tel.archives-ouvertes.fr/tel-00767152

tel-00767152
https://tel.archives-ouvertes.fr/tel-00767152
https://tel.archives-ouvertes.fr/tel-00767152/document
https://tel.archives-ouvertes.fr/tel-00767152/file/thesegaring.pdf

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
To evaluate the performance of downhole and surface geophysical monitoring methods, a series of shallow gas injection-monitoring experiments has been performed in a coastal saline aquifer at Maguelone, France. The recorded data include pressure measurements with a Westbay multilevel completion and CO2 saturation at an observation well derived from electrical resistivity with a modified Waxman-Smits (MWS) model. In this work, the aim is to develop a simulation model capturing the gas transport behavior and consistent with field data. For this purpose, the simulation of the CO2 injection experiment is carried out with two conceptual models, a homogeneous model and a heterogeneous model treated with multiple realization Monte Carlo simulations. Numerical simulator TOUGH2 with the equation of state module EOS7C is used for the simulations. Comparison of the model results with field data suggests that the pressure responses are captured with relatively good accuracy. Similarly, the model also provides an overall reasonable agreement and correct order of magnitude for predicted gas saturation values. However, as the heterogeneity pattern in the field data remains largely unknown, the model predictions can only be used to capture the mean behavior as well as to provide insights into how heterogeneity can influence the system behavior, by means of sensitivity analyses of the influence of heterogeneities on individual realizations.
ISSN: 1750-5836

hal-01467240
https://hal.archives-ouvertes.fr/hal-01467240
DOI : 10.1016/j.ijggc.2016.09.006

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
Continental rifts and passive continental margins show fundamental along-axis segmentation patterns that have been attributed to one or a number of different processes: extensional fault geometry, variable stretching along strike, preexisting lithospheric compositional and structural heterogeneities, oblique rifting, and the presence or absence of eruptive volcanic centers. The length and width scales of the rift stage fault-bounded basin systems change during the late evolution of the new plate boundary, and the role of magmatism may increase as rifting progresses to continental rupture. Along obliquely spreading ridges, first-order mid-ocean ridge geometries originate during the synrift stage, indicating an intimate relationship between magma production and transform fault spacing and location. The Gulf of Aden rift is a young ocean basin in which the earliest synrift to breakup structures are well exposed onshore and covered by thin sediment layers offshore. This obliquely spreading rift is considered magma-poor and has several large-offset transforms that originated during late stage rifting and control the first-order axial segmentation of the spreading ridge. Widely spaced geophysical transects of passive margins that produce only isolated 2-D images of crust and uppermost mantle structure are inadequate for evaluation of competing rift evolution models. Using closely spaced new geophysical and geological observations from the Gulf of Aden we show that rift sectors between transforms have a large internal variability over short distances (∼10 km): the ocean-continent transition (OCT) evolves from a narrow magmatic transition to wider zones where continental mantle is probably exhumed. We suggest that this small-scale variability may be explained (1) by the distribution of volcanism and (2) by the along-strike differences in time-averaged extension rate of the oblique rift system. The volcanism may be associated with (1) the long-offset Alula-Fartak Fracture Zone, which may enhance magma production on its younger side, or (2) channeled flow from the Afar plume material along the newly formed OCT and the spreading ridge. Oblique extension and/or hot spot interactions may thereby have a significant control on the styles of rifting and continental breakup and on the evolution of many magma-poor margins.
ISSN: 1525-2027

hal-00504866
https://hal.archives-ouvertes.fr/hal-00504866
https://hal.archives-ouvertes.fr/hal-00504866/document
https://hal.archives-ouvertes.fr/hal-00504866/file/G3-Leroy_al-2010.pdf
DOI : 10.1029/2009GC002963

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
We report measurements of long-period P wave polarization (P-pol) in Australia and Tahiti made by combining modeling of the polarization deviation and harmonic analysis. The analysis of the deviation of the horizontal polarization of the P wave as a function of event back azimuth may be used to obtain information about (1) sensor misorientation, (2) dipping discontinuities, (3) seismic anisotropy, and (4) velocity heterogeneities beneath a seismic station. The results from harmonic analysis and a grid search using Snell's law suggest the presence of a dipping seismic discontinuity beneath stations CTAO and CAN in Australia. These results are consistent with published receiver function studies for these stations. The P-pol fast axis orientation is close to the N-S absolute plate motion direction at station TAU (Tasmania), which may be due to plate-motion-driven alignment of olivine crystals in the asthenosphere. Interestingly, measurements of SKS splitting at Tahiti (French Polynesia) show an apparent isotropy, whereas an inversion of P-pol observations at PPTL seismic station located in Tahiti suggests the presence of two anisotropic layers. The fast axis azimuth is oriented E-Win the upper layer, and it is close to the NW-SE orientation in the lower layer. Since P-pol orientations are used for real-time earthquake locations, especially in poorly instrumented areas such as the South Pacific, we show that the bias from anisotropy and sensor misorientation determined here can be corrected to improve the location accuracy, which yields fundamental data for rapid location necessary for effective tsunami warning.
ISSN: 2169-9313

hal-00413141
https://hal.archives-ouvertes.fr/hal-00413141
https://hal.archives-ouvertes.fr/hal-00413141/document
https://hal.archives-ouvertes.fr/hal-00413141/file/FontaineetalJGR2008.pdf
DOI : 10.1029/2008JB005709