Geology Reports⌕ Search

Geology topics

J. Boatwright

Publications and source records attributed to J. Boatwright.

At least 19 recordsLinked to original sources
USGS2017

Geophysical characterization of seismic station sites in the United States – The importance of a flexible, multi-method approach

Noninvasive geophysical site characterization methods were used in two recent projects to obtain shear-wave velocity (VS) profiles to a minimum depth of 30 m and the time-averaged VS of the upper 30 meters (VS30) at seismic station sites. These projects include the 2009 American Recovery and Reinvestment Act (ARRA) funded U.S. Geological Survey site characterization project for 191 sites in California and the Central-eastern United States (CEUS), and the 2012 Electric Power Research Institute (EPRI) funded project for 33 additional CEUS sites. These sites are located in rural to urban settings with topographic conditions ranging from relatively flat sedimentary basins to mountaintop ridges. About 60 percent of the ARRA sites and 80 percent of the EPRI sites are located on rock or have thin sediment cover over rock, including Quaternary volcanic rock, Tertiary sediments and sedimentary rock, and Mesozoic (or older) crystalline or sedimentary rock. The remaining sites consist of thick sequences of Quaternary sediments overlying older sediments and rock. ARRA sites were characterized using non-invasive active and passive surface-wave methods, including the horizontal-tovertical spectral ratio (HVSR) method and one or more of the following: spectral analysis of surface waves (SASW), multichannel analysis of surface waves (MASW; Rayleigh and Love waves) and, occasionally, array microtremor (linear and 2-D arrays) methods. P-wave seismic refraction data were also acquired at rock and shallow-rock sites. S-wave seismic refraction and/or Love-wave MASW methods were applied at sites where characterization proved difficult with Rayleighwave methods. Based on our experience from the ARRA project, we acquired Rayleigh- and Love-wave based MASW and P- and S-wave refraction data for the EPRI project at CEUS sites. The HVSR method was found to be useful for identifying shallow-rock sites and for evaluating the relative variability of the depth-to-rock interface beneath the seismic station and the testing array(s). The fundamental mode modeling assumption was generally valid at most of these sites; nevertheless, multi-mode or effective-mode modeling routines were occasionally required, particularly in the case of shallow high-velocity layers. Deep sediment sites were characterized using active and, when appropriate, passive surface-wave based methods. Rock and shallow sediment sites were generally more challenging to characterize than deep sediment sites. About 10 percent of rock sites could not be characterized using surface wave methods, thus these sites were characterized using body-wave refraction methods. Love wave methods were found to be more effective than Rayleigh wave methods at some rock and shallow-rock sites (e.g., sites with shallow rock and sites with a thin low-velocity, highly attenuating surface layer). Lateral velocity variability was found to be very common at rock and shallow-rock sites, often causing significant scatter in the surface-wave dispersion data. Seismic refraction models have demonstrated that it may not be unusual for VS30 to vary by 20 percent, or more, over small distances (several tens of meters) at such sites. Based on these experiences, it is important to consider the application of combinations of methods when using noninvasive geophysical approaches to characterize seismic site conditions.

Conference Paper↗
USGS2012

A teleseismic analysis of the New Brunswick earthquake of January 9, 1982

The analysis of the New Brunswick earthquake of January 9, 1982, has important implications for the evaluation of seismic hazards in eastern North America. Although moderate in size ( m b 5.7), it was well-recorded teleseismically. Source characteristics of this earthquake have been determined from analysis of data that were digitally recorded by the Global Digital Seismograph Network. From broadband displacement and velocity records of P waves, we have obtained a dynamic description of the rupture process as well as conventional static properties of the source. The depth of the hypocenter is estimated to be 9 km from depth phases. The focal mechanism determined from the broadband data corresponds to predominantly thrust faulting. From the variation in the waveforms the direction of slip is inferred to be updip on a west dipping NNE striking fault plane. The steep dip of the inferred fault plane suggests that the earthquake occurred on a preexisting fault that was at one time a normal fault. From an inversion of bodywave pulse durations, the estimated rupture length is 5.5 km. Average properties of the rupture process were examined by a moment tensor analysis of long-period P and SH body waves. The long-period moment of this earthquake was 5.3 × 10 24 dyne cm. The static and dynamic stress drops are 41 and 65 bars, respectively, similar to those of many earthquakes with similar moment in regions that are more seismically active. The joint epicenter determination algorithm was used to locate, relative to the mainshock, the three teleseismically recorded aftershocks that occurred through March 31, 1982. The relocated hypocenters of the aftershocks are significantly different from each other and from that of the mainshock; they provide additional support for the source dimensions inferred from the waveform analysis.

