One of the current challenges in physics-based ground-motion simulations is to refine the modeling of local site effects. These effects require a finer spatial resolution in the material modeling than that generally considered in regional-scale simulations. Because of this, empirical amplification factors are typically applied to capture these unmodeled phenomena. The ergodic nature of this approach suggests that there is room for improvement. In this study, the predictive capability of simulations is evaluated using alternative methods for capturing local site effects. In addition to the conventional empirical approach, two methods are examined that allow for more site-specific information to be incorporated: the square-root impedance method and the 1D time domain site-response analysis. The three approaches are tested using 1000+ observed ground motions from 150+ small-magnitude events (3.5 ≤ Mw ≤ 5.0), recorded at 20 strong-motion stationsin the Canterbury, New Zealand, region. These 20 well-characterized sites represent a wide range of soil conditions, including stiff gravels with Vs30 values greater than 500 m/s, and sand and silt deposits with Vs30 valuesless than 200 m/s. Multiple intensity measures are computed and prediction residuals are partitioned using mixed-effects regression to rigorously assess the relative performance of the different approaches considered. The results indicate that the benefit of using more sophisticated methods is highly dependent on the characteristics of the site. Key site parameters and trends are identified and discussed in light of the assumptions and limitations of each approach.
This paper provides a summary of initial research results investigating systematic site effects from the prediction residuals of empirical- and physics-based ground-motion models (GMMs) for small magnitude (i.e., 3.5 ≤ MW ≤ 5) active shallow crustal earthquakes in New Zealand (NZ). Advancing ground-motion predictability through physics-based GMMs is an iterative process and requires addressing fundamental questions like: Is there salient physics which has been overlooked? Which geographic regions have predictions that significantly deviate from observations and why? Which sites exhibit systematic prediction residuals and how can the attributes influencing them be identified? This preliminary study examines these questions by classifying 171 sites from the Canterbury and Wellington regions into four geomorphic categories: basin, basin-edge, hill, and valley, following the categorisation by Nweke et al. (2022). Trends in the site-to-site residuals for each geomorphic category indicate apparent differences between the four categories, with residuals for valley sites illustrating a clear dependence with the inferred fundamental site period. Computed residuals from both empirical- and physics-based GMMs also provided insight into the role of site-specific attributes vs. the different prediction methods, assisting to understand the salient causes of these residuals.
The Mѡ=7.1 Darfield (Canterbury) earthquake struck on 4 September 2010, approximately 45 km west of Christchurch, New Zealand. It revealed a previously unknown fault (the Greendale fault) and caused billions of dollars of damage due to high peak ground velocities and extensive liquefaction. It also triggered the Mw=6.3 Christchurch earthquake on 22 February 2011, which caused further damage and the loss of 185 lives. The objective of this research was to determine the relationship between stress and seismic properties in a seismically active region using manually-picked P and S wave arrival times from the aftershock sequence between 8 September 2010-13 January 2011 to estimate shear-wave splitting (SWS) parameters, VP =VS-ratios, anisotropy (delay-time tomography), focal mechanisms, and tectonic stress on the Canterbury plains. The maximum horizontal stress direction was highly consistent in the plains, with an average value of SHmax=116 18 . However, the estimates showed variation in SHmax near the fault, with one estimate rotating by as much as 30° counter-clockwise. This suggests heterogeneity of stress at the fault, though the cause remains unclear. Orientations of the principal stresses predominantly indicate a strike-slip regime, but there are possible thrust regimes to the west and north/east of the fault. The SWS fast directions (ø) on the plains show alignment with SHmax at the majority of stations, indicating stress controlled anisotropy. However, structural effects appear more dominant in the neighbouring regions of the Southern Alps and Banks Peninsula.
