Research article 10 May 2021
Research article  10 May 2021
Reliabilitybased strength modification factor for seismic design spectra considering structural degradation
 ^{1}Instituto de Ingenieria, Universidad Nacional Autonoma de Mexico, Ciudad Universitaria, Coyoacan, C. P. 04510, Mexico City, CDMX, Mexico
 ^{2}Facultad de Ingeniería, Universidad Autonoma de Sinaloa, Calzada de las Americas y B. Universitarios s/n, C. P. 80040 Culiacán, Sinaloa, Mexico
 ^{1}Instituto de Ingenieria, Universidad Nacional Autonoma de Mexico, Ciudad Universitaria, Coyoacan, C. P. 04510, Mexico City, CDMX, Mexico
 ^{2}Facultad de Ingeniería, Universidad Autonoma de Sinaloa, Calzada de las Americas y B. Universitarios s/n, C. P. 80040 Culiacán, Sinaloa, Mexico
Correspondence: Sonia E. Ruiz (sruizg@iingen.unam.mx) and Edén Bojórquez (eden@uas.edu.mx)
Hide author detailsCorrespondence: Sonia E. Ruiz (sruizg@iingen.unam.mx) and Edén Bojórquez (eden@uas.edu.mx)
For earthquakeresistant design, structural degradation is considered using traditional strength modification factors, which are obtained via the ratio of the nonlinear seismic response of degrading and nondegrading structural singledegreeoffreedom (SDOF) systems. In this paper, with the aim to avoid the nonlinear seismic response to compute strength modification factors, a methodology based on probabilistic seismic hazard analyses (PSHAs), is proposed in order to obtain strength modification factors of design spectra which consider structural degradation through the spectralshape intensity measure ${I}_{{N}_{\mathrm{p}}}$. PSHAs using ${I}_{{N}_{\mathrm{p}}}$ to account for structural degradation and Sa(T_{1}), which represents the spectral acceleration associated with the fundamental period and does not consider such degradation, are performed. The ratio of the uniform hazard spectra in terms of ${I}_{{N}_{\mathrm{p}}}$ and Sa(T_{1}), which represent the response of degrading and nondegrading systems, provides new strength modification factors without the need to develop nonlinear time history analysis. A mathematical expression is fitted to the ratios that correspond to systems located in different soil types. The expression is validated by comparing the results with those derived from nonlinear time history analyses of structural systems.
Structures subjected to cyclic loading induced by intense ground motions can exhibit stiffness and/or strength degradation due to the inelastic nonlinear behavior of their structural elements, which can give place to lengthening of the structural fundamental vibration period T_{1}. The effect of such lengthening can be beneficial for structures whose fundamental period is in the descendant branch of the acceleration response spectrum, and their higher vibration modes have little influence on the structural response. Conversely, the effect can be detrimental for structures whose vibration period is in the ascendant branch of the response spectrum. In the latter case, the effect of “structural softening” can have severe consequences because the structure undergoes seismic loading greater than that assumed for its design (Akkar et al., 2004; Chenouda and Ashraf, 2008; Chopra and Chintanapakdee, 2004; TeránGilmore and Espinosa Johnson, 2008). For example, during the GuerreroMichoacán 19 September 1985 Mexican earthquake, many midrise buildings (5 to 10story buildings) with T_{1}=0.7–1.2 s approximately, located in soft soil of Mexico City, which has a vibration period around 2 s, suffered severe structural damage (including collapse) because of the degrading structural effect (Montiel and Ruiz, 2007).
Seismic design guidelines for building structures recommend modifying the responsespectra ordinates by a series of factors in order to include relevant structural behavior that affects the structural response. Those factors are related, for example, to seismic behavior, structural overstrength, structural irregularity, degrading behavior, etc. A common practice to derive those modification factors is by means of the ratio between specific response spectra of singledegreeoffreedom (SDOF) systems. Indeed, most current seismic code provisions implement simplified analyses based on these ratios. For example, the Federal Emergency Management Agency (FEMA) introduced the socalled coefficient method (FEMA273, 1997; FEMA356, 2000), which consists of multiplying the elastic design spectrum by several coefficients. One of them takes into account the hysteretic structuredegrading behavior. More recently, FEMA440 (2005) presented some improvements to current nonlinear analysis procedures. Accordingly, the coefficient method suffered slight adjustments, where the coefficient that incorporates the effect of degrading structural behavior was updated. At present, the simplified nonlinear approach is available in FEMA P58 (2012) methodology. Another example is the Manual for Civil Structures Design (MCSD, 2015), developed by the Federal Commission of Electricity of Mexico, which specifies a degrading factor that increases or decreases the design spectral ordinates, due to structural deterioration.
The hysteretic degrading behavior is particularly severe for structures located in soft soil, like that in the lake bed zone of Mexico City, where there is a highdensity population, and the site effects make it susceptible to severe earthquake damage (Singh et al., 1988, 2018). In spite of that, the current Mexico City Building Code (MCBC, 2017) does not specify any structuredegrading factor.
This study is aiming to propose a methodology for obtaining a mathematical expression corresponding to a structuredegrading factor for seismic design of buildings that exhibit period lengthening. The expression is a function of both the structural period and the dominant period of the soil. The methodology can be applied to any highseismichazard region of the world. Finally, notice that the variation in the vibration periods of a structure from the undamaged to the damaged state strongly depends of several parameters, and this is crucial to consider different design limit states. Although the procedure is not affected by these parameters, the variation in the structural period could be taken into account considering different values of T_{N} (see definition of N_{p} below); however, the assessment of this value accounts for the design limit state, structural type, interaction of the structural elements with the nonstructural ones requiring the study of specific structural systems such as reinforced concrete, momentresisting steel frames, masonry, structures with eccentrically buckling restrained braces, and posttensioned based isolators, which are out of the scope of the present study. On the other hand, soil–structure interaction (SSI) was not taken into account to compute the structuredegrading modification factors for seismic design spectra; nevertheless, notice that the effect of SSI is more important for stiff structures located on soft soil, in such a way that for this type of structure, the ordinates of the response spectra tend to increase while the opposite occurs for flexible structures (Avilés and PérezRocha, 2007). The results obtained in the present study could be modified to include the effect of SSI via the current Mexico City Building Code (MCBC, 2017) which provides recommendations about this issue.
First, it is necessary to perform probabilistic seismic hazard analyses (PSHAs) corresponding to a firmground site and then softsoil sites located in the seismic area of interest. PSHAs are associated with Sa(T_{1}) and alternatively with ${I}_{{N}_{\mathrm{p}}}$ intensity measures, where Sa(T_{1}) represents the spectral acceleration at the fundamental period of a structure, and I_{Np }is an intensity measure that accounts for the period lengthening due to structuredegrading behavior (${I}_{{N}_{\mathrm{p}}}$ is defined below). Although Sa(T_{1}) is the most used ground motion intensity measure (IM) around the world for PSHAs, it has some limitations. For example, it does not consider the effect of period lengthening of the structure due to its nonlinear behavior and mechanical property degradation (Baker and Cornell, 2005; Bojórquez et al., 2008, 2017a; Bojórquez and Iervolino, 2011; Kostinakis et al., 2018; Cordova et al., 2001; Shome et al., 1998; Tothong and Luco, 2007).
Second, uniform hazard spectra (UHS) of ${I}_{{N}_{\mathrm{p}}}$ and Sa(T_{1}), which represent the response of degrading and nondegrading systems, respectively, are obtained. The UHS are computed for several seismic recording stations located in different soil conditions. Subsequently, the effect of the structural degradation on the response of SDOF systems is characterized by the ratio between the uniform hazard spectra: ${I}_{{N}_{\mathrm{p}}}/\mathrm{Sa}\left({T}_{\mathrm{1}}\right)$.
Finally, a mathematical expression is adjusted to the spectral ratios. In order to verify that the mathematical expression leads to reasonable results, it is convenient to compare these with those obtained with other expressions found in the literature.
In what follows, a description of the methodology is presented (see Fig. 1).

