USGS ScienceSearch

Geology topics

A.D. Frankel

Publications and source records attributed to A.D. Frankel.

At least 19 recordsLinked to original sources

Challenges in making a seismic hazard map for Alaska and the Aleutians

We present a summary of the data and analyses leading to the revision of the time-independent probabilistic seismic hazard maps of Alaska and the Aleutians. These maps represent a revision of existing maps based on newly obtained data, and reflect best current judgments about methodology and approach. They have been prepared following the procedures and assumptions made in the preparation of the 2002 National Seismic Hazard Maps for the lower 48 States, and will be proposed for adoption in future revisions to the International Building Code. We present example maps for peak ground acceleration, 0.2 s spectral amplitude (SA), and 1.0 s SA at a probability level of 2% in 50 years (annual probability of 0.000404). In this summary, we emphasize issues encountered in preparation of the maps that motivate or require future investigation and research.

Alaska

Fault displacement hazard for strike-slip faults

In this paper we present a methodology, data, and regression equations for calculating the fault rupture hazard at sites near steeply dipping, strike-slip faults. We collected and digitized on-fault and off-fault displacement data for 9 global strikeslip earthquakes ranging from moment magnitude M 6.5 to M 7.6 and supplemented these with displacements from 13 global earthquakes compiled byWesnousky (2008), who considers events up to M 7.9. Displacements on the primary fault fall off at the rupture ends and are often measured in meters, while displacements on secondary (offfault) or distributed faults may measure a few centimeters up to more than a meter and decay with distance from the rupture. Probability of earthquake rupture is less than 15% for cells 200 m??200 m and is less than 2% for 25 m??25 m cells at distances greater than 200mfrom the primary-fault rupture. Therefore, the hazard for off-fault ruptures is much lower than the hazard near the fault. Our data indicate that rupture displacements up to 35cm can be triggered on adjacent faults at distances out to 10kmor more from the primary-fault rupture. An example calculation shows that, for an active fault which has repeated large earthquakes every few hundred years, fault rupture hazard analysis should be an important consideration in the design of structures or lifelines that are located near the principal fault, within about 150 m of well-mapped active faults with a simple trace and within 300 m of faults with poorly defined or complex traces.

Bulletin of the Seismological Society of America

Uniform California earthquake rupture forecast, version 2 (UCERF 2)

The 2007 Working Group on California Earthquake Probabilities (WGCEP, 2007) presents the Uniform California Earthquake Rupture Forecast, Version 2 (UCERF 2). This model comprises a time-independent (Poisson-process) earthquake rate model, developed jointly with the National Seismic Hazard Mapping Program and a time-dependent earthquake-probability model, based on recent earthquake rates and stress-renewal statistics conditioned on the date of last event. The models were developed from updated statewide earthquake catalogs and fault deformation databases using a uniform methodology across all regions and implemented in the modular, extensible Open Seismic Hazard Analysis framework. The rate model satisfies integrating measures of deformation across the plate-boundary zone and is consistent with historical seismicity data. An overprediction of earthquake rates found at intermediate magnitudes (6.5 ??? M ???7.0) in previous models has been reduced to within the 95% confidence bounds of the historical earthquake catalog. A logic tree with 480 branches represents the epistemic uncertainties of the full time-dependent model. The mean UCERF 2 time-dependent probability of one or more M ???6.7 earthquakes in the California region during the next 30 yr is 99.7%; this probability decreases to 46% for M ???7.5 and to 4.5% for M ???8.0. These probabilities do not include the Cascadia subduction zone, largely north of California, for which the estimated 30 yr, M ???8.0 time-dependent probability is 10%. The M ???6.7 probabilities on major strike-slip faults are consistent with the WGCEP (2003) study in the San Francisco Bay Area and the WGCEP (1995) study in southern California, except for significantly lower estimates along the San Jacinto and Elsinore faults, owing to provisions for larger multisegment ruptures. Important model limitations are discussed.

Bulletin of the Seismological Society of America

Earthquake hazard in the New Madrid Seismic Zone remains a concern

