Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
The Cretaceous paleogeography and the kinematic evolution of the Iberian plate are poorly constrained. Especially problematic is to reconcile Iberian paleomagnetic data with paleomagnetic data of the neighboring plates and with Euler poles derived from seafloor magnetic anomalies. The first limitation arises from the Cretaceous Normal Polarity Superchron where paleogeographic reconstruction using marine magnetic anomalies is handicapped. The second arises from the paucity of reliable paleomagnetic poles with satisfactory statistical criteria and age. In order to address these shortcomings and provide new high quality paleomagnetic poles for Iberia, we conducted a detailed rock magnetic and paleomagnetic study of two Cretaceous magmatic sills, the Paço de Ilhas (PI) and Foz da Fonte (FF) sills, from the Lusitanian Basin, Portugal, recently dated at about 88 and 94 Ma, respectively. Our results show that the magnetic mineralogy of the sills is primary, i.e., acquired during magma cooling, and essentially represented by titanomagnetite. The corresponding paleomagnetic poles match the synthetic APWP from the African plate at 80 and 100 Ma. On the basis of a rigorous selection of Iberian Cretaceous poles, we then calculated mean paleomagnetic poles for different time intervals and found that Iberian paleomagnetic data fit well the global APWP between 70 and 120 Ma, but move far away from the APWP at pre-rift times. Our approach shows that new and better constrained paleomagnetic poles can aide in solving part of the contradiction between Iberian and African APWPs.
ISSN: 0148-0227

hal-00758103
https://hal.archives-ouvertes.fr/hal-00758103
DOI : 10.1029/2011JB009067

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The origin of plate tectonics is one of the most fundamental issues in earth and planetary sciences. Laboratory experiments indicate that the viscosity of silicate rocks is so strongly temperature-dependent that the entire surface of the Earth should be one immobile rigid plate. The rheology of oceanic lithosphere is, however, still poorly understood, and there exist few constraints on the temperature dependency of viscosity on the field scale. Here we report a new kind of observational constraint based on the geoid along oceanic fracture zones. We identify a large number of conspicuous small-scale geoid anomalies, which cannot be explained by the standard evolution model of oceanic lithosphere, and estimate their source density perturbations using a new Bayesian inversion method. Our results suggest that they are caused most likely by small-scale convection involving temperature perturbations of ∼300 K±100 K∼300 K±100 K. Such thermal contrast requires the activation energy of mantle viscosity to be as low as View the MathML source100±50 kJmol−1 in case of diffusion creep, and View the MathML source225±112 kJmol−1 in case of dislocation creep, substantially reducing the thickness of the stiffest part of oceanic lithosphere. Oceanic lithosphere may thus be broken and bent much more easily than previously thought, facilitating the operation of plate tectonics.
ISSN: 0012-821X

hal-01467827
https://hal.archives-ouvertes.fr/hal-01467827
DOI : 10.1016/j.epsl.2016.05.027

É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 Oxford University Press (OUP)
Résumé : International audience
Several studies, both theoretical and experimental, show that sedimentary rocks have a fractal pore-grain interface. In this paper a computer simulated 3-D sedimentary rock structure generated by the Relaxed Ballistic Bidisperse Deposition Model (RBBDM), is investigated to characterize the micro structure of its pores. The pore volume and the rock-pore interface show the same fractal dimension indicating that the pore volume is a fractal. The two point density correlation is computed for the pore space and the results compare favourably with the range reported from experiments. An array of 2-D X-ray tomography micrograph sections of a real sedimentary rock, an oolitic limestone (pure calcite) from the Mondeville formation of Middle Jurassic age (Paris Basin, France), was used to generate a 3-D bitmap. The 3-D real rock sample generated in this manner, was analysed for similar studies as the simulated structure. The results were compared with those obtained from simulation. The simulation results agree qualitatively with the real rock sample. Diffusion through the connected pore space of the simulated structure was studied using a random walk algorithm and the results compared with the similar simulation study done on the 3-D oolitic limestone specimen. In both cases diffusion was found to be anomalous indicating that the sedimentary rock has a fractal geometry. The favourable comparability of results between the simulated and real rock supports the usefulness of the model of sedimentary rock generation which can be applicable to transport phenomena.
ISSN: 0956-540X

