Damage seismic scenarios for Angra do Heroísmo, Azores (Portugal)

Angra do Heroísmo, located in Terceira Island, is the main town of the Central Group of Azores Archipelago and it is classified as World Heritage by Unesco. The preservation of this patrimony, located inside a seismic region, includes the implementation of a policy to protect it against earthquakes. The development of credible damage scenarios could help this objective, as they can be used in prevention measures to minimize the impact of future earthquakes. Two methods were used for seismic damage estimation: the capacity spectrum method and a new method, based on the direct response of a 1D non-linear oscillator to simulate the building seismic performance. Both methods seem to overestimate the damages for the near sources, but this could be due to the empirical model used to estimate the input ground motion. Comparison with the damage observed after the January 1st, 1980 earthquake show a large discrepancy with an underestimation bias, which can be explained, mainly, by the difference of building performance and on the plausible hypothesis of 2D site effects.     

Instability structures, synsedimentary faults and turbidites, witnesses of a Liassic seismotectonic activity in the Dauphiné Zone (French Alps): A case example in the Lower Pliensbachian at Saint-Michel-en-Beaumont

Instability structures, synsedimentary faults and turbidites have been studied in the Lower Pliensbachian succession of Saint-Michel-en-Beaumont, belonging to the Taillefer block, an ancient half-graben emplaced during the Liassic Tethyan rifting. Geometrical and mechanical analyses demonstrate that the instability structures occurred thanks to movements along spineless synsedimentary normal faults, when the turbiditic and limestone layers were already case-hardened and partly fractured by tension gashes even when the mudstones were still unlithified. Both the tension gashes and the synsedimentary faults are homogeneous in strike with the major regional faults and are in good agreement with the regional direction of extension for this period. The characters of the turbiditic beds, with erosive base, graded bedding, and incomplete Bouma sequence, are in favour of a seismic origin. Instability structures, spineless synsedimentary faults and turbiditic inflows are thus considered as seismites and interpreted as the result of high seismicity periods including some events with M > 5 in the general extensive ambiance of the Liassic Tethyan rifting. The analysis of the geometrical relationships between all these sedimentary features allows to distinguish the successive stage of occurrence of an instability structure, from the sedimentation of alternating marls and limestones, and sudden turbiditic inflows, then early case-hardening of the turbidites, until the important seismotectonic event generating the spineless normal faults, themselves triggering the fall of indurated blocks and locally the forming of breccias. The Ornon Fault, which constitutes the border of the Taillefer block, 15 km eastward, played a major role during the Liassic sedimentation and may represent the major seismic fault related to the seismites occurrence in the Beaumont basin.     

Pb,Hf and Nd isotope compositions of the two Réunion volcanoes (Indian Ocean): A tale of two small-scale mantle “blobs”?

Pb, Hf, Nd and Sr isotopes of basaltic lavas from the two Réunion Island volcanoes are reported in order to examine the origin of the sources feeding these volcanoes and to detect possible changes through time. Samples, chosen to cover the whole lifetime of the two volcanoes (from 2 Ma to present), yield a chemically restricted (compared to OIB lavas) but complex distribution. Réunion plume isotopic characteristics have been defined on the basis of the composition of uncontaminated shield-building lavas from the Piton de la Fournaise volcano. The average ?_Nd, ?_Hf, ⁸⁷Sr/⁸⁶Sr and ²⁰⁶Pb/²⁰⁴Pb, ²⁰⁷Pb/²⁰⁴Pb, and ²⁰⁸Pb/²⁰⁴Pb isotope ratios calculated for this component are + 4.4, + 9.1, 0.70411, 18.97, 15.59 and 39.03, respectively. In Pb–Pb isotope space, each volcano defines a distinct linear trend but slight variations are also detected within the various volcanic sequences. The Piton des Neiges volcano yields a distinct and significantly more scattered isotopic distribution than Piton de la Fournaise for both Pb, Hf and Nd isotope tracers. A principal component analysis of the Pb isotope data from Piton de la Fournaise reveals a major contribution of the C and EM-1 components (with a clear Dupal flavor) as main components for the modern Réunion plume. The same components have been identified for Piton des Neiges but with a stronger participation of a depleted mantle component and a weaker EM-1 contribution. The compositional change of the lavas erupted by the Piton des Neiges and Piton de la Fournaise volcanoes is attributed to the impingement of two small-scale blobs of plume material at the base of the Réunion lithosphere. Compared to other hot-spots worldwide, in particular Hawaii and Kerguelen, magmas beneath Réunion are generated from a considerably more homogeneous, compositionally more primitive plume higher in ²⁰⁶Pb. Although shallow-level contamination processes have been locally detected they did not alter significantly the composition of the plume magmas. This is tentatively attributed to mantle dynamics producing small, high-velocity blobs that ascend rapidly through the lithosphere, and to the lack of a well-developed magma chamber at depth in the lithosphere.     

