Site Response Models Research Articles

Two-station Lg wave attenuation tomography in Eastern Asia

SUMMARY Utilizing over 31 000 Lg waveforms from 136 crustal earthquakes recorded at 346 regional stations, we conduct detailed tomographic mappings of the Lg Q structure across Eastern Asia in a frequency range from 0.5 to 4.0 Hz. By improving the standard two-station (TS) method, we effectively correct non-unity site response ratios using site responses estimated at individual stations. This innovative approach combines the flexible recording geometry of the TS method with the precision of reversed two-station (RTS) and reversed two-event (RTE) methods, producing a comprehensive data set devoid of source and site effects for Q tomography. To address unsolvable 3-D structural effects in the Lg spectral amplitude modelling, we justify these as modelling errors with a Gaussian distribution. This approach supports our SVD-based tomographic method, allowing for effective inversion of attenuation parameters and quantitative assessment of model resolution and errors. Our results reveal a complex relationship between Lg Q and the tectonic characteristics of Eastern Asia. In well-resolved regions, low Qo (1-Hz Q) values correspond to areas with high heat flow, partial melt, thick sediment and recent tectonic-thermal activities, in contrast to high Qo values in stable, ancient crusts lacking recent tectonic activity. Rift basins are characterized by low Lg Qo, whereas flexural basins generally have high Qo basements. We also note that post-formation factors, such as sedimentation and crustal flow intrusion, significantly impact Qo values. Furthermore, Lg Q shows a complex frequency relationship, though the power-law approximation with positive power η remains useful. The frequency dependence power η is inversely related to Qo: the regions with low Qo typically have high η and vice versa. This study provides reliable attenuation tomographic and relative site response models for Lg waves in Eastern Asia, pertinent for relative geophysical studies.

Geophysical Investigation and 3D Modeling of Bedrock Morphology in an Urban Sediment-Filled Basin: The Case of Bolzano (Northern Italy)

Bedrock mapping is essential for understanding seismic amplification, particularly in sediment-filled valleys or basins. However, this can be hard in urban environments. We conducted a geophysical investigation of the sediment-filled Bolzano basin in Northern Italy, where three valleys converge. This study uses low-impact, single-station geophysical methods suitable for urban areas to address the challenges of mapping in such environments. A dataset of 574 microtremor and gravity measurements, along with three seismic reflection lines, allows for the inference of the basin’s deep bedrock morphology, even without direct stratigraphic data. The dataset facilitates a detailed analysis of the spatial patterns of resonance frequencies and amplitudes, revealing 1D and 2D characteristics of the resonances. Notably, 2D resonances predominate along the Adige valley, i.e., the deepest part of the basin with depths up to 900 m. These 2D resonances, which cannot be interpreted through simple 1D frequency-depth relationships, are better understood by integrating gravity data to develop a depth model. The study identifies resonance frequencies ranging from 0.27 to over 3 Hz in Bolzano, affecting different building types during earthquakes. Maximum resonance amplitudes occur at lower frequencies, specifically at 2D resonance sites, therefore primarily impacting high structures. The 2D resonances are directional, with the most significant amplification occurring longitudinally along the valley axes. The resulting 3D bedrock model aids in seismic site response modeling, hydrogeological studies, and geothermal exploration and provides insights into the geological history of the basin, highlighting the role of the Adige Valley as a major drainage pathway.

Modeling of kappa factor using multivariate adaptive regression splines: application to the western Türkiye ground motion dataset

