# Urban Seismic Risk Assessment and Damage Estimation: Case of Rif Buildings (North of Morocco)

Seif-Eddine Cherif, Aboubakr Chaaraoui, Mimoun Chourak, Mohamed Oualid Mghazli, Abdelhay El Omari, Tiago Miguel Ferreira

► To cite this version:

Seif-Eddine Cherif, Aboubakr Chaaraoui, Mimoun Chourak, Mohamed Oualid Mghazli, Abdelhay El Omari, et al.. Urban Seismic Risk Assessment and Damage Estimation: Case of Rif Buildings (North of Morocco). Buildings, 2022, 12 (6), pp.742. 10.3390/buildings12060742. hal-04082657

HAL Id: hal-04082657

https://hal.science/hal-04082657v1

Submitted on 11 Mar 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Urban Seismic Risk Assessment and Damage Estimation: Case of Rif Buildings (North of Morocco)

Seif-Eddine Cherif <fnref n="1" /> ID, Aboubakr Chaaraoui <fnref n="2" />, Mimoun Chourak <fnref n="3" /> ID, Mohamed Oualid Mghazli <fnref n="4" />, Abdelhay EL Omari <fnref n="3" /> ID and Tiago Miguel Ferreira <fnref n="5" />,* ID

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Cherif, S.-E.; Chaaraoui, A.; Chourak, M.; Oualid Mghazli, M.; EL Omari, A.; Ferreira, T.M. Urban Seismic Risk Assessment and Damage Estimation: Case of Rif Buildings (North of Morocco).

Buildings **2022**, *12*, 742. https://doi.org/10.3390/buildings12060742

Academic Editor: Lucio Soibelman

Received: 9 May 2022

Accepted: 28 May 2022

Published: 30 May 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** The interest in assessing seismic risk in earthquake-prone regions in Morocco has been increasing over recent decades, to a large extent due to the substantial amount of damage imposed by recent events and because the population in these regions has grown significantly. In this context, the present study is aimed to contribute to the understanding of seismic risk in the urban areas of the Rif region, one of the most seismically active zones of Morocco, through the development and analysis of 36 building models representative of this area. Two earthquake scenarios were considered for the assessment of the seismic hazard, based on the national seismic code and the European code adapted to local seismic parameters. The performance points, determined following generated response and capacity spectra, made it possible to establish damage probability matrices. Obtained results corroborate those of previous reports, confirming that the damage is more significant in Imzouren due to the nature of the soil. It has also been shown that the credibility of the response spectra drawn from the national code is questioned, given the extreme damage estimated. The adapted European spectrum proved to be a more reliable probabilistic earthquake scenario for damage estimation.

**Keywords:** capacity spectrum method; capacity curves; fragility curves; response spectra; performance point; damage states

## 1. Introduction

In Morocco, the rural world is lagging in terms of development indicators such as weak infrastructure, isolation, and insecurity. These structural weaknesses, combined with the droughts that occasionally affect the country, have strongly contributed to the rural exodus. In 1960, 71% of Moroccans lived in rural areas [1]. In recent years, and mainly as a result of the rural exodus, less than 40% of the country's population still live in rural areas [2].

Across the country, urban planning has not been able to keep up with the significant increase in population following rural migration. The inadequate occupancy of the land contributes to increasing the damage due to seismic catastrophes. The high concentration of people, buildings, infrastructures, and exposed values turn these zones into high-risk areas. Districts that do not meet building standards are found in every city in Morocco. However, by adding moderate seismic activity to the equation, it becomes a recipe for disaster. This is the case for the urban and rural areas in the central Rif region.

<footnote n="1">¹ Georesources, Geoenvironment and Civil Engineering Laboratory (L3G), Faculty of Science and Technology, Cadi Ayyad University, Marrakech 40000, Morocco; s.cherif@uca.ma</footnote>

<footnote n="2">² Laboratory of Applied Geosciences, Department of Geology, Faculty of Sciences, Mohamed First University, Oujda 60000, Morocco; aboubakr.chaaraoui@gmail.com</footnote>

<footnote n="3">³ Mechanics and Applied Mathematics Department, Industrial and Seismic Engineering Research Team, National School of Applied Sciences of Oujda, Mohammed First University, Oujda 60000, Morocco; m.chourak@ump.ac.ma (M.C.); elomari_abdelhay1718@ump.ac.ma (A.E.O.)</footnote>

<footnote n="4">⁴ Tribology and System Dynamics Laboratory (LTDS UMR CNRS 5513), ENTPE, University of Lyon, 69120 Lyon, France; mohamed.oualid.mghazli@gmail.com</footnote>

<footnote n="5">5 Department of Geography and Environmental Management, University of the West of England, Bristol BS16 1QY, UK</footnote>

<footnote>* Correspondence: tiago.ferreira@uwe.ac.uk</footnote>
The recent past has demonstrated that seismic risk is significant in this area [3,4]. Successive earthquake sequences between 1992 and 1994 [5,6], the 2004 earthquake [7–9] and, more recently, the 2016 earthquake [10–12], have caused fatalities and significant damage. Moderate seismicity combined with a high vulnerability of existing buildings and local site effects were mainly believed to be the causes of the extensive damage. On a larger scale, the Mediterranean region is no stranger to earthquakes of moderate intensity, causing significant damage. On the 12 June 2017, an earthquake of Mw = 6.3 struck Lesvos Island, causing one human fatality and severe damage to the built environment [13]. An Mw = 6.4 earthquake hit the NW region of Albania on 26 November 2019, resulting in extensive damage to the civil structures in the broader area of Durrës city and its surroundings [14]. More recently, a severe earthquake hit Zagreb on 22 March 2020 (magnitude ML = 5.5). The event occurred during the COVID-19 lockdown and caused significant damage to the built environment and enormous disruption in everyday life [15].

Pre-emptive measures were eventually taken to reduce the damage. The introduction of a national seismic standard R.P.S. in 2002 [16] and its revision in 2011 [17] were one of the most important steps, along with multiple studies to estimate risk aiming to guide the decision making at the urban level. In the literature, several methods are suitable to solve this kind of problem on a large scale and can be categorized into three classes [18,19]: (i) Empirical methods, based on damage reports following earthquakes and building characteristics to estimate a seismic vulnerability index. These empirical methods were the only reasonable approaches that could be used in a large-scale seismic hazard analysis [20–26]; (ii) mechanical and analytical methods, which are more detailed but laborious for a large-scale study. They use slightly more detailed algorithms for vulnerability assessments, which allow more delicate and detailed studies [27–30]; (iii) hybrid methods, which can be particularly advantageous in the absence of damage data at certain levels of intensity for the geographical area considered and they also make it possible to adjust the analytical model. Furthermore, the use of observational data reduces the computational effort that would be required to produce a full set of analytical vulnerability curves [31,32].

In recent years, the use of the Georeferenced Information System (GIS) has gained momentum in civil engineering applications, especially in the field of risk assessment on a large scale. One cannot deny the increase in popularity of data collection for seismic vulnerability and risk assessment and some proposals have been developed to rapidly obtain a large amount of information processing images of wide sets of existing buildings. For instance, the Global Earthquake Model (GEM) has been developing an open global earthquake risk model [33]. These efforts have led to the development of a repository of probabilistic seismic hazard models, a global exposure dataset comprising structural and occupancy information regarding the residential, commercial, and industrial buildings, and a comprehensive set of fragility and vulnerability functions for the most common building classes. These components were used to estimate probabilistic earthquake risk globally using the OpenQuake Engine, an open-source software for seismic hazard and risk analysis [34,35].

In Morocco, the first adapted methods are the European macroseismic scale and vulnerability index methods [36–38]. These approaches are time consuming in gathering data, but they are easy to apply and show conclusive results. Other studies like the seismic index [39] tried to tackle the problem differently by adapting the Japanese Seismic Index Method, making it more time efficient and tuned to the context of Moroccan construction.

The present study follows the same strategy of estimating damage in the region by considering an advanced model to evaluate seismic risk in the cities of Al Hoceima and Imzouren. It is a predictive model obtained by applying the Capacity Spectrum Based Method (CSBM). The method considers four damage states in addition to a no-damage state [40]. Seismic action is defined in terms of response spectra and building vulnerability is evaluated through its capacity spectrum. A probabilistic seismic hazard scenario with a probability of exceedance of 10% in 50 years, in which soil effects are considered, is used herein.
Capacity curves are developed from building models representative of the building typologies present in the study area. These models have been designed in previous studies based on a database gathered over the years [37–39]. Since 2014, information on existing buildings in the region has been collected, catalogued, and analyzed. The gathered information is based on field investigations, design and reinforced plans of some buildings, and information from consulting civil engineering offices working on residential buildings in the region.

The obtained results are stored and displayed in a GIS model where the predictive damage of buildings in both cities is spatially represented. The detailed risk maps obtained offer an excellent opportunity to guide the decision making in the field of seismic risk prevention, preparedness, and mitigation in the central Rif region. The study would also constitute the first step towards a newer and more reliable version of the national seismic standard R.P.S. [17].

## 2. Capacity Spectrum Based Method

The Capacity Spectrum Based Method (CSBM) [41] is a performance-based seismic analysis technique and may be used for a variety of purposes and on different scales. The approach can potentially be a time-efficient evaluation of a large inventory of buildings on a large scale, or on a smaller scale, an assessment of the seismic performance of a new or an existing structure to estimate its damage state. The CSBM was developed by Freeman to be a rapid evaluation procedure for assessing the seismic vulnerability of buildings at the Puget Sound Naval Shipyard [42]. The procedure compares the capacity of a structure with the demand of earthquake motion on the structure.

The capacity of the structure is represented in terms of a force-displacement curve obtained from a non-linear static (pushover) analysis. The shear force and roof displacement are converted to a spectral acceleration and spectral displacement of an equivalent Single Degree of Freedom (SDOF) system. The earthquake demand curve is represented by response spectra. An Acceleration/Displacement Response Spectrum (ADRS) format [43], where spectral accelerations are plotted against spectral displacements, is used. The graphical superposition of the two curves makes possible a visual evaluation of how the structure will perform under earthquake ground motion.

The CSBM considers five damage states: slight, moderate, severe, and complete, including the none damage state. For reinforced concrete buildings, these damage states have the same meaning as in HAZUS [44]. The only exception is the complete damage state, which is subdivided into two states: very heavy and destruction, as per the 1998 European Macroseismic Scale [45]. A weighted average damage index, $DS_m$ can be calculated as:

$$DS_m = \sum_{k=0}^{4} k \cdot P[DS_k] \qquad (1)$$

where k takes the values 0, 1, 2, 3, and 4 for the damage states k considered in the analysis and $P[DS_k]$ are the corresponding probabilities. It can be considered that $DS_m$ is close to the most likely damage state of the structure. The damage index $DS_m$ will be useful for representing and analyzing damage distributions by using a single parameter.

Fragility curves define the probability that the expected damage d of a structure exceeds a given damage state $ds_i$, as a function of a parameter quantifying the severity of the seismic action. In this study, the quantifying parameter is the spectral displacement $Sd$. It is usually assumed that fragility curves are well described by the following lognormal probability density function:

$$P\left[\frac{ds_i}{Sd}\right] = \phi\left[\frac{1}{\beta_{ds_i}} \ln\left(\frac{Sd}{\bar{Sd}_{ds_i}}\right)\right] \qquad (2)$$
where $\overline{Sd}_{ds_i}$ is the threshold spectral displacement at which the probability of the damage state $ds_i$ is 50%, $\beta_{ds_i}$ is the standard deviation of the natural logarithm of this spectral displacement, $\phi$ is the standard cumulative distribution function, and $Sd$ is the spectral displacement.

Capacity spectra are approximated with bilinear curves using FEMA 273 guidelines [46]. Table 1 shows how the thresholds $\overline{Sd}_{ds_i}$ are obtained as a function of the yield displacement $D_y$ and the ultimate displacement $D_u$ of the structure. Concerning $\beta_{ds_i}$, it is well known that the expected seismic damage in buildings follows a binomial probability distribution [45]. $\beta_{ds_i}$ is calculated directly as a function of ultimate ductility $\mu_u$ [40] and is displayed in Table 1.