2D equivalent linear site effect simulation: example applications to two deep valleys

In the framework of the Sismovalp European project, an equivalent linear 2D code was developed to compute the response of a valley to SH waves, using the discrete wave-number method proposed by Aki and Larner (Aki K, Larner KL (1970) J Geophys Res 75:5). To overcome the frequency upper bound limitation, the Aki and Larner’s method is combined with a one-dimensional computation using a classical multi-layer method (Aki K, Richards PG (1980) Quantitative Seismology: Theory and Methods, vols. 1 & 2. W.H. Freeman & Co, San Francisco). The so-called “Aki–Larner extended method” is associated to an iterative algorithm, as proposed by Seed and Idriss (Seed HB, Idriss IM (1969) Report No. EERC 70–10, Earthquake Research Center, University of California, Berkeley, California) which accounts for the modulus and damping degradation using a linear visco-elastic model. A comparison of the results in the linear and the equivalent linear cases, for a magnitude 6.0 earthquake, shows that the account for the equivalent linear behaviour of the soil significantly reduces the amplification level, especially at frequencies higher than the fundamental resonance frequency of the site. In the case of site effects or microzonation studies devoted to produce design spectra for engineering structures, this can have a major impact on the associated results and costs, depending on the frequency of interest for the considered structure. As a first application of the developed technique, 2D equivalent linear Aki–Larner computations are used to perform the seismic microzonation study of the upper Rhone valley, in the Visp area (Switzerland), a typical 2D alpine valley. These investigations made it possible to determine site specific spectra, associated with different zones, to be used instead of the code spectra that do not take into account the local 2D amplification.     

Three-dimensional hydromechanical modeling of internal erosion in dike-on-foundation

Currently, numerical studies at the real scale of an entire engineering structure considering internal erosion are still rare. This paper presents a three-dimensional (3D) numerical simulation of the effects of internal erosion within a linear dike located on a foundation. A two-dimensional (2D) finite element code has been extended to 3D in order to analyze the impact of internal erosion under more realistic hydromechanical conditions. The saturated soil has been considered as a mixture of four interacting constituents: soil skeleton, erodible fines, fluidized fine particles, and fluid. The detachment and transport of the fine particles have been modeled with a mass exchange model between the solid and the fluid phases. An elastoplastic constitutive model for sand-silt mixtures has been developed to monitor the effect of the evolution of both the porosity and the fines content induced by internal erosion upon the behavior of the soil skeleton. An unsaturated flow condition has been implemented into this coupled hydromechanical model to describe more accurately the seepage within the dike and the foundation. A stabilized finite element method was used to eliminate spurious numerical oscillations in solving the convection-dominated transport of fluidized particles. This numerical tool was then applied to a specific dike-on-foundation case subjected to internal erosion induced by a leakage located at the bottom of the foundation. Different failure modes were observed and analyzed for different boundary conditions, including the significant influence of the leakage cavity size and the elevation of the water level at the upstream and downstream sides of the dike.     

A procedure to analyze nonlinear density waves in Saturn's rings using several occultation profiles

Cassini radio science experiments have provided multiple occultation optical depth profiles of Saturn's rings that can be used in combination to analyze density waves. This paper establishes an accurate procedure of inversion of the wave profiles to reconstruct the wave kinematic parameters as a function of semi-major axis, in the nonlinear regime. This procedure is established using simulated data in the presence of realistic noise perturbations, to control the reconstruction error. It is then applied to the Mimas 5:3 density wave. There are two important concepts at the basis of this procedure. The first one is that it uses the nonlinear representation of density waves, and the second one is that it relies on a combination of optical depth profiles instead of just one profile. A related method to analyze density waves was devised by Longaretti and Borderies [Longaretti, P.-Y., Borderies, N., 1986. Icarus 67, 211-223] to study the nonlinear density wave associated with the Mimas 5:3 resonance, but the single photopolarimetric profile provided limited constraints. Other studies of density waves analyzing Cassini data [Colwell, J.E., Esposito, L.W., 2007. Bull. Am. Astron. Soc. 39, 461; Tiscareno, M.S., Burns, J.A., Nicholson, P.D., Hedman, M.M., Porco, C.C., 2007. Icarus 189, 14-34] are based on the linear theory and find inconsistent results from profile to profile. Multiple cuts of the rings are helpful in a fundamental way to ensure the accuracy of the procedure by forcing consistency among the various optical depth profiles. By way of illustration we have applied our procedure to the Mimas 5:3 density wave. We were able to recover precisely the kinematic parameters from the radio experiment occultation data in most of the propagation region; a preliminary analysis of the pressure-corrected dispersion allowed us to determine new but still uncertain values for the opacity (K?0.02 cm²/g) and velocity dispersion of (c₀?0.6 cm/s) in the wave region. Our procedure constitutes the first step in our planned analysis of the density waves of Saturn's rings. It is very accurate and efficient in the far-wave region. However, improvements are required within the first wavelength. The ways in which this method can be used to establish diagnostics of ring physics are outlined.     

