Éditeur(s) : HAL CCSD
Résumé : International audience
The purpose of our study is to model a geothermal reservoir. When geothermal reservoir are assumed to be composed of pure water, the transfer of mass and energy is classically described by two balance equations: mass balance equation and the energy balance equation. In addition to those equations, fluid velocity ist classically given by the Darcy law while thermodynamic properties, inferred from theoretical or empirical equations of state, are used to close the mathematical system. Once this system is closed, there exist different solutions. The first one is to solve for pressure and temperature with a variable switch to saturation in the two-phase region (e.g. TOUGH2). The second one is to solve for pressure and enthalpy to increase stability of phase transition between single and two-phase states (e.g. Hydrotherm). We adopted the second option and chose te use a splitting method to get rid of the complexity of coupling equations and a finite volume method for the spatial discretization. Selecting object-oriented languages, we developed a multi-language framework, combining Python, Fortran and a C++ implementation of IAPWS (from the freesteam project) including the supercritical equations, in porous media velocity is given by Darcy law and to close the system physical properties are determined by the IAPWS-IF97 thermodynamic formulation. We resolve the equations in pressure and enthalpy instead of pressure and temperature in order to increase stability and to handle easier the passage from a single-phase to a two-phase system. We solve the system by using a splitting method - to get rid of the complexity of coupling equations - and a finite volume method. We offer some freedom to users thanks to the implementation of several methods like explicit or implicit Euler, Runge-Kutta or BDF2 for time solvers or GMRES and BICGSTAB for the linear solver. We can handle several boundary conditions like no-flow - describing a boundary which can not exchange matter with the exterior - or like a mixed-therm condition - a Dirichlet condition to the pressure and a Dirichlet or an outflow condition to the temperature in order to describe a recharge or a discharge zone - ... We're developing a multi-language framework, combining Python, Fortran and the C++ implementation of IAPWS (from the freesteam project). All these languages are object-oriented. We've applied this simulation model to the dogger in Paris, France, to several one-dimensional systems and a two-dimensional one made by Coumou with the CSMP++ platform. The dogger is a reservoir exploited to produce heat by pumping water at 70 and reinjecting it in the reservoir at 40. In the one-dimensional systems we wanted to observe the process of heat transfer from a higher temperature boundary to a smaller one in a high-energy domain. The last simulation shows the natural convection of water in a fault. For every simulation we compared the solutions we found with an other code (TOUGH2 or CSMP++) and they agreed. The next step will be to model the geothermal plant in Guadeloupe, West Indies. It's the only place in France - and in the West Indies so far - producing electricity with the earth power. The temperature can reach up to 1000 and the pressure range is around a few hundreds MPa. In some surface zones we can see two-phase water at atmospheric pressure. In the 1980s Bouillante was a laboratory for France. Since 1995 Bouillante has given 30GWh electricity a year to the Guadeloupeans.
PROCEEDINGS, Thirty-Ninth Workshop on Geothermal Reservoir Engineering
Stanford, United States

Droits : info:eu-repo/semantics/OpenAccess
hal-00944133
https://hal-brgm.archives-ouvertes.fr/hal-00944133
https://hal-brgm.archives-ouvertes.fr/hal-00944133/document
https://hal-brgm.archives-ouvertes.fr/hal-00944133/file/Copol.pdf

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
P>The Strait of Hormuz area is a transition zone between the continental collision of the Zagros (west) and the subduction of an oceanic part of the Arabian Plate beneath the Makran wedge (east). Geology and recent GPS measurements indicate that about 15 mm yr-1 of relative motion in N10 degrees E direction is accommodated by two major fault systems: (1) the NNW-trending Minab-Zendan-Palami (MZP) fault system that connects the Main Zagros Thrust (MZT) to the inner Makran thrust system and the Frontal subduction thrust and (2) the N-trending Sabzevaran-Kahnuj-Jiroft (SKJ) fault system that bounds the Jazmurian depression to the west. We use dense GPS measurements along four transects across these fault systems in order to determine the strains spatial distribution. The northern GPS transect confirms the total fault slip rates for both fault systems estimated by the tectonic analyses (about 10 and 7.3 mm yr-1 in N10 degrees direction across the MZP and SKJ fault systems, respectively). For both fault systems, the elastic deformation spreads over shear zones that are several tens of kilometres wide. However, transects located close to latitude 27 degrees N reveal a much narrower shear zone (similar to 10 km) for the MZP fault system. Moreover, we confirm that most of the present-day strain is transferred towards the frontal subduction thrust rather than towards the inner Makran thrusts. In order to complement this new GPS velocity field with spatially dense measurements, we processed a set of ERS radar images by the radar interferometry (InSAR) technique. We used both a 'stacking' and a 'persistant-scatterers' approach to differentiate the ground deformation signal which spatial gradient is expected to be very low, from the atmospheric signal. Results from these interferograms appear to be relatively in agreement with the GPS-determined strain distribution. Nevertheless, they confirm the absence of any superficial creep behaviour since no sharp discontinuity on interferometric phase can be noted on any interferogram. Finally, we use a purely kinematic 'block model' inversion process to calculate slip rates and locking depths for each fault system from our GPS measurements. These models suggest that the relative quiescence over the last 200 yr has certainly produced a slip deficit as high as 2 m. So, we may wonder if the MZP fault system is not late in the interseismic phase of its earthquake cycle.
ISSN: 0956-540X