<table><caption>**Table 1.** Damage state thresholds.</caption><thead><tr><td>Damage State Thresholds</td><td>Definition</td><td>β<sub>ds</sub></td><td>Damage State</td></tr></thead><tbody><tr><td>Sd<sub>1</sub></td><td>0.7 D<sub>y</sub></td><td>β<sub>Sd1</sub> = 0.25 + 0.07 ln(μ<sub>u</sub>)</td><td>Slight</td></tr><tr><td>Sd<sub>2</sub></td><td>D<sub>y</sub></td><td>β<sub>Sd2</sub> = 0.2 + 0.18 ln(μ<sub>u</sub>)</td><td>Moderate</td></tr><tr><td>Sd<sub>3</sub></td><td>D<sub>y</sub> + 0.25(D<sub>u</sub> - D<sub>y</sub>)</td><td>β<sub>Sd3</sub> = 0.1 + 0.4 ln(μ<sub>u</sub>)</td><td>Severe</td></tr><tr><td>Sd<sub>4</sub></td><td>D<sub>u</sub></td><td>β<sub>Sd4</sub> = 0.15 + 0.5 ln(μ<sub>u</sub>)</td><td>Complete</td></tr></tbody></table>

### 3. Description of the Study Area

#### 3.1. Seismicity in the Cities of Al Hoceima and Imzouren

The north region of Morocco is the most seismic area of the country. In recent decades, the region has experienced several earthquakes, some very violent, as was the case of the 1994 (MW = 6.0) and 2004 (MW = 6.4) earthquakes [47]. Recently, a 6.3 magnitude earthquake struck 50 km off the coast of Al Hoceima on Monday, January 26th 2016 [10–12] causing property damage. The central Rif area, in particular, has been affected by several recent seismic episodes, which proves the high level of seismicity of this zone [48–50]. This issue has been the subject of multiple studies related to seismology, seismic hazard, and seismic vulnerability [37–39,51,52].

Al Hoceima (35°25′ N, 3°93′ W) and Imzouren (35°09′ N, 3°52′ W) are cities located on the northern coast of Morocco, both in the same province (Figure 1). Being at the heart of the Rif chain, they have been subject to many seismic episodes. As a reminder, the earthquake of 24 February 2004 caused 629 deaths, 966 injuries, 2539 households damaged or destroyed, and left 15,600 people homeless. The difference in observed damage following the same earthquake was significant. Even though Imzouren and Al Hoceima are equidistant from the epicenter, Imzouren was damaged the most. It was reported that the damage differences were linked to the nature of the soil and construction quality [7]. While it is almost impossible to check the quality of construction for existing buildings, according to Cherif et al. [39], the nature of soil definitely played a role in the damage differences observed.

As previously stated, several studies have targeted the two cities given the consecutive seismic episodes, population growth, and the significant vulnerability of existing structures, especially those recently built due to rural exodus. The probability of significant damage due to earthquakes will only grow because of these factors.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** (a) Landsat satellite image of western Mediterranean (Google Earth) and location of the study area. (b) Satellite image (Google Earth) of the cities of Al Hoceima and Imzouren.</figcaption></figure>

## 3.2. Seismic Hazard

The seismic hazard is defined in terms of 5% damped elastic response spectra from a probabilistic point of view corresponding to a 475-year return period. Two response spectra are proposed: a response spectrum taken from the national seismic code [17] and another based on Eurocode 8 [53], considering local parameters. The second scenario is considered, because building conception and construction in Morocco are still profoundly influenced by European construction standards. To this day, the French construction code, BAEL91 [54], and the Eurocode are still considered in our construction practices.

RPS2000 does not contain specific site response spectra and shows many limitations in terms of the typology of the structures and their regularity. Any building project outside of those limitations would require a complementary study following the Eurocode 8. Thus, an additional probabilistic scenario has been considered: Moroccan seismic parameters used in the basic representation of the seismic action in Eurocode 8.

The national seismic standard [17] introduces two seismic zoning maps, one in terms of acceleration and the other in terms of velocity. Both maps contain five seismic zones (Figure 2), in which acceleration values range between 0.04 g and 0.18 g and seismic velocity from 5 cm/s to 17 cm/s. The expected elastic response spectrum $S_v(T)$ is defined in Table 2. Five soil types are considered: S1 rocky ground, S2 firm ground, S3 loose soil, S4 soft soil, and S5 special-case soil. The corresponding shear-wave velocity ranges are displayed in Table 3. According to Eurocode 8, two types of spectra are proposed: Type 1 and Type 2. Since the Rif region, where Imzouren and Al Hoceima are located, is subjected to moderate seismic activity, Type 2 response spectrum seems to be the most suitable option.

<table><caption>**Table 2.** Seismic acceleration values based on the national seismic standard [17].</caption><thead><tr><th>Z<sub>a</sub>/Z<sub>v</sub> Ratio *</th><th>0 ≤ T ≤ 0.25</th><th>0.25 ≤ T ≤ 0.5</th><th>0.5 ≤ T</th></tr></thead><tbody><tr><td>Z<sub>a</sub>/Z<sub>v</sub> > 1</td><td>1.9 S·a<sub>g</sub> **</td><td>1.9 S·a<sub>g</sub></td><td></td></tr><tr><td>Z<sub>a</sub>/Z<sub>v</sub> = 1</td><td>2.5 S·a<sub>g</sub></td><td>(-2.4 T + 3.1) S·a<sub>g</sub></td><td>1.20 a<sub>g</sub>·S/(T)<sup>2/3</sup></td></tr><tr><td>Z<sub>a</sub>/Z<sub>v</sub> < 1</td><td>3.5 S·a<sub>g</sub></td><td>(-6.4T + 5.1) S·a<sub>g</sub></td><td></td></tr></tbody></table>

<footnote>* Z_a is the acceleration value and Z_v is the velocity value corresponding to the seismic zone. ** S is the soil factor, having values displayed in Table 3.</footnote>
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Seismic zoning in terms of (a) acceleration and velocity and (b) for probabilities of 10% in 50 years according to the national seismic code (RPS2000) [17].</figcaption></figure>

<table><caption>**Table 3.** Ground types, according to Eurocode 8 [53] and the national seismic standard (RPS2000) [17].</caption><thead><tr><th rowspan="2">Description of Stratigraphic Profile</th><th colspan="2">EUROCODE 8</th><th colspan="3">National Seismic Code (RPS2000)</th></tr><tr><th>Ground Type</th><th>V<sub>s</sub> (m/s)</th><th>Soil Class</th><th>V<sub>s</sub> (m/s)</th><th>Soil Factor</th></tr></thead><tbody><tr><td>Rock or other rock-like geological formation</td><td>A</td><td>&gt;800</td><td>S1</td><td>V<sub>s</sub> ≥ 760</td><td>1</td></tr><tr><td>Deposits of very dense sand, gravel, or very stiff clay</td><td>B</td><td>360–800</td><td>S2</td><td>360–760</td><td>1.2</td></tr><tr><td>Deep deposits of dense or medium dense sand, gravel, or stiff clay</td><td>C</td><td>180–360</td><td>S3</td><td>180–360</td><td>1.4</td></tr><tr><td>Deposits of a loose-to-medium cohesionless soil or predominantly soft-to-firm cohesive soil</td><td>D</td><td>&lt;180</td><td>S4</td><td>&lt;180</td><td>1.8</td></tr><tr><td>A soil profile consisting of a surface alluvium layer</td><td>E</td><td>-</td><td>S5</td><td>-</td><td>-</td></tr></tbody></table>

For each town, two response spectra are proposed (one for each soil type). The first is based on the RPS2000 and the second on Eurocode 8. For rocky soil, the RPS assigns the same acceleration value of 0.18 g to the towns of Imzouren and Al Hoceima. For the second spectrum, and in keeping up with the most recent works, acceleration values are based on a conducted study on risk analysis by the RMSI (Risk Management Solutions Inc.) in 2012 [55], which provides the following $a_g$ results:

* Al Hoceima: a peak ground acceleration of 0.253 g for a return period of 475 years;

* Imzouren: a peak ground acceleration of 0.303 g for a return period of 475 years.

Figure 3 shows the obtained 5% damped elastic response spectra on rocky soil. As we can see, the RPS response spectrum is relatively flat compared with the two other spectra, which decreases the amplification values for the short return periods but dampens the decline effect of the accelerations values the more the return period increases.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** The 5% elastic response spectra for the probabilistic scenario expected on outcrop rock.</figcaption></figure>

As for the other types of soil, Table 3 shows that both seismic standards introduce the same soil classification. Figures 4 and 5 show soil distribution in the towns of Al Hoceima and Imzouren, respectively. The soil type will undoubtedly come into play in the evaluation of seismic risk between the two cities. The structures of Al Hoceima are mostly built on firm or rocky ground, while the structures of Imzouren are mostly built on relatively soft soils. All of the response spectra for each soil type are shown in Figure 6.

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Distribution of soil types in the city of Al Hoceima [37].</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Distribution of soil types in the city of Imzouren [56].</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>(b)</figcaption></figure>

**Figure 6.** *Cont.*
<figure><img src="image_10.png" /><figcaption>**Figure 6.** Response spectra of the cities of Al Hoceima and Imzouren for each soil type: (a) Eurocode 8 response spectra Al Hoceima; (b) Eurocode 8 response spectra Imzouren; (c) RPS response spectra Al Hoceima and Imzouren.</figcaption></figure>

### 3.3. *Building Typologies*

Urban communes within the central Rif area are modern twentieth-century establishments. According to the most recent census [2], the modern Moroccan house is the dominant building typology, see Figure 7. It is usually a low- to mid-rise structure with a reinforced concrete moment frame. The existing constructions in Al Hoceima and Imzouren are no exception and have similar features. Most of the buildings are regular in geometry and do not exceed six stories (Figure 8). They also have relatively small construction surfaces, ranging between 100 m² and 150 m² [39].

<figure><img src="image_11.png" /><figcaption>**Figure 7.** Distribution of building typologies in Al Hoceima and Imzouren of Central Rif [2].</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>**Figure 8.** Examples of residential buildings in the study region.</figcaption></figure>

The used structural models that allowed us to generate capacity and fragility curves have been obtained by considering the constructive peculiarities of the existing reinforced concrete buildings. Detailed information on their design has been obtained through years of collecting, georeferencing, arranging, classifying, and improving the database of existing buildings in both cities [37,38]. The database collection initially started with a preliminary survey in the cities to determine the best method for assessing the seismic performance of existing buildings. In addition to that, building design and conception, obtained from design offices working in the study area, helped us generate 3D building models designed according to the existent structural typologies [39].

The inspection focused mainly on residential buildings spread across the two cities representing the different districts and typologies. A total number of 2746 existing buildings between the two cities was considered (Figures 9 and 10). Information like the structural system, position, number of floors, observed irregularities, etc., were stored in a GIS database. The following observations highlighted the existence of several similarities, which made it possible to classify the buildings according to a reduced number of typologies:

* Most buildings are reinforced concrete moment frame structures;

* Low-rise structures built on relatively small surfaces (100–150 m²);

* Most buildings are regular structures having simple geometrical shapes.

<figure><img src="image_13.png" /><figcaption>**Figure 9.** *Cont.*</figcaption></figure>
<figure><img src="image_14.png" /><figcaption>**Figure 9.** Inspected buildings in the city of Al Hoceima based on (a) seismic code and (b) construction period.</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>**Figure 10.** *Cont.*</figcaption></figure>
<figure><img src="image_16.png" /><figcaption>**Figure 10.** Inspected buildings in the city of Imzouren based on (a) seismic code and (b) construction period.</figcaption></figure>

A total of 2746 buildings were inspected and classified into 36 typologies according to three factors: seismic code, construction period, and the number of floors. Each typology is represented by a 3D building model that can be used for a suitable and convenient illustration of these structures.

