SEARCH · Geology Reports

Results for “Bulletin of Seismological Society of America”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14Linked to original sources

USGS2005

Evidence for a global seismic-moment release sequence

Temporal clustering of the larger earthquakes (foreshock-mainshock-aftershock) followed by relative quiescence (stress shadow) are characteristic of seismic cycles along plate boundaries. A global seismic-moment release history, based on a little more than 100 years of instrumental earthquake data in an extended version of the catalog of Pacheco and Sykes (1992), illustrates similar behavior for Earth as a whole. Although the largest earthquakes have occurred in the circum-Pacific region, an analysis of moment release in the hemisphere antipodal to the Pacific plate shows a very similar pattern. Monte Carlo simulations confirm that the global temporal clustering of great shallow earthquakes during 1952-1964 at M ??? 9.0 is highly significant (4% random probability) as is the clustering of the events of M ??? 8.6 (0.2% random probability) during 1950-1965. We have extended the Pacheco and Sykes (1992) catalog from 1989 through 2001 using Harvard moment centroid data. Immediately after the 1950-1965 cluster, significant quiescence at and above M 8.4 begins and continues until 2001 (0.5% random probability). In alternative catalogs derived by correcting for possible random errors in magnitude estimates in the extended Pacheco-Sykes catalog, the clustering of M ??? 9 persists at a significant level. These observations indicate that, for great earthquakes, Earth behaves as a coherent seismotectonic system. A very-large-scale mechanism for global earthquake triggering and/or stress transfer is implied. There are several candidates, but so far only viscoelastic relaxation has been modeled on a global scale.

Bulletin of the Seismological Society of America↗

USGS2005

Hypocenter locations in finite-source rupture models

We use a database of more than 80 finite-source rupture models for more than 50 earthquakes (Mw, 4.1-8.1) with different faulting styles occurring in both tectonic and subduction environments to analyze the location of the hypocenter within the fault and to consider the correlation between hypocenter location and regions of large slip. Rupture in strike-slip and crustal dip-slip earthquakes tends to nucleate in the deeper sections of the fault; subduction earthquakes do not show this tendency. Ratios of the hypocentral slip to either the average or the maximum slip show that rupture can nucleate at locations with any level of relative displacement. Rupture nucleates in regions of very large slip (D ??? 2/3 Dmax,) in only 16% of the events, in regions of large slip (1/3 Dmax < D < 2/3 Dmax,) in 35% of the events, and in regions of low slip (D ??? 1/3 Dmax) in 48% of the events. These percentages significantly exceed the percentages of fault area with very large (???7%) and large (???28%) slip. Ruptures that nucleate in regions of low slip, however, tend to nucleate close to regions of large slip and encounter a zone of very large slip within half the total rupture length. Applying several statistical tests we conclude that hypocenters are not randomly located on a fault but are located either within or close to regions of large slip.

Bulletin of the Seismological Society of America↗

USGS2005

Neotectonics of the offshore Oak Ridge fault near Ventura, southern California

The Oak Ridge fault is a large-offset, south-dipping reverse fault that forms the south boundary of the Ventura Basin in southern California. Previous research indicates that the Oak Ridge fault south of the town of Ventura has been inactive since 200-400 ka ago and that the fault tip is buried by ??? 1 km of Quaternary sediment. However, very high-resolution and medium-resolution seismic reflection data presented here show a south-dipping fault, on strike with the Oak Ridge fault, that is truncated at 80 m depth by an unconformity that is probably at the base of late Pleistocene and Holocene sediment. Furthermore, if vertically aligned features in seismic reflection data are eroded remnants of fault scarps, then a subsidiary fault within the Oak Ridge system deforms the shallowest imaged sediment layers. We propose that this subsidiary fault has mainly left-slip offset. These observations of Holocene slip on the Oak Ridge fault system suggest that revision of the earthquake hazard for the densely populated Santa Clara River valley and the Oxnard coastal plain may be needed.

