On 8 August 2020, northwest North Carolina experienced an Mw 5.1 earthquake that caused damage to buildings and roads in the city of Sparta. A regional centroid moment tensor solution shows that the earthquake was the result of slip on a reverse fault with a minor strike-slip component. Interferometric Synthetic Aperture Radar (InSAR) data, from the Japan Aerospace Exploration Agency’s Advanced Land Observing Satellite #2 (ALOS2) satellite, reveal a deformation field that is more complex than expected from a single reverse fault earthquake. The data also reveal an apparent fault rupture at the Earth’s surface that caused damage to local roads. Modeling of the InSAR deformation field indicates the fault rupture is associated with a very shallow normal faulting event with an equivalent Mw of about 5.1 that overprinted the reverse fault deformation field and possibly occurred aseismically.

ABSTRACT The 2019 Ridgecrest, California, earthquake sequence involved predominantly right-lateral strike slip on a northwest–southeast-trending subvertical fault in the 6 July M 7.1 mainshock, preceded by left-lateral strike slip on a northeast–southwest-trending subvertical fault in the 4 July M 6.4 foreshock. To characterize the postseismic deformation, we assemble displacements measured by Global Positioning System (GPS) and Interferometric Synthetic Aperture Radar. The geodetic measurements illuminate vigorous postseismic deformation for at least 21 months following the earthquake sequence. The postseismic transient deformation is particularly well constrained from survey-mode GPS (sGPS) in the epicentral region carried out during the weeks after the mainshock. We interpret these observations with mechanical models including afterslip and viscoelastic relaxation of the lower crust and mantle asthenosphere. During the first 21 months, up to several centimeters of horizontal motions are measured at continuous GPS and sGPS sites, with amplitude that diminishes slowly with distance from the mainshock rupture, suggestive of deeper afterslip or viscoelastic relaxation. We find that although afterslip involving right-lateral strike slip along the mainshock fault traces and their deeper extensions reach a few decimeters, most postseismic deformation is attributable to viscoelastic relaxation of the lower crust and mantle. Within the Basin and Range crust and mantle, we infer a transient lower crust viscosity several times that of the mantle asthenosphere. The transient mantle asthenosphere viscosity is ∼1.3×1017 Pa s, and the adjacent Central Valley transient mantle asthenosphere viscosity is ∼7×1017 Pa s, about five times higher and consistent with an asymmetry in postseismic horizontal motions across the mainshock surface rupture.

Crustal inflation near the Three Sisters volcanic center documented since the mid‐1990s has persisted for more than two decades. We update past analyses of the event through 2020 by simultaneously inverting InSAR interferograms, GPS time series, and leveling data for time‐dependent volcanic deformation source parameters. We explore several source models to estimate how the deformation rate varied through time and to identify parameters that can reproduce measured deformation. Our preferred model is a Mogi source 4.1 km below sea level (5.9 km below the surface) about 5 km west of the summit of South Sister. Inflation started in late 1995 or 1996; the rate increased rapidly during 1998–1999, and peaked in late 1999, resulting in maximum surface uplift of about 30 cm by mid‐2020. Since 2000, the inflation rate generally declined exponentially with a time constant of about 6 years. Two source inflation scenarios fit the data equally well. In the first, the crust surrounding the source is elastic and the net source‐volume increase, which we attribute to persistent magma input, has been about 49 × 10 6 m 3 . The second scenario adds a viscoelastic shell surrounding the Mogi source. In that case, an injection of about 21 × 10 6 m 3 of magma prior to 2000, followed by continuing relaxation of the viscoelastic shell, can account for most of the observed surface deformation. In both scenarios, modeling reveals quasiperiodic increases in the inflation rate (pulses) with a recurrence interval of 3–4 years, both before and after 2000.

The 2020 M 6.5 Stanley, Idaho, earthquake produced rupture in the north of the active Sawtooth fault in the northern basin and range at depth, without any observable surface rupture. Global Positioning System (GPS) and Interferometric Synthetic Aperture Radar (InSAR) data yield several millimeters of static offsets out to ∼100  km from the rupture and up to ∼0.1  m of near-field crustal deformation. We combine the GPS and InSAR data with long-period regional seismic waveforms to derive models of kinematic slip and afterslip. We find that the coseismic rupture is complex, likely involving up to 2 m combined left-lateral strike slip and normal slip on a previously unidentified ∼south-southeast-striking fault. This slip is predominantly left-lateral strike slip, different from the dominant east-northeast–west-northwest normal faulting of the region. At least one ∼northeast-trending fault, likely associated with the Trans-Challis fault system, is inferred to have accommodated a few decimeters of right-lateral afterslip, consistent with vigorous aftershock activity at depth along northeast-trending lineations.