Journal of Geophysical Research Solid Earth↗
USGS2012

Acceleration spectra for subduction zone earthquakes

We estimate the source spectra of shallow earthquakes from digital recordings of teleseismic P wave groups, that is, P + p P + s P , by making frequency dependent corrections for the attenuation and for the interference of the free surface. The correction for the interference of the free surface assumes that the earthquake radiates energy from a range of depths. We apply this spectral analysis to a set of 12 subduction zone earthquakes which range in size from M S = 6.2 to 8.1, obtaining corrected P wave acceleration spectra on the frequency band from 0.01 to 2.0 Hz. Seismic moment estimates from surface waves and normal modes are used to extend these P wave spectra to the frequency band from 0.001 to 0.01 Hz. The acceleration spectra of moderate subduction zone earthquakes, that is, earthquakes whose seismic moments are less than 10 27 dyn cm, exhibit ω-square or Brune-type spectra, while the acceleration spectra of large subduction zone earthquakes, that is, earthquakes whose seismic moments are greater than 10 27 dyn cm, exhibit intermediate slopes where ü(ω) ∝ ω 5/4 for frequencies from 0.005 to 0.05 Hz. For this set of earthquakes, spectral shape appears to be a discontinuous function of seismic moment. Using reasonable assumptions for the phase characteristics, we transform the spectral shape observed for large earthquakes into the time domain to fit Ekström's (1987) moment rate functions for the M S = 8.1 Michoacan earthquake of September 19, 1985, and the M S = 7.6 Michoacan aftershock of September 21, 1985.

Journal of Geophysical Research Solid Earth↗
USGS2009

Differential energy radiation from two earthquakes in Japan with identical M w : The Kyushu 1996 and Tottori 2000 earthquakes

We examine two closely located earthquakes in Japan that had identical moment magnitudes M w but significantly different energy magnitudes M e . We use teleseismic data from the Global Seismograph Network and strong-motion data from the National Research Institute for Earth Science and Disaster Prevention's K-Net to analyze the 19 October 1996 Kyushu earthquake (M w 6.7, M e 6.6) and the 6 October 2000 Tottori earthquake (M w 6.7, M e 7.4). To obtain regional estimates of radiated energy E S we apply a spectral technique to regional (<200 km) waveforms that are dominated by S and Lg waves. For the thrust-fault Kyushu earthquake, we estimate an average regional attenuation Q(f) 230f 0:65 . For the strike-slip Tottori earthquake, the average regional attenuation is Q(f) 180f 0:6 . These attenuation functions are similar to those derived from studies of both California and Japan earthquakes. The regional estimate of E S for the Kyushu earthquake, 3:8 ?? 10 14 J, is significantly smaller than that for the Tottori earthquake, E S 1:3 ?? 10 15 J. These estimates correspond well with the teleseismic estimates of 3:9 ?? 10 14 J and 1:8 ?? 10 15 J, respectively. The apparent stress (T a = ??E s /M 0 with ?? equal to rigidity) for the Kyushu earthquake is 4 times smaller than the apparent stress for the Tottori earthquake. In terms of the fault maturity model, the significantly greater release of energy by the strike-slip Tottori earthquake can be related to strong deformation in an immature intraplate setting. The relatively lower energy release of the thrust-fault Kyushu earthquake can be related to rupture on mature faults at a subduction environment. The consistence between teleseismic and regional estimates of E S is particularly significant as teleseismic data for computing E S are routinely available for all large earthquakes whereas often there are no near-field data.

Bulletin of the Seismological Society of America↗
USGS2008

The distribution of modified mercalli intensity in the 18 April 1906 San Francisco earthquake