First, PSHAs are carried out for the firmground site of interest, corresponding to Sa(T_{1}) and, alternatively, to ${I}_{{N}_{\mathrm{p}}}$. With the purpose of performing the analyses, the seismic tectonic zones that contribute to the seismic hazard of the site are identified.

Then, the probability distribution for earthquake magnitude and sourcetosite distance are assumed. Additionally, it is necessary to define adequate ground motion prediction equations (GMPEs).

With the total probability theorem and the information previously defined, the mean annual rates of exceedance (seismic hazard curves) corresponding to the site located in firm ground are obtained.

Once the hazard curves for firm ground are available, the mean annual rates of exceedance of seismic recording stations located in different soil types of the seismic area of interest are estimated (using a technique described in the following sections). The stations are grouped in different zones, which depend on the dominant period of the soil, T_{s}.

For each recording station site, UHS associated with a given return period are computed for Sa(T_{1}), and alternatively, for ${I}_{{N}_{\mathrm{p}}}$.

Next, the spectral ratio ${I}_{{N}_{\mathrm{p}}}/\mathrm{Sa}\left({T}_{\mathrm{1}}\right)$ is estimated for each site. ${I}_{{N}_{\mathrm{p}}}/\mathrm{Sa}\left({T}_{\mathrm{1}}\right)$ represents the ratio of strength demands between systems with degrading and systems with nondegrading structural behavior.

Finally, a simplified mathematical expression is adjusted to the spectral ratio ${I}_{{N}_{\mathrm{p}}}/\mathrm{Sa}\left({T}_{\mathrm{1}}\right)$. The expression contains parameters that depend on the zone of interest.