hal-00823061
https://hal.archives-ouvertes.fr/hal-00823061
DOI : 10.1093/gji/ggs084

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : Gravity data and P-wave teleseismic traveltime residuals from 29 temporary broad-band stations spread over the northern margin of the Gulf of Aden (Dhofar region, Oman) were used to image lithospheric structure. We apply a linear relationship between density and velocity to provide consistent density and velocity models from mid-crust down to about 250 km depth. The accuracy of the resulting models is investigated through a series of synthetic tests. The analysis of our resulting models shows: (1) crustal heterogeneities that match the main geological features at the surface; (2) the gravity edge effect and disparity in anomaly depth locations for layers at 20 and 50 km; (3) two low-velocity anomalies along the continuation of Socotra-Hadbeen and Alula-Fartak fracture zones between 60 and 200 km depth; and (4) evidence for partial melting (3-6 per cent) within these two negative anomalies. We discuss the presence of partial melting in terms of interaction between the Sheba ridge melts and its along-axis segmentation.
ISSN: 0956-540X

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

Éditeur(s) : HAL CCSD
Résumé : Joint inversions are now commonly used in Earth Sciences. They have been developed to better understand Earth structure since they provide more constraints on the inverted parameters. We propose a new process to simultaneously invert several data sets in order to better image 3-D crustal and upper mantle structures. Our inversion uses 3 kind of data which present a good complementarity : (1) P receiver functions to provide Moho depth variations, (2) teleseismic delay times of P waves to retrieve velocity anomalies in the crust and upper mantle, (3) gravity anomalies to image density variations at lithospheric scale. In our scheme, receiver functions are first inverted. The resulting Moho depths are then interpolated and incorporated as a priori information into the joint inversion of teleseismic delay times and gravity anomalies process. In our new approach, we perform a model space search for Moho variations, P-velocity and density structure to find acceptable fit to the three data sets. In order to preferentially sample the good data fit region, we chose the Neighborhood algorithm of Sambridge to optimistically survey the model space. We model the delay times with a 3-D raytracing using velocity nodes evenly spaced and linked with density nodes via a linear relationship. This relationship can evolve through depth, simulating its pressure dependance. We test our scheme on synthetic examples, and apply the method in a continental rifting region where deformation processes are complex and badly known. This joint inversion thus provides more constraints on the inverted parameters, and allows us to better image crust - mantle interactions in the case of the Hangai dome (Central Mongolia). Indeed, in this area, we expect an asthenospheric upwelling to dynamically support the topography, but we image a density - velocity anomaly in the lower crust. The effects on topography and the origin (composition vs. thermal) of this anomaly are still under debate.
Les inversions conjointes sont maintenant communément utilisées en Sciences de la Terre. Elles ont été développées afin d'améliorer notre connaissance et compréhension de la structure interne de la Terre puisqu'elles apportent de plus en plus de contraintes sur les paramètres inversés. Dans cette thèse, nous proposons un nouveau processus d'inversion conjointe qui prend en compte trois paramètres différents et qui mène à l'obtention d'un modèle lithosphérique tridimensionnel. Notre méthode utilise trois types de données gravimétriques et sismologiques qui présentent une bonne complémentarité : (1) Les fonctions récepteur P afin d'obtenir les variations de profondeur du Moho, (2) les délais de temps d'arrivées P des téléséismes pour retrouver les anomalies de vitesse dans la croûte et le manteau supérieur, (3) les anomalies gravimétriques qui donnent accès aux variations de densité à l'échelle lithosphérique. Dans le schéma d'inversion proposé, nous inversons d'abord les fonctions récepteur. Les variations de profondeur du Moho qui en résultent sont alors interpolées puis incorporées comme information a priori dans le processus d'inversion conjointe gravimétrie - tomographie télésismique. Que ce soit pour les fonctions récepteur ou pour la partie conjointe, notre méthode d'inversion se base sur un algorithme de recherche afin de trouver la structure qui minimise l'écart aux trois types de données. Nous avons opté pour l'algorithme de voisinage qui présente l'avantage de concentrer sa recherche dans les régions qui minimisent l'écart aux données sans toutefois abandonner les autres régions afin de limiter les risques de convergence vers un minimum local. Nous utilisons un modèle 3-D constitué de noeuds de vitesse afin de modéliser les délais de temps. Ces noeuds de vitesse sont liés aux noeuds de densité grâce à une relation linéaire entre vitesse et densité. Cette relation peut évoluer avec la profondeur de manière à simuler les effets dus aux variations de pression. Nous avons réalisé plusieurs tests synthétiques afin d'évaluer le comportement de notre nouvelle méthode. Ils ont permis de montrer que les variations du Moho ainsi que les anomalies de vitesse et densité sont bien retrouvées. Ils ont également montré l'importance de la prise en compte de ces trois types de données dans un même schéma d'inversion, notamment en ce qui concerne la faculté de discerner les variations d'interfaces des anomalies volumiques. Nous avons ensuite appliqué la méthode dans une région continentale en phase de rifting où les processus de déformation sont complexes et mal connus. Notre nouvelle méthode a apporté davantage de contraintes sur les paramètres inversés et a permis d'obtenir un modèle correctement résolu de la croûte et du manteau supérieur sous le dôme de Hangai (Mongolie Centrale). Une anomalie asthénosphérique sous le dôme est souvent invoquée pour expliquer sa topographie élevée mais c'est une anomalie (en vitesse et densité) crustale (croûte inférieure) qui s'exprime dans nos modèles. Ses effets sur la topographie et son origine (compositionnelle et/ou thermique) sont encore débattus et nécessiteraient des études plus poussées (paramétriques, VP =VS, ...).
https://tel.archives-ouvertes.fr/tel-00783524

