Geology topics

David M. Boore

Publications and source records attributed to David M. Boore.

At least 19 recordsLinked to original sources

USGS2016

NGA-West2 equations for predicting vertical-component PGA, PGV, and 5%-damped PSA from shallow crustal earthquakes

We present ground motion prediction equations (GMPEs) for computing natural log means and standard deviations of vertical-component intensity measures (IMs) for shallow crustal earthquakes in active tectonic regions. The equations were derived from a global database with M 3.0–7.9 events. The functions are similar to those for our horizontal GMPEs. We derive equations for the primary M - and distance-dependence of peak acceleration, peak velocity, and 5%-damped pseudo-spectral accelerations at oscillator periods between 0.01–10 s. We observe pronounced M -dependent geometric spreading and region-dependent anelastic attenuation for high-frequency IMs. We do not observe significant region-dependence in site amplification. Aleatory uncertainty is found to decrease with increasing magnitude; within-event variability is independent of distance. Compared to our horizontal-component GMPEs, attenuation rates are broadly comparable (somewhat slower geometric spreading, faster apparent anelastic attenuation), V S 30 -scaling is reduced, nonlinear site response is much weaker, within-event variability is comparable, and between-event variability is greater.

Earthquake Spectra↗

USGS2014

A generalization of the double-corner-frequency source spectral model and its use in the SCEC BBP validation exercise

The stochastic method of simulating ground motions requires the specification of the shape and scaling with magnitude of the source spectrum. The spectral models commonly used are either single-corner-frequency or double-corner-frequency models, but the latter have no flexibility to vary the high-frequency spectral levels for a specified seismic moment. Two generalized double-corner-frequency ω 2 source spectral models are introduced, one in which two spectra are multiplied together, and another where they are added. Both models have a low-frequency dependence controlled by the seismic moment, and a high-frequency spectral level controlled by the seismic moment and a stress parameter. A wide range of spectral shapes can be obtained from these generalized spectral models, which makes them suitable for inversions of data to obtain spectral models that can be used in ground-motion simulations in situations where adequate data are not available for purely empirical determinations of ground motions, as in stable continental regions. As an example of the use of the generalized source spectral models, data from up to 40 stations from seven events, plus response spectra at two distances and two magnitudes from recent ground-motion prediction equations, were inverted to obtain the parameters controlling the spectral shapes, as well as a finite-fault factor that is used in point-source, stochastic-method simulations of ground motion. The fits to the data are comparable to or even better than those from finite-fault simulations, even for sites close to large earthquakes.

Bulletin of the Seismological Society of America↗

USGS2014

Path durations for use in the stochastic‐method simulation of ground motions

The stochastic method of ground‐motion simulation assumes that the energy in a target spectrum is spread over a duration D T . D T is generally decomposed into the duration due to source effects (D S ) and to path effects (D P ). For the most commonly used source, seismological theory directly relates D S to the source corner frequency, accounting for the magnitude scaling of D T . In contrast, D P is related to propagation effects that are more difficult to represent by analytic equations based on the physics of the process. We are primarily motivated to revisit D T because the function currently employed by many implementations of the stochastic method for active tectonic regions underpredicts observed durations, leading to an overprediction of ground motions for a given target spectrum. Further, there is some inconsistency in the literature regarding which empirical duration corresponds to D T . Thus, we begin by clarifying the relationship between empirical durations and D T as used in the first author’s implementation of the stochastic method, and then we develop a new D P relationship. The new D P function gives significantly longer durations than in the previous D P function, but the relative contribution of D P to D T still diminishes with increasing magnitude. Thus, this correction is more important for small events or subfaults of larger events modeled with the stochastic finite‐fault method.

Bulletin of the Seismological Society of America↗

USGS2014

NGA-West 2 Equations for predicting PGA, PGV, and 5%-Damped PSA for shallow crustal earthquakes

