Geology Reports⌕ Search

Geology topics

David A. Schmidt

Publications and source records attributed to David A. Schmidt.

6 recordsLinked to original sources

3D broadband ground motion simulations of Cascadia megathrust earthquakes: Full-margin scenarios incorporating variable magnitudes and frontal thrust splay faulting

Great megathrust earthquakes (∼ M w 8–9) that rupture along the Cascadia Subduction Zone pose significant shaking and tsunami hazards to the U.S. Pacific Northwest. Although these large earthquakes are documented in the geologic record, none have been recorded with modern instrumentation, complicating regional seismic hazard assessments. To improve constraints on seismic hazard estimates for future earthquakes, we develop 3D broadband (0–10 Hz) ground‐motion simulations for 36 full‐margin megathrust earthquake scenarios through the Cascadia Coastlines and Peoples Hazards Research Hub. These simulations expand on the M9 Project ( Frankel et al. , 2018 ; Wirth, Frankel, et al. , 2018 ) by incorporating a broader range of rupture parameters—variable magnitudes, splay faulting, and geodetically informed slip distributions. Benchmarking these simulations against the Next Generation Attenuation Subduction (NGA‐Sub) ground‐motion models shows reasonable fits for non‐basin sites, although simulations predict stronger shaking at short (0–2 s) and long (6–10 s) periods. NGA‐Sub basin terms improve fits to synthetic shaking within the Seattle Basin but underpredict amplitudes by a factor of 2 at periods of 2–5 s. We also evaluate how rupture characteristics influence regional inland shaking variability and basin response in the Puget Lowlands. Larger magnitude scenarios generally cause stronger inland shaking due to farther inland extents of slip and larger strong‐motion‐generating fault patches but show little difference in basin amplification within the Puget Lowlands compared to smaller magnitude events. Along‐dip hypocenter location strongly influences overall inland shaking, with offshore (updip) hypocenters generating on average ∼45% higher shaking at long periods (>6 s) compared to onshore (downdip) hypocenters due to rupture directivity effects. Offshore splay faults minimally affect inland shaking but alter offshore deformation important for coupled ground‐shaking and tsunami scenarios. Finally, we show that rupture directivity increases amplification within the Seattle Basin through the generation and interaction of higher order mode surface waves, an effect of 3D basin geometry that cannot be captured using 1D approximations.

British Columbia, California, Oregon, Washington↗

The impact of 3D structure on coseismic coastal land-level change and tsunami generation in the Cascadia Subduction Zone

Estimates of coseismic vertical displacements from past and potential future subduction zone earthquakes provide critical constraints on regional seismic and tsunami hazard. Many studies use elastic homogeneous half-space models to calculate vertical displacements from a specified earthquake slip distribution, neglecting complexities of 3D structure in subduction zone settings. Here, we use 3D ground motion simulations of potential Cascadia Subduction Zone earthquakes to investigate the impact of realistic 3D Earth structure on estimated vertical displacements. We show that offshore uplift increases when including 3D structure due to high fault slip in low-rigidity accretionary wedge sediments and coastal subsidence decreases due to higher rigidity material at depth. These larger offshore uplifts cause increased tsunami maximum wave heights, and coastal subsidence at paleoseismic sites decreases on average by ∼17 cm, or ∼60%, which can have important implications for the amplitude and extent of slip in earthquake scenarios developed based on paleoseismic data constraints.

California, Oregon, Washington↗

Development of a geodetic component for the U.S. West Coast Earthquake Early Warning System

An earthquake early warning (EEW) system, ShakeAlert, is under development for the West Coast of the United States. This system currently uses the first few seconds of waveforms recorded by seismic instrumentation to rapidly characterize earthquake magnitude, location, and origin time; ShakeAlert recently added a seismic line source algorithm. For large to great earthquakes, magnitudes estimated from the earliest seismic data alone generally saturate. Real‐time Global Navigation Satellite System (GNSS) data can directly measure large displacements, enabling accurate magnitude estimates for M w 7 + "> M w 7 + events, possibly before rupture termination. GNSS‐measured displacements also track evolving slip and, alone or in combination with seismic data, constrain finite‐fault models. Particularly for large‐magnitude, long‐rupture events, GNSS‐based magnitude and rupture extent estimates can improve updates to predicted shaking and thus alert accuracy. GNSS data processing centers at ShakeAlert partner institutions provide real‐time streams to the EEW system, and three geodetic EEW algorithms have been developed through the ShakeAlert collaboration. These algorithms will undergo initial testing within ShakeAlert’s computational architecture using a suite of input data that includes simulated real‐time displacements from synthetic earthquakes and GNSS recordings from recent earthquakes worldwide. Performance will be evaluated using metrics and standards consistent with those adopted for ShakeAlert overall. This initial assessment will guide method refinement and synthesis of the most successful features into a candidate geodetic algorithm for the ShakeAlert production system. In parallel, improvements to geodetic networks and streamlining approaches to data processing and exchange will ensure robust geodetic data availability in the event of an earthquake.

