SEISMIC DEMANDS AS A RESULT OF DIRECTIVITY EFFECTS FROM THE MW 5.1 2011 LORCA (SPAIN) EARTHQUAKE

Carlos GORDO-MONSO1, Eduardo MIRANDA2

ABSTRACT

The May 11th 2011, Lorca earthquake in Southeastern Spain was a damaging event, especially considering its relatively low magnitude (Mw 5.1), which produced 9 fatalities, more than 300 injuries, and more than 462 million Euros in economic loses. Recorded ground motions and corresponding spectral ordinates significantly exceed those in current Spanish probabilistic seismic hazard models, as well as those in the Eurocode and Spanish building codes. A pulse directivity effect has been postulated by several investigators to explain these differences, the objective of this paper is to assess such effects on ground motions recorded at the LOR station during the 2011 Lorca earthquake, and to evaluate the significance of these effects on earthquake resistant design in moderate seismic regions. First, we studied the likelihood of the presence of a directivity pulse, by comparing different parameters of the recorded ground motions to analytical pulses. Additionally, we present a novel technique to extract directivity pulse features, based on the match of the response spectrum of a simple analytical pulse and the spectrum of the recorded ground motion. In a second part, we relate the inelastic displacement spectra to some models that try to capture the displacement demand features of earthquakes presenting directivity-pulse characteristics. We conclude that a large part of the impact caused by this relatively small earthquake owes to directivity effects, an aspect which is typically ignored, both in probabilistic seismic hazard analysis and in most building codes.

Keywords: Directivity; Directionality; Near-Fault; Pulse-Like Ground Motions; Inelastic Spectra

1. INTRODUCTION

The Lorca, May 11th 2011, earthquake in South-Eastern Spain was a damaging event, especially considering its relatively low magnitude Mw 5.1 (Lopez-Comino et al, 2012) (Cabañas et al., 2013). The earthquake lead to a death-toll of 9 people, more than 300 people injured, and approximately €462 M in direct economic losses (Alvarez-Cabal et al., 2014). Damages included the collapse of a 4-story modern reinforced concrete building (Fig 1a) and the partial collapse of an 18th century church (Fig. 1b), along with extended damages in nonstructural elements.

Figure 1. (a) Collapsed modern reinforced concrete building. (b) Partial collapse of the Santiago Church, a masonry building dating from the XVIII century (from Cabañas-Rodriguez et al., 2011)

1Assistant Professor, ETSI Caminos. Universidad Politécnica de , Madrid, Spain, [email protected] 2Professor, Civil and Environmental Eng. Dept., Stanford University, Stanford, USA, [email protected]

The focal mechanism solution corresponds to a reverse and strike-strike slip faulting mechanism with a very shallow crustal depth of approximately 3 km (Cabañas-Rodriguez et al., 2011), to 4 km (Martinez- Solares et al., 2012). The maximum PGA in the event, which was recorded at the Lorca station, was 0.37·g in the N30W component and is more than three times larger than the one specified in the Spanish code for this region (based on a 10% probability of exceedance in 50 years) being also the largest ever recorded in Spain (Cabañas et al., 2013) (Moratto et al., 2014). A possible induced origin for this earthquake, related to groundwater crustal unloading, has been pointed as a cause by some researchers (González et al. 2012). However, the study presented in this paper does not address the specific causes (natural or induced stress build-up) that may have led to the earthquake triggering.

It is well known that rupture directivity effects can lead to strong pulse-like ground motions (Bertero et al., 1978), (Anderson and Bertero, 1987), (Hall et al., 1995), (Iwan, 1998), (Alavi and Krawinkler, 2001). The propagation of fault rupture toward a site, at a velocity close to the shear wave velocity, causes most of the seismic energy from the rupture to arrive in the form of a single large pulse of motion that occurs early in the seismic record (Somerville et al., 1997). Furthermore, the radiation pattern of the shear dislocation on the fault causes these large pulse-like ground motions that tend to be oriented in the direction perpendicular to the fault plane, leading to more intense ground motions in the normal component than in the parallel component.

Although near-fault pulse-like records have been included in the development of conventional ground motion prediction models (GMPE), they do not properly account for directivity effects because most of them are aimed at predicting the geometric median of the intensity of the two horizontal components, therefore systematically underestimating the intensity of ground motions affected by directivity in the maximum direction, which is typically the fault-normal component. Furthermore, the standard deviation of most ground motion prediction models is averaged over all distances leading to underestimation of the dispersion, especially near the fault (Abrahamson, 2000). Some more recent GMPM have explicitly incorporated directivity effects (e.g. Somerville et al., 1997) but modifications are only introduced for events with magnitudes larger than 6.5. Recently, such directivity effects have been observed in earthquakes with magnitudes between 6.0 and 6.5 such as the 2004 Mw 6.0 Parkfield, earthquake (Shakal et al., 2005), the 2009 L’Aquila Mw 6.3 earthquake (Chioccarelli and Iervolino, 2010) or the 2011 Mw 6.3 Christchurch earthquake (Bradley et al., 2014), pointing out that these damaging pulse-like ground motions can occur in earthquakes with magnitudes smaller than 6.5.

