Geology Reports⌕ Search

Geology topics

C.H. Thurber

Publications and source records attributed to C.H. Thurber.

11 recordsLinked to original sources
USGS2020

Relocated aftershocks and background seismicity in eastern Indonesia shed light on the 2018 Lombok and Palu earthquake sequences

High seismicity rates in eastern Indonesia occur due to the complex interaction of several tectonic plates which resulted in two deadly, destructive earthquake sequences that occurred in Lombok Island and the city of Palu, Sulawesi in 2018. The first sequence began in July with an M w 6.4 event near Lombok, culminating in an M w 7.0 event 8 d later. This was then followed by a nearby M w 6.9 event 12 d later. Approximately 1000 km to the northeast, a separate sequence began several weeks later near Palu where an M w 7.5 event occurred that triggered a tsunami. In this study, we present hypocentre relocations for both earthquake sequences as well as all other regional earthquakes in eastern Indonesia. The relocations were performed using a teleseismic double-difference relocation method and arrival times for P and S waves from stations at local, regional, and teleseismic distances. The catalogue and phase data were taken from the Agency for Meteorology, Climatology and Geophysics (BMKG) of Indonesia and the International Seismological Centre (ISC) for the period of April 2009 through November 2018. The relocated catalogue provides an improved view of seismicity in eastern Indonesia over the study period, sharpening locations and interpretations of seismogenic features throughout the region. In the Lombok area, the relocated earthquakes clearly show a backarc thrust to the north of the Sunda-Banda Arc transition zone. The relocated aftershocks show that the destructive M w 7.0 and M w 6.9 earthquakes of the Lombok sequence ruptured two different regions: The M w 7.0 earthquake propagated westward, whereas the M w 6.9 earthquake propagated eastward. The entire sequence of Lombok earthquakes was most likely started by the M w 6.4 event as the initial event or foreshock, which then triggered backarc thrusts on both sides. Several weeks later and far to the northeast, the M w 7.5 Palu earthquake occurred along the Palu-Koro Fault, filling a seismic gap that had not ruptured in an M w 6.0 event or larger since at least 1900. The distribution of aftershocks indicates that the northern part of the Palu-Koro Fault has lower relative seismicity rates than the southern part at shallow depths, and that off fault aftershocks are mostly located to the east of the Palu-Koro Fault.

Lombok, Palu↗
USGS2012

Determination and uncertainty of moment tensors for microearthquakes at Okmok Volcano, Alaska

Efforts to determine general moment tensors (MTs) for microearthquakes in volcanic areas are often hampered by small seismic networks, which can lead to poorly constrained hypocentres and inadequate modelling of seismic velocity heterogeneity. In addition, noisy seismic signals can make it difficult to identify phase arrivals correctly for small magnitude events. However, small volcanic earthquakes can have source mechanisms that deviate from brittle double-couple shear failure due to magmatic and/or hydrothermal processes. Thus, determining reliable MTs in such conditions is a challenging but potentially rewarding pursuit. We pursued such a goal at Okmok Volcano, Alaska, which erupted recently in 1997 and in 2008. The Alaska Volcano Observatory operates a seismic network of 12 stations at Okmok and routinely catalogues recorded seismicity. Using these data, we have determined general MTs for seven microearthquakes recorded between 2004 and 2007 by inverting peak amplitude measurements of P and S phases. We computed Green's functions using precisely relocated hypocentres and a 3-D velocity model. We thoroughly assessed the quality of the solutions by computing formal uncertainty estimates, conducting a variety of synthetic and sensitivity tests, and by comparing the MTs to solutions obtained using alternative methods. The results show that MTs are sensitive to station distribution and errors in the data, velocity model and hypocentral parameters. Although each of the seven MTs contains a significant non-shear component, we judge several of the solutions to be unreliable. However, several reliable MTs are obtained for a group of previously identified repeating events, and are interpreted as compensated linear-vector dipole events.

Alaska↗
USGS2011

The Augustine magmatic system as revealed by seismic tomography and relocated earthquake hypocenters from 1994 through 2009

We incorporate 14 years of earthquake data from the Alaska Volcano Observatory with data from a 1975 controlled‐source seismic experiment to obtain the three‐dimensional P and S wave velocity structure and the first high‐precision earthquake locations at Augustine Volcano to be calculated in a fully three‐dimensional velocity model. Velocity tomography shows two main features beneath Augustine: a narrow, high‐velocity column beneath the summit, extending from ∼2 km depth to the surface, and elevated velocities on the south flank. Our relocation results allow a thorough analysis of the spatio‐temoral patterns of seismicity and the relationship to the magmatic and eruptive activity. Background seismicity is centered beneath the summit at an average depth of 0.6 km above sea level. In the weeks leading to the January 2006 eruption of Augustine, seismicity focused on a NW‐SE line along the trend of an inflating dike. A series of drumbeat earthquakes occurred in the early weeks of the eruption, indicating further magma transport through the same dike system. During the six months following the onset of the eruption, the otherwise quiescent region 1 to 5 km below sea level centered beneath the summit became seismically active with two groups of earthquakes, differentiated by frequency content. The deep longer‐period earthquakes occurred during the eruption and are interpreted as resulting from the movement of magma toward the summit, and the post‐eruptive shorter‐period earthquakes may be due to the relaxation of an emptied magma tube. The seismicity subsequently returned to its normal background rates and patterns.

