Site Response Evaluation in the Trans-Mexican Volcanic Belt Based on HVSR from Ambient Noise and Regional Seismicity

Featured Application: The results obtained in this study allow an initial overview of the variation of the site response in an area for which low to moderate seismic risk is usually considered. The information presented may be used to analyze the seismic effects in the study zone associated with the geological characteristics. Abstract: The Trans-Mexican Volcanic Belt (TMVB), located in central Mexico, is an area for which low to moderate seismic risk is considered. This is based on the limited instrumental data available, even though large historical earthquakes have damaged some urban centers in the past. However, site effects is an aspect that must be considered in estimating risk, because there are some instances of important ampliﬁcations that have been documented with serious effects. In this work, ambient noise and earthquake records from 90 seismic permanent and temporary stations are used to analyze site response in the TMVB. The results obtained show a heterogeneous range in the value of the fundamental frequency. When possible, a comparison was made of the results obtained from ambient noise and earthquake records. In almost all these comparisons, no signiﬁcant differences were observed in terms of the fundamental frequency. However, there were some stations with a ﬂat average HVSR ambient noise curve that contradicted earthquake data results, which showed peaks at some frequencies. Our results are a ﬁrst step towards categorizing the different site responses in the TMVB but in order to provide ﬁner details, it is necessary to improve the actual monitoring conditions.


Introduction
Within the context of seismic engineering, the study of site response involves changes in variables related to the seismic intensity, in terms of amplitude, duration, and frequency content. These changes depend on the geological features at the measurement site, and usually, lead to larger amplitudes on soil sites than on hard rock.
Throughout history it is possible to find documented cases in which site effects have played a decisive role in observed damage after significant earthquakes, such as in San Francisco (1906), Mexico City (1985 and 2017), and Kobe (1995). Thus, in the seismic design of buildings and civil infrastructure it is important to estimate the maximum expected intensities at the site.
The evaluation of changes in amplitude of ground motion of a particular site is commonly made with respect to a reference position. The most popular and reliable way to do this is through the Standard Spectral Ratio (SSR) technique [1]. This method is based on the ratio of the Fourier amplitude spectra observed during an earthquake at soil A i ( f ) to reference site A j ( f ). The Fourier spectra of ground motion A( f ) can be expressed as the multiplication of source S( f ), path P( f ), and site H( f ) terms in the frequency domain: where f is the frequency.
In the SSR technique, the reference site is required to have the most homogeneous geological conditions as it is assumed to be free of anomalous amplifications. The two sites should be close to each other in comparison with the distance to the source. Considering that the analysis is performed for stations which register the same earthquake: However, it is not always possible to find records of the same event with a good signal-to-noise (S/R) ratio in two stations that meet these conditions. Thus, the Horizontalto-Vertical Spectral Ratio (HVSR) technique [2] is used as an alternative to studying the site response. In this method, the ratio between the mean of the horizontal components A Hi ( f ) and the vertical component A Vi ( f ), recorded in a single station, is calculated: This technique can be used to estimate the fundamental resonant frequency ( f 0 ) of a site, but is unreliable for determining its transfer function [3].
In the literature, several studies can be found in which the site response is analyzed in areas around the world, based on these empirical techniques and numerical models [4][5][6][7][8].
The study of this phenomenon remains current and it is an important issue to be considered in the estimation of seismic risk and damage prevention, even in areas of low seismicity [9].
The Trans-Mexican Volcanic Belt (TMVB), located in central Mexico (Figure 1), is an example of those areas where the site response might represent a potential risk despite the low frequency of earthquakes. It is a volcanic arc that covers Cretaceous and Cenozoic magmatic provinces [10]. Seismicity in this region is related to extensional faults in the crust and normal faults in the subducted Cocos plate. This zone is located between 18°30 and 21°30 , and extends from the coast of the Pacific Ocean to the Gulf of Mexico. It has an approximate length of 1000 km and a variable width of 90-230 km. This area presents low seismicity, compared to the Pacific Ocean coast, where the Cocos plate subducts under the North American plate. Nevertheless, there are important recent and historical occurrences of destructive events [11,12]. Examples of these is the earthquake that occurred in Acambay on 12 November 1912 (Mw 6.9, H = 33 km) [13], and more recently in Puebla on 19 September 2017 (Mw 7.1, H = 51.2 km) [14]. This aspect is particularly important, in view that the TMVB comprises some of the largest cities in Mexico with significant industrial development and population growth in the last years.
It may be difficult to find a zone in the TMVB with negligible site effects, due to the variability in geological properties [10]. Several articles have been published on this topic related to Mexico, mainly focused on the behavior of site effects in Mexico City.
Ordaz et al. [15] found that stations located in the hill zone in the Valley of Mexico present amplifications that are 10 times higher than those predicted from ground motion attenuation equations. In a later study, García et al. [16] compared amplification responses at some of the stations analyzed by [15]. The records analyzed corresponded to inslab earthquakes with an epicenter outside the TMVB, and they observed similar behavior despite the source type and location.
Singh et al. [17] briefly discuss HVSR curves for 2 stations crossing the TMVB, obtained from records of 9 shallow coastal earthquakes on the Pacific Ocean coast. Likewise, Lozano et al. [18] studied the influence of source characteristics on the site effects in the Valley of Mexico from 36 interplate and inslab Mexican earthquakes and 12 teleseismic events originated in South America. Their results agree with those obtained by [16] and concluded that observed site effects were independent of the characteristics and location of the source.
Clemente-Chávez et al. [19] presented the first study of site effects in the TMVB outside the Valley of Mexico, based on local shallow seismicity. They analyzed HVSR curves obtained from 22 earthquakes recorded by 25 seismological stations located in the study zone. The events studied have depths H < 10 km, and magnitudes between 3.6 and 4.3. Average values for the fundamental frequency f 0 and amplification factors obtained from the HVSR curves were reported. Their results were compared with previous studies and they attributed the differences found to the location of the source. However, due to the limited amount of data employed, some of their conclusions were supported only with one or two records. Furthermore, it should be considered that HVSR curves do not unequivocally represent the actual site amplification.
Considering the low density of permanent seismological stations in the TMVB [20], the number of records that can be analyzed with the SSR technique is very limited. So, the use of ambient noise records represents a feasible option to include additional data for the analysis of site response in terms of f 0 . In this paper, we analyze HVSR curves obtained from ambient noise and regional instrumental records related to 121 crustal and inslab events with an epicenter inside the TMVB. Data considered were recorded by 90 seismological stations belonging to seven seismic networks.
The results of our study may be useful in disaster prevention and in estimating the behavior of future buildings in the study area. The information presented may also be used to complement studies examining the damage observed to existing infrastructure after the occurrence of a seismic event [21,22]. However, it is important to consider that a more accurate and reliable analysis can be carried out throughout microzonation studies, so our results only give a general overview.