Several researchers (Lopez-Comino, 2012) (Rueda-Nunez et al., 2012) (Rueda-Nunez, 2014) (Alguacil et al., 2014) (Pro et al., 2014) have conjectured the presence of a directivity pulse in the 2011 Lorca earthquake despite its relative small magnitude. The objective of this paper is to evaluate directivity effects on ground motions recorded during the 2011 Lorca earthquake, and to evaluate the significance of these effects on earthquake resistant design in moderate seismic regions. In the first part of this paper, we study the likelihood of the presence of that pulse, by conducting a comparison of different parameters of recorded ground motions to analytical pulses. In the second part of the paper, we relate the recorded ground motion and its inelastic displacement spectra to some of the most recent statistical models that try to capture the displacement demand features of earthquakes presenting directivity-pulse characteristics.

2. LORCA 2011 EARTHQUAKE GROUND MOTION DATA

The accelerograph station in the city of Lorca recorded two horizontal ground motion signals in the N30W and N60E directions (Fig. 2a), which roughly match the fault-normal and fault-parallel directions of the causative Alhama de Fault (AMF). The corresponding elastic pseudo-acceleration spectra and their geometric mean are shown in Fig. 2b. It can be seen that spectral ordinates in the fault-normal direction over a wide range of periods are more than twice than those in the fault-parallel direction. This figure clearly illustrates how, if one were to design based on the geometric mean, one would significantly underestimate the intensity of the motion in the fault-normal direction and overestimate the motion in the fault-parallel direction. 2

[g] a Spectral Acceleration S

Figure 2. (a) Ground motion acceleration recorded at the Lorca station; (b) 5% damped spectra for the fault- normal and fault-parallel components, and their corresponding geometric mean, plotted along with the uniform hazard spectrum for TR=475yr. (Gaspar-Escribano et al., 2008), (Benito et al., 2012), the Spanish Seismic Building Code (NCSE-02) spectrum, and Eurocode 8 type 1 and 2 spectra for an equivalent soil classification

Several researchers (Martínez-Díaz et al. 2012) (Rueda et al., 2014) have been able to determine the approximate slip distribution of the fault area of rupture indicating that the rupture propagated from the hypocenter to an upward-southwest direction in a strike-reverse mechanism, toward the city of Lorca where the recording station is located. The Joyner and Boore distance is approximately Rjb=1.26 km, and no hanging-wall amplification effect is expected in the recorded ground motion, given that the recording site lies on the footwall side of the surface projection of the top edge of rupture.

According to the Spanish Earthquake-Resistant Building Code site classification, the recording site is located on a Class II, (Blazquez-Martinez et al., 2014), corresponding to medium to soft rock, with shear wave velocity 400 m/s ≤ Vs,30 ≤ 750 m/s. For the rest of the analysis, the site has been considered to have a Vs,30 = 575 m/s corresponding to the average of the upper and lower boundaries of site class II.

3. LORCA 2011 EARTHQUAKE NEAR-FAULT PULSE-LIKE GROUND MOTION

Following the methodology proposed by Shahi and Baker (2013) to search for pulse-like patterns in earthquake ground motions and extract its prominent features (pulse period and pulse velocity and acceleration components), based on a wavelet decomposition of the ground velocity time-history, it was found that a distinct pulse feature dominated the ground motion recordings (Fig. 3b). The extracted pulse signal shows a pulse period Tp=0.48s, a result identical to that obtained by (Rueda, 2014) when applying the algorithm by Baker (2007).

Figure 3. (a) Rupture surface projection (adapted from Martinez-Diaz et al. 2012). (b) Ground motion velocity time history of the record signal plotted along with the Db4 wavelet and with the M&P wavelet 3

Additionally to the method proposed by Shahi and Baker (2013), which is based on the use of the Daubechies 4 wavelet (Db4), a similar pulse extraction was performed using the Mavroeidis and Papageorgiou (2003) wavelet (M&P). The M&P wavelet is a modification to the Gabor wavelet, except that instead of having a Gaussian time modulation it has a harmonic modulation facilitating closed-form solutions to oscillators subjected to this pulse. This wavelet, has the advantage of being defined by a small number of input parameters which have an unambiguous physical meaning. The parameter Tp controls the period of the pulse, A controls the amplitude of the harmonic, while  controls the width/duration of the modulating function which, together with Tp control the number of zero crossings in the wavelet, and finally the parameter  defines the phase of the amplitude-modulated harmonic. The results of this M&P wavelet pulse extraction show a pulse similar in shape and coincident in time with the previous one, although with a pulse period somewhat larger of Tp=0.62s. The other M&P wavelet parameters were found to be =2 and =0 (Fig. 3b).