tel-00783524
https://tel.archives-ouvertes.fr/tel-00783524
https://tel.archives-ouvertes.fr/tel-00783524/document
https://tel.archives-ouvertes.fr/tel-00783524/file/Basuyau2011.pdf

Éditeur(s) : HAL CCSD American Geophysical Union - AGU
Résumé : IODP Expedition 335 "Superfast Spreading Rate Crust 4" returned to ODP Hole 1256D with the intent of deepening this reference penetration of intact ocean crust several hundred meters into cumulate gabbros. This was the fourth cruise of the superfast campaign to understand the formation of oceanic crust accreted at fast spreading ridges, by exploiting the inverse relationship between spreading rate and the depth to low velocity zones seismically imaged at active mid-ocean zones, thought to be magma chambers. Site 1256 is located on 15-million-year-old crust formed at the East Pacific Rise during an episode of superfast ocean spreading (>200 mm/yr full rate). Three earlier cruises to Hole 1256D have drilled through the sediments, lavas and dikes and 100 m into a complex dike-gabbro transition zone. The specific objectives of IODP Expedition 335 were to: (1) test models of magmatic accretion at fast spreading ocean ridges; (2) quantify the vigor of hydrothermal cooling of the lower crust; (3) establish the geological meaning of the seismic Layer 2-3 boundary at Site 1256; and (4) estimate the contribution of lower crustal gabbros to marine magnetic anomalies. It was anticipated that even a shortened IODP Expedition could deepen Hole 1256D a significant distance (300 m) into cumulate gabbros. Operations on IODP Expedition 335 proved challenging from the outset with almost three weeks spent re-opening and securing unstable sections of the Hole. When coring commenced, the destruction of a hard-formation C9 rotary coring bit at the bottom of the hole required further remedial operations to remove junk and huge volumes of accumulated drill cuttings. Hole-cleaning operations using junk baskets returned large samples of a contact-metamorphic aureole between the sheeted dikes and a major heat source below. These large (up to 3.5 kg) irregular samples preserve magmatic, hydrothermal and structural relationships hitherto unseen because of the narrow diameter of drill core and previous poor core recovery. Including the ~60 m-thick zone of granoblastic dikes overlying the uppermost gabbro, the dike-gabbro transition zone at Site 1256 is over 170 m thick, of which more than 100 m are recrystallized granoblastic basalts. This zone records a dynamically evolving thermal boundary layer between the principally hydrothermal domain of the upper crust and a deeper zone of intrusive magmatism. The recovered samples document a sequence of evolving geological conditions and the intimate coupling between temporally and spatially intercalated intrusive, hydrothermal, contact-metamorphic, partial melting and retrogressive processes. Despite the operational challenges, we achieved a minor depth advance to 1522 m, but this was insufficient penetration to complete any of the primary objectives. However, Hole 1256D has been thoroughly cleared of junk and drill cuttings that have hampered operations during this and previous Expeditions. At the end of Expedition 335, we briefly resumed coring and stabilized problematic intervals with cement. Hole 1256D is open to its full depth and ready for further deepening in the near future.
InterRidge News

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

