A Natural Seismic Isolating System: The Buried Mangrove Effects

The Belleplaine test site, located in the island of Guadeloupe (French Lesser Antilles) includes a three-accelerometer vertical array, designed for liquefac- tion studies. The seismic response of the soil column at the test site is computed using three methods: the spectral ratio method using the vertical array data, a numerical method using the geotechnical properties of the soil column, and an operative fre- quency domain decomposition (FDD) modal analysis method. The Belleplaine test site is characterized by a mangrove layer overlaid by a stiff sandy deposit. This con- figuration is widely found at the border coast of the Caribbean region, which is exposed to high seismic hazard. We show that the buried mangrove layer plays the role of an isolation system equivalent to those usually employed in earthquake engineering aimed at reducing the seismic shear forces by reducing the internal stress within the structure. In our case, the flexibility of the mangrove layer reduces the distortion and the stress in the sandy upper layer, and consequently reduces the potential of liquefaction of the site

Fault constitutive relations inferred from the 2009-2010 slow slip event in Guerrero, Mexico

The spatiotemporal evolution of stress state is analyzed during the 2009-2010 Slow Slip Event (SSE) of Guerrero, Mexico, based on the kinematic inversion results and using an integral expression for stress changes. A linear slip weakening behavior is generally observed during the SSE with an average slope of -0.5 +/- 0.2 MPa/m regardless the perturbation due to the 27 February 2010 M-w = 8.8 Maule, Chile earthquake. This slope remains unchanged before and after the Maule earthquake. However, for some area, the friction behavior changes from slip hardening to slip weakening following the Maule earthquake. The complex trajectory between shear stress and slip velocity is fitted with a rate-and state friction law through an inversion. The direct (rate) effect (parameter A) is found to be very small, lower by an order of magnitude than the evolutional (state) effect (parameter B). The characteristic length L is obtained as 5 cm on average

Adapted numerical modelling strategy developed to support EGS deployment.

International audienceThe exploitation of Enhanced/Engineered Geothermal Systems (EGS), for electricity and/or heat production, is a promising way to increase the amount of renewable energies contribution in the energetic mix in Europe. In regard to the required production characteristics (production temperature and flowrate) for the economical viability of EGS, the favourable targeted geological systems are deep and fractured. In order to reduce the risks and the prohibitive costs linked to the depth of such geothermal systems, numerical modelling is a useful tool to understand such deep fractured systems and to help in the construction and in the management of the deep infrastructures (wells architecture, stimulation of wells, implementation of adapted network of wells). Nevertheless, this forces to a change of paradigm in comparison to « classical » reservoir modelling based on mechanics of continuum media. Indeed 3D Discrete Fracture Network (DFN) approach looks fairly adapted to catch the mechanical and hydraulic phenomena in the fractured rock mass around wells and to understand the global systems in the network of wells. The conceptualisation of the fractured rock mass is a crucial step for such DFN models not only for the geometry but also to constrain the constitutive behaviour of singularities (fault zones, fractures etc.), depending on the tectonic context. We present some results illustrating how DFNs can be used to study the EGS behaviour at several scales

Multimethod Characterization of the French-Pyrenean Valley of Bagnères-de-Bigorre for Seismic-Hazard Evaluation: Observations and Models

International audienceA narrow rectilinear valley in the French Pyrenees, affected in the past by damaging earthquakes, has been chosen as a test site for soil response characteriza- tion. The main purpose of this initiative was to compare experimental and numerical approaches. A temporary network of 10 stations has been deployed along and across the valley during two years; parallel various experiments have been conducted, in particular ambient noise recording, and seismic profiles with active sources for struc- ture determination at the 10 sites. Classical observables have been measured for site amplification evaluation, such as spectral ratios of horizontal or vertical motions between site and reference stations using direct S waves and S coda, and spectral ratios between horizontal and vertical (H/V) motions at single stations using noise and S-coda records. Vertical shear-velocity profiles at the stations have first been obtained from a joint inversion of Rayleigh wave dispersion curves and ellipticity. They have subsequently been used to model the H/V spectral ratios of noise data from synthetic seismograms, the H/V ratio of S-coda waves based on equipartition theory, and the 3D seismic response of the basin using the spectral element method. General good agreement is found between simulations and observations. The 3D simulation reveals that topography has a much lower contribution to site effects than sedimentary filling, except at the narrow ridge crests. We find clear evidence of a basin edge effect, with an increase of the amplitude of ground motion at some distance from the edge inside the basin and a decrease immediately at the slope foot