The recent seismic activity on Türkiye’s west coast, especially in the Aegean Sea region, shows that this region requires further attention. The region has significant seismic hazards because of its location in an active tectonic regime of North–South extension with multiple basin structures on soft soil deposits. Recently, despite being 70 km from the earthquake source, the Samos event (with a moment magnitude of 7.0 on October 30, 2020) caused significant localized damage and collapse in the Izmir city center due to a combination of basin effects and structural susceptibility. Despite this activity, research on site characterization and site response modeling, such as local velocity models and kappa estimates, remains sparse in this region. Kappa values display regional characteristics, necessitating the use of local kappa estimations from previous earthquake data in region–specific applications. Kappa estimates are multivariate and incorporate several characteristics such as magnitude and distance. In this study, we assess and predict the trend in mean kappa values using three–component strong–ground motion data from accelerometer sites with known VS30 values throughout western Türkiye. Multiple linear regression (MLR) and multivariate adaptive regression splines (MARS) were used to build the prediction models. The effects of epicentral distance Repi, magnitude Mw, and site class (VS30) were investigated, and the contributions of each parameter were examined using a large dataset containing recent seismic activity. The models were evaluated using well–known statistical accuracy criteria for kappa assessment. In all performance measures, the MARS model outperforms the MLR model across the selected sites.

Investigation of Site Response in Xinjiang, Northwestern China, from S–Lg Fourier Amplitude Spectra

ABSTRACT Abundant seismic waveforms have been collected in Xinjiang since 2007, following the deployment of permanent stations by the China National Seismic Network and the Xinjiang Seismic Network. We investigated regional attenuation of the S–Lg phase in the crustal waveguide and site response in Xinjiang using broadband recordings at 44 stations from 118 earthquakes with magnitudes of 4.4–7.0 occurring in Xinjiang and adjacent areas between January 2009 and February 2022. We used linear regression analysis of the Fourier acceleration amplitude spectra in 18 frequency bands from 0.1 to 24.86 Hz to estimate site response terms and distance-dependent anelastic attenuation quality factor Q in the crustal waveguide at regional distances. The regression site terms of stations located in the Tarim and Junggar Basins show significant amplification at low frequencies and strong attenuation at higher frequencies due to the presence of thick sedimentary deposits. We quantified the linear behavior of the logarithms of the high-frequency site terms versus frequency using the zero distance attenuation parameter κ0 (Anderson and Hough, 1984). A positive correlation was observed between κ0 and the thickness of the sedimentary layers in Xinjiang. We observed some consistent patterns in the behavior of the site response terms as a function of sediment thickness over the frequency range of 0.1–24.86 Hz. Linear functions of sediment thickness were used to model the site response in Xinjiang at each of the 18 frequencies. An alternative approach that adopted the κ0 model as a function of sediment thickness can be used to estimate site responses at high frequencies (≥5.66 Hz). The results of this study demonstrate the nature of site effects on earthquake ground motions caused by the thick sedimentary deposits in Xinjiang and provide a preliminary site response model as a function of sediment thickness for this region.

Preliminary observations of an ergodic site response model in California conditioned on V<sub>s30</sub> and HVSR Parameters

Preliminary observations of an ergodic site response model in California conditioned on V<sub>s30</sub> and HVSR Parameters

A Site-Response Model for the Vertical Component of Ground-Motion Prediction Equation Using a New Site-Response Parameter TVH

ABSTRACT This study presents an alternate site-response model to the existing ground-motion prediction equation (GMPE) from our previous study with some calibration to the magnitude term. We used a new site-response proxy TVH=4×max(30,HB)/VS30 that combines VS30 (the travel-time-averaged shear-wave velocity to 30 m depth) with the engineering bedrock depth HB. VS30 is available for many strong-motion recording stations. However, GMPEs using VS30 do not directly account for the response of the soil layers between 30 m and bedrock depth. Site period TS (four times the shear-wave travel time to the engineering bedrock depth) has also been used in recent GMPEs, with TS being considered a theoretically better parameter than a pseudosite period TVS30 (four times the shear-wave travel time to a depth of 30 m). Obtaining the shear-wave velocity profiles between 30 m and HB for a deep soil site can be expensive, whereas obtaining bedrock depth based on the geotechnical description of borehole data may be relatively easy. We used velocity profiles from the Kyoshin net and Kiban–Kyoshin net strong-motion networks and found that a pseudosite period TVH and TS have an excellent correlation with a small standard deviation at all spectral periods, suggesting that TVH is a suitable site-effect parameter. The poor correlations between TVS30 and TS and between TVS30 and TVH for sites with TS>0.4 s led to poor model performance at long spectral periods. We modified a GMPE by replacing TVS30 with TVH, and we made a minor modification to the moment magnitude term of the GMPE that was necessary for the new site term. The response spectra predicted by models using TVH, TS, and TVS30 at short spectral periods up to 0.6 s are generally similar. At long spectral periods, the spectra predicted by the TVS30 are much smaller than those from the other two GMPEs. For all spectral periods, TVH is an excellent substitute for TVS30.

