Valley of Mexico site response from non-reference site approach

Abstract

We present a model for Valley of Mexico site response that is referenced to a 30 m time-averaged shear wave velocity (V_S₃₀) of 760 m/s. This reference condition ensures compatibility with regional ground motion models (GMM) for Central America and Mexico (CAM), which were developed from an expanded version of the NGA-Subduction (NGA-Sub) database as described in a companion paper. The site response model is derived using ground motion data from 89 sites within the Valley of Mexico and uses V_S₃₀ and soft soil depth as independent variables, the former of which was evaluated using shear wave velocity from literature for 56 sites (10 co-located with ground motion stations). As a consequence of the unique Lake Texcoco geology, the model predicts far greater levels of site response than those provided by global ergodic models, and notably, larger site response than prior models for the Valley of Mexico that had utilized a reference site approach using a medium-stiff soil site as the reference site. Linear amplification increases with decreasing V_S₃₀ in a period-dependent manner that captures strong resonance effects associated with variable depths of soft lacustrine soils. Site period is found to not appreciably improve predictions, due to strong negative correlation with V_S₃₀. Differential depth (depth minus V_S₃₀-conditioned mean depth) produces repeatable trends of reduced short period amplification and increased long period amplification for especially deep sites, the effects of which are modeled. We identify nonlinear site response for high-frequency ground motions using variable-amplitude reference site ground motions. The provided models improve upon current practice for Mexico because the reference GMM applies across CAM, V_S₃₀ and soft soil depth are used in lieu of site classes, and nonlinearity is incorporated.

Keywords

Soft soil site response ground motion model site period nonlinearity Mexico City

Introduction

Site response is fundamentally a ground motion difference between two site conditions. One is a reference condition, generally involving relatively competent materials, while the other is the site condition at a particular site of interest. Reference site and non-reference site approaches can be used to measure these ground motion differences (Field and Jacob, 1995). Site response in the Valley of Mexico, which contains Mexico City, has typically been evaluated in past work using reference site approaches (Borcherdt, 1970, 1994), in which ground motions on various site conditions are normalized by ground motions from reference sites. The site amplification that is derived depends on the attributes of the reference site or sites, which have usually been taken (for applications outside of the Valley of Mexico) as bedrock sites. Non-reference site approaches define the reference condition not for a single site, but for a parameterized site condition as derived from seismological inversions or as used in ground motion models (GMMs); large values of time-averaged 30 m shear wave velocity (V_S₃₀ = 760–1000 m/s) are typically selected.

Current practice is that ground motion estimations for sites in the Valley of Mexico are based on a reference site approach having two basic steps (e.g. Ordaz et al., 2024; Ordaz et al., 1994; Reinoso and Ordaz, 1999). The first step pertains to the estimation of reference site ground motions, which are for a specific location on the UNAM campus in the hill zone adjacent to the Valley of Mexico lakebed (referred to as the CU site). The second step estimates site response for various locations in the Valley of Mexico relative to the reference site. There are several drawbacks to this traditionally applied approach, including the need for a GMM for a single reference site, which necessarily produces data sampling problems (the available data for a single site are limited) and divorces the ground motion estimation for the Valley of Mexico from locations elsewhere in Mexico, which use GMMs derived from much larger databases.

We revisit the well-studied problem of Valley of Mexico site response using a non-reference site approach, which to our knowledge has not been applied previously for this important locale where site response has produced structural damage (e.g. Mayoral et al., 2019a; Singh et al., 2018). The benefit of this alternative approach is that it can produce site factors that are consistent with regional GMMs, thus avoiding the need for distinct ground motion estimation procedures in Mexico City and the remainder of Mexico. Such approaches have been applied elsewhere to derive regional (e.g. Parker and Stewart, 2022) and local (Wang et al., 2022) site response models in subduction zone regions. We first review the lakebed geology and zonation procedures used to characterize spatial variations in site conditions. We then synthesize site data for the purpose of developing site parameters for ground motion stations and for estimating values of site parameters in geotechnical zones. We describe current procedures used for ground motion estimation in the Valley of Mexico and provide example results. We then describe the application of a non-reference site approach with available ground motion data, from which a subregional site response model (i.e. model specific to the Valley of Mexico) is developed that is conditional on V_S₃₀. The model is compared to site response estimated using currently applied reference site procedures—large differences are observed that are related to appreciable site response at the reference site—and is validated against a testing data set from the 1985 M8.0 Michoacán earthquake independent of the training data. We conclude with recommendations on how the present results can be utilized to improve ground motion hazard estimation in Mexico City. A previous version of the material presented in this article appeared in the dissertation of the first author (Contreras, 2022) and a subsequent report based on that dissertation (Contreras et al., 2023a). Portions of the work presented here and in a companion paper (Contreras et al., 2025) have also been the subject of conference presentations (Contreras et al., 2023b; Contreras et al., 2024).

Geology and zonation

Mexico City is located within the eastern sector of the Trans-Mexican Volcanic Belt, north of the volcanic front (i.e. on the backarc side). Figure 1 depicts the geology of the subregion, as adapted from Flores-Estrella et al. (2007) and Arce et al. (2019). Zones delineating different geotechnical conditions are shown in Figure 2. These zones were originally derived from geologic maps and geotechnical logs (Marsal and Mazari, 1969; Romo et al., 1988). As indicated by Lermo et al. (2020), the zones were refined by Normas Técnicas Complementarias para el diseño por Sismos (NTCS) (2004) in consideration of site period maps by Lermo et al. (1988) and Lermo and Chávez-García (1993, 1994a, 1994b) and have been used for site classification purposes.

Figure 1.

Map showing Valley of Mexico (VM) region including surface geology (adapted from Flores-Estrella et al., 2007 and Arce et al., 2019) and locations of ground motion stations and seismic velocity profiles considered in this research. The location of the reference CU site is shown.

Figure 2.

Map of Valley of Mexico (VM) showing zones used for site classification per NTCS (2004). Also shown are locations of selected recording sites used to analyze site response. Locations of seven sites located in six different zones are indicated, which are subsequently used to provide example results.

Within the Valley of Mexico, the lakebed contains lacustrine sediments that are remnants of ancient Lake Texcoco. The lake occupied the entire Valley of Mexico between the Pleistocene epoch and the last glacial period (about 11,000 years ago). The lake zone is indicated in Figure 1 by the Ql unit and includes the Texcoco Lake (while this was once encompassing the entire Valley of Mexico, as used here it refers to a smaller lake in the post-Aztec era), Xochimilco Lake, and Chalco Lake (Romo et al., 1988). These smaller lakes existed as recently as 1245 (upon arrival of the Aztecs), after which the lakes were impacted by human activity (Alcocer and Williams, 1996). Geotechnical conditions in these lake zones consist of a surface layer of desiccated alluvial deposits overlying thick and soft lacustrine clay interbedded with thin seams of sands, silty sands, volcanic glass, and fossils, which is locally known as the Upper Clay formation. Underlying this formation are hard deposits of silty, weakly cemented sands, a deeper clay stratum known as the Lower Clay formation, and compact lacustrine cemented silty sands and gravels. The Xochimilco–Chalco Lake consists of a clay deposit (somewhat stiffer than the Upper Clay formation in Texcoco) with interbedded seams of silty sands, silts, and sands. This clayey deposit is underlain by a basalt layer (lava flow). The lakebed area has four subzones (IIIa–IIId in Figure 2) that are differentiated in part by the period of the soil column (1.0–1.5 s, 1.5–2.5 s, 2.5–3.5 s, >3.5 s, respectively) (NTCS, 2004).

As shown in Figure 1, Quaternary alluvium occurs between the lake deposits and hill areas to the west and south, producing a geotechnical transition zone (denoted Zone II by NTCS, 2004; shown in yellow in Figure 2). These alluvial deposits have variable sequences of firm soils, sands, silty sands, and soft clays. The hill areas west of Lake Texcoco consist of Tertiary age epiclastic rock and alluvial fans (Flores-Estrella et al., 2007). Geotechnical materials encountered near the surface generally consist of silty sands with gravels and cemented tuffs. To the west and south of Xochimilco–Chalco Lake, Quaternary basaltic, andesitic, and epiclastic deposits up to 20 m in thickness overlie these formations (Flores-Estrella et al., 2007; Arce et al., 2019). These volcanic deposits were derived from a series of volcanoes in the region, including Iztaccíhuatl, Popocatépetl, La Malinche, Ajusco, and Nevado de Toluca. As shown in Figure 2, the hill zone is denoted as Zone I by NTCS (2004).

Site parameters

Data compilation

We have compiled available information on shear wave velocity (V_S) profiles, site periods (T_S), and depths to a firm horizon (Lower Coarse Grain layer, $z_{LGC}$ ) in the Valley of Mexico. The V_S profiles and T_S data are from the following technical reports, research papers, and student dissertations: Cárdenas-Soto et al. (2016), Chamorro (2016), De La Rosa et al. (2013), Díaz-Rodríguez (2003), Flores-Guzmán et al. (2014), Luna (2012), Mayoral et al. (2011), Mayoral et al., (2016), Mayoral et al. (2017), Mayoral et al. (2019b), Pérez (2017), Ramos-Martínez et al. (1997), Romo and Garcia (2003), Seed et al. (1987), Stephenson and Lomnitz (2005), Vergara-Huerta and Aguirre (2020), Wood et al., (2019, 2024, 2023). The $z_{LGC}$ information is from a contour map provided by Juárez-Camarena et al. (2016).