Bulletin of the Seismological Society of America↗

USGS2005

Orientation of three-component geophones in the San Andreas Fault observatory at depth Pilot Hole, Parkfield, California

To identify and constrain the target zone for the planned SAFOD Main Hole through the San Andreas Fault (SAF) near Parkfield, California, a 32-level three-component (3C) geophone string was installed in the Pilot Hole (PH) to monitor and improve the locations of nearby earthquakes. The orientation of the 3C geophones is essential for this purpose, because ray directions from sources may be determined directly from the 3D particle motion for both P and S waves. Due to the complex local velocity structure, rays traced from explosions and earthquakes to the PH show strong ray bending. Observed azimuths are obtained from P-wave polarization analysis, and ray tracing provides theoretical estimates of the incoming wave field. The differences between the theoretical and the observed angles define the calibration azimuths. To investigate the process of orientation with respect to the assumed velocity model, we compare calibration azimuths derived from both a homogeneous and 3D velocity model. Uncertainties in the relative orientation between the geophone levels were also estimated for a cluster of 36 earthquakes that was not used in the orientation process. The comparison between the homogeneous and the 3D velocity model shows that there are only minor changes in these relative orientations. In contrast, the absolute orientations, with respect to global North, were significantly improved by application of the 3D model. The average data residual decreased from 13?? to 7??, supporting the importance of an accurate velocity model. We explain the remaining residuals by methodological uncertainties and noise and with errors in the velocity model.

Bulletin of the Seismological Society of America↗

USGS2005

The most recent large earthquake on the Rodgers Creek fault, San Francisco bay area

The Rodgers Creek fault (RCF) is a principal component of the San Andreas fault system north of San Francisco. No evidence appears in the historical record of a large earthquake on the RCF, implying that the most recent earthquake (MRE) occurred before 1824, when a Franciscan mission was built near the fault at Sonoma, and probably before 1776, when a mission and presidio were built in San Francisco. The first appearance of nonnative pollen in the stratigraphic record at the Triangle G Ranch study site on the south-central reach of the RCF confirms that the MRE occurred before local settlement and the beginning of livestock grazing. Chronological modeling of earthquake age using radiocarbon-dated charcoal from near the top of a faulted alluvial sequence at the site indicates that the MRE occurred no earlier than A.D. 1690 and most likely occurred after A.D. 1715. With these age constraints, we know that the elapsed time since the MRE on the RCF is more than 181 years and less than 315 years and is probably between 229 and 290 years. This elapsed time is similar to published recurrence-interval estimates of 131 to 370 years (preferred value of 230 years) and 136 to 345 years (mean of 205 years), calculated from geologic data and a regional earthquake model, respectively. Importantly, then, the elapsed time may have reached or exceeded the average recurrence time for the fault. The age of the MRE on the RCF is similar to the age of prehistoric surface rupture on the northern and southern sections of the Hayward fault to the south. This suggests possible rupture scenarios that involve simultaneous rupture of the Rodgers Creek and Hayward faults. A buried channel is offset 2.2 (+ 1.2, - 0.8) m along one side of a pressure ridge at the Triangle G Ranch site. This provides a minimum estimate of right-lateral slip during the MRE at this location. Total slip at the site may be similar to, but is probably greater than, the 2 (+ 0.3, - 0.2) m measured previously at the nearby Beebe Ranch site.

Bulletin of the Seismological Society of America↗

USGS2005

Calculation of broadband time histories of ground motion, Part II: Kinematic and dynamic modeling using theoretical Green's functions and comparison with the 1994 northridge earthquake

