Geology Reports⌕ Search

Geology topics

Robert B. Herrmann

Publications and source records attributed to Robert B. Herrmann.

17 recordsLinked to original sources

Estimating earthquake source depth using teleseismic broadband waveform modeling at the USGS National Earthquake Information Center

The U.S. Geologic Survey National Earthquake Information Center (NEIC) monitors global seismicity, producing a catalog of earthquake source parameters in near-real-time to provide information that can help mitigate the societal impact of earthquakes. The NEIC commonly relies on teleseismic observations to constrain earthquake source parameters (e.g., location, depth, magnitude, and mechanism) due to a lack of local and regional observations. For these ‘teleseismic’ events, depth phase (i.e., pP , sP ) arrival time observations provide the best estimate on source depth. However, depth phases are often difficult to accurately identify and/or pick. Therefore, NEIC relies on waveform modeling, such as those determined from W-phase ( M ww ), body wave ( M wb ), and regional ( M wr ) moment tensor estimations, to provide constraints on source depth. While depth estimates from these approaches are informative, higher frequency observations provide more precise estimates because depth phases are more prominently observed at higher frequencies. Here, we present NEIC’s relatively high-frequency (~0.04 to 1 Hz) teleseismic waveform modeling approach, termed Synthetic Depth Phase Modeling (SynDepth), for determining source depth. SynDepth was developed to provide NEIC with a tool that enables rapid, accurate, and quantifiable estimates of earthquake source depth in cases where locator depths are not reliable. This relatively simple and fast procedure searches over 1 km-incremented source depths and an expanding triangular source-time function to find the best-fitting solution. We compare automatic SynDepth solutions for a dataset of 1,216 earthquakes (M5.5-M7.6) between 2017 and 2021 to NEIC-derived depth estimates from other methods. Our approach provides a robust depth estimate for earthquakes lacking local arrival time data, and it minimizes the need for analyst review of depth-phase picks ( pP , sP ) or using predefined ‘fixed’ depths.

Seismological Research Letters↗

RMT focal plane sensitivity to seismic network geometry and faulting style

Modern tectonic studies often use regional moment tensors (RMTs) to interpret the seismotectonic framework of an earthquake or earthquake sequence; however, despite extensive use, little existing work addresses RMT parameter uncertainty. Here, we quantify how network geometry and faulting style affect RMT sensitivity. We examine how data-model fits change with fault plane geometry (strike and dip) for varying station configurations. We calculate the relative data fit for incrementally varying geometries about a best-fitting solution, applying our workflow to real and synthetic seismograms for both real and hypothetical station distributions and earthquakes. Initially, we conduct purely observational tests, computing RMTs from synthetic seismograms for hypothetical earthquakes and a series of well-behaved network geometries. We then incorporate real data and station distributions from the International Maule Aftershock Deployment (IMAD), which recorded aftershocks of the 2010 M W 8.8 Maule earthquake, and a set of regional stations capturing the ongoing earthquake sequence in Oklahoma and southern Kansas. We consider RMTs computed under three scenarios: (1) real seismic records selected for high data quality; (2) synthetic seismic records with noise computed for the observed source-station pairings and (3) synthetic seismic records with noise computed for all possible station-source pairings. To assess RMT sensitivity for each test, we observe the ‘fit falloff’, which portrays how relative fit changes when strike or dip varies incrementally; we then derive the ranges of acceptable strikes and dips by identifying the span of solutions with relative fits larger than 90 per cent of the best fit. For the azimuthally incomplete IMAD network, Scenario 3 best constrains fault geometry, with average ranges of 45° and 31° for strike and dip, respectively. In Oklahoma, Scenario 3 best constrains fault dip with an average range of 46°; however, strike is best constrained by Scenario 1, with a range of 26°. We draw two main conclusions from this study. (1) Station distribution impacts our ability to constrain RMTs using waveform time-series; however, in some tectonic settings, faulting style also plays a significant role and (2) increasing station density and data quantity (both the number of stations and the number of individual channels) does not necessarily improve RMT constraint. These results may be useful when organizing future seismic deployments (e.g. by concentrating stations in alignment with anticipated nodal planes), and in computing RMTs, either by guiding a more rigorous data selection process for input data or informing variable weighting among the selected data (e.g. by eliminating the transverse component when strike-slip mechanisms are expected).

Geophysical Journal International↗

Reactivated faulting near Cushing, Oklahoma: Increased potential for a triggered earthquake in an area of United States strategic infrastructure

In October 2014 two moderate-sized earthquakes ( Mw 4.0 and 4.3) struck south of Cushing, Oklahoma, below the largest crude oil storage facility in the world. Combined analysis of the spatial distribution of earthquakes and regional moment tensor focal mechanisms indicate reactivation of a subsurface unnamed and unmapped left-lateral strike-slip fault. Coulomb failure stress change calculations using the relocated seismicity and slip distribution determined from regional moment tensors, allow for the possibility that the Wilzetta-Whitetail fault zone south of Cushing, Oklahoma, could produce a large, damaging earthquake comparable to the 2011 Prague event. Resultant very strong shaking levels (MMI VII) in the epicentral region present the possibility of this potential earthquake causing moderate to heavy damage to national strategic infrastructure and local communities.

Oklahoma↗

Earthquake hypocenters and focal mechanisms in central Oklahoma reveal a complex system of reactivated subsurface strike-slip faulting

The sharp increase in seismicity over a broad region of central Oklahoma has raised concern regarding the source of the activity and its potential hazard to local communities and energy industry infrastructure. Since early 2010, numerous organizations have deployed temporary portable seismic stations in central Oklahoma in order to record the evolving seismicity. In this study, we apply a multiple-event relocation method to produce a catalog of 3,639 central Oklahoma earthquakes from late 2009 through 2014. RMT source parameters were determined for 195 of the largest and best-recorded earthquakes. Combining RMT results with relocated seismicity enabled us to determine the length, depth and style-of-faulting occurring on reactivated subsurface fault systems. Results show that the majority of earthquakes occur on near vertical, optimally oriented (NE-SW and NW-SE), strike-slip faults in the shallow crystalline basement. These are necessary first order observations required to assess the potential hazards of individual faults in Oklahoma.

Oklahoma↗

Efforts to monitor and characterize the recent increasing seismicity in central Oklahoma

The sharp increase in seismicity over a broad region of central Oklahoma has raised concerns regarding the source of the activity and its potential hazard to local communities and energy-industry infrastructure. Efforts to monitor and characterize the earthquake sequences in central Oklahoma are reviewed. Since early 2010, numerous organizations have deployed temporary portable seismic stations in central Oklahoma to record the evolving seismicity. A multiple-event relocation method is applied to produce a catalog of central Oklahoma earthquakes from late 2009 into early 2015. Regional moment tensor (RMT) source parameters were determined for the largest and best-recorded earthquakes. Combining RMT results with relocated seismicity enabled determination of the length, depth, and style of faulting occurring on reactivated subsurface fault systems. It was found that the majority of earthquakes occur on near-vertical, optimally oriented (northeast-southwest and northwest-southeast) strike-slip faults in the shallow crystalline basement. In 2014, 17 earthquakes occurred with magnitudes of 4 or larger. It is suggested that these recently reactivated fault systems pose the greatest potential hazard to the region.

Oklahoma↗

Using regional moment tensors to constrain the kinematics and stress evolution of the 2010–2013 Canterbury earthquake sequence, South Island, New Zealand

On September 3, 2010, a M W 7.0 (U.S. Geological Survey moment magnitude) earthquake ruptured across the Canterbury Plains in South Island, New Zealand. Since then, New Zealand GNS Science has recorded over 10,000 aftershocks M L 2.0 and larger, including three destructive ~ M W 6.0 earthquakes near Christchurch. We treat the Canterbury earthquake sequence as an intraplate earthquake sequence, and compare its kinematics to an Andersonian model for fault slip in a uniform stress field. We determined moment magnitudes and double couple solutions for 150 earthquakes having M W 3.7 and larger through the use of a waveform inversion technique using data from broadband seismic stations on South Island, New Zealand. The majority (126) of these double couple solutions have strike-slip focal mechanisms, with right-lateral slip on ENE fault planes or equivalently left-lateral slip on SSE fault planes. The remaining focal mechanisms indicate reverse faulting, except for two normal faulting events. The strike-slip segments have compatible orientations for slip in a stress field with a horizontal σ 1 oriented ~ N115°E, and horizontal σ 3 . The preference for right lateral strike-slip earthquakes suggests that these structures are inherited from previous stages of deformation. Reverse slip is interpreted to have occurred on previously existing structures in regions with an absence of existing structures optimally oriented for strike-slip deformation. Despite the variations in slip direction and faulting style, most aftershocks had nearly the same P-axis orientation, consistent with the regional σ 1 . There is no evidence for significant changes in these stress orientations throughout the Canterbury earthquake sequence.

South Island↗

The Mw 5.8 Mineral, Virginia, earthquake of August 2011 and aftershock sequence: constraints on earthquake source parameters and fault geometry

The M w 5.8 earthquake of 23 August 2011 (17:51:04 UTC) (moment, M0 5.7×10 17 N·m) occurred near Mineral, Virginia, within the central Virginia seismic zone and was felt by more people than any other earthquake in United States history. The U.S. Geological Survey (USGS) received 148,638 felt reports from 31 states and 4 Canadian provinces. The USGS PAGER system estimates as many as 120,000 people were exposed to shaking intensity levels of IV and greater, with approximately 10,000 exposed to shaking as high as intensity VIII. Both regional and teleseismic moment tensor solutions characterize the earthquake as a northeast‐striking reverse fault that nucleated at a depth of approximately 7±2 km. The distribution of reported macroseismic intensities is roughly ten times the area of a similarly sized earthquake in the western United States (Horton and Williams, 2012). Near‐source and far‐field damage reports, which extend as far away as Washington, D.C., (135 km away) and Baltimore, Maryland, (200 km away) are consistent with an earthquake of this size and depth in the eastern United States (EUS). Within the first few days following the earthquake, several government and academic institutions installed 36 portable seismograph stations in the epicentral region, making this among the best‐recorded aftershock sequences in the EUS. Based on modeling of these data, we provide a detailed description of the source parameters of the mainshock and analysis of the subsequent aftershock sequence for defining the fault geometry, area of rupture, and observations of the aftershock sequence magnitude–frequency and temporal distribution. The observed slope of the magnitude–frequency curve or b‐value for the aftershock sequence is consistent with previous EUS studies (b=0.75), suggesting that most of the accumulated strain was released by the mainshock. The aftershocks define a rupture that extends between approximately 2–8 km in depth and 8–10 km along the strike of the fault plane. Best‐fit modeling of the geometry of the aftershock sequence defines a rupture plane that strikes N36°E and dips to the east‐southeast at 49.5°. Moment tensor solutions of the mainshock and larger aftershocks are consistent with the distribution of aftershock locations, both indicating reverse slip along a northeast–southwest striking southeast‐dipping fault plane.

Virginia↗

Monitoring the Earthquake source process in North America

With the implementation of the USGS National Earthquake Information Center Prompt Assessment of Global Earthquakes for Response system (PAGER), rapid determination of earthquake moment magnitude is essential, especially for earthquakes that are felt within the contiguous United States. We report an implementation of moment tensor processing for application to broad, seismically active areas of North America. This effort focuses on the selection of regional crustal velocity models, codification of data quality tests, and the development of procedures for rapid computation of the seismic moment tensor. We systematically apply these techniques to earthquakes with reported magnitude greater than 3.5 in continental North America that are not associated with a tectonic plate boundary. Using the 0.02-0.10 Hz passband, we can usually determine, with few exceptions, moment tensor solutions for earthquakes with M w as small as 3.7. The threshold is significantly influenced by the density of stations, the location of the earthquake relative to the seismic stations and, of course, the signal-to-noise ratio. With the existing permanent broadband stations in North America operated for rapid earthquake response, the seismic moment tensor of most earthquakes that are M w 4 or larger can be routinely computed. As expected the nonuniform spatial pattern of these solutions reflects the seismicity pattern. However, the orientation of the direction of maximum compressive stress and the predominant style of faulting is spatially coherent across large regions of the continent.

Bulletin of the Seismological Society of America↗

Improved earthquake monitoring in the central and eastern United States in support of seismic assessments for critical facilities

Evaluation of seismic monitoring capabilities in the central and eastern United States for critical facilities - including nuclear powerplants - focused on specific improvements to understand better the seismic hazards in the region. The report is not an assessment of seismic safety at nuclear plants. To accomplish the evaluation and to provide suggestions for improvements using funding from the American Recovery and Reinvestment Act of 2009, the U.S. Geological Survey examined addition of new strong-motion seismic stations in areas of seismic activity and addition of new seismic stations near nuclear power-plant locations, along with integration of data from the Transportable Array of some 400 mobile seismic stations. Some 38 and 68 stations, respectively, were suggested for addition in active seismic zones and near-power-plant locations. Expansion of databases for strong-motion and other earthquake source-characterization data also was evaluated. Recognizing pragmatic limitations of station deployment, augmentation of existing deployments provides improvements in source characterization by quantification of near-source attenuation in regions where larger earthquakes are expected. That augmentation also supports systematic data collection from existing networks. The report further utilizes the application of modeling procedures and processing algorithms, with the additional stations and the improved seismic databases, to leverage the capabilities of existing and expanded seismic arrays.

Open-File Report↗

The April 18, 2008 Illinois earthquake: an ANSS monitoring success

The largest-magnitude earthquake in the past 20 years struck near Mt. Carmel in southeastern Illinois on Friday morning, 18 April 2008 at 09:36:59 UTC (04:37 CDT). The Mw 5.2 earthquake was felt over an area that spanned Chicago and Atlanta, with about 40,000 reports submitted to the U.S. Geological Survey (USGS) “Did You Feel It?” system. There were at least six felt aftershocks greater than magnitude 3 and 20 aftershocks with magnitudes greater than 2 located by regional and national seismic networks. Portable instrumentation was deployed by researchers of the University of Memphis and Indiana University (the first portable station was installed at about 23:00 UTC on 18 April). The portable seismographs were deployed both to capture near-source, high-frequency ground motions for significant aftershocks and to better understand structure along the active fault. The previous similar-size earthquake within the Wabash Valley seismic zone (WVSZ) of southeastern Illinois and southwestern Indiana was a magnitude 5.0 in June 1987. The seismicity associated with the WVSZ is thought to occur in a complex horst and graben system of Precambrian igneous and metamorphic units at depths between 12 and 20 km. Paleoliquefaction evidence suggests several major shaking events have occurred within the past 12,000 years (Munson et al. 1997).

Illinois↗

Shallow P- and S-wave velocities and site resonances in the St. Louis region, Missouri-Illinois

As part of the seismic hazard-mapping efforts in the St. Louis metropolitan area we determined the compressional and shear-wave velocities (Vp and Vs) to about a 40-m depth at 17 locations in this area. The Vs measurements were made using high-resolution seismic refraction and reflection methods. We find a clear difference in the Vs profiles between sites located on the river floodplains and those located in the upland urban areas of St. Louis. Vs30 (average Vs to 30-m depth) values in floodplain areas range from 200 to 290 m/s (NEHRP category D) and contrast with sites on the upland areas of St. Louis, which have Vs30 values ranging from 410 to 785 m/s (NEHRP categories C and B). The lower Vs30 values and earthquake recordings in the floodplains suggest a greater potential for stronger and more prolonged ground shaking in an earthquake. Spectral analysis of a M3.6 earthquake recorded on the St. Louis-area ANSS seismograph network indicates stronger shaking and potentially damaging S-wave resonant frequencies at NEHRP category D sites compared to ground motions at a rock site located on the Saint Louis University campus. ?? 2007, Earthquake Engineering Research Institute.

Earthquake Spectra↗

Faulting parameters of the January 16, 1994 Wyomissing Hills, Pennsylvania earthquakes

Two events dominated the January 1994, Wyomissing, PA earthquake sequence, an M w 4.0 foreshock, followed by an M w 4.6 mainshock. We modeled regional waveforms to estimate the event depth and the moment tensors for the two largest events in the sequence, and examine teleseismic waveforms recorded on the ARCESS short-period seismic array to estimate the depth and source time function of the mainshock. Our data constrain the depth of the events to be shallower than 5 km, and prefer a depth of 3–5 km. For an assumed depth of 3 km, the mainshock moment tensor is 75% double couple, with (the major double couple) planes striking at 135°N, 347°N, dips of 49°, 46°, and rakes of 68°, 114°. The estimated moment is 8.9 × 10 22 dyne-cm. The P axis strikes 241°N and plunges 2°, the Tension axis strikes 336°N and plunges 73°. The foreshock inversion results are virtually identical to the mainshock results; for a source depth of three km, we find a major double couple with a strike, dip, and rake of 121°N, 60°, and 66°, respectively. The seismic moment for the foreshock is 1.2 × 10 22 dyne-cm, which is approximately 13% of the mainshock moment release. These events did not excite high-frequency Lg waves as effectively as typical eastern North American events, and the mainshock had a stress drop in the range of 25–50 bars.

Pennsylvania↗

Broadband seismology and small regional seismic networks

In the winter of 1811-12, three of the largest historic earthquakes in the United States occurred near New Madrid, Missouri. Seismicity continues to the present day throughout a tightly clustered pattern of epicenters centered on the bootheel of Missouri, including parts of northeastern Arkansas, northwestern Tennessee, western Kentucky, and southern Illinois. In 1990, the New Madrid seismic zone/Central United States became the first seismically active region east of the Rocky Mountains to be designated a priority research area within the National Earthquake Hazards Reduction Program (NEHRP). This Professional Paper is a collection of papers, some published separately, presenting results of the newly intensified research program in this area. Major components of this research program include tectonic framework studies, seismicity and deformation monitoring and modeling, improved seismic hazard and risk assessments, and cooperative hazard mitigation studies.

Professional Paper↗

Earthquake hazard studies in southeastern Missouri

This report is intended to be a compact review of the current "state of the art" of operating a seismic telemetry network as applied to the New Madrid area, and incidentally, to the Adak, Alaska network because of the author's position as technical director of that network under a joint agreement with CIRES of the University of Colorado. The experiences and opinions expressed herein are necessarily those of the author and as such are limited to his acquaintance with the matter, and should not be construed to contradict or diminish the efforts of others operating telemetry networks under considerably different circumstances.

Missouri↗