Besides these two pulse-feature extraction methods based on the comparison between the ground motion record and the analytic pulse in the time domain, a different approach based on the comparison of the response spectra was used. In this complementary approach, the parameters are selected to minimize the differences between the pseudo-velocity response spectrum of the earthquake record and the pseudo- velocity response spectrum of an M&P analytic pulse as modified by Alonso and Miranda (Alonso- Rodriguez and Miranda, 2015). The likely pulse parameters are found by performing a least-squares fit between both spectra. With this approach, the most likely M&P pulse parameters were found to be Ap=0.041·g, fp=1.49 Hz, gp=1.00, a=0.18 (Fig. 4).

Table 1 summarizes the period values obtained by each extraction methodology. The differences in period, which can be attributed to the uncertainty intrinsic to each methodology, become relevant when assessing the significance of the pulse in terms of the intended intensity measure: if this were a response spectrum, the spectrum match would provide a better fit.

[g] a Spect. Acceleration S [cm/s] v Spect. Velocity S Velocity Spect. [g] g Acceleration a

Figure 4. (a) Pseudo-acceleration 5% damped spectra for the N30W component, and for the M&P wavelet with pulse parameters that match the pseudo-velocity spectrum. (b) Pseudo-velocity spectra for the fault-normal component, and for the M&P wavelet. (c) Acceleration record of the fault-normal component and of the velocity-spectrum fitted pulse

4

Fault Normal Component (N30W) | LOR Station 0.4 a) 12 0.3

0.2

0.1

0.0

-0.1

-0.2 Ratio of Squared Cumulative Velocity Ratio

-0.3

-0.4 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 Time t [s]

Figure 5. (a) Fault-Normal record component, where intervals 1 and 2 show significant zero-crossings of the pulse. (b) Normalized cumulative squared velocity for the fault-normal component and the M&P pulse, where intervals 1 and 2 corresponding to zero-crossings amount 79% of the CSV

Once that the residual ground motion has been extracted, the cumulative squared velocity (CSV) was computed, in order to apply the criterion by Baker (2007) and to identify if the build-up of pulse energy arrived before than a large build-up of the original signal energy occurs. The limits of obtaining a 10% of the CSV in the extracted pulse before obtaining a 20% of the CSV in the original ground motion are satisfied, as can be appreciated in Fig. 5.b for the M&P pulse.

Table 1. List of pulse periods obtained with each extraction method.

Extraction Method Pulse Period Tp [s] Tg/Tp Wavelet Db4 0.48 1.10 Wavelet M&P 0.62 0.85 Velocity Spectrum Match 0.67 0.79

The zero-crossings of the significant pulse are shown in intervals 1 and 2 of Fig 5a, while Fig 5b shows the contribution to the cumulative squared velocity (CSV) for this interval totaling a 79% of the CSV, along with the corresponding CSV for the extracted pulse. Comparing the previous results, we can observe a very good fit both for the Mw-Tp relationship obtained by Shahi and baker (2013) for the pulse- like records in the NGA West2 database, as well as for the Mavroeidis and Papageorgiou (2003) expression for their ground motion catalogue (Fig 7a).

Figure 6. (a) Elastic pseudo-acceleration 5% damped spectra for the fault-normal component, and Db4 and M&P pulses (b) Elastic pseudo-velocity spectra for the fault-normal component, and Db4 and M&P pulses.

5

T - M Relationship p w a) 101

100

NGA-West2 database M&P 2003 database Shahi & Baker 2013 M&P 2003 with self-similarity Lorca 2011 Earthquake (Db4) Lorca 2011 Earthquake (M&P) 10-1 5.0 5.5 6.0 6.5 7.0 7.5 8.0 Magnitude M w

Figure 7. (a) Pulse period versus earthquake magnitude for the NGA-West2 database records containing a pulse, the M&P database and Lorca 2011 earthquake, plotted along with the (Shahi and Baker, 2011) and (Mavroeidis and Papageorgiou, 2003) expressions for pulse-period expectation. (b) Spectral acceleration amplification factor Af due to directivity pulse (Shahi and Baker, 2011), and pulse to residual ratios for each wavelet.

Shahi and Baker (2013) proposed a local amplification factor Af=Sa,pulse/Sa,residual which amplifies spectral ordinates in the neighborhood of the pulse period. For this purpose, they adopted the functional form previously proposed by Ruiz-Garcia and Miranda (2005). This amplification factor is plotted in Fig. 7b for the Lorca pulse and residual motions, along with the Shahi and Baker expression, which captures well the main characteristic of the pulse, although the Lorca earthquake is somewhat larger.