Éditeur(s) : HAL CCSD Springer Verlag
Résumé : International audience
Polarized emittance measurements were acquired for synthetic forsterite, the pure magnesium end member of the olivines group, on the whole infrared spectral range and up to the melting point by using CO2 laser heating. The experimental data, fitted with a semi-quantum dielectric function model, allowed the retrieval of the temperature dependence of the absorption coefficient of forsterite both in the opaque and semi-transparent regions. The analysis of the phonon parameters indicates that the lattice dynamics evolve drastically with increasing temperature. The normal modes involving motions of the magnesium cations located in site 1 are the more impacted, and some of them vanish around 1,200 K. The results confirm that the enhancement of the lattice anharmonicity and the increasing mobility of the magnesium cations are closely linked and are at the origin of the anomalies observed in the evolution of the thermophysical properties. This complete set of spectroscopic data may be a step toward a more precise evaluation of the impact of thermal radiation heat transfer inside systems involving forsterite and quantification of their heat budget.
ISSN: 0342-1791

hal-00817079
https://hal.archives-ouvertes.fr/hal-00817079
DOI : 10.1007/s00269-013-0570-z

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
Estimating how topography is maintained provides insights into the different factors responsible for surface deformations and their relative roles. Here, we develop a new and simple approach to assess the degree of isostatic compensation of continental topography at regional scale from GOCE gravity gradients. We calculate the ratio between the radial gradient observed by GOCE and that calculated from topography only. From analytical and statistical formulations, simple relationships between this ratio and the degree of compensation are obtained under the Airy–Heiskanen isostasy hypothesis. Then, a value of degree of compensation at each point of study area can be easily deduced. We apply our method to the Alaska-Canada Cordillera and validate our results by comparison with a standard isostatic gravity anomaly model and additional geophysical information for this area. Both our GOCE-based results and the isostatic anomaly show that Airy–Heiskanen isostasy prevails for the Yukon Plateau whereas additional mechanisms are required to support topography below the Northwest Territories Craton and the Yakutat collision zone.
ISSN: 0956-540X

hal-01467540
https://hal.archives-ouvertes.fr/hal-01467540
https://hal.archives-ouvertes.fr/hal-01467540/document
https://hal.archives-ouvertes.fr/hal-01467540/file/cadioGJI2016.pdf
DOI : 10.1093/gji/ggw281

Éditeur(s) : HAL CCSD
Résumé : Guadeloupe Island (West French Indies) is one of the twenty islands that compose the Lesser Antilles Arc, which results from the subduction of the Atlantic Ocean plate beneath the Caribbean one. The island lies in a complex volcano-tectonic system and the need to understand its geological context has led to numerous on- and offshore geophysical investigations. This work presents the compilation and processing of available, on-land, airborne and marine, gravity and magnetic data acquired during the last 40 years on Guadeloupe Islands and at the scale of the Lesser Antilles Arc. The overall dataset provides new Bouguer and reduced to the pole magnetic anomaly maps at the highest achievable resolution. Regionally, the main central negative gravity trend of the arc allows defining two subsident areas. The first one is parallel to the arc direction (~N160°E) to the north, whereas the second unexpected southern one is oriented parallel to oceanic ridges (N130°E). Along the Outer Arc, the long wavelength positive anomaly is interpreted, at least along the Karukera Spur, as an up-rise of the volcanic basement in agreement with the seismic studies. To the NE of Guadeloupe, the detailed analysis of the geophysical anomalies outlines a series of structural discontinuities consistent with the main bathymetric morphologies, and in continuity of the main fault systems already reported in this area. Based on geophysical evidences, this large scale deformation and faulting of the Outer Arc presumably primarily affects the Atlantic subducting plate and secondarily deforms the upper Caribbean plate and the accretion prism. At the scale of Guadeloupe Island, joined gravity and magnetic modeling has been initiated based on existing interpretation of old seismic refraction profiles, with a general structure in three main layers. According to our geophysical anomalies, additional local structures are also modeled in agreement with geological observations: i) the gravity and magnetic signals confirm an up-rise of the volcanic basement below the limestone platforms outcropping on Grande-Terre Island ; ii) the ancient volcanic complexes of Basse-Terre Island are modeled with high density and reverse magnetized formations; iii) the recent volcanic centre is associated with formations consistent with the low measured density and the underlying hydrothermal system. The EW models coherently image a NNW-SSE depression structure in half-graben beneath Basse-Terre Island, its western scarp following the arc direction in agreement with bathymetric and seismic studies to the north of the island. The so-defined depressed area, and particularly its opening in half-graben toward the SW, is interpreted as the present-day front of deformation of the upper plate associated with the recent volcanic activity on and around Guadeloupe. Based on this regional deformation model, perspectives are given for further integrated investigation of key targets to address the internal structure and evolution of the Lesser Antilles Arc and Guadeloupe volcanic system.
Colloque 2012 du Comité National Français de Géodésie et Géophysique (CNFGG)
Clermont-Ferrand, France