In the evolution of methods for calculating synthetic time histories of ground motion for postulated earthquakes, kinematic source models have dominated to date because of their ease of application. Dynamic models, however, which incorporate a physical relationship between important faulting parameters of stress drop, slip, rupture velocity, and rise time, are becoming more accessible. This article compares a class of kinematic models based on the summation of a fractal distribution of subevent sizes with a dynamic model based on the slip-weakening friction law. Kinematic modeling is done for the frequency band 0.2 to 10.0. Hz, dynamic models are calculated from 0.2 to 2.0. Hz. The strong motion data set for the 1994 Northridge earthquake is used to evaluate and compare the synthetic time histories. Source models are propagated to the far field by convolution with 1D and 3D theoretical Green’s functions. In addition, the kinematic model is used to evaluate the importance of propagation path effects: velocity structure, scattering, and nonlinearity. At present, the kinematic model gives a better broadband fit to the Northridge ground motion than the simple slip-weakening dynamic model. In general, the dynamic model overpredicts rise times and produces insufficient shorter-period energy. Within the context of the slip-weakening model, the Northridge ground motion requires a short slip-weakening distance, on the order of 0.15 m or less. A more complex dynamic model including rate weakening or one that allows shorter rise times near the hypocenter may fit the data better.

Bulletin of the Seismological Society of America↗

USGS2005

On pads and filters: Processing strong-motion data

Processing of strong-motion data in many cases can be as straightforward as filtering the acceleration time series and integrating to obtain velocity and displacement. To avoid the introduction of spurious low-frequency noise in quantities derived from the filtered accelerations, however, care must be taken to append zero pads of adequate length to the beginning and end of the segment of recorded data. These padded sections of the filtered acceleration need to be retained when deriving velocities, displacements, Fourier spectra, and response spectra. In addition, these padded and filtered sections should also be included in the time series used in the dynamic analysis of structures and soils to ensure compatibility with the filtered accelerations.

Bulletin of the Seismological Society of America↗

USGS2006

Seismic attenuation structure of the Seattle Basin, Washington State from explosive-source refraction data

We used waveform data from the 1999 SHIPS (Seismic Hazard Investigation of Puget Sound) seismic refraction experiment to constrain the attenuation structure of the Seattle basin, Washington State. We inverted the spectral amplitudes of compressional- and shear-wave arrivals for source spectra, site responses, and one- and two-dimensional Q-1 models at frequencies between 1 and 40 Hz for P waves and 1 and 10 Hz for S waves. We also obtained Q-1 models from t* values calculated from the spectral slopes of P waves between 10 and 40 Hz. One-dimensional inversions show that Qp at the surface is 22 at 1 Hz, 130 at 5 Hz, and 390 at 20 Hz. The corresponding values at 18 km depth are 100, 440, and 1900. Qs at the surface is 16 and 160 at 1 Hz and 8 Hz, respectively, increasing to 80 and 500 at 18 km depth. The t* inversion yields a Qp model that is consistent with the amplitude inversions at 20 and 30 Hz. The basin geometry is clearly resolved in the t* inversion, but the amplitude inversions only imaged the basin structure after removing anomalously high-amplitude shots near Seattle. When these shots are removed, we infer that Q-1 values may be ???30% higher in the center of the basin than the one-dimensional models predict. We infer that seismic attenuation in the Seattle basin will significantly reduce ground motions at frequencies at and above 1 Hz, partially countering amplification effects within the basin.

Bulletin of the Seismological Society of America↗

USGS2005

Stress distribution along the Fairweather-Queen Charlotte transform fault system