We provide ground-motion prediction equations for computing medians and standard deviations of average horizontal component intensity measures (IMs) for shallow crustal earthquakes in active tectonic regions. The equations were derived from a global database with M 3.0–7.9 events. We derived equations for the primary M- and distance-dependence of the IMs after fixing the V S30 -based nonlinear site term from a parallel NGA-West 2 study. We then evaluated additional effects using mixed effects residuals analysis, which revealed no trends with source depth over the M range of interest, indistinct Class 1 and 2 event IMs, and basin depth effects that increase and decrease long-period IMs for depths larger and smaller, respectively, than means from regional V S30 -depth relations. Our aleatory variability model captures decreasing between-event variability with M, as well as within-event variability that increases or decreases with M depending on period, increases with distance, and decreases for soft sites.

Earthquake Spectra↗

USGS2014

NGA-West2 Research Project

The NGA-West2 project is a large multidisciplinary, multi-year research program on the Next Generation Attenuation (NGA) models for shallow crustal earthquakes in active tectonic regions. The research project has been coordinated by the Pacific Earthquake Engineering Research Center (PEER), with extensive technical interactions among many individuals and organizations. NGA-West2 addresses several key issues in ground-motion seismic hazard, including updating the NGA database for a magnitude range of 3.0–7.9; updating NGA ground-motion prediction equations (GMPEs) for the “average” horizontal component; scaling response spectra for damping values other than 5%; quantifying the effects of directivity and directionality for horizontal ground motion; resolving discrepancies between the NGA and the National Earthquake Hazards Reduction Program (NEHRP) site amplification factors; analysis of epistemic uncertainty for NGA GMPEs; and developing GMPEs for vertical ground motion. This paper presents an overview of the NGA-West2 research program and its subprojects.

Earthquake Spectra↗

USGS2013

The effect of sampling rate and anti-aliasing filters on high-frequency response spectra

The most commonly used intensity measure in ground-motion prediction equations is the pseudo-absolute response spectral acceleration (PSA), for response periods from 0.01 to 10 s (or frequencies from 0.1 to 100 Hz). PSAs are often derived from recorded ground motions, and these motions are usually filtered to remove high and low frequencies before the PSAs are computed. In this article we are only concerned with the removal of high frequencies. In modern digital recordings, this filtering corresponds at least to an anti-aliasing filter applied before conversion to digital values. Additional high-cut filtering is sometimes applied both to digital and to analog records to reduce high-frequency noise. Potential errors on the short-period (high-frequency) response spectral values are expected if the true ground motion has significant energy at frequencies above that of the anti-aliasing filter. This is especially important for areas where the instrumental sample rate and the associated anti-aliasing filter corner frequency (above which significant energy in the time series is removed) are low relative to the frequencies contained in the true ground motions. A ground-motion simulation study was conducted to investigate these effects and to develop guidance for defining the usable bandwidth for high-frequency PSA. The primary conclusion is that if the ratio of the maximum Fourier acceleration spectrum (FAS) to the FAS at a frequency f saa corresponding to the start of the anti-aliasing filter is more than about 10, then PSA for frequencies above f saa should be little affected by the recording process, because the ground-motion frequencies that control the response spectra will be less than f saa . A second topic of this article concerns the resampling of the digital acceleration time series to a higher sample rate often used in the computation of short-period PSA. We confirm previous findings that sinc-function interpolation is preferred to the standard practice of using linear time interpolation for the resamplin

Bulletin of Earthquake Engineering↗

USGS2013

What do data used to develop ground-motion prediction equations tell us about motions near faults?