IMBeR into the future Science Plan and Implementation Strategy 2016-2025

The Integrated Marine Biosphere Research (IMBeR) project, formerly the Integrated Marine Biogeochemistry and Ecosystem Research (IMBER1) project, is a global environmental change research initiative. Since its start in 2005, IMBeR has advanced understanding about potential marine environmental effects of global change, and the impacts and linkages to human systems at multiple scales. It is apparent that complex environmental issues and associated societal/sustainability choices operate at and across the interfaces of natural and social sciences and the humanities, and require both basic, curiosity-driven research and problem-driven, policy-relevant research. Collaborative, disciplinary, interdisciplinary, transdisciplinary and integrated research that addresses key ocean science issues generated by and/or impacting society is required to provide evidence-based knowledge and guidance, along with options for policy-makers, managers and marine-related communities, to help achieve sustainability of the marine realm under global change. This recognition underlies a new vision, “Ocean sustainability under global change for the benefit of society”, to guide IMBeR research for the next decade (2016-2025). This vision recognises that the evolution of marine ecosystems (including biogeochemical cycles and human systems) is linked to natural and anthropogenic drivers and stressors, as articulated in the new IMBeR research goal to, “Understand, quantify and compare historic and present structure and functioning of linked ocean and human systems to predict and project changes including developing scenarios and options for securing or transitioning towards ocean sustainability”. To implement its new vision and goal in the next decade, IMBeR’s mission is to, “Promote integrated marine research and enable capabilities for developing and implementing ocean sustainability options within and across the natural and social sciences, and communicate relevant information and knowledge needed by society to secure sustainable, productive and healthy oceans”. This Science Plan and Implementation Strategy provides a 10-year (2016-2025) marine research agenda for IMBeR. It is developed around three Grand Challenges (GC, see Graphical Executive Summary) focusing on climate variability, global change and drivers and stressors. The qualitative and quantitative understanding of historic and present ocean variability and change (Grand Challenge I) are the basis for scenarios, projections and predictions of the future (Grand Challenge II). These are linked in Grand Challenge III to understand how humans are causing the variability and changes, and how they, in turn, are impacted by these changes, including feedbacks between the human and ocean systems. Priority research areas with overarching and specific research questions are identified for each Grand Challenge. The Grand Challenges are supplemented with Innovation Challenges (IC, see graphical executive summary) that focus on new topics for IMBeR where research is needed and where it is believed that major achievements can be made within three to five years. The Innovation Challenges also provide a means for IMBeR to adjust its focus as major science discoveries are made and new priorities arise, especially regarding scientific innovations

Analyse du potentiel sismique d'un secteur lithosphérique au nord ouest des Alpes

The north-west of the Alps is an intraplate domain with very slow deformations. So, it seems difficult to determine the probability of occurrence of a lithospheric earthquake (magnitude greater than 7) from microseismic observations. Such observations are superficial processes with little relation to deeper and bigger ones. The aim is to determine the seismic potential of a lithospheric sector north-west of the Alps, studying the stress field generated by a gravity driven model. This model is 360 km by 400 km by 230 km deep, centered on the west alpine fossil subduction and going up to the north of Strasbourg. The study of the north-west alpine structures shows the importance of the alpine orogen which generates variations in depth of the lithospheric interfaces. A study of the stress field in the basement shows a variation of principal stress directions along the strike of the Alpine chain. Even if the absolute magnitude of stresses could not be determined a relative magnitude ratio is calculated. Results underline the importance of rheology for a gravity driven model. If an elastic rheology is modeled, calculated stress directions do not match observations. However, using an elasto-plastic rheology with a realistic geometry of the lithospheric interfaces, we can obtain stress directions coherent with the data.Le nord-ouest des Alpes est un domaine intraplaque présentant de très faibles déformations. C'est pourquoi il paraît délicat de déduire la probabilité d'occurrence d'un séisme de taille lithosphérique (magnitude supérieure à 7) à partir des observations de microsismicité. De telles observations sont en effet des processus superficiels et présentent peu ou pas de lien avec des processus profonds de plus grande ampleur. L'objectif est de déterminer le potentiel sismique d'un secteur au nord-ouest des Alpes en étudiant le champ de contrainte résultant d'un chargement gravitaire. Seuls les objets de taille lithosphérique, i.e. de l'ordre de la centaine de kilomètres sont pris en compte. Un modèle de contraintes à l'échelle 360 km par 400 km par 230 km d'épaisseur, centré sur la subduction fossile des Alpes de l'ouest et s'étendant jusqu'au nord de Strasbourg, est établi. L'étude des structures du nord-ouest alpin montre l'importance de l'orogène alpin qui se retrouve, enparticulier, dans les variations de profondeur des interfaces de la lithosphère. Une étude du champ de contrainte dans le socle a permis d'identifier une rotation des contraintes principales horizontales avec l'axe des Alpes. Bien que la valeur absolue des contraintes principales n'ait pas pu être déterminée, un rapport de valeur relative est calculé. Le résultat de la modélisation montre l'importance de la rhéologie dans le cas d'un chargement gravitaire. Si une rhéologie élastique est prise en compte, les directions de contrainte calculées sont totalement différentes des observations. Par contre, l'utilisation d'une rhéologie élasto-plastique combinée à l'utilisation d'une géométrie réaliste des interfaces lithosphériques permet d'obtenir des directions de contraintes cohérentes avec les données