Site effects in Mexico City eight years after the September 1985 Michoacan earthquakes

The purpose of this paper is to take a comprehensive look at site effects in Mexico City for the 1985 Michoacan earthquake. We examine, successively, 1D and 2D models. For the latter, we consider in detail both large scale and small scale heterogeneities, using extensively the Aki-Larner wave propagation method, in the version given by Bard and Gariel. In particular, we make a critical review of the different explanations proposed for the large duration of strong ground motion in the lake zone. Our purpose is two-sided. We first outline the difference between what is well established and what remains still unexplained regarding the seismic response of Mexico City basin. On the other hand, we wish to make explicit the conditions that the proposed models require to explain strong motion duration. Our results allow us to qualify the models proposed to date and to point out what could be the experiments and the new data required to find a truly satisfactory explanation of strong ground motion at Mexico City.     

Theory of metamorphic segregation and related processes

Metamorphic segregation is defined here as the formation and growth of bands or domains of different bulk compositions within an originally unbanded rock. It can result from an instability arising in some deforming rocks when diffusion transfer is significant. The nature of this instability is demonstrated separately for a differentiation associated with a crenulation (Type C) and for one without crenulation (Type L). However, the continuous gradation between the two types observed in nature is also expected from the models.A requirement for both types is that one chemical component of the rock, typically silica, diffuse more rapidly than others, e.g. phyllosilicate components, in response to stress-induced pressure gradients. In addition, Type L requires that, at least at an early stage, portions of the rock rich in this mobile component, e.g. rich in quartz, be more competent than those portions which are poor in quartz. By a process akin to the development of pressure shadows, silica diffuses toward domains which are already the richest in quartz. Alternatively, if the rock is chemically open to loss of silica, this silica preferentially originates from dissolution of quartz in the more mica-rich domains. In either case the compositional contrast between domains is increased and metamorphic segregation results. Type C differentiation is best explained if we accept the suggestion made by many petrographers that layer silicates catalyse the pressure-induced transfer of silica. This suggestion can in fact account for other features of metamorphic segregation bands.Metamorphic segregation should also proceed during the development of quartz segregations, quartz rods, slaty cleavage and mylonite banding. The stress distribution argued for Type L bands may also occur in many banded migmatites. In migmatites. however, the assemblage crystallizing in the leucosome may be derived from a circulating hydrothermal fluid as well as from diffusion transfer.     

Strong motion dataset of Turkey: data processing and site classification

Strong motion records of Turkey are studied in order to prepare a catalog to be used as a database for further studies (for instance empirical attenuation laws). The network started to be installed in 1973, and the first record was obtained in 1976. The instruments are of SMA-1 analog recorders and SIG SM-2 and GeoSys GSR-16 digital types. Out of a total of 426 records released on the web sites of the General Directorate of the Disaster Affairs and of Kandili Observatory, a set of 210 records was selected with a satisfactory quality, for which it was possible to associate correctly determined source parameters (source magnitudes and epicentral distances). Most of the records are obtained from around North and East Anatolian Fault zones, as well as from western and southwestern parts of Turkey. The main outcome of this paper is a strong motion catalog of Turkey, with the indication of site conditions, of the frequency band of the reliability of the records, peak values of acceleration, velocity and displacements, source parameters (magnitude, epicentral and macroseismic distances), intensity and finally the fault plane solutions whenever possible. The aim is to have, with other regional dataset, a homogenous and good quality dataset.     

Field temperature measurements at Erta'Ale Lava Lake, Ethiopia

The shield volcano Erta'Ale, situated in the Danakil Depression, Ethiopia, is known for its active lava lake. In February 2001, our team visited this lake, located inside an 80-m-deep pit, to perform field temperature measurements. The distribution and variation of temperature inside the lake were obtained on the basis of infrared radiation measurements performed from the rim of the pit and from the lake shores. The crust temperature was also determined from the lake shores with a thermocouple to calibrate the pyrometer. We estimated an emissivity of the basalt of 0.74 from this experiment. Through the application of the Stefan-Boltzmann law, we then obtained an estimate of the total radiative heat flux, constrained by pyrometer measurements of the pit, and visual observations of the lake activity. Taking into account the atmospheric convective heat flux, the convected magma mass flux needed to balance the energy budget was subsequently derived and found to represent between 510 and 580 kg s^-1. The surface circulation of this mass flux was also analyzed through motion processing techniques applied to video images of the lake. Electronic supplementary material to this paper can be obtained by using the Springer LINK server located at http://dx.doi.org/10.1007/s00445-002-0224-3.