A large database of ground motions from shallow earthquakes occurring in active tectonic regions around the world, recently developed in the Pacific Earthquake Engineering Center’s NGA-West2 project, has been used to investigate what such a database can say about the properties and processes of crustal fault zones. There are a relatively small number of near-rupture records, implying that few recordings in the database are within crustal fault zones, but the records that do exist emphasize the complexity of ground-motion amplitudes and polarization close to individual faults. On average over the whole data set, however, the scaling of ground motions with magnitude at a fixed distance, and the distance dependence of the ground motions, seem to be largely consistent with simple seismological models of source scaling, path propagation effects, and local site amplification. The data show that ground motions close to large faults, as measured by elastic response spectra, tend to saturate and become essentially constant for short periods. This saturation seems to be primarily a geometrical effect, due to the increasing size of the rupture surface with magnitude, and not due to a breakdown in self similarity.

Pure and Applied Geophysics↗

USGS2013

The attenuation of Fourier amplitudes for rock sites in eastern North America

We develop an empirical model of the decay of Fourier amplitudes for earthquakes of M 3–6 recorded on rock sites in eastern North America and discuss its implications for source parameters. Attenuation at distances from 10 to 500 km may be adequately described using a bilinear model with a geometric spreading of 1/R 1.3 to a transition distance of 50 km, with a geometric spreading of 1/R 0.5 at greater distances. For low frequencies and distances less than 50 km, the effective geometric spreading given by the model is perturbed using a frequency‐ and hypocentral depth‐dependent factor defined in such a way as to increase amplitudes at lower frequencies near the epicenter but leave the 1 km source amplitudes unchanged. The associated anelastic attenuation is determined for each event, with an average value being given by a regional quality factor of Q=525f 0.45 . This model provides a match, on average, between the known seismic moment of events and the inferred low‐frequency spectral amplitudes at R=1 km (obtained by correcting for the attenuation model). The inferred Brune stress parameters from the high‐frequency source terms are about 600 bars (60 MPa), on average, for events of M>4.5.

Bulletin of the Seismological Society of America↗

USGS2013

Extreme ground motions and Yucca Mountain

Yucca Mountain is the designated site of the underground repository for the United States' high-level radioactive waste (HLW), consisting of commercial and military spent nuclear fuel, HLW derived from reprocessing of uranium and plutonium, surplus plutonium, and other nuclear-weapons materials. Yucca Mountain straddles the western boundary of the Nevada Test Site, where the United States has tested nuclear devices since the 1950s, and is situated in an arid, remote, and thinly populated region of Nevada, ~100 miles northwest of Las Vegas. Yucca Mountain was originally considered as a potential underground repository of HLW because of its thick units of unsaturated rocks, with the repository horizon being not only ~300 m above the water table but also ~300 m below the Yucca Mountain crest. The fundamental rationale for a geologic (underground) repository for HLW is to securely isolate these materials from the environment and its inhabitants to the greatest extent possible and for very long periods of time. Given the present climate conditions and what is known about the current hydrologic system and conditions around and in the mountain itself, one would anticipate that the rates of infiltration, corrosion, and transport would be very low—except for the possibility that repository integrity might be compromised by low-probability disruptive events, which include earthquakes, strong ground motion, and (or) a repository-piercing volcanic intrusion/eruption. Extreme ground motions (ExGM), as we use the phrase in this report, refer to the extremely large amplitudes of earthquake ground motion that arise at extremely low probabilities of exceedance (hazard). They first came to our attention when the 1998 probabilistic seismic hazard analysis for Yucca Mountain was extended to a hazard level of 10 -8 /yr (a 10-4/yr probability for a 10 4 -year repository “lifetime”). The primary purpose of this report is to summarize the principal results of the ExGM research program as they have developed over the past 5 years; what follows will be focused on Yucca Mountain, but not restricted to it.

Nevada↗

USGS2012

Predominant-period site classification for response spectra prediction equations in Italy