# 4. Results and Discussion

## 4.1. *Structural Capacity of the Buildings in Al Hoceima and Imzouren*

The 36 building models are designed to represent a large part of the existing residential buildings in the urban municipalities of the central Rif. The observations collected as a result of the field investigations, the analysis of the reinforced concrete plans, and the consultations carried out with the design offices are all aspects considered in the development of the models. Even construction practices, often prioritizing the cost of building materials over build quality, were included in the analysis. The seismic performance of a building, represented in this study by a building model, can be characterized by its capacity spectrum obtained by means of a pushover analysis [57]. This capacity spectrum is usually modeled in its simplified bilinear form defined by the yielding (Dy, Ay) and ultimate capacity points (Du, Au).

As stated before, seismic code, construction period, and the number of stories are the main factors contemplated in the 36 building models. All existent buildings in the target area are low- to mid-rise constructions, so the considered number of floors ranges between one and six stories. Constructions were built with different types of materials (steel and concrete) depending on the construction period. Thus, three different construction periods are considered: before 1991, between 1991 and 2002, and after 2002 [39]. The third and final factor is the seismic code and there is only one national seismic code (RPS 2000), introduced in 2002 and revised in 2011. Buildings respecting the seismic code are designed to support
seismic loads and they are projected essentially in the cross-sectional area of the columns, giving more strength to vertical structural elements.

The floor area of the built models has a rectangular shape and is equal to 12 m by 10 m. The structural nature of the models is an RC frame since this is the structural system that is most frequently found in residential dwellings. Particular emphasis was given to the vertical structural elements because of their importance in resisting horizontal loads. The cross-sectional area of the vertical structural element and the disposition of steel rebar for the adopted 3D models were calculated similarly to the ones accomplished by design offices. The following mechanical properties have been assumed: concrete compression strength is equal to 16 MPa or 22 MPa and steel yield stress is equal to 235 MPa, 400 MPa, or 500 MPa, depending on the construction period (Table 4). Constructions were built with different types of materials (steel and concrete) depending on the period of construction. For instance, Fe E 235 steel rebar was commonly used before 1991, while Fe E 500 was more frequently used after the 2000s. In addition, concrete compression strength is taken equal to 22 MPa for buildings with seismic code (after RPS2000) and 16 MPa (before RPS2000) for those without seismic code. This assumption, of course, is only a proxy for real strength values. However, these two values (16 MPa and 22 MPa) are target values for both the construction company and the client. Any value below would not respect regulations and the values above would lead to unnecessary additional costs [39].

<table><caption>**Table 4.** Construction material characteristics for defined periods [39].</caption><thead><tr><td>Construction Period</td><td>Concrete Strength</td><td>Steel Rebar</td></tr></thead><tbody><tr><td>Period I (Pre 1991)</td><td>16 MPa</td><td>Fe E 235</td></tr><tr><td>Period II (1991–2002)</td><td>16 MPa</td><td>Fe E 400</td></tr><tr><td>Period III (Post 2002)</td><td>22 MPa</td><td>Fe E 500</td></tr></tbody></table>

Based on the results of this study, the longest fundamental period is equal to 0.74 s, found at the level of the highest buildings (six floors for our case). The existing buildings, like their corresponding models, are regular. We can therefore assume that the response of structures is represented by their fundamental mode of vibration. This makes pushover analysis, despite its limitations, suitable for assessing seismic risk scenarios in the towns of Al Hoceima and Imzouren.

For a developing country like Morocco, savings on construction materials is important and is of great concern to construction companies, to the point of altering and penalizing the behavior of structures. One example would be the cross-sectional area of columns, which can change from one story to another if the load is theoretically supported, which can potentially compromise the performance of the building due to weak linking. Even more, it could indirectly increase the price since the connection of the nodes is more complex, requiring skilled labor and taking more time to complete.

Considering these construction practices, the following models have been proposed [39]. The designed structures have a rectangular shape and a 12 m × 10 m surface (Figure 11). The construction area of most residential buildings in the study area varies between 100 m² and 150 m². Taking a surface of 120 m² is suitable and convenient for the representation of the inspected buildings. The structural nature of the models is an RC frame since it is the main type found in residential dwellings. Tables 5 and 6 show the cross-sectional area and disposition of steel rebar for the adopted 3D models. The loading and sizing calculations are carried out in a similar manner to that accomplished by design offices. Each building model was subjected to a static non-linear pushover analysis in the weakest direction in order to obtain the base shear versus roof displacement curves [57].
<figure><img src="image_17.png" /><figcaption>**Figure 11.** Distribution of column types for constructed models [39].</figcaption></figure>

<table><caption>**Table 5.** Cross-sectional area of columns and adopted steel rebar for buildings with seismic code [39].</caption><thead><tr><th>Cross-Sectional Area (cm²)</th><th>Floor 1</th><th>Floor 2</th><th>Floor 3</th><th>Floor 4</th><th>Floor 5</th><th>Floor 6</th></tr></thead><tbody><tr><td>P1</td><td>30 × 30</td><td>30 × 30</td><td>25 × 30</td><td>25 × 25</td><td>25 × 25</td><td>25 × 25</td></tr><tr><td>P2</td><td>30 × 40</td><td>25 × 40</td><td>25 × 40</td><td>25 × 30</td><td>25 × 30</td><td>25 × 25</td></tr><tr><td>P3</td><td>30 × 40</td><td>25 × 40</td><td>25 × 40</td><td>25 × 35</td><td>25 × 30</td><td>25 × 25</td></tr><tr><td>P4</td><td>45 × 45</td><td>40 × 40</td><td>35 × 35</td><td>35 × 35</td><td>30 × 30</td><td>25 × 25</td></tr><tr><td>Steel rebar</td><td>Floor 1</td><td>Floor 2</td><td>Floor 3</td><td>Floor 4</td><td>Floor 5</td><td>Floor 6</td></tr><tr><td>P1</td><td>6T14</td><td>6T14</td><td>4T14 + 2T12</td><td>6T12</td><td>6T12</td><td>6T12</td></tr><tr><td>P2</td><td>8T14</td><td>6T14</td><td>6T14</td><td>4T14 + 2T12</td><td>4T14 + 2T12</td><td>6T12</td></tr><tr><td>P3</td><td>8T14</td><td>6T14</td><td>6T14</td><td>6T14</td><td>4T14 + 2T12</td><td>6T12</td></tr><tr><td>P4</td><td>8T16</td><td>8T16</td><td>8T14</td><td>8T14</td><td>6T14</td><td>6T12</td></tr></tbody></table>

<table><caption>**Table 6.** Cross-sectional area of columns and adopted steel rebar for buildings without seismic code [39].</caption><thead><tr><th>Cross-Sectional Area (cm²)</th><th>Floor 1</th><th>Floor 2</th><th>Floor 3</th><th>Floor 4</th><th>Floor 5</th><th>Floor 6</th></tr></thead><tbody><tr><td>P1</td><td>25 × 25</td><td>25 × 25</td><td>25 × 25</td><td>20 × 20</td><td>20 × 20</td><td>20 × 20</td></tr><tr><td>P2</td><td>25 × 30</td><td>25 × 30</td><td>25 × 25</td><td>25 × 25</td><td>20 × 20</td><td>20 × 20</td></tr><tr><td>P3</td><td>25 × 30</td><td>25 × 30</td><td>25 × 25</td><td>25 × 25</td><td>20 × 20</td><td>20 × 20</td></tr><tr><td>P4</td><td>40 × 40</td><td>35 × 35</td><td>30 × 30</td><td>25 × 30</td><td>25 × 25</td><td>20 × 20</td></tr><tr><td>Steel rebar</td><td>Floor 1</td><td>Floor 2</td><td>Floor 3</td><td>Floor 4</td><td>Floor 5</td><td>Floor 6</td></tr><tr><td>P1</td><td>4T12</td><td>4T12</td><td>4T12</td><td>4T10</td><td>4T10</td><td>4T10</td></tr><tr><td>P2</td><td>6T12</td><td>6T12</td><td>4T12</td><td>4T12</td><td>4T12</td><td>4T12</td></tr><tr><td>P3</td><td>4T14 + 2T12</td><td>6T12</td><td>4T12</td><td>4T12</td><td>4T10</td><td>4T10</td></tr><tr><td>P4</td><td>8T14</td><td>6T14</td><td>6T14</td><td>6T12</td><td>4T12</td><td>4T10</td></tr></tbody></table>

Capacity curves were generated and allowed us to obtain capacity spectra, as can be shown in an example in Figure 12. One can clearly see from the results that the spectral displacement increases while the spectral acceleration decreases as the height of the building increases, as one would expect. The curves only represent six models, but the rest of the models show a similar pattern. Subsequently, the capacity spectra were approximated to bilinear curves using the FEMA 273 guidelines [46]. The set of points representing the yield limit and the ultimate limit of the 36 models is shown in Tables A1 and A2 (Appendix A).
<figure><img src="image_18.png" /><figcaption>**Figure 12.** Capacity spectra of buildings with seismic code built after 2002.</figcaption></figure>

## 4.2. Fragility Curves and Damage Index for the Buildings of Al Hoceima and Imzouren

Fragility curves of the building models were obtained after having calculated the $\overline{Sd}_i$ and $\beta_i$ parameters used in Equations (2) and (3). The complete results of the listed parameters are given in Table A3. As an example of buildings of the same height, we chose to display the fragility curves of four-story RC building models in Figure 13. We notice no significant difference in spectral displacements between buildings that comply with the seismic regulations and those that do not. This difference is mainly reflected at the level of spectral accelerations, as can be seen in Tables A1 and A2. However, there is a noticeable change in terms of the construction period; newly constructed buildings perform slightly better than older buildings in terms of spectral displacement. These results show the high vulnerability of the buildings in the target area. Even in the case of a moderate seismic event, the damage probability is significant.

Fragility curves enable us to estimate the seismic performance of the structure for a given spectral displacement. Therefore, it is important to determine a spectral displacement indicative of the most probable performance of the structure in the event of a seismic event. The spectral displacement resulting from the performance point is exactly the desired parameter, derived from the intersection between the capacity spectrum and the response spectrum as previously indicated. Table A4 represents said displacements for the probabilistic scenarios considered in the study.

The damage probability matrices are then generated from the limit damage states $\overline{Sd}_i$ with $i = 1$ to 4, presented above in Table A3, and the performance points are given in Table A4. For some cases, performance points could not be computed because of the relatively high RPS acceleration values in the higher fundamental period range (Figure 3).
<figure><img src="image_19.png" /></figure>

<figure><img src="image_20.png" /></figure>

<figure><img src="image_21.png" /></figure>

<figure><img src="image_22.png" /></figure>

<figure><img src="image_23.png" /></figure>

<figure><img src="image_24.png" /><figcaption>Figure 13. Fragility curves for 4-floor RC building models.</figcaption></figure>
Then, the weighted average damage index $DS_m$ is calculated, which is used to represent the damage distribution using one parameter only (Table A5). The values of the damage index generate some expected results. In fact, the buildings respecting the seismic code have a better performance and the newly built structures also perform better. The highest damage index values are mostly for structures built on soft soils S4. There are also significant differences when we compare the results of the probabilistic scenarios (RPS and Eurocode 8). According to the first scenario, the higher the rise of the building the more the damage increase. This increase is more significant for loose and soft soils. For the second scenario, the values are more balanced with respect to the elevation of the structure. In addition, we note that on average, the $DS_m$ values of the RPS scenario are a little higher than those from the Eurocode 8 scenario.

### 4.3. Probabilistic Risk Scenarios

Damage distribution according to probabilistic risk scenarios in Al Hoceima and Imzouren is discussed in this section. The average damage indices $DS_m$ calculated in the previous section will be represented for each of the 2746 buildings targeted in this study. The results are shown in Figure 14 for the town of Al Hoceima and in Figure 15 for Imzouren. The average value of the damage index in Al Hoceima is equal to 2.05 and 1.78 according to the RPS and Eurocode 8 scenarios, respectively, which corresponds to a "moderate" damage state. The damage index in the town of Imzouren is equal to 2.77 and 2.75 according to the RPS and Eurocode 8 risk scenarios, which corresponds to a "heavy" damage state.