Data and Methods
To perform this study, we compiled a catalog of earthquakes with epicenter inside the area of the TMVB, based on the seismicity catalogs of the Servicio Sismológico Nacional [23], the U.S. Geological Survey [24], and the International Seismological Centre (ISC) [25]. Avoiding duplications and magnitudes reported as "non-calculable", 2113 earthquakes with magnitude 1.5-7.8 and depth 1-191 km, were selected.
Based on this information, instrumental records were searched in seven permanent and temporary seismological networks in Mexico (Table 1), with broadband stations located inside the TMVB. It should be considered that over time there have been significant changes in instrumental density in this zone [20]. Additionally, the stations that belong to permanent networks are widely dispersed, and the availability of continuous data is limited to the last 30 years. Data were gathered from stations with velocity sensors belonging to Servicio Sismológico Nacional (SSN) [23], Centro de Geociencias de la Universidad Nacional Autónoma de México (CGEO), GEOSCOPE [26], and the temporary networks: The MesoAmerican Subduction Experiment (MASE) [27,28], Colima Deep Seismic Experiment (CODEX) [29], Mapping the Rivera Subduction Zone (MARS) [30], and Geometry of the Cocos Plate (GECO) [31].
The stations belonging to SSN and CGEO are equipped with STS-2 broadband seismometers with sampling rates of 80 and 100 Hz. Only the ACIG station has a Trillium 120 seismometer with a sampling rate of 20 Hz. In the case of the GEOSCOPE network, the UNM station has a velocity sensor STS-1 20 samples/s. During their period of operation, the GECO network stations were equipped with Reftek 151B-60 and Guralp 40T sensors. The temporary networks MASE and CODEX had CMG 3T, and 40T seismometers, respectively, with a sampling rate of 100 Hz. The MARS network was equipped with CMG 3T sensors of 40 sps.
Only records that met the following criteria were considered for the analyses: (1) Station and epicenter located inside the TMVB; (2) signal-to-noise ratio (S/R) > 2.0; and (3) clear arrival of P and S waves, based on visual inspection.
Considering these restrictions, a total of 352 records related to 24 inslab and 97 crustal earthquakes were chosen ( Figure 2). Such events have magnitudes between 2.0 and 7.1, a hypocentral depth of 1 ≤ H ≤ 105 km, and occurred between 1993 and 2019 in the TMVB. The main characteristics of the selected earthquakes are listed in Table 2.
In Figure 2 it can be observed that in the central zone there is the lowest number of records. As mentioned by Zúñiga et al. [20], the other regions present higher seismic activity due to active fault systems, the subduction zone on the Pacific Ocean coast, and the occurrence of inslab earthquakes in the central-eastern sector.
The number of records selected for each station was variable, and in most of them, the data with the characteristics sought was considerably limited. An HVSR analysis was performed at all stations with at least 10 earthquake records. This threshold was considered to include as many stations as possible, trying not to compromise the reliability of the results. It is important to note that despite the large number of seismological stations belonging to MASE, MARS, and CODEX networks, the suitable records were limited due to the characteristics of the earthquakes during their short period of operation.
The processing of records was performed by means of the package Geopsy [32]. For each analyzed record, the mean and the trend were removed. The corresponding Fourier Acceleration Spectra (FAS) were smoothed using the Konno-Omachi function [33], considering a b-value of 20. Window lengths that included 95% of the total energy were taken, beginning from the S-wave arrival. Due to the lack of earthquake records, an HVSR analysis was performed in all the stations using ambient noise data. Although the energy of ambient noise does not compare with that of an earthquake, this second analysis allowed us to increase the number of stations analyzed and to compare the results obtained with both procedures. For each station, a database of 30 records of one-hour duration was collected on random dates within its period of operation. In all ambient noise records, the mean and trend were removed and the same smoothing process was applied as in the earthquake records. For the HVSR analyses, the records were divided into 60 windows of 1 min and the root mean square of the horizontal components was calculated. SESAME criteria [34] were considered for the analyses and classification of the curves obtained. Finally, the results were correlated with the geological characteristics of the study area.