We propose a site‐classification scheme based on the predominant period of the site, as determined from the average horizontal‐to‐vertical (H/V) spectral ratios of ground motion. Our scheme extends Zhao et al. (2006) classifications by adding two classes, the most important of which is defined by flat H/V ratios with amplitudes less than 2. The proposed classification is investigated by using 5%‐damped response spectra from Italian earthquake records. We select a dataset of 602 three‐component analog and digital recordings from 120 earthquakes recorded at 214 seismic stations within a hypocentral distance of 200 km. Selected events are in the moment‐magnitude range 4.0≤M w ≤6.8 and focal depths from a few kilometers to 46 km. We computed H/V ratios for these data and used them to classify each site into one of six classes. We then investigate the impact of this classification scheme on empirical ground‐motion prediction equations (GMPEs) by comparing its performance with that of the conventional rock/soil classification. Although the adopted approach results in only a small reduction of the overall standard deviation, the use of H/V spectral ratios in site classification does capture the signature of sites with flat frequency‐response, as well as deep and shallow‐soil profiles, characterized by long‐ and short‐period resonance, respectively; in addition, the classification scheme is relatively quick and inexpensive, which is an advantage over schemes based on measurements of shear‐wave velocity.

BSSA↗

USGS2011

Regional correlations of V S 30 averaged over depths less than and greater than 30 meters

Using velocity profiles from sites in Japan, California, Turkey, and Europe, we find that the time-averaged shear-wave velocity to 30 m ( V S 30 ), used as a proxy for site amplification in recent ground-motion prediction equations (GMPEs) and building codes, is strongly correlated with average velocities to depths less than 30 m ( V Sz , with z being the averaging depth). The correlations for sites in Japan (corresponding to the KiK-net network) show that V Sz is systematically larger for a given V Sz than for profiles from the other regions. The difference largely results from the placement of the KiK-net station locations on rock and rocklike sites, whereas stations in the other regions are generally placed in urban areas underlain by sediments. Using the KiK-net velocity profiles, we provide equations relating V S 30 to V Sz for z ranging from 5 to 29 m in 1-m increments. These equations (and those for California velocity profiles given in Boore, 2004b) can be used to estimate V S 30 from V Sz for sites in which velocity profiles do not extend to 30 m. The scatter of the residuals decreases with depth, but, even for an averaging depth of 5 m, a variation in log V S 30 of ±1 standard deviation maps into less than a 20% uncertainty in ground motions given by recent GMPEs at short periods. The sensitivity of the ground motions to V S 30 uncertainty is considerably larger at long periods (but is less than a factor of 1.2 for averaging depths greater than about 20 m). We also find that V S 30 is correlated with V Sz for z as great as 400 m for sites of the KiK-net network, providing some justification for using V S 30 as a site-response variable for predicting ground motions at periods for which the wavelengths far exceed 30 m.

Bulletin of the Seismological Society of America↗

USGS2011

Estimating unknown input parameters when implementing the NGA ground-motion prediction equations in engineering practice

The ground-motion prediction equations (GMPEs) developed as part of the Next Generation Attenuation of Ground Motions (NGA-West) project in 2008 are becoming widely used in seismic hazard analyses. However, these new models are considerably more complicated than previous GMPEs, and they require several more input parameters. When employing the NGA models, users routinely face situations in which some of the required input parameters are unknown. In this paper, we present a framework for estimating the unknown source, path, and site parameters when implementing the NGA models in engineering practice, and we derive geometrically-based equations relating the three distance measures found in the NGA models. Our intent is for the content of this paper not only to make the NGA models more accessible, but also to help with the implementation of other present or future GMPEs.

Earthquake Spectra↗

USGS2011

Representation of bidirectional ground motions for design spectra in building codes

The 2009 NEHRP Provisions modified the definition of horizontal ground motion from the geometric mean of spectral accelerations for two components to the peak response of a single lumped mass oscillator regardless of direction. These maximum-direction (MD) ground motions operate under the assumption that the dynamic properties of the structure (e.g., stiffness, strength) are identical in all directions. This assumption may be true for some in-plan symmetric structures, however, the response of most structures is dominated by modes of vibration along specific axes (e.g., longitudinal and transverse axes in a building), and often the dynamic properties (especially stiffness) along those axes are distinct. In order to achieve structural designs consistent with the collapse risk level given in the NEHRP documents, we argue that design spectra should be compatible with expected levels of ground motion along those principal response axes. The use of MD ground motions effectively assumes that the azimuth of maximum ground motion coincides with the directions of principal structural response. Because this is unlikely, design ground motions have lower probability of occurrence than intended, with significant societal costs. We recommend adjustments to make design ground motions compatible with target risk levels.