We analyze Boatwright and Bundock's (2005) modified Mercalli intensity (MMI) map for the 18 April 1906 San Francisco earthquake, reviewing their interpretation of the MMI scale and testing their correlation of 1906 cemetery damage with MMI intensity. We consider in detail four areas of the intensity map where Boatwright and Bundock (2005) added significantly to the intensity descriptions compiled by Lawson (1908). We show that the distribution of off-fault damage in Sonoma County suggests that the rupture velocity approached the P-wave velocity along Tomales Bay. In contrast, the falloff of intensity with distance from the fault appears approximately constant throughout Mendocino County. The intensity in Humboldt County appears somewhat higher than the intensity in Mendocino County, suggesting that the rupture process at the northern end of the rupture was relatively energetic and that there was directivity consistent with a subsonic rupture velocity on the section of the fault south of Shelter Cove. Finally, we show that the intensity sites added in Santa Cruz County change the intensity distribution so that it decreases gradually along the southeastern section of rupture from Corralitos to San Juan Bautista and implies that the stress release on this section of rupture was relatively low.

Conference Paper↗
USGS2007

The energy radiated by the 26 December 2004 Sumatra-Andaman earthquake estimated from 10-minute P-wave windows

The rupture process of the Mw 9.1 Sumatra-Andaman earthquake lasted for approximately 500 sec, nearly twice as long as the teleseismic time windows between the P and PP arrival times generally used to compute radiated energy. In order to measure the P waves radiated by the entire earthquake, we analyze records that extend from the P-wave to the S-wave arrival times from stations at distances ?? >60??. These 8- to 10-min windows contain the PP, PPP, and ScP arrivals, along with other multiply reflected phases. To gauge the effect of including these additional phases, we form the spectral ratio of the source spectrum estimated from extended windows (between TP and TS) to the source spectrum estimated from normal windows (between TP and TPP). The extended windows are analyzed as though they contained only the P-pP-sP wave group. We analyze four smaller earthquakes that occurred in the vicinity of the Mw 9.1 mainshock, with similar depths and focal mechanisms. These smaller events range in magnitude from an Mw 6.0 aftershock of 9 January 2005 to the Mw 8.6 Nias earthquake that occurred to the south of the Sumatra-Andaman earthquake on 28 March 2005. We average the spectral ratios for these four events to obtain a frequency-dependent operator for the extended windows. We then correct the source spectrum estimated from the extended records of the 26 December 2004 mainshock to obtain a complete or corrected source spectrum for the entire rupture process (???600 sec) of the great Sumatra-Andaman earthquake. Our estimate of the total seismic energy radiated by this earthquake is 1.4 ?? 1017 J. When we compare the corrected source spectrum for the entire earthquake to the source spectrum from the first ???250 sec of the rupture process (obtained from normal teleseismic windows), we find that the mainshock radiated much more seismic energy in the first half of the rupture process than in the second half, especially over the period range from 3 sec to 40 sec.

Bulletin of the Seismological Society of America↗
USGS2007

The persistence of directivity in small earthquakes

We derive a simple inversion of peak ground acceleration (PGA) or peak ground velocity (PGV) for rupture direction and rupture velocity and then test this inversion on the peak motions obtained from seven 3.5≤ M ≤4.1 earthquakes that occurred in two clusters in November 2002 and February 2003 near San Ramon, California. These clusters were located on two orthogonal strike-slip faults so that the events share the same approximate focal mechanism but not the same fault plane. Three earthquakes exhibit strong directivity, but the other four earthquakes exhibit relatively weak directivity. We use the residual PGAs and PGVs from the other six events to determine station corrections for each earthquake. The inferred rupture directions unambiguously identify the fault plane for the three earthquakes with strong directivity and for three of the four earthquakes with weak directivity. The events with strong directivity have fast rupture velocities (0.63 β ≤ v ≤0.87 β ); the events with weak directivity either rupture more slowly (0.17 β ≤ v ≤0.35 β ) or bilaterally. The simple unilateral inversion cannot distinguish between slow and bilateral ruptures: adding a bilateral rupture component degrades the fit of the rupture directions to the fault planes. By comparing PGAs from the events with strong and weak directivity, we show how an up-dip rupture in small events can distort the attenuation of peak ground motion with distance. When we compare the rupture directions of the earthquakes to the location of aftershocks in the two clusters, we find than almost all the aftershocks of the three earthquakes with strong directivity occur within 70° of the direction of rupture.

Bulletin of the Seismological Society of America↗
USGS2006

Using Modified Mercalli Intensities to estimate acceleration response spectra for the 1906 San Francisco earthquake

