USGS Science⌕ Search

Geology topics

David M. Boore

Publications and source records attributed to David M. Boore.

At least 37 records · Page 2Linked to original sources
USGS2008

TSPP - A collection of FORTRAN programs for processing and manipulating time series

This report lists a number of FORTRAN programs that I have developed over the years for processing and manipulating strong-motion accelerograms. The collection is titled TSPP, which stands for Time Series Processing Programs. I have excluded 'strong-motion accelerograms' from the title, however, as the boundary between 'strong' and 'weak' motion has become blurred with the advent of broadband sensors and high-dynamic range dataloggers, and many of the programs can be used with any evenly spaced time series, not just acceleration time series. This version of the report is relatively brief, consisting primarily of an annotated list of the programs, with two examples of processing, and a few comments on usage. I do not include a parameter-by-parameter guide to the programs. Future versions might include more examples of processing, illustrating the various parameter choices in the programs. Although these programs have been used by the U.S. Geological Survey, no warranty, expressed or implied, is made by the USGS as to the accuracy or functioning of the programs and related program material, nor shall the fact of distribution constitute any such warranty, and no responsibility is assumed by the USGS in connection therewith. The programs are distributed on an 'as is' basis, with no warranty of support from me. These programs were written for my use and are being publicly distributed in the hope that others might find them as useful as I have. I would, however, appreciate being informed about bugs, and I always welcome suggestions for improvements to the codes. Please note that I have made little effort to optimize the coding of the programs or to include a user-friendly interface (many of the programs in this collection have been included in the software usdp (Utility Software for Data Processing), being developed by Akkar et al. (personal communication, 2008); usdp includes a graphical user interface). Speed of execution has been sacrificed in favor of a code that is intended to be easy to understand, although on modern computers speed of execution is rarely a problem. I will be pleased if users incorporate portions of my programs into their own applications; I only ask that reference be made to this report as the source of the programs.

Open-File Report↗
USGS2008

Comments on potential geologic and seismic hazards affecting proposed liquefied natural gas site in Santa Monica Bay, California

In a letter to the U.S. Geological Survey (USGS) dated March 25, 2008, Representative Jane Harman (California 36th district) requested advice on geologic hazards that should be considered in the review of a proposed liquefied natural gas (LNG) facility off the California coast in Santa Monica Bay. In 2004, the USGS responded to a similar request from Representative Lois Capps, regarding two proposed LNG facilities offshore Ventura County, Calif., with a report summarizing potential geologic and seismic hazards (Ross and others, 2004). The proposed LNG Deepwater Port (DWP) facility includes single point moorings (SPMs) and 35 miles of underwater pipelines. The DWP submersible buoys, manifolds, and risers would be situated on the floor of the southern Santa Monica Basin, in 3,000 feet of water, about 23 miles offshore of the Palos Verdes Peninsula. Twin 24-inch diameter pipelines would extend northeastward from the buoys across the basin floor, up the basin slope and across the continental shelf, skirting north around the Santa Monica submarine canyon. Figure 1 provides locations of the project and geologic features. Acronyms are defined in table 1. This facility is being proposed in a region of known geologic hazards that arise from both the potential for strong earthquakes and geologic processes related to sediment transport and accumulation in the offshore environment. The probability of a damaging earthquake (considered here as magnitude 6.5 or greater) in the next 30 years within about 30 miles (50 km) of the proposed pipeline ranges from 16% at the pipeline's offshore end to 48% where it nears land (Petersen, 2008). Earthquakes of this magnitude are capable of producing strong shaking, surface fault offsets, liquefaction phenomena, landslides, underwater turbidity currents and debris flow avalanches, and tsunamis. As part of the DWP license application for the Woodside Natural Gas proposal in Santa Monica Bay (known as the OceanWay Secure Energy Project), Fugro West, Inc., had already prepared a document discussing geologic hazards in the area, titled 'Exhibit B Topic Report 6 - Geological Resources' (Fugro West, Inc., 2007); hereafter, this will be called the 'Geological Resources document'. The USGS agreed to evaluate the information in the Geological Resources document regarding (1) proximity of active faults to the proposed project, (2) potential magnitude of seismic events from nearby faults, (3) thoroughness of the assessment of earthquake hazards in general, (4) potential hazards from ground rupture and strong shaking, (5) potential hazards from tsunamis, and (6) other geologic hazards including landslides and debris flows. Because two new earthquake probability reports were scheduled to be released in mid-April, 2008, by the USGS and the California Geological Survey (CGS), the USGS suggested a 6-month review period to enable a thorough incorporation of this new information. Twenty-seven scientists from the USGS and the CGS reviewed various sections of the Geological Resources document. This report outlines our major conclusions. The appendix is a longer list of comments by these reviewers, grouped by section of the Geological Resources document. Before discussing our reviews, we first provide a brief overview of geologic hazards in the proposed site area. This report is a snapshot in time and any future work in the area will need to take into account ongoing research efforts. For example, USGS scientists collected seismic reflection data in the spring of 2008 to study the structure and seismic potential of several faults in the area. Their interpretations (Conrad and others, 2008a and 2008b) are too preliminary to be included in this report, but their final results, along with other researchers' studies in the project area, should be considered in any future work on the Deepwater Port project.