Earthquake Spectra↗

USGS2010

Orientation-independent, nongeometric-mean measures of seismic intensity from two horizontal components of motion

New measures of spectral intensity based on the horizontal components of ground shaking are introduced. These new measures are independent of the in situ orientation of the recordings and encompass the full range of spectral amplitudes over all possible rotation angles. Unlike previously introduced measures that are also orientation independent, no geometric means are used in the computation of the new measures. The new measures based on fiftieth percentile values of the response spectra show small but systematic increases (to a factor of about 1.07 at a 10 sec period) compared to the comparable geometric-mean measure.

Bulletin of the Seismological Society of America↗

USGS2010

Implementation of the Next Generation Attenuation (NGA) ground-motion prediction equations in Fortran and R

This report presents two methods for implementing the earthquake ground-motion prediction equations released in 2008 as part of the Next Generation Attenuation of Ground Motions (NGA-West, or NGA) project coordinated by the Pacific Earthquake Engineering Research Center (PEER). These models were developed for predicting ground-motion parameters for shallow crustal earthquakes in active tectonic regions (such as California). Of the five ground-motion prediction equations (GMPEs) developed during the NGA project, four models are implemented: the GMPEs of Abrahamson and Silva (2008), Boore and Atkinson (2008), Campbell and Bozorgnia (2008), and Chiou and Youngs (2008a); these models are abbreviated as AS08, BA08, CB08, and CY08, respectively. Since site response is widely recognized as an important influence of ground motions, engineering applications typically require that such effects be modeled. The model of Idriss (2008) is not implemented in our programs because it does not explicitly include site response, whereas the other four models include site response and use the same variable to describe the site condition (VS30). We do not intend to discourage the use of the Idriss (2008) model, but we have chosen to implement the other four NGA models in our programs for those users who require ground-motion estimates for various site conditions. We have implemented the NGA models by using two separate programming languages: Fortran and R (R Development Core Team, 2010). Fortran, a compiled programming language, has been used in the scientific community for decades. R is an object-oriented language and environment for statistical computing that is gaining popularity in the statistical and scientific community. Derived from the S language and environment developed at Bell Laboratories, R is an open-source language that is freely available at http://www.r-project.org/ (last accessed 11 January 2011). In R, the functions for computing the NGA equations can be loaded as an add-on user-contributed code, which is referred to as a ?package? in R. The details of the nga package (Kaklamanos and Thompson, 2010) are presented in this report. In addition, differences between the R and Fortran implementations will be discussed later in this report. The NGA models have established a new baseline for seismic hazard assessments, and they have been incorporated into the most recent national seismic hazard maps published by the U.S. Geological Survey (Petersen and others, 2008). However, many of the new models are significantly more complicated than previous GMPEs and, therefore, require a substantial investment of time to implement and validate. We hope that the scientific and engineering communities find our implementations to be useful in research and practice. Our implementations may be considered as an alternate to the Microsoft Excel spreadsheet implementation available on the PEER NGA project Web site (http://peer.berkeley.edu/ngawest/, last accessed 11 January 2011). The implementations in Fortran and R are more appropriate for performing batch calculations than the implementation in Microsoft Excel. Spreadsheets and Fortran code for some of the individual models also are available on the PEER NGA project Web site; our programs implement the four GMPEs simultaneously. Our programs give the same results as the programs on the PEER NGA Web site, but we offer some additional flexibility of input, additional methods of estimating unknown input parameters, and additional options for output. Although these programs have been used by the U.S. Geological Survey (USGS), Tufts University, and others, no warranty, expressed or implied, is made by Tufts or the USGS as to the accuracy or functioning of the programs and related material, nor shall the fact of distribution constitute any such warranty, and no responsibility is assumed by Tufts or the USGS in connection therewith.