hal-00429386
https://hal.archives-ouvertes.fr/hal-00429386
DOI : 10.1111/j.1365-246X.2009.04321.x

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The Chalkidiki block is a major domain in the North Aegean that, contrary to other domains in the region, largely escaped thermal perturbations during Tertiary extension. As a result, the Chalkidiki block is an ideal candidate to glean information related to the timing of Mesozoic thermal events using appropriate geochronological techniques. We have undertaken a laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) study (U-Th-Pb on monazites and U-Pb on zircons) coupled with 40Ar/39Ar dating on nine samples from various structural levels within the thrust system of the Chalkidiki block. The eastern, and structurally lower part of the system revealed a complete isotopic reset of Carboniferous – Early Triassic monazites coeval with partial monazite destruction, REE-mobilisation and formation of apatite-allanite-epidote coronas at ~ 132 Ma, a reaction that is commonly observed in amphibolite-facies rocks. These coronas formed after crystallisation of garnet (i.e., at T > 580 °C) and, in all probability, either close to the peak-temperature conditions (~ 620 °C) on a prograde path or during retrogression between the peak-temperature and the low-temperature boundary of the amphibolite facies. Cooling of these rocks and arrival at mid-crustal levels occurred at 95–100 Ma. By contrast, the western, and structurally uppermost part of the system went through the same event by 120–125 Ma. Further structural considerations with respect to medium-temperature geochronology data imply that syn-metamorphic thrusting must have ceased by early Late Cretaceous. We emphasize that, with the sole exception of the Chalkidiki block, no pre-45 Ma medium-temperature geochronology data are preserved in other North Aegean domains, a feature that is clearly related to the extension-induced thermal perturbation of the region during the Tertiary.
ISSN: 0040-1951

insu-01329472
https://hal-insu.archives-ouvertes.fr/insu-01329472
https://hal-insu.archives-ouvertes.fr/insu-01329472/document
https://hal-insu.archives-ouvertes.fr/insu-01329472/file/kydonakis-tectono-2016.pdf
DOI : 10.1016/j.tecto.2016.06.006

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
We report trace element concentrations of silicate phases in chondrules from LL3 ordinary chondrites Bishunpur and Semarkona. Results are similar to previously reported data for carbonaceous chondrites, with rare earth element (REE) concentrations increasing in the sequence olivine < pyroxene < mesostasis, and heavy REE (HREE) being enriched by 1–2 orders of magnitude (CI-normalized) relative to light REE (LREE) in ferromagnesian silicates, although no single olivine with very large LREE/HREE fractionation has been found. On average, olivine in type II chondrules is poorer in refractory lithophile incompatible elements (such as REE) than its type I counterpart by a factor of ∼2. This suggests that olivine in type I and II chondrules formed by batch and fractional crystallization, respectively, implying that type II chondrules formed under faster cooling rates (>∼10 K/h) than type I chondrules. Appreciable Na concentrations (3–221 ppm) are measured in olivine from both chondrule types; type II chondrules seem to have behaved as closed systems, which may require chondrule formation in the vicinity of protoplanets or planetesimals. At any rate, higher solid concentrations in type II chondrule forming regions may explain the higher oxygen fugacities they record compared to type I chondrules. Type I and type II chondrules formed in different environments and the correlation between high solid concentrations and/or oxygen fugacities with rapid cooling rates is a key constraint that chondrule formation models must account for.
ISSN: 0016-7037

hal-01174084
https://hal.archives-ouvertes.fr/hal-01174084
DOI : 10.1016/j.gca.2015.02.005

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
In order to constrain the vertical and lateral extent of deformation and the interactions between lithosphere and asthenosphere in a context of a transpressional plate boundary, we performed teleseismic shear wave splitting measurements for 65 permanent and temporary broadband stations in central California. We present evidence for the presence of two anisotropic domains: (1) one with clear E–W trending fast directions and delay times in the range 1.5 to 2.0 s and (2) the other closely associated with the San Andreas Fault system with large azimuthal variations of the splitting parameters that can be modeled by two anisotropic layers. The upper of the two layers provides fast directions close to the strike of the main Californian faults and averaged delay times of 0.7 s; the lower layers show E–W directions and delay times in the range 1.5 to 2.5 s and thus can be compared to what is observed in stations that require a single layer. We propose the E–W trending anisotropic layer to be a 150 to 200 km thick asthenospheric layer explained by the shearing associated with the absolute plate motion of the North American lithosphere. The shallower anisotropic layer ought to be related to the dynamics of the San Andreas Fault system and thus characterized by a vertical foliation with lineation parallel to the strike of the faults localized in the lithosphere. We also propose that the anisotropic layer associated with each fault of the San Andreas Fault system is about 40 km wide at the base of the lithosphere.
ISSN: 0148-0227