We derive and test relations between the Modified Mercalli Intensity (MMI) and the pseudo-acceleration response spectra at 1.0 and 0.3 s - SA(1.0 s) and SA(0.3 s) - in order to map response spectral ordinates for the 1906 San Francisco earthquake. Recent analyses of intensity have shown that MMI ??? 6 correlates both with peak ground velocity and with response spectra for periods from 0.5 to 3.0 s. We use these recent results to derive a linear relation between MMI and log SA(1.0 s), and we refine this relation by comparing the SA(1.0 s) estimated from Boatwright and Bundock's (2005) MMI map for the 1906 earthquake to the SA(1.0 s) calculated from recordings of the 1989 Loma Prieta earthquake. South of San Jose, the intensity distributions for the 1906 and 1989 earthquakes are remarkably similar, despite the difference in magnitude and rupture extent between the two events. We use recent strong motion regressions to derive a relation between SA(1.0 s) and SA(0.3 s) for a M7.8 strike-slip earthquake that depends on soil type, acceleration level, and source distance. We test this relation by comparing SA(0.3 s) estimated for the 1906 earthquake to SA(0.3 s) calculated from recordings of both the 1989 Loma Prieta and 1994 Northridge earthquakes, as functions of distance from the fault. ?? 2006, Earthquake Engineering Research Institute.

Earthquake Spectra↗
USGS2006

A stochastic estimate of ground motion at Oceano, California, for the M 6.5 22 December 2003 San Simeon earthquake, derived from aftershock recordings

The U.S. Geological Survey deployed a digital seismic station in Oceano, California, in February 2004, to investigate the cause of damage and liquefaction from the 22 December 2003 M 6.5 San Simeon earthquake. This station recorded 11 M > 2.8 aftershocks in almost 8 weeks. We analyze these recordings, together with recordings of the mainshock and the same aftershocks obtained from nearby stations in Park Hill and San Luis Obispo, to estimate the mainshock ground motion in Oceano. We estimate the Fourier amplitude spectrum using generalized spectral ratio analysis. We test a set of aftershocks as Green's functions by comparing simulated and recorded acceleration amplitude spectra for the mainshock at San Luis Obispo and Park Hill. We convolve the aftershock accelerograms with a stochastic operator to simulate the duration and phase of the mainshock accelerograms. This approximation allows us to extend the range of aftershocks that can be used as Green's functions to events nearly three magnitude units smaller than the mainshock. Our realizations for the mainshock accelerogram at Oceano yield peak ground accelerations distributed as 28% ?? 4%g. We interpret these realizations as upper bounds for the actual ground motion, because our analysis assumes a linear response, whereas the presence of liquefaction indicates that the ground behaved nonlinearly in Oceano.

Bulletin of the Seismological Society of America↗
USGS2004

Radiated energy and the rupture process of the Denali fault earthquake sequence of 2002 from broadband teleseismic body waves

Displacement, velocity, and velocity-squared records of P and SH body waves recorded at teleseismic distances are analyzed to determine the rupture characteristics of the Denali fault, Alaska, earthquake of 3 November 2002 ( M W 7.9, M e 8.1 ). Three episodes of rupture can be identified from broadband ( ∼0.1–5.0 Hz ) waveforms. The Denali fault earthquake started as a M W 7.3 thrust event. Subsequent right-lateral strike-slip rupture events with centroid depths of 9 km occurred about 22 and 49 sec later. The teleseismic P waves are dominated by energy at intermediate frequencies (0.1-1 Hz) radiated by the thrust event, while the SH waves are dominated by energy at lower frequencies (0.05-0.2 Hz) radiated by the strike-slip events. The strike-slip events exhibit strong directivity in the teleseismic SH waves. Correcting the recorded P -wave acceleration spectra for the effect of the free surface yields an estimate of 2.8 × 10 15 N m for the energy radiated by the thrust event. Correcting the recorded SH -wave acceleration spectra similarly yields an estimate of 3.3 × 10 16 N m for the energy radiated by the two strike-slip events. The average rupture velocity for the strike-slip rupture process is 1.1 β –1.2 β . The strike-slip events were located 90 and 188 km east of the epicenter. The rupture length over which significant or resolvable energy is radiated is, thus, far shorter than the 340-km fault length over which surface displacements were observed. However, the seismic moment released by these three events, 4 × 10 20 N m, was approximately half the seismic moment determined from very low-frequency analyses of the earthquake. The difference in seismic moment can be reasonably attributed to slip on fault segments that did not radiate significant or coherent seismic energy. These results suggest that very large and great strike-slip earthquakes can generate stress pulses that rapidly produce substantial slip with negligible stress drop and little discernible radiated energy on fault segments distant from the initial point of nucleation. The existence of this energy-deficient rupture mode has important implications for the evaluation of the seismic hazard of very large strike-slip earthquakes.

