É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/

Résumé : Ce travail est l'expression d'une volonté de chercheurs, de la Carabe d'améliorer les connaissances scientifiques méso- et micro-météorologiques appliquées aux milieux insulaires sous influence des alizés et de développer la recherche dans ces domaines.On sait que le phénomène météorologique le plus remarquable impactant les îles de la Caraïbe reste le cyclone tropical. Mais d'autres phénomènes, a des échelles inferieures, telles que les pluies intenses, les houles, la dégradation de la qualité de l'air ont une importance aigüe en termes de risques naturels ou de risques sur la sante. Ces exemples attestent la nécessité d'utiliser des méthodes de descente d'échelle pour exploiter l'information météorologique et/ou climatique de grande échelle et dériver des scenarios locaux et régionaux appliques aux territoires insulaires. Ce défi est important car l'attente d'analyses scientifiques pertinentes dans ces domaines est grande.Les travaux présentés dans ce mémoire ont pour principaux objectifs la simulation numérique puis l'analyse des mécanismes méso- et micro-échelles qui induisent des circulations locales diurnes et nocturnes sur l'archipel de la Guadeloupe à l'aide de codes numériques météorologiques éprouves car largement utilises en recherche et en prévision opérationnelle.Ils constituent donc la première étude de modélisation numérique à haute résolution en basse atmosphère, par descente d'échelle dynamique, pour des intervalles d'espace compris entre 1 km et 111 m sur cet archipel.Le modelé atmosphérique Weather Research and Forecasting ARW V3.4 (WRF ARW), non-hydrostatique, a été utilise pour l'ensemble des simulations pour modéliser la troposphère depuis l'échelle globale à l'échelle de la turbulence.Trois situations météorologiques classiques réelles d'une durée de 48 heures, correspondant à 80 % des situations météorologiques observées dans la zone, alizes soutenus (AS), alizes moyens (AM), alizes faibles (AF) ont été examinées. Ces situations sont caractérisées par les valeurs du nombre de Froude local suivantes : 0,82 (AS), 0,41 (AM) et 0,21 (AF). Six domaines de modélisation ont été sélectionnés pour effectuer les descentes d'échelle dynamiques : D01 (maille de 27 km), D02 (maille de 9 km), D03 (maille de 3 km), D04 (maille de 1 km), D05 (maille de 333 m) et D06 (maille de 111 m) avec soixante-dix niveaux verticaux. Les quatre premiers domaines (D01 à D04) couvrent l'archipel de la Guadeloupe et sont utilisés en mode méso-échelle à l'aide d'un schéma d'ensemble de couche limite planétaire YSU. Les domaines D05 (couvrant l'île de la Basse-Terre et le centre de l'archipel) et D06 (couvrant la zone littorale et rurale du Petit Cul-de-Sac Marin et la zone urbaine de l'agglomération pointoise) sont utilisés en mode Large Eddy Simulation avec une fermeture de la turbulence 1,5 TKE 3D. Le modèle WRF a été forcé toutes les six heures par l'assimilation des champs d'analyses globales du modèle NCEP FNL (1° de résolution). Les simulations effectuées ont permis d'obtenir des champs de variables météorologiques 10-minutes à très hautes résolutions spatiales.Les résultats des simulations méso- et micro-échelles ont été confrontés aux valeurs expérimentales obtenues à l'aide de capteurs places sur des mâts météorologiques (campagne Atmo-Mangrov et réseau de mesure Météo-France). Il s'agissait d'optimiser l'utilisation des couplages de codes numériques tout en conservant la possibilité de les confronter aux observations expérimentales.Les résultats des simulations numériques micro-échelles, des différents cas (AS, AM, AF) sont utilisées pour forcer (c'est-à-dire définir les conditions limites) un modèle lagrangien de dispersion de particules : FLEXPART. Le système couple FLEXPART-WRF a été employé dans le domaine D06 pour étudier la dispersion du panache d'oxydes d'azote émis par la principale usine de production d'électricité de l'île.
The present work expresses the will of Caribbean researchers to improve the meso- and micro-meteorological scientific knowledge of the trade winds influenced island areas, and to develop research in these domains.It is well known that tropical hurricanes remain the most remarkable meteorological phenomena that affect the Caribbean islands. But some other phenomena, of smaller scale, such as intense rainfall events, swells, or air quality degradation, are of extreme importance for natural or health hazards. These examples show the need to use downscaling methods to exploit large scale meteorological or climatic information, and to obtain local and regional scenarios for the island areas. This is an important challenge, as sound scientific studies in these matters are eagerly expected. The aim of the research works exposed in the present dissertation is numerical simulation and analysis of the meso- and microscale mechanisms that induce diurnal and nocturnal local circulations in the Guadeloupe archipelago, using numerical meteorological models that are widely used in research and in operational forecasting.These works represent the first high-resolution (1 km to 111 m) numerical study of the lower atmosphere over the Guadeloupe archipelago. The Weather Research and Forecasting ARW 3.4 (WRF-ARW) model is used to simulate the troposphere from global scale to turbulence scale. Real cases of three typical weather types (80% of cases during a year) are examined during 48 hours : strong trade winds (STW), medium trade winds (MTW) and weak trade winds (WTW). These weather types are characterized by values of the local Froude number : 0.82 (STW), 0.41 (MTW) and 0.21 (WTW). Six domains have been selected for the dynamical downscaling : D01 (grid spacing of 27 km), D02 (grid spacing of 9 km), D03 (grid spacing of 3 km), D04 (grid spacing of 1 km), D05 (grid spacing of 333 m) and D06 (grid spacing of 111 m), including 70 vertical levels. The first four domains (D01 to D04) cover the Guadeloupe archipelago and are used in the meso-scale simulations with the planetary boundary layer scheme YSU (ensemble mean). Domain D05 (covering the Basse-Terre island and the middle of the archipelago) and domain D06 (covering the coastal and rural area of Le Petit Cul-de-Sac Marin and the urban area of Pointe-à-Pitre), are employed in the micro-scale simulation (LES) with the 3D TKE 1.5 order closure scheme. WRF has been 6 hourly reinitialized with the NCEP FNL global analyses (resolution of 1°). These simulations permitted to obtain 10-minutes meteorological variable fields with a very high resolution (111 m).Meso-scale and micro-scale model results have been evaluated with observational data from meteorological stations (field campaign Atmo-Mangrov, French Met Office).Once validated, the micro-scale model outputs have been used for the assimilation of the lagrangian particle dispersion model : FLEXPART. The coupling FLEXPART-WRF has been employed in domain D06, in order to analyze the dispersion of the nitrogen oxide plume emitted by the main power plant of the archipelago.

http://www.theses.fr/2014AGUY0735/document

Éditeur(s) : American Geophysical Union
Résumé : [1] Extreme weather events such as tropical cyclones are difficult to observe with conventional means. Satellite-based observations provide essential measurements of key parameters governing tropical cyclones. They are critical for short-term forecasting. Radiometers onboard the Defense Meteorological Satellite Program satellite series, WindSat and Tropical Rainfall Measuring Mission satellites, scatterometers onboard the ERS, ADEOS, and QuikScat satellites offer unprecedented synoptic observations of surface wind and atmospheric liquid water content, revealing the storm structures with good accuracy. However, satellite estimates do not provide direct measurements of geophysical parameters and can suffer from limitations linked to the sensors characteristics, such as the signal wavelength and polarization or the measurement incidence angle. For example, measurements at Ku band are strongly affected by rain. Still, each observing system can offer specific information that can be combined with the others. In particular, we highlight the capabilities of dual-frequency altimeter to provide very high resolution measurements of rain rate, surface wind speed, and wave characteristics. A method is proposed to obtain continuous along-track 5 km resolution measurements of these parameters in the tropical cyclone Isabel. The results shows that dual-frequency altimeters can provide useful information to complement and validate the operational fields provided by the atmospheric numerical models and by NOAA observing systems.
Journal of Geophysical Union - Research C - Oceans (0148-0227) (American Geophysical Union), 2006 , Vol. 111 , N. C1 , P. NIL_38-NIL_50

Droits : 2006 American Geophysical Union
http://archimer.ifremer.fr/doc/2006/publication-1033.pdf
DOI:10.1029/2005JC003068
http://archimer.ifremer.fr/doc/00000/1033/

É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 Elsevier
Résumé : International audience
Tropical islands are characterized by thermal and orographical forcings which may generate microscale air mass circulations. The Lesser Antilles Arc includes small tropical islands (width lower than 50 km) where a total of one-and-a-half million people live. Air quality over this region is affected by anthropogenic and volcanic emissions, or saharan dust. To reduce risks for the population health, the atmospheric dispersion of emitted pollutants must be predicted. In this study, the dispersion of anthropogenic nitrogen oxides (NOx) is numerically modelled over the densely populated area of the Guadeloupe archipelago under weak trade winds, during a typical case of severe pollution. The main goal is to analyze how microscale resolutions affect air pollution in a small tropical island. Three resolutions of domain grid are selected: 1 km, 333 m and 111 m. The Weather Research and Forecasting model (WRF) is used to produce real nested microscale meteorological fields. Then the weather outputs initialize the Lagrangian Particle Dispersion Model (FLEXPART). The forward simulations of a power plant plume showed good ability to reproduce nocturnal peaks recorded by an urban air quality station. The increase in resolution resulted in an improvement of model sensitivity. The nesting to subkilometer grids helped to reduce an overestimation bias mainly because the LES domains better simulate the turbulent motions governing nocturnal flows. For peaks observed at two air quality stations, the backward sensitivity outputs identified realistic sources of NOx in the area. The increase in resolution produced a sharper inverse plume with a more accurate source area. This study showed the first application of the FLEXPART-WRF model to microscale resolutions. Overall, the coupling model WRF-LES-FLEXPART is useful to simulate the pollutant dispersion during a real case of calm wind regime over a complex terrain area. The forward and backward simulation results showed clearly that the subkilometer resolution of 333 m is necessary to reproduce realistic air pollution patterns in this case of short-range transport over a complex terrain area. Globally, this work contributes to enrich the sparsely documented domain of real nested microscale air pollution modelling. This study dealing with the determination of the proper resolution grid and proper turbulence scheme, is of significant interest to the near-source and complex terrain air quality research community.
ISSN: 1352-2310

hal-01349208
https://hal.archives-ouvertes.fr/hal-01349208
DOI : 10.1016/j.atmosenv.2016.05.028

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
Knowing the slip amplitudes that large earthquakes produced in prehistorical times is one key to anticipate the magnitude of large forthcoming events. It is long known that the morphology is preserving remnants of paleoearthquake slips in the form of fault-offset landforms. However, the measured offsets that can be attributed to the most recent paleoearthquakes are generally few along a fault, so that they rarely allow recovering the slip distributions and largest slips of these earthquakes. We acquired ~1 m resolution airborne lidar data on a 30 km stretch of a fast-slipping strike-slip fault (eastern Hope Fault, New Zealand) located in a region of high alluvial dynamics where landforms are rapidly evolving. Data analysis reveals >200 offset landforms; only 30% allow a good to moderate quality offset measurement. From these good to moderate quality measures, we recover the slip-length distributions and largest slips of the four most recent large paleoearthquakes and find evidence of 4–6 prior events. The record suggests that large earthquake slip recurred in multiples of about 4 m along the 30 km stretch. Although they have larger uncertainties, the more numerous lower-quality offsets that we also measured reveal a similar earthquake slip record. This shows that, although offset landforms are partly degraded in dynamically active landscapes, they store valuable information on paleoearthquake slips. This information might be recovered provided that the morphology is analyzed at high resolution and " continuously " over a significant fault length. Remote lidar data are powerful to perform such analyses.
ISSN: 2169-9313

insu-01469219
https://hal-insu.archives-ouvertes.fr/insu-01469219
https://hal-insu.archives-ouvertes.fr/insu-01469219/document
https://hal-insu.archives-ouvertes.fr/insu-01469219/file/Manighetti_et_al-2015-Journalsur Geophysical_Research__Solid_Earth.pdf
DOI : 10.1002/2014JB011787

Éditeur(s) : HAL CCSD Ocean Drilling Programm
Résumé : This paper provides a summary of postcruise scientific results from Ocean Drilling Program (ODP) Leg 209 available to date, building upon shipboard observations and syntheses summarized in the Leg 209 Initial Results volume. During Leg 209, 19 holes were drilled at 8 sites along the Mid-Atlantic Ridge from 14°43´ to 15°44´N, mainly in residual mantle peridotite intruded by gabbroic rocks, in order to understand the tectonic and structural processes responsible for formation of oceanic lithosphere with abundant residual peridotite exposed on the seafloor coupled with a relatively low proportion of volcanic rocks. Based on proportions of recovered lithologies, the entire area may be underlain by mantle peridotite with ~20%–40% gabbroic intrusions and impregnations. Impregnated peridotites with olivine + two pyroxenes + plagioclase + spinel that apparently formed in equilibrium probably record crystallization from primitive mid-ocean-ridge basalt at pressures of 0.5–0.6 GPa. Metamorphic equilibria record isobaric cooling to ~1100°C at this pressure. Thus, the conductively cooled thermal boundary layer beneath the Mid-Atlantic Ridge in this region is >15 km thick. Combined crystallization and reaction with residual peridotite formed a series of impregnated peridotites recording increasing Na content at nearly constant Mg#; this process could explain some of the variation in fractionation-corrected Na (e.g., Na = 8.0) observed in mid-ocean-ridge basalts. Clinopyroxene textures and compositions record such impregnation processes, and they are particularly well documented for Site 1274. Other Leg 209 gabbroic rocks formed from extensive crystallization of highly evolved melts, indicating that a substantial proportion of melt entering the thermal boundary layer crystallizes entirely beneath the seafloor, with no volcanic equivalent. Alteration of peridotites occurred over a range of temperatures and is the result of three distinct processes: rock-dominated serpentinization with formation of brucite in olivine-rich lithologies, fluid-dominated serpentinization with formation of magnetite and no brucite, and fluid-dominated talc alteration with addition of SiO2 as well as H2O and oxygen. The latter two processes also exhibit detectable trace element metasomatism that is distinct in its character from the igneous impregnation described in the previous paragraph. Microstructures show that most residual peridotites were not ductilely deformed at temperatures less than ~1200°C. Structural and paleomagnetic data require tectonic rotations of relatively undeformed blocks; some rotations probably exceeded 60° around nearly horizontal axes parallel to the rift axis. Rotations occurred along several generations of high-temperature mylonitic shear zones extending deeper than 15 km depth and numerous faults at lower temperature. Early formed shear zones and faults were passively rotated around later features; such a process could have produced low-angle fault surfaces without slip on low-angle faults. This region provides end-member examples of processes that are common at many or most slow-spreading ridges. Osmium isotope ratios indicate an ancient history of depletion for residual peridotites from the 14°–16°N region along the Mid-Atlantic Ridge. Though depleted Os isotope ratios in peridotite have been reported elsewhere along the global ridge system, the values from this region are among the most depleted. In general, Os isotope ratios from mid-ocean-ridge basalts are systematically more radiogenic than Os isotope ratios from ridge peridotite samples, suggesting a polygenetic heterogeneous source for mid-ocean-ridge basalts. Geochemical studies of zircons from Leg 209 gabbroic rocks and impregnated peridotites, together with other ridge and arc-related zircons, indicate that ridge zircons have systematically lower fractionation-corrected U and Th concentrations compared to arc zircons. This observation provides a tool for interpreting the tectonic provenance of ancient detrital zircons and indicates an arclike provenance for Hadean detrital zircons. Geobiological studies and aerobiological studies were also undertaken during Leg 209. The geobiological work found no measurable microbial enhancement of olivine dissolution rate, possibly because the samples from Leg 209 were sterile. The aerobiological study determined that dust from North Africa, collected from the derrick of the JOIDES Resolution during Leg 209, contains a variety of abundant microorganisms.
Proceedings ODP, scientific results, 209

hal-00407966
https://hal.archives-ouvertes.fr/hal-00407966
DOI : 10.2973/odp.proc.sr.209.001.2007

Éditeur(s) : IEEE
Résumé : This paper aims at studying the quality of the sea surface salinity (SSS) retrieved from soil moisture and ocean salinity (SMOS) data in coastal areas. These areas are characterized by strong and variable SSS gradients [several practical salinity units (psu) on relatively small scales: the extent of river plumes is highly variable, typically at kilometric and daily scales. Monitoring this variability from SMOS measurements is particularly challenging because of their resolution (typically 30-100 km) and because of the contamination by the nearby land. A set of academic tests was conducted with a linear coastline and constant geophysical parameters, and more realistic tests were conducted over the Bay of Biscay. The bias of the retrieved SSS has been analyzed, as well as the root mean square (rms) of the bias, and the retrieved SSS compared to a numerical hydrodynamic model in the semirealistic case. The academic study showed that the Blackman apodization window provides the best compromise in terms of magnitude and fluctuations of the bias of the retrieved SSS. Whatever the type of vegetation cover, a strong negative bias, greater than 1 psu, was found when nearer than 36 km from the coast. Between 44 and 80 km, the type of vegetation cover has an impact of less than a factor 2 on the bias, and no influence further than 80 km from the coast. The semirealistic study conducted in the Bay of Biscay showed a bias over ten days lower than 0.2 psu for distances greater than 47 km, due to an averaging over various geometries (coastline orientation, swath orientation, etc.). The bias showed a weak dependence on the location of the grid point within the swath. Despite the noise on the retrieved SSS, contrasts due to the plume of the Loire River and the Gironde estuary remained detectable on ten-day averaged maps with an rms of 0.57 psu. Finally, imposing thresholds on the major axis of the measurements brought little improvement to the bias, whereas it increased the rms and- could lead to strong swath restriction: a 49-km threshold on the major axis resulted in an effective swath of 800-900 km instead of 1200 km. NOT CONTROLLED OCR
Transactions on geoscience an remote sensing IEEE (0196-2892) (IEEE), 2007-07 , Vol. 45 , N. 7 , P. 2061-2072

Droits : 2007 IEEE
http://archimer.ifremer.fr/doc/2007/publication-3643.pdf
DOI:10.1109/TGRS.2007.894934
http://archimer.ifremer.fr/doc/00000/3643/

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The Gulf of Cadiz spans the plate boundary between Africa and Eurasia west of the Betic–Rif mountain belt. A narrow east dipping subduction zone descends beneath the Gulf of Cadiz and the straits of Gibraltar. The deep crustal structure of the Gulf and the adjacent SW Iberian and Moroccan margins is constrained by numerous multi-channel seismic reflection and wide-angle seismic surveys. A compilation of these existing studies is presented in the form of depth to basement, sediment thickness, depth to Moho and crustal thickness maps. These structural maps image an E–W trending trough, with thin (< 10 km) crust beneath the Gulf of Cadiz. This trough is filled by an eastward thickening wedge of sediments, reaching a thickness of 10–15 km in the eastern Gulf. These sediments are tectonically deformed, primarily along a series of westward-vergent thrust faults and represent a 200–250 km wide accretionary wedge. The northern and especially the southern limits of the accretionary wedge are marked by sharp morphological lineaments showing evidence of recent deformation. These tectonic limits are situated in an internal position with respect to the Miocene deformation front (external Betic and Rif allocthons), which has been abandoned. At the western boundary of the accretionary wedge, near the adjacent Seine and Horseshoe abyssal plains, an E–W trending basement high (Coral Patch Ridge) can be seen indenting the deformation front in an asymmetric manner. Analog modeling is performed using granular materials accreted against a semicircular backstop (representing the basement of the Rif and Betic mountain belts). The modeling initially produces a symmetric, arcuate accretionary wedge. The ensuing collision of an oblique rigid indenter retards accretion on one side, resulting in an embayment and a locally steeper deformation front. The deformation pattern observed in morphology and high-resolution seismic profiles suggests the accretionary wedge and underlying subduction system is still active. The implications of active subduction for the source region of the 1755 Lisbon earthquake and the regional seismic hazard assessment are discussed.
ISSN: 0040-1951

hal-00424922
http://hal.univ-brest.fr/hal-00424922
http://hal.univ-brest.fr/hal-00424922/document
http://hal.univ-brest.fr/hal-00424922/file/Gutscher_09Tphys_rev.pdf
DOI : 10.1016/j.tecto.2008.11.031

Éditeur(s) : HAL CCSD Lambert Academic Publishing
Résumé : The Himalaya-Tibet system is considered as the classical case of continental collision. Nonetheless, some fundamental questions concerning its structure, as well as physical processes influencing its evolution are still pending. This book presents several advances in the subject matter. Based on new high-resolution geophysical data from the field, combined with robust imaging techniques, the structure of the collision is revealed in detail: the underplating of India beneath half of Tibet, as well as the geometry of the surface along which this happens. Various numerical modelling techniques allow to discuss physical processes at work in the collision zone, such as the rheology of the India plate whose strong mantle supports the Himalayas, or the metamorphic reactions at depth emphasizing the role of water. We recommend this book to Earth Scientists in a broader sense. While specialists in seismology or modelling can thumb technical details, all Readers can get an interdisciplinary overview of mountain building in the Himalayas, with a retrospective outlook on early day explorers and the dawn of some classical geological principles.
https://hal.archives-ouvertes.fr/hal-00854592

ISBN : 978-3838323534
hal-00854592
https://hal.archives-ouvertes.fr/hal-00854592

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

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

Éditeur(s) : HAL CCSD Geological Society of London
Résumé : We have applied transmission electron microscopy (TEM) analyses coupled with viscoplastic self-consistent (VPSC) numerical modelling to identify the active slip systems and to better understand the crystal preferred orientation (CPO) development of the Torridon quartz mylonite (NW Scotland). TEM analyses showed evidence of activation of 1/3〈a〉{π′}, 1/3〈a〉{z} and possible 〈a〉(c) slip systems, as well as dislocation climb and dynamic recrystallization. All the CPOs generated by VPSC models share common characteristics with the Torridon quartz mylonite, but only Models 2 and 3 reproduce the [c]-axes maxima at low angle (<20°) to the foliation pole along the YZ plane, as observed in the mylonite. In Model 2, this concentration only occurs at γ≥2.6, whereas in Model 3 this maxima occurs at lower shear strains. The models that start with a previous preferred orientation acquire very strong CPOs after small-imposed strains, followed by the rapid rotation of the fabric in relation to the new imposed finite strain axes. The combined activation of 〈a〉{π′}, 〈a〉{z} and possibly 〈a〉(c) slip systems, as demonstrated by TEM analyses, suggests that the VPSC model that best predicts CPO development in the Torridon quartz mylonite is Model 2, where the critical resolved shear stress (CRSS) of 〈a〉{π/π′} is assumed to be slightly stronger than 〈a〉(c).
Deformation Mechanisms, Rheology and Tectonics: Microstructures, Mechanics and Anisotropy

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

Éditeur(s) : Springer
Résumé : Series of sensitivity tests were performed with a z-coordinate, global eddy-permitting (1/4 degrees) ocean/sea-ice model (the ORCA-R025 model configuration developed for the DRAKKAR project) to carefully evaluate the impact of recent state-of-the-art numerical schemes on model solutions. The combination of an energy-enstrophy conserving (EEN) scheme for momentum advection with a partial step (PS) representation of the bottom topography yields significant improvements in the mean circulation. Well known biases in the representation of western boundary currents, such as in the Atlantic the detachment of the Gulf Stream, the path of the North Atlantic Current, the location of the Confluence, and the strength of the Zapiola Eddy in the south Atlantic, are partly corrected. Similar improvements are found in the Pacific, Indian, and Southern Oceans, and characteristics of the mean flow are generally much closer to observations. Comparisons with other state-of-the-art models show that the ORCA-R025 configuration generally performs better at similar resolution. In addition, the model solution is often comparable to solutions obtained at 1/6 or 1/10 degrees resolution in some aspects concerning mean flow patterns and distribution of eddy kinetic energy. Although the reasons for these improvements are not analyzed in detail in this paper, evidence is shown that the combination of EEN with PS reduces numerical noise near the bottom, which is likely to affect current-topography interactions in a systematic way. We conclude that significant corrections of the mean biases presently seen in general circulation model solutions at eddy-permitting resolution can still be expected from the development of numerical methods, which represent an alternative to increasing resolution.
Ocean Dynamics (1616-7341) (Springer), 2006-12 , Vol. 56 , N. 5-6 , P. 543-567

Droits : © Springer 2006. Part of Springer Science+Business Media
http://archimer.ifremer.fr/doc/2006/publication-3514.pdf
DOI:10.1007/s10236-006-0082-1
http://archimer.ifremer.fr/doc/00000/3514/

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
We propose a new thermomechanical numerical model of mid-ocean ridge accretion aimed at investigating asymmetric spreading rates, diverse configurations of lens and sill magma injections, crystallization and depth, and on- and off-axis patterns of hydrothermal cooling. The numerical algorithm iteratively resolves temperature and motion equations until it reaches a stationary solution. The motion equation was written in a vorticity-stream function formalism, with boundary and internal conditions applied to the stream function to impose the style of magma injection. Unlike in previous models, our model does not assume an a priori shape for the temperature field, which is initiated by an initial half-space cooling according to the left and right spreading rates. Complex patterns of hydrothermal cooling are simulated by enhanced thermal diffusivity. The model succeeds in describing the dynamic and thermal effects of spreading rates, the style of magma intrusion, and the hydrothermal cooling. Accurate descriptions of these are essential to study the cooling histories of crustal rocks and geophysical observables.
ISSN: 1525-2027

hal-00413118
https://hal.archives-ouvertes.fr/hal-00413118
DOI : 10.1029/2008GC002270

Éditeur(s) : HAL CCSD
Résumé : Among erosion processes, river incision is classicaly described as a key process controlling erosion of landscapes. Incision efficiency is mainly influenced by climate and erodibility. This latter is not only dependent on the nature of the bedrock, but also on its past deformation, which affects its rheological effective properties, such as fracture density. The main objectives of this thesis are: (1) to better constrain the relationship between effective properties and erodibility, and (2) to quantify the influence of erodibility and erosion on both the temporal and spatial building or decay of the topography. Several numerical tools are developed. A 1D formalism of landscape evolution is introduced, including river incision with stochastic distribution of water discharge and hillslope landsliding. A new remeshing algorithm called Surface Lagrangian Remeshing (SLR) is developed as a complement to remeshing algorithms dealing with internal elements. It allows one to take into account long-term erosion into 2D Lagrangian numerical codes based on triangular finite elements. Then the potentiality of measuring erodibility in-situ using a Schmidt hammer (R) is assessed for the active orogen of Taiwan, the diagenetic Annot sandstones and St Clement fault zone. Results suggest a strong control of R by effective properties. A linear model based on effective medium theory is applied to a fault zone with an unmacthed resolution (750 measures, 25 measures per square meter). The model successfully correlates R to fracture density. These results demonstrate that effective elasticity as well as erodibility are sensitive to the density and type of fractures. Finally I focus on the erosional and rheological conditions that allows reproducing post-orogenic evolution of mountain belts. A model coupling surface erosion and regional isostatic uplift is consistent with observations. The topographic decay and decrease of the ratio of surface elevation over crustal root thickness is at first order controlled by the initial geometry of the mountain belt and erosion efficiency. This new model highlights the control ofclimate and erodibility on the topographic decay and of lithospheric rheology on the perseverance of crustal roots.
Parmi les processus d'érosion, l'incision des rivières est classiquement décrit comme un processus clé contrôlant l'érosion des paysages. L'efficacité de l'incision est principalement influencée par le climat et par l'érodabilité. Ce dernier paramètre ne dépend pas seulement de la nature du substratum rocheux, mais aussi de sa déformation passée qui affecte ses propriétés rhéologiques équivalentes. Les principaux objectifs de cette thèse sont: (1) de mieux contraindre les relations entre propriétés équivalentes et érodabilité, et (2) de quantifier l'influence de l'érosion et de l'érodabilité sur la formation ou la décadence, spatiale et temporelle, de la topographie. Plusieurs outils numériques sont développés. Un formalisme 1D d'évolution des paysages est proposé, prenant en compte simultanément l'incision des rivières avec une distribution stochastique des débits en eau et l'érosion des versants par glissements de terrain. Un nouvel algorithme de remaillage appelé Surface Lagrangian Remeshing (SLR) est développé. Il permet de prendre en compte l'érosion à long terme dans les codes numériques 2D Lagrangiens basés sur des éléments finis triangulaires. Ensuite, la possibilité de mesurer in-situ l'érodabilité avec un marteau de Schmidt (R) est évaluée pour: l'orogène actif de Taiwan, les grès diagénétiques d'Annot et la zone de faille de St Clément. Les résultats suggèrent un fort contrôle de R par les propriétés équivalentes des roches testées. Un modèle linéaire, basé sur la théorie des milieux équivalents, est appliqué à une zone de faille avec une résolution inégalée (750 mesures, 25 mesures par mètre carré). Le modèle permet de corréler avec succès la densité de fracture et R. Ces résultats démontrent empiriquement que l'élasticité équivalente ainsi que l'érodabilité sont sensibles à la densité et à la nature des fractures. Enfin j'étudie les conditions d'érosion et de rhéologie qui permettent de reproduire l'évolution des chaînes de montagnes post-orogéniques. Un modèle couplant érosion en surface et soulèvement isostatique régional est compatible avec les observations. Les taux de décroissance topographique et de diminution du rapport de l'élévation de surface sur l'épaisseur de racine crustale sont contrôlés au premier ordre par la géométrie initiale de la chaîne de montagne et par l'efficacité de l'érosion. Ce nouveau modèle met en évidence le contrôle du climat et de l'érodabilité sur la décroissance topographique et de la rhéologie lithosphérique sur la persévérance des racines crustales.
https://tel.archives-ouvertes.fr/tel-00602341

tel-00602341
https://tel.archives-ouvertes.fr/tel-00602341
https://tel.archives-ouvertes.fr/tel-00602341/document
https://tel.archives-ouvertes.fr/tel-00602341/file/thesephilippesteer.pdf

Résumé : La thèse est consacrée à la recherche des catastrophes marines (tsunamis, ondes de tempête) dans les Antilles françaises, en utilisant des modèles analytiques et numériques de mécanique des fluides. L'accent est mis sur le développement de modèles de mouvement de glissements de terrain et des vagues causées par ces glissements. Le modèle le plus simple du glissement de terrain «solide block» est utilisé pour évaluer les caractéristiques des flux pyroclastiques du volcan Soufrière-Hills (Montserrat). Un modèle plus complexe de modélisation de glissement de terrain (modèle de Savage-Hutter) a été étudié analytiquement, donnant ainsi une nouvelle famille de solutions exactes décrivant le mouvement de l'écoulement par gravité non linéaire dans une vallée ou dans un canyon sous-marin. Le modèle comprend des ondes simples (Riemann waves),le cas d'un barrage qui cède (dam-break problem), des solutions auto-similaires dans la M-fonne et « chapeaux paraboliques ». Grâce à la théorie linéaire de l'eau peu profonde nous avons étudié le processus de génération de tsunamis par des glissements de terrain de volume variable, se déplaçant à une vitesse variable dans un bassin de profondeur variable. Dans le cas d'un fond marin particulier (cas sans "réflexion"), les phénomènes de résonnance ont été étudiés dans un bassin à profondeur variable. Nous avons utilisé des méthodes numériques pour la résolution non-linéaire des équations des eaux peu profondes afin d'analyser des catastrophes marines réelles (les ondes de tempête causées par le cyclone Lili en 2002, le tsunami volcanique de 2003 à Montserrat) et probables (un tsunami prés des côtes de la Martinique). Des données sur les catastrophes
The dissertation is devoted to research in the field of marine natural hazards (tsunamis, storm surges) in the French West Indies, using analytical and numerical models of fluid mechanics. Emphasis is placed on the development of models of landslide motion and generated tsunami waves. The simple "solid block" model is used to evaluate the characteristics of pyroclastic flow Soufriere Hills volcano (Montserrat). The "fluid model" of a landslide (so called Savage-Hutter model) is studied analytically; within this model a new family of exact solutions that describe the motion of nonlinear gravity flow in a valley or underwater canyon is found: nonlinear Riemann wave, dan break problem, self-similar solutions (M - wave and parabolic cap). In the framework of the linear shallow water theory the process of generation of tsunami waves by landslides of variable volume moving with variable velocity above the basin of variable depth is studied. For the specific bottom profile ("reflectionless" one) the resonant phenomena is investigated in the basin of variable depth. Numerical methods are used to analyze marine hazards: historical (storm surges, caused by Cyclone LILI in 2002; volcanic tsunami 2003 on Montserrat) and possible events (possible tsunami of the coast of Martinique). Various data on marine natural disasters are obtained during field surveys (volcanic tsunami in 2003, stonn surges caused by Hurricane Dean in 2007). Designed catalogs of tsunamis and storm surges are created based on results of numerical modeling and field studies; some statistical analysis is perfomed

http://www.theses.fr/2011AGUY0435/document

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
The South Pacific region is characterized by a broadly elevated seafloor known as the South Pacific superswell. This region has a concentration of midplate volcanoes that experienced massive eruptions in the mid-Cretaceous period (90-120 Ma). These characteristics suggest the presence of a large-scale mantle plume beneath the South Pacific, called the South Pacific superplume. The geometry, origin depth, temperature, and composition of the superplume remain controversial, however, mainly due to the lack of seismological data that documents the mantle structure beneath the South Pacific. Seismic stations are sparse in the area due to its remote ocean environment. To obtain a better seismic image of the superplume, we deployed temporary broadband seismographs on oceanic islands and the seafloor in the South Pacific, which made possible the highest spatial resolution that has ever been achieved for the mantle structure beneath the region. The seismic image obtained from this new seismic data indicates that large-scale low-velocity anomalies (on the order of 1000 km in diameter), indicative of the superplume, are located from the bottom of the mantle to a depth of 1000 km, and small-scale low-velocity anomalies (on the order of 100 km in diameter) are present above it. A comparison of the seismic image with recent mantle convection studies based upon laboratory and numerical experiments suggests that the superplume may be a hot and chemically distinct mantle dome, and that the small-scale anomalies may be narrow plumes generated from the top of the dome. This model may explain various characteristics of hot spots in the South Pacific, such as the seafloor swell, short-lived hot spot chains, and the periodicity of massive eruptions.
ISSN: 1525-2027

hal-00443845
https://hal.archives-ouvertes.fr/hal-00443845
https://hal.archives-ouvertes.fr/hal-00443845/document
https://hal.archives-ouvertes.fr/hal-00443845/file/Suetsugu_PLUME_tomo_G32009_hal.pdf
DOI : 10.1029/2009GC002533

Éditeur(s) : HAL CCSD Geological Society of America
Résumé : International audience
Lakes are major depositional systems for which the related depositional processes have long been considered relatively simple. Breaking this statement, this study presents a detailed analysis of deposits in Lake Saint-Jean, the third largest natural lake in Québec. In addition to postglacial deltaic and coastal depositional systems fringing the lake, current-controlled features such as a large subaqueous prograding wedge and three sediment drifts have been identifi ed in its central portion based on two-dimensional (2-D) acoustic high-resolution subbottom profiles. The large subaqueous prograding wedge is a 4-km-long and up to 15-m-thick heterolithic shelf-like construction in the southeastern part of the lake. The three sediment drifts are 0.1–0.5-km-long and 2–5-mthick mud mounds distributed on the lake floor in the central portion of the lake. Diatom analyses and radiocarbon dating show that the development of these current-controlled features occurred during the lacustrine phase, after the disconnection with the postglacial marine Laflamme Gulf at 8.5 cal. k.y. B.P. Depositional facies show evidence of recurrent bottom-current activity. Related deposits alternate with pelagic sedimentation stages characterized by the settling of mud and biogenic accumulations. We investigated the origin of bottom currents using a numerical simulation (SYMPHONIE, an oceanographic model), with the aim of modeling wind-induced lake-scale water circulation. Simulations suggest that the subaqueous prograding wedge and the three sediment drifts result from wind-induced bottom currents generated by storm events having wind speed greater than 10 m s–1. Such strong winds are able to significantly affect sedimentation in the central portion of Lake Saint-Jean. Theresulting wind-induced sedimentary features were integrated into a refi ned lacustrine depositional model that summarizes the evolution of a group of water bodies referred to as “wind-driven water bodies.” This study applies a new tool for lake strata characterization and highlights the potential diffi cultyin differentiating them from marine deposits in the geological record.
ISSN: 0016-7606

insu-01184724
https://hal-insu.archives-ouvertes.fr/insu-01184724
DOI : 10.1130/B31145.1

Éditeur(s) : HAL CCSD Springer-Verlag
Résumé : International audience
High-resolution borehole strainmeters are usually installed in tectonically active regions in order to detect slow-slip events, and to estimate slow transients related to earthquake swarms. However, they are also sensitive to other numerous influences, internal or external. Furthermore, the quality of their coupling to the rock through cementation, and the mechanical properties of the rock mass around them, have a critical influence on their records. Many of the existing strainmeters present such problems, and the correction for these effects often remains a challenge. In this paper, we present the analysis of the records of a high-resolution borehole dilatometer (Sacks–Evertson), located in the seismically active rift of Corinth (Greece) (station TRZ in the Trizonia island). We show that the instrument suffers from an only partial-solid coupling, and that the nearby sea tides have a direct (through elastic response) and indirect (through pore-pressure diffusion) effect on the dilatation signal, which adds up to the solid tidal strain source. We propose a methodology that allows, in a first step, to better separate the internal (solid tide) from the external (air pressure, sea level) influences, by calculating a frequency-dependent transfer function outside the range of the tidal periods. We then extrapolate this function, in particular at the tidal periods. In a second step, the resulting variation with frequency of the coupling coefficients with sea level led us to estimate the proportion of instrument not solidly cemented to rock (thus in contact with water pore pressure), which is about 90 % of the total height. Despite the small proportion of solid coupling, the sensor resolution remains very good up to a few tens of hours of a time period, thanks to the confining effects of the rocks on the local pore pressure. These results allow us to correct for the external effects, and reduce the associated variance by 80–90 % (in the period range of minutes to days). The empirical correction of the sea level effect could be explained using a simple Boussinesq’s approximation and 1D pore-pressure diffusion model, which contributed to better constraint of some of the poro-elastic parameters in the vicinity of the instrument. After correction, the solid tidal signal at the 24-h period is almost anti-correlated with those of the theoretical solid tide. This surprising result is consistent with a similar anti-correlation observed for the longest period surface waves (200 s) comparing the TRZ dilatometer signals to the strain measured by a nearby borehole strainmeter (MOK, 15 km). This could be related to the presence of a shallow fault close to the instrument, which would creep in response to seismic wave-related stress.
ISSN: 0033-4533

hal-01122853
https://hal.archives-ouvertes.fr/hal-01122853
DOI : 10.1007/s00024-013-0742-2

Éditeur(s) : Amer Geophysical Union
Résumé : Bolide impacts and large-scale volcanic eruptions have been proposed as possible causes of the massive extinction of life that has occurred episodically in Earth's history. In spite of the catastrophic disruption of the local environment that accompanies bolide impacts and volcanic eruptions, it has been difficult to explain why these events sometimes lead to global extinction of species. We propose that, in some cases, the missing link may be provided by hypercanes, runaway hurricanes that are capable of injecting massive amounts of water and aerosols into the middle and upper stratosphere, where they may have profound effects on atmospheric chemistry and radiative transfer. Hypercanes are theorized to occur when the sea surface temperature exceeds a critical threshold, which may occur when sea water is locally heated by bolide impact, shallow-sea volcanism, or possibly, by overturning of superheated brine pools formed by underwater volcanic activity. Simulations using a convection-resolving nonhydrostatic, axisymmetric numerical model show that hypercanes can indeed develop when the sea surface temperature is high, and that they inject large amounts of mass into the stratosphere.
Journal Of Geophysical Research-atmospheres (2169-897X) (Amer Geophysical Union), 1995-07 , Vol. 100 , N. D7 , P. 13755-13765

Droits : 1995 by the American Geophysical Union.
http://archimer.ifremer.fr/doc/00258/36966/35611.pdf
DOI:10.1029/95JD01368
http://archimer.ifremer.fr/doc/00258/36966/