The 2020 Magna, Utah, earthquake produced observable crustal deformation over an similar to 100 km(2) area around the southeast margin of Great Salt Lake, but it did not produce any surface rupture. To obtain a detailed picture of the fault slip, we combine strong-motion seismic waveforms with Global Positioning System static offsets and Interferometric Synthetic Aperture Radar observations to obtain kinematic and static slip models of the event. We sample the regional seismic wavefield with three-component records from 68 stations of the University of Utah Seismograph Stations network. We find that coseismic slip and afterslip, with predominantly normal slip, distributed on a shallowly west-dipping plane, possibly augmented by afterslip on a steeply northeast-dipping plane, best fits the joint dataset. The west-dipping plane locates near previously inferred sources of interseismic creep at depth. Hence, the earthquake may have occurred on the down-dip extension of the Wasatch fault and activated further slip (afterslip) at shallow depth east of the hypocenter. This inferred afterslip may have driven the vigorous aftershock activity that was concentrated east of the hypocenter.

On March 26, 2020, a M 5.0 earthquake occurred in the Delaware Basin, Texas, near the border between Reeves and Culberson Counties. This was the third largest earthquake recorded in Texas and the largest earthquake in the Central and Eastern United States since the three M 5.0-5.8 induced events in Oklahoma during 2016. Using multistation waveform template matching, we detect 3,940 earthquakes in the sequence with the first event in the area occurring in May 2018. The M 5.0 earthquake sequence occurred on a ENE (similar to 082 degrees) normal fault dipping similar to 37 degrees toward the south. The earthquake caused 6 mm of oblique surface deformation, and geodetic slip inversion suggests slip was isolated above 6 km depth. We find that the sequence was most likely induced by nearby wastewater disposal operations, and seismicity rates in the region surrounding the M 5.0 will likely continue to increase in the future if disposal operations continue unaltered.

ABSTRACT The 2019 Ridgecrest, California, earthquake sequence produced observable crustal deformation over much of central and southern California, as well as surface rupture over several tens of kilometers. To obtain a detailed picture of the fault slip involved in the 4 July M 6.4 foreshock and 6 July M 7.1 mainshock, we combine strong-motion seismic waveforms with crustal deformation observations to obtain kinematic and static slip models of both events. We sample the regional seismic wavefield for both the foreshock and mainshock with three-component records from 31 stations of the California Integrated Seismic Network. The deformation observations include Global Positioning System (GPS), Interferometric Synthetic Aperture Radar (InSAR), and borehole strainmeter recordings of the dynamic strain field. These data collectively constrain the kinematic coseismic slip distributions of the events, with measurements variously observing coseismic slip from one event (e.g., seismic waveforms, kinematic solutions from continuous GPS, and strainmeter time series) or coseismic slip from both events combined (InSAR). We find that the foreshock ruptured two separate faults, one with left-lateral strike slip on a northeast–southwest-trending fault and the other with right-lateral strike slip on an orthogonal fault, with unilateral rupture propagation along both. The mainshock ruptured a series of northwest–southeast-trending faults with right-lateral strike slip concentrated in the uppermost 6 km with exceptionally low-rupture velocity averaging 1.0–1.5 km/s. A possible explanation for the low-rupture velocity is that the mainshock rupture expended relatively high energy, generating secondary fractures in off-fault deformation, which is consistent with field and seismic evidence of plastic deformation on small fault strands adjacent to the main rupture trace.

Recent activity has provided new insights into the causes of surface deformation in and around the Yellowstone Caldera, a topic that has been debated since the discovery of caldera floor uplift more than four decades ago. An episode of unusually rapid uplift (>15 cm/yr) centered near Norris Geyser Basin along the north caldera rim began in late 2013 and continued until a M w 4.9 earthquake on 30 March 2014; thereafter, uplift abruptly switched to subsidence. Uplift at rates of several centimeters per year resumed in 2016 and continued at least through the end of 2018. Modeling of Global Positioning System and interferometric synthetic aperture radar data suggests an evolving process of deep magma intrusion during 1996–2001 followed by volatile ascent and accumulation at shallow levels, perhaps as shallow as a few hundred meters depth. The depth of shallow volatile accumulation appears to have shallowed from the 2014 to the 2016 deformation episode, and frequent eruptions of Steamboat Geyser since March 2018 are likely a surface manifestation of this ongoing process. Hydrothermal explosion features are prominent in the Norris Geyser Basin area, and the apparent shallow nature of the volatile accumulation implies an increased risk of hydrothermal explosions.