4. LORCA 2011 EARTHQUAKE FORWARD DIRECTIONALITY

The difference in terms of spectral acceleration that can be present in a given pair of orthogonal ground motion records of the same event, and their relation to the spectral acceleration predicted by some GMPE model, has been pointed out by different authors (e.g., Baker and Cornell, 2006; Boore, 2010). The usual intensity measure of many attenuation models is the geometric mean of two spectral ordinates computed from as-recorded orthogonal records. This geometric mean value can depart significantly from the spectral acceleration in an orientation where a maximum would occur, especially in those cases where large polarization of the ground motion can be present. Boore (2010) developed direction-independent measures spectral ordinates. SaRotD50 corresponds to the median of spectral ordinates at all possible orientations, while SaRotD100 and SaRotD00 are intensity measures that provide information on what the maximum and minimum spectral acceleration are in any possible direction, respectively. Fig. 8 shows polarization plots of the Lorca 2011 earthquake for three chosen periods, one being the pulse period Tp=0.48s, where the magenta discontinuous line shows spectral acceleration amplitude at any given direction for that particular period.

Lorca 2011 Displacement Response Trace (T=0.48s) 90 b) 1 120 60 0.8

0.6 150 30 0.4

0.2

180 0

210 330

Sa/Sa 240 300 RotD100 / re sp max 270 Fault trace

Displacement along fault-parallel axis / X max

Figure 8. Normalized displacement polarization plots, where resp is the elastic displacement response and max is the maximum displacement in each direction for a SDOF oscillator with vibration period of (a) T=0.20s, showing maximum polarization at an angle fairly close to the fault-normal direction (N30W recorded ground motion). (b) T=0.48s, the pulse period, showing maximum polarization at an angle 30º East of fault-normal direction. (c) T=1.00s, showing maximum polarization at an angle 25º East of fault-normal direction

6

Lorca 2011 Pseudo-Acceleration Spectra | LOR Station 1.2 Fault-normal component (N30W) Fault-parallel component (N60E) 1.0 Geometric mean

Sa RotD50 Drot00

[g] Sa /Sa

a Drot50 0.8 Sa Drot100 RotD100

0.6

0.4 Spectral Acceleration S Acceleration Spectral

0.2 Amplification factor Sa a) 0.0 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 Period T [s]

Figure 9. (a) Acceleration 5% damped spectra for both the fault-normal and fault-parallel components, their geometric mean, and orientation independent spectra for the 0%, 50%, and 100% percentiles. (b) Spectral amplification factor SaDRot100/SaDRot50 due to polarization of the response for the Lorca 2011 earthquake, plotted along with the median and ± of the Shahi and Baker 2014, model, and the theoretical maximum value of

SaDRot100/SaDRot50=√2 corresponding to a completely polarized response

From this graph, one can observe a strong polarization of the ground motion, especially for orientations close to 25º East of the fault-normal direction. The orientation independent spectra corresponding to the minimum, median, and maximum accelerations are plotted in Fig. 9a, along with the as-recorded N30W and N60E components and their respective geometric mean. It can be seen that the fault-normal component practically corresponds to the maximum spectral ordinates in any direction SaRotD100.

The degree of directionality or polarization of the ground motion can be analyzed by the ratio SaDRot100/SaDRot50, which is a measure of how large the strongest component is respect to the median value obtained for any arbitrary direction. This directionality amplification factor can be shown to be, at maximum √2=1.41 for a perfectly polarized ground motion. Fig. 9b shows the SaDRot100/SaDRot50 amplification factor computed for the recorded Lorca 2011 earthquake, along with the expression proposed by (Shahi and Baker, 2014) and ± one standard deviation (±) ranges. From the plot we can observe how the Lorca 2011 earthquake presented strong directionality effects, larger than most of the ones derived empirically from the NGA-West2 database, and for many periods of engineering interest very close to a perfect polarization.

5. LORCA 2011 EARTHQUAKE AS PREDICTED BY NGA-WEST2 GMPE

Shai and Baker (2011, 2014) proposed a statistical model to take into account both the polarization and pulse directivity effects (directionality and forward-directivity), departing from the Boore and Atkinson (2008) GMPE (B&A) and including adjusting terms. Figure 13 shows the results computed for a reverse- slip fault mechanism earthquake, Mw 5.1, Joyner and Boore distance Rjb=1.26km, and a Vs,30=575m/s, which likely correspond to the characteristics of the recorded ground motion at the Lorca accelerometer. The graphs shown in Figure 10 correspond to the B&A GMPE without any directionality or directivity adjustments, with each of these adjustments, and with both adjustments. From Figure 10 we can note that, although the value of PGA and spectral ordinates for very short periods (e.g. between 0 and 0.2 s) show differences from the expected median value, these still fall within the uncertainty range of this particular GMPE. However, the results show how the recorded N30W ground motion departs significantly from the B&A 2008 GMPE, with many spectral ordinates (e.g., between 0.2 and 0.8s) being more than two standard deviations above the median values, and how the adjustments proposed by Shahi and Baker present reasonably good agreement with the recorded ground motion if both directionality and directivity are taken into account.