Results
Due to the large extension of the study area, the stations analyzed were grouped into four zones. This division was based on the distribution of stations and on geologic information given by [10]. Most of the average curves shown correspond to the HVSR analyses based on ambient noise data, between 0.1 and 25 Hz. Some of the considered stations have a Nyquist frequency of 10 and 20 Hz, so not all curves are displayed in the same frequency range. When possible, a comparison is made with results from earthquake records. Regarding geology, information from the Geologic Map of North America [35] was used.

Western TMVB
In the westernmost part of the TMVB, there are three main fault systems: Tepic-Zacoalco, Colima, and Chapala [10]. Regarding significant earthquakes, there is a historical account of the occurrence of a damaging event during the 16th century near the town of Ameca, in the state of Jalisco. For this event, a magnitude Mw 7.2 ± 0.3 was estimated from a rupture vs. magnitude scaling relation [36,37]. Another earthquake occurred on 2 October 1847, with a mean magnitude of 5.7 ± 0.4 which has been associated to the faults of the Chapala Graben. It affected several towns and produced important damages to buildings. An intensity of IX has been estimated for this event [38].
Later, on 11 February 1875 an earthquake near the city of Guadalajara caused great destruction in nearby areas including dozens of deaths [39].   In this sector, there are a variety of geological formations ranging from the Lower Cretaceous to the Quaternary [35]. Most of the stations analyzed are located in the vicinity of the Colima volcanic complex. Figure 4 shows some of the average HVSR curves obtained in this zone, in which only the ANIG station could be analyzed with earthquake records.