California↗
USGS2006

On low-frequency errors of uniformly modulated filtered white-noise models for ground motions

Low-frequency errors of a commonly used non-stationary stochastic model (uniformly modulated filtered white-noise model) for earthquake ground motions are investigated. It is shown both analytically and by numerical simulation that uniformly modulated filter white-noise-type models systematically overestimate the spectral response for periods longer than the effective duration of the earthquake, because of the built-in low-frequency errors in the model. The errors, which are significant for low-magnitude short-duration earthquakes, can be eliminated by using the filtered shot-noise-type models (i. e. white noise, modulated by the envelope first, and then filtered).

Earthquake Engineering and Structural Dynamics↗
USGS2005

An evaluation of the strong ground motion recorded during the May 1, 2003 Bingol Turkey, earthquake

An important record of ground motion from a M 6.4 earthquake occurring on May 1, 2003, at epicentral and fault distances of about 12 and 9 km, respectively, was obtained at a station near the city of Bingöl, Turkey. The maximum peak ground values of 0.55g and 36 cm/s are among the largest ground-motion amplitudes recorded in Turkey. From simulations and comparisons with ground motions from other earthquakes of comparable magnitude, we conclude that the ground motion over a range of frequencies is unusually high. Site response may be responsible for the elevated ground motion, as suggested from analysis of numerous aftershock recordings from the same station. The mainshock motions have some interesting seismological features, including ramps between the P- and S-wave that are probably due to near- and intermediate-field elastic motions and strong polarisation oriented at about 39 degrees to the fault (and therefore not in the fault-normal direction). Simulations of motions from an extended rupture explain these features. The N10E component shows a high-amplitude spectral acceleration at a period of 0.15 seconds resulting in a site specific design spectrum that significantly overestimates the actual strength and displacement demands of the record. The pulse signal in the N10E component affects the inelastic spectral displacement and increases the inelastic displacement demand with respect to elastic demand for very long periods.

Journal of Earthquake Engineering↗
USGS2004

Estimating V̄s(30) (or NEHRP site classes) from shallow velocity models (depths < 30 m)

The average velocity to 30 m [V??s(30)] is a widely used parameter for classifying sites to predict their potential to amplify seismic shaking. In many cases, however, models of shallow shear-wave velocities, from which V??s(30) can be computed, do not extend to 30 m. If the data for these cases are to be used, some method of extrapolating the velocities must be devised. Four methods for doing this are described here and are illustrated using data from 135 boreholes in California for which the velocity model extends to at least 30 m. Methods using correlations between shallow velocity and V??s(30) result in significantly less bias for shallow models than the simplest method of assuming that the lowermost velocity extends to 30 m. In addition, for all methods the percent of sites misclassified is generally less than 10% and falls to negligible values for velocity models extending to at least 25 m. Although the methods using correlations do a better job on average of estimating V??s(30), the simplest method will generally result in a lower value of V??s(30) and thus yield a more conservative estimate of ground motion [which generally increases as V??s(30) decreases].

