USGS Science⌕ Search

Geology topics

G. L. Choy

Publications and source records attributed to G. L. Choy.

15 recordsLinked to original sources
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↗
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↗
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↗
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↗
USGS1999

Seismic subduction of the Nazca Ridge as shown by the 1996-97 Peru earthquakes

By rupturing more than half of the shallow subduction interface of the Nazca Ridge, the great November 12, 1996 Peruvian earthquake contradicts the hypothesis that oceanic ridges subduct aseismically. The mainshock's rupture has a length of about 200 km and has an average slip of about 1.4 m. Its moment is 1.5 x 1028 dyne-cm and the corresponding M(w) is 8.0. The mainshock registered three major episodes of moment release as shown by a finite fault inversion of teleseismically recorded broadband body waves. About 55% of the mainshock's total moment release occurred south of the Nazca Ridge, and the remaining moment release occurred at the southern half of the subduction interface of the Nazca Ridge. The rupture south of the Nazca Ridge was elongated parallel to the ridge axis and extended from a shallow depth to about 65 km depth. Because the axis of the Nazca Ridge is at a high angle to the plate convergence direction, the subducting Nazca Ridge has a large southwards component of motion, 5 cm/yr parallel to the coast. The 900-1200 m relief of the southwards sweeping Nazca Ridge is interpreted to act as a 'rigid indenter,' causing the greatest coupling south of the ridge's leading edge and leading to the large observed slip. The mainshock and aftershock hypocenters were relocated using a new procedure that simultaneously inverts local and teleseismic data. Most aftershocks were within the outline of the Nazca Ridge. A three-month delayed aftershock cluster' occurred at the northern part of the subducting Nazca Ridge. Aftershocks were notably lacking at the zone of greatest moment release, to the south of the Nazca Ridge. However, a lone foreshock at the southern end of this zone, some 140 km downstrike of the mainshock's epicenter, implies that conditions existed for rupture into that zone. The 1996 earthquake ruptured much of the inferred source zone of the M(w) 7.9-8.2 earthquake of 1942, although the latter was a slightly larger earthquake. The rupture zone of the 1996 earthquake is immediately north of the seismic gap left by the great earthquakes (M(w) ~8.8-9.1) of 1868 and 1877. The M(w) 8.0 Antofagasta earthquake of 1995 occurred at the southern end of this great seismic gap. The M(w) 8.2 deep-focus Bolivian earthquake of 1994 occurred directly downdip of the 1868 portion of that gap. The recent occurrence of three significant earthquakes on the periphery of the great seismic gap of the 1868 and 1877 events, among other factors, may signal an increased seismic potential for that zone.

Pure and Applied Geophysics↗
USGS1997

The rupture process of the Manjil, Iran earthquake of 20 June 1990 and implications for intraplate strike-slip earthquakes

In terms of seismically radiated energy or moment release, the earthquake of 20 January 1990 in the Manjil Basin-Alborz Mountain region of Iran is the second largest strike-slip earthquake to have occurred in an intracontinental setting in the past decade. It caused enormous loss of life and the virtual destruction of several cities. Despite a very large meizoseismal area, the identification of the causative faults has been hampered by the lack of reliable earthquake locations and conflicting field reports of surface displacement. Using broadband data from global networks of digitally recording seismographs, we analyse broadband seismic waveforms to derive characteristics of the rupture process. Complexities in waveforms generated by the earthquake indicate that the main shock consisted of a tiny precursory subevent followed in the next 20 seconds by a series of four major subevents with depths ranging from 10 to 15 km. The focal mechanisms of the major subevents, which are predominantly strike-slip, have a common nodal plane striking about 285°–295°. Based on the coincidence of this strike with the dominant tectonic fabric of the region we presume that the EW striking planes are the fault planes. The first major subevent nucleated slightly south of the initial precursor. The second subevent occurred northwest of the initial precursor. The last two subevents moved progressively southeastward of the first subevent in a direction collinear with the predominant strike of the fault planes. The offsets in the relative locations and the temporal delays of the rupture subevents indicate heterogeneous distribution of fracture strength and the involvement of multiple faults. The spatial distribution of teleseismic aftershocks, which at first appears uncorrelated with meizoseismal contours, can be decomposed into stages. The initial activity, being within and on the periphery of the rupture zone, correlates in shape and length with meizoseismal lines. In the second stage of activity the aftershock zone expands and appears to cluster about the geomorphic and geologic features several tens of kilometres from the rupture zone. The activity is interpreted as a regional response to quasistatic stress migration along zones of tectonic weakness. The radiated energy of the main shock and the estimate of seismic moment yields an apparent stress of 20 bars. High apparent stress may be typical of strike slip earthquakes occurring in intracontinental environments undergoing continental collision.

Studia Geophysica et Geodaetica↗
USGS1996

Inelastic models of lithospheric stress - I. Theory and application to outer-rise plate deformation

Outer-rise stress distributions determined in the manner that mechanical engineers evaluate inelastic stress distributions within conventional materials are contrasted with those predicted using simple elastic-plate models that are frequently encountered in studies of outer-rise seismicity. This comparison indicates that the latter are inherently inappropriate for studies of intraplate earthquakes, which are a direct manifestation of lithospheric inelasticity. We demonstrate that the common practice of truncating elastically superimposed stress profiles so that they are not permitted to exceed laboratory-based estimates of lithospheric yield strength will result in an accurate characterization of lithospheric stress only under relatively restrictive circumstances. In contrast to elastic-plate models, which predict that lithospheric stress distributions depend exclusively upon the current load, inelastic plate models predict that stress distributions are also significantly influenced by the plate-loading history, and, in many cases, this influence is the dominant factor in determining the style of potential seismicity (e.g. thrust versus normal faulting). Numerous ‘intuitive’ interpretations of outer-rise earthquakes have been founded upon the implicit assumption that a unique relationship exists between a specified combination of plate curvature and in-plane force, and the resulting lithospheric stress distribution. We demonstrate that the profound influence of deformation history often invalidates such interpretations. Finally, we examine the reliability of ‘yield envelope’ representations of lithospheric strength that are constructed on the basis of empirically determined frictional sliding relationships and silicate plastic-flow laws. Although representations of this nature underestimate the strength of some major interplate faults, such as the San Andreas, they appear to represent a reliable characterization of the strength of intraplate oceanic lithosphere.

Geophysical Journal International↗
USGS1996

Inelastic models of lithospheric stress - II. Implications for outer-rise seismicity and dynamics

Outer-rise seismicity and dynamics are examined using inelastic models of lithospheric deformation, which allow a more realistic characterization of stress distributions and failure behaviour. We conclude that thrust- and normal-faulting outer-rise earthquakes represent substantially different states of stress within the oceanic lithosphere. Specifically, the normal-faulting events occur in response to downward plate bending, which establishes the ‘standard’, bending-dominated state of outer-rise stress, and the thrust-faulting events occur in response to an elevated level of in-plane compression, which develops only in response to exceptional circumstances. This interpretation accounts for the observation that normal-faulting outer-rise earthquakes occur more frequently and are more widely distributed than their thrust-faulting counterparts, an observation for which the simple bending model offers no explanation. In addition, attributing both thrust- and normal-faulting outer-rise earthquakes to plate bending implies that both classes of events should occur within relatively close lateral proximity to one another because both are allegedly a manifestation of the same bendingdominated stress distribution, whereas, in reality, this is not observed. We propose that the tendency for thrust-faulting outer-rise earthquakes to exhibit greater source depths than their normal-faulting counterparts (an observation that is frequently cited in support of the bending interpretation of the former) is merely a consequence of the fact that bending-induced tension is confined to the upper lithosphere. Our model predicts that outer-rise in-plane-force variations may promote thrust-faulting outerrise activity prior to an underthrusting interplate subduction earthquake and normalfaulting outer-rise activity following such an earthquake, but that both forms of outerrise activity are unlikely to be associated with the same subduction earthquake. A corollary implication of our model is t...

Geophysical Journal International↗
USGS1989

Measuring the size of an earthquake

Earthquakes range broadly in size. A rock-burst in an Idaho silver mine may involve the fracture of 1 meter of rock; the 1965 Rat Island earthquake in the Aleutian arc involved a 650-kilometer length of the Earth's crust. Earthquakes can be even smaller and even larger. If an earthquake is felt or causes perceptible surface damage, then its intensity of shaking can be subjectively estimated. But many large earthquakes occur in oceanic areas or at great focal depths and are either simply not felt or their felt pattern does not really indicate their true size. Today, state-of-the-art seismic systems transmit data from the seismograph via telephone line and satellite directly to a central digital computer. A preliminary location, depth-of-focus, and magntidue can now be obtained within minutes of the onset of an earthquake. The only limiting factor is how long the seismic wave stake to travel from the epicenter to the stations-usually less than 10 minutes.

Earthquakes & Volcanoes (USGS)↗
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↗
USGS1987

Analysis of broadband seismograms from selected IASPEI events

Broadband seismograms of body waves that are flat to displacement and velocity in the frequency range from 0.01 to 5.0 Hz can now be routinely obtained for most earthquakes of magnitude greater than about 5.5. These records are obtained either directly or through multichannel deconvolution of waveforms from digitally recording seismograph stations. In contrast to data from conventional narrowband seismographs, broadband records have sufficient frequency content to define the source-time functions of body waves, even for shallow events for which the source functions of direct and surface-reflected phases may overlap. Broadband seismograms for selected IASPEI events are systematically analysed to identify depth phases and the presence of subevents. The procedure results in improved estimates of focal depth, identification of subevents in complex earthquakes, and better resolution of focal mechanisms. We propose that it is now possible for reporting agencies, such as the National Earthquake Information Center, to use broadband digital waveforms routinely in the processing of earthquake data. ?? 1987.

Physics of the Earth and Planetary Interiors↗