Alaska↗
USGS2003

Strong ground-motion prediction from Stochastic-dynamic source models

In the absence of sufficient data in the very near source, predictions of the intensity and variability of ground motions from future large earthquakes depend strongly on our ability to develop realistic models of the earthquake source. In this article we simulate near-fault strong ground motion using dynamic source models. We use a boundary integral method to simulate dynamic rupture of earthquakes by specifying dynamic source parameters (fracture energy and stress drop) as spatial random fields. We choose these quantities such that they are consistent with the statistical properties of slip heterogeneity found in finite-source models of past earthquakes. From these rupture models we compute theoretical strong-motion seismograms up to a frequency of 2 Hz for several realizations of a scenario strike-slip Mw 7.0 earthquake and compare empirical response spectra, spectra obtained from our dynamic models, and spectra determined from corresponding kinematic simulations. We find that spatial and temporal variations in slip, slip rise time, and rupture propagation consistent with dynamic rupture models exert a strong influence on near-source ground motion. Our results lead to a feasible approach to specify the variability in the rupture time distribution in kinematic models through a generalization of Andrews' (1976) result relating rupture speed to apparent fracture energy, stress drop, and crack length to 3D dynamic models. This suggests that a simplified representation of dynamic rupture may be obtained to approximate the effects of dynamic rupture without having to do full dynamic simulations.

Bulletin of the Seismological Society of America↗
USGS2003

The dependence of PGA and PGV on distance and magnitude inferred from Northern California ShakeMap data

We analyze peak ground velocity (PGV) and peak ground acceleration (PGA) data from 95 moderate (3.5 ??? M < 5.5) and 9 large (5.5 ??? M ??? 7.1) earthquakes in northern California. The 95 moderate earthquakes occurred from August 1998 through December 2002, and their peak motions were compiled and mapped by ShakeMap. The nine large earthquakes include the M 6.2 Morgan Hill earthquake in 1984, the M 6.9 Loma Prieta earthquake in 1989, and the M 7.1 Petrolia earthquake in 1992. For r > 100 km, the peak motions attenuate more rapidly than a simple power law (that is, r-??) can fit. Instead, we use an attenuation function that combines a fixed power law (r-0.7) with a fitted exponential dependence on distance, which is estimated as expt(-0.0063r) and exp(-0.0073r) for PGV and PGA, respectively, for moderate earthquakes. We regress log(PGV) and log(PGA) as functions of distance and magnitude. We assume that the scaling of log(PGV) and log(PGA) with magnitude can differ for moderate and large earthquakes, but must be continuous. Because the frequencies that carry PGV and PGA can vary with earthquake size for large earthquakes, the regression for large earthquakes incorporates a magnitude dependence in the exponential attenuation function. We fix the scaling break between moderate and large earthquakes at M 5.5; log(PGV) and log(PGA) scale as 1.06M and 1.00M, respectively, for moderate earthquakes and 0.58M and 0.31M for large earthquakes.

California↗
USGS2002

Regional estimates of radiated seismic energy

We revise the spectral technique for estimating radiated energy from recordings of large earthquakes at regional distances (Δ < 200 km) by correcting for geometric spreading and for site amplification as explicit functions of frequency. We analyze 65 recordings of the 1999 Hector Mine earthquake as functions of frequency, distance, and azimuth. For r > 27.5 km from the source, we model the geometrical spreading of the regional wavefield as r – γ where γ = 0.5 for f ≤ 0.2 Hz and γ = 0.7 for f ≥ 0.25 Hz. We fit the spectral falloff with distance using a frequency-dependent attenuation Q = 400( f /1.5) 0.6 , where Q = 400 for f ≤ 1.5 Hz. There is little directivity apparent in the corrected velocity spectra: the velocity spectra observed to the northwest along strike are amplified by a factor of 2.5 from 0.3 to 1.0 Hz and those to the southeast are amplified by a factor of 1.6 from 0.3 to 0.7 Hz. We group the stations in NEHRP site classes, using average 1-D velocity structures to estimate site amplification as a function of frequency and assuming 0.40 ≤ κ ≤ 0.55 sec for the near-surface attenuation. We increase the amplification of the soft-soil sites from 0.1 to 1.0 Hz by a factor that reaches 1.7 at 0.3 Hz because they are more strongly amplified than the NEHRP-D velocity structure predicts. We combine the 65 single-station estimates of radiated energy using an equal-azimuth weighting scheme that compensates for station distribution and incorporates the observed directivity, yielding a regional estimate of E s = 3.4 ± 0.7 × 10 22 dyne cm. This regional estimate of radiated energy corresponds closely to the teleseismic estimate of E s = 3.2 × 10 22 dyne cm.