7

[g] [g] a a Spectral Acceleration S Spectral Acceleration S

Shahi & Baker 2011 and 2014 Corrections - Yes Pulse, Yes Directionality 1.8 d) Shahi & Baker 2011, 2014 1.6 Median Shahi & Baker 2011, 2014 Lorca 2011 earthquake

1.4 [g]

a 1.2

1.0

0.8

0.6 Spectral Acceleration S

0.4

0.2

0.0 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 Period T [s]

Figure 10. (a) Boore and Atkinson (2008) 5% damped median spectral acceleration and ± intervals. Response spectra for the Lorca N30W component, (b) The Shahi and Baker (2014) modified model of spectral acceleration to the B&A (2008) model, considering the amplification due to directionality of the strongest component, (c) The Shahi and Baker (2011) modified model, considering the amplification due to the near-fault forward directivity pulse. (d) The Shahi and Baker (2011, 2014) modified for both directionality and directivity pulse

Some other nearby stations records (in particular AM2, MUL, ZAR, VLR, OLU, and VER, in a distance range of 30 km to 50 km), and corresponding response spectra, have been carefully examined in the same terms as for the LOR station. Unfortunately, and as indicated by Moratto et al. (2014), the location of all other stations are much farther away from the source which, in combination with the rather small magnitude of this event, make the study of possible directivity effects at the other stations far more difficult. In this sense, for these other stations a clear pulse signature has not been found as for the LOR station case.

6. LORCA 2011 EARTHQUAKE INELASTIC SPECTRA

Baez and Miranda (2000) showed that modifying the elastic spectral ordinates to take into account directivity effects in near-fault pulse-like ground motions was not enough because most structures are not designed to remain elastic and because the ratio of inelastic to elastic spectral ordinates of these types of ground motions can differ from those recorded away from the rupture. Various authors (e.g., Baez and Miranda 2000, Akkar et al. 2004; Ruiz-García and Miranda 2005; Ruiz-García 2011; Iervolino et al. 2012) have proposed expressions for describing inelastic to elastic displacement ratios for near- fault pulse-like ground motions. Figure 11 shows the displacement time-histories corresponding to systems with periods ranging from 0.01s to 2.00s, force reduction factors (R=1, 2, 4, 5, 6), and a post- yield stiffness ratio of 3% for the Lorca 2011 N30W component for different period ratios T/Tp (short period range, pulse period range, and long period range). In all cases, the displacement time-histories seem dominated by a large abrupt inelastic displacement increment soon after the 6s mark, rather than a steady increase in inelastic displacement produced by consecutive pulses.

8

Displacement history T/T =0.4 Displacement history T/T =1.0 p p 30 80 R=1.0 R=1.0 a) R=2.0 b) R=2.0 60 20 R=4.0 R=4.0 R=5.0 R=5.0 R=6.0 R=6.0 40 10

20 0 0

-10 -20

-20 -40

-30 -60 0 2 4 6 8 10 12 14 16 0 2 4 6 8 10 12 14 16 Time t [s] Time t [s]

Figure 11. (a)(b)(c) Displacement time-histories for increasing force reduction factors for systems with natural period T=0.20s (T/Tp=0.40), T=0.48s (T/Tp=1.0), and T=1.00s (T/Tp=2.10)

Ruiz-García and Miranda (2003) proposed estimating the ratio of inelastic to elastic spectral ordinates with an expression, which formed the basis of the simplified expression (1) included in ASCE 41-06, where R is the ratio of lateral strength required to maintain the system elastic to the lateral strength of the structure, T is the fundamental period of vibration of the structure, and a is a parameter that depends on site conditions..

11 ∙ (1) ∙

Inelastic to elastic displacement ratio C for R=6 R 4.5 Lorca 2011 earthquake 4.0 Median ASCE 41-06 Median Ruiz-Garcia & Miranda, 2003 3.5 Median Ruiz-Garcia, 2011 Median Iervolino et al. 2012 + Iervolino et al. 2012

elast 3.0

/ 2.5

inelast,R 2.0 =

R 1.5 C

1.0 0.5 b) 0.0 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 Period T [s] elast elast / / inelast,R inelast,R = = R R C C