Applications of Nonergodic Site Response Models to ShakeAlert Case Studies in the Los Angeles Area

ABSTRACT In this study, we explore whether the Parker and Baltay (2022) site response models for the Los Angeles (LA) basin region can improve ground-motion forecasts in the U.S. Geological Survey ShakeAlert earthquake early warning system (hereafter ShakeAlert). We implement the peak ground acceleration and peak ground velocity site response models of Parker and Baltay (2022) in ShakeAlert via the earthquake information to ground-motion (hereafter eqinfo2GM) module, which predicts ground motions from the estimated earthquake parameters of magnitude, rupture length, and location. The nonergodic site response models for the greater LA area were developed using ground motions from 414 M 3–7.3 earthquakes in southern California. We test nonergodic ground-motion forecasts for five earthquakes in the LA area: the 1994 M 6.7 Northridge earthquake, the 2008 M 5.4 Chino Hills earthquake, the 2019 M 7.1 Ridgecrest earthquake, the 2020 M 4.5 South El Monte earthquake, and a synthetic M 7.8 earthquake on the southern San Andreas fault from the ShakeOut scenario, which was the basis of a statewide emergency response exercise. From the test results, we find that with the nonergodic site response applied, ShakeAlert not only alerts larger areas but can also result in longer warning times in LA region. In addition, the modified Mercalli intensity (MMI) ground-motion predictions generated by the ShakeAlert eqinfo2GM module are improved in accuracy when compared with the corresponding ShakeMap ground-truth MMI when the nonergodic site response model is applied.

Identification Protocols for Horizontal-to-Vertical Spectral Ratio Peaks

ABSTRACTPeaks in horizontal-to-vertical spectral ratios (HVSRs) of Fourier amplitudes from three-component recordings of ground vibrations without undue local anthropogenic influence are used to identify site resonances, which are an important component of site response. We address two topics: (1) how should HVSR peaks be identified and (2) are there appreciable differences in HVSR derived using different instruments recording microtremors and seismic strong ground motions? We propose identifying peaks by considering peak amplitudes relative to neighboring ordinates and peak width. The procedure incorporates a regression tree algorithm that can be tuned to conform with user preferences toward relatively “conservative” or “liberal” peak identification (producing few or many sites with peaks, respectively). We then investigate the consistency of microtremor-based HVSRs (mHVSRs) derived from seismometers and accelerometers, which show a high rate of false negatives (missed peaks) from accelerometers with a full scale of ± 2g or greater. In contrast, mHVSRs derived from collocated temporary and permanent seismometers (optimized to record teleseismic signals) have about 60%–80% consistency (with no apparent bias in peak attributes). This indicates that mHVSRs from accelerometers having a broad full scale are unreliable but that mHVSRs can be reliably obtained from temporary or permanent seismometers. Finally, we compare seismometer-based HVSR from microtremor and earthquake (eHVSRs) sources. Results are consistent for 60%–70% of sites (i.e., both either do or do not have significant peaks, and when peaks are present, they occur at similar frequencies, <20% change). For sites with an mHVSR peak, the rate of corresponding eHVSR peaks is nearly 50%, whereas for sites without an mHVSR peak the eHVSR peak rate is low (about 20%). The mismatch rate for mHVSR peak sites is sufficiently high that the use of eHVSR to derive site response models is likely too optimistic (overestimates model effectiveness); mHVSR is preferred for consistency with information available in forward applications.

A study of site response in the Longmen Shan and adjacent regions and site response models for the Sichuan Basin