hal-01236161
http://hal.univ-reunion.fr/hal-01236161
http://hal.univ-reunion.fr/hal-01236161/document
http://hal.univ-reunion.fr/hal-01236161/file/Bonnin_SKS_California_JGR2010.pdf
DOI : 10.1029/2009JB006438

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : Geochemistry, Geophysics, Geosystems, v. 9, p. Q12001, 2008. http://dx.doi.org/10.1029/2008GC002098
International audience
Rb/Sr, 40Ar/39Ar, and Sm/Nd isotopic data are reported in Caledonian eclogites from the Lindås Nappe, Bergen Arcs, Norway, in order to investigate processes controlling isotopic equilibrium at mineral scale in polymetamorphic rocks. The Bergen Arcs exposes Sveconorwegian ∼950–930 Ma granulites, partially overprinted by Caledonian eclogite-facies metamorphism at ∼425 Ma and amphibolite-facies metamorphism at ∼410 Ma. Geochemical and Rb-Sr data from more than 10 phengite fractions separated from one sample reflect the composition of the microdomain (a few hundred microns in size) in which phengite crystallized. Phengite crystallized after garnet or plagioclase by dissolution-precipitation processes yield apparent age between 700 and 600 Ma. At the time of their crystallization, these phengites inherited the isotopic composition of their precursor minerals, at a microdomain scale. Phengite from quartz veins, which crystallized from elements mobilized by the circulating fluid, yield an age closer to the eclogite-facies metamorphic age. The closed system evolution of the eclogitizing fluid, the segregation of textural and reactional microdomains, the high Sr content of the studied phengite, and the short duration of the recrystallization processes (<1Ma) are interpreted as the main factors responsible for the lack of a Rb/Sr isotopic equilibrium, at the scale of hand samples. Such equilibrium is nevertheless reached in quartz veins where the crystallization of minerals implies that the fluid circulation acted as a factor of isotopic homogenization. The in situ single mineral Ar/Ar datings revealed that both eclogite- and amphibolite-facies minerals are characterized by an excess of radiogenic argon (ages between 425 and 520 Ma). The excess of argon has been inherited from the previous granulite and has been only partially evacuated from the system by the circulating fluid. The apparent Sm/Nd ages from garnets inherited from the granulite-facies metamorphism (c. 930 Ma) are in agreement with previous estimates. This result confirms that the diffusion of Nd in garnet does not occur at temperatures lower than 700°C. This study highlights the complexity of radiochronometer behavior during HP metamorphism and demonstrates that coupling different radiochronometers, such as Ar/Ar and Rb/Sr, does not always guarantee the validity of the geochronological results.
ISSN: 1525-2027

hal-00407885
https://hal.archives-ouvertes.fr/hal-00407885
DOI : 10.1029/2008GC002098

Éditeur(s) : HAL CCSD Springer
Résumé : International audience
One of the most important topics that occupy public health problems is the air quality. That is the reason why mechanical ventilation and air handling units (AHU) were imposed by the different governments in the collective or individual buildings. Many buildings create an artificial climate using heating, ventilation and airconditioning (HVAC) systems. Among the existing aerosols in the indoor air, we can distinguish the bioaerosol with biological nature such as bacteria, viruses, fungi. Respiratory viral infections are a major public health issue because they are usually highly infective. We spend about 90% of our time in closed environments such as homes, workplaces or transport. Some studies have shown that AHU contribute to the spread and transport of viral particles within buildings. The aim of this work is to study the characterization of viral bioaerosols in indoor environments and to understand the fate of mengovirus eukaryote RNA virus on glass fiber filter F7 used in AHU. In this study, a setup close to reality of AHU system was used. The mengovirus aerosolized was characterized and measured with the Electrical Low Pressure Impact (ELPI) and the Scanner Mobility Particle Size (SMPS) and detected with RT-qPCR. The results about quantification and the level of infectivity of mengovirus on the filter and in the biosampler showed that mengovirus can pass through the filter and remain infectious upstream and downstream the system. Regarding the virus effectiveness on the filter under a constant air flow, and mengovirus was remained infectious during 10 hours after aerosolization.
ISSN: 1867-0334

hal-01636427
https://hal-univ-rennes1.archives-ouvertes.fr/hal-01636427
https://hal-univ-rennes1.archives-ouvertes.fr/hal-01636427/document
https://hal-univ-rennes1.archives-ouvertes.fr/hal-01636427/file/Fate%20of%20mengovirus%20on%20fiberglass%20filter_accepted.pdf
DOI : 10.1007/s12560-017-9310-8