The distribution of damage in the city of Al Hoceima and Imzouren is quite similar to the results found in previous studies [37–39]. The disparity of results in Al Hoceima comes mainly from the vulnerability of the structures. The seismic hazard uniformly affects the whole city and site effects are minimal since the soil is mostly of S1 or S2 types. Buildings with high vulnerabilities are structures that have been built without considering seismic regulations, often having three or four stories. These structures are being built on the outskirts of Al Hoceima, indicating that it is generally linked to the rural exodus. There are also old vulnerable structures located randomly in the city center, which have not been reinforced since the 2004 earthquake. As for the town of Imzouren, the distribution of damage is more linked to the type of soil. We can clearly see that the distribution of damage follows the pattern of soil typology shown in Figure 5. The buildings with the highest damage indices are structures that do not respect the seismic code and are built on loose soils.

In terms of risk scenarios, the damage in the city of Al Hoceima for the RPS2000 scenario is more significant than the one obtained for Eurocode 8. In Imzouren, the damage values are more spread out according to the first scenario, ranging from 1.61 to 3.76, compared with the values of the Eurocode 8 scenario, which are between 1.94 and 3.25. According to the first scenario, the percentage of buildings having a "collapse" damage state is equal to 25% of the buildings inspected, all built on S4 soft soils. These results call into question the reliability of the spectrum taken from the national seismic code [17], which, as the conclusions of this study suggest, should be reviewed. The Eurocode 8 risk scenario results are closer to the results and observations reported by Cherif et al. [37–39] for these two cities.
<figure><img src="image_25.png" /><figcaption>**Figure 14.** Damage distribution in the city of Al Hoceima following (a) the RPS2000 risk scenario and (b) the Eurocode 8 risk scenario.</figcaption></figure>
<figure><img src="image_26.png" /><figcaption>**Figure 15.** Damage distribution in the city of Imzouren following (a) the RPS2000 risk scenario and (b) the Eurocode 8 risk scenario.</figcaption></figure>
## 5. Final Remarks

Seismic risk assessment based on the Capacity Spectrum method was used in this study and applied to the buildings of Al Hoceima and Imzouren, the most important towns of the central Rif region in Morocco. The region has experienced devastating earthquakes in the recent past and multiple studies have been carried out to assess the vulnerability and seismic risk of buildings in the region. This study falls within this framework, where importance is given to the development of building capacity spectra for the estimation of seismic performance. Two probabilistic risk scenarios were considered, using 5% damped elastic response spectra. The first scenario is based on the response spectrum taken from the national seismic code RPS2000, while the second scenario is based on the response spectrum model of Eurocode 8, using the seismic parameters and properties of the central Rif region. Response spectra were generated for four types of soil, ranging from rocky soil to soft soil.

The capacity spectra were developed for 36 building models, which were designed to represent the two cities' most common types of structures. The model determination is related to the seismic code, the construction period, and the number of floors. Capacity spectra were generated by performing a non-linear pushover analysis. Fragility curves were estimated for the building models, following capacity spectra. They are used to estimate the damage induced following a chosen seismic action. By specifying the spectral displacement of the structure, one can build damage probability matrices and calculate the average damage index. Performance points for each model and each suggested risk scenario are calculated afterwards. The adopted method was applied to the buildings of the city of Al Hoceima and Imzouren. A database of 2746 existing buildings in the region was used to carry out this study. The damage estimate was made by calculating the average damage index.

The most important results show that the damage is more significant in Imzouren than in Al Hoceima, mainly because of the seismic hazard and the soil typology. The distribution of damage in the two cities is quite similar to previous studies carried out in the same cities. Buildings with high damage indices are mainly high vulnerability buildings in the city of Al Hoceima and buildings that do not comply with earthquake regulations and are built on loose soil in Imzouren. By reflecting on the risk scenarios considered, the credibility of the response spectra drawn from the RPS2000 is questioned, given the extreme damage estimated, especially in Imzouren, and which do not reflect the results found in studies carried out previously. The models of the response spectra derived from the Eurocode and adapted to the seismic properties of the region prove to be a reliable model for future studies in the region.

The results related to the capacity spectra and fragility curves, as well as the generated seismic risk maps, constitute a valuable source of information that can be used to support effective seismic risk mitigation and management actions in the central Rif region. However, it is part of a broader vision for a better understanding of the seismic risk in the Rif region. The limitations of the work will be explored in future projects that will revolve around the following points:

* Expanding the existing databases and planning investigations in the rural areas of the Rif chain.

* Developing response spectra of target regions using seismograph recordings.

* Estimating site effects using shear wave velocity vs30.

* Considering new building models such as shear wall structures or mixed structures, which became widespread after the seismic regulations.

**Author Contributions:** Conceptualization, S.-E.C., M.C. and T.M.F.; methodology, S.-E.C.; validation A.C., M.O.M. and A.E.O.; formal analysis, A.C., M.O.M. and A.E.O.; investigation, S.-E.C.; resources, S.-E.C.; data curation, M.O.M.; writing—original draft preparation, S.-E.C.; writing—review and editing, M.C. and T.M.F.; visualization, S.-E.C.; supervision, M.C. and T.M.F.; project administration, M.C. All authors have read and agreed to the published version of the manuscript.
**Funding:** This research received no external funding.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Not applicable.

**Conflicts of Interest:** The authors declare no conflict of interest.

## Appendix A

<table><caption>Table A1. Yield and ultimate capacity for buildings with seismic code.</caption><thead><tr><th rowspan="2">Buildings</th><th colspan="4">Construction Period I</th><th colspan="4">Construction Period II</th><th colspan="4">Construction Period III</th></tr><tr><th colspan="2">Yield Capacity</th><th colspan="2">Ultimate Capacity</th><th colspan="2">Yield Capacity</th><th colspan="2">Ultimate Capacity</th><th colspan="2">Yield Capacity</th><th colspan="2">Ultimate Capacity</th></tr><tr><th></th><th>Dy (cm)</th><th>Ay (g)</th><th>Du (cm)</th><th>Au (g)</th><th>Dy (cm)</th><th>Ay (g)</th><th>Du (cm)</th><th>Au(g)</th><th>Dy (cm)</th><th>Ay (g)</th><th>Du (cm)</th><th>Au(g)</th></tr></thead><tbody><tr><td>1 floor</td><td>0.409</td><td>0.297</td><td>8.934</td><td>0.422</td><td>0.578</td><td>0.442</td><td>8.856</td><td>0.521</td><td>0.63</td><td>0.484</td><td>9.25</td><td>0.969</td></tr><tr><td>2 floors</td><td>0.689</td><td>0.234</td><td>9.078</td><td>0.274</td><td>0.874</td><td>0.298</td><td>9.145</td><td>0.366</td><td>1.02</td><td>0.352</td><td>9.82</td><td>0.588</td></tr><tr><td>3 floors</td><td>1.403</td><td>0.250</td><td>10.178</td><td>0.356</td><td>1.494</td><td>0.275</td><td>12.833</td><td>0.501</td><td>1.56</td><td>0.303</td><td>11.64</td><td>0.609</td></tr><tr><td>4 floors</td><td>1.193</td><td>0.155</td><td>12.361</td><td>0.277</td><td>2.186</td><td>0.243</td><td>16.165</td><td>0.366</td><td>2.03</td><td>0.245</td><td>15.01</td><td>0.430</td></tr><tr><td>5 floors</td><td>1.536</td><td>0.123</td><td>16.673</td><td>0.233</td><td>1.773</td><td>0.159</td><td>17.725</td><td>0.303</td><td>2.04</td><td>0.192</td><td>17.56</td><td>0.333</td></tr><tr><td>6 floors</td><td>1.507</td><td>0.104</td><td>20.829</td><td>0.196</td><td>2.404</td><td>0.132</td><td>20.601</td><td>0.237</td><td>2.49</td><td>0.147</td><td>20.69</td><td>0.270</td></tr></tbody></table>

<table><caption>Table A2. Yield and ultimate capacity for buildings without seismic code.</caption><thead><tr><th rowspan="2">Buildings</th><th colspan="4">Construction Period I</th><th colspan="4">Construction Period II</th><th colspan="4">Construction Period III</th></tr><tr><th colspan="2">Yield Capacity</th><th colspan="2">Ultimate Capacity</th><th colspan="2">Yield Capacity</th><th colspan="2">Ultimate Capacity</th><th colspan="2">Yield Capacity</th><th colspan="2">Ultimate Capacity</th></tr><tr><th></th><th>Dy (cm)</th><th>Ay (g)</th><th>Du (cm)</th><th>Au (g)</th><th>Dy (cm)</th><th>Ay (g)</th><th>Du (cm)</th><th>Au(g)</th><th>Dy (cm)</th><th>Ay (g)</th><th>Du (cm)</th><th>Au(g)</th></tr></thead><tbody><tr><td>1 floor</td><td>0.439</td><td>0.145</td><td>9.100</td><td>0.164</td><td>0.559</td><td>0.188</td><td>9.057</td><td>0.280</td><td>0.64</td><td>0.217</td><td>9.14</td><td>0.260</td></tr><tr><td>2 floors</td><td>0.650</td><td>0.100</td><td>8.689</td><td>0.115</td><td>0.869</td><td>0.134</td><td>9.025</td><td>0.174</td><td>0.97</td><td>0.150</td><td>8.90</td><td>0.175</td></tr><tr><td>3 floors</td><td>1.446</td><td>0.143</td><td>9.259</td><td>0.161</td><td>1.598</td><td>0.165</td><td>10.052</td><td>0.227</td><td>1.48</td><td>0.152</td><td>11.00</td><td>0.284</td></tr><tr><td>4 floors</td><td>1.832</td><td>0.128</td><td>11.064</td><td>0.138</td><td>2.138</td><td>0.154</td><td>11.325</td><td>0.174</td><td>2.28</td><td>0.166</td><td>10.88</td><td>0.223</td></tr><tr><td>5 floors</td><td>1.923</td><td>0.101</td><td>14.336</td><td>0.142</td><td>2.608</td><td>0.139</td><td>13.251</td><td>0.159</td><td>2.56</td><td>0.143</td><td>13.52</td><td>0.178</td></tr><tr><td>6 floors</td><td>2.279</td><td>0.086</td><td>18.493</td><td>0.129</td><td>2.880</td><td>0.117</td><td>18.435</td><td>0.159</td><td>2.94</td><td>0.128</td><td>20.80</td><td>0.186</td></tr></tbody></table>