Mafic rocks (Neogene and Quaternary)
In some curves there can be observed flat responses or complex shapes with low amplitudes at different frequencies: MA47, MORA, MAZE, CUAT, ZAPO, MA24, OLOT, SNID, SINN, EMBG ( Figure 4), SCRI, and COMA ( Figure 5). Station MA47 is located on sedimentary rocks from the Lower Cretaceous. Its HVSR curve has a flat shape and no significant impedance contrast is identified. In contrast, stations MORA, MAZE, CUAT, and ZAPO were deployed on Oligocene felsic rocks and near their location there are volcanic and sedimentary rocks from the Neogene and Quaternary. The site response could be attributed to lateral heterogeneity and the presence of sedimentary rocks from the Cretaceous. Stations MA24, OLOT, SINN, and EBMG were located on mafic and intermediate rocks from the Neogene and Quaternary. MA24 and OLOT were relatively close to each other, but a resonant frequency can not be identified in either of them. As for SINN and EBMG, the corresponding curves have a flat shape, indicating hard rock sites. Stations SCRI and COMA ( Figure 5) were deployed on Quaternary sedimentary rocks, and in both of them a predominantly flat shape is also observed.
There is another group of stations in which the HVSR curves present multiple peaks, with amplitudes similar to what could be identified as the clearest peak. Station MA54 were located on sedimentary rocks from the Upper Cretaceous. The peak around 6 Hz could be attributed to a thin layer of sediments and the amplitude between 0.6 and 2 Hz to the contact with formations from the Lower Cretaceous. Meanwhile, MA46 were deployed on Oligocene felsic rocks. The observed peaks above 2 Hz could indicate multiple layers of sediments of varying thickness above the bedrock. Station MA42 was on mafic and intermediate rocks from the Neogene and Quaternary. The shape of the HVSR curve corresponding to this station could be attributed to variations in the topography or a deep interface with Oligocene formations. SNID ( Figure 4) and MA48 ( Figure 5) were deployed on Quaternary sedimentary rocks. The amplitudes in SNID could be associated to lateral heterogeneities by means of its location near Neogene formations. Likewise, some of the multiple peaks in MA48 could be spurious or related to multiple layers of sediments at variable depth.  Figure 5). This is also indicative of possible amplifications in ground motion during the occurrence of seismic events. In the case of stations BAVA and JANU, located on Oligocene felsic rocks, the observed peaks at low frequencies could be associated to thick sedimentary deposits. High clear peaks at 3 and 4 Hz are observed at stations MA44 and MA45. These may be produced by strong impedance contrasts with interlayered sedimentary and volcanic rocks from the Lower Cretaceous.  (Figure 4). In general, a clear peak is observed above 5 Hz, which could be attributed to impedance contrast with thin layers of sediment. It should be noted that the average HVSR curve obtained with earthquake data at ANIG station does not vary significantly in its shape to that obtained from ambient noise. There is a peak above 10 Hz that could have been produced by the presence of thin volcanic layers.
Regarding the average HVSR curves of stations on Quaternary sedimentary rocks (ESPN, CDGZ, COLM, and CANO) there are varied shapes ( Figure 5). In ESPN, there is a peak above 10 Hz which may require the analysis of a wider frequency bandwidth. The clear peak in CANO could be associated to contact with Neogene and Quaternary rocks. As for CDGZ, there is a high peak near 1 Hz, attributable to impedance contrast with layers of the Lower Cretaceous or thick layers of sediments. Likewise, in COLM there is a small peak around 10 Hz, probably due to the presence of shallow unconsolidated sediments.

Central TMVB
The Central sector is the one with the lowest instrumental density. There is a group of nine temporary stations concentrated on the southwest of the sector, and only two permanent stations that are considerably separated in the east ( Figure 6). There were few earthquake records that met the established restrictions. For this reason, the results presented here only correspond to the analysis of ambient noise data.
Most of the stations in this sector were deployed on Neogene and Quaternary mafic rocks. In MA15, MA18, MA27, MA29, and MOIG, there are HVSR curves with flat shapes and no significant impedance contrast is observed, while stations MA20 and MA28 do not have clear peaks (Figure 7). Station MA21 has a clear peak at 1 Hz and additional unclear peaks above this frequency. In this case, the shape of the HVSR curve could be related to the calderas that predominate in the area or lateral variations with sedimentary rocks from the Lower Cretaceous. Furthermore, the curve corresponding to the IGIG station shows peaks at frequencies beyond 10 Hz. This permanent station is located on Felsic rocks from Neogene, and its response could be related to thin layers of sediments.   In stations MA16 and MA17, there are clear peaks at 1.5 and 5 Hz, respectively. This could be attributed to impedance contrasts with layers of sediments of different thickness.