Éditeur(s) : HAL CCSD Elsevier
Résumé : Different types of fractures and fault rocks occur within the major normal fault zones bordering the NW of the Neogene Penedes basin (NE Spain). In the southwest part of the basin, the main fault is the Baix Penedes fault that places the Upper Jurassic-Lower Cretaceous dolomicrites of the footwall in contact with the Miocene siliciclastic sediments of the hangingwall. Within the footwall damage zone, different fracture systems and their respective associated breccias and cements have been recognized. Using a combination of structural, petrologic and geochemical analyses, we show that these fracture systems evolved from mainly tectonic to mainly karstic features during the progressive development and upward propagation of the fault. The earlier fractures and associated breccias represent the initial stages of the fault propagation, when the process zone formed around the fault tip was dominated by dilatant fracturing (fractures 1 and dilatant breccias B1). The fluid involved in this early deformation precipitated the dolomite cement Cl, in equilibrium with the host-rock in a closed hydrogeological system. Localization of frictional processes generated cataclasites and ultracataclasites along distinct second-order faults (fractures 2). In these fractures, the dolomite cement C2 precipitated in relation to the incipient opening of the hydrological system to percolating meteoric waters. The later fractures and associated breccias represent the late stages of the fault propagation. These structures (fractures 3, fractures 4 and the dilatant breccias B2) were formed when the extensional deformation reached the surface, resulting in an increasingly open hydrogeological system. At these stages, deposition of the dolomitic sediment DS and precipitation of the calcite cements C3, C4 and C5 show the circulation of freshwaters in a vadose meteoric environment, i.e. in karstic conditions.
ISSN: 0375-6742

hal-00512690
https://hal.archives-ouvertes.fr/hal-00512690
DOI : 10.1016/j.gexplo.2009.11.010

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
Advances in groundwater storage monitoring are crucial for water resource management and hydrological processes understanding. The evaluation of water storage changes (WSC) often involve point measurements (observation wells, moisture probes, etc.), which may be inappropriate in heterogeneous media. Over the past few years, there has been an increasing interest in the use of gravimetry for hydrological studies. In the framework of the GHYRAF (Gravity and Hydrology inAfrica) project, 3 yr of repeated absolute gravity measurements using a FG5-type gravimeter have been undertaken at Nalohou, a Sudanian site in northern Benin. Hydrological data are collected within the long-term observing system AMMA-Catch. Once corrected for solid earth tides, ocean loading, air pressure effects, polar motion contribution and non-local hydrology, seasonal gravity variations reach up to 11 μGal, equivalent to a WSC of 260-mm thick infinite layer of water. Absolute temporal gravity data are compared to WSC deduced from neutron probe and water-table variations through a direct modelling approach. First, we use neutronic measurements available for the whole vertical profile where WSC occur (the vadose zone and a shallow unconfined aquifer). The RMSD between observed and modelled gravity variations is 1.61 μGal, which falls within the error bars of the absolute gravity data. Second, to acknowledge for the spatial variability of aquifer properties, we use a 2-D model for specific yield (Sy) derived from resistivity mapping and Magnetic Resonance Soundings (MRS). The latter provides a water content (θMRS) known to be higher than the specific yield. Hence, we scaled the 2-D model of θMRS with a single factor (α). WSC are calculated from water-table monitoring in the aquifer layer and neutronic measurements in the vadose layer. The value of α is obtained with aMonte-Carlo sampling approach, minimizing the RMSD between modelled and observed gravity variations. This leads to α = Sy/θMRS = 0.63 ± 0.15, close to what is found in the literature on the basis of pumping tests experiments, with a RMSD value of 0.94 μGal. This hydrogeophysical experiment is a first step towards the use of time-lapse gravity data as an integrative tool to monitor interannual WSC even in complicated subsurface distribution.
ISSN: 0956-540X

hal-00821181
https://hal.archives-ouvertes.fr/hal-00821181
DOI : 10.1093/gji/ggt146

Résumé : In the French West Indies islands, the transfer of chlordecone (kepone) from treated grounds towards the marine environment has been proved since the beginning of 2000s. In 2008, a maximal residue limit (MRL) was set at 20 µg / kg wet weight in fish and seafood, and therefore intensified the concern about the risks of high human exposure to the pesticide due to the consumption of contaminated marine products. To answer this concern, several sampling cruises were carried by public authorities and by Ifremer in order to analyze the contamination of the fish fauna around the Guadeloupe island. The data collected within these surveys improved our knowledge on the characteristics of the contamination of fish fauna and allowed to assess the distribution of this contamination in relation with optional geographical plans of the marine area around the Guadeloupe island. Further, this study was completed by an assessment of the fishery activity in areas which are the more exposed to the chlordecone contamination, taking into account data from a fishery observatory which has been set up in the archipelago since 2007. The results allow to precise two main characteristics of distribution of the contamination: its spatial and specific dimensions. About the spatial aspect, contamination of the fish fauna is centred on marine areas which are closely downstream contaminated watersheds. Concerning the marine fauna, two species groups are particularly exposed: the very coastal ones, particularly those living at or close to the mouth of rivers, and crustaceans wherever they live. Among fishes, the most contaminated species belong to many various groups which are present in the very coastal belt, more or less regularly. Lastly, this study provided orders of magnitude on the relative part of the fishery activity around the archipelagos for which the contamination is of great concern.
Aux Antilles, l’existence d’un transfert de la molécule de chlordécone des sols traités vers le milieu marin est avérée depuis le début des années 2000. L’établissement d’une limite maximale de résidus à 20 µg/kg pf pour tous les produits marins en 2008 y a intensifié la préoccupation vis-à-vis des risques d’exposition des populations humaines par les produits de la pêche. Pour répondre à cette préoccupation, depuis 2008, plusieurs campagnes d’échantillonnage ont été réalisées par les services de l’Etat (DAAF) et par l’Ifremer pour analyser la contamination de la faune halieutique autour de la Guadeloupe. Les données collectées de 2008 à 2011 par ces plans d’observation ont été utilisées pour actualiser les connaissances sur les caractéristiques de la contamination de la faune halieutique autour de la Guadeloupe, et analyser cette distribution en référence à différents schémas de segmentation du domaine marin. L’étude a été complétée par une analyse de l’activité halieutique dans les secteurs les plus sensibles à la contamination, à partir des données du système d’informations halieutiques déployé dans l’archipel depuis 2007. L’analyse a permis de préciser deux traits majeurs de cette distribution : sa dimension spatiale et son caractère spécifique. Sur le plan spatial, la contamination de la faune halieutique est centrée sur les espaces maritimes situées immédiatement en aval des bassins versants contenant des sols contaminés. Sur le plan faunistique, deux groupes apparaissent particulièrement sensibles, les espèces très côtières, en particulier celles fréquentant les embouchures des cours d’eau, et les crustacés, quelque soit leur milieu de vie. Chez les poissons, parmi les espèces contaminées, on trouve un large cortège d’espèces fréquentant plus ou moins régulièrement la bande littorale. Enfin, l’étude a permis de fournir des ordres de grandeur sur la part de l’activité halieutique de l’archipel directement concernée par cette pollution.

