Geology Reports⌕ Search

Geology topics

J. Boatwright

Publications and source records attributed to J. Boatwright.

26 records · Page 2Linked to original sources
USGS2001

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

Mapping the rupture process of moderate earthquakes by inverting accelerograms

We present a waveform inversion method that uses recordings of small events as Green's functions to map the rupture growth of moderate earthquakes. The method fits P and S waveforms from many stations simultaneously in an iterative procedure to estimate the subevent rupture time and amplitude relative to the Green's function event. We invert the accelerograms written by two moderate Parkfield earthquakes using smaller events as Green's functions. The first earthquake ( M = 4.6) occurred on November 14, 1993, at a depth of 11 km under Middle Mountain, in the assumed preparation zone for the next Parkfield main shock. The second earthquake ( M = 4.7) occurred on December 20, 1994, some 6 km to the southeast, at a depth of 9 km on a section of the San Andreas fault with no previous microseismicity and little inferred coseismic slip in the 1966 Parkfield earthquake. The inversion results are strikingly different for the two events. The average stress release in the 1993 event was 50 bars, distributed over a geometrically complex area of 0.9 km 2 . The average stress release in the 1994 event was only 6 bars, distributed over a roughly elliptical area of 20 km 2 . The ruptures of both events appear to grow spasmodically into relatively complex shapes: the inversion only constrains the ruptures to grow more slowly than the S wave velocity but does not use smoothness constraints.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Using safety inspection data to estimate shaking intensity for the 1994 Northridge earthquake

We map the shaking intensity suffered in Los Angeles County during the 17 January 1994, Northridge earthquake using municipal safety inspection data. The intensity is estimated from the number of buildings given red, yellow, or green tags, aggregated by census tract. Census tracts contain from 200 to 4000 residential buildings and have an average area of 6 km 2 but are as small as 2 and 1 km 2 in the most densely populated areas of the San Fernando Valley and downtown Los Angeles, respectively. In comparison, the zip code areas on which standard MMI intensity estimates are based are six times larger, on average, than the census tracts. We group the buildings by age (before and after 1940 and 1976), by number of housing units (one, two to four, and five or more), and by construction type, and we normalize the tags by the total number of similar buildings in each census tract. We analyze the seven most abundant building categories. The fragilities (the fraction of buildings in each category tagged within each intensity level) for these seven building categories are adjusted so that the intensity estimates agree. We calibrate the shaking intensity to correspond with the modified Mercalli intensities ( MMI ) estimated and compiled by Dewey et al. (1995); the shapes of the resulting isoseismals are similar, although we underestimate the extent of the MMI = 6 and 7 areas. The fragility varies significantly between different building categories (by factors of 10 to 20) and building ages (by factors of 2 to 6). The post-1940 wood-frame multi-family (≧5 units) dwellings make up the most fragile building category, and the post-1940 woodframe single-family dwellings make up the most resistant building category.

California↗
USGS1996

Frictional constraints on crustal faulting

We consider how variations in fault frictional properties affect the phenomenology of earthquake faulting. In particular, we propose that lateral variations in fault friction produce the marked heterogeneity of slip observed in large earthquakes. We model these variations using a rate- and state-dependent friction law, where we differentiate velocity-weakening behavior into two fields: the strong seismic field is very velocity weakening and the weak seismic field is slightly velocity weakening. Similarly, we differentiate velocity-strengthening behavior into two fields: the compliant field is slightly velocity strengthening and the viscous field is very velocity strengthening. The strong seismic field comprises the seismic slip concentrations, or asperities. The two “intermediate” fields, weak seismic and compliant, have frictional velocity dependences that are close to velocity neutral: these fields modulate both the tectonic loading and the dynamic rupture process. During the interseismic period, the weak seismic and compliant regions slip aseismically, while the strong seismic regions remain locked, evolving into stress concentrations that fail only in main shocks. The weak seismic areas exhibit most of the interseismic activity and aftershocks but can also creep seismically. This “mixed” frictional behavior can be obtained from a sufficiently heterogeneous distribution of the critical slip distance. The model also provides a mechanism for rupture arrest: dynamic rupture fronts decelerate as they penetrate into unloaded complaint or weak seismic areas, producing broad areas of accelerated afterslip. Aftershocks occur on both the weak seismic and compliant areas around a fault, but most of the stress is diffused through aseismic slip. Rapid afterslip on these peripheral areas can also produce aftershocks within the main shock rupture area by reloading weak fault areas that slipped in the main shock and then healed. We test this frictional model by comparing the seismicity and the coseismic slip for the 1966 Parkfield, 1979 Coyote Lake, and 1984 Morgan Hill earthquakes. The interevent seismicity and aftershocks appear to occur on fault areas outside the regions of significant slip: these regions are interpreted as either weak seismic or compliant, depending on whether or not they manifest interevent seismicity.

Journal of Geophysical Research B: Solid Earth↗
USGS1989

Inverting measurements of surface slip on the Superstition Hills fault