Central-Eastern TMVB
In the central-eastern zone of the TMVB is the Acambay fault system, covering some faults considered as active [40]. In the same area, there are the Taxco-San Miguel de Allende and Chapala-Tula systems, with evidence of activity in recent years. In this sector, a large number of relevant earthquakes can be mentioned over time, highlighting the following.
In the 18th century, there are reports of crustal earthquakes related to the Venta de Bravo fault. The first occurred in November 1734, followed by more than 30 strong and small events between November 1734 and March 1735 [39]. The largest crustal earthquake in the TMVB during the 20th century took place on 12 November 1912, in Acambay [13]. It had a magnitude of Mw 6.9 with a maximum intensity of IX in the epicentral area. Then, a series of 90 earthquakes were registered between February and June 1979 in a region comprising Maravatio and Mexico state. The mainshock occurred on 22 February 1979 with a magnitude Mw 5.3 and maximum intensity of VIII [40].  In this sector, there are some stations located on Neogene and Quaternary geological units (Figure 8). In JRQG, ECID, SNLU, COAC, PTRP, ESTA, and MIXC there are mostly flat HVSR curves, and no clear peaks are identified. In DHIG, CUIG, UNM, SABI, and PSIQ, the shape of the average curves obtained from ambient noise and earthquake records do not vary significantly ( Figure 9). Small peaks can be observed at different frequencies that show no significant impedance contrasts or possible amplifications. As a result of this, the locations of these stations and those mentioned above may be considered as hard rock sites. In the HVSR curve corresponding to station TEPE (Figure 10), it is not possible to point out some value as the resonance frequency. This shape could be produced by lateral heterogeneity due to the vicinity of some Quaternary mafic rocks.
In some stations located on Quaternary volcanic rocks, there are almost HVSR curves of a flat shape (ACIG, PACH, MIMO, SAPE, SAPA, and TOSU). On the other hand, stations CUCE and CUNO have HVSR curves in which it is not possible to clearly identify any resonant frequency and could be considered as flat (Figure 10).
In station SALU, deployed on Quaternary sedimentary rocks, there are multiple peaks below 1 Hz and between 3 and 10 Hz (Figure 9). This could be due to the influence of topography or underlying geological units. The curve corresponding to station CHIC, located on Quaternary volcanic rocks, has multiple peaks at high frequencies, but only the one at 20 Hz can be considered as clear, indicating shallow layers of sediments. Additionally, in stations VEGU, ARBO, and VLAD, there are not very well-defined peaks between 1 and 10 Hz (Figure 10). These could be attributed to small impedance contrasts at different depths or lateral heterogeneity due to the presence of Neogene sedimentary rocks.  There are some stations in which clear peaks of varied amplitude can be identified. In the stations TECA, TIZA, and PASU, located on Quaternary sedimentary rocks, there are small peaks of a wide shape but low amplitude ( Figure 9). This could indicate the presence of deposits of different stiffness and width, which by consolidation time have no significant contrasts with the underlying bedrock. The MULU and CIRE stations are special cases since the average HVSR curves present very high peaks at frequencies below 1 Hz. These responses can be explained by the fact that these temporary stations were located in the lakebed zone of Mexico City. In this area, there are reports of devastating effects during the earthquakes of 1985 and 2017.
At some of the remaining stations deployed on Quaternary volcanic rocks, it was possible to perform the analyses with ambient noise and earthquake data (NOGA, AGBE, KM67, and PPIG). The shape of the average HVSR curves obtained did not differ significantly and the amplitude only increased when earthquake records were used ( Figure 10). A particularly special case is that of station KM67, where the HVSR curve obtained using ambient noise data is flat. However, when doing analysis with earthquake data, a clear peak is observed, revealing a significant impedance contrast. At this site, the energy associated with the seismic movement may have produced reverberations in underlying Cretaceous layers. In the case of the PPIG station, located at the Popocatépetl volcano, there is a peak above 10 Hz using ambient noise and earthquake data. Likewise, a peak below 1 Hz can be observed in the analysis with earthquake records. In ATOT, SUPA, and PTCU there are clear peaks between 1 and 3 Hz that may be due to impedance contrasts with Cretaceous layers.

