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

: 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 signiﬁcantly. 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, conﬁrming that the damage is more signiﬁcant 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.


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.
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][8][9] and, more recently, the 2016 earthquake [10][11][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][21][22][23][24][25][26]; (ii) mechanical and analytical methods, which are more detailed but laborious for a largescale study.They use slightly more detailed algorithms for vulnerability assessments, which allow more delicate and detailed studies [27][28][29][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][37][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][38][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].

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: 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: where Sd ds i is the threshold spectral displacement at which the probability of the damage state ds i is 50%, β ds i is the standard deviation of the natural logarithm of this spectral displacement, φ 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 Sd ds i are obtained as a function of the yield displacement D y and the ultimate displacement D u of the structure.Concerning β ds , it is well known that the expected seismic damage in buildings follows a binomial probability distribution [45].β ds is calculated directly as a function of ultimate ductility µ u [40] and is displayed in Table 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][11][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][49][50].This issue has been the subject of multiple studies related to seismology, seismic hazard, and seismic vulnerability [37][38][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.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.

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   () 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.

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 e (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 2. Seismic acceleration values based on the national seismic standard [17].3. Table 3. Ground types, according to Eurocode 8 [53] and the national seismic standard (RPS2000) [17].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   results: Table 3. Ground types, according to Eurocode 8 [53] and the national seismic standard (RPS2000) [17].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.


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.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.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.


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.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.

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 2 and 150 m 2 [39].
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].

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 2 and 150 m 2 [39].

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 2 and 150 m 2 [39].
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 2 );  Most buildings are regular structures having simple geometrical shapes.
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.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;  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 2 );  Most buildings are regular structures having simple geometrical shapes.
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.

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.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.

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].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 2 and 150 m 2 .Taking a surface of 120 m 2 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].Table 5. Cross-sectional area of columns and adopted steel rebar for buildings with seismic code [39].

Cross-Sectional Area (cm 2 )
Floor 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 Table 5. Cross-sectional area of columns and adopted steel rebar for buildings with seismic code [39].Table 6.Cross-sectional area of columns and adopted steel rebar for buildings without seismic code [39].

Cross-Sectional Area (cm
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).
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).

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  ̅̅̅̅  and   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.

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 Sd i and β 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 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).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.

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][38][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][38][39] for these two cities.

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 vs 30.

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

Figure 1 .
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.

Figure 1 .
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.

Figure 2 .
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].

Figure 2 .
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].

Figure 3 .
Figure 3.The 5% elastic response spectra for the probabilistic scenario expected on outcrop rock.

Figure 4 .Figure 3 .
Figure 4. Distribution of soil types in the city of Al Hoceima [37].

Figure 3 .
Figure 3.The 5% elastic response spectra for the probabilistic scenario expected on outcrop rock.

Figure 4 .
Figure 4. Distribution of soil types in the city of Al Hoceima [37].

Figure 4 .
Figure 4. Distribution of soil types in the city of Al Hoceima [37].

Figure 7 .
Figure 7. Distribution of building typologies in Al Hoceima and Imzouren of Central Rif [2].

Figure 6 .
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.

Figure 8 .
Figure 8. Examples of residential buildings in the study region.

Figure 8 .
Figure 8. Examples of residential buildings in the study region.

Figure 8 .
Figure 8. Examples of residential buildings in the study region.

Figure 10 .
Figure 10.Inspected buildings in the city of Imzouren based on (a) seismic code and (b) construction period.

Figure 12 .
Figure 12.Capacity spectra of buildings with seismic code built after 2002.

Figure 12 .
Figure 12.Capacity spectra of buildings with seismic code built after 2002.

Figure 14 .
Figure 14.Damage distribution in the city of Al Hoceima following (a) the RPS2000 risk scenario and (b) the Eurocode 8 risk scenario.

Figure 14 .
Figure 14.Damage distribution in the city of Al Hoceima following (a) the RPS2000 risk scenario and (b) the Eurocode 8 risk scenario.

Figure 15 .
Figure 15.Damage distribution in the city of Imzouren following (a) the RPS2000 risk scenario and (b) the Eurocode 8 risk scenario.Figure 15.Damage distribution in the city of Imzouren following (a) the RPS2000 risk scenario and (b) the Eurocode 8 risk scenario.

Figure 15 .
Figure 15.Damage distribution in the city of Imzouren following (a) the RPS2000 risk scenario and (b) the Eurocode 8 risk scenario.Figure 15.Damage distribution in the city of Imzouren following (a) the RPS2000 risk scenario and (b) the Eurocode 8 risk scenario.
y + 0.25 D u − D y β S d3 = 0.1 + 0.4 ln(µ u ) d4 = 0.15 + 0.5 ln(µ u ) Complete * 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 Table3.

Table 2 .
[17]mic acceleration values based on the national seismic standard[17]..  *   is the acceleration value and   is the velocity value corresponding to the seismic zone.**  is the soil factor, having values displayed in Table

Table 6 .
[39]s-sectional area of columns and adopted steel rebar for buildings without seismic code[39].

Table A2 .
Yield and ultimate capacity for buildings without seismic code.

Table A3 .
Parameters characterizing the fragility curves for the building models.

Table A4 .
Spectral displacement (cm) of the performance points for the RC building models according to both probabilistic scenarios (RPS2000 and Eurocode 8).

Table A5 .
Damage index DS m values for the RC building models according to both probabilistic scenarios (RPS2000 and Eurocode 8).