There is broad agreement in the scientific community that a continuing concern exists for a major destructive earthquake in the New Madrid seismic zone. Many structures in Memphis, Tenn., St. Louis, Mo., and other communities in the central Mississippi River Valley region are vulnerable and at risk from severe ground shaking. This assessment is based on decades of research on New Madrid earthquakes and related phenomena by dozens of Federal, university, State, and consulting earth scientists. Considerable interest has developed recently from media reports that the New Madrid seismic zone may be shutting down. These reports stem from published research using global positioning system (GPS) instruments with results of geodetic measurements of strain in the Earth's crust. Because of a lack of measurable strain at the surface in some areas of the seismic zone over the past 14 years, arguments have been advanced that there is no buildup of stress at depth within the New Madrid seismic zone and that the zone may no longer pose a significant hazard. As part of the consensus-building process used to develop the national seismic hazard maps, the U.S. Geological Survey (USGS) convened a workshop of experts in 2006 to evaluate the latest findings in earthquake hazards in the Eastern United States. These experts considered the GPS data from New Madrid available at that time that also showed little to no ground movement at the surface. The experts did not find the GPS data to be a convincing reason to lower the assessment of earthquake hazard in the New Madrid region, especially in light of the many other types of data that are used to construct the hazard assessment, several of which are described here.

Arkansas, Illinois, Kentucky, Missouri, Tennessee

Model uncertainties of the 2002 update of California seismic hazard maps

In this article we present and explore the source and ground-motion model uncertainty and parametric sensitivity for the 2002 update of the California probabilistic seismic hazard maps. Our approach is to implement a Monte Carlo simulation that allows for independent sampling from fault to fault in each simulation. The source-distance dependent characteristics of the uncertainty maps of seismic hazard are explained by the fundamental uncertainty patterns from four basic test cases, in which the uncertainties from one-fault and two-fault systems are studied in detail. The California coefficient of variation (COV, ratio of the standard deviation to the mean) map for peak ground acceleration (10% of exceedance in 50 years) shows lower values (0.1-0.15) along the San Andreas fault system and other class A faults than along class B faults (0.2-0.3). High COV values (0.4-0.6) are found around the Garlock, Anacapa-Dume, and Palos Verdes faults in southern California and around the Maacama fault and Cascadia subduction zone in northern California.

Bulletin of the Seismological Society of America

Probabilistic seismic hazard analysis for Sumatra, Indonesia and across the Southern Malaysian Peninsula

The ground motion hazard for Sumatra and the Malaysian peninsula is calculated in a probabilistic framework, using procedures developed for the US National Seismic Hazard Maps. We constructed regional earthquake source models and used standard published and modified attenuation equations to calculate peak ground acceleration at 2% and 10% probability of exceedance in 50 years for rock site conditions. We developed or modified earthquake catalogs and declustered these catalogs to include only independent earthquakes. The resulting catalogs were used to define four source zones that characterize earthquakes in four tectonic environments: subduction zone interface earthquakes, subduction zone deep intraslab earthquakes, strike-slip transform earthquakes, and intraplate earthquakes. The recurrence rates and sizes of historical earthquakes on known faults and across zones were also determined from this modified catalog. In addition to the source zones, our seismic source model considers two major faults that are known historically to generate large earthquakes: the Sumatran subduction zone and the Sumatran transform fault. Several published studies were used to describe earthquakes along these faults during historical and pre-historical time, as well as to identify segmentation models of faults. Peak horizontal ground accelerations were calculated using ground motion prediction relations that were developed from seismic data obtained from the crustal interplate environment, crustal intraplate environment, along the subduction zone interface, and from deep intraslab earthquakes. Most of these relations, however, have not been developed for large distances that are needed for calculating the hazard across the Malaysian peninsula, and none were developed for earthquake ground motions generated in an interplate tectonic environment that are propagated into an intraplate tectonic environment. For the interplate and intraplate crustal earthquakes, we have applied ground-motion prediction relations that are consistent with California (interplate) and India (intraplate) strong motion data that we collected for distances beyond 200 km. For the subduction zone equations, we recognized that the published relationships at large distances were not consistent with global earthquake data that we collected and modified the relations to be compatible with the global subduction zone ground motions. In this analysis, we have used alternative source and attenuation models and weighted them to account for our uncertainty in which model is most appropriate for Sumatra or for the Malaysian peninsula. The resulting peak horizontal ground accelerations for 2% probability of exceedance in 50 years range from over 100% g to about 10% g across Sumatra and generally less than 20% g across most of the Malaysian peninsula. The ground motions at 10% probability of exceedance in 50 years are typically about 60% of the ground motions derived for a hazard level at 2% probability of exceedance in 50 years. The largest contributors to hazard are from the Sumatran faults.

