Simulations of strong ground motion in SW Iberia for the 1969 February 28 (Ms= 8.0) and the 1755 November 1 (M∼ 8.5) earthquakes – II. Strong ground motion simulations

This is the second paper of a series of two concerning strong ground motion in SW Iberia due to earthquakes originating from the adjacent Atlantic area. The aim of this paper is to use the velocity model that was proposed and validated in the companion paper for seismic intensity modelling of the 1969 ( M _s= 8.0) and 1755 ( M = 8.5–8.7) earthquakes.
First, we propose a regression to convert simulated values of Peak Ground Velocity (PGV) into Modified Mercalli Intensity (MMI) in SW Iberia, and using this regression, we build synthetic isoseismal maps for a large ( M _s= 8.0) earthquake that occurred in 1969. Based on information on the seismic source provided by various authors, we show that the velocity model effectively reproduces macroseismic observations in the whole region. We also confirm that seismic intensity distribution is very sensitive to a small number of source parameters: rupture directivity, fault strike and fault dimensions. Then, we extrapolate the method to the case of the great ( M = 8.5–8.7) 1755 earthquake, for a series of hypotheses recently proposed by three authors about the location of the epicentral region. The model involving a subduction-related rupture in the Gulf of Cádiz results in excessive ground motion in northern Morocco, suggesting that the source of the 1755 earthquake should be located further west. A rupture along the western coast of Portugal, compatible with an activation of the passive western Iberian margin, would imply a relatively low average slip, which, alone, would could not account for the large tsunami observed in the whole northern Atlantic ocean. A seismic source located below the Gorringe Bank seems the most likely since it is more efficient in reproducing the distribution of high intensities in SW Iberia due to the 1755 earthquake.     

Spatial variability and non-linearity of strong ground motion near a fault

We present observations of ground accelerations recorded at a small array close to the fault during the Düzce earthquake and its early aftershocks. The records show the strong spatial variability of ground acceleration over distances of only a few hundred metres. During the main shock, the peak horizontal acceleration values ranged from 0.3 to about 1.0 g at stations distant of 1.5 km only. We attribute this spatial variability to a fault zone site effect as peak ground acceleration steadily increases as the distance to the fault trace decreases. The spectral ratio between the ground motion recorded near the fault and the one outside the fault zone shows a shift of the spectral peak to lower frequencies with increasing peak accelerations. Such an observation suggests a non-linear behaviour of the fault zone due to the strong ground shaking. As much as a 45 per cent reduction in the shear wave velocity is necessary for the observed shifts. The opening of pre-existing cracks throughout the fault zone is the proposed mechanism to account for the observed shear wave reductions. The observation that elastic fault zone properties are soon recovered following episodes of large strains shows that cracks and fissures close rapidly after the strong shaking is over.     

Ground motion characterization and seismic hazard assessment in Spain: context,problems and recent developments

This article points out some particular features conditioning seismic hazard assessments (SHA) in Spain, a region with low–moderate seismicity. Although sized earthquakes occurred in the past, as evidenced by historical documents and neotectonic studies, no large events occurred during the last decades. The absence of strong motion records corresponding to earthquakes with magnitude larger than 5.5 is an important obstacle for the development of ground motion models constrained by local data, with the consequent difficulty in SHA studies. In this paper, some recent developments aiming at providing solutions to these difficulties are presented. Specifically, a strong motion databank containing a massive collection of accelerograms and response spectra from different configurations source-path-site corresponding to earthquakes all over the world is introduced, together with software utilities for its management. A first application of this databank is the development of specific ground motion models for Spain and for the Mediterranean region that predict peak ground accelerations as a function of several definitions of magnitude, distance and soil class. The predictive power of these ground motion models is tested by contrasting their estimates with recently recorded ground motions. The comparison between our ground-motion models with others proposed in the literature for other areas reveals a regular overestimation of the expected ground motions at Spanish sites by the non-local models. Consequently, SHA studies based in external models may overestimate the predicted hazard at the Iberian sites. In the last part of the paper a method for checking whether the response spectra proposed in the Spanish Building Code (NCSE-02) are consistent with actual accelerometric data from recent low magnitude earthquakes is applied. The spectral shapes of the Spanish Building Code NCSE-02 are compared with the response spectral shapes deduced from the available accelerograms by normalising the response spectra with the recorded PGA. It is appreciated that the NCSE-02 spectral shapes are exceeded by a large number of actual spectral shapes for short periods (around 0.2 s), a result to be taken into account in further revisions of the NCSE-02 code. The issues tackled in this work constitute not only an improvement for ground-motion characterisation in Spain, but also provide guidelines of general interest for potential applications in other regions with similar seismicity.     

Physics of basic motions in asymmetric continuum

We consider the advantages of a formalism based on concept of the asymmetric continuum and we present some equivalence theorems relating it to the asymmetric elasticity and to micropolar and micromorphic theories as founded by Nowacki, Cosserats and Eringen. We consider the basic processes in an asymmetric continuum which could be reduced to the point basic motions/deformations. The co-action of spin and shear motions is assumed to play the main role in fracturing process, while the constitutive relation between the antisymmetric stresses and rotations replaces the friction constitutive law.     

Inferred ground motions on Guadeloupe during the 2004 Les Saintes earthquake

Accurate estimates of the ground motions that occurred during damaging earthquakes are a vital part of many aspects of earthquake engineering, such as the study of the size and cause of the uncertainties within earthquake risk assessments. This article compares a number of methods to estimate the ground shaking that occurred on Guadeloupe (French Antilles) during the 21st November 2004 (M _w 6.3) Les Saintes earthquake, with the aim of providing more accurate shaking estimates for the investigation of the sources of uncertainties within loss evaluations, based on damage data from this event. The various techniques make differing use of the available ground-motion recordings of this earthquake and by consequence the estimates obtained by the different approaches are associated with differing uncertainties. Ground motions on the French Antilles are affected by strong local site effects, which have been extensively investigated in previous studies. In this article, use is made of these studies in order to improve the shaking estimates. It is shown that the simple methods neglecting the spatial correlation of earthquake shaking lead to uncertainties similar to those predicted by empirical ground-motion models and that these are uniform across the whole of Guadeloupe. In contrast, methods (such as the ShakeMap approach) that take account of the spatial correlation in motions demonstrate that shaking within roughly 10 km of a recording station (covering a significant portion of the investigated area) can be defined with reasonable accuracy but that motions at more distant points are not well constrained.     

Approaches for modeling longitudinal dispersion in pore-networks

The purpose of this study is to quantify the dispersivity in the longitudinal direction by upscaling pore scale mixing over a network domain and to verify the dispersivity with that obtained through the more rigorous upscaling technique, the Brownian particle tracking model (BPTM). We model a porous medium with a network of pore-units that are comprised of pore bodies and bonds of finite volume. Such a pore-unit is assumed to be a mixing cell with the steady state flow condition for a single fluid. Dispersivity can be obtained by solving the mixing cell model (MCM) for the concentration in each pore-unit and by averaging the concentrations for a large number of pore units (as a function of time and space). A minimal size of network that ascertains an asymptotic value of dispersivity was determined and verified with large size pore networks. This numerically computed dispersivity is compared with the results from the BPTM for the same porous medium and flow conditions. We show that the dispersivity obtained from the MCM is equally reliable for the heterogeneous pore-networks and can be estimated as a function of pore size heterogeneity. For homogeneous networks with the MCM, the iteration time step plays an important role. On the other hand, for networks with the BPTM, the assumption of intra-bond velocity profile affects the results.