Droits : 2013 Ifremer
http://archimer.ifremer.fr/doc/00136/24762/22840.pdf
http://archimer.ifremer.fr/doc/00136/24762/

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
CK chondrites are highly oxidized meteorites containing abundant magnetite and trace amounts of Fe,Ni metal. Although the group is predominately composed of equilibrated meteorites (types 4-6), in recent years a significant number of new samples have been classified as being either CK3 or CK3-anomalous. These unequilibrated CKs often display a close affinity with members of the CV oxidized subgroup. CKs and CVs (oxidized subgroup) may therefore form a continuum and by implication could be derived from a single common parent body. To investigate the relationship between these two groups a detailed study of the oxygen isotope composition, opaque mineralogy and major and trace element geochemistry of a suite of CV and CK chondrites has been undertaken. The results of oxygen isotope analysis confirm the close affinity between CV and CK chondrites, while excluding the possibility of a linkage between the CO and CK groups. Magnetites in both CV and CK chondrites show significant compositional similarities, but high Ti contents are a diagnostic feature of the latter group. The results of major and trace element analysis demonstrate that both CV and CK chondrites show overlapping variation. Supporting evidence for a single common source for both groups comes from their similar cosmic-ray exposure age distributions. Recent reflectance spectral analysis is consistent with both the CVs and CKs being derived from Eos family asteroids, which are believed to have formed by the catastrophic disruption of a single large asteroid. Thus, a range of evidence appears to be consistent with CV and CK chondrites representing samples from a single thermally stratified parent body. In view of the close similarity between CV and CK chondrites some modification of the present classification scheme may be warranted, possibly involving integration of the two groups. One means of achieving this would be to reassigned CK chondrites to a subgroup of the oxidized CVs. It is recognized that a full evaluation of this proposal may require further study of the still poorly understood CK3 chondrites.
ISSN: 0016-7037

hal-00474633
https://hal.archives-ouvertes.fr/hal-00474633
DOI : 10.1016/j.gca.2009.11.029

Éditeur(s) : HAL CCSD Taylor & Francis
Résumé : The Latest Cretaceous/Paleocene paleokarsts from Bas-Languedoc (South of France) are characterized by : 1) exokarstic paleosurfaces (sink holes, pinnacles, canyons) which are mostly superimposed onto the Late Jurassic limestones and partly filled up with breccias and sandy clays containing Paleocene planktonic foraminifera; 2) endokarstic cavities filled up by sandy pelites and laminated mudstones with similar micropaleontological assemblage. All these cavities constitute three cut-and-filled paleokarstic systems. A general model of formation for this polyphase system is proposed. These paleosurfaces are generated by a Latest Cretaceous tectonic phase and Paleocene movements associated with eustatic rises inducing rapid floodings. During Upper Danian-Lower Selandian times three successive base levels drops, whose cumulated duration could be as long as 3 Ma, have induced several hundred metres of incisions within the exposed Jurassic carbonate series. A paleogeographic reconstruction shows three main marine gulfs which extend towards NE the EW Pyrenean thrust belt ("Pyrenean Paleocene Trough", PPT) located on the axis of the future Pyrenees. The similarity of the facies and micropaleontological content in the two domains suggests connexions between the marine Paleocene deposits of Bas-Languedoc and the PPT. It could thus explain the succession of the same tectono-eustatic events. We propose, as an hypothesis, that the rapid sea-level changes, as recorded by the imbricated karstic paleomorphologies, could be induced by the closing and the opening of a strait in the Eastern part of the PPT acting as a sill. This would be controlled by tectonic movements along the active orogenic axis of the Pyrenean Range and eustatic variations of the Paleocene World Ocean. The karstic systems developed during these low-stand episodes may have been later reactivated after the Early Selandian and more specifically during the Messinian desiccation event. This late evolution may have generated deep flooded karsts, now situated below the present sea level. Such karsts, partly inherited from Paleocene, could be important aquifers of economic interest.
ISSN: 0985-3111