Figure 12. Inelastic displacement ratio CR for different force reduction factors for the Lorca earthquake, plotted along with the median prediction by Ruiz-García & Miranda (2003), ASCE/SEI (2007), Ruiz-García (2011), and with the median and + prediction by Iervolino et al. (2012) (a) for R=4. (b) for R=6, (c) for R=4 with normalized abscissa T/Tp, (d) for R=6 with normalized abscissa T/Tp

Ruiz-Garcia and Miranda (2005); Ruiz-García, (2011), proposed modifications to the previous equation to estimate inelastic displacements of structures near major faults as follows: 9

11 ∙ ∙ ∙ ∙ 0.08 (2) ∙

Similarly, Iervolino et al. (2012) proposed a similar functional form as a function of the period normalized to the pulse period T/Tp but added an extra term to capture that local maximum in displacement demand. All the four expressions (ASCE/SEI 2007, Ruiz-García and Miranda 2003, Ruiz- García 2011, and Iervolino et al. 2012) have been evaluated (with values Tg=0.56s, and Tp=0.48s) for force reduction factors of R=[4, 6], and then compared to the inelastic to elastic displacement ratios CR obtained for the Lorca 2011 earthquake. The apparent trends show that the Lorca 2011 earthquake presented larger inelastic displacements than what could be expected at all range of periods and force reduction factors for the Iervolino et al. (2012) expression (in the order of +), while the Ruiz-García and Miranda (2005) expression shows a better fit in predicting the inelastic displacements. The Ruiz- García and Miranda, (2005) and Ruiz-García (2011) expression shows a more conservative trend, predicting larger displacements than observed for T≥0.25 s. However none of the aforementioned expressions is able to capture the large local maxima occurring at T=0.20 s, or conversely T/Tp=0.40.

7. CONCLUSIONS

The Lorca, May 11th 2011, earthquake recorded ground motions have been analyzed both in terms of spectral accelerations and inelastic displacement spectra. Response spectra ordinates were significantly larger than those expected using probabilistic seismic hazard analyses, and those specified in the Spanish building code and the Eurocode. We conclude that these unexpectedly large accelerations were caused by a strong forward-directivity pulse with a period close to Tp=0.6s. This finding has important consequences for earthquake resistant design because directivity effects are currently not explicitly addressed in most seismic codes. Furthermore, this earthquake illustrates that important directivity effects can occur even in Mw 5.1 events.

Analyses conducted in this study show that two relatively simple wavelets are capable of reproducing the main features of the time series of the fault-normal component of the recorded ground motions as well as the main characteristics of response spectral ordinates. Furthermore, modifications recently proposed for considering directionality and directivity effects based on NGA2 ground motions, despite currently being recommended only for seismic events with Mw>6.0, together with expressions of inelastic displacements ratios specific for near-fault pulse-like ground motions, if they are applied to these smaller magnitude events, are capable of capturing the main features of both elastic and inelastic spectra of the motion recorded in the fault-normal direction.

8. REFERENCES

Abrahamson, N. A. (2000) Effects of rupture directivity on probabilistic seismic hazard analysis, in Proceedings of the 6th International Conference on Seismic Zonation (Vol. 1, pp. 151-156). CA: Palm Springs. Akkar, S. D., Yazgan, U. Gulkan, P., (2004) Deformation limits for simple non-degrading Systems subjected to near-fault ground motions, in Proceedings 13th World Conference on Earthquake Engineering, Vancouver, Canada. Alavi, B., Krawinkler, H. (2001) Effects of near-fault ground motions on frame structures, Technical Report 138, John A. Blume Earthquake Engineering Center, Stanford, CA. Alonso-Rodríguez, A. Miranda, E. (2015) Assessment of building behavior under near-fault pulse-like ground motions through simplified models. Soil Dynamics and Earthquake Engineering 79, 45-58. Alvarez-Cabal, R., Diaz-Pavón, E., Rodríguez-Escribano, R., (2013) El Terremoto de Lorca, Efectos en los Edificios, Consorcio de Compensación de Seguros, Ministerio de Economía y Competitividad, Madrid, Spain. (in Spanish).

10