<table><caption>Table A3. Parameters characterizing the fragility curves for the building models.</caption><thead><tr><th rowspan="2"></th><th colspan="8">Construction Period I</th></tr><tr><th>Number of floors</th><th>Sd1 (cm)</th><th>βSD1</th><th>Sd2 (cm)</th><th>βSD2</th><th>Sd3 (cm)</th><th>βSD3</th><th>Sd4 (cm)</th><th>βSD4</th></tr></thead><tbody><tr><td rowspan="6">Buildings without seismic code</td><td>1</td><td>0.31</td><td>0.46</td><td>0.44</td><td>0.75</td><td>2.60</td><td>1.31</td><td>9.10</td><td>1.67</td></tr><tr><td>2</td><td>0.46</td><td>0.43</td><td>0.65</td><td>0.67</td><td>2.66</td><td>1.14</td><td>8.69</td><td>1.45</td></tr><tr><td>3</td><td>1.01</td><td>0.38</td><td>1.45</td><td>0.53</td><td>3.40</td><td>0.84</td><td>9.26</td><td>1.08</td></tr><tr><td>4</td><td>1.28</td><td>0.38</td><td>1.83</td><td>0.52</td><td>4.14</td><td>0.82</td><td>11.06</td><td>1.05</td></tr><tr><td>5</td><td>1.35</td><td>0.39</td><td>1.92</td><td>0.56</td><td>5.03</td><td>0.90</td><td>14.34</td><td>1.15</td></tr><tr><td>6</td><td>1.60</td><td>0.40</td><td>2.28</td><td>0.58</td><td>6.33</td><td>0.94</td><td>18.49</td><td>1.20</td></tr><tr><td rowspan="6"></td><th colspan="8">Construction Period II</th></tr><tr><th>Number of floors</th><th>Sd1 (cm)</th><th>βSD1</th><th>Sd2 (cm)</th><th>βSD2</th><th>Sd3 (cm)</th><th>βSD3</th><th>Sd4 (cm)</th><th>βSD4</th></tr><tr><td>1</td><td>0.39</td><td>0.44</td><td>0.56</td><td>0.70</td><td>2.68</td><td>1.21</td><td>9.06</td><td>1.54</td></tr><tr><td>2</td><td>0.61</td><td>0.41</td><td>0.87</td><td>0.62</td><td>2.91</td><td>1.04</td><td>9.03</td><td>1.32</td></tr><tr><td>3</td><td>1.12</td><td>0.38</td><td>1.60</td><td>0.53</td><td>3.71</td><td>0.84</td><td>10.05</td><td>1.07</td></tr><tr><td>4</td><td>1.50</td><td>0.37</td><td>2.14</td><td>0.50</td><td>4.43</td><td>0.77</td><td>11.32</td><td>0.98</td></tr><tr><td>5</td><td>1.83</td><td>0.36</td><td>2.61</td><td>0.49</td><td>5.27</td><td>0.75</td><td>13.25</td><td>0.96</td></tr><tr><td></td><td>6</td><td>2.02</td><td>0.38</td><td>2.88</td><td>0.53</td><td>6.77</td><td>0.84</td><td>18.43</td><td>1.08</td></tr></tbody></table>
<table><caption>Table A3. Cont.</caption><thead><tr><th colspan="9">Construction Period III</th></tr><tr><th>Number of floors</th><th>Sd1 (cm)</th><th>βSD1</th><th>Sd2 (cm)</th><th>βSD2</th><th>Sd3 (cm)</th><th>βSD3</th><th>Sd4 (cm)</th><th>βSD4</th></tr></thead><tbody><tr><td>1</td><td>0.45</td><td>0.44</td><td>0.64</td><td>0.68</td><td>2.77</td><td>1.16</td><td>9.14</td><td>1.48</td></tr><tr><td>2</td><td>0.68</td><td>0.41</td><td>0.97</td><td>0.60</td><td>2.95</td><td>0.99</td><td>8.90</td><td>1.26</td></tr><tr><td>3</td><td>1.03</td><td>0.39</td><td>1.48</td><td>0.56</td><td>3.86</td><td>0.90</td><td>11.00</td><td>1.15</td></tr><tr><td>4</td><td>1.60</td><td>0.36</td><td>2.28</td><td>0.48</td><td>4.43</td><td>0.72</td><td>10.88</td><td>0.93</td></tr><tr><td>5</td><td>1.80</td><td>0.37</td><td>2.56</td><td>0.50</td><td>5.30</td><td>0.76</td><td>13.52</td><td>0.98</td></tr><tr><td>6</td><td>2.06</td><td>0.39</td><td>2.94</td><td>0.55</td><td>7.41</td><td>0.88</td><td>20.80</td><td>1.13</td></tr></tbody><thead><tr><th colspan="9">Construction Period I</th></tr><tr><th>Number of floors</th><th>Sd1 (cm)</th><th>βSD1</th><th>Sd2 (cm)</th><th>βSD2</th><th>Sd3 (cm)</th><th>βSD3</th><th>Sd4 (cm)</th><th>βSD4</th></tr></thead><tbody><tr><td>1</td><td>0.29</td><td>0.47</td><td>0.41</td><td>0.75</td><td>2.54</td><td>1.33</td><td>8.93</td><td>1.69</td></tr><tr><td>2</td><td>0.48</td><td>0.43</td><td>0.69</td><td>0.66</td><td>2.79</td><td>1.13</td><td>9.08</td><td>1.44</td></tr><tr><td>3</td><td>0.98</td><td>0.39</td><td>1.40</td><td>0.56</td><td>3.60</td><td>0.89</td><td>10.18</td><td>1.14</td></tr><tr><td>4</td><td>0.84</td><td>0.41</td><td>1.19</td><td>0.62</td><td>3.99</td><td>1.04</td><td>12.36</td><td>1.32</td></tr><tr><td>5</td><td>1.08</td><td>0.42</td><td>1.54</td><td>0.63</td><td>5.32</td><td>1.05</td><td>16.67</td><td>1.34</td></tr><tr><td>6</td><td>1.05</td><td>0.43</td><td>1.51</td><td>0.67</td><td>6.34</td><td>1.15</td><td>20.83</td><td>1.46</td></tr></tbody><thead><tr><th colspan="9">Construction Period II</th></tr><tr><th>Number of floors</th><th>Sd1 (cm)</th><th>βSD1</th><th>Sd2 (cm)</th><th>βSD2</th><th>Sd3 (cm)</th><th>βSD3</th><th>Sd4 (cm)</th><th>βSD4</th></tr></thead><tbody><tr><td>1</td><td>0.40</td><td>0.44</td><td>0.58</td><td>0.69</td><td>2.65</td><td>1.19</td><td>8.86</td><td>1.51</td></tr><tr><td>2</td><td>0.61</td><td>0.41</td><td>0.87</td><td>0.62</td><td>2.94</td><td>1.04</td><td>9.14</td><td>1.32</td></tr><tr><td>3</td><td>1.05</td><td>0.40</td><td>1.49</td><td>0.59</td><td>4.33</td><td>0.96</td><td>12.83</td><td>1.23</td></tr><tr><td>4</td><td>1.53</td><td>0.39</td><td>2.19</td><td>0.56</td><td>5.68</td><td>0.90</td><td>16.17</td><td>1.15</td></tr><tr><td>5</td><td>1.24</td><td>0.41</td><td>1.77</td><td>0.61</td><td>5.76</td><td>1.02</td><td>17.72</td><td>1.30</td></tr><tr><td>6</td><td>1.68</td><td>0.40</td><td>2.40</td><td>0.59</td><td>6.95</td><td>0.96</td><td>20.60</td><td>1.22</td></tr></tbody><thead><tr><th colspan="9">Construction Period III</th></tr><tr><th>Number of floors</th><th>Sd1 (cm)</th><th>βSD1</th><th>Sd2 (cm)</th><th>βSD2</th><th>Sd3 (cm)</th><th>βSD3</th><th>Sd4 (cm)</th><th>βSD4</th></tr></thead><tbody><tr><td>1</td><td>0.44</td><td>0.44</td><td>0.63</td><td>0.68</td><td>2.79</td><td>1.17</td><td>9.25</td><td>1.49</td></tr><tr><td>2</td><td>0.72</td><td>0.41</td><td>1.02</td><td>0.61</td><td>3.22</td><td>1.00</td><td>9.82</td><td>1.28</td></tr><tr><td>3</td><td>1.09</td><td>0.39</td><td>1.56</td><td>0.56</td><td>4.08</td><td>0.90</td><td>11.64</td><td>1.15</td></tr><tr><td>4</td><td>1.42</td><td>0.39</td><td>2.03</td><td>0.56</td><td>5.28</td><td>0.90</td><td>15.01</td><td>1.15</td></tr><tr><td>5</td><td>1.43</td><td>0.40</td><td>2.04</td><td>0.59</td><td>5.92</td><td>0.96</td><td>17.56</td><td>1.23</td></tr><tr><td>6</td><td>1.75</td><td>0.40</td><td>2.49</td><td>0.58</td><td>7.04</td><td>0.95</td><td>20.69</td><td>1.21</td></tr></tbody></table>