Tectonic loading and Coulomb stress transfer are modeled along the right-lateral Fairweather-Queen Charlotte transform fault system using a threedimensional boundary element program. The loading model includes slip below 12 km along the transform as well as motion of the Pacific plate, and it is consistent with most available Global Positioning System (GPS) displacement rate data. Coulomb stress transfer is shown to have been a weak contributing factor in the failure of the southeastern (Sitka) segment of the Fairweather fault in 1972, hastening the occurrence of the earthquake by only about 8 months. Failure of the Sitka segment was enhanced by a combination of cumulative loading from below (95%) by slip of about 5 cm/yr since 1848, by stress transfer (about 1%) from major earthquakes on straddling segments of the Queen Charlotte fault (M 8.1 in 1949) and the Fairweather fault (M 7.8 in 1958), and by viscoelastic relaxation (about 4%) following the great 1964 Alaska earthquake, modeled by Pollitz et al. (1998). Cumulative stress increases in excess of 7 MPa at a depth of 8 km are projected prior to the M 7.6 earthquake. Coulomb stress transferred by the rupture of the great M 9.2 Alaska earthquake in 1964 (Bufe, 2004a) also hastened the occurrence of the 1972 event, but only by a month or two. Continued tectonic loading over the last half century and stress transfer from the M 7.6 Sitka event has resulted in restressing of the adjacent segments by about 3 MPa at 8 km depth. The occurrence of a M 6.8 earthquake on the northwestern part of the Queen Charlotte fault on 28 June 2004, the largest since 1949, also suggests increased stress. The Cape St. James segment of the fault immediately southeast of the 1949 Queen Charlotte rupture has accumulated about 6 MPa at 8 km through loading since 1900 and stress transfer in 1949. A continued rise in earthquake hazard is indicated for the Alaska panhandle and Queen Charlotte Islands region in the decades ahead as the potential for damaging earthquakes increases.

Bulletin of the Seismological Society of America↗

USGS2005

Erratum: Empirical evidence for acceleration-dependent amplification factors

Incorrect versions of Figures 5 and 6 containing normalization errors were accidentally published by Borcherdt (2002) . They should be replaced with the figures shown here. The text and tabulated regression values published in Borcherdt (2002) are correct and refer to the figures shown here.

Bulletin of the Seismological Society of America↗

USGS2006

A stochastic estimate of ground motion at Oceano, California, for the M 6.5 22 December 2003 San Simeon earthquake, derived from aftershock recordings

The U.S. Geological Survey deployed a digital seismic station in Oceano, California, in February 2004, to investigate the cause of damage and liquefaction from the 22 December 2003 M 6.5 San Simeon earthquake. This station recorded 11 M > 2.8 aftershocks in almost 8 weeks. We analyze these recordings, together with recordings of the mainshock and the same aftershocks obtained from nearby stations in Park Hill and San Luis Obispo, to estimate the mainshock ground motion in Oceano. We estimate the Fourier amplitude spectrum using generalized spectral ratio analysis. We test a set of aftershocks as Green's functions by comparing simulated and recorded acceleration amplitude spectra for the mainshock at San Luis Obispo and Park Hill. We convolve the aftershock accelerograms with a stochastic operator to simulate the duration and phase of the mainshock accelerograms. This approximation allows us to extend the range of aftershocks that can be used as Green's functions to events nearly three magnitude units smaller than the mainshock. Our realizations for the mainshock accelerogram at Oceano yield peak ground accelerations distributed as 28% ?? 4%g. We interpret these realizations as upper bounds for the actual ground motion, because our analysis assumes a linear response, whereas the presence of liquefaction indicates that the ground behaved nonlinearly in Oceano.

Bulletin of the Seismological Society of America↗

USGS2006

Site response and attenuation in the Puget Lowland, Washington State

Simple spectral ratio (SSR) and horizontal-to-vertical (HN) site-response estimates at 47 sites in the Puget Lowland of Washington State document significant attenuation of 1.5- to 20-Hz shear waves within sedimentary basins there. Amplitudes of the horizontal components of shear-wave arrivals from three local earthquakes were used to compute SSRs with respect to the average of two bedrock sites and H/V spectral ratios with respect to the vertical component of the shear-wave arrivals at each site. SSR site-response curves at thick basin sites show peak amplifications of 2 to 6 at frequencies of 3 to 6 Hz, and decreasing spectra amplification with increasing frequency above 6 Hz. SSRs at nonbasin sites show a variety of shapes and larger resonance peaks. We attribute the spectral decay at frequencies above the amplification peak at basin sites to attenuation within the basin strata. Computing the frequency-independent, depth-dependent attenuation factor (Qs,int) from the SSR spectral decay between 2 and 20 Hz gives values of 5 to 40 for shallow sedimentary deposits and about 250 for the deepest sedimentary strata (7 km depth). H/V site responses show less spectral decay than the SSR responses but contain many of the same resonance peaks. We hypothesize that the H/V method yields a flatter response across the frequency spectrum than SSRs because the H/V reference signal (vertical component of the shear-wave arrivals) has undergone a degree of attenuation similar to the horizontal component recordings. Correcting the SSR site responses for attenuation within the basins by removing the spectral decay improves agreement between SSR and H/V estimates.

Bulletin of the Seismological Society of America↗

USGS2005

Holocene and latest Pleistocene oblique dextral faulting on the southern Inyo Mountains fault, Owens Lake basin, California

The Inyo Mountains fault (IMF) is a more or less continuous range-front fault system, with discontinuous late Quaternary activity, at the western base of the Inyo Mountains in Owens Valley, California. The southern section of the IMF trends ???N20??-40?? W for at least 12 km at the base of and within the range front near Keeler in Owens Lake basin. The southern IMF cuts across a relict early Pliocene alluvial fan complex, which has formed shutter ridges and northeast-facing scarps, and which has dextrally offset, well-developed drainages indicating long-term activity. Numerous fault scarps along the mapped trace are northeast-facing, mountain-side down, and developed in both bedrock and younger alluvium, indicating latest Quaternary activity. Latest Quaternary multiple- and single-event scarps that cut alluvium range in height from 0.5 to 3.0 m. The penultimate event on the southern IMF is bracketed between 13,310 and 10,590 cal years B.P., based on radiocarbon dates from faulted alluvium and fissure-fill stratigraphy exposed in a natural wash cut. Evidence of the most recent event is found at many sites along the mapped fault, and, in particular, is seen in an ???0.5-m northeast-facing scarp and several right-stepping en echelon ???0.5-m-deep depressions that pond fine sediment on a younger than 13,310 cal years B.P. alluvial fan. A channel that crosses transverse to this scarp is dextrally offset 2.3 ?? 0.8 m, providing a poorly constrained oblique slip rate of 0.1-0. 3 m/ k.y. The identified tectonic geomorphology and sense of displacement demonstrate that the southern IMF accommodates predominately dextral slip and should be integrated into kinematic fault models of strain distribution in Owens Valley.

Bulletin of the Seismological Society of America↗

USGS2006

Fault parameter constraints using relocated earthquakes: A validation of first-motion focal-mechanism data

We estimate the strike and dip of three California fault segments (Calaveras, Sargent, and a portion of the San Andreas near San Jaun Bautistia) based on principle component analysis of accurately located microearthquakes. We compare these fault orientations with two different first-motion focal mechanism catalogs: the Northern California Earthquake Data Center ( NCEDC ) catalog, calculated using the FPFIT algorithm (Reasenberg and Oppenheimer, 1985), and a catalog created using the HASH algorithm that tests mechanism stability relative to seismic velocity model variations and earthquake location (Hardebeck and Shearer, 2002). We assume any disagreement (misfit >30° in strike, dip, or rake) indicates inaccurate focal mechanisms in the catalogs. With this assumption, we can quantify the parameters that identify the most optimally constrained focal mechanisms. For the NCEDC/FPFIT catalogs, we find that the best quantitative discriminator of quality focal mechanisms is the station distribution ratio ( STDR ) parameter, an indicator of how the stations are distributed about the focal sphere. Requiring STDR > 0.65 increases the acceptable mechanisms from 34%–37% to 63%–68%. This suggests stations should be uniformly distributed surrounding, rather than aligning, known fault traces. For the HASH catalogs, the fault plane uncertainty ( FPU ) parameter is the best discriminator, increasing the percent of acceptable mechanisms from 63%–78% to 81%–83% when FPU ≤ 35°. The overall higher percentage of acceptable mechanisms and the usefulness of the formal uncertainty in identifying quality mechanisms validate the HASH approach of testing for mechanism stability.

Bulletin of the Seismological Society of America↗

USGS2005

Blind shear-wave velocity comparison of ReMi and MASW results with boreholes to 200 m in Santa Clara Valley: Implications for earthquake ground-motion assessment

Multichannel analysis of surface waves (MASW) and refraction microtremor (ReMi) are two of the most recently developed surface acquisition techniques for determining shallow shear-wave velocity. We conducted a blind comparison of MASW and ReMi results with four boreholes logged to at least 260 m for shear velocity in Santa Clara Valley, California, to determine how closely these surface methods match the downhole measurements. Average shear-wave velocity estimates to depths of 30, 50, and 100 m demonstrate that the surface methods as implemented in this study can generally match borehole results to within 15% to these depths. At two of the boreholes, the average to 100 m depth was within 3%. Spectral amplifications predicted from the respective borehole velocity profiles similarly compare to within 15 % or better from 1 to 10 Hz with both the MASW and ReMi surface-method velocity profiles. Overall, neither surface method was consistently better at matching the borehole velocity profiles or amplifications. Our results suggest MASW and ReMi surface acquisition methods can both be appropriate choices for estimating shearwave velocity and can be complementary to each other in urban settings for hazards assessment.

Bulletin of the Seismological Society of America↗

USGS2004

Triggered deformation and seismic activity under Mammoth Mountain in Long Valley caldera by the 3 November 2002 Mw 7.9 Denali fault earthquake

The 3 November 2002 Mw 7.9 Denali fault earthquake triggered deformational offsets and microseismicity under Mammoth Mountain (MM) on the rim of Long Valley caldera, California, some 3460 km from the earthquake. Such strain offsets and microseismicity were not recorded at other borehole strain sites along the San Andreas fault system in California. The Long Valley offsets were recorded on borehole strainmeters at three sites around the western part of the caldera that includes Mammoth Mountain - a young volcano on the southwestern rim of the caldera. The largest recorded strain offsets were -0.1 microstrain at PO on the west side of MM, 0.05 microstrain at MX to the southeast of MM, and -0.025 microstrain at BS to the northeast of MM with negative strain extensional. High sample rate strain data show initial triggering of the offsets began at 22:30 UTC during the arrival of the first Rayleigh waves from the Alaskan earthquake with peak-to-peak dynamic strain amplitudes of about 2 microstrain corresponding to a stress amplitude of about 0.06 MPa. The strain offsets grew to their final values in the next 10 min. The associated triggered seismicity occurred beneath the south flank of MM and also began at 22:30 UTC and died away over the next 15 min. This relatively weak seismicity burst included some 60 small events with magnitude all less than M = 1. While poorly constrained, these strain observations are consistent with triggered slip and intrusive opening on a north-striking normal fault centered at a depth of 8 km with a moment of l016 N m, or the equivalent of a M 4.3 earthquake. The cumulative seismic moment for the associated seismicity burst was more than three orders of magnitude smaller. These observations and this model resemble those for the triggered deformation and slip that occurred beneath the north side of MM following the 16 October 1999 M 7.1 Hector Mine, California, earthquake. However, in this case, we see little post-event slip decay reflected in the strain data after the Rayleigh-wave arrivals from the Denali fault earthquake and onset of triggered seismicity did not lag the triggered deformation by 20 min. These observations are also distinctly different from the more widespread and energetic seismicity and deformation triggered by the 1992 M 7.3 Landers earthquake in the Long Valley caldera. Thus, each of the three instances of remotely triggered unrest in Long Valley caldera recorded to date differ. In each case, however, the deformation moment inferred from the strain meter data was more than an order of magnitude larger than the cumulative moment for the associated triggered seismicity.