Open-File Report↗

USGS2009

Along-Arc and Back-Arc Attenuation, Site Response, and Source Spectrum for the Intermediate-Depth 8 January 2006 M 6.7 Kythera, Greece, Earthquake

An M 6.7 intermediate-depth (66 km), in-slab earthquake occurring near the island of Kythera in Greece on 8 January 2006 was well recorded on networks of stations equipped with acceleration sensors and with broadband velocity sensors. All data were recorded digitally using recording instruments with resolutions ranging from almost 11 to 24 bits. We use data from these networks to study the distance dependence of the horizontal-component Fourier acceleration spectra (FAS) and horizontal-component pseudoabsolute response spectral acceleration (PSA). For purposes of simulating motions in the future, we parameterize the distance decay using several forms of the geometrical-spreading function, for each of which we derive Q as a function of frequency. By extrapolating the distance decay back to 1 km, we obtain a reference spectrum that can be used in future simulations. This spectrum requires a more complicated spectral shape than the classic single-corner-frequency model; in particular, there appears to be an enhancement of motion around 0.2-0.3 Hz that may be due to the radiation of a 3-5 sec pulse from the source. We infer a ?? 0 value of about 0.055 sec for rock stations and a stress parameter in the range of 400-600 bars. We also find distinctive differences in the site response of stations on soft soil and soil; both the FAS and the 5% damped PSA amplifications have similar peak amplitudes (about 2 and 4 for soil and soft-soil sites, respectively, relative to the rock sites) at similar frequencies (between about 0.4 and 2.0 Hz, with the soft-soil amplifications peaking at somewhat lower frequencies than the soil amplifications). One of the most distinctive features of the data is the clear difference in the motions for along-arc and back-arc stations, with the former being significantly higher than the latter over a broad range of frequencies at distances beyond about 250 km. The motions from the Kythera earthquake are roughly comparable to those from intermediate-depth earthquakes elsewhere, but they appear to be significantly higher than those from recordings of shallow earthquakes in Greece of comparable magnitude and hypocentral distance.

Bulletin of the Seismological Society of America↗

USGS2009

On baseline corrections and uncertainty in response spectrafor baseline variations commonly encountered in digital accelerograph records

Most digital accelerograph recordings are plagued by long-period drifts, best seen in the velocity and displacement time series obtained from integration of the acceleration time series. These drifts often result in velocity values that are nonzero near the end of the record. This is clearly unphysical and can lead to inaccurate estimates of peak ground displacement and long-period spectral response. The source of the long-period noise seems to be variations in the acceleration baseline in many cases. These variations could be due to true ground motion (tilting and rotation, as well as local permanent ground deformation), instrumental effects, or analog-to-digital conversion. Very often the trends in velocity are well approximated by a linear trend after the strong shaking subsides. The linearity of the trend in velocity implies that no variations in the baseline could have occurred after the onset of linearity in the velocity time series. This observation, combined with the lack of any trends in the pre-event motion, allows us to compute the time interval in which any baseline variations could occur. We then use several models of the variations in a Monte Carlo procedure to derive a suite of baseline-corrected accelerations for each noise model using records from the 1999 Chi-Chi earthquake and several earthquakes in Turkey. Comparisons of the mean values of the peak ground displacements, spectral displacements, and residual displacements computed from these corrected accelerations for the different noise models can be used as a guide to the accuracy of the baseline corrections. For many of the records considered here the mean values are similar for each noise model, giving confidence in the estimation of the mean values. The dispersion of the ground-motion measures increases with period and is noise-model dependent. The dispersion of inelastic spectra is greater than the elastic spectra at short periods but approaches that of the elastic spectra at longer periods. The elastic spectra from the most basic processing, in which only the pre-event mean is removed from the acceleration time series, do not diverge from the baseline-corrected spectra until periods of 10-20 sec or more for the records studied here, implying that for many engineering purposes elastic spectra can be used from records with no baseline correction or filtering.