West coast↗

Deep long-period earthquakes west of the volcanic arc in Oregon: evidence of serpentine dehydration in the fore-arc mantle wedge

Here we report on deep long-period earthquakes (DLPs) newly observed in four places in western Oregon. The DLPs are noteworthy for their location within the subduction fore arc: 40–80 km west of the volcanic arc, well above the slab, and near the Moho. These “offset DLPs” occur near the top of the inferred stagnant mantle wedge, which is likely to be serpentinized and cold. The lack of fore-arc DLPs elsewhere along the arc suggests that localized heating may be dehydrating the serpentinized mantle wedge at these latitudes and causing DLPs by dehydration embrittlement. Higher heat flow in this region could be introduced by anomalously hot mantle, associated with the western migration of volcanism across the High Lava Plains of eastern Oregon, entrained in the corner flow proximal to the mantle wedge. Alternatively, fluids rising from the subducting slab through the mantle wedge may be the source of offset DLPs. As far as we know, these are among the first DLPs to be observed in the fore arc of a subduction-zone system.

Oregon↗

Compilation of surface creep on California faults and comparison of WGCEP 2007 deformation model to Pacific-North American plate Mmtion

This Appendix contains 3 sections that 1) documents published observations of surface creep on California faults, 2) constructs line integrals across the WG-07 deformation model to compare to the Pacific - North America plate motion, and 3) constructs strain tensors of volumes across the WG-07 deformation model to compare to the Pacific - North America plate motion. Observation of creep on faults is a critical part of our earthquake rupture model because if a fault is observed to creep the moment released as earthquakes is reduced from what would be inferred directly from the fault's slip rate. There is considerable debate about how representative creep measured at the surface during a short time period is of the whole fault surface through the entire seismic cycle (e.g. Hudnut and Clark, 1989). Observationally, it is clear that the amount of creep varies spatially and temporally on a fault. However, from a practical point of view a single creep rate is associated with a fault section and the reduction in seismic moment generated by the fault is accommodated in seismic hazard models by reducing the surface area that generates earthquakes or by reducing the slip rate that is converted into seismic energy. WG-07 decided to follow the practice of past Working Groups and the National Seismic Hazard Map and used creep rate (where it was judged to be interseismic, see Table P1) to reduce the area of the fault surface that generates seismic events. In addition to following past practice, this decision allowed the Working Group to use a reduction of slip rate as a separate factor to accommodate aftershocks, post seismic slip, possible aseismic permanent deformation along fault zones and other processes that are inferred to affect the entire surface area of a fault, and thus are better modeled as a reduction in slip rate. C-zones are also handled by a reduction in slip rate, because they are inferred to include regions of widely distributed shear that is not completely expressed as earthquakes large enough to model. Because the ratio of the rate of creep relative to the total slip rate is often used to infer the average depth of creep, the depth of creep can be calculated and used to reduce the surface area of a fault that generates earthquakes in our model. This reduction of surface area of rupture is described by an aseismicity factor, assigned to each creeping fault in Appendix A. An aseismicity factor of less than 1 is only assigned to faults that are inferred to creep during the entire interseismic period. A single aseismicity factor was chosen for each section of the fault that creeps by expert opinion from the observations documented here. Uncertainties were not determined for the aseismicity factor, and thus it represents an unmodeled (and difficult to model) source of error. This Appendix simply provides the documentation of known creep, the type and precision of its measurement, and attempts to characterize the creep as interseismic, afterslip, transient or triggered. Parts 2 and 3 of this Appendix compare the WG-07 deformation model and the seismic source model it generates to the strain generated by the Pacific - North American plate motion. The concept is that plate motion generates essentially all of the elastic strain in the vicinity of the plate boundary that can be released as earthquakes. Adding up the slip rates on faults and all others sources of deformation (such as C-zones and distributed background seismicity) should approximately yield the plate motion. This addition is usually accomplished by one of four approaches: 1) line integrals that sum deformation along discrete paths through the deforming zone between the two plates, 2) seismic moment tensors that add up seismic moment of a representative set of earthquakes generated by a crustal volume spanning the plate boundary, 3) strain tensors generated by adding up the strain associated with all of the faults in a crustal volume spanning the plate