Anderson, J. C., Bertero, V. V. (1988) Uncertainties in establishing design earthquakes. Journal of Structural Engineering. 113(8), 1709-1724. ASCE/SEI Seismic Rehabilitation Standards Committee (2007) Seismic Rehabilitation of Existing Buildings (ASCE/SEI 41-06). American Society of Civil Engineers, Reston, VA, US. Baez, J. I., Miranda, E. (2000) Amplification Factors to Estimate Inelastic Displacement Demands for the Design of Structures in the Near Field, Proceedings of 12th World Conference on Earthquake Engineering, Auckland, NZ, 1561. Baker, J. W., Cornell, C. A. (2006) Which Spectral Acceleration Are You Using? Earthquake Spectra. 22(2), 293- 312. Baker, J. W. (2007) Quantitative Classification of Near-Fault Ground Motions Using Wavelet Analysis. Bulletin of the Seismological Society of America. 97(5), 1486-1501. Benito-Oterino, B., Rivas-Medina, A., Gaspar-Escribano, J. M., Murphy, P. (2012) El terremoto de Lorca (2011) en el contexto de la peligrosidad y el riesgo sísmico en Murcia, Física de la Tierra, 24, pp. 255-287. (in Spanish). Bertero, V. V., Mahin, S. A., & Herrera, R. A. (1978) Aseismic design implications of near‐fault San Fernando earthquake records. Earthquake Engineering & Structural Dynamics, 6(1), 31-42. Blazquez-Martinez, R., Rodriguez-Estrella, T. Ibargüen-Soler, J., Martinez-Pagan, P. Navarro-Bernal, M., Trigueros-Tornero, E., Alhama-Manteca, I., Tomas-Espin, A., Parra-Costa, C., Garcia-Ayllon, S., Alcaraz- Aparicio, M., Hausen-Vargas, M. (2014). Estudio de Microzonacion Sismica de Lorca. Universidad Politecnica de Cartagena – Consejeria de Fomento, Obras Publicas y Ordenacion del Territorio (Comunidad Autonoma Region de Murcia). Accessed Nov 2017. (in Spanish). Boore, D. M., Watson-Lamprey, J., Abrahamson, N. A., (2006) Orientation-Independent Measures of Ground Motion, Bulletin of the Seismological Society of America, 96(4A), pp. 1502-1511. Boore, D. M., Atkinson, G. M., (2008) Ground-Motion Prediction Equations for the Average Horizontal Component of PGA, PGV, and 5%-Damped PSA at Spectral Periods between 0.01 s and 10.0 s, Earthquake Spectra, 24(1), pp. 99-138. Boore, D., (2010) Orientation-Independent, Nongeometric-Mean Measures of Seismic Intensity from Two Horizontal Components of Motion, Bulletin of the Seismological Society of America, 100(4), pp. 1830-1835. Bradley, B. A., Quigley, M. C., Van Dissen, R. J., Litchfield, N. J., (2014) Ground Motion and Seismic Source Aspects of the Canterbury Earthquake Sequence. Earthquake Spectra. 30(1), pp. 1-15. Cabañas-Rodriguez, L., Carreño-Herrero, E., Izquierdo-Alvarez, A., Martinez-Solares, J. M., Capote-Villar, R., Martinez-Diaz, J., Benito-Oterino, B., Gaspar-Escribano, J., Garcia-Mayordomo, J., Perez-Lopez, R., Rodriguez.Pascua, M. A., Murphy-Corella, P., (2011) Informe del Sismo de Lorca del 11 de mayo de 2011. Instituto Geografico Nacional. Ministerio de Fomento. Madrid, Spain. Accessed Nov 2016. (in Spanish). Cabañas, L. Alcalde, J. M., Carreño, E. Bravo, J. B. (2013). Characteristics of observed strong ground motion accelerograms from the 2011 Lorca (Spain) Earthquake. Bulletin of Earthquake Engineering. 12:1909-1932. doi: 10.1007/s10518-013-9501-0 Chioccarelli, E., Iervolino, I., (2010) Near-source seismic demand and pulse-like records: A discussion for L’Aquila earthquake, Earthquake Engineering & Structural Dynamics, 39(2), pp. 1039-1062. Gaspar-Escribano, J. M., Benito, B., & García-Mayordomo, J. (2008) Hazard-consistent response spectra in the (Southeast Spain): comparison to earthquake-resistant provisions. Bulletin of Earthquake Engineering, 6(2), pp. 179-196. García-Mayordomo, J., Gaspar-Escribano, J. M., & Benito, B. (2007) Seismic hazard assessment of the Province of Murcia (SE Spain): analysis of source contribution to hazard. Journal of Seismology, 11(4), pp. 453-471. González, P. J., Tiampo, K. F., Palano, M. Cannavó, F. Fernández, J. (2012) The 2011 Lorca earthquake slip distribution controlled by groundwater cristal unloading. Nature Geoscience. 5. 821–825. Hall, J. J., Heaton, T. H., Halling, M. W., Wald, D. J. (1995) Near-source ground-motion and its effects in flexible buildings. Earthquake Spectra. 13(4). 569-605. Iwan. W. D. (1997) Drift spectrum: Measure of demand for earthquake ground motions. Journal of Structural Engineering. 132(4). 397-404.

11