We investigated the regional attenuation and site responses in the Sichuan Basin and adjacent Songpan-Ganze terrane of the Tibetan Plateau using seismic data recorded at 41 stations from regional earthquakes occurring between January 2009 and October 2020. Fourier amplitude spectra of Lg waves were computed and binned into 18 frequency bins with center frequencies ranging from 0.1 Hz to 20.4 Hz. The quality factor is estimated as Q(f)=313f0.74 for the Sichuan Basin and Q(f)=568f0.34 for the Songpan-Ganze terrane, reflecting significant differences in the crustal structure beneath these two regions. Relative to the Songpan-Ganze terrane, site responses in the Sichuan Basin are characterized by strong amplification effects at frequencies lower than 6 Hz and obvious attenuation at higher frequencies (>10 Hz). κ0 of stations in the Sichuan Basin show clearly geographical dependence with an average value of 0.045 s, whereas stations in the Songpan-Ganze terrane generally have smaller κ0 values with an average value of 0.028 s. In particular, site response and κ0 of stations in the Sichuan Basin are found to be dependent on the geographically variable thickness of the sedimentary deposits (sediment thickness). These units are comprised of sedimentary rock and semi-consolidated sediments, with a maximum thickness reaching approximately 10 km. Site response terms in the Sichuan Basin derived from the Lg Fourier spectra exhibit consistent patterns versus sediment thickness as frequency increases. We developed site response models as functions of sediment thickness for stations in the Sichuan Basin. The site response model derived from Lg site terms is consistent with that based on site response terms from coda amplitude spectra and horizontal to vertical (H/V) spectral ratios. The models were then incorporated in the stochastic method of ground motion predictions in the Sichuan Basin for six earthquakes occurring between October 2020 and June 2022. Residual analysis suggests that incorporating the site response models as functions of sediment thickness can improve the ground motion prediction model for the Sichuan Basin from moderate earthquakes.

Impact of local site effects on seismic risk assessment of reinforced concrete bridges

In most practical applications of performance-based earthquake engineering, the local site characteristics are considered simplistically, and soil-structure interaction effects are ignored. Within such a conventional approach, the ground motion record selection uses the time-averaged shear-wave velocity in the top 30 m of the soil profile (VS30) as the main proxy to represent the site where the structure is located. This study aims to assess the impact of local site effects on the risk assessment of reinforced concrete (RC) bridges by incorporating the local-site response and soil-structure interaction effects more realistically within a more detailed approach. For this purpose, seismic risk assessment of RC bridge structures with various configurations located on varying site characteristics was carried out using different levels of refinement in the numerical modelling. Local site amplification was incorporated using site response analysis, and p-y models were used to consider foundation flexibility. Cloud analyses were then carried out on both modelling approaches to obtain and compare the corresponding scenario-based seismic loss estimates for the structures. Moreover, the results of numerical analyses were scrutinized in terms of peak demand to capacity ratio values, and statistical hypothesis testing was used to quantify the similarity between the two approaches further. In some cases of site response modeling, especially when the fundamental periods of the bridge structure and site are sufficiently close, higher displacement demands on bridge piers were obtained, likely leading to higher direct seismic losses.

Empirical Map-Based Nonergodic Models of Site Response in the Greater Los Angeles Area

ABSTRACTWe develop empirical estimates of site response at seismic stations in the Los Angeles area using recorded ground motions from 414 M 3–7.3 earthquakes in southern California. The data are from a combination of the Next Generation Attenuation-West2 project, the 2019 Ridgecrest earthquakes, and about 10,000 newly processed records. We estimate site response using an iterative mixed-effects residuals partitioning approach, accounting for azimuthal variations in anelastic attenuation and potential bias due to spatial clusters of colocated earthquakes. This process yields site response for peak ground acceleration, peak ground velocity, and pseudospectral acceleration relative to a 760 m/s shear-wave velocity (VS) reference condition. We employ regression kriging to generate a spatially continuous site response model, using the linear site and basin terms from Boore et al. (2014) as the background model, which depend on VS30 and depth to the 1 km/s VS isosurface. This is different from past approaches to nonergodic models, in which spatially varying coefficients are regressed. We validate the model using stations in the Community Seismic Network (CSN) that are in the middle of our model spatial domain but were not considered in model development, finding strong agreement between the interpolated model and CSN data for long periods. Our model could be implemented in regional seismic hazard analyses, which would lead to improvements especially at long return periods. Our site response model also has potential to improve both ground-motion accuracy and warning times for the U.S. Geological Survey ShakeAlert earthquake early warning (EEW) system. For a point-source EEW simulation of the 1994 M 6.7 Northridge earthquake, our model produces ground motions more consistent with the ground-truth ShakeMap and would alert areas with high population density such as downtown Los Angeles at lower estimated magnitudes (i.e., sooner) than an ergodic model for a modified Mercalli intensity 4.5 alerting threshold.