Tectonophysics

The 2002 Denali fault earthquake, Alaska: A large magnitude, slip-partitioned event

The MW (moment magnitude) 7.9 Denali fault earthquake on 3 November 2002 was associated with 340 kilometers of surface rupture and was the largest strike-slip earthquake in North America in almost 150 years. It illuminates earthquake mechanics and hazards of large strike-slip faults. It began with thrusting on the previously unrecognized Susitna Glacier fault, continued with right-slip on the Denali fault, then took a right step and continued with right-slip on the Totschunda fault. There is good correlation between geologically observed and geophysically inferred moment release. The earthquake produced unusually strong distal effects in the rupture propagation direction, including triggered seismicity.

Science

Simulations of seismic hazard for the Pacific Northwest of the United States from earthquakes associated with the Cascadia subduction zone

We investigate the impact of different rupture and attenuation models for the Cascadia subduction zone by simulating seismic hazard models for the Pacific Northwest of the U.S. at 2% probability of exceedance in 50 years. We calculate the sensitivity of hazard (probabilistic ground motions) to the source parameters and the attenuation relations for both intraslab and interface earthquakes and present these in the framework of the standard USGS hazard model that includes crustal earthquakes. Our results indicate that allowing the deep intraslab earthquakes to occur anywhere along the subduction zone increases the peak ground acceleration hazard near Portland, Oregon by about 20%. Alternative attenuation relations for deep earthquakes can result in ground motions that differ by a factor of two. The hazard uncertainty for the plate interface and intraslab earthquakes is analyzed through a Monte-Carlo logic tree approach and indicates a seismic hazard exceeding 1 g (0.2 s spectral acceleration) consistent with the U.S. National Seismic Hazard Maps in western Washington, Oregon, and California and an overall coefficient of variation that ranges from 0.1 to 0.4. Sensitivity studies indicate that the paleoseismic chronology and the magnitude of great plate interface earthquakes contribute significantly to the hazard uncertainty estimates for this region. Paleoseismic data indicate that the mean earthquake recurrence interval for great earthquakes is about 500 years and that it has been 300 years since the last great earthquake. We calculate the probability of such a great earthquake along the Cascadia plate interface to be about 14% when considering a time-dependent model and about 10% when considering a time-independent Poisson model during the next 50-year interval.

Pure and Applied Geophysics

Nonlinear and linear site response and basin effects in Seattle for the M 6.8 Nisqually, Washington, earthquake