<table><caption>Table A4. Spectral displacement (cm) of the performance points for the RC building models according to both probabilistic scenarios (RPS2000 and Eurocode 8).</caption><thead><tr><th rowspan="2">Number of stories</th><th rowspan="2">Soil Type</th><th colspan="3">Construction Period I</th><th colspan="3">Construction Period II</th><th colspan="3">Construction Period III</th></tr><tr><th>RPS</th><th>E8 ALH</th><th>E8 IMZ</th><th>RPS</th><th>E8 ALH</th><th>E8 IMZ</th><th>RPS</th><th>E8 ALH</th><th>E8 IMZ</th></tr></thead><tbody><tr><td rowspan="4">1</td><td>S1</td><td>2.80</td><td>1.70</td><td>2.40</td><td>2.00</td><td>1.60</td><td>2.10</td><td>1.80</td><td>1.50</td><td>2.00</td></tr><tr><td>S2</td><td>4.50</td><td>3.00</td><td>4.20</td><td>2.90</td><td>2.50</td><td>3.30</td><td>2.50</td><td>2.30</td><td>3.10</td></tr><tr><td>S3</td><td>6.70</td><td>5.10</td><td>6.40</td><td>4.10</td><td>4.00</td><td>5.40</td><td>3.50</td><td>3.70</td><td>5.10</td></tr><tr><td>S4</td><td></td><td>5.60</td><td>6.90</td><td>6.80</td><td>4.50</td><td>6.00</td><td>6.30</td><td>4.10</td><td>5.60</td></tr><tr><td rowspan="4">2</td><td>S1</td><td>5.60</td><td>2.50</td><td>3.40</td><td>3.40</td><td>2.20</td><td>2.90</td><td>3.10</td><td>2.10</td><td>2.70</td></tr><tr><td>S2</td><td>8.50</td><td>4.00</td><td>4.80</td><td>5.20</td><td>3.40</td><td>4.70</td><td>4.70</td><td>3.30</td><td>4.40</td></tr><tr><td>S3</td><td></td><td>5.30</td><td>6.60</td><td>7.30</td><td>5.50</td><td>6.80</td><td>6.70</td><td>5.30</td><td>6.60</td></tr><tr><td>S4</td><td></td><td>5.70</td><td>7.00</td><td>6.00</td><td>7.30</td><td>6.00</td><td>7.30</td><td>5.80</td><td>7.10</td></tr><tr><td rowspan="4">3</td><td>S1</td><td>4.20</td><td>2.70</td><td>3.40</td><td>3.40</td><td>2.40</td><td>3.00</td><td>3.10</td><td>2.40</td><td>2.90</td></tr><tr><td>S2</td><td>5.80</td><td>4.00</td><td>5.20</td><td>4.60</td><td>3.60</td><td>4.50</td><td>4.10</td><td>3.40</td><td>4.30</td></tr><tr><td>S3</td><td>7.90</td><td>5.90</td><td>7.00</td><td>6.10</td><td>5.30</td><td>7.00</td><td>5.40</td><td>5.00</td><td>6.50</td></tr><tr><td>S4</td><td></td><td>6.30</td><td>7.40</td><td>9.80</td><td>5.90</td><td>7.60</td><td>8.80</td><td>5.40</td><td>7.20</td></tr></tbody></table>
<table><caption>Table A4. Cont.</caption><thead><tr><th rowspan="2"></th><th rowspan="2">Number of stories</th><th rowspan="2">Soil Type</th><th colspan="3">Construction Period I</th><th colspan="3">Construction Period II</th><th colspan="3">Construction Period III</th></tr><tr><th>RPS</th><th>E8 ALH</th><th>E8 IMZ</th><th>RPS</th><th>E8 ALH</th><th>E8 IMZ</th><th>RPS</th><th>E8 ALH</th><th>E8 IMZ</th></tr></thead><tbody><tr><td rowspan="12">Buildings without seismic code</td><td rowspan="4">4</td><td>S1</td><td>5.20</td><td>3.20</td><td>4.00</td><td>4.30</td><td>2.90</td><td>3.60</td><td>4.00</td><td>2.90</td><td>3.50</td></tr><tr><td>S2</td><td>7.20</td><td>4.70</td><td>5.60</td><td>5.80</td><td>4.20</td><td>5.30</td><td>5.30</td><td>4.00</td><td>5.00</td></tr><tr><td>S3</td><td>9.90</td><td>6.20</td><td>7.10</td><td>7.70</td><td>6.10</td><td>7.20</td><td>6.80</td><td>5.80</td><td>7.40</td></tr><tr><td>S4</td><td></td><td>6.50</td><td>7.60</td><td></td><td>6.40</td><td>7.50</td><td>10.60</td><td>6.40</td><td>7.90</td></tr><tr><td rowspan="4">5</td><td>S1</td><td>6.90</td><td>3.70</td><td>4.40</td><td>5.40</td><td>3.40</td><td>4.20</td><td>4.90</td><td>3.30</td><td>4.00</td></tr><tr><td>S2</td><td>9.40</td><td>4.90</td><td>5.80</td><td>7.20</td><td>4.90</td><td>5.80</td><td>6.50</td><td>4.70</td><td>5.80</td></tr><tr><td>S3</td><td>12.50</td><td>6.40</td><td>7.60</td><td>9.60</td><td>6.40</td><td>7.50</td><td>8.60</td><td>6.40</td><td>7.40</td></tr><tr><td>S4</td><td></td><td>6.90</td><td>8.10</td><td></td><td>6.80</td><td>7.90</td><td></td><td>6.70</td><td>7.70</td></tr><tr><td rowspan="4">6</td><td>S1</td><td>8.70</td><td>3.60</td><td>4.40</td><td>6.90</td><td>3.80</td><td>4.50</td><td>6.10</td><td>3.80</td><td>4.60</td></tr><tr><td>S2</td><td>12.00</td><td>5.00</td><td>5.90</td><td>9.10</td><td>5.10</td><td>5.90</td><td>8.10</td><td>5.10</td><td>6.00</td></tr><tr><td>S3</td><td></td><td>6.60</td><td>7.70</td><td>11.90</td><td>6.60</td><td>7.90</td><td>10.50</td><td>6.60</td><td>7.80</td></tr><tr><td>S4</td><td></td><td>6.90</td><td>8.20</td><td></td><td>7.00</td><td>8.30</td><td>16.70</td><td>7.00</td><td>8.30</td></tr></tbody></table><table><thead><tr><th rowspan="2"></th><th rowspan="2">Number of stories</th><th rowspan="2">Soil Type</th><th colspan="3">Construction Period I</th><th colspan="3">Construction Period II</th><th colspan="3">Construction Period III</th></tr><tr><th>RPS</th><th>E8 ALH</th><th>E8 IMZ</th><th>RPS</th><th>E8 ALH</th><th>E8 IMZ</th><th>RPS</th><th>E8 ALH</th><th>E8 IMZ</th></tr></thead><tbody><tr><td rowspan="12">Buildings with seismic code</td><td rowspan="4">1</td><td>S1</td><td>0.62</td><td>0.76</td><td>0.87</td><td>0.60</td><td>0.78</td><td>0.92</td><td>0.60</td><td>0.79</td><td>0.94</td></tr><tr><td>S2</td><td>0.75</td><td>1.00</td><td>1.30</td><td>0.73</td><td>1.00</td><td>1.20</td><td>0.72</td><td>1.00</td><td>1.20</td></tr><tr><td>S3</td><td>0.88</td><td>1.60</td><td>2.10</td><td>0.86</td><td>1.40</td><td>1.80</td><td>0.85</td><td>1.40</td><td>1.70</td></tr><tr><td>S4</td><td>1.70</td><td>1.80</td><td>2.40</td><td>1.20</td><td>1.50</td><td>1.90</td><td>1.10</td><td>1.50</td><td>1.80</td></tr><tr><td rowspan="4">2</td><td>S1</td><td>1.10</td><td>1.10</td><td>1.30</td><td>1.10</td><td>1.10</td><td>1.30</td><td>1.10</td><td>1.20</td><td>1.30</td></tr><tr><td>S2</td><td>1.40</td><td>1.50</td><td>1.80</td><td>1.20</td><td>1.40</td><td>1.70</td><td>1.30</td><td>1.40</td><td>1.60</td></tr><tr><td>S3</td><td>1.70</td><td>2.10</td><td>2.80</td><td>1.50</td><td>1.90</td><td>2.40</td><td>1.40</td><td>1.80</td><td>2.30</td></tr><tr><td>S4</td><td>2.70</td><td>2.30</td><td>3.20</td><td>2.20</td><td>2.10</td><td>2.70</td><td>1.90</td><td>2.00</td><td>2.50</td></tr><tr><td rowspan="4">3</td><td>S1</td><td>1.60</td><td>1.40</td><td>1.70</td><td>1.60</td><td>1.50</td><td>1.70</td><td>1.60</td><td>1.50</td><td>1.80</td></tr><tr><td>S2</td><td>2.00</td><td>2.00</td><td>2.50</td><td>1.90</td><td>2.00</td><td>2.40</td><td>1.90</td><td>2.00</td><td>2.30</td></tr><tr><td>S3</td><td>2.50</td><td>2.80</td><td>3.60</td><td>2.30</td><td>2.70</td><td>3.30</td><td>2.20</td><td>2.60</td><td>3.20</td></tr><tr><td>S4</td><td>3.70</td><td>3.10</td><td>3.90</td><td>3.30</td><td>2.90</td><td>3.60</td><td>3.00</td><td>2.80</td><td>3.40</td></tr><tr><td rowspan="4">4</td><td>S1</td><td>2.40</td><td>1.80</td><td>2.20</td><td>2.10</td><td>1.70</td><td>2.10</td><td>2.00</td><td>1.80</td><td>2.10</td></tr><tr><td>S2</td><td>3.20</td><td>2.60</td><td>3.40</td><td>2.60</td><td>2.40</td><td>2.90</td><td>2.50</td><td>2.30</td><td>2.80</td></tr><tr><td>S3</td><td>4.00</td><td>3.90</td><td>4.90</td><td>3.30</td><td>3.40</td><td>4.20</td><td>3.00</td><td>3.20</td><td>4.00</td></tr><tr><td>S4</td><td>5.90</td><td>4.20</td><td>5.30</td><td>4.80</td><td>3.60</td><td>4.70</td><td>4.20</td><td>3.50</td><td>4.30</td></tr><tr><td rowspan="4">5</td><td>S1</td><td>3.20</td><td>2.10</td><td>2.60</td><td>2.70</td><td>2.00</td><td>2.50</td><td>2.50</td><td>2.00</td><td>2.40</td></tr><tr><td>S2</td><td>4.40</td><td>3.10</td><td>4.10</td><td>3.40</td><td>2.80</td><td>3.50</td><td>3.10</td><td>2.70</td><td>3.30</td></tr><tr><td>S3</td><td>5.60</td><td>4.70</td><td>5.90</td><td>4.30</td><td>4.00</td><td>5.20</td><td>3.90</td><td>3.80</td><td>4.70</td></tr><tr><td>S4</td><td>8.20</td><td>5.20</td><td>6.40</td><td>6.40</td><td>4.40</td><td>5.60</td><td>5.70</td><td>4.10</td><td>5.20</td></tr><tr><td rowspan="4">6</td><td>S1</td><td>4.10</td><td>2.40</td><td>3.10</td><td>3.80</td><td>2.50</td><td>3.00</td><td>3.50</td><td>2.40</td><td>3.00</td></tr><tr><td>S2</td><td>5.50</td><td>3.60</td><td>4.30</td><td>4.90</td><td>3.50</td><td>4.40</td><td>4.50</td><td>3.40</td><td>4.20</td></tr><tr><td>S3</td><td>7.10</td><td>4.90</td><td>5.80</td><td>6.30</td><td>5.00</td><td>6.00</td><td>5.60</td><td>4.80</td><td>6.00</td></tr><tr><td>S4</td><td>10.50</td><td>5.20</td><td>6.20</td><td>9.30</td><td>5.30</td><td>6.40</td><td>8.20</td><td>5.20</td><td>6.40</td></tr></tbody></table>
<table><caption>**Table A5.** Damage index $DS_m$ values for the RC building models according to both probabilistic scenarios (RPS2000 and Eurocode 8).</caption><thead><tr><th rowspan="2">Number of stories</th><th rowspan="2">Soil Type</th><th colspan="3">Construction Period I</th><th colspan="3">Construction Period II</th><th colspan="3">Construction Period III</th></tr><tr><th>RPS</th><th>E8 ALH</th><th>E8 IMZ</th><th>RPS</th><th>E8 ALH</th><th>E8 IMZ</th><th>RPS</th><th>E8 ALH</th><th>E8 IMZ</th></tr></thead><tbody><tr><td rowspan="4">1</td><td>S1</td><td>2.76</td><td>2.49</td><td>2.68</td><td>2.53</td><td>2.40</td><td>2.56</td><td>2.43</td><td>2.30</td><td>2.49</td></tr><tr><td>S2</td><td>3.00</td><td>2.79</td><td>2.96</td><td>2.75</td><td>2.66</td><td>2.82</td><td>2.63</td><td>2.58</td><td>2.76</td></tr><tr><td>S3</td><td>3.19</td><td>3.06</td><td>3.17</td><td>2.94</td><td>2.92</td><td>3.09</td><td>2.83</td><td>2.86</td><td>3.05</td></tr><tr><td>S4</td><td>3.39</td><td>3.11</td><td>3.20</td><td>3.20</td><td>2.99</td><td>3.14</td><td>3.16</td><td>2.92</td><td>3.10</td></tr><tr><td rowspan="4">2</td><td>S1</td><td>3.12</td><td>2.65</td><td>2.84</td><td>2.78</td><td>2.47</td><td>2.67</td><td>2.69</td><td>2.39</td><td>2.59</td></tr><tr><td>S2</td><td>3.34</td><td>2.93</td><td>3.04</td><td>3.05</td><td>2.78</td><td>2.99</td><td>2.98</td><td>2.74</td><td>2.94</td></tr><tr><td>S3</td><td>3.51</td><td>3.09</td><td>3.21</td><td>3.25</td><td>3.08</td><td>3.21</td><td>3.21</td><td>3.06</td><td>3.20</td></tr><tr><td>S4</td><td>3.60</td><td>3.13</td><td>3.24</td><td>3.53</td><td>3.14</td><td>3.25</td><td>3.51</td><td>3.12</td><td>3.24</td></tr><tr><td rowspan="4">3</td><td>S1</td><td>2.81</td><td>2.39</td><td>2.62</td><td>2.53</td><td>2.15</td><td>2.41</td><td>2.45</td><td>2.18</td><td>2.38</td></tr><tr><td>S2</td><td>3.06</td><td>2.77</td><td>2.98</td><td>2.81</td><td>2.59</td><td>2.79</td><td>2.69</td><td>2.53</td><td>2.73</td></tr><tr><td>S3</td><td>3.28</td><td>3.08</td><td>3.20</td><td>3.04</td><td>2.93</td><td>3.14</td><td>2.90</td><td>2.85</td><td>3.04</td></tr><tr><td>S4</td><td>3.58</td><td>3.13</td><td>3.24</td><td>3.37</td><td>3.01</td><td>3.20</td><td>3.24</td><td>2.90</td><td>3.11</td></tr><tr><td rowspan="4">4</td><td>S1</td><td>2.82</td><td>2.34</td><td>2.58</td><td>2.56</td><td>2.07</td><td>2.36</td><td>2.46</td><td>2.00</td><td>2.28</td></tr><tr><td>S2</td><td>3.09</td><td>2.73</td><td>2.89</td><td>2.86</td><td>2.54</td><td>2.78</td><td>2.78</td><td>2.46</td><td>2.72</td></tr><tr><td>S3</td><td>3.31</td><td>2.97</td><td>3.08</td><td>3.11</td><td>2.91</td><td>3.05</td><td>3.02</td><td>2.87</td><td>3.09</td></tr><tr><td>S4</td><td>3.63</td><td>3.01</td><td>3.13</td><td>3.49</td><td>2.95</td><td>3.08</td><td>3.37</td><td>2.96</td><td>3.15</td></tr><tr><td rowspan="4">5</td><td>S1</td><td>2.89</td><td>2.36</td><td>2.52</td><td>2.62</td><td>2.02</td><td>2.32</td><td>2.51</td><td>1.99</td><td>2.26</td></tr><tr><td>S2</td><td>3.11</td><td>2.62</td><td>2.75</td><td>2.90</td><td>2.51</td><td>2.69</td><td>2.80</td><td>2.46</td><td>2.69</td></tr><tr><td>S3</td><td>3.30</td><td>2.83</td><td>2.96</td><td>3.15</td><td>2.79</td><td>2.94</td><td>3.05</td><td>2.79</td><td>2.92</td></tr><tr><td>S4</td><td>3.52</td><td>2.89</td><td>3.01</td><td>3.50</td><td>2.85</td><td>2.99</td><td>3.40</td><td>2.83</td><td>2.96</td></tr><tr><td rowspan="4">6</td><td>S1</td><td>2.89</td><td>2.13</td><td>2.33</td><td>2.64</td><td>1.97</td><td>2.19</td><td>2.46</td><td>1.91</td><td>2.16</td></tr><tr><td>S2</td><td>3.11</td><td>2.45</td><td>2.59</td><td>2.88</td><td>2.34</td><td>2.49</td><td>2.71</td><td>2.27</td><td>2.44</td></tr><tr><td>S3</td><td>3.45</td><td>2.68</td><td>2.80</td><td>3.09</td><td>2.60</td><td>2.76</td><td>2.92</td><td>2.53</td><td>2.68</td></tr><tr><td>S4</td><td>3.76</td><td>2.71</td><td>2.84</td><td>3.51</td><td>2.65</td><td>2.80</td><td>3.24</td><td>2.58</td><td>2.73</td></tr><tr><th colspan="11" style="text-align: center;">Construction Period I</th></tr><tr><th colspan="11" style="text-align: center;">Construction Period II</th></tr><tr><th colspan="11" style="text-align: center;">Construction Period III</th></tr><tr><th>Number of stories</th><th>Soil Type</th><th>RPS</th><th>E8 ALH</th><th>E8 IMZ</th><th>RPS</th><th>E8 ALH</th><th>E8 IMZ</th><th>RPS</th><th>E8 ALH</th><th>E8 IMZ</th></tr><tr><td rowspan="4">1</td><td>S1</td><td>1.87</td><td>2.03</td><td>2.13</td><td>1.48</td><td>1.81</td><td>1.97</td><td>1.35</td><td>1.73</td><td>1.91</td></tr><tr><td>S2</td><td>2.02</td><td>2.22</td><td>2.37</td><td>1.73</td><td>2.05</td><td>2.19</td><td>1.61</td><td>1.98</td><td>2.14</td></tr><tr><td>S3</td><td>2.13</td><td>2.48</td><td>2.62</td><td>1.90</td><td>2.31</td><td>2.47</td><td>1.81</td><td>2.25</td><td>2.39</td></tr><tr><td>S4</td><td>2.52</td><td>2.54</td><td>2.69</td><td>2.19</td><td>2.35</td><td>2.50</td><td>2.06</td><td>2.30</td><td>2.43</td></tr><tr><td rowspan="4">2</td><td>S1</td><td>2.01</td><td>2.01</td><td>2.16</td><td>1.79</td><td>1.79</td><td>1.99</td><td>1.59</td><td>1.71</td><td>1.82</td></tr><tr><td>S2</td><td>2.22</td><td>2.27</td><td>2.41</td><td>1.90</td><td>2.07</td><td>2.25</td><td>1.82</td><td>1.91</td><td>2.06</td></tr><tr><td>S3</td><td>2.37</td><td>2.51</td><td>2.69</td><td>2.14</td><td>2.35</td><td>2.53</td><td>1.91</td><td>2.18</td><td>2.40</td></tr><tr><td>S4</td><td>2.67</td><td>2.57</td><td>2.77</td><td>2.46</td><td>2.42</td><td>2.61</td><td>2.24</td><td>2.28</td><td>2.47</td></tr><tr><td rowspan="4">3</td><td>S1</td><td>1.72</td><td>1.50</td><td>1.81</td><td>1.60</td><td>1.49</td><td>1.69</td><td>1.55</td><td>1.44</td><td>1.74</td></tr><tr><td>S2</td><td>2.04</td><td>2.04</td><td>2.29</td><td>1.85</td><td>1.91</td><td>2.13</td><td>1.82</td><td>1.89</td><td>2.07</td></tr><tr><td>S3</td><td>2.29</td><td>2.41</td><td>2.64</td><td>2.08</td><td>2.25</td><td>2.43</td><td>2.02</td><td>2.21</td><td>2.42</td></tr><tr><td>S4</td><td>2.66</td><td>2.50</td><td>2.70</td><td>2.43</td><td>2.32</td><td>2.51</td><td>2.36</td><td>2.29</td><td>2.48</td></tr><tr><td rowspan="4">4</td><td>S1</td><td>2.28</td><td>2.01</td><td>2.21</td><td>1.44</td><td>1.05</td><td>1.44</td><td>1.48</td><td>1.29</td><td>1.56</td></tr><tr><td>S2</td><td>2.51</td><td>2.35</td><td>2.56</td><td>1.78</td><td>1.66</td><td>1.94</td><td>1.83</td><td>1.71</td><td>1.99</td></tr><tr><td>S3</td><td>2.67</td><td>2.65</td><td>2.81</td><td>2.10</td><td>2.14</td><td>2.36</td><td>2.08</td><td>2.15</td><td>2.39</td></tr><tr><td>S4</td><td>2.93</td><td>2.71</td><td>2.86</td><td>2.49</td><td>2.20</td><td>2.47</td><td>2.43</td><td>2.25</td><td>2.46</td></tr><tr><td rowspan="4">5</td><td>S1</td><td>2.30</td><td>1.89</td><td>2.11</td><td>2.03</td><td>1.65</td><td>1.94</td><td>1.79</td><td>1.45</td><td>1.73</td></tr><tr><td>S2</td><td>2.54</td><td>2.27</td><td>2.49</td><td>2.25</td><td>2.07</td><td>2.28</td><td>2.06</td><td>1.90</td><td>2.13</td></tr><tr><td>S3</td><td>2.71</td><td>2.59</td><td>2.74</td><td>2.45</td><td>2.39</td><td>2.59</td><td>2.30</td><td>2.27</td><td>2.47</td></tr><tr><td>S4</td><td>2.95</td><td>2.66</td><td>2.80</td><td>2.74</td><td>2.47</td><td>2.65</td><td>2.62</td><td>2.35</td><td>2.55</td></tr><tr><td rowspan="4">6</td><td>S1</td><td>2.42</td><td>2.00</td><td>2.22</td><td>2.11</td><td>1.55</td><td>1.82</td><td>1.98</td><td>1.43</td><td>1.78</td></tr><tr><td>S2</td><td>2.61</td><td>2.33</td><td>2.45</td><td>2.36</td><td>2.02</td><td>2.26</td><td>2.26</td><td>1.94</td><td>2.19</td></tr><tr><td>S3</td><td>2.76</td><td>2.53</td><td>2.64</td><td>2.57</td><td>2.38</td><td>2.54</td><td>2.46</td><td>2.32</td><td>2.52</td></tr><tr><td>S4</td><td>2.99</td><td>2.57</td><td>2.68</td><td>2.87</td><td>2.43</td><td>2.59</td><td>2.77</td><td>2.39</td><td>2.57</td></tr></tbody></table>
## References