Convergence Study of Variational Space‐Time Coupled Least‐Squares Frameworks in Simulation of Wave Propagation in Viscoelastic Medium

There are many novel applications of space‐time decoupled least squares and Galerkin formulations that simulate wave propagation through different types of media. Numerical simulation of stress wave propagation through viscoelastic medium is commonly carried out using the space‐time decoupled Galerkin weak form in site response problem, etc. In a recent investigation into accuracy of this formulation in simulating elastic wave propagation, it was shown that the diffusive and dispersive errors are greatly reduced when space‐time coupled least squares formulation is used instead in variational form. This paper investigates convergence characteristics of both formulations. To this end, two test cases, which are site response and impact models for viscoelastic medium, are employed, one dominated by dispersive and the other by diffusive numerical error. Convergence studies reveal that, compared to the commonly used space‐time decoupled Galerkin and the coupled least squares formulation has much lower numerical errors, higher rates of convergence, and ability to take far larger time increments in the evolution process. In solving such models, coefficient matrices would require updating after each time step, a process that can be very costly. However large time steps allowed by cLs are expected to be a significant feature in reducing the overall computational cost.

Relational Database for Horizontal-to-Vertical Spectral Ratios

Abstract Frequency-dependent horizontal-to-vertical spectral ratios (HVSRs) of Fourier amplitudes from three-component recordings can provide useful information for site response modeling. However, such information is not incorporated into most ground-motion models, including those from Next-Generation Attenuation projects, which instead use the time-averaged shear-wave velocity (VS) in the upper 30 m of the site and sediment depth terms. To facilitate utilization of HVSR, we developed a publicly accessible relational database. This database is adapted from a similar repository for VS data and provides microtremor-based HVSR data (mHVSR) and supporting metadata, but not parameters derived from the data. Users can interact with the data directly within a web portal that contains a graphical user interface (GUI) or through external tools that perform cloud-based computations. Within the database GUI, the median horizontal-component mHVSR can be plotted against frequency, with the mean and mean ± one standard deviation (representing variability across time windows) provided. Using external interactive tools (provided as a Jupyter Notebook and an R script), users can replot mHVSR (as in the database) or create polar plots. These tools can also derive parameters of potential interest for modeling purposes, including a binary variable indicating whether an mHVSR plot contains peaks, as well as the fitted properties of those peaks (frequencies, amplitudes, and widths). Metadata are also accessible, which includes site location, details about the instruments used to make the measurements, and data processing information related to windowing, antitrigger routines, and filtering.

Ergodic site response model for subduction zone regions

We present an ergodic site response model with regional adjustments for use with subduction zone ground-motion models. The model predicts site amplification of peak ground acceleration, peak ground velocity, and 5% damped pseudo-spectral accelerations of the orientation-independent horizonal component for oscillator periods from 0.01 to 10 s. The model depends on the time-averaged shear-wave velocity in the upper 30 m ( VS30), basin depth, and region and is independent of subduction earthquake type. It has three components: a linear site-amplification term in the form of VS30-scaling, a nonlinear term that depends on VS30 and shaking intensity parameterized by peak ground acceleration at the reference-rock velocity condition of 760 m/s, and a basin sediment-depth term for Japan and Cascadia conditioned on the depth to the 2.5 km/s shear-wave velocity isosurface ( Z2.5). A global VS30-scaling model is provided along with regional adjustments for Japan, Taiwan, South America, Alaska, and Cascadia. The nonlinear model is global, with a functional form that has often been used to fit nonlinear responses inferred from simulations, but here we calibrate it empirically. Relative to a prior model for shallow earthquakes in active tectonic regions, our subduction zone global VS30-scaling is comparable at short periods (<1.0 s) but weaker at long periods, while the nonlinear site response is generally less pronounced but extends to lower levels of shaking. Basin depth models are conditioned on the difference of the actual Z2.5 and a VS30-conditioned mean Z2.5. Sites with positive differential depths have increased long-period site responses and decreased short-period responses, with the opposite occurring for negative differential depths.

