Next Article in Journal
The Clean Label Concept: Novel Approaches in Gluten-Free Breadmaking
Next Article in Special Issue
A Case Study of a Large Unstable Mass Stabilization: “El Portalet” Pass at the Central Spanish Pyrenees
Previous Article in Journal
Advanced Crystal Plasticity Modeling of Multi-Phase Steels: Work-Hardening, Strain Rate Sensitivity and Formability
Previous Article in Special Issue
Detecting Areas Vulnerable to Flooding Using Hydrological-Topographic Factors and Logistic Regression
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
L. Francisco Pérez-Moreno
1,
Quetzalcoatl Rodríguez-Pérez
2,*,
F. Ramón Zúñiga
3,
Jaime Horta-Rangel
1,
M. de la Luz Pérez-Rea
1 and
Miguel A. Pérez-Lara
1
1
Facultad de Ingeniería, Universidad Autónoma de Querétaro, Querétaro 76010, Mexico
2
Consejo Nacional de Ciencia y Tecnología, Dirección Adjunta de Desarrollo Científico, Ciudad de México 03940, Mexico
3
Centro de Geociencias, Universidad Nacional Autónoma de México, Querétaro 76230, Mexico
*
Author to whom correspondence should be addressed.
Appl. Sci. 2021, 11(13), 6126; https://doi.org/10.3390/app11136126
Submission received: 28 April 2021 / Revised: 25 June 2021 / Accepted: 28 June 2021 / Published: 30 June 2021
(This article belongs to the Special Issue Geohazards: Risk Assessment, Mitigation and Prevention)

Abstract

:

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 amplifications 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 significant differences were observed in terms of the fundamental frequency. However, there were some stations with a flat average HVSR ambient noise curve that contradicted earthquake data results, which showed peaks at some frequencies. Our results are a first step towards categorizing the different site responses in the TMVB but in order to provide finer details, it is necessary to improve the actual monitoring conditions.

1. 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:
A ( f ) = S ( f ) × P ( f ) × H ( f )
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:
S S R ( f ) = | A i ( f ) | | A j ( f ) | = | H i ( f ) | | H j ( f ) | .
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 Horizontal-to-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 H i ( f ) and the vertical component A V i ( f ) , recorded in a single station, is calculated:
H V S R ( f ) = | A H i ( f ) | | A V i ( f ) | .
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.

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

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