The results of the expression proposed are compared with those obtained from other expressions found in the literature, which were obtained from time history analyses.
For illustrative purposes, in the following sections, the methodology proposed above is applied in order to find mathematical expressions of structuredegrading factors of the design spectra specified in MCBC; however, the approach can be applied to any seismic region in the world.
3.1 Earthquake sources
The evaluation of a probabilistic seismic hazard analysis for a particular site requires identification of all possible earthquake sources capable of producing a significant seismic event. For this purpose, Zúñiga et al. (2017) proposed a seismic regionalization for Mexico, which is used in the present study. Figure 2a shows the shallowdepth seismic zones where interplate earthquakes occur due to the subduction of the Rivera and Cocos plates (SUB1–SUB4). Figure 2b illustrates the intermediatedepth seismic zones. This region corresponds to intraslab events that take place inside the subducted Rivera and Cocos plates below southcentral Mexico (IN1 to IN3). Additionally, Fig. 2c displays the seismic zones for characteristic seismic events (C1 to C14) proposed by Ordaz and Reyes (1999). Seismic zones in Fig. 2c are also included in the present study to compute PSHA.
3.2 Magnitude probability distribution
Seismic sources are capable of producing different earthquake sizes. Therefore, it is crucial to define the probability distribution of the earthquake magnitudes and corresponding rates of occurrence for each source. In this regard, the distribution of earthquake sizes is commonly described by the bounded Gutenberg–Richter recurrence law (Eq. 1).
where λ_{m} is the mean annual rate of exceedance for earthquakes between a minimum magnitude M_{min} and a maximum magnitude M_{max}, and $\mathit{\nu}=\mathrm{exp}(\mathit{\alpha}\mathit{\beta}{M}_{\mathrm{min}})$ is the mean annual number of earthquakes of magnitude M_{w}≥M_{min}, where α=2.303p and β=2.303q. The values of p and q are indicated in Fig. 2a and b, according to Zúñiga et al. (2017).
For the seismic sources related to characteristic earthquakes (Fig. 2c), the bounded Gutenberg–Richter recurrence law does not accurately describe the magnitude exceedance rates. Accordingly, for M_{w}>7, we employ a Gaussian probability distribution function (pdf) of magnitudes to account for the characteristic events in the Mexican subduction zones (see Eq. 2) (Ordaz and Reyes, 1999).
where ν_{7} is the mean annual number of earthquakes of magnitude M_{w}>7; ${E}_{{M}_{\mathrm{w}}}$ and ${\mathit{\sigma}}_{{M}_{\mathrm{w}}}$ are the mean and standard deviation of the magnitude, respectively, and Φ (.) is the normal distribution function. The corresponding parameters to evaluate the distribution are shown in Fig. 2c.
The present study assumes M_{min}=4.5 and M_{max}=6.9 for the interplate shallowdepth seismic zones SUB1–SUB3 (see Fig. 2a). In contrast, M_{min}=4.5 and M_{max}=7.2, 7.8 and 7.9 are assumed for IN1–IN3, respectively (intermediatedepth seismic zones, Fig. 2b). Finally, M_{min}=7.0 and M_{max}=8.1 are assumed for the 14 earthquake sources shown in Fig. 2c.
3.3 Sourcetodistance distribution
Once the earthquake magnitude distribution is established, the pdf of distances from the earthquake location to the site of interest must be characterized. A uniform pdf is generally assigned to any point in the seismic zone (McGuire, 1995; Kramer, 1996). Since the area sources, where earthquakes can occur, are welldelimited (Fig. 2a–c), it is straightforward to determine the sourcetodistance distribution.
3.4 Ground motion prediction equations
Attenuation relationships are fundamental for PSHA. They are commonly developed to predict the peak ground acceleration, PGA, or the spectral acceleration, Sa(T_{1}). Unfortunately, attenuation models have not yet been devised to provide ${I}_{{N}_{\mathrm{p}}}$ as a function of the vibration period (as is done with existing GMPEs); however, with GMPEs for Sa(T_{1}) currently available, it is possible to perform PSHA using ${I}_{{N}_{\mathrm{p}}}$. Here we employ the GMPEs proposed by Reyes et al. (2002) and Jaimes et al. (2015) for interplate and intraslab events, respectively. They were developed using accelerometric data recorded at Ciudad Universitaria station (CU), which is located at the hill zone (firm ground) of Mexico City, basically conformed by a surface layer of lava flows and volcanic tuffs with a shear wave velocity in the upper 30 m of 750 m s^{−1} (Ordaz and Singh, 1992; Singh et al., 2018).
3.5 Seismic hazard curves
The final product of a PSHA can be expressed in different forms. Seismic hazard curves are used frequently to represent the seismic hazard. They indicate the annual rate of exceeding a variety of intensity levels of a ground motion parameter at a site of interest. The procedure to compute a ground motion hazard curve is based on the total probability theorem (Baker, 2013; Cornell, 1968; Esteva, 1968; McGuire, 1995; Kramer, 1996).
In order to overcome the limitations of traditional IMs (e.g., PGA, Sa(T_{1})), advanced seismic IMs have been proposed. Some researchers suggest using vectorvalued ground motion IMs. By including two or more representative parameters of the ground motion, accurate evaluations of seismic performance can be achieved (Baker and Cornell, 2005; Bojórquez et al., 2008, 2017a; Bojórquez and Iervolino, 2011; Kostinakis et al., 2018; Cordova et al., 2001; Tothong and Luco, 2007). Accordingly, Bojórquez et al. (2008) developed the vectorvalued intensity measure <Sa(T_{1}), N_{p}>, where N_{p} is a parameter proxy for the spectral shape. This IM is an advancement in predicting the seismic response in comparison with other IMs. However, the evaluation of PSHA using vectorvalued IMs is a complicated and impractical task; therefore, Bojórquez and Iervolino (2011) introduced a scalar IM based on Sa(T_{1}) and N_{p}, called ${I}_{{N}_{\mathrm{p}}}$; both N_{p} scalar and vectorvalued intensity measures have been effectively used (Bojórquez et al., 2012, 2017b).
Accordingly, Buratti (2012) made an exhaustive comparison of the most influential scalar IMs available in the literature with respect to efficiency and sufficiency. The study concluded that the most effective intensity measure was ${I}_{{N}_{\mathrm{p}}}$. Additionally, De Biasio et al. (2014), based on a comparative study of structures with nonlinear behavior, showed the good performance of ${I}_{{N}_{\mathrm{p}}}$ to predict maximum interstory drift and maximum ductility demands. Moreover, Modica and Stafford (2014) using <Sa(T_{1}) and N_{p}>, estimated the fragility and performance of buildings with higher efficiency with respect to different IMs. In this context, Minas and Galasso (2019) showed the advantages of ${I}_{{N}_{\mathrm{p}}}$ comparing Sa(T_{1}) fragility curves, for different damage states. Additionally, Yakhchalian et al. (2015) demonstrated the efficiency of the parameter N_{p}. They showed that the parameter N_{p} works appropriately, particularly in performance levels related to moderate levels of nonlinearity. Similarly, Kostinakis and Athanatopoulou (2016) proved the adequate efficiency of ${I}_{{N}_{\mathrm{p}}}$ to reduce the uncertainty in the prediction of the response of reinforced concrete buildings. In addition, Jamshidiha et al. (2018) examined the ability of different IMs to predict the seismic collapse capacity of steel momentresisting frames with fluid viscous dampers. They concluded that the scalar IM that resulted from the combination with the parameter N_{p} was most efficient.
Based on the literature mentioned above, the authors of the present study concluded that ${I}_{{N}_{\mathrm{p}}}$ is a promising tool to perform PSHA.
4.1 Methodology to perform a PSHA using ${I}_{{N}_{\mathrm{p}}}$
In this section a methodology to perform PSHA using ${I}_{{N}_{\mathrm{p}}}$ is proposed. First, ${I}_{{N}_{\mathrm{p}}}$ is defined as follows (Bojórquez and Iervolino, 2011):
where ${I}_{{N}_{\mathrm{p}}}$ is the scalar intensity measure, α is a parameter that should be calibrated according to the structure and the earthquake demand parameter selected (in this study α=0.5 is adopted, as recommended in Bojórquez and Iervolino, 2011), and Sa_{avg}(T_{1} … T_{N}) is the geometric mean of the spectral acceleration at N numbers of structural vibration periods considered. Sa_{avg}(T_{1} …T_{N}) takes into account the vibration period lengthening due to structural damage and is expressed as
Substituting Eqs. (4) and (5) into Eq. (3), applying the natural logarithm, results in
Then, the expected value and the variance of $\mathrm{ln}\left({I}_{{N}_{\mathrm{p}}}\right)$ can be expressed as in Eqs. (7) and (8), respectively.
The values of ln [Sa(T_{i})] are obtained from existing attenuation models (e.g., the GMPEs described in Sect. 3.4). On the other hand, ln [Sa(T_{i})] terms are commonly assumed to have a joint Gaussian pdf; consequently, the summation also has Gaussian distribution. Therefore, the variance Var{ln [Sa_{avg}(T_{1} … T_{N})]} and the correlation coefficient ρln [Sa_{avg}(T_{1} … T_{N})]ln [Sa(T_{i})] can be obtained by Eqs. (9) and (10), respectively:
where the term ρln [Sa(T_{i})], ln [Sa(T_{j})] represents the correlation between spectral acceleration values at periods T_{i} and T_{j}. The correlation coefficients have been obtained by the authors of the present study (RodríguezCastellanos et al., 2019, 2020).
4.2 Values of T_{N}
Among the parameters that define the intensity measure ${I}_{{N}_{\mathrm{p}}}$, the geometric mean, Sa_{avg}(T_{1} … T_{N}), has a crucial role when computing the uniform hazard spectra (UHS). The T_{N} value (Nth structural vibration period) takes into account the level of nonlinearity developed by the structure. Bojórquez et al. (2008) and Bojórquez and Iervolino (2011) recommend using T_{N}=2.0T1. Nevertheless, we consider that there is no optimal period range for Sa_{avg}(T_{1} … T_{N}) that meets the entire range of structural vibration periods; therefore, here we propose that T_{N} should depend on the structural vibration period, which is in agreement with Tsantaki et al. (2012, 2017).
It has been pointed out that the stiffer the structure, the larger the period lengthening. Accordingly, for structures with short vibration periods, we adopt T_{N}=2.0T_{1}, which agrees with recommendations made by Bianchini et al. (2009), Katsanos and Sextos (2015), and Tsantaki et al. (2017), for relatively stiff structures, and assuming a ductility demand between 2 and 3.
At short to moderate vibration periods, the structural period lengthening diminishes somewhat linearly until it reaches a semiconstant behavior (which is independent of the level of nonlinearity developed by the structure) (Katsanos and Sextos, 2015). In this regard, Di Sarno and Amiri (2019) quantified the fundamental period lengthening of structures by the ratio of response spectra corresponding to the lengthened and the elastic structural vibration period (${T}_{\mathrm{in}}/{T}_{\mathrm{el}}$). They suggested dividing the response spectra into two main regions: the first associated with short to moderate period structures, whose period shift ratio ${T}_{\mathrm{in}}/{T}_{\mathrm{el}}$ decreases with increasing the elastic period, and the second region related to longperiod structures, where the ratio period ${T}_{\mathrm{in}}/{T}_{\mathrm{el}}$ behaves practically constant. Consequently, there must be a certain bound where the period shift ratio switches to remain constant; therefore, we propose T_{N}=T_{s} as that bound from which the lengthening of the structural vibration period remains almost constant. In this context, Miranda and RuizGarcia (2002), RuizGarcia and Miranda (2003), and independently TeránGilmore and Espinosa Johnson (2008), found that strength demands between degrading and nondegrading systems are similar when the structural period and dominant soil period are comparable, which means that the mean ratio value should be approximate to one when T_{n}≈T_{s}.
For vibration periods longer than the dominant soil period, it is assumed T_{N}=1.25T_{1}, which is, on average, the period shift ratio value for structures with a short to moderate nonlinearity level, that is, with a ductility ratio around 2 to 3 (Katsanos and Sextos, 2015; Di Sarno and Amiri, 2019).
Summarizing, we used T_{N}=2.0T_{1} in this study for structural systems with a short fundamental period; T_{N}=T_{s} for those with an intermediate period and T_{N}=1.25T_{1} for systems with a long fundamental period. It is possible to get a better approximation of T_{N} bounds, by means of a parametric study of the ratios of the equivalent period of SDOF degraded systems and that of the elastic systems (${T}_{\mathrm{in}}/{T}_{\mathrm{el}}$), as a function of T_{el}, for a given ductility. Such a study can consider both ground motion characteristics and structural properties (such as degrading stiffness ratio, pinching factor, accumulated damage factor, etc.), as was done by Di Sarno and Amiri (2019). They proposed a mathematical expression for estimating the lengthening of the fundamental period as a function of the structural elastic period and the significant structural parameters, which is applicable to systems in site classes D and C according to ASCE 7–10 (2010), with shear wave velocities $\mathrm{182.88}<{V}_{{\mathrm{s}}_{\mathrm{30}}}<\mathrm{365.76}$ and $\mathrm{365.76}<{V}_{{\mathrm{S}}_{\mathrm{30}}}<\mathrm{762}$ m s^{−1}, respectively. However, the T_{N} bounds used here lead to reasonable results, as is verified below.
5.1 Uniform hazard spectra corresponding to firm ground
The uniform hazard spectra are computed, first, for the CU site, which is in firm ground. Figure 3a shows the UHS if only interplate, or alternatively intraslab, earthquakes occur. It also displays when both types of events are considered simultaneously (Total). Figure 3b shows the total UHS of Sa(T_{1}) and ${I}_{{N}_{\mathrm{p}}}$, both associated with a 250year return period. It can be seen that the spectra are quite similar; practically, they reach the same acceleration levels, and slight differences occur at long periods.
5.2 Uniform hazard spectra corresponding to softsoil sites
Estimating the seismic hazard at firm ground allows us to proceed with a technique to assess the seismic hazard at softsoil sites. In this regard, Esteva (1970) presented a formulation in which through a known hazard curve at a reference site it is feasible to estimate a hazard curve at a recipient site. In this study, we used CU station as the reference site because, since 1964, it has recorded all the significant ground motions that have struck Mexico City. In addition, different studies have taken CU as a reference site (Ordaz et al., 1988; Reinoso and Ordaz, 1999; Singh et al., 1988). Therefore, it is viable to perform a hazard analysis for CU station and then to compute the annual rate of exceedance at other sites located in soft or medium soils, as follows:
where ν_{Y}(y) is the mean annual rate of exceedance of a seismic IM, for the recipient site. ${\mathit{\nu}}_{x}(y/\mathit{\tau})$ is the mean annual rate of exceedance of a seismic IM for the reference site, divided by the variable τ. τ represents the response spectral ratios between the response spectra corresponding to the recipient site and the reference site ($Y/X$). fτ(τ) is the pdf of τ.
Therefore, to evaluate the previous function, firstly, the spectral ratios are estimated, and then they are coupled with the seismic hazard curves via Eq. (11). In this respect, Fig. 4a–f show the mean response of the spectral ratios for Sa(T_{1}) (solid line) and ${I}_{{N}_{\mathrm{p}}}$ (dashed line) for one representative station located in each of the zones listed in Table 1. In this sense, the spectral ratios roughly represent the spectral amplification of soft soil with respect to firm ground. It is observed how the peak values shift towards increasingly longer periods, which, approximately, match with the dominant soil period (see Table 1). For this analysis, more than 1100 ground motion records corresponding to the different recording stations were used. The stations are grouped depending on the dominant soil period where these are located, as follows: zone A: T_{s}<0.5 s; zone B: 0.5 s < T_{s} < 1.0 s; zone C: 1.0 s < T_{s} < 1.5 s; zone D: 1.5 s < T_{s} < 2.0 s; zone E: 2.0 s < T_{s} < 2.5 s; and zone F: 2.5 s < T_{s} < 3.0 s. Additionally, Fig. 5 shows the location of the recording stations in Mexico City, which are represented with circles of different colors associated with each of the proposed zones (see Table 1).
Next, in order to compute the mean annual rate of exceedance of Sa(T_{1}) and ${I}_{{N}_{\mathrm{p}}}$, the seismic hazard curves corresponding to CU station are coupled with the response spectral ratios, using Eq. (11). Figure 6a to f show the hazard curves (λ) of Sa(T_{1}) and ${I}_{{N}_{\mathrm{p}}}$, associated with different vibration periods, corresponding to CU and the same recording stations of Fig. 4a to f. First, as expected, the rates of exceedance for all the recording stations analyzed are higher than the corresponding ones of CU (up and down). Additionally, concerning the CU site, the hazard curves for both intensity measures ${I}_{{N}_{\mathrm{p}}}$ and Sa(T_{1}) are very similar, and differences are barely visible at long return periods. Now, for the rest of the recording stations, Fig. 4c and d show noticeable variations between exceedance rates of Sa(T_{1}) and ${I}_{{N}_{\mathrm{p}}}$; nevertheless, Fig. 4e and f display almost no contrast between the rates of exceedance of the two intensity measures. The previous is relative, because to fully characterize the variations between exceedance rates of Sa(T_{1}) and ${I}_{{N}_{\mathrm{p}}}$, a wide range of periods need to be covered; for this reason, we estimate the UHS in the following.
Then, having the mean rates of exceedance for each recording station site (see Table 1), the UHS are estimated for a given return interval. Figure 7a to f show the UHS of Sa(T_{1}) and I_{Np} for the same stations of Fig. 6a to f, for a 250year return period. It is observed that, at vibration periods shorter than the dominant soil period, the spectral ordinates corresponding to firm ground (zones A–C) are comparable for both IMs. However, at soft soil (zones D–F), the spectral ordinates of ${I}_{{N}_{\mathrm{p}}}$ are notably higher than those of Sa(T_{1}) (up to 30 %). In contrast, at vibration periods longer than T_{s}, they are smaller than those corresponding to Sa(T_{1}) (5 % to 20 %, depending on the soil type). The same can be appreciated for different sites of the city in the maps shown in Fig. 8a to d, which corresponds to Sa(T_{1}) (left side) and ${I}_{{N}_{\mathrm{p}}}$ (right side), for T_{1}=0.5 s (up side) and T_{1}=1.0 s (down side), for a return interval of T_{r}=250 years.
Once the uniform hazard spectra of Sa(T_{1}) and ${I}_{{N}_{\mathrm{p}}}$ were estimated, the degrading structural behavior effect is evaluated by means of the ratio ${I}_{{N}_{\mathrm{p}}}/\mathrm{Sa}\left({T}_{\mathrm{1}}\right)$. It represents the ratio of strength demands between a system with degrading and the same system with nondegrading structural behavior. The ratios are obtained for each station of the zones listed in Table 1. Figure 9a to f show the ${I}_{{N}_{\mathrm{p}}}/\mathrm{Sa}\left({T}_{\mathrm{1}}\right)$ ratios (thin gray lines) as a function of the normalized periods ${T}_{n}/{T}_{\mathrm{s}}$, for zones A to F, respectively.
Based on these ratios, the following spectral modification function (SMF) was proposed, which is a variation of that specified by MCSD (2015):
where the values of a–e are shown in Table 2. It is noticed that the values of the parameters depend on the type of soil where the structure is located; conversely, those in the MCSD function are constant values. In addition, such a function is restricted only to soft soils.
Figure 9a to f show the equation proposed here (Eq. 12) (thick dashed line), as well as the MCSD (2015) function (thick solid line). In the figures, the horizontal and vertical dotted lines, aligned at ${I}_{{N}_{\mathrm{p}}}/\mathrm{Sa}\left({T}_{\mathrm{1}}\right)=\mathrm{1}$ and ${T}_{n}/{T}_{\mathrm{s}}=\mathrm{1}$, approximately delimitate the increase or decrease in the spectral amplification.
The figures show the following.
 a.