Reduction of Bias and Uncertainty in Regional Seismic Site Amplification Factors for Seismic Hazard and Risk Analysis

Site amplification factors in National Building Codes are typically specified as a function of the average shear wave velocity over the first 30 m (Vs30) or site class (A, B, C, D and E) for defined ranges of Vs30 and/or ranges of depth to bedrock. However, a single set of amplification factors may not be representative of site conditions across the country, introducing a bias in seismic hazard and seismic risk analyses. This is exemplified by significant differences in geological settings between East and West coast locations in North America. Western sites are typically characterized by lower impedance contrasts between recent surface deposits and bedrock in comparison to Eastern sites. In North America, site amplification factors have been derived from a combination of field data on ground motions recorded during West Coast earthquakes and numerical models of site responses that are meant to be representative of a wide variety of soil profiles and ground motions. The bias on amplifications and their impact on seismic hazards is investigated for the Montreal area, which ranks second for seismic risks in Canada in terms of population and hazard (PGA of 0.25 g for a 2475 years return period). Representative soil profiles at several locations in Montreal are analyzed with 1-D site response models for natural and synthetic ground motions scaled between 0.1 to 0.5 g. Since bedrock depths are typically shallow (<30 m) across the island, bedrock shear wave velocities have a significant influence on the impedance contrast and amplifications. Bedrock shear wave velocity is usually very variable due to the differences in rock formations, level of weathering and fracturing. The level of this uncertainty is shown to be greatly decreased when rock quality designation (RQD) data, common information when bore hole data are logged, is available since it is highly correlated with both shear and compression wave velocities. The results are used to derive region-specific site amplification factors as a function of both Vs30 and site fundamental frequency and compared to those of the National Building Code of Canada (2015). The results of the study indicate that there are large uncertainties associated with these parameters due to variability in soil profiles, soil properties and input seismic ground motions. Average and confidence intervals for the mean and for predictions of amplification factors are calculated for each site class to quantify this uncertainty. Amplifications normalized relative to class C are obtained by accounting for the correlation between site class amplifications for given ground motions. Non-linearity in the analysis of equivalent linear 1-D site response is taken into account by introducing the non-linear G/Gmax and damping ratios curves. In this method, it is assumed that the shear strain compatible shear modulus and damping ratio values remains constant throughout the duration of the seismic excitation. This assumption is not fully applicable to a case when loose saturated soil profile undergo heavy shaking (PGA > 0.3 g). In this study, all simulations with input motion PGA >0.3 g have been performed by using the EL method instead of the NL method considering that cohesive soils (clay and silt) at Montreal sites are stiff and cohesionless soils (sand and gravel) are considerably dense. In addition, the field and laboratory data required to perform NL analyses are not currently available and may be investigated in future works.

NGA-subduction global ground motion models with regional adjustment factors