Velocity data come from invasive (downhole, crosshole, suspension logging, seismic cone penetration) and noninvasive methods (Multi-channel Analysis of Surface Waves—MASW, Microtremor Array Measurements—MAM, Modified microtremor seismic method using SPatial AutoCorrelation—MSPAC, seismic interferometry—SI, seismic refraction). Site period data are derived from microtremor-based measurements of Horizontal-to-Vertical Spectral Ratios (mHVSR), the lowest frequency peaks of which reveal the fundamental mode site period (Nogoshi and Igarashi, 1970, 1971; Field and Jacob, 1993, 1995; Bonilla et al., 1997, 2002; Cadet et al., 2012; Satoh et al., 2001; Theodulidis et al., 1996). With the exception of 23 sites from Wood et al. (2019, 2024, 2023), site periods are provided without the source mHVSR curves, which means that the procedures used to select site periods may be inconsistent and other parameters (such as the height of peaks) are not available. For sites with V_S profiles that reach a large impedance contrast (i.e. below the lower clay formation in Texcoco Lake or the basalt layer in Xochimilco–Chalco Lake), site period was also derived from V_S profiles using the Simplified Rayleigh Method as given in Urzúa et al. (2017) using soil layers generally down to $z_{LGC}$ .

Table A1 in the electronic supplement materials presents the compiled information for sites with a V_S profile, including the measurement locations, zones (from Figure 2), V_S₃₀ values, derived site periods co-located with the V_S profile where available, testing methods, $z_{LGC}$ values, and references. Two columns for site period are provided, one for results of mHVSR testing and one for period derived from the V_S profile. Figure 1 shows the locations of these sites.

Figure 3 shows box and whisker plots of V_S₃₀ and site period by zone. Both sources of site period (measured and derived from the V_S profile) are used; for sites with both types of site period, the mHVSR result was used. Where an mHVSR measurement is not co-located with the ground motion site, we used values of site period from the large compilation of Lermo and Chávez-García (1994a) (409 sites within the Valley of Mexico). Median and mean values of V_S₃₀ and site period are indicated by the blue and red lines in Figure 3, respectively. The lower and upper quartiles are indicated by the boxes and the extremes are indicated by the whiskers (extreme horizontal lines). Results in Figure 3 show that, on average, the three main zones (i.e. Zone I-Hills, Zone II-Transition, and Zone III-Lake) have statistically distinct values of V_S₃₀, as judged by non-overlapping boxes. The subdivision of Zone III shows significant differences between V_S₃₀ values in Zone IIIa and in the other zones (IIIb, c, d); however, the differences between these latter zones are lower and not necessarily significant. In the case of site periods, the zones have statistically distinct site periods. An exception is Zones IIIc–d, which present a similar median. Zone IIIc shows a large variability in comparison with the other zones, which may reflect natural variability or may be a consequence of data quality; further research would be needed to assess and potentially improve site period characterization in this area.

Figure 3.

Box and whisker plots of (a) V_S₃₀ and (b) site period ( $T_{S}$ ) for characterized sites in the Valley of Mexico (using data from Table A1). The number of sites within each zone are indicated.

Figure 4a shows the relationship between V_S₃₀ and site period. The negative correlation is strong, with short-period sites having high V_S₃₀, and V_S₃₀ decreasing with increasing $T_{S}$ . This result was expected given that about 57% of the data shown in Figure 4a are from measured V_S profiles that were constrained using the measured mHVSR site periods. Figure 4b shows a similarly strong negative correlation between V_S₃₀ and $z_{LGC}$ , which is considered below in model development.

Figure 4.

Correlation between V_S₃₀ and (a) site period ( $T_{S}$ ) and (b) depth to lower coarse grained layer ( $z_{LCG}$ ) for sites in the Valley of Mexico. Parameter ρ is the correlation coefficient.

V_S30 assignments

We apply protocols for characterizing V_S₃₀ for ground motion sites as described by Ahdi et al. (2022). That study classified different methods of V_S₃₀ prediction for sites without seismic velocity measurements. The data compilation in Table A1 was not available at the time the NGA-Sub database was assembled. Hence, while the general approach is the same as given by Ahdi et al. (2022), the assigned site parameters are different for sites that were already in the database, and many new sites have been added to the database.

The database contains 97 ground motion sites in the Valley of Mexico. Of these, 17 are Code 0 sites per Ahdi et al. (2022; their Table 5), meaning that V_S₃₀ was computed from a V_S profile, which is taken as the exponent of the mean of the log values ( $μ_{lnV}$ ) and the natural log standard deviation ( $σ_{lnV}$ ) is taken as 0.1. To assign a V_S profile to a ground motion site, we require separation distances ≤300 m and both locations to be in the same zone.

For the 80 sites without a V_S profile, we assigned $μ_{lnV}$ based on the zone the site is located in. This is an example of a region-specific prediction model, which is indicated as Code 2MC (Code 2 is as given in Table 5 of Ahdi et al. (2022); the letters “MC” are added for this Mexico City–specific model, which was not available when the NGA-Sub database was developed). The natural log means and standard deviation for each zone are given in Table 1, which are taken from the results in Figure 3.

Table 1.

Summary of V_S₃₀ parameters developed for Valley of Mexico stations without a V_S profile.

Zone	# sites with a V_S profile	$μ_{lnV}$ (m/s)	$σ_{lnV}$	# strong motion stations
I	7	331.3	0.296	20
II	10	172.3	0.394	13
IIIa	19	110.5	0.214	10
IIIb	15	67.1	0.234	21
IIIc	5	56.9	0.311	13
IIId	7	51.5	0.157	12

State of practice for ground motion prediction

Current seismic design in Mexico City uses response spectra that are derived using a two-step procedure:

Fourier amplitude spectra (FAS) are derived for a reference site in the hill zone (Ciudad Universitaria, CU site).

Reference site FAS are then modified using transfer functions for different locations and the resulting site FAS are converted to response spectra using random vibration theory.