Seismic potential analysis of a lithospheric sector north-west of the Alps

Le nord-ouest des Alpes est un domaine intraplaque présentant de très faibles déformations. C'est pourquoi il paraît délicat de déduire la probabilité d'occurrence d'un séisme de taille lithosphérique (magnitude supérieure à 7) à partir des observations de microsismicité. De telles observations sont en effet des processus superficiels et présentent peu ou pas de lien avec des processus profonds de plus grande ampleur. L'objectif est de déterminer le potentiel sismique d'un secteur au nord-ouest des Alpes en étudiant le champ de contrainte résultant d'un chargement gravitaire. Seuls les objets de taille lithosphérique, i.e. de l'ordre de la centaine de kilomètres sont pris en compte. Un modèle de contraintes à l'échelle 360 km par 400 km par 230 km d'épaisseur, centré sur la subduction fossile des Alpes de l'ouest et s'étendant jusqu'au nord de Strasbourg, est établi. L'étude des structures du nord-ouest alpin montre l'importance de l'orogène alpin qui se retrouve, enparticulier, dans les variations de profondeur des interfaces de la lithosphère. Une étude du champ de contrainte dans le socle a permis d'identifier une rotation des contraintes principales horizontales avec l'axe des Alpes. Bien que la valeur absolue des contraintes principales n'ait pas pu être déterminée, un rapport de valeur relative est calculé. Le résultat de la modélisation montre l'importance de la rhéologie dans le cas d'un chargement gravitaire. Si une rhéologie élastique est prise en compte, les directions de contrainte calculées sont totalement différentes des observations. Par contre, l'utilisation d'une rhéologie élasto-plastique combinée à l'utilisation d'une géométrie réaliste des interfaces lithosphériques permet d'obtenir des directions de contraintes cohérentes avec les données.The north-west of the Alps is an intraplate domain with very slow deformations. So, it seems difficult to determine the probability of occurrence of a lithospheric earthquake (magnitude greater than 7) from microseismic observations. Such observations are superficial processes with little relation to deeper and bigger ones. The aim is to determine the seismic potential of a lithospheric sector north-west of the Alps, studying the stress field generated by a gravity driven model. This model is 360 km by 400 km by 230 km deep, centered on the west alpine fossil subduction and going up to the north of Strasbourg. The study of the north-west alpine structures shows the importance of the alpine orogen which generates variations in depth of the lithospheric interfaces. A study of the stress field in the basement shows a variation of principal stress directions along the strike of the Alpine chain. Even if the absolute magnitude of stresses could not be determined a relative magnitude ratio is calculated. Results underline the importance of rheology for a gravity driven model. If an elastic rheology is modeled, calculated stress directions do not match observations. However, using an elasto-plastic rheology with a realistic geometry of the lithospheric interfaces, we can obtain stress directions coherent with the data

Seismic potential analysis of a lithospheric sector north-west of the Alps

Le nord-ouest des Alpes est un domaine intraplaque présentant de très faibles déformations. C'est pourquoi il paraît délicat de déduire la probabilité d'occurrence d'un séisme de taille lithosphérique (magnitude supérieure à 7) à partir des observations de microsismicité. De telles observations sont en effet des processus superficiels et présentent peu ou pas de lien avec des processus profonds de plus grande ampleur. L'objectif est de déterminer le potentiel sismique d'un secteur au nord-ouest des Alpes en étudiant le champ de contrainte résultant d'un chargement gravitaire. Seuls les objets de taille lithosphérique, i.e. de l'ordre de la centaine de kilomètres sont pris en compte. Un modèle de contraintes à l'échelle 360 km par 400 km par 230 km d'épaisseur, centré sur la subduction fossile des Alpes de l'ouest et s'étendant jusqu'au nord de Strasbourg, est établi. L'étude des structures du nord-ouest alpin montre l'importance de l'orogène alpin qui se retrouve, enparticulier, dans les variations de profondeur des interfaces de la lithosphère. Une étude du champ de contrainte dans le socle a permis d'identifier une rotation des contraintes principales horizontales avec l'axe des Alpes. Bien que la valeur absolue des contraintes principales n'ait pas pu être déterminée, un rapport de valeur relative est calculé. Le résultat de la modélisation montre l'importance de la rhéologie dans le cas d'un chargement gravitaire. Si une rhéologie élastique est prise en compte, les directions de contrainte calculées sont totalement différentes des observations. Par contre, l'utilisation d'une rhéologie élasto-plastique combinée à l'utilisation d'une géométrie réaliste des interfaces lithosphériques permet d'obtenir des directions de contraintes cohérentes avec les données.The north-west of the Alps is an intraplate domain with very slow deformations. So, it seems difficult to determine the probability of occurrence of a lithospheric earthquake (magnitude greater than 7) from microseismic observations. Such observations are superficial processes with little relation to deeper and bigger ones. The aim is to determine the seismic potential of a lithospheric sector north-west of the Alps, studying the stress field generated by a gravity driven model. This model is 360 km by 400 km by 230 km deep, centered on the west alpine fossil subduction and going up to the north of Strasbourg. The study of the north-west alpine structures shows the importance of the alpine orogen which generates variations in depth of the lithospheric interfaces. A study of the stress field in the basement shows a variation of principal stress directions along the strike of the Alpine chain. Even if the absolute magnitude of stresses could not be determined a relative magnitude ratio is calculated. Results underline the importance of rheology for a gravity driven model. If an elastic rheology is modeled, calculated stress directions do not match observations. However, using an elasto-plastic rheology with a realistic geometry of the lithospheric interfaces, we can obtain stress directions coherent with the data

Analysis of the slow slip events of Guerrero, Mexico: implications for numerical modeling.

Guerrero, in Mexico, is one of the subduction zones where long term slow slip events (SSEs) have been observed recurrently. Understanding the mechanics of these events is important to determine their role in the seismic cycle. SSEs in Guerrero have been found to have the same characteristics, along the interface of subduction, as classical earthquakes but with much longer slip time (around a year) and lower stress drop (0.1 MPa). We investigate the slip models of the Guerrero SSEs of 2006 and 2009 (Radiguet et al., JGR 2012). The kinematic slip models have been determined by inversion of GPS time series using two different methods. From these slip histories, the constitutive relation between stress and slip (or slip rate) on each subfault is determined. Analytical Green functions are used to calculate the shear stress in a homogeneous, elastic, isotropic medium. Whatever the kinematic slip modeling method used, a clear slip weakening law can be retrieved over the whole slipping area. While some spatial variation in the parameters of the slip weakening law is observed, a mean value of about 0.1 m for the slip weakening distance and 2.5 kJ/m2 for the fracture energy can be extracted on each subfault. Moreover the slip-weakening rate seems quite homogeneous (around 1 MPa/m), and this is roughly the same as the value found in coseismic processes. The yield stress is of the order of 0.01 MPa, a low value compared to a stress drop of 0.1 MPa. The stress-slip rate relationship presents a loop trajectory coherent with the one observed in classical earthquakes. The results of these analyses are used to numerically model the Guerrero SSEs. The aim is to reproduce the slip pattern using the mechanical laws determined in the study of the slip model. If a simple slip weakening law, with parameters found above, is used, we observe a rapid progress of the crack-like slip area. This is different from the observation of the migration of localized slip. So a slowing mechanism (healing) must be introduced in addition to the slip weakening law. A pseudo-dynamic model is developed, supposing a fully plastic fault strengthening with a prefixed slip distance