We develop semi-empirical ground motion models (GMMs) for peak ground acceleration, peak ground velocity, and 5%-damped pseudo-spectral accelerations for periods from 0.01 to 10 s, for the median orientation-independent horizontal component of subduction earthquake ground motion. The GMMs are applicable to interface and intraslab subduction earthquakes in Japan, Taiwan, Mexico, Central America, South America, Alaska, the Aleutian Islands, and Cascadia. The GMMs are developed using a combination of data inspection, data regression with respect to physics-informed functions, ground-motion simulations, and geometrical constraints for certain model components. The GMMs capture observed differences in source and path effects for interface and intraslab events, conditioned on moment magnitude, rupture distance, and hypocentral depth. Site effect and aleatory variability models are shared between event types. Regionalized GMM components include the model constant (that controls ground motion amplitude), anelastic attenuation, magnitude-scaling break point, linear site response, and sediment depth terms. We develop models for the aleatory between-event variability [Formula: see text], within-event variability [Formula: see text], single-station within-event variability [Formula: see text], and site-to-site variability [Formula: see text]. Ergodic analyses should use the median GMM and aleatory variability computed using the between-event and within-event variability models. An analysis incorporating non-ergodic site response should use the median GMM at the reference shear-wave velocity condition, a site-specific site response model, and aleatory variability computed using the between-event and single-station within-event variability models. Epistemic uncertainty in the median model is represented by standard deviations on the regional model constants, which facilitates scaled-backbone representations of model uncertainty in hazard analyses.

Liquefaction analysis of sandy soil during strong earthquake in Northern Thailand

Northern Thailand has experienced several earthquakes which led to soil liquefaction in the past few decades. Traditional methods of evaluating liquefaction potential involve standard penetration test (SPT) or cone penetration tests. This research augmented experimental results with numerical methods to evaluate the liquefaction potential of Mae Lao Sand in Chiang Rai province of northern Thailand. SPT and downhole seismic test data collected during a field investigation at the Mae Lao site were compared to a 1D site response model analysis of the site. A series of undrained monotonic and cyclic triaxial tests was conducted on Mae Lao Sand specimens with different initial void ratios and confining pressures. Cyclic triaxial test results with varying deviator stress amplitudes were used to draw liquefaction resistance curves. Results from numerical simulation of sand liquefaction were used to characterise the stress–strain-pore water pressure response of Mae Lao Sand. 1D site response analysis determined seismic responses with different geological and groundwater conditions. All put together, the results showed that pore water pressure ratio decreases with increasing sand stiffness, the thickness of a soil layer significantly increases its liquefaction potential, and in-situ conditions and groundwater depths have major influences on the liquefaction potential of sand layers.

Simplified methods to estimate mean hazard due to updated ground motion model: Application to nuclear power plants in CEUS

Two methods were developed to estimate updated mean seismic hazard for existing probabilistic seismic hazard analyses (PSHAs) due to a change in the ground motion model (GMM). Both methods were used to estimate updated hazard at nuclear power plant (NPP) sites in the Central and Eastern United States (CEUS) for a change from the Electric Power Research Institute (EPRI) 2013 GMM to the Next Generation Attenuation (NGA)-East GMM. These methods present efficient tools to inform decisions on whether to perform a full PSHA revision or other detailed evaluations, especially when a large number of sites must be analyzed. A Simplified Hazard (SiHaz) method was developed to estimate mean hazard explicitly using a reduced PSHA logic tree that incorporates the updated GMM and potential changes in the site response model. An alternative scaling method was independently developed to be applied directly to current CEUS NPP hazard. Both methods were validated using updated PSHA results at several sites. Estimates at 46 NPP sites using both methods showed good agreement for mean annual frequencies of exceedance between 1E-4 and 1E-5/yr.

Evaluation of the Influence of Shear Strength Correction through a Comparative Study of Nonlinear Site Response Models

Evaluation of the Influence of Shear Strength Correction through a Comparative Study of Nonlinear Site Response Models

Basin effects and limitations of 1D site response analysis from 2D numerical models of the Thorndon basin

Plane strain (2D) finite element models are used to examine factors contributing to basin effects observed for multiple seismic events at sites in the Thorndon basin of Wellington, New Zealand. The models consider linear elastic soil and rock response when subjected to vertically-propagating shear waves. Depth-dependent shear wave velocities are considered in the soil layers, and the effects of random variations of soil velocity within layers are modelled. Various rock shear wave velocity configurations are considered to evaluate their effect on the modelled surficial response. It is shown that these simple 2D models are able to capture basin reverberations and compare more favourably to observations from strong motion recordings than conventional 1D site response models. It is also shown that consideration of a horizontal impedance contrast across the Wellington Fault affects spectral response and amplification at longer periods, suggesting the importance of this feature in future ground motion modelling studies in the Wellington region.