We employ strong motion seismograms and static offsets from the Global Positioning System, Interferometric Synthetic Aperture Radar, and other measurements in order to derive a coseismic slip and afterslip model of the M6.0 24 August 2014 South Napa earthquake. This earthquake ruptured an similar to 13-km-long portion of the West Napa fault with predominantly right-lateral strike slip. In the kinematic seismic slip inversions, we couple the coseismic slip and afterslip distributions by requiring both distributions to involve right-lateral strike slip with positive amplitude, with the net static slip being the sum of the two. We consider several candidate fault geometries: a first involving two steeply east dipping fault planes that reach Earth's surface at the western surface trace (STW), where most surface rupture was observed, a second involving a steeply west dipping plane that also reaches Earth's surface at the STW, and a third involving a combination of two variably west dipping planes constrained to pass through the locus of postseismic seismicity located similar to 1 km west of the STW. The data are best fit using the model of two east dipping fault planes, with coseismic slip up to similar to 1.2 m on a dominant shallow asperity about 10 km north of the hypocenter and on deeper asperities on the southern part of the rupture. Afterslip up to 1 m is concentrated along the southern part of the rupture at depths less than or similar to 5 km, consistent with surface observations of afterslip. Seismic moments associated with coseismic slip and afterslip are 1.13 x 10(18) N m (Mw 6.00) and 3.64 x 10(17) N m, respectively.

The 2 September 2017 M 5.3 Sulphur Peak, Idaho, earthquake is one of the largest earthquakes in southern Idaho since the 1983 M 6.9 Borah Peak earthquake. It was followed by a vigorous aftershock sequence for nearly two weeks that included five events above M 4.5. The coseismic and early postseismic deformation was measured with both Interferometric Synthetic Aperture Radar and Global Positioning System (GPS), yielding up to 3 cm subsidence southwest of the mainshock epicenter and horizontal motions of similar to 1 mm at sites similar to 40 km east and west of the epicenter. We derive dislocation models of the net slip during the similar to 14-day swarm from Sentinel 1A interferograms and GPS offsets, allowing for both fault-zone collapse and normal faulting to account for the observed geodetic motions. Slip inversions yield several decimeters of normal slip on one or more normal faults near the mainshock hypocenter. Distributed normal slip on a moderately (55 degrees) east-dipping fault, normal slip on one or more shallowly west-dipping faults, or a combination thereof explain the data equally well and are difficult to distinguish from one another on the basis of geodetic data alone. Previously mapped regional Sevier-age thrust structures and later normal faults dip westward, suggesting that the sequence reactivated one or more ancient thrust structures with normal slip. If a moderately east-dipping fault accommodated substantial slip, it would imply a nascent fault structure that cuts across the reactivated ancient thrust structures. The inferred geodetic moment of 3.02-4.39 x 10(-17) N . m (M-w 5.62-5.73) greatly exceeds the 1.15 x 10(-17) N . m (M-w 5.34) seismic moment of the 2 September mainshock, showing that most of the moment release occurred during the aftershock sequence. Up to similar to 0.2 m of fault-zone collapse may have occurred on a shallow west-dipping fault, suggesting possible large-scale expulsion of fluids from the fault zone at depth.

The arrival of space-based imaging radar as a revolutionary land-surface mapping and monitoring tool little more than a quarter century ago enabled a spate of innovative volcano research worldwide. Soon after launch of European Space Agency’s ERS-1 spacecraft in 1991, the U.S. Geological Survey began SAR and InSAR studies of volcanoes in the Aleutian and Cascades arcs, in Hawaiʻi, and elsewhere in the western U.S. including the Yellowstone and Long Valley calderas. This paper summarizes results of that effort and presents new findings concerning: (1) prevalence of volcano deformation in the Aleutian and Cascade arcs; (2) surface-change detection and hazard assessment during eruptions at Aleutian and Hawaiian volcanoes; (3) geodetic imaging of magma storage and transport systems in Hawaiʻi; and (4) deformation sources and processes at the Yellowstone and Long Valley calderas. Surface deformation caused by a variety of processes is common in arc settings and could easily escape detection without systematic InSAR surveillance. Space-based SAR imaging of active lava flows and domes in remote or heavily vegetated settings, including during periods of bad weather and darkness, extends land-based monitoring capabilities and improves hazards assessments. At Kīlauea Volcano, comprehensive SAR and InSAR observations identify multiple magma storage zones beneath the summit area and along the East Rift Zone, and illuminate magma transport pathways. The same approach at Yellowstone tracks the ascent of magmatic volatiles from a mid-crustal intrusion to shallow depth and relates that process to increased hydrothermal activity at the surface. Together with recent and planned launches of highly capable imaging-radar satellites, these findings support an optimistic outlook for near-real time surveillance of volcanoes at global scale in the coming decade.