hal-00407601
https://hal.archives-ouvertes.fr/hal-00407601
DOI : 10.3166/ga.20.301-326

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The full numerical simulation of turbulent flows in the context of industrial applications remains a very challenging problem. One of the difficulties is the presence of a large number of interacting scales of different orders of magnitude, ranging from the macroscopic scale to the Kolmogorov dissipation scale. In order to better understand and simulate these interactions, new algorithms of incremental type have recently been introduced, stemming from the theory of infinite-dimensional dynamical systems (see, e.g., the algorithms of Foias, Jolly et al. (1988), Foias, Manley and Temam (1988), Laminie et al. (1993, 1994), Marion and Temam (1989, 1990), and Temam (1990)). These algorithms are based on decompositions of the unknown functions into a large- and a small-scale component, one of the underlying ideas being to approximate the corresponding attractor. In the context of finite element methods, the decomposition of a solution into small- and large-scale components appears when we consider hierarchical bases (Yserentant, 1986; Zienkiewicz et al., 1982). In the present article we derive several estimates of the size of the structures for the linear and nonlinear terms which correspond to interactions of different hierarchical components of the velocity field, and also their time variation. The one-step time variation of these terms can be smaller than the expected accuracy of the computation. Using this remark, we implement an adaptive spatial and temporal multilevel algorithm which treats differently the small- and large-scale components of the flow. We derive several a priori estimates in order to study the perturbation introduced into the approximated equations. All the interaction terms between small- and large-scale structures are frozen during several time steps. Finally, we implement the multilevel method in order to simulate a bidimensional flow described by the Burgers equations. We perform a parametric study of the procedure and its efficiency. The gain in CPU time is significant and the trajectories computed by our multi-scale method remain close to the trajectories obtained with the classical Galerkin method.
ISSN: 0168-9274

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

Éditeur(s) : Université de Bretagne Occidentale
Résumé : Study of mineral suspended particulate matter (SPM) is essential to understand ecosystem dynamics. SPM influences light penetration which is closely related to primary production and carries pollutants and nutrients within the water masses. This work aims to study near-surface mineral SPM dynamics in the English Channel using two complementary methods: remote-sensing by satellites and numerical modelling. Firstly, MODIS and MERIS satellite images processed by the IFREMER semi-analytical algorithm and describing non-algal SPM, are analysed to build three statistical models which are able to estimate near-surface SPM using simple variables, e.g., tidal coefficient, significant wave height and mean chlorophyll-a concentration. Secondly, the three-dimensional hydrosedimentary model ROMS (Regional Ocean Modeling System) is applied in the English Channel. It considers tidal forcing along open boundaries, wind stress at the sea-surface, wave-current interactions in the bottom boundary layer and realistic heterogeneous bottom sediments. Results of the more sophisticated statistical model is also used to impose SPM concentrations at open boundaries. Sensitivity tests on SPM boundary conditions show the importance of prescribing realistic sediment concentrations to improve the model predictions. Biological influences on mineral SPM are also tested using settling velocity and critical shear stress for erosion which vary with seasons. Both statistical and numerical models highlight and locate respective influences of tide, waves and biology on mineral SPM. Numerical predictions indicate that near-surface SPM is mainly composed of silts (>70%). Near-surface suspended silts come largely from seabed resuspension alongshore in the Eastern Channel and from the open boundaries of our limited-area coastal model at the West of Cotentin Peninsula and offshore. Near-surface SPM variations are investigated at semi-diurnal, neap-spring tidal cycle time scales and during a wave event. Finally, dynamics of sediments in the turbidity maximum zone around the Isle of Wight is examined depending on hydrodynamical conditions in order to distinguish the role played by advection and local erosion/settling processes in this area. Keywords: Suspended particulate matter, satellite, remote-sensing, modelling, English Channel, sediments, turbidity, tide, wave
L'étude des matières en suspension (MES) minérales est essentielle pour comprendre le fonctionnement des écosystèmes en Manche : les MES influencent la pénétration de la lumière dans la colonne d'eau, paramètre clé pour la production biologique, et sont susceptibles de transporter des polluants et nutriments. Lors de ce travail, la dynamique des MES minérales de surface en Manche est étudiée à l'aide de deux méthodes complémentaires : l'observation par satellite et la modélisation numérique tridimensionnelle. Dans un premier temps, les images satellite MODIS et MERIS, traitées par l'algorithme semi-analytique développé par l'IFREMER et représentant les MES non-algales de surface, sont analysées afin d'élaborer trois modèles statistiques permettant d'estimer la MES de surface en fonction de variables basiques telles que le coefficient de marée, la hauteur significative des vagues et la concentration moyenne en chlorophylle-a. Dans un deuxième temps, le modèle hydrosédimentaire multiclasse tridimensionnel ROMS (Regional Ocean Modeling System) est implémenté en Manche. Il prend en considération le forçage par la marée aux frontières ouvertes, la contrainte exercée par le vent à la surface de la mer, l'interaction houle-courant en couche limite de fond et l’hétérogénéité spatiale réaliste du lit sédimentaire. Les résultats du modèle statistique le plus sophistiqué sont également utilisés pour forcer la concentration en MES aux frontières ouvertes du domaine de calcul. Des tests de sensibilité sur les conditions limites montrent l'importance de considérer de manière réaliste les sédiments entrant en Manche par les frontières ouvertes lors de la modélisation numérique. L'influence de la biologie sur les MES minérales est testée en utilisant une formulation de la vitesse de chute et de la contrainte critique d'érosion variant saisonnièrement. Grâce à ces deux outils associés, les influences respectives de la marée, de la houle et de la biologie sur les MES minérales sont mises en avant et localisées en Manche. Le modèle numérique prédit que les MES de surface sont composées principalement de silts (>70%). Les silts en suspension en surface proviennent majoritairement du lit sédimentaire près des côtes en Manche orientale et des frontières ouvertes du domaine de calcul à l'ouest de la presqu'île du Cotentin et au large en Manche centrale et orientale. Les variations des MES de surface sont ensuite étudiées aux échelles des cycles de marée semi-diurnes, vive-eau/morte-eau et autour d'un épisode de houles. Enfin, l'évolution de la concentration de sédiments dans la zone de forte turbidité autour de l'île de Wight peut être analysée suivant les conditions hydrodynamiques afin d'évaluer les rôles respectifs des processus locaux d'érosion/dépot et d'advection dans la génération de ces turbidités remarquables. Mots-clés : Matière en suspension, télédétection, satellite, modélisation, Manche, sédiments, turbidité, marée, houle