Open-File Report↗

Instrumentation recommendations for volcano monitoring at U.S. volcanoes under the national volcano early warning system

As magma moves toward the surface, it interacts with anything in its path: hydrothermal systems, cooling magma bodies from previous eruptions, and (or) the surrounding “country rock.” Magma also undergoes significant changes in its physical properties as pressure and temperature conditions change along its path. These interactions and changes lead to a range of geophysical and geochemical phenomena. The goal of volcano monitoring is to detect and correctly interpret such phenomena in order to provide early and accurate warnings of impending eruptions. Given the well-documented hazards posed by volcanoes to both ground-based populations (for example, Blong, 1984; Scott, 1989) and aviation (for example, Neal and others, 1997; Miller and Casadevall, 2000), volcano monitoring is critical for public safety and hazard mitigation. Only with adequate monitoring systems in place can volcano observatories provide accurate and timely forecasts and alerts of possible eruptive activity. At most U.S. volcanoes, observatories traditionally have employed a two-component approach to volcano monitoring: (1) install instrumentation sufficient to detect unrest at volcanic systems likely to erupt in the not-too-distant future; and (2) once unrest is detected, install any instrumentation needed for eruption prediction and monitoring. This reactive approach is problematic, however, for two reasons. 1. At many volcanoes, rapid installation of new ground-1. based instruments is difficult or impossible. Factors that complicate rapid response include (a) eruptions that are preceded by short (hours to days) precursory sequences of geophysical and (or) geochemical activity, as occurred at Mount Redoubt (Alaska) in 1989 (24 hours), Anatahan (Mariana Islands) in 2003 (6 hours), and Mount St. Helens (Washington) in 1980 and 2004 (7 and 8 days, respectively); (b) inclement weather conditions, which may prohibit installation of new equipment for days, weeks, or even months, particularly at midlatitude or high-latitude volcanoes; (c) safety factors during unrest, which can limit where new instrumentation can safely be installed (particularly at near-vent sites that can be critical for precursor detection and eruption forecasting); and (d) the remoteness of many U.S. volcanoes (particularly those in the Aleutians and the Marianas Islands), where access is difficult or impossible most of the year. Given these difficulties, it is reasonable to anticipate that ground-based monitoring of eruptions at U.S. volcanoes will likely be performed primarily with instruments installed before unrest begins. 2. Given a growing awareness of previously undetected 2. phenomena that may occur before an eruption begins, at present the types and (or) density of instruments in use at most U.S. volcanoes is insufficient to provide reliable early warning of volcanic eruptions. As shown by the gap analysis of Ewert and others (2005), a number of U.S. volcanoes lack even rudimentary monitoring. At those volcanic systems with monitoring instrumentation in place, only a few types of phenomena can be tracked in near-real time, principally changes in seismicity, deformation, and large-scale changes in thermal flux (through satellite-based remote sensing). Furthermore, researchers employing technologically advanced instrumentation at volcanoes around the world starting in the 1990s have shown that subtle and previously undetectable phenomena can precede or accompany eruptions. Detection of such phenomena would greatly improve the ability of U.S. volcano observatories to provide accurate early warnings of impending eruptions, and is a critical capability particularly at the very high-threat volcanoes identified by Ewert and others (2005). For these two reasons, change from a reactive to a proactive volcano-monitoring strategy is clearly needed at U.S. volcanoes. Monitoring capabilities need to be expanded at virtually every volcanic center, regardless of its current state of unrest, with particular emphases on real-time data transmission and increasing the diversity, quality, and quantity of instrumentation at U.S. volcano observatories. In this report, we present recommendations for the types and numbers of instruments that should be deployed to monitor U.S. volcanoes. These recommendations are the result of discussions among members of a panel of government- and university-based scientists about the status and future directions of volcano monitoring, with discussions framed by the U.S. Geological Survey (USGS)’s National Volcano Early Warning System (NVEWS) initiative (Ewert and others, 2005).

Scientific Investigations Report↗