Bulletin of the Seismological Society of America↗
USGS2004

Comments on potential geologic and seismic hazards affecting coastal Ventura County, California

This report examines the regional seismic and geologic hazards that could affect proposed liquefied natural gas (LNG) facilities in coastal Ventura County, California. Faults throughout this area are thought to be capable of producing earthquakes of magnitude 6.5 to 7.5, which could produce surface fault offsets of as much as 15 feet. Many of these faults are sufficiently well understood to be included in the current generation of the National Seismic Hazard Maps; others may become candidates for inclusion in future revisions as research proceeds. Strong shaking is the primary hazard that causes damage from earthquakes and this area is zoned with a high level of shaking hazard. The estimated probability of a magnitude 6.5 or larger earthquake (comparable in size to the 2003 San Simeon quake) occurring in the next 30 years within 30 miles of Platform Grace is 50-60%; for Cabrillo Port, the estimate is a 35% likelihood. Combining these probabilities of earthquake occurrence with relationships that give expected ground motions yields the estimated seismic-shaking hazard. In parts of the project area, the estimated shaking hazard is as high as along the San Andreas Fault. The combination of long-period basin waves and LNG installations with large long-period resonances potentially increases this hazard.

California↗
USGS2003

A compendium of P- and S-wave velocities from surface-to-borehole logging; summary and reanalysis of previously published data and analysis of unpublished data

For over 28 years, the U.S. Geological Survey (USGS) has been acquiring seismic velocity and geologic data at a number of locations in California, many of which were chosen because strong ground motions from earthquakes were recorded at the sites. The method for all measurements involves picking first arrivals of P- and S-waves from a surface source recorded at various depths in a borehole (as opposed to noninvasive methods, such as the SASW method [e.g., Brown et al., 2002]). The results from most of the sites are contained in a series of U.S. Geological Survey Open-File Reports (see References). Until now, none of the results have been available as computer files, and before 1992 the interpretation of the arrival times was in terms of piecemeal interval velocities, with no attempt to derive a layered model that would fit the travel times in an overall sense (the one exception is Porcella, 1984). In this report I reanalyze all of the arrival times in terms of layered models for P- and for S-wave velocities at each site, and I provide the results as computer files. In addition to the measurements reported in the open-file reports, I also include some borehole results from other reports, as well as some results never before published. I include data for 277 boreholes (at the time of this writing; more will be added to the web site as they are obtained), all in California (I have data from boreholes in Washington and Utah, but these will be published separately). I am also in the process of interpreting travel time data obtained using a seismic cone penetrometer at hundreds of sites; these data can be interpreted in the same way of those obtained from surface-to-borehole logging. When available, the data will be added to the web site (see below for information on obtaining data from the World Wide Web (WWW)). In addition to the basic borehole data and results, I provide information concerning strong-motion stations that I judge to be close enough to the boreholes that the borehole velocity models can be used as the velocity models beneath the stations.

California↗
USGS2002

Borehole velocity measurements at five sites that recorded the Cape Mendocino, California earthquake of 25 April, 1992

The U.S. Geological Survey (USGS), as part of an ongoing program to acquire seismic velocity and geologic data at locations that recorded strong-ground motions during earthquakes, has investigated five sites in the Fortuna, California region (Figure 1). We selected drill sites at strong-motion stations that recorded high accelerations (Table 1) from the Cape Mendocino earthquake (M 7.0) of 25 April 1992 (Oppenheimer et al., 1993). The boreholes were drilled to a nominal depth of 95 meters (310 ft) and cased with schedule 80 pvc-casing grouted in place at each location. S-wave and P-wave data were acquired at each site using a surface source and a borehole three-component geophone. This report contains the velocity models interpreted from the borehole data and gives reference to locations and peak accelerations at the selected strong-motion stations.