Alaska↗
USGS2010

A California statewide three-dimensional seismic velocity model from both absolute and differential times

We obtain a seismic velocity model of the California crust and uppermost mantle using a regional-scale double-difference tomography algorithm. We begin by using absolute arrival-time picks to solve for a coarse three-dimensional (3D) P velocity (V P ) model with a uniform 30 km horizontal node spacing, which we then use as the starting model for a finer-scale inversion using double-difference tomography applied to absolute and differential pick times. For computational reasons, we split the state into 5 subregions with a grid spacing of 10 to 20 km and assemble our final statewide V P model by stitching together these local models. We also solve for a statewide S-wave model using S picks from both the Southern California Seismic Network and USArray, assuming a starting model based on the V P results and a V P =V S ratio of 1.732. Our new model has improved areal coverage compared with previous models, extending 570 km in the SW-NE directionand 1320 km in the NW-SE direction. It also extends to greater depth due to the inclusion of substantial data at large epicentral distances. Our V P model generally agrees with previous separate regional models for northern and southern California, but we also observe some new features, such as high-velocity anomalies at shallow depths in the Klamath Mountains and Mount Shasta area, somewhat slow velocities in the northern Coast Ranges, and slow anomalies beneath the Sierra Nevada at midcrustal and greater depths. This model can be applied to a variety of regional-scale studies in California, such as developing a unified statewide earthquake location catalog and performing regional waveform modeling.

Bulletin of the Seismological Society of America↗
USGS2007

Three-dimensional P wave velocity model for the San Francisco Bay region, California

[1] A new three-dimensional P wave velocity model for the greater San Francisco Bay region has been derived using the double-difference seismic tomography method, using data from about 5,500 chemical explosions or air gun blasts and approximately 6,000 earthquakes. The model region covers 140 km NE-SW by 240 km NW-SE, extending from 20 km south of Monterey to Santa Rosa and reaching from the Pacific coast to the edge of the Great Valley. Our model provides the first regional view of a number of basement highs that are imaged in the uppermost few kilometers of the model, and images a number of velocity anomaly lows associated with known Mesozoic and Cenozoic basins in the study area. High velocity ( V p > 6.5 km/s) features at ∼15-km depth beneath part of the edge of the Great Valley and along the San Francisco peninsula are interpreted as ophiolite bodies. The relocated earthquakes provide a clear picture of the geometry of the major faults in the region, illuminating fault dips that are generally consistent with previous studies. Ninety-five percent of the earthquakes have depths between 2.3 and 15.2 km, and the corresponding seismic velocities at the hypocenters range from 4.8 km/s (presumably corresponding to Franciscan basement or Mesozoic sedimentary rocks of the Great Valley Sequence) to 6.8 km/s. The top of the seismogenic zone is thus largely controlled by basement depth, but the base of the seismogenic zone is not restricted to seismic velocities of ≤6.3 km/s in this region, as had been previously proposed.

California↗
USGS2004

Precise relocation of earthquakes following the 15 June 1991 eruption of Mount Pinatubo (Philippines)

The 15 June 1991 climactic eruption of Mount Pinatubo (Philippines) was followed by intense seismicity that remained at a high level for several months. We located 10,839 events recorded between 1 July and mid-December 1991. In contrast to the preeruptive seismicity which was focused in two groups below the summit area, posteruptive events were widely distributed below and around the volcano. The classification of the events indicates the presence of several large multiplets, and the application of relative relocation techniques to the similar events by calculating high-precision delays between traces outlines a number of clear seismogenic structures. We used different methods to confirm the validity of our results; these tests indicate that reliable features can be detected with a small monitoring network. While the main cluster of activity can be attributed to an intrusive process starting from below the 15 June crater, the volcanic origin of the seismic activity in the other areas is more difficult to establish. Away from the summit, relocations define streaks or planes which are oriented predominantly southwest-northeast, with in several cases the presence of northwest-southeast conjugate structures. Most of the composite focal mechanisms that we could determine indicate predominantly strike-slip, right-lateral faulting. Our results indicate that most of the seismicity that occurred after the 15 June eruption is related to the east-west regional compressional stress field related to the subduction. We suggest that the regional stress field induces seismicity along new or preexisting faults in the medium surrounding the volcano where the stress field was locally disturbed by the volcanic eruption.

Mount Pinatubo↗
USGS2004

Dome growth behavior at Soufriere Hills Volcano, Montserrat, revealed by relocation of volcanic event swarms, 1995-1996

We have relocated a subset of events from the digital waveform catalogue of ???17,000 volcanic microearthquakes recorded between July 1995 and February 1996 at Soufriere Hills Volcano (SHV), Montserrat, using a cross-correlation-based phase repicking technique with a joint location method. Hypocenters were estimated for 3914 earthquakes having five or more corrected P-wave picks. The seismic source region collapsed to a volume of ???1 km3 from an initial ???100 km3. Relocated events represent 36 swarms, each containing nearly identical waveforms, having source dimensions of 10 to 100 m in diameter and spatial separations on the order of 500 m or less. Each swarm occurred over a span of several hours to a few days.Triggered data appear to miss between 65% and 98% of the events that occur within these swarms, based on review of helicorder records. Visual estimates of summit dome growth show a rough correspondence between episodes of intense swarming and increases in extruded magma, although dome observations are too sparse to make a direct comparison for this time period. The limited depth range over which dome-growth-related events occur is consistent with a dynamic model of cyclic plug extrusion behavior in the shallow conduit, governed by magma supply rate, overpressure buildup and physical properties of the magma and conduit geometry. Seismic sources may occur in locally overpressured regions that result from microlite formation in a zone of rapid decompression; we propose that this zone exists in the vicinity of a detachment plane associated with the cyclic plug extrusion. ?? 2004 Elsevier B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗
USGS2004

Upper crustal structure from the Santa Monica Mountains to the Sierra Nevada, Southern California: Tomographic results from the Los Angeles Regional Seismic Experiment, Phase II (LARSE II)

In 1999, the U.S. Geological Survey and the Southern California Earthquake Center (SCEC) collected refraction and low-fold reflection data along a 150-km-long corridor extending from the Santa Monica Mountains northward to the Sierra Nevada. This profile was part of the second phase of the Los Angeles Region Seismic Experiment (LARSE II). Chief imaging targets included sedimentary basins beneath the San Fernando and Santa Clarita Valleys and the deep structure of major faults along the transect, including causative faults for the 1971 M 6.7 San Fernando and 1994 M 6.7 Northridge earthquakes, the San Gabriel Fault, and the San Andreas Fault. Tomographic modeling of first arrivals using the methods of Hole (1992) and Lutter et al. (1999) produces velocity models that are similar to each other and are well resolved to depths of 5-7.5 km. These models, together with oil-test well data and independent forward modeling of LARSE II refraction data, suggest that regions of relatively low velocity and high velocity gradient in the San Fernando Valley and the northern Santa Clarita Valley (north of the San Gabriel Fault) correspond to Cenozoic sedimentary basin fill and reach maximum depths along the profile of ???4.3 km and >3 km , respectively. The Antelope Valley, within the western Mojave Desert, is also underlain by low-velocity, high-gradient sedimentary fill to an interpreted maximum depth of ???2.4 km. Below depths of ???2 km, velocities of basement rocks in the Santa Monica Mountains and the central Transverse Ranges vary between 5.5 and 6.0 km/sec, but in the Mojave Desert, basement rocks vary in velocity between 5.25 and 6.25 km/sec. The San Andreas Fault separates differing velocity structures of the central Transverse Ranges and Mojave Desert. A weak low-velocity zone is centered approximately on the north-dipping aftershock zone of the 1971 San Fernando earthquake and possibly along the deep projection of the San Gabriel Fault. Modeling of gravity data, using densities inferred from the velocity model, indicates that different velocity-density relationships hold for both sedimentary and basement rocks as one crosses the San Andreas Fault. The LARSE II velocity model can now be used to improve the SCEC Community Velocity Model, which is used to calculate seismic amplitudes for large scenario earthquakes.

Bulletin of the Seismological Society of America↗
USGS1997

Fault-zone guided waves from explosions in the San Andreas fault at Parkfield and Cienega Valley, California

Fault-zone guided waves were successfully excited by near-surface explosions in the San Andreas fault zone both at Parkfield and Cienega Valley, central California. The guided waves were observed on linear, three-component seismic arrays deployed across the fault trace. These waves were not excited by explosions located outside the fault zone. The amplitude spectra of guided waves show a maximum peak at 2 Hz at Parkfield and 3 Hz at Cienega Valley. The guided wave amplitude decays sharply with observation distance from the fault trace. The explosion-excited fault-zone guided waves are similar to those generated by earthquakes at Parkfield but have lower frequencies and travel more slowly. These observations suggest that the fault-zone wave guide has lower seismic velocities as it approaches the surface at Parkfield. We have modeled the waveforms as S waves trapped in a low-velocity wave guide sandwiched between high-velocity wall rocks, resulting in Love-type fault-zone guided waves. While the results are nonunique, the Parkfield data are adequately fit by a shallow wave guide 170 m wide with an S velocity 0.85 km/sec and an apparent Q ∼ 30 to 40. At Cienega Valley, the fault-zone wave guide appears to be about 120 m wide with an S velocity 0.7 km/sec and a Q ∼ 30.

California↗