Bulletin of the Seismological Society of America↗
USGS2002

Triggered surface slips in the Salton Trough associated with the 1999 Hector Mine, California, earthquake

Surface fracturing occurred along the southern San Andreas, Superstition Hills, and Imperial faults in association with the 16 October 1999 (Mw 7.1) Hector Mine earthquake, making this at least the eighth time in the past 31 years that a regional earthquake has triggered slip along faults in the Salton Trough. Fractures associated with the event formed discontinuous breaks over a 39-km-long stretch of the San Andreas fault, from the Mecca Hills southeastward to Salt Creek and Durmid Hill, a distance from the epicenter of 107 to 139 km. Sense of slip was right lateral; only locally was there a minor (~1 mm) vertical component of slip. Dextral slip ranged from 1 to 13 mm. Maximum slip values in 1999 and earlier triggered slips are most common in the central Mecca Hills. Field evidence indicates a transient opening as the Hector Mine seismic waves passed the southern San Andreas fault. Comparison of nearby strong-motion records indicates several periods of relative opening with passage of the Hector Mine seismic wave-a similar process may have contributed to the field evidence of a transient opening. Slip on the Superstition Hills fault extended at least 9 km, at a distance from the Hector Mine epicenter of about 188 to 196 km. This length of slip is a minimum value, because we saw fresh surface breakage extending farther northwest than our measurement sites. Sense of slip was right lateral; locally there was a minor (~1 mm) vertical component of slip. Dextral slip ranged from 1 to 18 mm, with the largest amounts found distributed (or skewed) away from the Hector Mine earthquake source. Slip triggered on the Superstition Hills fault commonly is skewed away from the earthquake source, most notably in 1968, 1979, and 1999. Surface slip on the Imperial fault and within the Imperial Valley extended about 22 km, representing a distance from the Hector Mine epicenter of about 204 to 226 km. Sense of slip dominantly was right lateral; the right-lateral component of slip ranged from 1 to 19 mm. Locally there was a minor (~1-2 mm) vertical component of slip; larger proportions of vertical slip (up to 10 mm) occurred in Mesquite basin, where scarps indicate long-term oblique-slip motion for this part of the Imperial fault. Slip triggered on the Imperial fault appears randomly distributed relative to location along the fault and source direction. Multiple surface slips, both primary and triggered slip, indicate that slip repeatedly is small at locations of structural complexity.

Bulletin of the Seismological Society of America↗
USGS2002

Comments on potential geologic and seismic hazards affecting Mare Island, Solano County, California

This report was prepared in response to a written request from the City of Vallejo, California, to the U.S. Geological Survey (USGS). By letter of October 4, 2002, the City requested that the USGS "provide advice to the City’s LNG Health and Safety Committee on its review of a potential liquid natural gas project" on the southern portion of Mare Island. The City specifically requested that the USGS advise the committee on potential hazards including fault rupture, earthquake ground motion, soil failure during earthquakes, tsunami and seiche, and landslides. The City requested that the USGS: (1) comment on these hazards, (2) describe its degree of confidence in its opinions, and (3) describe the scope of additional studies that will be needed if the City enters into an agreement with project sponsors. Advice was also requested on the selection of the safe shutdown and operating basis earthquakes as specified in the NFPA 59A standard (NFPA, 2001). This review of published reports and other publicly available information indicates that all of the hazards on which the USGS was asked to comment should be considered for the proposed project on the southern portion of Mare Island. Available information differs greatly for each of these potential hazards, and adequate understanding for design will require detailed site-specific investigations.

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