3.1. 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].
Figure 3 shows the seismological stations analyzed in the Western sector of the TMVB. As can be observed, there is only one station belonging to a permanent seismological network. The other temporary stations are mostly concentrated in the south of this sector.
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.
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.
In the remaining stations, clear peaks of different amplitude can be observed, revealing different impedance contrasts between sediment layers and the bedrock: BAVA, JANU, MA44, MA45, ANIG, MA40, MA41, PAVE, GARC, ALPI, SANM (Figure 4), CDGZ, COLM, CANO, and ESPN (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.
The stations located on mafic and intermediate rocks from the Neogene and Quaternary have a heterogeneous response (ANIG, MA40, MA41, PAVE, GARC, ALPI, and SANM) (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.

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

3.3. 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].
Figure 8 shows the seismological stations analyzed in this sector. In this case, it was possible to analyze nine stations with earthquake data and to compare the results with the use of ambient noise. Particularly noteworthy are the stations belonging to the MASE temporary network, deployed between 2004 and 2007. Despite being numerous, due to the characteristics of the events that occurred during this period, only a few stations had earthquake records with the characteristics sought.
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.

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

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

Author Contributions

Conceptualization, L.F.P.-M., Q.R.-P. and F.R.Z.; methodology, L.F.P.-M.; formal analysis, L.F.P.-M.; investigation, L.F.P.-M. and M.A.P.-L.; writing—original draft preparation, L.F.P.-M.; writing—review and editing, Q.R.-P., F.R.Z., J.H.-R. and M.d.l.L.P.-R.; supervision, J.H.-R., M.d.l.L.P.-R. and M.A.P.-L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding. The APC was funded by L.F.P.-M.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

Special thanks to the Mexican National Council for Science and Technology (CONACYT) for the support provided to L. Francisco Pérez-Moreno during the realization of his Ph.D. studies and to Quetzalcoatl Rodríguez-Pérez in the Cátedras program (project 1126). SSN data was obtained by the Servicio Sismológico Nacional (México). Station maintenance, data acquisition, and distribution is thanks to its personnel. The authors thank Antonio de J. Mendoza Carvajal and Jorge A. Real Pérez for the maintenance and operation of the Geometry of the Cocos Plate (GECO) array. This network was financed by PAPIIT (projects IN110913, IN105816, and IN106119) and CONACYT (project 177676). We are also particularly grateful to Xyoli Pérez-Campos for providing access to the network data. We appreciate the contributions of the four anonymous reviewers of the initial version of this paper. Their valuable comments allowed us to improve this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Borcherdt, R.D. Effects of local geology on ground motion near San Francisco bay. Bull. Seismol. Soc. Am. 1970, 60, 29–61. [Google Scholar]
  2. Nakamura, Y. A method for dynamic characteristics estimation of subsurface using microtremor on the ground surface. Railw. Tech. Res. Inst. Q. Rep. 1989, 30, 25–33. [Google Scholar]
  3. Perron, V.; Gélis, C.; Froment, B.; Hollender, F.; Bard, P.Y.; Cultrera, G.; Cushing, E.M. Can broad-band earthquake site responses be predicted by the ambient noise spectral ratio? Insight from observations at two sedimentary basins. Geophys. J. Int. 2018, 215, 1442–1454. [Google Scholar] [CrossRef]
  4. Gosar, A. Study on the applicability of the microtremor HVSR method to support seismic microzonation in the town of Idrija (W Slovenia). Nat. Hazards Earth Syst. Sci. 2017, 17, 925–937. [Google Scholar] [CrossRef] [Green Version]
  5. Madiai, C.; Facciorusso, J.; Gargini, E. Numerical Modeling of Seismic Site Effects in a Shallow Alluvial Basin of the Northern Apennines (Italy). Bull. Seismol. Soc. Am. 2017, 107, 2094–2105. [Google Scholar] [CrossRef]
  6. Stanko, D.; Gülerce, Z.; Markušić, S.; Šalić, R. Evaluation of the site amplification factors estimated by equivalent linear site response analysis using time series and random vibration theory based approaches. Soil Dyn. Earthq. Eng. 2019, 117, 16–29. [Google Scholar] [CrossRef]
  7. Sun, Q.; Guo, X.; Dias, D. Evaluation of the seismic site response in randomized velocity profiles using a statistical model with Monte Carlo simulations. Comput. Geotech. 2020, 120, 103442. [Google Scholar] [CrossRef]
  8. Falcone, G.; Romagnoli, G.; Naso, G.; Mori, F.; Peronace, E.; Moscatelli, M. Effect of bedrock stiffness and thickness on numerical simulation of seismic site response. Italian case studies. Soil Dyn. Earthq. Eng. 2020, 139, 106361. [Google Scholar] [CrossRef]
  9. Singh, S.K.; Quintanar-Robles, L.; Arroyo, D.; Cruz-Atienza, V.C.; Espíndola, V.H.; Bello-Segura, D.I.; Ordaz, M. Lessons from a Small Local Earthquake (Mw 3.2) That Produced the Highest Acceleration Ever Recorded in Mexico City. Seismol. Res. Lett. 2020, 91, 3391–3406. [Google Scholar] [CrossRef]
  10. Ferrari, L.; Orozco-Esquivel, T.; Vlad, C.; Marina, M. The dynamic history of the Trans-Mexican Volcanic Belt and the Mexico subduction zone. Tectonophysics 2012, 522–523, 122–149. [Google Scholar] [CrossRef]
  11. Suárez, G.; Caballero-Jiménez, G.V.; Novelo-Casanova, D.A. Active Crustal Deformation in the Trans-Mexican Volcanic Belt as Evidenced by Historical Earthquakes During the Last 450 Years. Tectonics 2019, 38, 3544–3562. [Google Scholar] [CrossRef]
  12. Suárez, G.; Ruiz-Barón, D.; Chico-Hernández, C.; Zúñiga, F.R. Catalog of Preinstrumental Earthquakes in Central Mexico: Epicentral and Magnitude Estimations Based on Macroseismic Data. Bull. Seismol. Soc. Am. 2020, 110, 3021–3036. [Google Scholar] [CrossRef]
  13. Suter, M.; Carrillo-Martínez, M.; Quintero, O. Macroseismic Study of Shallow Earthquakes in the Central and Eastern Parts of the Trans-Mexican Volcanic Belt, Mexico. Bull. Seismol. Soc. Am. 1996, 86, 1952–1963. [Google Scholar]
  14. Sahakian, V.J.; Melgar, D.; Quintanar, L.; Ramírez-Guzmán, L.; Pérez-Campos, X.; Baltay, A. Ground motions from the 7 and 19 September 2017 Tehuantepec and Puebla-Morelos, Mexico, earthquakes. Bull. Seismol. Soc. Am. 2018, 108, 3300–3312. [Google Scholar] [CrossRef] [Green Version]
  15. Ordaz, M.; Singh, S.K. Source spectra and spectral attenuation of seismic waves from Mexican earthquakes, and evidence of amplification in the hill zone of Mexico City. Bull. Seismol. Soc. Am. 1992, 82, 24–43. [Google Scholar]
  16. García, D.; Singh, S.K.; Herráiz, M.; Pacheco, J.F.; Ordaz, M. Inslab earthquakes of Central Mexico: Q, source spectra, and stress drop. Bull. Seismol. Soc. Am. 2004, 94, 789–802. [Google Scholar] [CrossRef]
  17. Singh, S.K.; Iglesias-Mendoza, A.; García, D.; Pacheco, J.F.; Ordaz, M. Q of Lg waves in the central Mexican Volcanic Belt. Bull. Seismol. Soc. Am. 2007, 97, 1259–1266. [Google Scholar] [CrossRef] [Green Version]
  18. Lozano, L.; Herráiz, M.; Singh, S.K. Site effect study in central Mexico using H/V and SSR techniques: Independence of seismic site effects on source characteristics. Soil Dyn. Earthq. Eng. 2009, 29, 504–516. [Google Scholar] [CrossRef]
  19. Clemente-Chávez, A.; Zúñiga, F.R.; Lermo, J.; Figueroa-Soto, A.; Valdés-González, C.; Montiel, M.; Chávez, O.; Arroyo, M. On the behavior of site effects in central Mexico (the Mexican volcanic belt-MVB), based on records of shallow earthquakes that occurred in the zone between 1998 and 2011. Nat. Hazards Earth Syst. Sci. 2014, 14, 1391–1406. [Google Scholar] [CrossRef] [Green Version]
  20. Zúñiga, F.R.; Lacan, P.; Rodríguez-Pérez, Q.; Márquez-Ramírez, V.H. Temporal and spatial evolution of instrumented seismicity in the Trans-Mexican Volcanic Belt. J. S. Am. Earth Sci. 2020, 98, 102390. [Google Scholar] [CrossRef]
  21. Ruggieri, S.; Tosto, C.; Rosati, G.; Uva, G.; Ferro, G.A. Seismic Vulnerability Analysis of Masonry Churches in Piemonte after 2003 Valle Scrivia Earthquake: Post-event Screening and Situation 17 Years Later. Int. J. Archit. Herit. 2020, 1–29. [Google Scholar] [CrossRef]
  22. Godínez-Domínguez, E.A.; Tena-Colunga, A.; Pérez-Rocha, L.E.; Archundia-Aranda, H.I.; Gómez-Bernal, A.; Ruiz-Torres, R.P.; Escamilla-Cruz, J.L. The September 7, 2017 Tehuantepec, Mexico, earthquake: Damage assessment in masonry structures for housing. Int. J. Disaster Risk Reduct. 2021, 56, 1–29. [Google Scholar] [CrossRef]
  23. Catálogo de Sismos. Available online: http://www2.ssn.unam.mx:8080/catalogo/ (accessed on 16 October 2020).
  24. Search Earthquake Catalog. Available online: https://earthquake.usgs.gov/earthquakes/search/ (accessed on 16 October 2020).
  25. ISC Bulletin: Event Catalogue Search. Available online: http://www.isc.ac.uk/iscbulletin/search/catalogue/ (accessed on 16 October 2020).
  26. GEOSCOPE. French Global Network of Broad Band Seismic Stations; Institut de Physique du Globe de Paris (IPGP) and Ecole Et Observatoire Des Sciences De La Terre De Strasbourg (EOST). [CrossRef]
  27. Meso America Subduction Experiment. Caltech. Dataset. [CrossRef]
  28. Pérez-Campos, X. MASE: Undergraduate Research and Outreach as Part of a Large Project. Seismol. Res. Lett. 2008, 92, 232–236. [Google Scholar] [CrossRef] [Green Version]
  29. West, M. The Colima Deep Seismic Experiment: Imaging the Magmatic Root of Colima Volcano. Int. Fed. Digit. Seismogr. Netw. 2006. [Google Scholar]
  30. Yang, T.; Grand, S.P.; Wilson, D.; Guzman-Speziale, M.; Gomez-Gonzalez, J.M.; Dominguez-Reyes, T. Seismic structure beneath the Rivera subduction zone from finite-frequency seismic tomography. J. Geophys. Res. Solid Earth 2009, 114, B01302. [Google Scholar] [CrossRef]
  31. Rodríguez-Domínguez, M.; Pérez-Campos, X.; Montealegre-Cázares, C.; Clayton, R.W.; Cabral-Cano, E. Crustal structure variations in south-central Mexico from receiver functions. Geophys. J. Int. 2019, 219, 2174–2186. [Google Scholar] [CrossRef] [Green Version]
  32. Wathelet, M.; Chatelain, J.L.; Cornou, C.; Giulio, G.D.; Guillier, B.; Ohrnberger, M.; Savvaidis, A. Geopsy: A user-friendly open-source tool set for ambient vibration processing. Seismol. Res. Lett. 2020, 91, 1878–1889. [Google Scholar] [CrossRef]
  33. Konno, K.; Ohmachi, T. Ground-motion characteristics estimated from spectral ratio between horizontal and vertical components of microtremor. Bull. Seismol. Soc. Am. 1998, 88, 228–241. [Google Scholar]
  34. SESAME Team. Guidelines for the Implementation of the H/V Spectral Ratio Technique on Ambient Vibrations: Measurements, Processing and Interpretation; SESAME European Research Project. 2004. Available online: http://sesame.geopsy.org/Papers/HV_User_Guidelines.pdf (accessed on 15 August 2020).
  35. Reed, J.C.; Wheeler, J.O.; Tucholke, B.E. Decade of North American Geology Geologic Map of North America—Perspectives and Explanation; The Geological Society of America: Boulder, CO, USA, 2015. [Google Scholar]
  36. Suter, M. The A.D. 1567 Mw 7.2 Ameca, Jalisco, Earthquake (Western Trans-Mexican Volcanic Belt): Surface Rupture Parameters, Seismogeological Effects, and Macroseismic Intensities from Historical Source. Bull. Seismol. Soc. Am. 2015, 105, 646–656. [Google Scholar] [CrossRef]
  37. Suter, M. The 1563 MI 8 Puerto de la Navidad Subduction-Zone and 1567 Mw 7.2 Ameca Crustal Earthquakes (Western Mexico): New Insights from Sixteenth-Century Sources. Seismol. Res. Lett. 2019, 90, 366–375. [Google Scholar] [CrossRef]
  38. Suter, M. The 2 October 1847 MI 5.7 Chapala Graben Triggered Earthquake (Trans-Mexican Volcanic Belt, West-Central Mexico): Macroseismic Observations and Hazard Implications. Seismol. Res. Lett. 2017, 89, 35–46. [Google Scholar] [CrossRef]
  39. Lacan, P.; Aguirre-Díaz, G.; Audin, L.; Perea, H.; Baize, S.; Zúñiga, F.R.; Ortuño, M. Sedimentary evidence of historical and prehistorical earthquakes along the Venta de Bravo Fault System, Acambay Graben (Central Mexico). Sediment. Geol. 2018, 365, 62–77. [Google Scholar] [CrossRef]
  40. Rodríguez-Pérez, Q.; Zúñiga, F.R. Seismicity characterization of the Maravatío-Acambay and Actopan regions, central Mexico. J. S. Am. Earth Sci. 2017, 76, 264–275. [Google Scholar] [CrossRef]
  41. Suárez, G.; Novelo-Casanova, D.A. A Pioneering Aftershock Study of the Destructive 4 January 1920 Jalapa, Mexico, Earthquake. Seismol. Res. Lett. 2018, 89, 1894–1899. [Google Scholar] [CrossRef]
Figure 1. Location of the Trans-Mexican Volcanic Belt and interaction of tectonic plates in Mexico.
Figure 1. Location of the Trans-Mexican Volcanic Belt and interaction of tectonic plates in Mexico.
Applsci 11 06126 g001
Figure 2. Location of epicenters and seismological stations analyzed in this study. Cyan circles: Crustal earthquakes; orange circles: Inslab earthquakes. Red triangles: MARS temporary network stations; blue triangles: CODEX temporary network stations; green triangles: MASE temporary network stations; purple triangles: GECO temporary network stations; black triangles: SSN permanent network stations; yellow triangle: CGEO permanent station; and white triangle: GEOSCOPE permanent station. Dashed line: TMVB.
Figure 2. Location of epicenters and seismological stations analyzed in this study. Cyan circles: Crustal earthquakes; orange circles: Inslab earthquakes. Red triangles: MARS temporary network stations; blue triangles: CODEX temporary network stations; green triangles: MASE temporary network stations; purple triangles: GECO temporary network stations; black triangles: SSN permanent network stations; yellow triangle: CGEO permanent station; and white triangle: GEOSCOPE permanent station. Dashed line: TMVB.
Applsci 11 06126 g002
Figure 3. Geological units in the western sector of the TMVB [35]. Red triangles: MARS temporary network stations; blue triangles: CODEX temporary network stations; and black triangle: SSN permanent network station. Continuous line: Main fault systems in the area; dashed line: TMVB.
Figure 3. Geological units in the western sector of the TMVB [35]. Red triangles: MARS temporary network stations; blue triangles: CODEX temporary network stations; and black triangle: SSN permanent network station. Continuous line: Main fault systems in the area; dashed line: TMVB.
Applsci 11 06126 g003
Figure 4. Average HVSR curves for stations located in the western sector of the TMVB. Blackline: HVSR results using ambient noise records; redline: HVSR results from earthquake records.
Figure 4. Average HVSR curves for stations located in the western sector of the TMVB. Blackline: HVSR results using ambient noise records; redline: HVSR results from earthquake records.
Applsci 11 06126 g004
Figure 5. Average HVSR curves for stations located in the western sector of the TMVB.
Figure 5. Average HVSR curves for stations located in the western sector of the TMVB.
Applsci 11 06126 g005
Figure 6. Geological units in the central sector of the TMVB [35]. Red triangles: MARS temporary network stations; black triangles: SSN permanent network stations. Continuous line: Main fault systems in the area; dashed line: TMVB.
Figure 6. Geological units in the central sector of the TMVB [35]. Red triangles: MARS temporary network stations; black triangles: SSN permanent network stations. Continuous line: Main fault systems in the area; dashed line: TMVB.
Applsci 11 06126 g006
Figure 7. Average HVSR curves for stations located in the central sector of the TMVB.
Figure 7. Average HVSR curves for stations located in the central sector of the TMVB.
Applsci 11 06126 g007
Figure 8. Geological units in the central-eastern sector of the TMVB [35]. Green triangles: MASE temporary network stations; black triangles: SSN permanent network stations; yellow triangle: CGEO permanent station; and white triangle: GEOSCOPE permanent station. Continuous line: Main fault systems in the area; dashed line: TMVB.
Figure 8. Geological units in the central-eastern sector of the TMVB [35]. Green triangles: MASE temporary network stations; black triangles: SSN permanent network stations; yellow triangle: CGEO permanent station; and white triangle: GEOSCOPE permanent station. Continuous line: Main fault systems in the area; dashed line: TMVB.
Applsci 11 06126 g008
Figure 9. Average HVSR curves for some stations located in the central-eastern sector of the TMVB. Blackline: HVSR results using ambient noise records; redline: HVSR results using earthquake records.
Figure 9. Average HVSR curves for some stations located in the central-eastern sector of the TMVB. Blackline: HVSR results using ambient noise records; redline: HVSR results using earthquake records.
Applsci 11 06126 g009
Figure 10. Average HVSR curves for some stations located in the central-eastern sector of the TMVB. Blackline: HVSR results using ambient noise records; redline: HVSR results using earthquake records.
Figure 10. Average HVSR curves for some stations located in the central-eastern sector of the TMVB. Blackline: HVSR results using ambient noise records; redline: HVSR results using earthquake records.
Applsci 11 06126 g010
Figure 11. Geological units in the eastern sector of the TMVB [35]. Purple triangles: GECO temporary network stations; black triangles: SSN permanent network stations. Dashed line: TMVB.
Figure 11. Geological units in the eastern sector of the TMVB [35]. Purple triangles: GECO temporary network stations; black triangles: SSN permanent network stations. Dashed line: TMVB.
Applsci 11 06126 g011
Figure 12. Average HVSR curves for stations located in the Eastern sector of the TMVB. Blackline: HVSR results using ambient noise records; redline: HVSR results using earthquake records.
Figure 12. Average HVSR curves for stations located in the Eastern sector of the TMVB. Blackline: HVSR results using ambient noise records; redline: HVSR results using earthquake records.
Applsci 11 06126 g012
Table 1. Seismological networks that were analyzed in this study. The continuous data availability is variable depending on each station.
Table 1. Seismological networks that were analyzed in this study. The continuous data availability is variable depending on each station.
NetworkPeriod
SSN2002–present
CGEO2003–present
GEOSCOPE1990–present
MASE2004–2007
CODEX2006–2008
MARS2006–2007
GECO2013–2018
Table 2. A database of crustal and inslab earthquakes with records that met the considered criteria. The value of magnitude is the one reported by SSN.
Table 2. A database of crustal and inslab earthquakes with records that met the considered criteria. The value of magnitude is the one reported by SSN.
DateOrigin TimeMagnitudeDepthLatitudeLongitude
dd/mm/yyyyhh:mm:ss (km)( )( )
08/03/199314:37:413.31119.25−98.93
02/11/199415:31:213.7719.41−98.90
27/03/199608:05:153.4819.77−99.01
17/04/199611:19:103.51019.21−98.94
06/09/199606:55:004.0519.21−98.36
03/07/199719:39:322.8719.25−99.44
25/12/199701:29:104.11118.97−98.59
15/03/199802:01:494.2519.50−100.26
26/10/199802:35:103.0619.06−99.17
26/10/199807:50:063.3419.06−99.20
31/12/199806:44:053.8218.98−98.65
04/01/199923:21:593.5320.11−98.98
19/05/199920:23:063.1519.19−98.96
22/05/199906:57:263.1419.18−98.97
07/06/199900:24:143.3519.17−98.96
15/06/199920:42:047.06318.13−97.54
20/06/199910:25:214.0319.31−98.48
07/12/199913:53:493.4319.22−98.92
04/03/200020:55:234.1418.84−98.57
12/03/200003:44:484.1520.10−99.29
06/07/200008:36:243.3819.18−98.93
06/07/200013:15:343.4519.19−98.96
21/07/200006:13:396.04818.09−98.97
14/10/200002:44:483.11219.36−99.19
15/10/200007:50:432.7519.39−99.11
01/03/200116:26:423.54618.98−99.82
13/06/200103:57:433.21219.27−99.42
14/09/200117:13:152.91419.31−99.30
15/11/200122:18:223.5419.56−99.16
09/05/200214:25:433.71919.49−99.01
04/02/200310:59:034.1218.92−98.51
15/12/200310:39:384.0420.35−99.07
13/09/200420:58:343.11119.42−99.17
16/04/200522:55:253.81619.43−99.00
12/05/200508:06:253.43719.20−98.97
07/08/200503:25:094.0719.74−98.61
16/10/200514:12:363.51419.30−99.20
23/11/200523:11:263.62019.35−98.94
24/01/200612:59:383.7520.34−99.22
28/02/200623:58:493.4219.35−99.02
12/03/200601:41:323.6219.17−98.95
12/03/200601:47:213.4519.18−98.96
05/05/200615:24:063.5119.17−98.95
27/06/200615:40:104.01619.26−102.32
14/07/200605:15:283.77220.19−98.87
18/12/200603:38:313.94420.92−104.47
25/01/200714:38:223.73019.18−99.16
30/01/200719:00:103.7920.35−103.87
14/05/200708:23:274.41221.36−104.75
30/05/200720:42:284.0319.21−99.45
22/06/200713:36:524.22519.16−96.96
20/10/200719:25:543.77519.47−102.20
20/10/200702:37:273.92019.95−101.97
21/10/200711:29:533.79319.53−102.17
23/03/200902:53:363.94618.92−97.01
22/05/200919:24:185.76218.11−98.46
29/11/200906:34:124.0519.35−103.76
17/04/201007:03:294.0220.45−99.04
18/04/201018:26:254.07919.27−97.46
18/05/201005:29:094.0520.35−98.92
18/05/201005:32:373.41020.37−98.96
03/10/201023:07:114.0619.48−103.52
19/12/201002:46:454.11619.03−97.21
28/09/201110:51:184.12319.76−96.64
18/05/201203:07:584.41.620.26−103.44
05/10/201200:14:574.216.119.45−102.26
08/02/201301:24:583.86.520.12−100.48
21/06/201309:41:334.0218.52−98.74
18/09/201300:26:153.3519.62−96.79
14/12/201322:33:404.1104.619.45−103.71
24/12/201315:27:073.6319.52−103.80
16/01/201407:23:573.4919.27−97.26
09/05/201407:50:493.7318.98−97.27
19/07/201411:57:233.5819.17−98.97
04/08/201421:10:034.059.218.18−97.83
08/11/201402:46:503.6519.51−96.78
01/12/201408:50:063.4219.35−99.22
06/05/201504:06:463.71019.16−97.17
09/08/201515:17:173.7518.96−98.21
16/08/201514:29:173.42019.65−97.55
25/10/201502:52:433.777.319.26−97.30
07/11/201502:29:504.1519.73−96.70
26/12/201514:08:364.0519.55−102.11
16/01/201606:28:463.142.619.65−96.80
27/01/201621:22:444.169.518.28−97.33
28/01/201611:41:493.44419.45−96.72
08/02/201621:16:074.66.919.66−97.35
13/03/201619:13:503.755.319.31−97.03
01/04/201622:34:333.758.618.27−97.78
01/04/201605:06:513.1519.63−96.64
01/04/201605:39:393.44.119.68−96.73
11/05/201622:35:204.8820.81−103.52
15/05/201601:43:483.9520.82−103.47
11/06/201619:06:143.912.518.85−97.19
16/06/201600:56:293.31019.95−96.77
24/06/201600:22:543.42019.23−97.40
27/06/201619:12:253.41019.95−96.76
12/07/201617:33:293.5520.77−103.45
19/07/201604:34:433.956.718.29−97.38
28/07/201618:21:063.9520.73−103.45
08/08/201615:57:083.82019.44−96.72
29/08/201616:24:173.8318.95−98.58
26/09/201603:44:443.610.119.51−96.45
15/10/201601:51:484.151.318.87−96.99
26/02/201702:11:294.023.719.88−96.85
01/06/201713:36:593.8519.24−97.04
09/07/201709:03:333.530.719.69−96.81
04/08/201701:07:363.51019.50−96.43
10/09/201702:54:122.79.419.31−99.18
19/09/201718:14:407.151.218.33−98.68
29/09/201717:02:233.850.718.35−98.66
01/11/201720:48:044.352.818.25−98.60
06/11/201709:38:093.1319.19−98.44
09/11/201710:12:072.08.519.38−99.19
18/11/201707:21:543.632.418.64−98.48
24/12/201712:25:123.43.118.82−98.60
08/06/201818:09:343.911.519.95−96.75
30/07/201822:15:543.924.619.27−97.34
30/03/201922:09:223.9519.44−96.71
14/07/201919:37:063.83719.21−97.23
23/07/201919:03:204.03219.27−97.43
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pérez-Moreno, L.F.; Rodríguez-Pérez, Q.; Zúñiga, F.R.; Horta-Rangel, J.; Pérez-Rea, M.d.l.L.; Pérez-Lara, M.A. Site Response Evaluation in the Trans-Mexican Volcanic Belt Based on HVSR from Ambient Noise and Regional Seismicity. Appl. Sci. 2021, 11, 6126. https://doi.org/10.3390/app11136126

AMA Style

Pérez-Moreno LF, Rodríguez-Pérez Q, Zúñiga FR, Horta-Rangel J, Pérez-Rea MdlL, Pérez-Lara MA. Site Response Evaluation in the Trans-Mexican Volcanic Belt Based on HVSR from Ambient Noise and Regional Seismicity. Applied Sciences. 2021; 11(13):6126. https://doi.org/10.3390/app11136126

Chicago/Turabian Style

Pérez-Moreno, L. Francisco, Quetzalcoatl Rodríguez-Pérez, F. Ramón Zúñiga, Jaime Horta-Rangel, M. de la Luz Pérez-Rea, and Miguel A. Pérez-Lara. 2021. "Site Response Evaluation in the Trans-Mexican Volcanic Belt Based on HVSR from Ambient Noise and Regional Seismicity" Applied Sciences 11, no. 13: 6126. https://doi.org/10.3390/app11136126

APA Style

Pérez-Moreno, L. F., Rodríguez-Pérez, Q., Zúñiga, F. R., Horta-Rangel, J., Pérez-Rea, M. d. l. L., & Pérez-Lara, M. A. (2021). Site Response Evaluation in the Trans-Mexican Volcanic Belt Based on HVSR from Ambient Noise and Regional Seismicity. Applied Sciences, 11(13), 6126. https://doi.org/10.3390/app11136126

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop