Search USGS⌕ Search

Geology topics

C.H. Thurber

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

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