The highest ${I}_{{N}_{\mathrm{p}}}/\mathrm{Sa}\left({T}_{\mathrm{1}}\right)$ ratios are reached for structures with vibration periods shorter than the dominant soil period (approximately ${T}_{s}/\mathrm{2}$), which indicates that the lateral strength demand for degrading systems is higher than the strength demand for nondegrading systems.
 b.
When the vibration period of the system is close to the dominant soil period (${T}_{n}/{T}_{\mathrm{s}}\approx \mathrm{1}$), the strength demands for degrading and nondegrading systems are similar.
 c.
When ${T}_{n}/{T}_{\mathrm{s}}>\mathrm{1}$, the demands of the degrading systems decrease with respect to those of the nondegrading systems. This means that for structural vibration periods longer than T_{s}, the degrading behavior provides a beneficial effect.
 d.
It is noticed that for zone D (Fig. 9d), the MCSD function predicts spectral modification values which are similar to the function proposed in the present study (Eq. 12). This happens because the MCSD function was calibrated using ground motion data recorded at a station located in that zone (SCT station in zone D); however, it does not happen the same for other soil conditions, especially for ${T}_{n}/{T}_{\mathrm{s}}>\mathrm{1}$.
 e.
Equation (12) predicts values closer to unity at sites in zones A–C (firm ground and transition soil) than at zones D–F, which means that the structural softening is not as significant as it is for zones D–F. In this respect, several studies have observed that the degradation of the stiffness has little effect on the strength demands for structures located on firm sites (Akkar et al., 2004; Chenouda and Ashraf, 2008; Chopra and Chintanapakdee, 2004). Moreover, it is noticed that at very short vibration period systems (${T}_{n}/{T}_{\mathrm{s}}<\mathrm{0.1}$), the SMF proposed here predicts amplification values very close to unity, which is consistent for extremely stiff structures.
 f.
Finally, the reduction of strength demand according to Eq. (12) fits the observed data (thin gray lines) better for each type of soil (zones A to F) than that recommended by MCSD guidelines.
With the aim of verifying the validity of the proposed expression, Fig. 10a and b compare the results of Eq. (12) with those obtained from time history analysis of SDOF systems. The figures show the mean ratio of strength demands of degrading and of nondegrading systems (elastoplastic behavior) corresponding to a ductility value, μ_{u} (thin gray lines), using firmground and softsoil records, respectively (Miranda and RuizGarcia, 2002; TeránGilmore and Espinosa Johnson, 2008). The ground motions at firm ground (Fig.10a) correspond to synthetic accelerograms (T_{s}=1.0 s) (TeránGilmore and Espinosa Johnson, 2008) and ground motions recorded in the San Francisco Bay area during the 1989 the Loma Prieta earthquake (T_{s}≈1.1 s) (Miranda and RuizGarcia, 2002). In contrast, the ground motions in soft soil were recorded in the lake bed zone of Mexico City (T_{s}≈2.0 s) (Fig. 10b).
Figure 10a and b also include the ${I}_{{N}_{\mathrm{p}}}/\mathrm{Sa}\left({T}_{\mathrm{1}}\right)$ ratios, corresponding to the stations D11 and C3, estimated from the uniform hazard spectra normalization (thick red dotted lines). It can be observed that the ${I}_{{N}_{\mathrm{p}}}/\mathrm{Sa}\left({T}_{\mathrm{1}}\right)$ ratio agrees with the results of Miranda and RuizGarcia (2002), and of TeránGilmore and Espinosa Johnson (2008). The figures also show that the function given by Eq. (12) is in agreement with both the observed data obtained from the time history analyses and the ${I}_{{N}_{\mathrm{p}}}/\mathrm{Sa}\left({T}_{\mathrm{1}}\right)$ ratio calculated from the study based on seismic hazard analyses.
A methodology based on probabilistic seismic hazard analysis is proposed to evaluate the effect of degrading behavior on the strength demands of SDOF systems. For this aim, uniform hazard spectra are obtained for two alternative intensity measures: ${I}_{{N}_{\mathrm{p}}}$ and Sa(T_{1}), which represent the response of degrading and nondegrading systems, respectively. Thus, the ratio of the hazard spectra ${I}_{{N}_{\mathrm{p}}}/\mathrm{Sa}\left({T}_{\mathrm{1}}\right)$ characterizes the strength demands of systems with degrading behavior to those of systems with nondegrading behavior. Based on the ${I}_{{N}_{\mathrm{p}}}/\mathrm{Sa}\left({T}_{\mathrm{1}}\right)$ ratios, which correspond to systems located at different sites, grouped in different seismic zones (depending on the type of soil where the structures are located), a mathematical expression is proposed. The methodology is applied here to structural systems located in Mexico City, but it can be applied to any seismic region of the world.
From the study the following is concluded.

For structures with vibration periods shorter than the dominant soil period (${T}_{n}/{T}_{\mathrm{s}}<\mathrm{1}$), degrading systems exhibit strength demands up to 30 % higher than systems with nondegrading behavior.

For structures with vibration periods close to the dominant soil period (${T}_{n}/{T}_{\mathrm{s}}\approx \mathrm{1}$), the strength demands for degrading and nondegrading systems are similar.

For systems with vibration periods longer than the dominant soil period (${T}_{n}/{T}_{\mathrm{s}}>\mathrm{1}$), the strength demands for structures with degrading behavior are lower, approximately 5 % to 20 %, than structures with nondegrading behavior. That reduction highly depends on the dominant soil period at the site, and it is larger for systems with longer dominant soil periods. For these cases, the structuredegrading behavior produces a beneficial effect, reducing the lateral strength requirement of the structures.

A strength modification factor was proposed (Eq. 12). The expression was fitted according to the spectral ratios ${I}_{{N}_{\mathrm{p}}}/\mathrm{Sa}\left({T}_{\mathrm{1}}\right)$ corresponding to different soil conditions. The value of the parameters included in the equation depends on the type of soil where the structure is located.

The expression proposed (Eq. 12) is a useful tool for simplified nonlinear modal analyses, to explicitly incorporate the effect of degrading behavior according to the type of soil where the structure is located. It was verified that the mathematical expression proposed leads to results that are comparable to those obtained from time history analyses of SDOF systems located in soft soil.

