## 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 Sismologico 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 Sismologico Nacional (SSN) [23], Centro de Geociencias de la Universidad Nacional Autonoma de Mexico (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 \(\leq\) H \(\leq\) 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 [PERSON] 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], \begin{table}
\begin{tabular}{c c} \hline \hline
**Network** & **Period** \\ \hline SSN & 2002–present \\ CGEO & 2003–present \\ GEOSCOPE & 1990–present \\ MASE & 2004–2007 \\ CODEX & 2006–2008 \\ MARS & 2006–2007 \\ GECO & 2013–2018 \\ \hline \hline \end{tabular}
\end{table} 
Table 1: Seismological networks that were analyzed in this study.