hal-00742653
https://hal-brgm.archives-ouvertes.fr/hal-00742653
https://hal-brgm.archives-ouvertes.fr/hal-00742653/document
https://hal-brgm.archives-ouvertes.fr/hal-00742653/file/CNFGG_Caraibes.pdf

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
Both flow field heterogeneity and mass transfer between mobile and immobile domains have been studied separately for explaining observed anomalous transport. Here we investigate non-Fickian transport using high-resolution 3-D X-ray microtomographic images of Berea sandstone containing microporous cement with pore size below the setup resolution. Transport is computed for a set of representative elementary volumes and results from advection and diffusion in the resolved macroporosity (mobile domain) and diffusion in the microporous phase (immobile domain) where the effective diffusion coefficient is calculated from the measured local porosity using a phenomenological model that includes a porosity threshold ( inline image) below which diffusion is null and the exponent n that characterizes tortuosity-porosity power-law relationship. We show that both flow field heterogeneity and microporosity trigger anomalous transport. Breakthrough curve (BTC) tailing is positively correlated to microporosity volume and mobile-immobile interface area. The sensitivity analysis showed that the BTC tailing increases with the value of inline image, due to the increase of the diffusion path tortuosity until the volume of the microporosity becomes negligible. Furthermore, increasing the value of n leads to an increase in the standard deviation of the distribution of effective diffusion coefficients, which in turn results in an increase of the BTC tailing. Finally, we propose a continuous time random walk upscaled model where the transition time is the sum of independently distributed random variables characterized by specific distributions. It allows modeling a 1-D equivalent macroscopic transport honoring both the control of the flow field heterogeneity and the multirate mass transfer between mobile and immobile domains.
ISSN: 0043-1397

hal-01278285
https://hal.archives-ouvertes.fr/hal-01278285
DOI : 10.1002/2015WR017645

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
The nature and extent of the regional lithosphere-asthenosphere interaction beneath Ireland and Britain remains unclear. Although it has been established that ancient Caledonian signatures pervade the lithosphere, tertiary structure related to the Iceland plume has been inferred to dominate the asthenosphere. To address this apparent contradiction in the literature, we image the 3-D lithospheric and deeper upper-mantle structure beneath Ireland via non-linear, iterative joint teleseismic-gravity inversion using data from the ISLE (Irish Seismic Lithospheric Experiment), ISUME (Irish Seismic Upper Mantle Experiment) and GRACE (Gravity Recovery and Climate Experiment) experiments. The inversion combines teleseismic relative arrival time residuals with the GRACE long wavelength satellite derived gravity anomaly by assuming a depth-dependent quasilinear velocity-density relationship. We argue that anomalies imaged at lithospheric depths probably reflect compositional contrasts, either due to terrane accretion associated with Iapetus Ocean closure, frozen decompressional melt that was generated by plate stretching during the opening of the north Atlantic Ocean, frozen Iceland plume related magmatic intrusions, or a combination thereof. The continuation of the anomalous structure across the lithosphere-asthenosphere boundary is interpreted as possibly reflecting sub-lithospheric small-scale convection initiated by the lithospheric compositional contrasts. Our hypothesis thus reconciles the disparity which exists between lithospheric and asthenospheric structure beneath this region of the north Atlantic rifted margin.
ISSN: 0956-540X

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

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
Our knowledge and understanding of the 3-D lithospheric structure of the Himalayas and the Tibetan Plateau is still challenging although numerous geophysical studies have been performed in the region. The GOCE satellite mission has the ambitious goal of mapping Earth's gravity field with unprecedented precision (i.e. an accuracy of 1-2 mGal for a spatial resolution of 100 km) to observe the lithosphere and upper-mantle structure. Consequently, it gives new insights in the lithospheric structure beneath the Himalayas and the Tibetan Plateau. Indeed, the GOCE gravity data now allow us to develop a new strategy for joint gravimetry-seismology inversion. Combined with teleseismic data over a large region in a joint inversion scheme, they will lead to lithospheric velocity-density models constrained in two complementary ways. We apply this joint inversion scheme to the Hi-CLIMB (Himalayan-Tibetan Continental Lithosphere during Mountain Building) seismological network which was deployed in South Tibet and the Himalayas for a 3-yr period. The large size of the network, the high quality of the seismological data and the new GOCE gravity data set allow us to image the entire lithosphere of this active area in an innovative way. We image 3-D low velocity and density structures in the middle crust that fit the location of discontinuous low S-velocity zones revealed by receiver functions in previous geophysical studies. In the deeper parts of our velocity model we image a positive anomaly interpreted to be the heterogenous Indian lithosphere vertically descending beneath the centre of the Tibetan Plateau.
ISSN: 0956-540X

hal-00858024
https://hal.archives-ouvertes.fr/hal-00858024
DOI : 10.1093/gji/ggs110

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
3D seismic data located in the Gjallar Ridge (Vøring Basin, offshore Norway) reveals a closely-spaced polygonal fault system affecting more than 800 m of homogeneous mud-dominated Quaternary and Tertiary sequences. As some faults reach the modern seafloor, they represent an active polygonal fault system at present day. Even if the processes remain unclear and are still under debate, it is generally agreed that the initiation of polygonal faults is the result of shallow burial dewatering of fine-grained unconsolidated sediments by volumetric compaction. 3D seismic data are commonly interpreted by propagating horizons automatically and by picking faults manually. However, in the case of polygonal fault intervals, this approach is time consuming due to the huge number of faults and because automatic propagation can be misleading. In this study, we applied a new technique of 3D seismic interpretation based on a sequential stratigraphy analysis, using the new PaleoScan© software (Eliis Company). It allowed us to build a 3D geological model computing more than 300 horizons within the faulted intervals. We then used the coherency attribute, depicting anomalies in the shape of seismic waveform like faults, in order to constrain a possible link between fault distribution and stratigraphic levels. Our approach allows fault throws to be calculated in milliseconds on any polygonal fault plane. The result shows that fault segments have been reactivated by dip-linkage. Distribution of faults depends on mechanical units, intervals characterized by different petrophysical properties, which are independent from lithological and diagenetic changes. According to these results, we propose a model showing the evolution of polygonal fault intervals in which faulting stages are separated by a quiescence phase during burial. A first tier of polygonal faults is initiated at a specific depth, according to the Cam-clay model. Then, following a period of quiescence during which mud-rich sediments continued to accumulate, new fault segments are initiated above the first mechanical unit and within this undeformed interval. New nucleated faults then connect downward to pre-existing underlying polygonal fault system, thus progressively increasing the thickness of the faulted interval.
ISSN: 0025-3227

hal-00772470
https://hal.archives-ouvertes.fr/hal-00772470
DOI : 10.1016/j.margeo.2012.07.016

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Several studies, both theoretical and experimental, show that sedimentary rocks have fractal pore-grain interface. The real rocks are 3-d structures with highly tortuous and often fractal pore spaces. Before attempting simulation of this daunting geometry, we present here, as a preliminary study, a simpler 2-d version. In this paper, a computer simulated 2-d sedimentary rock structure is generated by the relaxed bidisperse ballistic deposition model. Grains of two different sizes are dropped ballistically on a linear substrate. By changing the fraction of the two types of particles, the porosity of the rock structure can be tuned. The structure undergoes compaction through the relaxation of possible unstable overhangs. The micro structure of the pore space is investigated. The pore mass and the rock-pore interface show a fractal behaviour with the same fractal dimension indicating that the pore volume is a fractal. Our simulation results indicate that the process of compaction of grains during the deposition process seems to erase the dependency of the fractal dimension on the grain size distribution. The two point density correlation is measured for the pore space. It shows anisotropy which is an outcome of the growth rule. X-ray tomography of two-dimensional sections of real sedimentary rocks obtained from Mallorca Island is subjected to the same study and the results compared with those obtained from simulation. The simulation results agree qualitatively with the real rock sample. We also study diffusion on the pore space. Diffusion is found to be anomalous as is expected in fractal spaces. It also bears the signature of anisotropy of the structure. Diffusion studies on the real rock sample could not yield conclusive results as the system size is not large enough.
ISSN: 0926-9851

hal-00761818
https://hal.archives-ouvertes.fr/hal-00761818
DOI : 10.1016/j.jappgeo.2012.09.002

Éditeur(s) : HAL CCSD Royal Meteorological Society
Résumé : Assessing the impact/adaptation of human activities on/to climate change is a key issue, especially in the tropics that concentrate major anthropogenic dynamics such as deforestation and nearly two-thirds of the planetary rainfall. However, this task is often made tough because human activities such as agricultural dynamics are usually analysed at local or regional scale whereas climate related studies are led at large to global scales due to a lack of reliable data, especially in the tropics. In this article we argue that the increased spatial resolution of remote sensing-based rainfall estimates enables assessing the spatiotemporal variability of rainfall regimes at regional and local scales, thus allowing fine analysis of the interactions with human activities. We processed Tropical Rainfall Measuring Mission (TRMM) 3B42 daily rainfall estimates over the state of Mato Grosso (southern Brazilian Amazon) for the 1998-2012 study period in order to compute rainfall metrics such as annual rainfall and duration, onset and end dates of the rainy season based on the Anomalous Accumulation methodology (at a 0.25◦ spatial resolution). We then crossed these metrics with agricultural maps (produced at a 250m spatial resolution) and proved that the adoption of intensive agricultural practices such as double cropping systems is partly the result of a strategy to adapt practices to local climatic conditions. Finally, we discuss how such results raise important issues regarding the sustainability of the agricultural development model in the Southern Amazon.
International Journal of Climatology

halshs-00910825
https://halshs.archives-ouvertes.fr/halshs-00910825
DOI : 10.1002/joc.3863

Éditeur(s) : HAL CCSD Nature Publishing Group
Résumé : International audience
Temperature and fluid pressure conditions control rock deformation and mineralization on geological faults, and hence the distribution of earthquakes1. Typical intraplate continental crust has hydrostatic fluid pressure and a near-surface thermal gradient of 31 ± 15 degrees Celsius per kilometre2, 3. At temperatures above 300–450 degrees Celsius, usually found at depths greater than 10–15 kilometres, the intra-crystalline plasticity of quartz and feldspar relieves stress by aseismic creep and earthquakes are infrequent. Hydrothermal conditions control the stability of mineral phases and hence frictional–mechanical processes associated with earthquake rupture cycles, but there are few temperature and fluid pressure data from active plate-bounding faults. Here we report results from a borehole drilled into the upper part of the Alpine Fault, which is late in its cycle of stress accumulation and expected to rupture in a magnitude 8 earthquake in the coming decades4, 5. The borehole (depth 893 metres) revealed a pore fluid pressure gradient exceeding 9 ± 1 per cent above hydrostatic levels and an average geothermal gradient of 125 ± 55 degrees Celsius per kilometre within the hanging wall of the fault. These extreme hydrothermal conditions result from rapid fault movement, which transports rock and heat from depth, and topographically driven fluid movement that concentrates heat into valleys. Shear heating may occur within the fault but is not required to explain our observations. Our data and models show that highly anomalous fluid pressure and temperature gradients in the upper part of the seismogenic zone can be created by positive feedbacks between processes of fault slip, rock fracturing and alteration, and landscape development at plate-bounding faults.
ISSN: 0028-0836