We used recordings of the M 6.8 Nisqually earthquake and its ML 3.4 aftershock to study site response and basin effects for 35 locations in Seattle, Washington. We determined site amplification from Fourier spectral ratios of the recorded horizontal ground motions, referenced to a soft-rock site. Soft-soil sites (generally National Earthquake Hazard Reduction Program [NEHRP] class E) on artificial fill and young alluvium have the largest 1-Hz amplifications (factors of 3-7) for both the mainshock and aftershock. These amplifications are correlated with areas of higher damage from the mainshock to major buildings and liquefaction. There are several indications of nonlinear response at the soft-soil sites for the mainshock ground motions, despite relatively modest peak accelerations in the S waves of 15%-22%g. First, the mainshock spectral ratios do not show amplification at 2-8 Hz as do the aftershock spectral ratios. Spectral peaks at frequencies below 2 Hz generally occur at lower frequencies for the mainshock spectral ratios than for the aftershock ratios. At one soft-soil site, there is a clear shift of the resonant frequency to a lower frequency for the mainshock compared with the aftershock. The frequency of this resonance increases in the coda of the mainshock record, indicating that the site response during the weaker motions of the coda is more linear than that of the initial S wave. Three of the soft-soil sites display cusped, one-sided mainshock accelerograms after the S wave. These soft-soil sites also show amplification at 10-20 Hz in the S wave, relative to the rock site, that is not observed for the aftershock. The cusped waveforms and 10-20-Hz amplification are symptomatic of nonlinear response at the soft-soil sites. These sites had nearby liquefaction. The largest amplifications for 0.5 Hz occur at soft-soil sites on the southern portion of the Seattle Basin. Stiff-soil sites (NEHRP classes D and C) on Pleistocene-age glacial deposits display similar spectral amplification for the mainshock and aftershock, indicating approximately linear response. The stiff-soil sites generally have moderate amplification (factors of 1.1-2.4) at 0.5 and 1 Hz. Amplifications at 1 and 5 Hz for all sites generally increase with decreasing shear-wave velocity measured in the top 30 m (Vs 30). However, larger amplifications at 0.5 and 1 Hz for sites with similar Vs 30 values are observed for sites in the Seattle Basin, illustrating the amplification from the deeper (>30 m) sediments and the contribution from basin surface waves. Record sections for the mainshock and aftershock show that basin surface waves produce the peak velocities for many of the sites in the Seattle Basin and often dominate the amplitude at 1 Hz and lower frequencies.

Bulletin of the Seismological Society of America

Seismic hazard in Hawaii: High rate of large earthquakes and probabilistics ground-motion maps

The seismic hazard and earthquake occurrence rates in Hawaii are locally as high as that near the most hazardous faults elsewhere in the United States. We have generated maps of peak ground acceleration (PGA) and spectral acceleration (SA) (at 0.2, 0.3 and 1.0 sec, 5% critical damping) at 2% and 10% exceedance probabilities in 50 years. The highest hazard is on the south side of Hawaii Island, as indicated by the MI 7.0, MS 7.2, and MI 7.9 earthquakes, which occurred there since 1868. Probabilistic values of horizontal PGA (2% in 50 years) on Hawaii's south coast exceed 1.75g. Because some large earthquake aftershock zones and the geometry of flank blocks slipping on subhorizontal decollement faults are known, we use a combination of spatially uniform sources in active flank blocks and smoothed seismicity in other areas to model seismicity. Rates of earthquakes are derived from magnitude distributions of the modem (1959-1997) catalog of the Hawaiian Volcano Observatory's seismic network supplemented by the historic (1868-1959) catalog. Modern magnitudes are ML measured on a Wood-Anderson seismograph or MS. Historic magnitudes may add ML measured on a Milne-Shaw or Bosch-Omori seismograph or MI derived from calibrated areas of MM intensities. Active flank areas, which by far account for the highest hazard, are characterized by distributions with b slopes of about 1.0 below M 5.0 and about 0.6 above M 5.0. The kinked distribution means that large earthquake rates would be grossly under-estimated by extrapolating small earthquake rates, and that longer catalogs are essential for estimating or verifying the rates of large earthquakes. Flank earthquakes thus follow a semicharacteristic model, which is a combination of background seismicity and an excess number of large earthquakes. Flank earthquakes are geometrically confined to rupture zones on the volcano flanks by barriers such as rift zones and the seaward edge of the volcano, which may be expressed by a magnitude distribution similar to that including characteristic earthquakes. The island chain northwest of Hawaii Island is seismically and volcanically much less active. We model its seismic hazard with a combination of a linearly decaying ramp fit to the cataloged seismicity and spatially smoothed seismicity with a smoothing half-width of 10 km. We use a combination of up to four attenuation relations for each map because for either PGA or SA, there is no single relation that represents ground motion for all distance and magnitude ranges. Great slumps and landslides visible on the ocean floor correspond to catastrophes with effective energy magnitudes ME above 8.0. A crude estimate of their frequency suggests that the probabilistic earthquake hazard is at least an order of magnitude higher for flank earthquakes than that from submarine slumps.

Bulletin of the Seismological Society of America

Configuration of the Seattle Urban Seismic Array for the February 28, 2001, M6.8, Nisqually earthquake and its aftershocks