These procedures are implemented for applications in an online software package, SASID (System of Seismic Actions for Design; https://sasid.unam.mx, last accessed October 2022, registration is required).

Four GMMs have been developed for the CU reference site (Arroyo et al., 2024; Jaimes et al., 2015; Jaimes et al., 2006; Ordaz et al., 1994). The number of recordings that are considered are equal to the number events and range from 21 to 23 for the first three models and 67 for the fourth (40 for “coastal” events, mainly interface, and 27 for intermediate depth events). These GMMs are referred to here as “single-site GMMs” to differentiate them from conventional GMMs that are derived considering data from many sites. It is currently unclear whether the GMM that is applied in the SASID software is a modified version of Ordaz et al. (1994) or the newer Arroyo et al. (2024) model. Contreras et al. (2023a) provided further details on the single-site GMMs and how they are applied to derive reference site spectra.

The transfer functions applied to modify the reference site FAS use data from recording sites in the Valley of Mexico and the CU reference site. The precise manner in which these calculations are performed in the current version of SASID is not well documented. The available documentation is provided by Reinoso and Ordaz (1999) and Ordaz et al. (2024), who describe procedures, whereby transfer functions are developed for a series of sites within each of the zones defined in Figure 2 (hill, transition, lakebed). As described by Reinoso and Ordaz (1999) and Ordaz et al. (2024), transfer functions are computed for each instrument location, and for locations between instruments, an interpolation procedure is applied.

Issues with these procedures that are addressed in the remainder of this article include:

The amount of data used to constrain the single-site GMM is small, which does not allow for the evaluation of important effects such as magnitude–scaling relationships that include region-specific break points, near-source saturation effects, region-specific anelastic attenuation effects that distinguish forearc and backarc attenuation, and differences between event types (interface vs intraslab).

The single-site GMM produces different ground motion estimates for the CU site than would be provided by regional models (e.g. Contreras et al., 2025; García-Soto and Jaimes, 2017; Jaimes and García-Soto, 2020; Kuehn et al., 2023), producing an artificial discontinuity in seismic hazard around the margins of Mexico City.

The site transfer function does not use site information (other than location) from the site of interest. As a result, the site information is assumed to be adequately represented by the averaging of transfer functions from neighboring sites. This does not allow for a local deviation of site condition, for example, from a filled stream channel or similar.

Nonlinearity in site response is not considered.

Due to limited documentation, the procedure is not reproducible. Analysts can apply the procedure through the SASID software but cannot independently apply or check it.

We illustrate site responses derived using the reference site approach as provided by the SASID software for 65 sites in the Valley of Mexico at the locations shown in Figure 2 that recorded three large events in 2017–2018. The SASID program requires only location as input. Effects of site period are implicit to the interpolation procedure. Nonlinearity is not considered, so parameters related to the strength of shaking are not required. As shown in Figure 5, response spectra amplification for sites in each of the five zones relative to Zone I follow a consistent pattern, whereby amplification decreases slightly from 0.1 s to about 0.2–0.3 s, increases to a peak that is variable in amplitude and peak period by zone, and then decreases at longer periods. The peak amplification values range, on average, from 1.0 in ln units (∼2.7) for Zone II to approximately 1.9–2.5 in ln units for lakebed zones (∼7–12). The period of the peaks increases with zone number, being about 0.8 s for Zone II and ranging from 1.5 to 3.0 s for Zone III. Within-zone dispersion across sites is modest.

Figure 5.

Response spectra amplification in natural log units for various locations in five Valley of Mexico zones using the methodology in SASID.

Analysis of Valley of Mexico site response

In this section, we first apply a non-reference site approach to estimate site response at ground motion recording sites in the Valley of Mexico. The estimated site responses are then used to develop a sub-regional site response model.

Ground motion analysis

The ground motion data used in this study are free-field acceleration recordings that include both body and surface waves, each of which contributes to the spectral accelerations used to infer site response. The analysis of ground motion data to evaluate site response effects follows a non-reference site approach (adapted from Field and Jacob (1995)). In this approach, residuals are computed between recorded data and a GMM, which are then partitioned to evaluate site effects. By using a GMM to develop the comparison spectra, our application of the non-reference site approach differs from Field and Jacob (1995), who used an inverted Fourier amplitude spectrum, although at a conceptual level, the approaches align. The GMM-based application of the non-reference site approach used here has extensive precedent (e.g. Parker and Stewart, 2022; Seyhan and Stewart, 2014). The data considered in this investigation are an expanded version of the NGA-Sub database—the expansion is to include additional events that had not originally been considered, including three large events in 2017–2018. Information on the expanded database is provided in Contreras et al. (2023a) and in a companion paper (Contreras et al., 2025) and is not repeated here. The GMM that is applied is the Parker et al. (2022) model developed in the NGA-Sub project with adjustments to account for additional anelastic attenuation in the backarc region (Contreras et al., 2023a, 2025).

The total (unpartitioned) residuals are computed as the difference between data and model,

δ_{i j}^{v} = \ln (Y_{ij}) - [μ_{\ln, ij}^{r} (M_{i}, F_{S}, {(R_{rup})}_{ij}) + F_{V, j}],

(1)

where ( $Y_{ij}$ ) is a ground motion intensity measure for event $i$ recorded at site $j$ , $μ_{\ln, i j}^{r}$ is the mean ground motion prediction for reference rock site conditions in natural log units from the adjusted GMM (with the arguments of moment magnitude, $M$ , event-type parameter, $F_{S}$ , and rupture distance $R_{rup}$ ), and $F_{V, j}$ is an ergodic site amplification model (Parker and Stewart, 2022) with linear and nonlinear components,

F_{V} = F_{lin} + F_{nl} .

(2)

The superscript $v$ in Equation 1 indicates that the residuals were computed using a site amplification model. Total residuals are partitioned using mixed effects analyses to identify systematic event terms ( $η_{E, i}$ ), site terms ( $η_{S, j}^{v}$ ), and remaining, essentially random, residual components ( $ε_{ij}$ ),

δ_{ij}^{v} = c + η_{E, i} + η_{S, j}^{v} + ε_{ij} .

(3)

The mixed-effects analyses were performed in MATLAB using the fitlme command (documentation available at https://www.mathworks.com/help/stats/linear-mixed-effects-models.html).

For the investigation of site response nonlinearity in the next section, it is useful to consider within-event residuals relative to a reference rock condition, $δ W_{ij}^{r}$ . These can be computed as,

δ W_{ij}^{r} = δ W_{ij}^{v} + {(F_{lin})}_{j} + {(F_{nl})}_{ij}

(4)

where $F_{lin}$ and $F_{nl}$ are the two components of the ergodic site response model ( $F_{V}$ ); note that while both components depend on the site, the nonlinear component also depends on the event because different events produced different shaking intensities and hence have different levels of nonlinearity.

The addition of the ergodic site amplification model in the $δ W_{ij}^{r}$ calculation (equation 4) offsets the subtraction of this same model in the original residuals calculation (equation 1). This is desirable because (1) the global model may be anticipated to be ineffective in the Valley of Mexico because of its unusual geologic conditions and because applications require extrapolation of the model to V_S₃₀ values below the operable range (<150 m/s) and (2) it provides residuals referenced to a 760 m/s site condition, which is needed to quantify site response effects as used in model development.

Site response nonlinearity

Figure 6 shows example plots of $δ W^{r}$ (representing site amplification) against the strength of shaking on reference rock (denoted $PG A_{r}$ ) for two intensity measures (PGA and PSA(2.0 s)) and sites within Zones I and IIIa. The results in Figure 6 are from three well-recorded earthquakes (EQ1: 2017 M8.3 Tehuantepec, EQ2: 2017 M7.2 Puebla-Morelos, EQ3: 2018 M7.2 Pinotepa Nacional). Two of these events (Tehuantepec and Pinotepa Nacional) produced weak shaking conditions ( $PG A_{r} \approx 0.0015 - 0.0025 g$ ) and one (Puebla-Morelos) produced relatively strong shaking ( $PG A_{r} \approx 0.06 - 0.09 g$ ). The results suggest appreciably lower levels of PGA site response for the Puebla-Morelos event than for the weaker shaking events, with the level of reduction being greater for the softer site condition (Zone IIIa) than for the firmer condition (Zone I). For PSA(2.0 s), the site responses for the three events are largely similar. These features are consistent with expectation when site response is nonlinear, which is known to reduce short-period site response, has little effect for long-period site response, and is of increasing importance as site conditions soften. These site response features, while consistent with many prior observations worldwide (e.g. Parker and Stewart, 2022; Seyhan and Stewart, 2014), are not incorporated into current design procedures (SASID) for the Valley of Mexico, which assume site response to be linear. This study is not the first to identify nonlinear site response in Mexico City; Singh et al. (1988b) and Sánchez-Sesma et al. (1988) identified nonlinear site response during the Michoacán earthquake at specific sites (Central de Abastos Oficina, CDAO and Central de Abastos Frigorífico, CDAF), although other investigators have disagreed (Chávez-García and Bard, 1994) and the topic has remained controversial (Flores-Estrella et al., 2007).

Figure 6.

Reference rock within-event residuals for sites in Valley of Mexico Zones I and IIIa for the intensity measures of PGA and PSA (2.0 s).

The evidence for nonlinearity shown in Figure 6 was derived from response spectral ordinates, which have some dependence on spectral shape. As shown by Stafford et al. (2017), spectral shape effects are important for small magnitudes (generally < 4) for which linear site response in Fourier amplitudes can appear to be nonlinear in PSA due to sensitivities of earthquake corner frequencies to magnitude within the frequency range of oscillator response saturation. However, given the M range of the data considered (7.2–8.3), such effects are very unlikely to explain the nonlinearity identified in this study.

For Zone I (hills), most of the recording stations (15) are located to the southwest of Mexico City (Figures 1 and 2). However, there are five stations in other areas of the city mapped as hills, three of which presented significant differences in site response. One of these stations is located on the slopes of an isolated hill within the lake zone (Cerro del Peñón) and two are located on the slopes of hills located in the north of the city (Cerro Tepeyac and Cerro de Guerrero). Site responses for these three stations are shown in Figure 6 with hollow triangles. The Cerro del Peñón station shows higher within-event residuals than the average for Zone I (possibly influenced by the lakebed), whereas Cerro Tepeyac and Cerro de Guerrero present lower within-event residuals than average for Zone I. Due to this anomalous behavior, the data from these three stations were not used to evaluate the overall nonlinear site response of Zone I.

To model the nonlinear components of site response, we fit the data in Figure 6 using Equation 5,

{(F_{V})}_{z} = f_{1 z} + f_{2 z} \ln (\frac{PG A_{r} + f_{3 z}}{f_{3 z}})

(5)

where f_1,z to f_3,z are coefficients for zone $z$ , and $PG A_{r}$ is as defined previously. Parameter $f_{1 z}$ represents the linear component of the site response for each zone z, whereas the second term in the sum is $F_{nl}$ . These fits were performed setting $f_{3 z} = 0.01 g$ and then regressing for $f_{1 z}$ and $f_{2 z}$ using least-squared procedures. Almost any reasonable value of $f_{3 z}$ could have been selected for the present application. The reason we selected 0.01 g was because it produced values of $f_{2 z}$ that are generally compatible with prior experience from subduction zone site responses (for similar V_S₃₀ values and IMs). This is shown in Figure 7, where we plot the regressed values of $f_{2 z}$ as a function of period for the six zones. Predictions from Parker and Stewart (2022) and smoothed values of the Valley of Mexico coefficients are also shown. The Parker and Stewart (2022) model uses a value of $f_{3 z} = 0.05 g$ , which is constrained from global data generally from firmer site conditions than in the Valley of Mexico. Parameter $f_{3 z}$ represents the level of reference site shaking level where nonlinear effects begin to become appreciable. Given the very soft soil conditions in the Valley of Mexico, larger strains could be expected for modest shaking conditions than for relatively firm sites; as a result, a reduction of $f_{3 z}$ to account for this effect has a physical basis.

Figure 7.

Values of parameter f₂ estimated using the Valley of Mexico data for each zone along with a smoothed fit and model predictions from Parker and Stewart (2022) using the median V_S₃₀ for each zone.

The results in Figure 7 show strongly negative values of $f_{2 z}$ at short periods for all zones. For softer sites (e.g. Zones III), the values are more negative than for relatively firm sites (Zone I) and the negative values extend to longer periods (e.g. at 1.0 s, $f_{2 z} \approx 0$ for Zone I but remains negative for all Zone III groups). A smooth fit to regressed coefficients was provided by eye as shown in the figure. The resulting period-dependent values of $f_{2 z}$ are given in Table 2.

Table 2.

Values of the $f_{2 z}$ parameter for use in nonlinear model (Equation 5).

Intensitymeasure	Oscillator period (s)	VM seismic zone
		I	II	IIIa	IIIb	IIIc	IIId
PGV	–	−0.2297	−0.1764	−0.1812	−0.2664	−0.4255	−0.4499
PGA	–	−0.2659	−0.4117	−0.4918	−0.6046	−0.6928	−0.6607
PSA	0.010	−0.2274	−0.3694	−0.4553	−0.5619	−0.6539	−0.6210
	0.020	−0.2410	−0.3875	−0.4720	−0.5800	−0.6660	−0.6380
	0.025	−0.2585	−0.4015	−0.4864	−0.5952	−0.6845	−0.6562
	0.030	−0.2756	−0.4187	−0.5021	−0.6130	−0.7017	−0.6753
	0.040	−0.2930	−0.4383	−0.5200	−0.6350	−0.7215	−0.7004
	0.050	−0.3100	−0.4600	−0.5370	−0.6597	−0.7456	−0.7202
	0.075	−0.3063	−0.4630	−0.5439	−0.6728	−0.7596	−0.7318
	0.100	−0.2894	−0.4450	−0.5310	−0.6670	−0.7633	−0.7229
	0.150	−0.2013	−0.3855	−0.4613	−0.6170	−0.7310	−0.6857
	0.200	−0.1500	−0.3375	−0.4100	−0.5484	−0.6760	−0.6160
	0.250	−0.1150	−0.3014	−0.3600	−0.5000	−0.6176	−0.5556
	0.300	−0.0896	−0.2700	−0.3283	−0.4600	−0.5511	−0.4878
	0.400	−0.0650	−0.2250	−0.2800	−0.3950	−0.4460	−0.3960
	0.500	−0.0460	−0.1920	−0.2430	−0.3500	−0.3760	−0.3330
	0.750	−0.0200	−0.1280	−0.1800	−0.2600	−0.2760	−0.2330
	1.000	−0.0070	−0.0800	−0.1300	−0.2050	−0.2160	−0.1730
	1.500	0.0000	−0.0250	−0.0700	−0.1300	−0.1420	−0.1000
	2.000	0.0000	0.0000	−0.0320	−0.0750	−0.0975	−0.0659
	2.500	0.0000	0.0000	−0.0100	−0.0350	−0.0720	−0.0500
	3.000	0.0000	0.0000	0.0000	−0.0120	−0.0520	−0.0370
	4.000	0.0000	0.0000	0.0000	0.0000	−0.0320	−0.0200
	5.000	0.0000	0.0000	0.0000	0.0000	−0.0175	−0.0120
	7.500	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
	10.000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000

Linear site responses at recordings sites

To model the linear component of site response, it is necessary to correct within-event residuals for nonlinearity. This adjustment is made using Equation 6,

δ W_{ij}^{r, lin} = δ W_{ij}^{r} - {(F_{nl})}_{ijz},

(6)

where ${(F_{nl})}_{ijz}$ is dependent on event $i$ and site $j$ , because of their effect on $PG A_{r}$ , and on zone $z$ , which determines the coefficients that are applied. As shown by Contreras et al. (2023a), when Equation 6 is applied, the adjusted within-event residuals are effectively independent of $PG A_{r}$ . These within-event residuals are then used in mixed-effects analyses to compute site terms referenced to a site condition of V_S₃₀ = 760 m/s, which are denoted $η_{S}^{r}$ ,

δ W_{ij}^{r, lin} = η_{S, j}^{r} + ε_{ij}

(7)

Equation 7 does not include a constant term, which forces all of the site response bias into the $η_{S}^{r}$ term and $ε_{ij}$ is the remaining residual. The physical meaning of $η_{S}^{r}$ is the mean difference between linearized ground motions for a given site and what would be expected for the site if its site condition was V_S₃₀ = 760 m/s. Subsequently, $η_{S}^{r}$ is referred to as the “linear site response.”

Typical characteristics of Valley of Mexico site responses are illustrated with example results for seven sites from different zones (site locations in Figure 2). The linear site responses for these seven sites are shown in Figure 8. The Zone I site (CUP5; this is the modern station code associated with the traditional CU site) has appreciable amplification, ranging from 0.5 to 1.7 (natural log units), with the maximum amplification being essentially constant for T = 1.5–5s. The Zone II site (AO24) has higher amplification for PGA and PSA, with a peaked response at 0.8 s that is much stronger than for CUP5 (∼2.5), which then decays sharply for T > 1 s. The Zone III sites (IIIa, IB22; IIIb, SCT2; IIIc, XP06; IIId, AE02) also have larger amplification than CUP5 for PGA and PSA. At the longer periods, the site responses are peaked at amplitudes ranging from 3.2 to 3.9 (factors of 25–50 approximately) at periods of 1.5 s for IIIa, 1.5–2 s for IIIb, 2.5 s for IIIc, and 4–5 s for IIId.

Figure 8.

Period-dependent linear site response as derived from non-reference site approach for seven sites in Mexico City (locations in Figure 2).

Results in Figure 8 for the CUP5 site are particularly notable, because that site has been widely used as a reference site in prior studie of Valley of Mexico site response (e.g. Reinoso and Ordaz, 1999). The V_S₃₀ for this site is 295 m/s (Table A1, see ID#7) and its site response is appreciable. Strong site response at the CUP5 site was previously recognized by Singh et al. (1995) and attributed to complex shallow structure. Reference sites used in other regions have generally had much stiffer geologic conditions (e.g. 750–1000 m/s; Borcherdt, 1994, 2002).

Site responses relative to a 760 m/s condition are compared to those derived using SASID in Figure 9. The patterns of site responses between zones and across periods are broadly consistent, although the present site responses are larger by amounts ranging from 0.7 to 1.4 in natural log units (factors of 2.0–4.0) due to the different reference site conditions. These substantial differences are largely a consequence of the significant site response at the CU site. This is illustrated by the “adjusted median” lines in Figure 9, which are computed as the sum of the SASID medians and the CU site response. These adjusted medians approximately convert the SASID site response from the CU reference condition to a 760 m/s reference condition. The results are much closer to the median non-reference site response, but are generally somewhat higher (∼0.1–0.5, or factors of ∼1.1–1.6) for T > 1.0 s and lower for T < 1.0 s.

Figure 9.

Period-dependent linear site response as derived from reference and non-reference site approaches for six zones in the Valley of Mexico (Figure 2). The different levels of site amplification within a zone are related to different reference conditions.

Linear amplification model

Due to the strong correlation of V_S₃₀ with depth and site period T_S (Figure 4), our model development approach begins by modeling amplification trends based solely on V_S₃₀ and subsequently checking dependencies on $z_{LCG}$ and T_S using residuals analysis. Figure 10 shows site response for all Valley of Mexico sites as a function of V_S₃₀ for the intensity measures of PGA, PGV, PSA (0.3 s), and PSA (3.0 s). At short periods (e.g. PGA), the site response increases as V_S₃₀ decreases from about 400 to 150 m/s. For softer sites(< 150 m/s), the site response does not scale with V_S₃₀. This general pattern is retained up to a period of 2.0 s, but for longer periods, the saturation at low V_S₃₀ gradually disappears, instead continuing to increase as V_S₃₀ decreases to its minimum value of 50 m/s. The data trends in Figure 10 are distinct from the Parker and Stewart (2022) ergodic model when that model is extrapolated below 150 m/s. This is notable because that model is generally compatible with site response trends elsewhere in Central America and Mexico (CAM; Contreras et al., 2025).

Figure 10.

Fit of Valley of Mexico model to site amplification data for peak acceleration, peak velocity, PSA (0.3 s), and PSA (3.0 s). Model coefficients are indicated.

To capture the subregional response, the following model is applied:

F_{lin}^{VM} = {\begin{matrix} s_{1 vm} \ln (\frac{V_{S 30}}{V_{1 vm}}) + s_{2 vm} \ln (\frac{V_{1 vm}}{V_{ref}}) + f_{VM} & V_{S 30} \leq V_{1 vm}, \\ s_{2 vm} \ln (\frac{V_{S 30}}{V_{ref}}) + f_{VM} & V_{1 vm} < V_{S 30} < V_{ref}, \end{matrix}

(8)

where the VM superscript indicates it is a Valley of Mexico model, the $vm$ subscripts indicate that the coefficients are specific to that subregion (otherwise the same variable names are used as in the Parker and Stewart (2022) global model), V_ref = 760 m/s, $V_{1 vm}$ is a break velocity, $s_{1 vm}$ is the V_S₃₀-scaling for ${V_{S}}_{30} < V_{1 vm}$ , and $s_{2 vm}$ is the V_S₃₀-scaling for ${V_{S}}_{30} > V_{1 vm}$ . This model is similar to that used by Parker and Stewart (2022), with the following changes: (1) site response is only modeled for V_S₃₀ < V_ref (higher velocities are not applicable in the Valley of Mexico) and (2) an amplification shift parameter f_VM is introduced that allows for higher site response across all V_S₃₀, which is needed across the full period range. Figure 10 shows the fit of the model in Equation 8 to the data, with model coefficients indicated. The selected function fits the data well.

Figure 11 shows the period-variations of model coefficients. Negative values of slope parameters indicate that amplification increases as V_S₃₀ decreases. For short periods (0.01–0.2 s), the slope parameters have relatively little change with period, with $s_{1 vm} ~ - 0.2$ and $s_{2 vm}$ ranging from −1.25 to −0.9. The strongest negative slopes are at long periods (>2.0 s for $s_{1 vm}$ , >0.7 s for $s_{2 vm}$ ), which are needed to model the large long-period amplification. Positive values of $s_{1 vm}$ occur for periods of 0.5–2.0 s, indicating a decrease of amplification as V_S₃₀ decreases in this range. Break velocity $V_{1 vm}$ ranges from 120 m/s at short periods to 70 m/s for T > 2 s. Amplification shift parameter f_VM ranges from −0.5 to 1, presenting negative values at short periods (0.01–0.2 s) and large, positive values at long periods (>2.0 s). The adopted period-dependent values of these parameters are given in Table 3.

Figure 11.

Period dependence of coefficients for Valley of Mexico subregional linear site response model in Equation 8: (a) V_S₃₀-scaling gradients, (b) reference velocity V_ref and corner velocity V_1vm, and(c) amplification shift parameter f_VM.

Table 3.

Values of the parameters of the linear Valley of Mexico site response model.

Intensitymeasure	Oscillatorperiod (s)	V _S ₃₀-scaling gradients		Cornervelocity	Amplificationshift parameter	$δ_{z}$ slope parameters
		$s_{1 vm}$	$s_{2 vm}$	$V_{1 vm}$	$f_{VM}$	$d_{1}$	$d_{2}$
PGV	–	−0.9889	−0.7682	120	0.6105	−0.0251	0.0036
PGA	–	−0.1911	−1.2318	120	−0.2107	−0.0224	−0.0077
PSA	0.010	−0.1902	−1.2467	120	−0.3080	−0.0220	−0.0084
	0.020	−0.1859	−1.2438	120	−0.3248	−0.0220	−0.0079
	0.025	−0.1965	−1.2258	120	−0.3300	−0.0220	−0.0077
	0.030	−0.1999	−1.2159	120	−0.3388	−0.0220	−0.0075
	0.040	−0.2077	−1.1935	120	−0.3586	−0.0220	−0.0068
	0.050	−0.2066	−1.1701	120	−0.3876	−0.0223	−0.0064
	0.075	−0.2075	−1.1237	120	−0.4202	−0.0235	−0.0059
	0.100	−0.2090	−1.0406	120	−0.3593	−0.0250	−0.0057
	0.150	−0.1884	−0.9579	120	−0.2206	−0.0260	−0.0058
	0.200	−0.1711	−0.8762	120	−0.0433	−0.0250	−0.0040
	0.250	−0.1450	−0.8178	120	0.1182	−0.0205	−0.0060
	0.300	−0.1393	−0.8388	120	0.1500	−0.0175	−0.0100
	0.400	0.0138	−0.9225	120	0.2458	−0.0165	−0.0117
	0.500	0.1153	−0.9932	120	0.3407	−0.0160	−0.0120
	0.750	0.1845	−1.1995	120	0.2776	−0.0150	−0.0115
	1.000	0.3848	−1.4444	120	0.0801	−0.0135	−0.0095
	1.500	1.0491	−1.5348	90	0.0586	−0.0120	−0.0075
	2.000	0.0364	−1.3148	80	0.1575	−0.0200	−0.0308
	2.500	−1.8641	−0.9172	80	0.5922	−0.0439	−0.0489
	3.000	−3.2350	−0.7303	70	0.7914	−0.0300	−0.0350
	4.000	−4.4205	−0.4244	70	1.0622	−0.0010	0.0132
	5.000	−3.7641	−0.3361	70	1.0711	0.0000	0.0340
	7.500	−1.8551	−0.2423	70	1.0375	0.0000	0.0367
	10.000	−1.2880	−0.2272	70	1.0242	0.0000	0.0367

Residuals analyses

In this subsection, we evaluate the performance of the V_S₃₀-scaling model proposed for the linear site response of the Valley of Mexico, extend the model to consider effects of sediment depth based on trends observed in residuals analyses, check the sensitivity of residuals to site period T_S, and develop models for within-event components of aleatory variability.

To evaluate model performance, we take the linearized within-event rock residuals from Equation 6 as “data” and subtract the proposed model prediction (Equation 8) to evaluate within-event residual $δ W^{v, lin}$ , where the “ $v$ ” superscript indicates that V_S30-dependent site response is considered:

δ W_{ij}^{v, lin} = δ W_{ij}^{r, lin} - {(F_{lin}^{VM})}_{j} .

(9)

Figure 12 shows that the model removes the principle trends of within-event residuals with V_S₃₀ for PGA, PGV and PSA at 0.3, 1.0, 2.0, and 3.0 s. This is evident from the binned means of these residuals generally being nearly zero and not exhibiting a trend with V_S₃₀. Some non-zero binned means ( $< \pm 0.5$ ) occur at periods between 1.5 and 3.0 s (e.g. the sites in Zone IIIc with V_S₃₀∼57 m/s have a positive binned mean for PSA at 3.0 s). This is caused by these sites having a natural site period of T_S∼2.5–3.5 s that approximately matches the oscillator period. Figure 12 also shows that Zone I (V_S₃₀∼300 m/s) sites have larger dispersion than other zones (although the binned means are centered at zero). This indicates a relatively large variability of site response in the Hills Zone, likely reflecting heterogeneity of geotechnical and topographic conditions.

Figure 12.

Trend of linearized within-event residuals with V_S₃₀ for sites in the Valley of Mexico for PGA, PGV, and PSA at T = 0.3, 1.0, 2.0, and 3.0 s. Error bars are 95% confidence intervals.

Sediment depth has been found for multiple regions to have predictive capacity for site response beyond that provided by V_S₃₀ (e.g. Day et al., 2008; Nweke et al., 2022). Figure 4b shows that depth to the Lower Coarse Grain layer is negatively correlated with V_S₃₀, which makes it poorly suited as an additional predictive parameter. To avoid this problem, we define a model for mean depth conditioned on V_S₃₀ ( ${\bar{z}}_{LCG}$ ) as shown in Figure 4b and use depth differential ( $δ z = z_{LCG} - {\bar{z}}_{LCG}$ ) as a parameter against which to evaluate model performance. Figure 13a and b shows trends of residuals ( $δ W_{ij}^{v, lin}$ ) with $δ z$ , which are significant for periods T > 1s for $δ z < - 5$ m and $δ z > 5$ m (means of residuals are zero for $δ z = - 5 to 5 m$ ). We fit the data using Equation 10,

F_{z} = {\begin{matrix} d_{1} (δ z + 5) δ z < - 5 m \\ 0 δ z = - 5 to 5 m \\ d_{2} (δ z - 5) δ z > 5 m \end{matrix}

(10)

Figure 13.

Trend of linearized within-event residuals with differential depth $δ z$ for sites in the Valley of Mexico for (a) PSA at T = 2.5 s and (b) PSA at T = 5 s. (c) Variation of slope coefficients with period T.

Figure 13c shows the trends with period of $d_{1}$ and $d_{2}$ , which we smooth as shown in the figure.

Although V_S₃₀ is used as the primary site parameter in the proposed model (equation 8), the site parameter traditionally associated with the different Valley of Mexico zones is site period, $T_{S}$ (Lermo et al., 2020). Accordingly, we investigate whether the proposed model, including the depth correction in Equation 10, is able to capture the known period dependency of Valley of Mexico site response. This is evaluated by computing linearized within-event residuals from Equation 11

δ W_{ij}^{vz, lin} = δ W_{ij}^{r, lin} - [{(F_{lin}^{VM})}_{j} + F_{z, j}] .

(11)

Residuals from Equation 11 are plotted against $T_{S}$ in Figure 14 (plots without adjustment for depth effects are provided by Contreras et al. (2023a)). The results for PGA, PGV, and PSA for T≤ 1.5 s show no trend with $T_{S}$ . This likely occurs because the intensity measures shown either do not correspond with a particular oscillator period (PGV) or their period is lower than those for Valley of Mexico sites, even in the case of Zone I.

Figure 14.

Trend of linearized within-event residuals with T_S for sites in the Valley of Mexico for PGA, PGV, and PSA at T = 0.3, 1.0, 2.0, and 3.0 s. Error bars are 95% confidence intervals.

At longer periods, the residuals patterns in Figure 14 illustrate site resonance effects that are not captured by the V_S₃₀-scaling model. Before discussing the results, it is useful to recall the median values of site period in the different zones (e.g. ∼1 s for Zone II and IIIa, ∼2 s for IIIb, and ∼3 s for IIIc). Residuals plots for a given oscillator period often have peaks when that period aligns with the zone median site period; for example, a peak occurs at $T_{S}$ = 2 s in the plots for PSA (2 s) for the Zone IIIb sites, and a peak occurs at $T_{S}$ = 3 s in the plots for PSA (3 s) for the Zone IIIc sites. These peaks indicate that there are resonance features in the site response that are not captured by the V_S₃₀-scaling model. We considered supplementing the V_S₃₀-scaling model in Equation 8 with a pulse model that would introduce peaks at a site period (e.g. as done by Buckreis et al. (2024) for a different region). While such a feature would improve data fit, the level of improvement would be marginal. Consider, for example, Zone IIIb—the overall level of amplification at the 2 s peak is approximately 3.5 in natural log units (Figure 9), whereas the height of the binned means in Figure 14 is approximately 0.25–0.3. Accordingly, the overwhelming majority of the site response is captured by the V_S₃₀-scaling model and a resonance model would introduce a relatively minor modification (as expected given the parameter correlation, Figure 4). Moreover, the development of a model conditioned on site period, which is derived from mHVSR, is impractical at the present time due to the limited availability of mHVSR curves (23 of 89 sites), as indicated previously. For these reasons, we decided against adding this component.

Residuals analyses using the Valley of Mexico data along with the recommended model (i.e. reference site regional GMM from Contreras et al. (2025), which is modified from Parker et al. (2022); linear site response from Equation 8; nonlinear site response from Equation 5) provides within-event residuals that can be partitioned using mixed-effects analysis into site terms ( $η_{S, j}^{v}$ ) and remaining residuals ( $ε_{ij}$ ). This final set of residuals includes updated estimates of event terms that account for the effects of the Valley of Mexico site response model. Mixed effects analyses also provide corresponding standard deviation terms, which are the site-to-site standard deviation ( $ϕ_{S 2 S}$ ) and single-station within-event standard deviation ( $ϕ_{SS}$ ). These are important components of the overall aleatory variability model. Their values could be expected to differ from values from global models as a result of the specific location and geology of the Valley of Mexico.

Figure 15 shows subregional $ϕ_{S 2 S}$ and $ϕ_{SS}$ terms derived for Valley of Mexico sites for the sampled oscillator periods along with a smoothed model. Also shown for comparison purposes are the corresponding standard deviation terms from the Parker et al. (2022) global model exercised for conditions generally appropriate for the Valley of Mexico data (V_S30 = 150 m/s, R_rup = 100 km). The subregional dispersion terms are lower than those from the global model, with the differences being particularly large for $ϕ_{S 2 S}$ . These lower dispersions are expected because the variations in path and site response in the Valley of Mexico subregion are clearly less than those present globally.

Figure 15.

Site-to-site and single-station dispersion terms derived for Valley of Mexico sites. Larger values from global model of Parker et al. (2022) shown for comparison.

Validation of the site response model

In this subsection, we validate the site response model for the Valley of Mexico presented in the previous section. The validation uses a data set from the 1985 M8.0 Michoacán earthquake. Ground motion data from 26 stations that recorded this event are included in the NGA-Sub database, 10 of which are located in Mexico City, which is known to have experienced substantial site response (Ordaz et al., 1988; Seed et al., 1987; Singh et al., 1988a, 1988b). However, the 10 Mexico City recordings were not considered in the development of the site response model because of data screening criteria that were applied (matching those from Parker et al. (2022)). The criteria that led to the Mexico City sites being screened mainly relate to maximum usable distance (i.e. the maximum distance for which recordings can be considered without introducing sampling bias, which is denoted $R_{\max}$ ) and lack of necessary site metadata (V_S₃₀). The $R_{\max}$ value for this event was 137 km and the Mexico City sites were at distances of 233–246 km. The V_S₃₀ screen no longer applies because in this study V_S₃₀ values have been provided for all stations using the procedures given in the Site Parameters section.

While the elimination of the Mexico City data was necessary in order to maintain consistent data selection criteria, the strong subregional site response likely causes the amplitudes to exceed noise levels to a sufficient degree that the recordings are representative of the ground shaking that was experienced. Instead of making a “one-time” exception to the screening criteria, we chose to use the data as an independent testing data set for model validation. Because data from the Michoacán event at closer distances (not in the Valley of Mexico) were used by Contreras et al. (2023a, 2025) in the development of the Parker et al. (2022) GMM adjustments, event terms ( $η_{E, i}$ ) for the modified model are available (Figure 16). The events terms are nearly 0 up to 1.5 s and then increase, indicating that the adjusted GMM is unbiased for T < 1.5 s.

Figure 16.

Period dependence of event terms for the 1985 M8.0 Michoacán earthquake overlaid on the mean $\pm τ_{\ln}$ of all CAM events (from Contreras et al. (2025)).

The validation is performed by summing (in natural log units) the adjusted reference rock GMM, the event terms from Figure 16, and the site amplification models (V_S₃₀-scaling only and combined model that includes differential depth effects). $PG A_{r}$ values for Valley of Mexico sites range from 0.013 to 0.015 g (using the event-adjusted Parker et al. (2022) model), which is large enough for modest effects of nonlinearity to appear. The results of these calculations are compared to data in Figure 17 for six sites (of the 10 available) from all of the zones except IIIa, for which no data are available. To examine the impact of the nonlinear model, predicted spectra are shown in Figure 17 both for the linear models and the linear models adjusted by the nonlinear model (Equation 5). The results shown for the V_S₃₀-scaling model applied with the nonlinear model are bracketed by bands of $\pm ϕ_{SS}$ . The depth adjustments only affect the SCT and Tlahuac Bombas sites because the others have $| δ z | < 5 m$ . Overall, the proposed site response model performs well. The misfit of the model to the data is reduced in 60% of cases by the use of the nonlinear model. This exercise provides confidence that the proposed model is consistent with observations from the landmark 1985 event. SASID results are not provided in Figure 17 because it is not configured to provide results for scenario events (it provides hazard results across multiple events).

Figure 17.

Comparison of observed PSA at six stations that recorded the 1985 event (red lines) with predicted PSA for engineering-rock conditions, and after applying the Valley of Mexico site response model (full and linear-only).

Summary and discussion

This research examined site response within the Valley of Mexico, which contains Mexico City, using a non-reference site approach. The work presented here on site response is part of a broader study of ground motions from subduction earthquakes in CAM. The extension of the NGA-Sub project database (Mazzoni et al., 2022) to include additional events and stations in Mexico is presented in a companion paper (Contreras et al., 2025). The supplemented database provides multiple ground motion recordings at sites across the Valley of Mexico that enabled the present application of non-reference site analyses to measure site response and its uncertainty.

We have summarized the unique geotechnical conditions in the Valley of Mexico and compiled available site parameters (V_S₃₀, T_S, $z_{LCG}$ ) from site-specific investigations and other literature. These data support the zonation of the lakebed and surrounding regions (Figures 1 and 2) as currently applied in forward applications of site response using a reference site approach (Ordaz et al., 2024; Reinoso and Ordaz, 1999). This approach relies on normalization of ground motion Fourier amplitudes within the lakebed relative to those at a hill zone site (CU).

We have shown that the CU site, commonly used as a reference, has an appreciable level of site response across a broad period range (based on data from the CUP5 station). The site response results presented in this article are broadly similar to previous findings with regard to relative amplification levels between lakebed zones and the period ranges where these amplifications are maximized. However, there are two notable aspects of our findings: (1) amplification levels are higher due to the referencing of amplification to V_ref = 760 m/s in this work, and (2) by examining data from multiple events, we find short period site amplifications exhibit nonlinearity. The use of the non-reference site approach for the derivation of amplifications is particularly significant from a ground motion modeling perspective, because it removes the need to derive single-station reference GMMs for the CU site as has been done in prior work (Arroyo et al., 2024; Jaimes et al., 2015; Jaimes et al., 2006; Ordaz et al., 1994). Doing so substantially increases the level of rigor with which the reference site GMM can be defined.

A nonlinear model conditioned on V_S₃₀ is presented for site response in the Valley of Mexico. The linear portion of the model is able to capture the main trends of site response variations across the lakebed. The linear model includes the option of applying a differential depth parameter derived from $z_{LCG}$ (this option can be turned off by setting δz = 0). While some resonance effects are seen in the data that are not captured, these effects are small compared to what is captured from the V_S₃₀-based model and are not included herein. We acknowledge that models conditioned solely on mHVSR parameters, which were not developed in this research, could be similarly effective as the proposed V_S₃₀-based model. However, a challenge related to developing such models is that mHVSR curves are only available for 23 of the 89 sites considered in this research, and 19 of those 23 sites are in only two zones (II and IIIa). The nonlinear portion of the proposed model is admittedly based upon limited data, but those data demonstrate nonlinearity unambiguously. The development of a more rigorous nonlinear model would require calibrated simulations in which the nonlinear soil behavior is considered, which was beyond the scope of this investigation.

The proposed modeling approach for sites in the Valley of Mexico involves the use of a regionally customized GMM and a subregional site amplification model. The GMM is an adaptation of the Parker et al. (2022) model for CAM, with the modifications being solely related to backarc attenuation effects. The site amplification model has linear and nonlinear components (Equation 2). The V_S₃₀ portion of the linear component is given by Equation 8 with the coefficients in Table 3. The optional modification based on δz is given by Equation 10 with coefficients in Table 3. The nonlinear component is given by Equation 5 with the coefficients in Table 2. Corresponding values for the two components of within-event variability ( $ϕ_{S 2 S}$ and $ϕ_{SS}$ ) are given in Figure 15. Coded versions of the subregional site response model presented in this article and the regionally customized GMM presented by Contreras et al. (2025) are provided on the NGA-Subduction project web portal (https://www.risksciences.ucla.edu/nhr3/nga-subduction).

The findings and recommendations presented here provide an opportunity to improve ground motion hazard analyses in Mexico due to the following factors:

Ground motion estimation in the Valley of Mexico (Mexico City) and elsewhere in Mexico is developed using consistent approaches.

The models are transparent, thus allowing for future research where individual components can more readily be interrogated and improved.

Site response is conditioned upon the site-specific V_S₃₀ parameter—this encourages site characterization and avoids ambiguities that can arise for sites near zone boundaries.

Nonlinearity in site response is considered.

The next steps are to evaluate the impact of these models on seismic hazards across Mexico and see whether consensus can be developed regarding their introduction into national and regional seismic design standards. Such analyses have not been undertaken in the present work.

Supplemental Material

sj-docx-1-eqs-10.1177_87552930251316816 – Supplemental material for Valley of Mexico site response from non-reference site approach

Supplemental material, sj-docx-1-eqs-10.1177_87552930251316816 for Valley of Mexico site response from non-reference site approach by Victor Contreras, Jonathan P Stewart, Juan M Mayoral and Xyoli Pérez-Campos in Earthquake Spectra

Supplemental Material

sj-docx-2-eqs-10.1177_87552930251316816 – Supplemental material for Valley of Mexico site response from non-reference site approach

Supplemental material, sj-docx-2-eqs-10.1177_87552930251316816 for Valley of Mexico site response from non-reference site approach by Victor Contreras, Jonathan P Stewart, Juan M Mayoral and Xyoli Pérez-Campos in Earthquake Spectra

Footnotes

Acknowledgements

The authors are grateful to the NGA-Sub project for providing the database and models from which this project was developed. SSN data was obtained by the Servicio Sismológico Nacional (México); station maintenance, data acquisition, and distribution are possible thanks to its personnel. The views expressed herein are those of the authors and do not necessarily reflect the views of the CTBTO Preparatory Commission. Daniel De La Rosa assisted with data collection for sites in Mexico, particularly geologic maps. We appreciate the helpful feedback from the Associate Editor and two anonymous reviewers, which improved the paper.

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: Financial support for the first author was provided by the Chilean National Agency for Research and Development (ANID) through the grant ANID BECAS CHILE-DOCTORADO EN EL EXTRANJERO 72180625 and the Civil and Environmental Engineering Department at UCLA.

ORCID iDs

Victor Contreras

Jonathan P Stewart

Data resources

Raw ground-motion recordings from the RAII-UNAM network, managed by the Institute of Engineering at UNAM, was accessed through the RAII-UNAM network website: http://aplicaciones.iingen.unam.mx/AcelerogramasRSM/Inicio.aspx (last accessed April 2022). Raw ground motions from the RACM network, managed by CIRES (Instrumentation and Seismic Recording Center), were requested from the CIRES website: http://www.cires.org.mx/racm_historico_es.php (last accessed April 2022). Our requests made in September 2017 and April 2020 were granted. Raw data from the IG (Pérez-Campos et al., 2018; SSN, 2020), VM (Quintanar et al., 2018), and UV (Córdoba-Montiel et al., 2018) networks, which are distributed by the National Seismological Service of Mexico (SSN), are accessible through a client called SSNstp upon request (). Our request made on 7 June 2020 was granted.

Supplemental material

Supplemental material for this article is available online.

References

Ahdi

Kwak

Ancheta

Contreras

Kishida

Kwok

Mazzoni

Ruz

Stewart

(2022) Site parameters applied in NGA-Sub database. Earthquake Spectra 38(1): 494–520.

Alcocer

Williams

(1996) Historical and recent changes in Lake Texcoco, a saline lake in Mexico. International Journal of Salt Lake Research 5: 45–61.

Arce

Layer

Macías

Morales-Casique

García-Palomo

Jiménez-Domínguez

Benowitz

Vásquez-Serrano

(2019) Geology and stratigraphy of the Mexico Basin (Mexico City), central Trans-Mexican volcanic belt. Journal of Maps 15(2): 320–332.

Arroyo

Ordaz

Singh

(2024) Prediction of Fourier amplitude spectrum of ground motion in Mexico City from subduction thrust earthquakes. Geofísica Internacional 62(2): 783–801.

Bonilla

Steidl

Gariel

Archuleta

(2002) Borehole response studies at the Garner Valley downhole array, Southern California. Bulletin of the Seismological Society of America 92: 3165–3179.

Bonilla

Steidl

Lindley

Tumarkin

Archuleta

(1997) Site amplification in the San Fernando Valley, California: Variability of site-effect estimation using the S-wave, coda, and H/V methods. Bulletin of the Seismological Society of America 87: 710–730.

Borcherdt

(1970) Effects of local geology on ground motion near San Francisco Bay. Bulletin of the Seismological Society of America 60: 29–61.

Borcherdt

(1994) Estimates of site-dependent response spectra for design (methodology and justification). Earthquake Spectra 10: 617–653.

Borcherdt

(2002) Empirical evidence for acceleration-dependent amplification factors. Bulletin of the Seismological Society of America 92: 761–782.

10.

Buckreis

Stewart

Brandenberg

Wang

(2024) Small-strain site response of soft soils in the Sacramento-San Joaquin Delta region of California conditioned on VS30 and mHVSR. Earthquake Spectra 40(1): 230–260.

11.

Cadet

Bard

P-Y

Duval

A-M

Bertrand

(2012) Site effect assessment using KiK-net data: Part 2—site amplification prediction equation based on f0 and Vsz. Bulletin of Earthquake Engineering 10: 451–489.

12.

Cárdenas-Soto

Ramos-Saldaña

Vidal-García

(2016) Interferometría de ruido sísmico para la caracterización de la estructura de velocidad 3D de un talud en la 3ª Sección del Bosque de Chapultepec, Ciudad de México. Boletín De La Sociedad Geológica Mexicana 68(2): 173–186 (in Spanish).

13.

Chamorro

(2016) Medición de los valores de velocidad de onda cortante (Vs) en suelos provenientes del ex lago Texcoco. MS Thesis, Institute of Engineering, National Autonomous University of Mexico (UNAM), Mexico (in Spanish).

14.

Chávez-García

Bard

(1994) Site effects in Mexico City eight years after the September 1985 Michoacan earthquakes. Soil Dynamics and Earthquake Engineering 13: 229–247.

15.

Contreras

(2022) Characteristics of subduction zone ground motions with an emphasis on Latin America. PhD Dissertation, Civil & Environmental Engineering, University of California, Los Angeles, CA.

16.

Contreras

Stewart

DeLaRosa

Mayoral

Pérez-Campos

(2023a) Regional Customization of Global Subduction Ground Motion Models for Mexico. ReportGIRS-2023-09, Natural Hazards Risk and Resiliency Research Center, B. John Garrick Institute for the Risk Sciences, University of California, LosAngeles, CA.

17.

Contreras

Stewart

De La Rosa

Mayoral

Pérez-Campos

(2023b) Empirical site response of Mexico city developed from regional customization of global subduction ground motion models. In: XIII Congreso Chileno De Sismología E Ingeniería Sísmica ACHISINA 2023, Viña Del Mar, Chile, 24–26 October.

18.

Contreras

Stewart

Pérez-Campos

Mayoral

De La Rosa

(2024) Empirical site response of Mexico city through regionalization of global subduction GMMs. In: 18th World Conference on Earthquake Engineering, Milan, 30 June–5 July.

19.

Contreras

Stewart

Mayoral

Pérez-Campos

(2025) Regionalization of global subduction ground motion model for application in Central America and Mexico. Earthquake Spectra, https://doi.org/10.1177/87552930251322789

20.

Córdoba-Montiel

Iglesias

Pérez-Campos

Sieron

Grupo de Trabajo del Servicio Sismológico Nacional, Valdés-González

Singh

Pacheco

(2018) The broadband seismological network of Veracruz, Mexico: Toward a regional seismotectonic interpretation. Seismological Research Letters 89(2A): 345–355.

21.

Day

Graves

Bielak

Dreger

Larsen

Olsen

Pitarka

Ramirez-Guzman

(2008) Model for basin effects on long-period response spectra in southern California. Earthquake Spectra 24(1): 257–277.

22.

De La Rosa

Flores

Guasch

(2013) Caracterización geo-sísmica de una zona del ex-lago Texcoco para el diseño de un túnel de desagüe. In: XIX National Conference on Earthquake Engineering, Boca Del Río, Mexico, 12–15 June (in Spanish).

23.

Díaz-Rodríguez

(2003) Characterization and engineering properties of Mexico City lacustrine soils. In: Tan

Phoon

Hight

, et al. (eds) Characterisation and Engineering Properties of Natural Soils. Lisse: Swets & Zeitlinger, pp. 593–643.

24.

Field

Jacob

(1993) The theoretical response of sedimentary layers to ambient seismic noise. Geophysical Research Letters 20: 2925–2928.

25.

Field

Jacob

(1995) A comparison and test of various site response estimation techniques, including three that are not reference site dependent. Bulletin of the Seismological Society of America 85(4): 1127–1143.

26.

Flores-Estrella

Yussim

Lomnitz

(2007) Seismic response of the Mexico City Basin: A review of twenty years of research. Natural Hazards 40: 357–372.

27.

Flores-Guzmán

Ovando-Shelley

Valle-Molina

(2014) Small-strain dynamic characterization of clayey soil from the Texcoco Lake, Mexico. Soil Dynamics and Earthquake Engineering 63: 1–7.

28.

García-Soto

Jaimes

(2017) Ground-motion prediction model for vertical response spectra from Mexican interplate earthquakes. Bulletin of the Seismological Society of America 107(2): 887–900.

29.

Jaimes

García-Soto

(2020) Updated ground motion prediction model for Mexican intermediate-depth intraslab earthquakes including V/H ratios. Earthquake Spectra 36(3): 1298–1330.

30.

Jaimes

Ramirez-Gaytan

Reinoso

(2015) Ground-motion prediction model from intermediate-depth intraslab earthquakes at the hill and lake-bed zones of Mexico city. Journal of Earthquake Engineering 19(8): 1260–1278.

31.

Jaimes

Reinoso

Ordaz

(2006) Comparison of methods to predict response spectra at instrumented sites given the magnitude and distance of an earthquake. Journal of Earthquake Engineering 10(6): 887–902.

32.

Juárez-Camarena

Auvinet-Guichard

Méndez-Sánchez

(2016) Geotechnical zoning of Mexico Valley subsoil. Ingeniería Investigación Y Tecnología 17(3): 297–308.

33.

Kuehn

Bozorgnia

Campbell

Gregor

(2023) A regionalized partially nonergodic ground-motion model for subduction earthquakes using the NGA-Sub database. Earthquake Spectra 39(3): 1625–1657.

34.

Lermo

Chávez-García

(1993) Site effect evaluation using spectral ratios with only one station. Bulletin of the Seismological Society of America 83(5): 1574–1594.

35.

Lermo

Chávez-García

(1994a) Site effect evaluation at Mexico City: Dominant period and relative amplification from strong motion and microtremor records. Soil Dynamics and Earthquake Engineering 13(6): 413–423.

36.

Lermo

Chávez-García

(1994b) Are microtremors useful in site response evaluation? Bulletin of the Seismological Society of America 84(5): 1350–1364.

37.

Lermo

Rodríguez

Singh

(1988) The Mexico earthquake of September 19, 1985—Natural period of sites in the valley of Mexico from microtremor measurements and strong motion data. Earthquake Spectra 4(4): 805–814.

38.

Lermo

Sánchez-Sesma

Ramos

Álvarez

Jiménez

Torres

Machado

(2020) Actualización de la zonificación sísmica de la Ciudad de México y áreas aledañas-Parte Norte, Informe Técnico No. De Convenio ISCDF/CEC-04/2020-13 (in Spanish).

39.

Luna

(2012) Análisis de movimiento sísmico en la porción sur de la ciudad de México, una evaluación de la amenaza sismológica. MS Thesis, School of Engineering and Architecture, National Polytechnic Institute of Mexico, Mexico (in Spanish).

40.

Marsal

Mazari

(1969) El subsuelo de la Ciudad de México. Reg. 63, Facultad De Ingeniería, UNAM, Mexico, p. 614.

41.

Mayoral

Asimaki

Tepalcapa

Wood

Roman-de la Sancha

Hutchinson

Franke

Montalva

(2019a) Site effects in Mexico City basin: Past and present. Soil Dynamics and Earthquake Engineering 121: 369–382.

42.

Mayoral

Castañon

Alcantara

Tepalcapa

(2016) Seismic response characterization of high plasticity clays. Soil Dynamics and Earthquake Engineering 84: 174–189.

43.

Mayoral

Castañon

Albarran

(2017) Regional subsidence effects on seismic soil-structure interaction in soft clay. Soil Dynamics and Earthquake Engineering 103: 123–140.

44.

Mayoral

De la Rosa

Tepalcapa

(2019b) Topographic effects during the September 19, 2017 Mexico city earthquake. Soil Dynamics and Earthquake Engineering 125: 105732.

45.

Mayoral

Flores

Romo

(2011) Seismic response evaluation of an urban overpass. Earthquake Engineering & Structural Dynamics 40: 827–845.

46.

Mazzoni

Kishida

Stewart

Contreras

Darragh

Ancheta

Chiou

BSJ

Silva

Bozorgnia

(2022) NGA-Sub relational database. Earthquake Spectra 38(2): 1529–1548.

47.

Nogoshi

Igarashi

(1970) On the amplitude characteristics of ambient noise (part 1). Journal of Atmospheric and Oceanic Technology 24: 282–308.

48.

Nogoshi

Igarashi

(1971) On the amplitude characteristics of ambient noise (part 2). Journal of Atmospheric and Oceanic Technology 24: 26–40.

49.

Normas Técnicas Complementarias para el diseño por Sismos (NTCS) (2004) Gaceta oficial del Distrito Federal (Ed.), México, Tomo II, No. 103-BIS, octubre, pp. 55–77 (in Spanish).

50.

Nweke

Stewart

Wang

Brandenberg

(2022) Site response of sedimentary basins and other geomorphic provinces in southern California. Earthquake Spectra 38(4): 2341–2370.

51.

Ordaz

Arroyo

Singh

Salgado-Gálvez

(2024) A PSHA for Mexico City based solely in Fourier-based GMM of the response spectra. Soil Dynamics and Earthquake Engineering 187: 109025.

52.

Ordaz

Singh

Arciniega

(1994) Bayesian attenuation regressions: An application to Mexico city. Geophysical Journal International 117(2): 335–344.

53.

Ordaz

Singh

Reinoso

Lermo

Espinosa

Domínguez

(1988) The Mexico Earthquake of September 19, 1985—Estimation of Response Spectra in the Lake Bed Zone of the Valley of Mexico. Earthquake Spectra 4(4): 815–834.

54.

Parker

Stewart

(2022) Ergodic site response model for subduction zone regions. Earthquake Spectra 38(2): 841–864.

55.

Parker

Stewart

Boore

Atkinson

Hassani

(2022) NGA-subduction global ground motion models with regional adjustment factors. Earthquake Spectra 38(1): 456–493.

56.

Pérez-Campos

Espíndola

Pérez

Estrada

Monroy

Bello

González-López

González Ávila

Contreras Ruiz Esparza

Maldonado

Tan

Rodríguez Rasilla

Vela Rosas

Cruz

Cárdenas

Navarro Estrada

Hurtado

Mendoza Carvajal

Montoya-Quintanar

Pérez-Velázquez

(2018) The Mexican National Seismological Service: An overview. Seismological Research Letters 89(2A): 318–323.

57.

Pérez

(2017) Estudio geofísico de refracción sísmica para la caracterización del terreno, en el conjunto urbano Tlalpan, ubicado en la delegación Coyoacán. CDMX. Technical report, AXDSCA Consultants.

58.

Quintanar

Cárdenas-Ramírez

Bello-Segura

Espíndola

Pérez-Santana

Cárdenas-Monroy

Carmona-Gallegos

Rodríguez-Rasilla

(2018) A seismic network for the Valley of Mexico: Present status and perspectives. Seismological Research Letters 89(2A): 356–362.

59.

Ramos-Martínez

Chávez-Garcia

Romero-Jiménez

Rodríguez-Zúñiga

Gómez-González

(1997) Site effects in Mexico City: Constraints from surface wave inversion of shallow refraction data. Journal of Applied Geophysics 36: 157–165.

60.

Reinoso

Ordaz

(1999) Spectral ratios for Mexico City from free-field recordings. Earthquake Spectra 15(2): 273–295.

61.

Romo

García

(2003) Neurofuzzy mapping of CPT values into soil dynamic properties. Soil Dynamics and Earthquake Engineering 23: 473–482.

62.

Romo

Jaime

Reséndiz

(1988) The Mexico earthquake of September 19, 1985—General soil conditions and clay properties in the Valley of Mexico. Earthquake Spectra 4(4): 731–752.

63.

Sánchez-Sesma

Chávez-Pérez

Suárez

Bravo

Pérez-Rocha

(1988) The Mexico Earthquake of September 19, 1985, on the seismic response of the valley of Mexico. Earthquake Spectra 4: 569–589.

64.

Satoh

Kawase

Matsushima

(2001) Differences between site characteristics obtained from microtremors, S-waves, P-waves, and codas. Bulletin of the Seismological Society of America 91: 313–334.

65.

Seed

Romo

Sun

Jaime

Lysmer

(1987) Relationships between soil conditions and earthquake ground motions in Mexico City in the earthquake of September 19, 1985, Report No. UCB/EERC-87/15, Earthquake Engineering Research Center, UC Berkeley.

66.

Seyhan

Stewart

(2014) Semi-empirical nonlinear site amplification from NGA-West2 data and simulations. Earthquake Spectra 30(3): 1241–1256.

67.

Singh

Lermo

Dominguez

Ordaz

Espinoza

Mena

Quaas

(1988b) A study of amplification of seismic waves in the Valley of Mexico with respect to a hill zone site (CU). The Mexico Earthquake of September 19, 1985. Earthquake Spectra 4: 653–673.

68.

Singh

Mena

Castro

(1988a) Some aspects of the source characteristics and ground motion amplification in and near Mexico City from acceleration data of the September, 1985, Michoacan, Mexico Earthquakes. Bulletin of the Seismological Society of America 78: 451–477.

69.

Singh

Quass

Ordaz

Mooser

Almora

Torres

Vásquez

(1995) Is there truly a “hard” rock site in the Valley of Mexico? Geophysical Research Letters 22(4): 481–484.

70.

Singh

Reinoso

Arroyo

Ordaz

Cruz-Atienza

Pérez-Campos

Iglesias

Hjörleifsdóttir

(2018) Deadly Intraslab Mexico earthquake of 19 September 2017 (Mw 7.1): Ground motion and damage pattern in Mexico City. Seismological Research Letters 89(6): 2193–2203.

71.

SSN (2020) Servicio Sismológico Nacional. Instituto De Geofísica, Universidad Nacional Autónoma De México, México. Available at: https://doi.org/10.21766/SSNMX/SN/MX. Accessed February 2025.

72.

Stafford

Rodriguez-Marek

Edwards

Kruiver

Bommer

(2017) Scenario dependence of linear site-effect factors for short-period response spectral ordinates. Bulletin of the Seismological Society of America 107(6): 2859–2872.

73.

Stephenson

Lomnitz

(2005) Shear-wave velocity profile at the Texcoco strong-motion array site, Valley of Mexico. Geofísica Internacional 44(1): 3–10.

74.

Theodulidis

Bard

Archuleta

Bouchon

(1996) Horizontal-to-vertical spectral ratio and geological conditions: The case of Garner Valley downhole array in southern California. Bulletin of the Seismological Society of America 86(2): 306–319.

75.

Urzúa

Dobry

Christian

(2017) Is harmonic averaging of shear wave velocity or the simplified Rayleigh method appropriate to estimate the period of a soil profile? Earthquake Spectra 33: 895–915.

76.

Vergara-Huerta

Aguirre

(2020) One-dimensional seismic velocity model of the sub-basin of Chalco, Mexico. Physics of the Earth and Planetary Interiors 300: 106426.

77.

Wang

Tsai

Stewart

Mikami

Brandenberg

(2022) Region-specific linear site amplification model for peaty organic soil sites in Hokkaido, Japan. Earthquake Spectra 38(3): 2207–2234.

78.

Wood

Deschenes

Ledezma

Meneses

Montalva

Morales-Velez

(2019) Dynamic site characterization of areas affected by the 2017 Puebla-Mexico city earthquake. Soil Dynamics and Earthquake Engineering 125: 105704.

79.

Wood

Rieth

Rahimi

(2024) Shallow three-dimensional shear wave velocity model for the Mexico City Basin. Earthquake Spectra 40(1): 314–345.

80.

Wood

Woodfield

Rahimi

Rosado-Fuentes

Sánchez-Sesma

Cruz-Jiménez

Mayoral

de la Rosa

(2023) Shear wave velocity and site period measurements for the western portion of the Mexico City Basin following the Mw7.1 2017 Puebla–Morelos, Mexico, earthquake. Earthquake Spectra 39(1): 505–527.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

0.03 MB

Valley of Mexico site response from non-reference site approach

Abstract

Keywords

Introduction

Geology and zonation

Site parameters

Data compilation

VS30 assignments

State of practice for ground motion prediction

Analysis of Valley of Mexico site response

Ground motion analysis

Site response nonlinearity

Linear site responses at recordings sites

Linear amplification model

Residuals analyses

Validation of the site response model

Summary and discussion

Supplemental Material

sj-docx-1-eqs-10.1177_87552930251316816 – Supplemental material for Valley of Mexico site response from non-reference site approach

Supplemental Material

sj-docx-2-eqs-10.1177_87552930251316816 – Supplemental material for Valley of Mexico site response from non-reference site approach

Footnotes

Acknowledgements

Declaration of conflicting interests

Funding

ORCID iDs

Data resources

Supplemental material

References

Supplementary Material

V_S30 assignments