In addition, the study presents a methodology to elaborate seismic hazard maps in terms of the intensity measure ${I}_{{N}_{\mathrm{p}}}$. Based on that methodology, the first seismic hazard map of Mexico City is presented, in terms of ${I}_{{N}_{\mathrm{p}}}$.
The data are available for free at http://cires.org.mx/registro_es.php (last access: 30 March 2020) (Centro de Instrumentación y Registro Sísmico, 2020).
ARC developed the theoretical framework, performed the computations and wrote the original manuscript. SER and EB conceived the study and contributed to the development and design of the methodology. MAO and ARS contributed to the sample preparation and review of the final manuscript.
We wish to confirm that there are no known conflicts of interest associated with this publication, and there has been no significant financial support for this work that could have influenced its outcome.
We confirm that the paper has been read and approved by all named authors and that there are no other persons who satisfied the criteria for authorship but are not listed. We further confirm that the order of authors listed in the paper has been approved by all of us. We confirm that we have given due consideration to the protection of intellectual property associated with this work and that there are no impediments to publication, including the timing of publication, with respect to intellectual property. In so doing we confirm that we have followed the regulations of our institution concerning intellectual property.
We understand that the first and the second corresponding authors are the contacts for the editorial process (including editorial manager and direct communications with the office). They are responsible for communicating with the other authors about progress, submissions of revisions and final approval of proofs. We confirm that we have provided current, correct email addresses which are accessible by the corresponding authors and which have been configured to accept email from sruizg@iingen.unam.mx and eden@uas.edu.mx.
Thanks are given to DGAPAUNAM (project PAPIIT IN100320) and to Instituto para la Seguridad de las Construcciones de la Ciudad de Mexico, for their support. The first and fourth authors acknowledge the scholarship given by Consejo Nacional de Ciencia y Tecnología (CONACyT) during their graduate studies. The authors wish to thank CIRES for having provided the seismic records used in this study.
This research has been supported by the Dirección General de Asuntos del Personal Académico, Universidad Nacional Autónoma de México (grant no. PAPIITIN100320).
This paper was edited by Filippos Vallianatos and reviewed by two anonymous referees.
Akkar, S., Yazgan, U., and Gülkan, P.: Deformation limits for simple nondegrading systems subjected to nearfault ground motions, in: 13th World Conference on Earthquake Engineering, Vancouver, Canada, 2004.
ASCE 710: ASCE/SEI 710: minimum design loads for buildings and other structures, Reston, VA, USA, 2010.
Avilés, J. and PérezRocha, L. E.: Damage analysis of structures on elastic foundation, J. Struct. Eng., 133, 1453–1461, https://doi.org/10.1061/(ASCE)07339445(2007)133:10(1453), 2007.
Baker, J. W.: Probabilistic Seismic Hazard Analysis, White Paper Version 2.0.1, 1–79, 2013.
Baker, J. W. and Cornell, C. A.: A vectorvalued ground motion intensity measure consisting of spectral acceleration and epsilon, Earthq. Eng. Struct. Dynam., 34, 1193–1217, https://doi.org/10.1002/eqe.474, 2005.
Bianchini, M., Diotallevi, P. P., and Baker, J. W.: Prediction of inelastic structural response using an average of spectral accelerations, in: 10th International Conference on Structural Safety and Reliability (ICOSSAR09), 13–17 September 2009, Osaka, Japan, 2009.
Bojórquez, E., Iervolino, I., and Manfredi, G.: Evaluating a new proxy for spectral shape to be used as an intensity measure, AIP Conf. Proc., 1020, 1599–1606, https://doi.org/10.1063/1.2963788, 2008.
Bojórquez, E. and Iervolino, I.: Spectral shape proxies and nonlinear structural response, Soil Dynam. Earthq. Eng., 31, 996–1008, https://doi.org/10.1016/j.soildyn.2011.03.006, 2011.
Bojórquez, E., Iervolino, I., ReyesSalazar, A., and Ruiz, S. E.: Comparing vectorvalued intensity measures for fragility analysis of steel frames in the case of narrowband ground motions, Eng. Struct., 45, 472–480, https://doi.org/10.1016/j.engstruct.2012.07.002, 2012.
Bojórquez, E., Chávez, R., ReyesSalazar, A., Ruiz, S. E., and Bojórquez, J.: A new ground motion intensity measure IB, Soil Dynam. Earthq. Eng., 99, 97–107, https://doi.org/10.1016/j.soildyn.2017.05.011, 2017a.
Bojórquez, E., Baca, V., Bojórquez, J., ReyesSalazar, A., Chávez, R., and Barraza, M.: A simplified procedure to estimate peak drift demands for midrise steel and R/C frames under narrowband motions in terms of the spectralshapebased intensity measure INp, Eng. Struct., 150, 334–345, https://doi.org/10.1016/j.engstruct.2017.07.046, 2017b.
Buratti, N.: A comparison of the performances of various groundmotion intensity measures, in: Proc. 15th World Conference on Earthquake Engineering, 24–28 September 2012, Lisbon, Portugal, 2012.
Centro de Instrumentación y Registro Sísmico, A. C.: Búsqueda de registros, available at: http://cires.org.mx/registro_es.php, last access: 30 March 2020.
Chenouda, M. and Ashraf, A.: Inelastic displacement ratios of degrading systems, J. Earthq. Eng., 134, 1030–1045, https://doi.org/10.1061/(ASCE)07339445(2008)134:6(1030), 2008.
Chopra, A. K. and Chintanapakdee, C.: Inelastic deformation ratios for design and evaluation of structures: Singledegreeoffreedom bilinear systems, J. Struct. Eng., 130, 1309–1319, https://doi.org/10.1061/(ASCE)07339445(2004)130:9(1309), 2004.
Cordova, P. P., Deierlein, G. G., Sameh, S. F., and Mehanny, C. A. C.: Development of a twoparameter seismic intensity measure and probabilistic design procedure, J. Eng. Appl. Sci., 51, 233–252, 2001.
Cornell, C. A.: Engineering seismic risk analysis, Bull. Seismol. Soc. Am., 58, 1583–1606, https://doi.org/10.1016/01676105(83)901435, 1968.
De Biasio, M., Grange, S., Dufour, F., Allain, F., and PetreLazar, I.: A simple and efficient intensity measure to account for nonlinear structural behavior, Earthq. Spectra, 30, 1403–1426, https://doi.org/10.1193/010614EQS006M, 2014.
Di Sarno, L. and Amiri, S.: Period elongation of deteriorating structures under mainshockaftershock sequences, Eng. Struct., 196, 109341, https://doi.org/10.1016/j.engstruct.2019.109341, 2019.
Esteva, L.: Bases para la formulación de decisiones de diseño sísmico, PhD Thesis, Universidad Nacional Autonóma de México, Ciudad de Mexico, Mexico, 1968.
Esteva, L.: Regionalización sísmica de México para fines de ingeniería, series No. 246, Instituto de ingeniería, Ciudad de Mexico, Mexico, 1970.
FEMA273: NEHRP guidelines for the seismic rehabilitation of buildings, Rep. No. FEMA273, Washington, D.C., 1997.
FEMA356: Prestandard and commentary for the seismic rehabilitation of buildings, Rep. No. FEMA356, Washington, D.C., 2000.
FEMA440: Improvement of nonlinear static seismic analysis procedures, Rep. No. FEMA440, Washington, D.C., 2005.
FEMA P58: Seismic performance assessment of buildings volume 1 – methodology, Rep. No. FEMA P581, Washington, D.C., 2012.
Jaimes, M. A., RamirezGaytan, A., and Reinoso, E.: GroundMotion Prediction Model from IntermediateDepth Intraslab Earthquakes at the Hill and LakeBed Zones of Mexico City, J. Earthq. Eng., 19, 1260–1278, https://doi.org/10.1080/13632469.2015.1025926, 2015.
Jamshidiha, H. R., Yakhchalian, M., and Mohebi, B.: Advanced scalar intensity measures for collapse capacity prediction of steel moment resisting frames with fluid viscous dampers, Soil Dynam. Earthq. Eng., 109, 102–118, https://doi.org/10.1016/j.soildyn.2018.01.009, 2018.
Katsanos, E. I. and Sextos, A. G.: Inelastic spectra to predict period elongation of structures under earthquake loading, Earthq. Eng. Struct. Dynam., 44, 1765–1782, https://doi.org/10.1002/eqe.2554, 2015.
Kostinakis, K. and Athanatopoulou, A.: Incremental dynamic analysis applied to assessment of structurespecific earthquake IMs in 3D R/C buildings, Eng. Struct., 125, 300–312, https://doi.org/10.1016/j.engstruct.2016.07.007, 2016.
Kostinakis, K., Fontara, I. K., and Athanatopoulou, A. M.: Scalar structurespecific ground Motion intensity measures for assessing the seismic performance of structures: A review, J. Earthq. Eng., 22, 630–665, https://doi.org/10.1080/13632469.2016.1264323, 2018.
Kramer, S. L.: Geotechnical Earthquake Engineering, PrenticeHall, Upper Saddle River, New Jersey, 1996.
MCBC: Normas técnicas complementarias para diseño por sismo de la xiudad de México, CDMX, México, 2017.
McGuire, R. K.: Probabilistic seismic hazard analysis and design earthquakes: closing the loop, Bull. Seismol. Soc. Am., 85, 1275–1284, https://doi.org/10.1016/01489062(96)833559, 1995.
MCSD – Manual of Civil Structures Design: Manual de Diseño de Obras Civiles Diseño por Sismo, Instituto de Investigaciones Eléctricas, Comisión Federal de Electricidad, Ciudad de Mexico, Mexico, 2015.
Minas, S. and Galasso, C.: Accounting for spectral shape in simplified fragility analysis of casestudy reinforced concrete frames, Soil Dynam. Earthq. Eng., 119, 91–103, https://doi.org/10.1016/j.soildyn.2018.12.025, 2019.
Miranda, E. and RuizGarcia, J.: Influence of stiffness degradation on strength demands of structures built on soft soil sites, Eng. Struct., 24, 1271–1281, https://doi.org/10.1016/S01410296(02)000524, 2002.
Modica, A. and Stafford, P. J.: Vector fragility surfaces for reinforced concrete frames in Europe, Bull. Earthq. Eng., 12, 1725–1753, https://doi.org/10.1007/s105180139571z, 2014.
Montiel, M. A. and Ruiz, S. E.: Influence of structural capacity uncertainty on seismic reliability of buildings under narrowband motions, Earthq. Eng. Struct. Dynam., 36, 1915–1934, https://doi.org/10.1002/eqe.711, 2007.
Ordaz, M. and Reyes, C.: Earthquake hazard in Mexico City: observations versus computations, Bull. Seismol. Soc. Am., 89, 1379–1383, 1999.
Ordaz, M. and Singh, S. K.: Source spectra and spectral attenuation of seismic waves from Mexican earthquakes, and evidence of amplification in the hill zone of Mexico City, Bull. Seismol. Soc. Am., 82, 24–43, 1992.
Ordaz, M., Singh, S. K., Reinoso, E., Lermo, J., Espinosa, J. M., and Dominguez, T.: Mexico earthquake of September 19, 1985 – estimation of response spectra in the lake bed zone of the Valley of Mexico, Earthq. Spectra, 4, 815–834, https://doi.org/10.1193/1.1585504, 1988.
Reinoso, E. and Ordaz, M.: Spectral ratios for Mexico City from freefield recordings, Earthq. Spectra, 15, 273–295, https://doi.org/10.1193/1.1586041, 1999.
Reyes, C., Miranda, E., Ordaz, M., and Meli, R.: Estimación de espectros de aceleraciones correspondientes a diferentes periodos de retorno para las distintas zonas sísmicas de la ciudad de méxico, Rev. Ing. Sísmica, 66, 95–121, https://doi.org/10.18867/ris.66.198, 2002.
RodríguezCastellanos, A., Bojórquez, E., and Ruiz, S. E.: Probabilistic seismic hazard analysis using an advanced intensity measure accounting for structural degradation, in: 12th Canadian conference on earthquake engineering, Quebec, Canada, 1–7, 2019.
RodríguezCastellanos, A., Ruiz, S. E., Bojórquez, E., and ReyesSalazar, A.: Influence of spectral acceleration correlation models on conditional mean spectra and probabilistic seismic hazard analysis, Earthq. Eng. Struct. Dynam., https://doi.org/10.1002/eqe.3331, in press, 2020.
RuizGarcía, J. and Miranda, E.: Inelastic displacement ratios for evaluation of existing structures, Earthq. Eng. Struct. Dynam., 32, 1237–1258, https://doi.org/10.1002/eqe.271, 2003.
Shome, N., Cornell, C. A., Bazzurro, P., and Carballo, J. E.: Earthquakes, records, and nonlinear responses, Earthq. Spectra, 14, 469–500, https://doi.org/10.1193/1.1586011, 1998.
Singh, S. K., Lermo, J., Dominguez, T., Ordaz, M., Espinosa, J. M., Mena, E., and Quaas, R.: Mexico earthquake of September 19, 1985 – a study of amplification of seismic waves in the valley of Mexico with respect to a hill zone site, Earthq. Spectra, 4, 653–673, https://doi.org/10.1193/1.1585496, 1988.
Singh, S. K., CruzAtienza, V., PérezCampos, X., Iglesias, A., Hjörleifsdóttir, V., Reinoso, E., Ordaz, M., and Arroyo, D.: Deadly intraslab Mexico earthquake of 19 September 2017 (Mw 7.1): Ground motion and damage pattern in Mexico City, Seismol. Res. Lett., 89, 2193–2203, https://doi.org/10.1785/0220180159, 2018.
TeránGilmore, A. and Espinosa Johnson, M.: Diseño por desempeño de estructuras dúctiles de concreto reforzado ubicadas en la zona del lago del Distrito Federal: La resistencia lateral de diseño, Rev. Ing. Sísmica, 78, 23–46, https://doi.org/10.18867/ris.78.20, 2008.
Tothong, P. and Luco, N.: Probabilistic seismic demand analysis using advanced ground motion intensity measures, Earthq. Eng. Struct. Dynam., 36, 1837–1860, https://doi.org/10.1002/eqe.696, 2007.
Tsantaki, S., Jäger, C., and Adam, C.: Improved seismic collapse prediction of inelastic simple systems vulnerable to the Pdelta effect based on average spectral acceleration, in: Proc. 15th World Conference on Earthquake Engineering, 24–28 September 2012, Lisbon, Portugal, 2012.
Tsantaki, S., Adam, C., and Ibarra, L. F.: Intensity measures that reduce collapse capacity dispersion of Pdelta vulnerable simple systems, Bull. Earthq. Eng., 15, 1085–1109, https://doi.org/10.1007/s1051801699944, 2017.
Yakhchalian, M., Nicknam, A., and Amiri, G. G.: Optimal vectorvalued intensity measure for seismic collapse assessment of structures, Earthq. Eng. Eng. Vib., 14, 37–54, https://doi.org/10.1007/s1180301500056, 2015.
Zúñiga, F. R., Suárez, G., FigueroaSoto, Á., and Mendoza, A.: A firstorder seismotectonic regionalization of Mexico for seismic hazard and risk estimation, J. Seismol., 21, 1295–1322, https://doi.org/10.1007/s1095001796660, 2017.
 Abstract
 Introduction
 Methodology proposed
 Probabilistic seismic hazard analysis (PSHA)
 ${I}_{{N}_{\mathrm{p}}}$ intensity measure
 Probabilistic seismic hazard analysis using ${I}_{{N}_{\mathrm{p}}}$
 Degrading structural behavior effect
 Conclusions
 Data availability
 Author contributions
 Competing interests
 Acknowledgements
 Financial support
 Review statement
 References
 Abstract
 Introduction
 Methodology proposed
 Probabilistic seismic hazard analysis (PSHA)
 ${I}_{{N}_{\mathrm{p}}}$ intensity measure
 Probabilistic seismic hazard analysis using ${I}_{{N}_{\mathrm{p}}}$
 Degrading structural behavior effect
 Conclusions
 Data availability
 Author contributions
 Competing interests
 Acknowledgements
 Financial support
 Review statement
 References