This is a description of the configuration of the Seattle Urban Seismic Array (click on thumbnail map to see full size map) at the time of the Nisqually earthquake and its aftershocks (from April 28 through March 3, 2001). The array is operated by the Geologic Hazards Team of the U. S Geological Survey. It has been recording small local earthquakes since its initial installation immediately after the Duvall earthquake in May 1996. Since then, we have been able to convert the array to all Kinemetrics K2s after retiring the aging DR-200s in June of 1998. The array stations have been moved around to suit our research interests. We had 26 instruments in urban Seattle at the time of the Nisqually earthquake. The primary research goal of the array is to record earthquake ground motion at sites scattered around the city that represent the major types of surficial geologic materials found in the city. With this information we plan to draw seismic hazard maps that will show the expected variability of ground shaking throughout the city. Our strategy to achieve our goal includes: 1) sampling ground motion at sites that represent the major geologic materials including man-made fill and alluvial surficial materials, 2) recording high-resolution ground motions from a tripartite array (SDN, SDS, and SDW) with about 200 meter spacing located on man-made fill in the SoDo (South of Downtown) industrial district, 3) recording ground motion in areas with man-made fill and alluvial surficial materials which have a high potential for liquifaction, 4) sampling ground motion at two sites on Capitol Hill (THO, and SEU) that are located near drill holes that have been logged for seismic velocity, 5) Providing long-term continuity at a few sites for comparison and reference, especially SEW and ALK, which are located on Tertiary bedrock. 6) recording ground motion at sites with similar surficial geology both inside and outside the Seattle Basin The configuration of the array is summarized in the accompanying table which is further described below.

Washington

Seismic design parameters - A user guide

The 1997 NEHRP Recommended Provisions for Seismic Regulations for New Buildings (1997 NEHRP Provisions ) introduced seismic design procedure that is based on the explicit use of spectral response acceleration rather than the traditional peak ground acceleration and/or peak ground velocity or zone factors. The spectral response accelerations are obtained from spectral response acceleration maps accompanying the report. Maps are available for the United States and a number of U.S. territories. Since 1997 additional codes and standards have also adopted seismic design approaches based on the same procedure used in the NEHRP Provisions and the accompanying maps. The design documents using the 1997 NEHRP Provisions procedure may be divided into three categories -(1) Design of New Construction, (2) Design and Evaluation of Existing Construction, and (3) Design of Residential Construction. A CD-ROM has been prepared for use in conjunction with the design documents in each of these three categories. The spectral accelerations obtained using the software on the CD are the same as those that would be obtained by using the maps accompanying the design documents. The software has been prepared to operate on a personal computer using a Windows (Microsoft Corporation) operating environment and a point and click type of interface. The user can obtain the spectral acceleration values that would be obtained by use of the maps accompanying the design documents, include site factors appropriate for the Site Class provided by the user, calculate a response spectrum that includes the site factor, and plot a response spectrum. Sites may be located by providing the latitude-longitude or zip code for all areas covered by the maps. All of the maps used in the various documents are also included on the CDROM

Open-File Report

Development of Maximum Considered Earthquake Ground Motion Maps

The 1997 NEHRP Recommended Provisions for Seismic Regulations for New Buildings use a design procedure that is based on spectral response acceleration rather than the traditional peak ground acceleration, peak ground velocity, or zone factors. The spectral response accelerations are obtained from maps prepared following the recommendations of the Building Seismic Safety Council's (BSSC) Seismic Design Procedures Group (SDPG). The SDPG-recommended maps, the Maximum Considered Earthquake (MCE) Ground Motion Maps, are based on the U.S. Geological Survey (USGS) probabilistic hazard maps with additional modifications incorporating deterministic ground motions in selected areas and the application of engineering judgement. The MCE ground motion maps included with the 1997 NEHRP Provisions also serve as the basis for the ground motion maps used in the seismic design portions of the 2000 International Building Code and the 2000 International Residential Code . Additionally the design maps prepared for the 1997 NEHRP Provisions , combined with selected USGS probabilistic maps, are used with the 1997 NEHRP Guidelines for the Seismic Rehabilitation of Buildings .

Earthquake Spectra