Droits : 2013 UBO
http://archimer.ifremer.fr/doc/00157/26783/24884.pdf
http://archimer.ifremer.fr/doc/00157/26783/

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The Roussillon sedimentary Basin (South France) is a complex multi-layered aquifer, close to the Mediterranean Sea facing seasonally increases of water abstraction and salinization issues. We report geochemical and isotopic vertical variability in this aquifer using groundwater sampled with a Westbay System® at two coastal monitoring sites: Barcarès and Canet. The Westbay sampling allows pointing out and explaining the variation of water quality along vertical profiles, both in productive layers and in the less permeable ones where most of the chemical processes are susceptible to take place.The aquifer layers are not equally impacted by salinization, with electrical conductivity ranging from 460 to 43,000 μS·cm− 1. The δ2H–δ18O signatures show mixing between seawater and freshwater components with long water residence time as evidenced by the lack of contribution from modern water using 3H, 14C and CFCs/SF6. S(SO4) isotopes also evidence seawater contribution but some signatures can be related to oxidation of pyrite and/or organically bounded S. In the upper layers 87Sr/86Sr ratios are close to that of seawater and then increase with depth, reflecting water–rock interaction with argillaceous formations while punctual low values reflect interaction with carbonate. Boron isotopes highlight secondary processes such as adsorption/desorption onto clays in addition to mixings. At the Barcarès site (120 m deep), the high salinity in some layers appear to be related neither to present day seawater intrusion, nor to Salses-Leucate lagoonwater intrusion. Groundwater chemical composition thus highlights binary mixing between fresh groundwater and inherited salty water together with cation exchange processes, water–rock interactions and, locally, sedimentary organic matter mineralisation probably enhanced by pyrite oxidation.Finally, combining the results of this study and those of Caballero and Ladouche (2015), we discuss the possible future evolution of this aquifer system under global change, as well as the potential management strategies needed to preserve quantitatively and qualitatively this water resource.
ISSN: 0048-9697

hal-01685219
https://hal-brgm.archives-ouvertes.fr/hal-01685219
DOI : 10.1016/j.scitotenv.2016.05.016

Éditeur(s) : Ieee-inst Electrical Electronics Engineers Inc
Résumé : The synthetic aperture radar (SAR) Doppler centroid has been used to estimate the scatter line-of-sight radar velocity. In weak to moderate ocean surface current environment, the SAR Doppler centroid is dominated by the directionality and strength of wave-induced ocean surface displacements. In this paper, we show how this sea state signature can be used to improve surface wind retrieval from SAR. Doppler shifts of C-band radar return signals from the ocean are thoroughly investigated by colocating wind measurements from the ASCAT scatterometer with Doppler centroid anomalies retrieved from Envisat ASAR. An empirical geophysical model function (CDOP) is derived, predicting Doppler shifts at both VV and HH polarization as function of wind speed, radar incidence angle, and wind direction with respect to radar look direction. This function is used into a Bayesian inversion scheme in combination with wind from a priori forecast model and the normalized radar cross section (NRCS). The benefit of Doppler for SAR wind retrieval is shown in complex meteorological situations such as atmospheric fronts or low pressure systems. Using in situ information, validation reveals that this method helps to improve the wind direction retrieval. Uncertainty of the calibration of Doppler shift from Envisat ASAR hampers the inversion scheme in cases where NRCS and model wind are accurate and in close agreement. The method is however very promising with respect of future SAR missions, in particular Sentinel-1, where the Doppler centroid anomaly will be more robustly retrieved.
Ieee Transactions On Geoscience And Remote Sensing (0196-2892) (Ieee-inst Electrical Electronics Engineers Inc), 2012-07 , Vol. 50 , N. 7 , P. 2901-2909

