Geology ReportsSearch

SEARCH · Geology Reports

Results for “Bulletin Seismological Society America”

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

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

At least 73 records · Page 4Linked to original sources

A satellite-based digital data system for low-frequency geophysical data

A reliable method for collection, display, and analysis of low-frequency geophysical data from isolated sites, which can be throughout North and South America and the Pacific Rim, has been developed for use with the Geostationary Operational Environmental Satellite (GOES) system. Geophysical data primarily intended for earthquake hazard and crustal deformation monitoring are digitized with either 12-bit or 16-bit resolution and transmitted every 10 min through a satellite link to a bank of UNIX-based computers in Menlo Park, California. There the data are available for analysis and display within a few seconds of their transmit time. This system provides real-time monitoring of crustal deformation parameters such as tilt, strain, fault displacement, local magnetic field, crustal geochemistry, and water levels, as well as meteorological and other parameters, along faults in California and Alaska, and in volcanic regions in the western United States, Rabaul, and other locations in the New Britain region of the South Pacific. Various mathematical, statistical, and graphical algorithms process the incoming data to detect changes in crustal deformation and fault slip that may indicate the first stages of catastrophic fault failure. Alert trigger levels based on physical models, signal resolution, and previous history have been defined for particular instrument types. Computer-driven remote paging and mail systems are used to notify appropriate personnel when alarm status is reached. The system supports continuous historical records of low-frequency geophysical data, software for extensive analysis of these data, and programs for modeling fault rupture with and without seismic radiation, as well as providing an environment for real-time attempts at earthquake prediction.

Bulletin of the Seismological Society of America

Mid- to late Quaternary geomorphic and paleoseismic event history, Cheraw fault, southeastern Colorado

Despite its subdued expression and isolated location within the Great Plains of southeastern Colorado, the 80 km long Cheraw fault may be one of the most active faults in North America east of the Rocky Mountains. We present geomorphic analyses, geochronology, and paleoseismic trenching data to (1) document the rupture history of the ∼45 km long southwestern section of the Cheraw fault over the past ∼19 ka, and (2) evaluate slip rate changes for the entire fault over the past ∼200 ka. Results from new trenches excavated at the Old Ranch site show evidence of four surface rupture events since ∼19 ka, each with an average vertical displacement of 0.75 ± 0.15 m. An additional event is likely only slightly older than ∼19 ka. Evidence for relatively small displacements at and near the Old Ranch site indicates that most of these earthquakes were magnitude ( M ) 7 or less and likely did not rupture the full length of the Cheraw fault. Since ∼19 ka, the average slip rate is ∼0.16 ± 0.3 mm/yr near the Old Ranch site with an average interevent time of 3–5 ka. New geochronologic data for mid‐ to late‐Quaternary geomorphic surfaces cut by the Cheraw fault imply rapid incision by local Arkansas River tributaries from ∼145 to ∼100 ka. Maximum vertical offsets of 7–9 m for these surfaces indicate that from ∼19 to >200 ka the average slip rate was no greater than ∼0.03 mm/yr. The accelerated slip rate since ∼19 ka indicates a possible response to rapid erosional unloading and/or a limited late Cenozoic, <40 ka, paleoseismic history for the Cheraw fault.

Colorado

Determination of stress parameters for eight well-recorded earthquakes in eastern North America

We determined the stress parameter, Δσ , for the eight earthquakes studied by Atkinson and Boore (2006), using an updated dataset and a revised point-source stochastic model that captures the effect of a finite fault. We consider four geometrical-spreading functions, ranging from 1/ R at all distances to two- or three-part functions. The Δσ values are sensitive to the rate of geometrical spreading at close distances, with 1/ R 1.3 spreading implying much higher Δσ than models with 1/ R spreading. The important difference in ground motions of most engineering concern, however, arises not from whether the geometrical spreading is 1/ R 1.3 or 1/ R at close distances, but from whether a region of flat or increasing geometrical spreading at intermediate distances is present, as long as Δσ is constrained by data that are largely at distances of 100 km–800 km. The simple 1/ R model fits the sparse data for the eight events as well as do more complex models determined from larger datasets (where the larger datasets were used in our previous ground-motion prediction equations); this suggests that uncertainty in attenuation rates is an important component of epistemic uncertainty in ground-motion modeling. For the attenuation model used by Atkinson and Boore (2006), the average value of Δσ from the point-source model ranges from 180 bars to 250 bars, depending on whether or not the stress parameter from the 1988 Saguenay earthquake is included in the average. We also find that Δσ for a given earthquake is sensitive to its moment magnitude M , with a change of 0.1 magnitude units producing a factor of 1.3 change in the derived Δσ .

Bulletin of the Seismological Society of America

Earthquake ground-motion prediction equations for eastern North America

New earthquake ground-motion relations for hard-rock and soil sites in eastern North America (ENA), including estimates of their aleatory uncertainty (variability) have been developed based on a stochastic finite-fault model. The model incorporates new information obtained from ENA seismographic data gathered over the past 10 years, including three-component broadband data that provide new information on ENA source and path effects. Our new prediction equations are similar to the previous ground-motion prediction equations of Atkinson and Boore (1995), which were based on a stochastic point-source model. The main difference is that high-frequency amplitudes (f ??? 5 Hz) are less than previously predicted (by about a factor of 1.6 within 100 km), because of a slightly lower average stress parameter (140 bars versus 180 bars) and a steeper near-source attenuation. At frequencies less than 5 Hz, the predicted ground motions from the new equations are generally within 25% of those predicted by Atkinson and Boore (1995). The prediction equations agree well with available ENA ground-motion data as evidenced by near-zero average residuals (within a factor of 1.2) for all frequencies, and the lack of any significant residual trends with distance. However, there is a tendency to positive residuals for moderate events at high frequencies in the distance range from 30 to 100 km (by as much as a factor of 2). This indicates epistemic uncertainty in the prediction model. The positive residuals for moderate events at < 100 km could be eliminated by an increased stress parameter, at the cost of producing negative residuals in other magnitude-distance ranges; adjustment factors to the equations are provided that may be used to model this effect.

Bulletin of the Seismological Society of America

Influence of lithostatic stress on earthquake stress drops in North America

We estimate stress drops for earthquakes in and near the continental United States using the method of spectral ratios. The ratio of acceleration spectra between collocated earthquakes recorded at a given station removes the effects of path and recording site and yields source parameters including corner frequency for, and the ratio of seismic moment between, the two earthquakes. We determine stress drop from these parameters for 1121 earthquakes greater than M ∼3 in 60 earthquake clusters. We find that the average Brune stress drop for the few eastern United States (EUS) tectonic mainshocks studied (2.6–36 MPa) is about three times greater than that of tectonic mainshocks in the western United States (WUS, 1.0–7.9 MPa) and five times greater than mainshocks potentially induced by wastewater injection in the central United States (CUS, 0.6–5.6 MPa). EUS events tend to be deeper thrusting events, whereas WUS events tend to be shallower but have a wide range of focal mechanisms. CUS events tend to be shallow with strike‐slip to normal‐faulting mechanisms. With the possible exception of CUS aftershocks, we find that differences in stress drop among all events can be taken into account, within one standard deviation of significance, by differences in the shear failure stress as outlined by Mohr–Coulomb theory. The shear failure stress is a function of vertical stress (or depth), the fault style (normal, strike slip, or reverse), and coefficient of friction (estimated here to be, on average, 0.64). After accounting for faulting style and depth dependence, we find that the average Brune stress drop is about 3% of the failure stress. These results suggest that high‐frequency shaking hazard (>∼1 Hz) from shallow induced events and aftershocks is reduced to some extent by lower stress drop. However, the shallow hypocenters will increase hazard within several kilometers of the source.

Bulletin of the Seismological Society of America

Neotectonics and paleoseismology of the Limón and pedro miguel faults in Panamá: earthquake hazard to the Panamá canal

We present new geologic, tectonic geomorphic, and geochronologic data on the slip rate, timing, and size of past surface ruptures for the right-lateral Lim??n and Pedro Miguel faults in central Panam??. These faults are part of a system of conjugate faults that accommodate the internal deformation of Panam?? resulting from the ongoing collision of Central and South America. There have been at least three surface ruptures on the Lim??n fault in the past 950-1400 years, with the most recent during the past 365 years. Displacement in this young event is at least 1.2 m (based on trenching) and may be 1.6-2 m (based on small channel offsets). Awell-preserved 4.2 m offset suggests that the penultimate event also sustained significant displacement. The Holocene slip rate has averaged about 6 mm=yr, based on a 30-m offset terrace riser incised into a 5-ka abandoned channel. The Pedro Miguel fault has sustained three surface ruptures in the past 1600 years, the most recent being the 2 May 1621 earthquake that partially destroyed Panam?? Viejo. At least 2.1 m of slip occurred in this event near the Canal, with geomorphic offsets suggesting 2.5-3 m. The historic Camino de Cruces is offset 2.8 m, indicating multimeter displacement over at least 20 km of fault length. Channel offsets of 100-400 m, together with a climate-induced incision model, suggest a Late Quaternary slip rate of about 5 mm=yr, which is consistent with the paleoseismic results. Comparison of the timing of surface ruptures between the Lim??n and Pedro Miguel faults suggests that large earthquakes may rupture both faults with 2-3 m of displacement for over 40 km, such as is likely in earthquakes in the M 7 range. Altogether, our observations indicate that the Lim??n and Pedro Miguel faults represent a significant seismic hazard to central Panam?? and, specifically, to the Canal and Panam?? City.

Bulletin of the Seismological Society of America

Neotectonic and paleoseismic analysis of the northwest extent of Holocene surface deformation along the Meers Fault, Oklahoma

TheMeers fault (Oklahoma) is one of fewseismogenic structures with evidence for Holocene surface rupture in the stable continental region of North America. The 37-kilometer-long southeast section of the full 54-kilometer-long Meers fault is interpreted to be Holocene active. The 17-kilometer-long northwest section is considered Quaternary active, but not Holocene active.We reevaluate surface expression and earthquake timing of the northwest Meers fault to improve seismic source characterization.We use airborne light detection and ranging and historical stereopaired aerial photos to evaluate the fault scarp and local faultzone geomorphology. In the northwest, complex surface deformation includes fault splays, subtle monoclinal warping, and a minor change in fault strike. We interpret that the alongstrike transition from surface faulting on the southeastMeers fault to surface folding on the northwest Meers fault occurs at the lithologic contact between Permian Post Oak conglomerate and Hennessey shale. We excavated a paleoseismic trench to evaluate the timing of surface-deforming earthquakes on the northwest section of the fault. The excavation revealed weathered Permian Hennessey shale and an ∼1–2-meter-thick veneer of Holocene alluvial deposits that were progressively deformed during two surface-folding earthquakes likely related to blind fault rupture beneath the site. Repeated onlapping to overlapping stratigraphic sequences and associated unconformities are intimately related to folding events along the monocline. OxCal paleoearthquake age modeling indicates that earthquakes occurred 4704–3109 yr B.P. and 5955–4744 yr B.P., and that part of the northwest section of the Meers fault is Holocene active. We find the Holocene-active section of the Meers fault should be lengthened 6.1 km to the northwest, to a total Holocene-active fault length of 43 km. Empirical scaling relationships between surface rupture length and magnitude reveal that the fault could generate an Mw 7.0 earthquake.

Oklahoma

Regional propagation characteristics and source parameters of earthquakes in northeastern North America

The vertical components of the S wave trains recorded on the Eastern Canadian Telemetered Network ( ECTN ) from 1980 through 1990 have been spectrally analyzed for source, site, and propagation characteristics. The data set comprises some 1033 recordings of 97 earthquakes whose magnitudes range from M ≈ 3 to 6. The epicentral distances range from 15 to 1000 km, with most of the data set recorded at distances from 200 to 800 km. The recorded S wave trains contain the phases S , S m S , S n , and L g and are sampled using windows that increase with distance; the acceleration spectra were analyzed from 1.0 to 10 Hz. To separate the source, site, and propagation characteristics, an inversion for the earthquake corner frequencies, low-frequency levels, and average attenuation parameters is alternated with a regression of residuals onto the set of stations and a grid of 14 distances ranging from 25 to 1000 km. The iteration between these two parts of the inversion converges in about 60 steps. The average attenuation parameters obtained from the inversion were Q = 1997 ± 10 and γ = 0.998 ± 0.003. The most pronounced variation from this average attenuation is a marked deamplification of more than a factor of 2 at 63 km and 2 Hz, which shallows with increasing frequency and increasing distance out to 200 km. The site-response spectra obtained for the ECTN stations are generally flat. The source spectral shape assumed in this inversion provides an adequate spectral model for the smaller events ( M o < 3 × 10 21 dyne-cm) in the data set, whose Brune stress drops range from 5 to 150 bars. For the five events in the data set with M o ≧ 10 23 dyne-cm, however, the source spectra obtained by regressing the residuals suggest that an ω 2 spectrum is an inadequate model for the spectral shape. In particular, the corner frequencies for most of these large events appear to be split, so that the spectra exhibit an intermediate behavior (where | ü ( ω )| is roughly proportional to ω ).

Bulletin of the Seismological Society of America

A time-dependent probabilistic seismic-hazard model for California

For the purpose of sensitivity testing and illuminating nonconsensus components of time-dependent models, the California Department of Conservation, Division of Mines and Geology (CDMG) has assembled a time-dependent version of its statewide probabilistic seismic hazard (PSH) model for California. The model incorporates available consensus information from within the earth-science community, except for a few faults or fault segments where consensus information is not available. For these latter faults, published information has been incorporated into the model. As in the 1996 CDMG/U.S. Geological Survey (USGS) model, the time-dependent models incorporate three multisegment ruptures: a 1906, an 1857, and a southern San Andreas earthquake. Sensitivity tests are presented to show the effect on hazard and expected damage estimates of (1) intrinsic (aleatory) sigma, (2) multisegment (cascade) vs. independent segment (no cascade) ruptures, and (3) time-dependence vs. time-independence. Results indicate that (1) differences in hazard and expected damage estimates between time-dependent and independent models increase with decreasing intrinsic sigma, (2) differences in hazard and expected damage estimates between full cascading and not cascading are insensitive to intrinsic sigma, (3) differences in hazard increase with increasing return period (decreasing probability of occurrence), and (4) differences in moment-rate budgets increase with decreasing intrinsic sigma and with the degree of cascading, but are within the expected uncertainty in PSH time-dependent modeling and do not always significantly affect hazard and expected damage estimates.

California

Comparison of phase velocities from array measurements of Rayleigh waves associated with microtremor and results calculated from borehole shear-wave velocity profiles

Shear-wave velocities (V S ) are widely used for earthquake ground-motion site characterization. V S data are now largely obtained using borehole methods. Drilling holes, however, is expensive. Nonintrusive surface methods are inexpensive for obtaining V S information, but not many comparisons with direct borehole measurements have been published. Because different assumptions are used in data interpretation of each surface method and public safety is involved in site characterization for engineering structures, it is important to validate the surface methods by additional comparisons with borehole measurements. We compare results obtained from a particular surface method (array measurement of surface waves associated with microtremor) with results obtained from borehole methods. Using a 10-element nested-triangular array of 100-m aperture, we measured surface-wave phase velocities at two California sites, Garner Valley near Hemet and Hollister Municipal Airport. The Garner Valley site is located at an ancient lake bed where water-saturated sediment overlies decomposed granite on top of granite bedrock. Our array was deployed at a location where seismic velocities had been determined to a depth of 500 m by borehole methods. At Hollister, where the near-surface sediment consists of clay, sand, and gravel, we determined phase velocities using an array located close to a 60-m deep borehole where downhole velocity logs already exist. Because we want to assess the measurements uncomplicated by uncertainties introduced by the inversion process, we compare our phase-velocity results with the borehole V S depth profile by calculating fundamental-mode Rayleigh-wave phase velocities from an earth model constructed from the borehole data. For wavelengths less than ~2 times of the array aperture at Garner Valley, phase-velocity results from array measurements agree with the calculated Rayleigh-wave velocities to better than 11%. Measurement errors become larger for wavelengths 2 times greater than the array aperture. At Hollister, the measured phase velocity at 3.9 Hz (near the upper edge of the microtremor frequency band) is within 20% of the calculated Rayleigh-wave velocity. Because shear-wave velocity is the predominant factor controlling Rayleigh-wave phase velocities, the comparisons suggest that this nonintrusive method can provide V S information adequate for ground-motion estimation.

California

Discrepancy between earthquake rates implied by historic earthquakes and a consensus geologic source model for California

We examine the difference between expected earthquake rates inferred from the historical earthquake catalog and the geologic data that was used to develop the consensus seismic source characterization for the state of California [California Department of Conservation, Division of Mines and Geology (CDMG) and U.S. Geological Survey (USGS) Petersen et al., 1996; Frankel et al., 1996]. On average the historic earthquake catalog and the seismic source model both indicate about one M 6 or greater earthquake per year in the state of California. However, the overall earthquake rates of earthquakes with magnitudes (M) between 6 and 7 in this seismic source model are higher, by at least a factor of 2, than the mean historic earthquake rates for both southern and northern California. The earthquake rate discrepancy results from a seismic source model that includes earthquakes with characteristic (maximum) magnitudes that are primarily between M 6.4 and 7.1. Many of these faults are interpreted to accommodate high strain rates from geologic and geodetic data but have not ruptured in large earthquakes during historic time. Our sensitivity study indicates that the rate differences between magnitudes 6 and 7 can be reduced by adjusting the magnitude-frequency distribution of the source model to reflect more characteristic behavior, by decreasing the moment rate available for seismogenic slip along faults, by increasing the maximum magnitude of the earthquake on a fault, or by decreasing the maximum magnitude of the background seismicity. However, no single parameter can be adjusted, consistent with scientific consensus, to eliminate the earthquake rate discrepancy. Applying a combination of these parametric adjustments yields an alternative earthquake source model that is more compatible with the historic data. The 475-year return period hazard for peak ground and 1-sec spectral acceleration resulting from this alternative source model differs from the hazard resulting from the standard CDMG-USGS model by less than 10% across most of California but is higher (generally about 10% to 30%) within 20 km from some faults.

Bulletin of the Seismological Society of America

Correlation of ground motion and intensity for the 17 January 1994 Northridge, California, earthquake

We analyze the correlations between intensity and a set of groundmotion parameters obtained from 66 free-field stations in Los Angeles County that recorded the 1994 Northridge earthquake. We use the tagging intensities from Thywissen and Boatwright (1998) because these intensities are determined independently on census tracts, rather than interpolated from zip codes, as are the modified Mercalli isoseismals from Dewey et al. (1995). The ground-motion parameters we consider are the peak ground acceleration (PGA), the peak ground velocity (PGV), the 5% damped pseudovelocity response spectral (PSV) ordinates at 14 periods from 0.1 to 7.5 sec, and the rms average of these spectral ordinates from 0.3 to 3 sec. Visual comparisons of the distribution of tagging intensity with contours of PGA, PGV, and the average PSV suggest that PGV and the average PSV are better correlated with the intensity than PGA. The correlation coefficients between the intensity and the ground-motion parameters bear this out: r = 0.75 for PGA, 0.85 for PGV, and 0.85 for the average PSV. Correlations between the intensity and the PSV ordinates, as a function of period, are strongest at 1.5 sec (r = 0.83) and weakest at 0.2 sec (r = 0.66). Regressing the intensity on the logarithms of these ground-motion parameters yields relations I ?? mlog?? with 3.0 ??? m ??? 5.2 for the parameters analyzed, where m = 4.4 ?? 0.7 for PGA, 3.4 ?? 0.4 for PGV, and 3.6 ?? 0.5 for the average PSV.

Bulletin of the Seismological Society of America

The 1998 earthquake sequence south of Long Valley Caldera, California: Hints of magmatic involvement

A significant episode of seismic and geodetic unrest took place at Long Valley Caldera, California, beginning in the summer of 1997. Activity through late May of 1998 was concentrated in and around the south moat and the south margin of the resurgent dome. The Sierran Nevada block (SNB) region to the south/southeast remained relatively quiet until a M 5.1 event occurred there on 9 June 1998 (UT). A second M 5.1 event followed on 15 July (UT); both events were followed by appreciable aftershock sequences. An additional, distinct burst of activity began on 1 August 1998. The number of events in the August sequence (over the first week or two) was similar to the aftershock sequence of the 15 July 1998 M 5.1 event, but the later sequence was not associated with any events larger than M 4.3. All of the summer 1998 SNB activity was considered tectonic rather than magmatic; in general the SNB is considered an unlikely location for future eruptions. However, the August sequence-an 'aftershock sequence without a mainshock'- is suggestive of a strain event larger than the cumulative seismotectonic strain release. Moreover, a careful examination of waveforms from the August sequence reveals a small handful of events whose spectral signature is strikingly harmonic. We investigate the waveforms of these events using spectral, autocorrelation, and empirical Green's function techniques and conclude that they were most likely associated with a fluid-controlled source. Our observations suggest that there may have been some degree of magma or magma-derived fluid involvement in the 1998 SNB sequence.

California

High-resolution seismic reflection surveys and modeling across an area of high damage from the 1994 Northridge earthquake, Sherman Oaks, California

Approximately 3.6 km of P-wave seismic-reflection data were acquired along two orthogonal profiles in Sherman Oaks, California to determine whether shallow (less than 1-km depth) geologic structures contributed to the dramatic localized damage resulting from the 1994 Northridge earthquake. Both lines, one along Matilija Avenue and one along Milbank Street, crossed areas of both high and low damage. We believe these data reveal a geologic structure in the upper 600 m that contributed to the increased earthquake ground shaking in the high-damage areas south of and along the Los Angeles River. Of interest in these data is a reflection interpreted to be from bedrock that can be traced to the north along the Matilija Avenue profile. This reflecting interface, dipping northward at 15°–22°, may be an important impedance boundary because it is the lower boundary of a wedge of overlying low-velocity sediments. The wedge thins and terminates in the area where we interpret down-warped reflections as evidence of a shallow subbasin. The low-velocity subbasin sediments (V s of 200 m/sec V p of 500 m/sec) may be up to 150 m thick beneath the channelized Los Angeles River. The area across the subbasin experienced greater earthquake damage from possible geometric focusing effects. Three-dimensional basin effects may be responsible for the variable damage pattern, but from these seismic profiles it is not possible to determine the regional structural trends. Two-dimensional elastic and SH-mode finite-difference modeling of the imaged structural geometry along Matilija Avenue suggests that a peak horizontal-velocity amplification factor of two-and-over can be explained in the high-damage area above the shallow subbasin and sediment wedge. Amplification factors up to 5 were previously observed in aftershock data, at frequencies of 2 to 6 Hz. Amplification in the elastic simulation at the Santa Monica Mountains range-front on the southern end of the Matilija profile, with the geologic layering and geometry interpreted from the seismic data, is also consistent with aftershock observations.

Bulletin of the Seismological Society of America

A geometricla error in some Computer Programs based on the Aki-Christofferson-Husebye (ACH) Method of Teleseismic Tomography

Some computer programs based on the Aki-Christofferson-Husebye (ACH) method of teleseismic tomography contain an error caused by identifying local grid directions with azimuths on the spherical Earth. This error, which is most severe in high latitudes, introduces systematic errors into computed ray paths and distorts inferred Earth models. It is best dealt with by explicity correcting for the difference between true and grid directions. Methods for computing these directions are presented in this article and are likely to be useful in many other kinds of regional geophysical studies that use Cartesian coordinates and flat-earth approximations.

Bulletin of the Seismological Society of America

A Double-difference Earthquake location algorithm: Method and application to the Northern Hayward Fault, California

We have developed an efficient method to determine high-resolution hypocenter locations over large distances. The location method incorporates ordinary absolute travel-time measurements and/or cross-correlation P-and S-wave differential travel-time measurements. Residuals between observed and theoretical travel-time differences (or double-differences) are minimized for pairs of earthquakes at each station while linking together all observed event-station pairs. A least-squares solution is found by iteratively adjusting the vector difference between hypocentral pairs. The double-difference algorithm minimizes errors due to unmodeled velocity structure without the use of station corrections. Because catalog and cross-correlation data are combined into one system of equations, interevent distances within multiplets are determined to the accuracy of the cross-correlation data, while the relative locations between multiplets and uncorrelated events are simultaneously determined to the accuracy of the absolute travel-time data. Statistical resampling methods are used to estimate data accuracy and location errors. Uncertainties in double-difference locations are improved by more than an order of magnitude compared to catalog locations. The algorithm is tested, and its performance is demonstrated on two clusters of earthquakes located on the northern Hayward fault, California. There it colapses the diffuse catalog locations into sharp images of seismicity and reveals horizontal lineations of hypocenter that define the narrow regions on the fault where stress is released by brittle failure.

Bulletin of the Seismological Society of America

Geographic deaggregation of seismic hazard in the United States

The seismic hazard calculations for the 1996 national seismic hazard maps have been geographically deaggregated to assist in the understanding of the relative contributions of sources. These deaggregations are exhibited as maps with vertical bars whose heights are proportional to the contribution that each geographical cell makes to the ground-motion exceedance hazard. Bar colors correspond to average source magnitudes. We also extend the deaggregation analysis reported in Harmsen et al. (1999) to the western conterminous United States. In contrast to the central and eastern United States (CEUS); the influence of specific faults or characteristic events can be clearly identified. Geographic deaggregation for 0.2-sec and 1.0-sec pseudo spectral acceleration (SA) is performed for 10% probability of exceedance (PE) in 50 yr (475-yr mean return period) and 2% PE in 50 yr (2475-yr mean return period) for four western U.S. cities, Los Angeles, Salt Lake City, San Francisco, and Seattle, and for three central and eastern U.S. cities, Atlanta, Boston, and Saint Louis. In general, as the PE is lowered, the sources of hazard closer to the site dominate. Larger, more distant earthquakes contribute more significantly to hazard for 1.0-sec SA than for 0.2-sec SA. Additional maps of geographically deaggregated seismic hazard are available on the Internet for 120 cities in the conterminous United States (http://geohazards. cr.usgs.gov/eq/) for 1-sec SA and for 0.2-sec SA with a 2% PE in 50 yr. Examination of these maps of hazard contributions enables the investigator to determine the distance and azimuth to predominant sources, and their magnitudes. This information can be used to generate scenario earthquakes and corresponding time histories for seismic design and retrofit. Where fault density is lower than deaggregation cell dimensions, we can identify specific faults that contribute significantly to the seismic hazard at a given site. Detailed fault information enables investigators to include rupture information such as source directivity, radiation pattern, and basin-edge effects into their scenario earthquakes used in engineering analyses.

Bulletin of the Seismological Society of America

Correlation of 1- to 10-Hz earthquake resonances with surface measurements of S-wave reflections and refractions in the upper 50 m

Resonances observed in earthquake seismograms recorded in Seattle, Washington, the central United States and Sherman Oaks, California, are correlated with each site's respective near-surface seismic velocity profile and reflectivity determined from shallow seismic-reflection/refraction surveys. In all of these cases the resonance accounts for the highest amplitude shaking at the site above 1 Hz. These results show that imaging near-surface reflections from the ground surface can locate impedance structures that are important contributors to earthquake ground shaking. A high-amplitude S-wave reflection, recorded 250-m northeast and 300-m east of the Seattle Kingdome earthquake-recording station, with a two-way travel time of about 0.23 to 0.27 sec (about 18- to 22-m depth) marks the boundary between overlying alluvium (VS < 180 m/sec) and a higher velocity material (VS about 400 m/sec). This reflector probably causes a strong 2-Hz resonance that is observed in the earthquake data for the site near the Kingdome. In the central United States, S-wave reflections from a high-impedance boundary (an S-wave velocity increase from about 200 m/sec to 2000 m/sec) at about 40-m depth corresponds to a strong fundamental resonance at about 1.5 Hz. In Sherman Oaks, strong resonances at about 1.0 and 4 Hz are consistently observed on earthquake seismograms. A strong S-wave reflector at about 40-m depth may cause the 1.0 Hz resonance. The 4.0-Hz resonance is possibly explained by constructive interference between the first overtone of the 1.0-Hz resonance and a 3.25- to 3.9-Hz resonance calculated from an areally consistent impedance boundary at about 10-m depth as determined by S-wave refraction data.

Bulletin of the Seismological Society of America