California↗

USGS2004

Radiated energy and the rupture process of the Denali fault earthquake sequence of 2002 from broadband teleseismic body waves

Displacement, velocity, and velocity-squared records of P and SH body waves recorded at teleseismic distances are analyzed to determine the rupture characteristics of the Denali fault, Alaska, earthquake of 3 November 2002 ( M W 7.9, M e 8.1 ). Three episodes of rupture can be identified from broadband ( ∼0.1–5.0 Hz ) waveforms. The Denali fault earthquake started as a M W 7.3 thrust event. Subsequent right-lateral strike-slip rupture events with centroid depths of 9 km occurred about 22 and 49 sec later. The teleseismic P waves are dominated by energy at intermediate frequencies (0.1-1 Hz) radiated by the thrust event, while the SH waves are dominated by energy at lower frequencies (0.05-0.2 Hz) radiated by the strike-slip events. The strike-slip events exhibit strong directivity in the teleseismic SH waves. Correcting the recorded P -wave acceleration spectra for the effect of the free surface yields an estimate of 2.8 × 10 15 N m for the energy radiated by the thrust event. Correcting the recorded SH -wave acceleration spectra similarly yields an estimate of 3.3 × 10 16 N m for the energy radiated by the two strike-slip events. The average rupture velocity for the strike-slip rupture process is 1.1 β –1.2 β . The strike-slip events were located 90 and 188 km east of the epicenter. The rupture length over which significant or resolvable energy is radiated is, thus, far shorter than the 340-km fault length over which surface displacements were observed. However, the seismic moment released by these three events, 4 × 10 20 N m, was approximately half the seismic moment determined from very low-frequency analyses of the earthquake. The difference in seismic moment can be reasonably attributed to slip on fault segments that did not radiate significant or coherent seismic energy. These results suggest that very large and great strike-slip earthquakes can generate stress pulses that rapidly produce substantial slip with negligible stress drop and little discernible radiated energy on fault segments distant from the initial point of nucleation. The existence of this energy-deficient rupture mode has important implications for the evaluation of the seismic hazard of very large strike-slip earthquakes.

Alaska↗

USGS2004

Ground motion in Anchorage, Alaska, from the 2002 Denali fault earthquake: Site response and displacement pulses

Data from the 2002 Denali fault earthquake recorded at 26 sites in and near Anchorage, Alaska, show a number of systematic features important in studies of site response and in constructing long-period spectra for use in earthquake engineering. The data demonstrate that National Earthquake Hazards Reduction Program (NEHRP) site classes are a useful way of grouping stations according to site amplification. In general, the sites underlain by lower shear-wave velocities have higher amplification. The amplification on NEHRP class D sites exceeds a factor of 2 relative to an average of motions on class C sites. The amplifications are period dependent. They are in rough agreement with those from previous studies, but the new data show that the amplifications extend to at least 10 sec, periods longer than considered in previous studies. At periods longer than about 14 sec, all sites have motion of similar amplitude, and the ground displacements are similar in shape, polarization, and amplitude for all stations. The displacement ground motion is dominated by a series of four pulses, which are associated with the three subevents identified in inversion studies (the first pulse is composed of P waves from the first subevent). Most of the high-frequency ground motion is associated with the S waves from subevent 1. The pulses from subevents 1 and 2, with moment releases corresponding to M 7.1 and 7.0, are similar to the pulse of displacement radiated by the M 7.1 Hector Mine earthquake. The signature from the largest subevent (M 7.6) is more subdued than those from the first two subevents. The two largest pulses produce response spectra with peaks at a period of about 15 sec. The spectral shape at long periods is in good agreement with the recent 2003 NEHRP code spectra but is in poor agreement with the shape obtained from Eurocode 8.

Alaska↗