Bulletin of the Seismological Society of America↗

USGS2009

Implementation of the SSHAC Guidelines for Level 3 and 4 PSHAs - Experience gained from actual applications

In April 1997, after four years of deliberations, the Senior Seismic Hazard Analysis Committee released its report 'Recommendations for Probabilistic Seismic Hazard Analysis: Guidance on Uncertainty and Use of Experts' through the U.S. Nuclear Regulatory Commission as NUREG/CR-6372, hereafter SSHAC (1997). Known informally ever since as the 'SSHAC Guidelines', SSHAC (1997) addresses why and how multiple expert opinions - and the intrinsic uncertainties that attend them - should be used in Probabilistic Seismic Hazard Analyses (PSHA) for critical facilities such as commercial nuclear power plants. Ten years later, in September 2007, the U.S. Geological Survey (USGS) entered into a 13-month agreement with the U.S. Nuclear Regulatory Commission (NRC) titled 'Practical Procedures for Implementation of the SSHAC Guidelines and for Updating PSHAs'. The NRC was interested in understanding and documenting lessons learned from recent PSHAs conducted at the higher SSHAC Levels (3 and 4) and in gaining input from the seismic community for updating PSHAs as new information became available. This study increased in importance in anticipation of new applications for nuclear power facilities at both existing and new sites. The intent of this project was not to replace the SSHAC Guidelines but to supplement them with the experience gained from putting the SSHAC Guidelines to work in practical applications. During the course of this project, we also learned that updating PSHAs for existing nuclear power facilities involves very different issues from the implementation of the SSHAC Guidelines for new facilities. As such, we report our findings and recommendations from this study in two separate documents, this being the first. The SSHAC Guidelines were written without regard to whether the PSHAs to which they would be applied were site-specific or regional in scope. Most of the experience gained to date from high-level SSHAC studies has been for site-specific cases, although three ongoing (as of this writing) studies are regional in scope. Updating existing PSHAs will depend more critically on the differences between site-specific and regional studies, and we will also address these differences in more detail in the companion report. Most of what we report here and in the second report on updating PSHAs emanates from three workshops held by the USGS at their Menlo Park facility: 'Lessons Learned from SSHAC Level 3 and 4 PSHAs' on January 30-31, 2008; 'Updates to Existing PSHAs' on May 6-7, 2008; and 'Draft Recommendations, SSHAC Implementation Guidance' on June 4-5, 2009. These workshops were attended by approximately 40 scientists and engineers familiar with hazard studies for nuclear facilities. This company included four of the authors of SSHAC (1997) and four other experts whose contributions to this document are mentioned in the Acknowledgments section; numerous scientists and engineers who in one role or another have participated in one or more high-level SSHAC PSHAs summarized later in this report; and representatives of the nuclear industry, the consulting world, the regulatory community, and academia with a keen interest and expertise in hazard analysis. This report is a community-based set of recommendations to NRC for improved practical procedures for implementation of the SSHAC Guidelines. In an early publication specifically addressing the SSHAC Guidelines, Hanks (1997) noted that the SSHAC Guidelines were likely to evolve for some time to come, and this remains true today. While the broad philosophical and theoretical dimensions of the SSHAC Guidelines will not change, much has been learned during the past decade from various applications of the SSHAC Guidelines to real PSHAs in terms of how they are implemented. We anticipate that, in their practical applications, the SSHAC Guidelines will continue to evolve as more experience is gained from future SSHAC applications. Indeed, to the extent that every PSHA has its

Open-File Report↗