Final Report Site Response Models
Total Page:16
File Type:pdf, Size:1020Kb
Final Report Award Number G18AC00026 Site Response Models for the Atlantic and Gulf Coastal Plain Martin C. Chapman and Zhen Guo Department of Geosciences Virginia Polytechnic Institute and State University 4044 Derring Hall Blacksburg, Virginia, 24061 email: [email protected], telephone: (540) 231-5036 September 12, 2019 Report Period April 1 2018 – March 31, 2019 This material is based upon work supported by the U.S. Geological Survey under Grant no. G18AC00026. The views and conclusions contained in this document are those of the authors and should not be interpreted as representing the opinions or policies of the U.S. Geological Survey. 1 Abstract The Atlantic and Gulf Coastal Plain in the southern and southeastern United States contains extensive Cretaceous and Cenozoic sedimentary sequences of variable thickness. We investigated the difference in response of sites in the Coastal Plain relative to sites outside that region using Fourier spectral ratios from 17 regional earthquakes occurring in 2010-2018 recorded by the EARTHSCOPE transportable array and other stations. We used mean coda and Lg spectra for sites outside the Coastal Plain as a reference. We found that Coastal Plain sites experience amplification of low-frequency ground motions and attenuation at high-frequencies relative to average site conditions outside the Coastal Plain. The spectral ratios at high frequencies gave estimates of the difference between kappa at Coastal Plain sites and the reference condition. Differential kappa values determined from the coda are correlated with the thickness of the sediment section and agree with previous estimates determined from Lg-waves. Averaged estimates of kappa reach ~ 120 ms at Gulf coast stations overlying ~12 km of sediments. Relations between Lg spectral ratio amplitudes versus sediment thickness in successive frequency bins exhibit consistent patterns, which were modeled using piecewise linear functions at frequencies ranging from 0.1 to 2.8 Hz. For sediment thickness greater than ~ 0.5 km, the spectral amplitude ratio at frequencies higher than approximately ~3 Hz is controlled by the value of kappa. The peak frequency and maximum relative amplification at frequencies less than ~1.0 Hz depend on sediment thickness. At 0.1 Hz, the mean Fourier amplitude ratio (Coastal Plain/ reference) is about 2.7 for sediment of 12 km thickness. Analysis of residuals between observed and predicted ground motions suggests that incorporating the amplification and attenuation as functions of sediment thickness may improve ground motion prediction models for the Coastal Plain region. 2 Investigations Undertaken We explored some of the differences in ground motion propagation for sites located within, versus outside, the Atlantic and Gulf Coastal Plain. Currently, ground motion prediction models are lacking for the Coastal Plain regions of the central and eastern United States. This study was a first step toward developing such a model. Figure 1 shows the study area and the seismic stations contributing data. Until recently, the relatively low levels of seismic activity and a lack of long-term operating seismic stations outside the New Madrid seismic zone have limited wave propagation studies in most parts of the Atlantic and Gulf Coastal Plain. Previous studies documented high-frequency attenuation in the Gulf Coastal Plain. Gupta et al. (1989) found lower Lg Q in the Gulf region than elsewhere in eastern North America. The operation of the Earthscope USArray Transportable Array (TA network) in the central United States during 2010-2012 when a series of moderate earthquakes occurred in Oklahoma, Arkansas and Texas resulted in an important data set. More recently, a few earthquakes have been recorded in the Atlantic Coastal Plain by the currently operational Central and Eastern United States Network (N4 network), the United States National Seismic Network (US network) and some other stations, including temporary array deployments. Pasyanos (2013), using Q tomography, found lower Q for crustal S waves in the Gulf coastal region than in regions to the north. Chapman and Conn (2016) observed geographic variation of the attenuation parameter kappa, 0, (Anderson and Hough, 1984) in the Gulf Coastal Plain, noting a clear positive correlation of 0 and the thickness of post-Jurassic sediments in the region. Incorporating a thickness-dependent Lg kappa model in stochastic ground motion simulations resulted in improved high-frequency ground motion prediction (Chapman and Conn, 2016). Figure 10 of Chapman and Conn (2016) shows that the area with largest kappa (and thickest Coastal Plain sediment) largely corresponds with the low Q Gulf Coastal Plain area resolved by Cramer (2018) using USArray (TA network) data. It also corresponds with the area that has experienced continental crustal thinning (Salvador, 1991a, Sawyer et al., 1991, Thomas, 2010). Lg propagation is known to be sensitive to changes in crustal structure (Kennett, 1986). Both Lg blockage due to crustal thinning, and absorption due to the increase in thickness of sediments may operate to increase attenuation in parts of the Gulf Coastal Plain. Chapman and Conn (2016) jointly estimated shear wave crustal Q associated with distance dependent attenuation and site terms for Lg wave Fourier amplitude spectra. They used the site 0.62 terms to estimate 0 in the Gulf Region. They found Q =365f , where f is frequency in Hz. Chapman and Conn concluded that the bulk of the attenuation in the Gulf Coastal Plain is not strongly related to crustal waveguide Q, but instead is dominated by near-receiver attenuation reflected by kappa values that are correlated with local sediment thickness. Recently, Cramer (2018) estimated Q = 259f0.72 for the Gulf coastal region. Relative to the results of Chapman and Conn (2016) and Cramer (2018), representative estimates of Q outside the Coastal Plain in eastern North America show higher values at 1 Hz by a factor of approximately 1.4 – 2.0, but significantly less frequency dependence. For example, Atkinson and Boore (2014) found Q= 525f0.45 for rock sites in eastern North America. These models predict lower Q in the Coastal Plain at frequencies less than approximately 8-14 Hz, but higher Q at higher frequencies. Purely on the basis of these crustal Q estimates, one might expect lower amplitudes in the Gulf at low frequencies, and similar or larger amplitudes at frequencies of approximately 12 Hz, relative to the average of sites outside the Coastal Plain. In this study we observed that Coastal Plain sites 3 Figure 1. Geologic map of the central and eastern United States. Locations and station codes of the Earthscope Transportable Array (TA) stations (triangles), the United States National Seismic Network (US) stations (hexagons), the Central and Eastern US Network (N4) stations (circles), the Lamont-Doherty Cooperative Seismographic Network (LD) stations (stars) and the Southeastern Suture of the Appalachian Margin Experiment (Z9) stations (diamonds) used in this study are indicated. The thick solid curve shows the boundary of the Atlantic and Gulf Coastal Plain. Adapted from Garrity and Soller (2009). exhibit smaller high frequency amplitudes and larger low-frequency amplitudes than average site conditions outside the Coastal Plain. The origin of the strong frequency dependence of the reported estimates of crustal Q(f) for the Gulf Coastal Plain may represent complex trade-offs between site terms, source terms and distance dependent attenuation parameters in the regression models used to invert for crustal Q. It is our view that ground motion prediction models for Coastal Plain sites will require information in addition to the estimated value of Q for the crustal waveguide. The higher frequency (greater than 1 or 2 Hz) attenuation as well as the amplitude and frequency range of low frequency amplification we observe in the Coastal Plain is geographically variable and is dependent on the thickness of sediments (Chapman and Conn, 2016). 4 The motivation for this study was simple. We attempted to quantify, in a straightforward way, the relative difference between site response in the Coastal Plain (Atlantic and Gulf) and the region outside the Coastal Plain in terms of the Fourier amplitude spectra. We focused on spectral ratios because we wanted to establish a basis for modifying existing or future ground motion prediction models established for rock-like conditions for application in the Coastal Plain. The existing ground motion prediction models are to a large degree founded on results derived from the stochastic method of ground motion simulation, and our approach here is amenable to the development of target spectra for stochastic simulation. We expanded the dataset used by Chapman and Conn (2016) by adding broadband stations in addition to the TA network and data from a few more recent earthquakes including some occurring in the Appalachian region. Selecting reference sites is an important step in the spectral ratio method (e.g., Borcherdt, 1970). Our study is handicapped by a lack of information on shallow geologic conditions and near-surface velocity at the great majority of recording locations. Most of the stations outside the Coastal Plain are not sited on hard rock outcrop, but instead have site conditions ranging from thin residual soil over hard crystalline rock (many sites in the Appalachian Piedmont), to sites on thick sequences of Paleozoic sedimentary rock (e.g., stations in the Appalachian Valley and Ridge, and many stations in the mid-continent area). We used mean coda and Lg spectra derived from large numbers of stations outside the Coastal Plain as the reference condition. This approach is simple, but it lacks rigor and introduces some ambiguity. Most of our data are from recent (post-2009) shocks occurring outside the Coastal Plain region. However, we find evidence that shocks occurring within the Coastal Plain produce motions outside the Coastal Plain that have reduced amplitudes at high frequency, an observation that suggests that Lg waves experience appreciable high-frequency attenuation near the source if in the Coastal Plain.