Droits : 2012 IEEE
http://archimer.ifremer.fr/doc/00088/19896/17554.pdf
DOI:10.1109/TGRS.2011.2174998
http://archimer.ifremer.fr/doc/00088/19896/

É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 Elsevier
Résumé : International audience
The Paleozoic-early Mesozoic geology of Tibet was controlled by the rift-drift, seafloor spreading and subduction zone tectonics of a Paleo-Tethyan realm, which evolved between the West Cathaysides (WC) and the East Cimmerides (EC). Different suture zones with ophiolites and ophiolitic mélanges, high-pressure metamorphic belts, magmatic arcs and accretionary prism complexes separating different terranes mark multiple subduction-accretion systems within this Paleo-Tethyan domain, reminiscent of the modern Western Pacific Ocean. Discrete basins separated by different continental blocks and magmatic arcs constituted a complex paleogeography of Paleo-Tethys, and these oceanic strands were closed as a result of subduction with different polarities during the late Paleozoic-Triassic. The Longmu Tso Shuanghu-Changning Menglian Suture zone (LS-CMS) in Tibet represents the main tectonic boundary between the WC and EC that developed in the Devonian. The East Kunlun-A’nyemaqen oceanic slab was subducted northward beneath the East Kunlun terrane in northern Tibet, whereas the Sumdo oceanic slab was subducted northward beneath the South Qiangtang-North Lhasa terrane in southern Tibet. The Longmu Tso Shuanghu-Changning Menglian and Jinshajiang-Ailaoshan-Song Ma ophiolites were developed and emplaced in subduction-accretion systems with opposite polarities (westward and eastward) beneath the North Qiangtang-Qamdo-Simao-Indochina terrane in central Tibet. The Jinshajiang-Ailaoshan-Song Ma oceanic slab was subducted westward beneath the North Qiangtang-Simao-Indochina terrane along the Jinshajiang and Ailaoshan-Song Ma trenches in a trench-ridge-trench triple junction plate configuration. The Emeishan mantle plume produced a large Permian basaltic terrane, developed on the western passive margin of the South China block. The final closure of the Paleo-Tethyan oceanic branches resulted in continental collisions and development of a vast Indosinian orogenic collage in the latest Triassic-Jurassic.
ISSN: 0743-9547

insu-01119604
https://hal-insu.archives-ouvertes.fr/insu-01119604
https://hal-insu.archives-ouvertes.fr/insu-01119604/document
https://hal-insu.archives-ouvertes.fr/insu-01119604/file/Xu-JAES-2015.pdf
DOI : 10.1016/j.jseaes.2015.01.021

Éditeur(s) : HAL CCSD Kluwer
Résumé : This paper presents some considerations about hybrid optimization algorithms useful to optimize real complex system. The given indication could be used by readers to conceive hybrid algorithms. These considerations have been deducted from a concrete application case: the satellites constellations design problem. But each of the advanced techniques proposed in this paper are considered in a more general way to solve other problems. This problem is used to illustrate the techniques along the paper because it is grouping many characteristics (di culties) of contemporary real complex systems: the size and the characteristics of the search space engendered by a combinatorial problem; The irregularity of the criterions; The mathematical and physical heterogeneity of parameters forbids the use of classical algorithms; The evaluation of a solution, which uses a time consuming simulation; A need of accurate values. More details are available in previous papers ([?], [?], [?], [?]). For these reasons, we could learn a lot from this experiment in order to detach hybrid techniques usable for problems having close characteristics. This paper presents the historic leading to the current algorithm, the modeling of the complex system and the sophisticated algorithm proposed to optimize it. Application cases and ways to built signi cant tests of hybrid algorithm are also given.
International Federation for Information Processing

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

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
In petroleum industry, the difference between pore pressure (Pp) and minimum horizontal stress Sh (termed the seal or retention capacity) is of major consideration because it is often assumed to represent how close a system is to hydraulic failure and thus the maximum hydrocarbon column height that can be maintained. While Sh and Pp are often considered to be independent parameters, several studies in the last decade have demonstrated that Sh and Pp are in fact coupled. However, the nature of this coupling relationship remains poorly understood. In this paper, we explore the influences of the spatial pore pressure distribution on Sh/Pp coupling and then on failure pressure predictions and trap integrity evaluation. With analytical models, we predict the fluid pressure sustainable within a reservoir before failure of its overpressured shale cover. We verify our analytical predictions with experiments involving analogue materials and fluids. We show that hydraulic fracturing and seal breach occur for fluid pressure greater than it would be expected from conventional retention capacity. This can be explained by the impact of the fluid overpressure field in the overburden and the pressure diffusion around the reservoir on the principal stresses. We calculate that supralithostatic pressure could locally be reached in overpressured covers. We also define the retention capacity of a cover (RC) surrounding a fluid source or reservoir as the difference between the failure pressure and the fluid overpressure prevailing in shale at the same depth. In response to a localized fluid pressure rise, we show that the retention capacity does not only depend on the pore fluid overpressure of the overburden but also on the tensile strength of the cover, its Poisson's ratio, and the depth and width of the fluid source.
ISSN: 0264-8172

hal-00590497
https://hal.archives-ouvertes.fr/hal-00590497
DOI : 10.1016/j.marpetgeo.2010.08.007