Iwan, W. D., Huang, C.T., Guyader, A. C. (1998) Evaluation of the effects of near-source ground motions. Final Report. California Institute of Technology. Pasadena, CA. Iervolino, I., Chioccarelli, E., Baltzopoulos, G. (2012) Inelastic displacement ratio of near source pulselike ground motions. Earthquake Engineering & Structural Dynamics, 41, pp. 2351-2357. Lopez-Comino, J. A., Mancilla, F., Morales, J., Stich, D., (2012) Rupture directivity of the 2011, Mw 5.2 Lorca earthquake (Spain). Geophysical Research Letters, 39(3), L03301. doi:10.1029/2011GL050498. Makris, N., & Black, C. J. (2004). Evaluation of peak ground velocity as a “good” intensity measure for near- source ground motions. Journal of Engineering Mechanics, 130(9), 1032-1044 Martínez-Díaz, J. J., Álvarez-Gómez, J. A., García-Mayordomo, J., Insua-Arévalo, J. M., Martín-González, F., & Rodríguez-Peces, M. J. (2012) Interpretación tectónica de la fuente del terremoto de Lorca de 2011 (Mw 5, 2) y sus efectos superficiales. Boletín Geológico y Minero, 123(4), 441-458. (in Spanish). Martínez-Díaz, J. J., Bejar-Pizarro, M., Álvarez-Gómez, J. A., Mancilla, F. D. L., Stich, D., Herrera, G., Morales, J. (2012) Tectonic and seismic implications of an intersegment rupture: The damaging May 11th 2011 Mw 5.2 Lorca, Spain, earthquake. Tectonophysics, 546, 28-37. Martinez-Solares, J.M. Cantavella-Nadal, J.V. Cabañas-Rodríguez, L. Valero-Zornosa. J.F. (2012) El terremoto de Lorca de 11 de mayo de 2011 y la sismicidad de la región. Física de la Tierra, 24, pp. 17-40. (in Spanish). Mavroeidis, G. P., Papageorgiou, A. S., (2003), A Mathematical Representation of Near-Fault Ground Motions. Bulletin of the Seismological Society of America, 93(3), 1099-1131. Moratto L., Saraò A., Vuan A., Mucciarelli M., Jimenez M., Garcia Fernandez M. (2017). The 2011 MW 5.2 Lorca earthquake as a case study to investigate the ground motion variability related to the source model. Bulletin of Earthquake Engineering. 15, 3463-3482. Pro, C. Buforn, E. Cesca, S. Sanz de Galdeano, C. Udías, A. (2014) Rupture process of the Lorca (southeast Spain) 11 May 2011 (Mw=5.1) earthquake. Journal of Seismology, 18, 481-495. Rueda-Nuñez, J., Mezcua-Rodríguez, J. García-Blanco, R. M., (2012) Directividad de la fuente sísmica en el terremoto de Lorca del 11 de mayo de 2011. Física de la Tierra, 24, 83-111 (in Spanish). Rueda, J., Dreger, D., Blanco, R. M. G., Mezcua, J., (2014) Directivity Detection and Source Properties of the 11 May 2011 Mw 5.2 Lorca, Spain, Earthquake. Bulletin of the Seismological Society of America, 104(4), 1735-1749. Ruiz-García, J., Miranda, E., (2003) Inelastic displacement ratios for evaluation of existing structures. Earthquake Engineering & Structural Dynamics, 32(8), 1237-1258. Ruiz-García, J., & Miranda, E., (2005) Performance-based assessment of existing structures accounting for residual displacements, Report No. 153, John A. Blume Earthquake Center, Stanford University, https://purl.stanford.edu/tz123dt3431. Ruiz-García, J., (2011) Inelastic displacement ratios for seismic assessment of structures subjected to forward- directivity near-fault ground motions. Journal of Earthquake Engineering, 15(3), 449-468. Shahi, S. K., Baker, J. W., 2011, An Empirically Calibrated Framework for Including the Effects of Near-Fault Directivity in Probabilistic Seismic Hazard Analysis, Bulletin of the Seismological Society of America, 101(2), pp. 742-755. Shahi, S. K., Baker, J. W., (2013) A Probabilistic Framework to Include the Effects of Near-fault Directivity in Seismic Hazard Assessment, Report No. 180, The John A. Blume Earthquake Engineering Research Center, Stanford University. Shahi, S. K., Baker, J. W., (2014) NGA-West2 models for ground-motion directionality, Earthquake Spectra, 30(3), pp. 1285-1300. Shakal, A., Graizer, V., Huang, M., Borcherdt, R., Haddadi, H., Lin, K-W., Stephens, C., Roffers, P., (2005) Preliminary analysis of strong-motion recordings from the 28 September 2004 Parkfield, California earthquake. Seismological Research Letters. 76(1). pp. 27-39. Somerville, P. G., Smith, N. F., Graves, R. W., Abrahamson, N. A., (1997) Modification of Empirical Strong Ground Motion Attenuation Relations to Include the Amplitude and Duration Effects of Rupture Directivity. Seismological Research Letters. 68(1). pp. 199-222.

12