We derive and test a set of inversions of surface-slip measurements based on the empirical relation u(t) = uf/(1 + T/t)c proposed by Sharp and Saxton (1989) to estimate the final slip uf, the power-law exponent c, and the power-law duration T. At short times, Sharp's relation behaves like the simple power law, u(t) ∼ u1tc, where u1 is the initial slip, that is, the slip at 1 day after the earthquake. At long times, the slip approaches the final slip asymptotically. The inversions are designed in part to exploit the accuracy of measurements of differential slip; that is, measurements of surface slip which are made relative to a set of nails or stakes emplaced after the earthquake. We apply the inversions to slip measurements made at 53 sites along the Superstition Hills fault for the 11 months following the M = 6.2 and 6.6 earthquakes of 24 November 1987. In general, estimates of the initial slip and the power law exponent are well resolved, while estimates of the power-law durations and the final slip are less well resolved because the durations of the surface slip measurements are often less than the derived power-law durations. The slip on the three fault strands is a relatively smooth function of position; the initial slip and final slip are well correlated. The time dependence of surface slip at the 53 sites is roughly similar along the entire fault, where the power-law exponents are distributed as c = 0.14 ± 0.04 and the power-law durations range from 100 < T < 1000 days.

California↗
USGS1988

The seismic radiation from composite models of faulting

The failure of an asperity, i.e., the dynamic rupture of a small fault area with finite stress drop surrounded by a broken or weak fault area which has no stress drop but which slips after the asperity fails, is proposed as a model for the rupture process of a subevent in a composite earthquake. The rupture area of the composite earthquake surrounding the subevent is modeled by the weak fault area surrounding the asperity in the subevent model. The resulting seismic moment of the subevent is proportional to the stress drop and the rupture area of the subevent, as well as the radius of the composite earthaquake. By setting the stress drops of the asperity models equal to the dynamic stress drops of the subevents, the composite earthquake can be modeled as the sum of a set of subevents which cover the rupture area of the composite earthquake. The scaling of the high- and low-frequency radiation from composite earthquakes composed of asperities is commensurate with generally observed spectral scaling laws, in contrast to composite earthquakes composed of cracks, or smaller earthquakes. A simple filtering strategy is proposed for filtering the waveforms radiated by cracks to approximate the waveforms radiated by asperities. The P and S waves radiated by an ML = 5.2 earthquake which occurred on 9 May 1983, at Coalinga, California, are simulated using the P and S waves radiated by an ML = 3.6 aftershock. The aftershock waveforms are first filtered to approximate the radiation from asperities with the appropriate rupture areas, and then the waveforms of 12 asperity subevents are summed together to simulate the waveforms and spectra of the composite earthquake.

Bulletin of the Seismological Society of America↗
USGS1988

Teleseismic and near-field analysis of the Nahanni earthquakes in the Northwest Territories, Canada

The analysis of the Nahanni earthquakes of October 5, 1985 (MS 6.6), and December 23, 1985 (MS 6.9), will have important implications for the assessment of seismic hazards in intraplate environments. To maximize the information available to seismic engineers, broadband data recorded teleseismically are analyzed jointly with strong-motion data recorded in the near field. The time-domain analysis of teleseismic data yields the source mechanisms, depths, and complexities of rupture of each earthquake. Both earthquakes occurred as shallow thrusts with centroid depths (6 to 7 km) and shallowly dipping fault planes that correspond well with the aftershock distributions obtained from a local survey run by the Canadian Geological Survey. The shallow nodal plane for the October 5 earthquake dips 30° to the WSW, while the shallow nodal planes of the subevents for the December 23 earthquake dip an average of 23° to the WSW. The October 5 earthquake has an impulsive initial rupture, followed by a weak subevent of longer duration but smaller moment release. The December 23 earthquake exhibits more complexity, being comprised of three subevents of similar size. The subevent delays derived from the teleseismic analysis are used to help interpret arrivals in records of ground velocity recorded in the near field of the December 23 earthquake. The rupture geometries inferred from the joint near- and far-field analysis suggest that the rupture processes were unusually complicated and that the 2g peak that occurs late in one of the near-field records could be a localized phenomenon. Spectral analyses of the teleseismic P waves yield the following source parameters for the October 5 and December 23 earthquakes, respectively: the seismic moments are 1.2 and 1.8 × 1026 dyne-cm, the radiated energies are 1.8 and 2.8 × 1021 dyne-cm, and the dynamic stress drops are 65 and 50 bar. The acceleration source spectra of both earthquakes exhibit an intermediate slope (| üα(ω) | ∝ω) from 0.03 to 0.3 Hz, suggesting that the earthquakes represent the failure of asperities. Extrapolating the teleseismic P-wave spectra to estimate the near-field S-wave spectra yields good fits to the acceleration spectra from two strong motion records, but underestimates the spectra from a third strong motion record with the strongest, but possibly localized, accelerations.

Bulletin of the Seismological Society of America↗