California↗
USGS2001

Borehole P- and S-wave velocity at thirteen stations in Southern California

The U.S. Geological Survey (USGS), as part of a program to acquire seismic velocity data at locations of strong-ground motion in earthquakes (e.g., Gibbs et al., 2000), has investigated thirteen additional sites in the Southern California region. Of the thirteen sites, twelve are in the vicinity of Whittier, California, and one is located in San Bernardino, California. Several deployments of temporary seismographs were made after the Whittier Narrows, California earthquake of 1 October 1987 (Mueller et al., 1988). A deployment, between 2 October and 9 November 1987, was the motivation for selection of six of the drill sites. Temporary portable seismographs at Hoover School (HOO), Lincoln School (LIN), Corps of Engineers Station (NAR), Olive Junior High School (OLV), Santa Anita Golf Course (SAG), and Southwestern Academy (SWA) recorded significant aftershock data. These portable sites, with the exception of Santa Anita Golf Course, were co-sited with strong-motion recorders. Stations at HOO, Lincoln School Whittier (WLB), Saint Paul High School (STP), Alisos Adult School (EXC), Cerritos College Gymnasium (CGM), Cerritos College Physical Science Building (CPS), and Cerritos College Police Building (CPB) were part of an array of digital strong-motion stations deployed from "bedrock" in Whittier to near the deepest part of the Los Angeles basin in Norwalk. Although development and siting of this new array (partially installed at the time of this writing) was generally motivated by the Whittier Narrows earthquake, these new sites (with the exception of HOO) were not part of any Whittier Narrows aftershock deployments. A similar new digital strong-motion site was installed at the San Bernardino Fire Station during the same time frame. Velocity data were obtained to depths of about 90 meters at two sites, 30 meters at seven sites, and 18 to 25 meters at four sites. Lithology data from the analysis of cuttings and samples was obtained from the two 90-meter deep holes and from five of the shallower holes to supplement the velocity interpretation. The two 90-meter boreholes (SB1, CPB) have been instrumented with borehole seismometers for continuous monitoring of earthquake activity (Rogers et al., 1998). No drill samples or cuttings were obtained from the other six sites, but driller's logs were scanned for major changes noted there. The velocity models at those sites were interpreted using only the measured data and major changes in the driller's log as noted above. The sites are shown in Figure 1 and listed in Table 1, which gives references to information regarding the strong-motion data. Several hundred strong-motion records of the main-shock were written by this moderate size earthquake (ML = 5.9), making it important from a scientific and engineering prospective (Brady et al., 1988; Shakal et al., 1988).

California↗
USGS2000

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↗
USGS2000

SMSIM--Fortran programs for simulating ground motions from earthquakes: Version 2.0.--a revision of OFR 96-80-A

A simple and powerful method for simulating ground motions is based on the assumption that the amplitude of ground motion at a site can be specified in a deterministic way, with a random phase spectrum modified such that the motion is distributed over a duration related to the earthquake magnitude and to distance from the source. This method of simulating ground motions often goes by the name "the stochastic method." It is particularly useful for simulating the higher-frequency ground motions of most interest to engineers, and it is widely used to predict ground motions for regions of the world in which recordings of motion from damaging earthquakes are not available. This simple method has been successful in matching a variety of ground-motion measures for earthquakes with seismic moments spanning more than 12 orders of magnitude. One of the essential characteristics of the method is that it distills what is known about the various factors affecting ground motions (source, path, and site) into simple functional forms that can be used to predict ground motions. SMSIM is a set of programs for simulating ground motions based on the stochastic method. This Open-File Report is a revision of an earlier report (Boore, 1996) describing a set of programs for simulating ground motions from earthquakes. The programs are based on modifications I have made to the stochastic method first introduced by Hanks and McGuire (1981). The report contains source codes, written in Fortran, and executables that can be used on a PC. Programs are included both for time-domain and for random vibration simulations. In addition, programs are included to produce Fourier amplitude spectra for the models used in the simulations and to convert shear velocity vs. depth into frequency-dependent amplification. The revision to the previous report is needed because the input and output files have changed significantly, and a number of new programs have been included in the set.