1. Abdesslam, B. Evolution of rural-urban health gaps in Morocco: 1992–2011. *BMC Res. Notes* **2012**, *5*, 381. [CrossRef]

2. Institut National de la Statistique. *Recensement Général de la Population et de l'Habitat 2014 (RGPH 2014)*; Haut Commissariat au Plan (HCP): Casablanca, Morocco, 2014.

3. Cherkaoui, T.E.; Hatzfeld, D.; Jebli, H.; Medina, F.; Caillot, V. Etude microsismique de la région d'Al Hoceima. *Bull. Inct. Sci. Rabat.* **1990**, *12*, 25–34.

4. Kariche, J.; Meghraoui, M.; Timoulali, Y.; Cetin, E.; Toussaint, R. The Al Hoceima earthquake sequence of 1994, 2004 and 2016: Stress transfer and poroelasticity in the Rif and Alboran Sea region. *Geophys. J. Int.* **2018**, *212*, 42–53. [CrossRef]

5. Bezzeghoud, M.; Buforn, E. Source parameters of the 1992 Melilla (Spain, MW = 4.8), 1994 Alhoceima (Morocco, MW = 5.8), and 1994 Mascara (Algeria, MW = 5.7) earthquakes and seismotectonic implications. *Bull. Seismol. Soc. Am.* **1999**, *89*, 359–372. [CrossRef]

6. El Alami, S.O.; Tadili, B.A.; Cherkaoui, T.E.; Medina, F.; Ramdani, M.; Brahim, L.A.; Harnafi, M. The Al Hoceima earthquake of May 26, 1994 and its aftershocks: A seismotectonic study. *Ann. Geophys.* **1998**, *41*, 519–537. [CrossRef]

7. Goula, X.; Gonzalez, M. Visite Technique à la zone Endommagée par le séisme d'Al-hoceima du 24 février 2004. *Generalitat de Catalunya. Institut Cartogràfic de Catalunya*; ICC-GS-192/2004 (fr): Barcelona, Spain, 2004.

8. Jabour, N.; Kasmi, M.; Menzhi, M.; Birouk, A.; Hni, L.; Hahou, Y.; Timoulali, Y.; Bedrane, S. *The February 24th, 2004 Al Hoceima Earthquake*; Newsletter of the European-Mediterranean Seismological Centre: Essonne, France, 2004.

9. Cherkaoui, T.E.; Harnafi, M. Le séisme d'Al Hoceima du 24 février 2004. Trav. Col. Nat. Sur: Région de Taza-Al Hoceima-Taounate, ressources et stratégies de développement. *Pub. Fac. Polydiscip. Taza* **2005**, 27–38. Available online: http://www.afps-seisme.org/content/download/34357/740092/file/AFPS_Rapport_mission_2004_Seisme%20Al%20Hoceima%20(Maroc)%202.pdf (accessed on 8 May 2022).

10. Buforn, E.; Pro, C.; Sanz de Galdeano, C.; Cantavella, J.V.; Cesca, S.; Caldeira, B.; Udías, A.; Mattesini, M. The 2016 south Alboran earthquake (Mw=6.4): A reactivation of the Ibero-Maghrebian region? *Tectonophysics* **2017**, 712–713, 704–715. [CrossRef]

11. Medina, F.; Cherkaoui, T. The South-Western Alboran Earthquake Sequence of January-March 2016 and Its Associated Coulomb Stress Changes. *Open J. Earthq. Res.* **2016**, *6*, 35–54. [CrossRef]

12. Galindo-Zaldivar, J.; Ercilla, G.; Estrada, F.; Catalán, M.; d'Acremont, E.; Azzouz, O.; Casas, D.; Chourak, M.; Vazquez, J.T.; Chalouan, A.; et al. Imaging the growth of recent faults: The case of 2016–2017 seismic sequence sea bottom deformation in the Alboran Sea (Western Mediterranean). *Tectonics* **2018**, *37*, 2513–2530. [CrossRef]

13. Vlachakis, G.; Vlachaki, E.; Lourenço, P.B. Learning from failure: Damage and failure of masonry structures, after the 2017 Lesvos earthquake (Greece). *Eng. Fail. Anal.* **2020**, *117*, 104803. [CrossRef]

14. Bilgin, H.; Shkodrani, N.; Hysenlliu, M.; Ozmen, H.B.; Isik, E.; Harirchian, E. Damage and performance evaluation of masonry buildings constructed in 1970s during the 2019 Albania earthquakes. *Eng. Fail. Anal.* **2022**, *131*, 105824. [CrossRef]