hal-01685757
https://hal.archives-ouvertes.fr/hal-01685757
DOI : 10.1038/nature22355

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
We report paleomagnetic results from four overlapping stratigraphic sections (55 flows) through the lava pile in Tjörnes peninsula (North Iceland). The initial aim of our study was to check the existence of intermediate directions previously interpreted as belonging to the Matuyama-Brunhes geomagnetic reversal (MBR), and complete this record with Thellier-Thellier paleointensity determinations and 40Ar/39Ar radioisotopic dating. The directional results corroborate the findings by Kristjánsson et al. (1988): a sequence of reversed, transitional, normal, and transitional polarity (R-T-N-T) is found in each section. The polarity change is characterized by a jump from reversed virtual geomagnetic poles (VGPs) to a quasi-cluster of transitional VGPs located over China, followed by a second jump to a cluster of normal polarity before moving again towards intermediate polarities. Reversely magnetized flows exhibit a magnetic mineralogy well suited for paleointensity determination. Of these 25 flows, 20 yielded paleointensity estimates of good quality. Unfortunately, the systematic presence of coarse multidomain (MD) titanomagnetite in the transitional and normal polarity flows hampered paleointensity estimates for these flows. The Virtual Dipole Moments (VDM) calculated for the reversed polarity flows vary from 4.9 to 7.0 × 1022 Am2 with an arithmetic mean value of 5.5 ± 0.8 × 1022 Am2. This value is identical to the mean VDM obtained for the 0.3-5 Ma time window and thus strengthens the conclusion that the recent geomagnetic field strength is anomalously high compared to that older than 0.3 Ma. 40Ar/39Ar isochron ages were obtained with the incremental heating technique on groundmass separates or phenocryst-poor whole rock samples. Six of the transitionally magnetized lavas yielded isochron ages that are indistinguishable from one another at the 95% confidence level. The weighted average age is 862 ± 51 ka. In the upper part of this volcanic sequence, the normally magnetized flows yielded a weighted average age of 787 ± 40 ka. This finding supports the conclusion that a geomagnetic excursion within the Matuyama chron is recorded in the lower part of the Tjörnes volcanic sequence, whereas only the latest portion of the MBR is recorded by flows in the upper part of the lava pile. This conclusion is supported via comparison of the VGPs obtained from the lavas with those predicted by a geomagnetic field model.
ISSN: 0012-821X

hal-00618515
https://hal.archives-ouvertes.fr/hal-00618515
https://hal.archives-ouvertes.fr/hal-00618515/document
https://hal.archives-ouvertes.fr/hal-00618515/file/Camps-epsl2011.pdf
DOI : 10.1016/j.epsl.2011.07.026

Éditeur(s) : HAL CCSD
Résumé : National audience
New experimental data on electrical and anelastic properties of oceanic rocks at high pressures and temperatures are compared with the field data of seismic and electrical anomalies observed at oceanic spreading centers. The laboratory results of electrical conductivity and internal friction on gabbro and peridotites permit to estimate a degree of partial melting and maximum temperature in axial magma chambers (12-15% of melting at 1155-1165°C) and at the upper boundary of a source zone at the depth of 40-90 km (less than 5% of melting in the absence of pyroxene crystals and up to 11-14% in the pyroxene reach rocks at 1360-1370°C) at mid-oceanic ridges.
ISSN: 0934-6554

hal-01322920
http://hal.univ-reunion.fr/hal-01322920

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
An analytical solution in the Laplace domain is derived for modeling anomalous pressure diffusion during pumping tests in aquifer displaying hierarchical fractal fracture networks. The proposed solution generalizes all of the analytical models for fractal flow published previously by combining multifractal diffusion and nested multiporosity with transient exchanges, interface skin effects, and well storage effect. Solutions are derived for fracture-delimited blocks with planar, cylindrical, and spherical shapes, as well as with any fractional dimensional shapes. Any combinations of these shapes can be defined in order to model a large range of situations. Within each permeability level, the fractal properties of the fracture network can be specified and the fractal dimension can be distinct from the shape dimension of the block.
ISSN: 0043-1397

hal-00411549
https://hal.archives-ouvertes.fr/hal-00411549
DOI : 10.1029/2008WR007125