The 3 September 2016 M-w 5.8 Pawnee earthquake in northern Oklahoma is the largest earthquake ever recorded in Oklahoma. The coseismic deformation was measured with both Interferometric Synthetic Aperture Radar and Global Positioning System (GPS), with measureable signals of order 1 cm and 1 mm, respectively. We derive a coseismic slip model from Sentinel-1A and Radarsat 2 interferograms and GPS static offsets, dominated by distributed left-lateral strike slip on a primary west-northwest-east-southeast-trending subvertical plane, whereas strike slip is concentrated near the hypocenter (5.6 km depth), with maximum slip of similar to 1 m located slightly east and down-dip of the hypocenter. Based on systematic misfits of observed interferogram line-of-sight (LoS) displacements, with LoS based on shear-dislocation models, a few decimeters of fault-zone collapse are inferred in the hypocentral region where coseismic slip was the largest. This may represent the postseismic migration of large volumes of fluid away from the high-slip areas, made possible by the creation of a temporary high-permeability damage zone around the fault.

Semipermanent GPS (SPGPS) is an alternative to conventional campaign or survey-mode GPS (SGPS) and to continuous GPS (CGPS) that offers several advantages for monitoring ground deformation. Unlike CGPS installations, SPGPS stations can be deployed quickly in response to changing volcanic conditions or earthquake activity such as a swarm or aftershock sequence. SPGPS networks can be more focused or more extensive than CGPS installations, because SPGPS equipment can be moved from station to station quickly to increase the total number of stations observed in a given time period. SPGPS networks are less intrusive on the landscape than CGPS installations, which makes it easier to satisfy land-use restrictions in ecologically sensitive areas. SPGPS observations are preferred over SGPS measurements because they provide better precision with only a modest increase in the amount of time, equipment, and personnel required in the field. We describe three applications of the SPGPS method that demonstrate its utility and flexibility. At the Yellowstone caldera, Wyoming, a 9-station SPGPS network serves to densify larger preexisting networks of CGPS and SGPS stations. At the Three Sisters volcanic center, Oregon, a 14-station SPGPS network complements an SGPS network and extends the geographic coverage provided by 3 CGPS stations permitted under wilderness land-use restrictions. In the Basin and Range province in northwest Nevada, a 6-station SPGPS network has been established in response to a prolonged earthquake swarm in an area with only sparse preexisting geodetic coverage. At Three Sisters, the estimated precision of station velocities based on annual ~3month summertime SPGPS occupations from 2009 to 2015 is approximately half that for nearby CGPS stations. Conversely, SPGPS-derived station velocities are about twice as precise as those based on annual ~1week SGPS measurements. After 5years of SPGPS observations at Three Sisters, the precision of velocity determinations is estimated to be 0.5mm/yr in longitude, 0.6mm/yr in latitude, and 0.8mm/yr in height. We conclude that an optimal approach to monitoring volcano deformation includes complementary CGPS and SPGPS networks, periodic InSAR observations, and measurements from in situ borehole sensors such as tiltmeters or strainmeters. This comprehensive approach provides the spatial and temporal detail necessary to adequately characterize a complex and evolving deformation pattern. Such information is essential to multi-parameter models of magmatic or tectonic processes that can help to guide research efforts, and also to inform hazards assessments and land-use planning decisions.

In 2001, a rare swarm of small, shallow earthquakes beneath the city of Spokane, Washington, caused ground shaking as well as audible booms over a five-month period. Subsequent Interferometric Synthetic Aperture Radar (InSAR) data analysis revealed an area of surface uplift in the vicinity of the earthquake swarm. To investigate the potential faults that may have caused both the earthquakes and the topographic uplift, we collected similar to 3 km of high-resolution seismic-reflection profiles to image the upper-source region of the swarm. The two profiles reveal a complex deformational pattern within Quaternary alluvial, fluvial, and flood deposits, underlain by Tertiary basalts and basin sediments. At least 100 m of arching on a basalt surface in the upper 500 m is interpreted from both the seismic profiles and magnetic modeling. Two west-dipping faults deform Quaternary sediments and project to the surface near the location of the Spokane fault defined from modeling of the InSAR data.