Open-File Report↗
USGS1995

The Los Angeles Dam Story

In 1971, the near-failure of a dam during a magnitude 6.7 earthquake forced 80,000 people to evacuate their residences. In 1994, the replacement dam survived an almost identical earthquake with little damage. Underlying this progress in designing critical structures are years of research on the powerful shaking during large earthquakes.

California↗
USGS1994

The attenuation of seismic shear waves in quaternary alluvium in Santa Clara Valley, California

We used shear waves, generated by an air-powered source at the ground surface and recorded in a borehole, to estimate the shear-wave quality factor at strong-motion station Gilroy no. 2. We find similar values of Q using both the decay of the spectra with depth and the slope of the spectral ratio at two depths; we find no evidence of a frequency dependence of Q . The mean value of Q over the depth range 10 to 115 m is close to 10. The use of this value over the depth of the borehole and the observed travel time of 0.358 sec gives a cumulative attenuation factor t * of 0.036 sec for the upper 180 m of the Quaternary alluvium. This is comparable to the differential decay between Gilroy no. 2 and a rock site 1.9 km away (Gilroy no. 1), as measured from the decay of the high-frequency spectra of accelerograms from large earthquakes, plotted on a log-linear scale: t * = 0.05, 0.04, and 0.03 sec for the 1979 Coyote Lake, 1984 Morgan Hill, and 1989 Loma Prieta earthquakes, respectively. The similarity between the attenuations measured from the low-strain surface source and those from the larger amplitude earthquake sources suggests that increases of damping due to nonlinear wave propagation effects are limited.

California↗
USGS1989

The Richter scale: its development and use for determining earthquake source parameters

The M L scale, introduced by Richter in 1935, is the antecedent of every magnitude scale in use today. The scale is defined such that a magnitude-3 earthquake recorded on a Wood-Anderson torsion seismometer at a distance of 100 km would write a record with a peak excursion of 1 mm. To be useful, some means are needed to correct recordings to the standard distance of 100 km. Richter provides a table of correction values, which he terms -log A o , the latest of which is contained in his 1958 textbook. A new analysis of over 9000 readings from almost 1000 earthquakes in the southern California region was recently completed to redetermine the -log A o values. Although some systematic differences were found between this analysis and Richter's values (such that using Richter's values would lead to under and overestimates of M L at distances less than 40 km and greater than 200 km, respectively), the accuracy of his values is remarkable in view of the small number of data used in their determination. Richter's corrections for the distance attenuation of the peak amplitudes on Wood-Anderson seismographs apply only to the southern California region, of course, and should not be used in other areas without first checking to make sure that they are applicable. Often in the past this has not been done, but recently a number of papers have been published determining the corrections for other areas. If there are significant differences in the attenuation within 100 km between regions, then the definition of the magnitude at 100 km could lead to difficulty in comparing the sizes of earthquakes in various parts of the world. To alleviate this, it is proposed that the scale be defined such that a magnitude 3 corresponds to 10 mm of motion at 17 km. This is consistent both with Richter's definition of M L at 100 km and with the newly determined distance corrections in the southern California region. Aside from the obvious (and original) use as a means of cataloguing earthquakes according to size, ML has been used in predictions of ground shaking as a function of distance and magnitude; it has also been used in estimating energy and seismic moment. There is a good correlation of peak ground velocity and the peak motion on a Wood-Anderson instrument at the same location, as well as an observationally defined (and theoretically predicted) nonlinear relation between M L and seismic moment. An important byproduct of the establishment of the M L scale is the continuous operation of the network of Wood-Anderson seismographs on which the scale is based. The records from these instruments can be used to make relative comparisons of amplitudes and waveforms of recent and historic earthquakes; furthermore, because of the moderate gain, the instruments can write onscale records from great earthquakes at teleseismic distances and thus can provide important information about the energy radiated from such earthquakes at frequencies where many instruments have saturated.

California↗