15. Stepinac, M.; Lourenço, P.B.; Atalić, J.; Kišiček, T.; Uroš, M.; Baniček, M.; Novak, M.Š. Damage classification of residential buildings in historical downtown after the ML5. 5 earthquake in Zagreb, Croatia in 2020. *Int. J. Disaster Risk Reduct.* **2021**, *56*, 102140. [CrossRef]

16. RPS. *Règlement de Construction Parasismique. Ministère de l'Aménagement du Territoire, de l'Urbanisme, de l'Habitat et de l'Environnement*; Secrétariat d'Etat à l'Habitat: Rabat, Morocco, 2000.

17. RPS. *Règlement de Construction Parasismique*; Version 2011; Ministère de l'Habitat, de l'Urbanisme et de l'Aménagement de l'Espace: Rabat, Morocco, 2000.

18. Ferreira, T.M.; Mendes, N.; Silva, R. Multiscale Seismic Vulnerability Assessment and Retrofit of Existing Masonry Buildings. *Buildings* **2019**, *9*, 91. [CrossRef]

19. Shabani, A.; Kioumarsi, M.; Zucconi, M. State of the art of simplified analytical methods for seismic vulnerability assessment of unreinforced masonry buildings. *Eng. Struct.* **2021**, *239*, 112280. [CrossRef]

20. Benedetti, D.; Petrini, V. On the seismic vulnerability of masonry buildings: Proposal of year procedure evaluation. *Ind. Constr.* **1984**, *18*, 66–74.

21. Cosenza, E.; Manfredi, G.; Polese, M.; Verderame, G.M. A multi-level approach to the capacity assessment of existing RC buildings. *J. Earthq. Eng.* **2005**, *9*, 1–22. [CrossRef]

22. Lagomarsino, S.; Giovinazzi, S. Macroseismic and mechanical models for the vulnerability and damage assessment of current buildings. *Bull. Earthq. Eng.* **2006**, *4*, 415–443. [CrossRef]

23. Rosti, A.; Rota, M.; Penna, A. Damage classification and derivation of damage probability matrices from L'Aquila (2009) post-earthquake survey data. *Bull. Earthq. Eng.* **2018**, *16*, 3687–3720. [CrossRef]

24. Diana, L.; Lestuzzi, P.; Podestá, L.; Luchini, G. Improved Urban Seismic Vulnerability Assessment Using Typological Curves and Accurate Displacement Demand Prediction. *J. Earthq. Eng.* **2019**, *25*, 1709–1731. [CrossRef]

25. Sangiorgio, V.; Uva, G.; Fatiguso, F.; Adam, J.M. A new index to evaluate exposure and potential damage to RC building structures incoastal areas. *Eng. Fail. Anal.* **2019**, *100*, 439–455. [CrossRef]

26. Rosti, A.; Del Gaudio, C.; Di Ludovico, M.; Magenes, G.; Penna, A.; Polese, M.; Prota, A.; Ricci, P.; Rota, M.; Verderame, G.M. Empirical vulnerability curves for italian residential buildings. *Boll. Di. Geofis. Teor. Ed. Appl.* **2020**, *61*, 357–374.

27. Del Gaudio, C.; Ricci, P.; Verderame, G.M.; Manfredi, G. Development and urban-scale application of a simplified method for seismic fragility assessment of RC buildings. *Eng. Struct.* **2015**, *91*, 40–57. [CrossRef]
28. Gaetani d'Aragona, M.; Polese, M.; Prota, A. Stick-IT: A simplified model for rapid estimation of IDR and PFA for existing low-rise symmetric infilled RC building typologies. *Eng. Struct.* **2020**, *223*, 111182. [CrossRef]

29. Ruggieri, S.; Porco, F.; Uva, G.; Vamvatsikos, D. Two frugal options to assess class fragility and seismic safety for low-rise reinforced concrete school buildings in Southern Italy. *Bull. Earthq. Eng.* **2021**, *19*, 1415–1439. [CrossRef]

30. Leggieri, V.; Ruggieri, S.; Zagari, G.; Uva, G. Appraising seismic vulnerability of masonry aggregates through an automated mechanical-typological approach. *Autom. Constr.* **2021**, *132*, 103972. [CrossRef]

31. Kappos, A.J.; Panagopoulos, G.; Panagiotopoulos, C.; Penelis, C.G. A hybrid method for the vulnerability assessment of R/C and URM buildings. *Bull. Earthq. Eng.* **2006**, *4*, 391–413. [CrossRef]

32. Barbat, A.H.; Moya, F.Y.; Canas, J.A. Damage scenarios simulation for seismic risk assessment in urban zones. *Earthq. Spectra* **1996**, *12*, 371–394. [CrossRef]

33. Pagani, M.; Monelli, D.; Weatherill, G.; Danciu, L.; Crowley, H.; Silva, V.; Henshaw, P.; Butler, L.; Nastasi, M.; Panzeri, L.; et al. OpenQuake Engine: An Open Hazard (and Risk) Software for the Global Earthquake Model. *Seismol. Res. Lett.* **2014**, *85*, 692–702. [CrossRef]

34. Crowley, H.; Despotaki, V.; Rodrigues, D.; Silva, V.; Toma-Danila, D.; Riga, E.; Karatzetzou, A.; Fotopoulou, S.; Zugic, Z.; Sousa, L.; et al. Exposure model for European seismic risk assessment. *Earthq. Spectra* **2020**, *36* (Suppl. S1), 252–273. [CrossRef]

35. Silva, V.; Amo-Oduro, D.; Calderon, A.; Costa, C.; Dabbeek, J.; Despotaki, V.; Martins, L.; Pagani, M.; Rao, A.; Simionato, M.; et al. Development of a global seismic risk model. *Earthq. Spectra* **2020**, *36* (Suppl. S1), 372–394. [CrossRef]

36. Louhibi, S.; Chourak, M.; El Kihel, E.; El Hammoumi, A.; Douiri, A. Evaluation of the seismic vulnerability of buildings in Al Hoceima (Morocco). *Int. J. Earth Sci. Eng.* **2013**, *6*, 1505–1511.

37. Cherif, S.; Chourak, M.; Abed, M.; Pujades, L. Seismic risk in the city of Al Hoceima (north of Morocco) using the vulnerability index method, applied in Risk-UE project. *Nat. Hazards* **2017**, *85*, 329–347. [CrossRef]

38. Cherif, S.; Chourak, M.; Abed, M.; Douiri, A. Potential Seismic Damage Assessment of Residential Buildings in Imzouren City (Northern Morocco). *Buildings* **2018**, *8*, 179. [CrossRef]

39. Cherif, S.; Chourak, M.; Abed, M.; Pujades, L. Seismic evaluation method for existing reinforced concrete buildings in North of Morocco. *Bull. Earthquake Eng.* **2019**, *17*, 3873–3894. [CrossRef]

40. Milutinovic, Z.V.; Trendafiloski, G.S. *Risk-UE an Advanced Approach to Earthquake Risk Scenarios with Applications to Different European Towns*; Contract: EVK4-CT-2000-00014, WP4; Vulnerability of Current Buildings; 2003; pp. 1–111. Available online: http://www.civil.ist.utl.pt/~{}mlopes/conteudos/DamageStates/Risk%20UE%20WP04_Vulnerability.pdf (accessed on 8 May 2022).

41. Freeman, S.A. Review of the development of the capacity spectrum method. *ISET J. Earthq. Technol.* **2004**, *41*, 1–13.

42. Freeman, S.A. Evaluations of existing buildings for seismic risk-A case study of Puget Sound Naval Shipyard. In Proceedings of the U.S. National Conference on Earthquake Engineering, Bremerton, WA, USA, 18–20 June 1975; pp. 113–122.

43. Mahaney, J.A.; Paret, T.F.; Kehoe, B.E.; Freeman, S.A. *The Capacity Spectrum Method for Evaluating Structural Response during the Loma Prieta Earthquake*; Mitigation and damage to the built environment; U.S. Central United States Earthquake Consortium (CUSEC); Memphis, Tennessee, 1993; pp. 501–510.

44. HAZUS. *Earthquake Loss Estimation Methodology Technical Manual*; National Institute of Building Sciences for Federal Emergency Management Agency (FEMA); Washington, DC, USA, 1999.

45. Grünthal, G. *European Macroseismic Scale 1998*; European Seismological Commission (ESC); Luxembourg, 1998.

46. Seminar, C. *NEHRP Guidelines for the Seismic Rehabilitation of Buildings (FEMA 273)*; Building Seismic Safety Council: Washington, DC, USA, 1997.

47. van der Woerd, J.; Dorbath, C.; Ousadou, F.; Dorbath, L.; Delouis, B.; Jacques, E.; Tapponnier, P.; Hahou, Y.; Menzhi, M.; Frogneux, M.; et al. The Al Hoceima Mw 6.4 earthquake of 24 February 2004 and its aftershocks sequence. *J. Geodyn.* **2014**, *77*, 89–109. [CrossRef]

48. Cherkaoui, T.E.; El Hassani, A. Seismicity and Seismic Hazard in Morocco 1901–2010. *Bull. L'institut Sci. Sect. Sci. Terre* **2012**, *34*, 45–55.

49. Mourabit, T.; Abou Elenean, K.M.; Ayadi, A.; Benouar, D.; Ben Suleman, A.; Bezzeghoud, M.; Cheddadi, A.; Chourak, M.; ElGabry, M.N.; Harbi, A.; et al. Neo-deterministic seismic hazard assessment in North Africa. *J. Seismol.* **2013**, *18*, 301–318. [CrossRef]

50. Valerio, P.; Garcia-Peláez, J.; Styron, R.; Pagani, M.; Gee, R. A probabilistic seismic hazard model for North Africa. *Bull. Earthq. Eng.* **2020**, *18*, 2917–2951. [CrossRef]

51. d'Acremont, E.; Gutscher, M.A.; Rabaute, A.; Mercier de Lépinay, B.; Lafosse, M.; Poort, J.; Ammar, A.; Tahayt, A.; Le Roy, P.; Smit, J.; et al. High-resolution imagery of active faulting offshore Al Hoceima, Northern Morocco. *Tectonophysics* **2014**, *632*, 160–166. [CrossRef]

52. Tahayt, A.; Feigl, K.L.; Mourabit, T.; Rigo, A.; Reilinger, R.; McClusky, S.; Fadil, A.; Berthier, E.; Dorbath, L.; Serroukh, M.; et al. The Al Hoceima (Morocco) earthquake of 24 February 2004, analysis and interpretation of data from ENVISAT ASAR and SPOT5 validated by ground-based observations. *Remote Sens. Environ.* **2009**, *113*, 306–316. [CrossRef]

53. European Committee for Standardization. *Eurocode 8: Design Provisions for Earthquake Resistance of Structures*; Central Secretariat, CEN: Bruxelles, Belgium, 1994.

54. BAEL 91. *Règles Techniques de Conception et de Calcul des Ouvrages et Constructions en Béton Armé suivant la Méthode des Etats Limites*; Edition Eyrolles: Paris, France, 1992; ISBN 221210023X.
55. RMSI. Morocco natural hazards probabilistic analysis and national strategy development. In *Earthquake Hazard Report*; Department of Economic and General Affairs: Rabat, Morocco, 2012.

56. Douiri, A.; Mourabit, T.; Cheddadi, A.; Chourak, M.; Louhibi, S. Coïncidences spectrales du bruit de fond entre sol et bâti existant dans la ville d'Imzouren, Maroc (Spectral coincidences of background noise between soil and existing buildings in Imzouren city, Morocco). *J. Mater. Environ. Sci.* **2015**, *6*, 366-376.

57. Council, A.T. *Seismic Evaluation and Retrofit of Concrete Buildings*; Report No. SSC 96-01: ATC-40, 1; Applied Technology Council: Redwood City, CA, USA, 1996.