Eastern TMVB
Regarding relevant historic earthquakes in this sector, it can be mentioned the one occurred in 1546 that caused important damage in the city of Jalapa and nearby towns. There is no specific date reported for this event nor wider descriptions in the historical records, but it is known that the first Catholic church built in America was destroyed [41]. Later, on 4 January 1920 an earthquake with magnitude Mw 6.4 occurred in the city of Jalapa. This event caused structural damage in nearby towns and it seconds the 1985 Michoacán earthquake in terms of fatalities. Figure 11 shows the seismological stations analyzed in this sector. Most of them belong to the temporary network GECO and only two to the permanent network of the SSN. In this sector, it was possible to analyze almost half of the stations with earthquake data ( Figure 12). In the case of the LVIG station, no significant differences are found between the curves obtained using ambient noise and the ones obtained through earthquake data. As well as in this station, in AYAH and TPIG, no contrast of impedances is observed, and the shape of the curves is almost flat, indicating hard rock sites.
Moreover, TEPY was located on Quaternary mafic rocks, while QUEC and LUPE on Neogene and Quaternary sedimentary rocks. In these stations, there are multiple peaks. For station TEPY it is not possible to identify a value for the resonant frequency. As for QUEC and LUPE, there is a peak near 10 Hz, which may indicate the presence of thin deposits of sediments. Likewise, the peaks at lower frequencies could be due to the influence of sedimentary rocks from the Neogene and Upper Cretaceous. Stations TATA, NAOL, HUEY, and HUAT are on Quaternary mafic rocks and have varied fundamental frequencies. The observed clear peaks could be attributed to impedance contrasts related to layers of sediments with varied thickness. In TATA and HUAT, there are visible peaks between 1 and 10 Hz. The shape of these curves does not vary significantly when using earthquake and ambient noise data. In the case of the NAOL station, the amplitude of the peak is increased but its location in the considered frequency range does not change.

Conclusions
From the results obtained in this study, it is observed that the site response in terms of the fundamental frequency f 0 had a wide variety throughout the TMVB and there was no visible correlation between the shape of the HVSR curves and extent of any particular geological unit. In almost 46% of the stations analyzed, the response was mostly flat or it was not possible to identify clear peaks. In most of these cases, no significant impedance contrasts were observed, so they could be considered as hard rock sites. There was another small group of stations distributed in the four sectors, in which multiple peaks were observed in the HVSR curves. Some of these peaks could be identified as spurious or attributable to the influence of underlying layers of varying stiffness and thickness.
Additionally, in approximately 36% of the analyzed sites there were clear peaks of varied amplitude in the HVSR curves. In these cases, it was possible to identify the value of the resonance frequency, which was mostly between 1 and 10 Hz. Similarly, in some sites there were peaks in frequencies lower than 1 Hz, which should not be underestimated considering the modern trend of tall buildings in the populated areas of the TMVB.
To date, most of the seismicity studies in Mexico have focused on earthquakes originated on the Pacific Ocean coast. Except for Mexico City, due to its geotechnical properties, the amplitude of these events usually reaches the TMVB considerably attenuated. For this reason, this study focused on analyzing what happens in this zone with regional earthquakes. In this regard, it is important to note that there were sites with high peaks in the four sectors analyzed, not only in the central-eastern sector where Mexico City is located.
In some stations, it was possible to make a comparison of the results obtained from ambient noise and earthquake records. In most of them, there were no significant differences in the shape of the average HVSR curves. In the analyses carried out with seismic records, no evidence of non-linear behavior was identified in the site response with respect to the ambient noise results. This could be attributed to the fact that most of the events considered are of low magnitude. In general, there are few major earthquakes with both epicenter and instrumental records inside the TMVB. However, there were some sites with a flat HVSR ambient noise curve in which new peaks appeared when analyzing earthquake data. Considering that the energy associated with ambient noise is not comparable to that associated with a seismic event, it is possible that several sites considered as hard rock may have significant amplifications during the occurrence of an earthquake.
Based on these observations, we consider it is necessary to increase the density of the permanent seismic instrumentation in the TMVB. This will allow to have greater coverage, better quality records, and the possibility of carrying out more analyses using earthquake data, including the SSR technique. A greater permanent seismic instrumentation would allow obtaining isofrequency maps, which could be used in soil-structure interaction analyses and as a basis for estimating seismic hazards and risk assessment in the TMVB. This would also allow to better delimit the risk of populated zones based on local seismicity, considering the significant growth that the cities in the area have had in recent years.