Geology topics

Charles Wicks

Publications and source records attributed to Charles Wicks.

At least 19 recordsLinked to original sources

USGS2022

Surface rupture on a secondary fault associated with the August 8, 2020, Mw 5.1 Sparta North Carolina Earthquake

On August 8, 2020 northwest North Carolina experienced a M w 5.1 earthquake that caused damage to buildings and roads in the city of Sparta. A regional centroid moment tensor solution shows the earthquake was the result of slip on a reverse fault with a minor strike-slip component. InSAR data, from the Japan Aerospace Exploration Agency’s 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 M w of about 5.1, that overprinted the reverse fault deformation field and possibly occurred aseismically.

North Carolina↗

USGS2022

Kinematic slip model of the July 8, 2021 M6.0 Antelope Valley, California, earthquake

We present a kinematic slip model of the July 8, 2021 Antelope Valley earthquake from a finite-source inversion based on regional seismic waveforms and static offsets from GPS and InSAR. Seismic waveforms are employed at 6s dominant period out to 100 km from the epicenter, and the combined GPS and InSAR datasets cover the near field and far field out to ∼ 100 km and constrain the overall rupture size. The aftershock pattern defines a nearly north-striking, 50◦ east-dipping fault plane. We find a unilateral rupture along this fault plane propagating southward and updip with predominantly normal slip up to ∼ 1.5m. The estimated seismic moment of 8.47 × 10 22 17 Nm is equivalent to Mw 5.92. A finite-source inversion that retains seismic waveforms and GPS static offsets but omits InSAR range changes yields a seismic moment of 1.08 × 10 25 18 Nm (Mw 5.99). Despite vigorous aftershock activity between 10 km and Earth’s surface, coseismic slip is concentrated in the depth interval 7 - 10 km.

California↗

USGS2021

Geodetic constraints on a 25-year magmatic inflation episode near Three Sisters, central Oregon

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.

Oregon↗

USGS2021

Post-seismic relaxation following the 2019 Ridgecrest, California, earthquake sequence

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 × 10 17    Pa   s ">

California↗

USGS2021

Coseismic fault slip and afterslip associated with the M5.7 March 18, 2020 Magna, Utah, earthquake

The 2020 Magna, Utah, earthquake produced observable crustal deformation over a ∼ 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 GPS static oﬀsets and Interferometric Synthetic Aperture Radar (InSAR) observations to obtain kinematic and static slip models of the event. We sample the regional seismic waveﬁeld with 3-component records from 68 stations of the University of Utah Seismograph Stations network. We ﬁnd 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 ﬁts 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 downdip ex-tension 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.

Utah↗

USGS2020

The induced Mw 5.0 March 2020 west Texas seismic sequence

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 (∼082°) normal fault dipping ∼37° 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.

Texas↗

USGS2020

Rupture process of the M6.5 Stanley, Idaho, earthquake inferred from seismic waveform and geodetic data

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.

Idaho↗

USGS2020

Kinematics of fault slip associated with the July 4-6 2019 Ridgecrest, Californai earthquakes sequence

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 ">

California↗

USGS2020

Magma intrusion and volatile ascent beneath Norris Geyser Basin, Yellowstone National Park

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.

Wyoming↗

USGS2019

Fault slip associated with the 2 September 2017 M 5.3 Sulphur Peak, Idaho, earthquake and aftershock sequence

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 ∼ 1    mm "> ∼ 1 mm at sites ∼ 40    km "> ∼ 40 km east and west of the epicenter. We derive dislocation models of the net slip during the ∼ 14 ‐ day "> ∼ 14 ‐ day s warm 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°) 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 × 10 17    N · m "> 3.02 – 4.39 × 10 17 N ⋅ m ( ⁠ M w "> M w 5.62–5.73) greatly exceeds the 1.15 × 10 17    N · m "> 1.15×10 17 N·m ( ⁠ M w "> 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 ∼ 0.2    m "> ∼ 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.

Idaho↗

USGS2011

InSAR observations of aseismic slip associated with an earthquake swarm in the Columbia River flood basalts

In 2009 a swarm of small shallow earthquakes occurred within the basalt flows of the Columbia River Basalt Group (CRBG). The swarm occurred within a dense seismic network in the U.S. Department of Energys Hanford Site. Data from the seismic network along with interferometric synthetic aperture radar (InSAR) data from the European Space Agencys (ESA) ENVISAT satellite provide insight into the nature of the swarm. By modeling the InSAR deformation data we constructed a model that consists of a shallow thrust fault and a near horizontal fault. We suggest that the near horizontal lying fault is a bedding-plane fault located between basalt flows. The geodetic moment of the modeled fault system is about eight times the cumulative seismic moment of the swarm. Precise location estimates of the swarm earthquakes indicate that the area of highest slip on the thrust fault, ???70mm of slip less than ???0.5km depth, was not located within the swarm cluster. Most of the slip on the faults appears to have progressed aseismically and we suggest that interbed sediments play a central role in the slip process. Copyright 2011 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2011

The role of dyking and fault control in the rapid onset of eruption at Chaitén Volcano, Chile

Rhyolite is the most viscous of liquid magmas, so it was surprising that on 2 May 2008 at Chaitén Volcano, located in Chile’s southern Andean volcanic zone, rhyolitic magma migrated from more than 5 km depth in less than 4 hours and erupted explosively with only two days of detected precursory seismic activity . The last major rhyolite eruption before that at Chaitén was the largest volcanic eruption in the twentieth century, at Novarupta volcano, Alaska, in 1912. Because of the historically rare and explosive nature of rhyolite eruptions and because of the surprisingly short warning before the eruption of the Chaitén volcano, any information about the workings of the magmatic system at Chaitén, and rhyolitic systems in general, is important from both the scientific and hazard perspectives. Here we present surface deformation data related to the Chaitén eruption based on radar interferometry observations from the Japan Aerospace Exploration Agency (JAXA) DAICHI (ALOS) satellite. The data on this explosive rhyolite eruption indicate that the rapid ascent of rhyolite occurred through dyking and that melt segregation and magma storage were controlled by existing faults.

Chaitén Volcano↗

USGS2010

Ground surface deformation patterns, magma supply, and magma storage at Okmok volcano, Alaska, from InSAR analysis: 1. Intereruption deformation, 1997–2008

Starting soon after the 1997 eruption at Okmok volcano and continuing until the start of the 2008 eruption, magma accumulated in a storage zone centered ~3.5 km beneath the caldera floor at a rate that varied with time. A Mogi-type point pressure source or finite sphere with a radius of 1 km provides an adequate fit to the deformation field portrayed in time-sequential interferometric synthetic aperture radar images. From the end of the 1997 eruption through summer 2004, magma storage increased by 3.2–4.5 × 10 7 m 3 , which corresponds to 75–85% of the magma volume erupted in 1997. Thereafter, the average magma supply rate decreased such that by 10 July 2008, 2 days before the start of the 2008 eruption, magma storage had increased by 3.7–5.2 × 10 7 m 3 or 85–100% of the 1997 eruption volume. We propose that the supply rate decreased in response to the diminishing pressure gradient between the shallow storage zone and a deeper magma source region. Eventually the effects of continuing magma supply and vesiculation of stored magma caused a critical pressure threshold to be exceeded, triggering the 2008 eruption. A similar pattern of initially rapid inflation followed by oscillatory but generally slowing inflation was observed prior to the 1997 eruption. In both cases, withdrawal of magma during the eruptions depressurized the shallow storage zone, causing significant volcano-wide subsidence and initiating a new intereruption deformation cycle.

Alaska↗

USGS2010

Characterizing 6 August 2007 Crandall Canyon mine collapse from ALOS PALSAR InSAR

We used ALOS InSAR images to study land surface deformation over the Crandall Canyon mine in Utah, which collapsed on 6 August 2007 and killed six miners. The collapse was registered as a ML 3.9 seismic event. An InSAR image spanning the time of the collapse shows 25–30 cm surface subsidence over the mine. We used distributed dislocation sources to model the deformation ﬁeld, and found that a collapse source model alone does not adequately ﬁt the deformation ﬁeld. Normal faulting is also required, such that the event is best characterized as a ‘trapdoor’ collapse. The calculated moment of the normal fault is about the same as the moment of the collapse source, with each larger than the seismically computed moment. Our InSAR results, including the location of the event, the extent of the collapsed area, and constraints on the shearing component of the deformation source, all conﬁrm and extend recent seismic studies of the 6 August 2007 event.

Utah↗

USGS2009

Rupture parameters of the 2003 Zemmouri (M w 6.8), Algeria, earthquake from joint inversion of interferometric synthetic aperture radar, coastal uplift, and GPS

We study the surface deformation associated with the 21 May 2003 (M w = 6.8) Zemmouri (Algeria) earthquake, the strongest seismic event felt in the Algiers region since 1716. The thrust earthquake mechanism and related surface deformation revealed an average 0.50 m coastal uplift along ??55-km-long coastline. We obtain coseismic interferograms using Envisat advanced synthetic aperture radar (ASAR) (IS2) and RADARSAT standard beam (ST4) data from both the ascending and descending orbits of Envisat satellite, whereas the RADARSAT data proved useful only in the descending mode. While the two RADARSAT interferograms cover the earthquake area, Envisat data cover only the western half of the rupture zone. Although the interferometric synthetic aperture radar (InSAR) coherence in the epicenter area is poor, deformation fringes are observed along the coast in different patches. In the Boumerdes area, the maximum coseismic deformation is indicated by the high gradient of fringes visible in all interferograms in agreement with field measurements (tape, differential GPS, leveling, and GPS). To constrain the earthquake rupture parameters, we model the interferograms and uplift measurements using elastic dislocations on triangular fault patches in an elastic and homogeneous half-space. We invert the coseismic slip using first, a planar surface and second, a curved fault, both constructed from triangular elements using Poly3Dinv program that uses a damped least square minimization. The best fit of InSAR, coastal uplift, and GPS data corresponds to a 65-km-long fault rupture dipping 40?? to 50?? SE, located at 8 to 13 km offshore with a change in strike west of Boumerdes from N60??-65?? to N95??-105??. The inferred rupture geometry at depth correlates well with the seismological results and may have critical implications for the seismic hazard assessment of the Algiers region. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗

USGS2007

Accelerated uplift and magmatic intrusion of the Yellowstone caldera, 2004 to 2006

The Yellowstone caldera began a rapid episode of ground uplift in mid-2004, revealed by Global Positioning System and interferometric synthetic aperture radar measurements, at rates up to 7 centimeters per year, which is over three times faster than previously observed inflation rates. Source modeling of the deformation data suggests an expanding volcanic sill of ???1200 square kilometers at a 10-kilometer depth beneath the caldera, coincident with the top of a seismically imaged crustal magma chamber. The modeled rate of source volume increase is 0.1 cubic kilometer per year, similar to the amount of magma intrusion required to supply the observed high heat flow of the caldera. This evidence suggests magma recharge as the main mechanism for the accelerated uplift, although pressurization of magmatic fluids cannot be ruled out.

Wyoming↗

USGS2004

Constraining the slip distribution and fault geometry of the M w 7.9, 3 November 2002, Denali fault earthquake with Interferometric Synthetic Aperture Radar and Global Positioning System data

The M w 7.9, Denali fault earthquake ( DFE ) is the largest continental strike-slip earthquake to occur since the development of Interferometric Synthetic Aperture Radar (In SAR ). We use five interferograms, constructed using radar images from the Canadian Radarsat-1 satellite, to map the surface deformation at the western end of the fault rupture. Additional geodetic data are provided by displacements observed at 40 campaign and continuous Global Positioning System ( GPS ) sites. We use the data to determine the geometry of the Susitna Glacier fault, thrusting on which initiated the DFE , and to determine a slip model for the entire event that is consistent with both the In SAR and GPS data. We find there was an average of 7.3 ± 0.4 m slip on the Susitna Glacier fault, between 1 and 9.5 km depth on a 29 km long fault that dips north at 41 ± 0.7° and has a surface projection close to the mapped rupture. On the Denali fault, a simple model with large slip patches finds a maximum of 8.7 ± 0.7 m of slip between the surface and 14.3 ± 0.2 km depth. A more complex distributed slip model finds a peak of 12.5 ± 0.8 m in the upper 4 km, significantly higher than the observed surface slip. We estimate a geodetic moment of 670 ± 10 × 10 18 N m ( M w 7.9), consistent with seismic estimates. Lack of preseismic data resulted in an absence of In SAR coverage for the eastern half of the DFE rupture. A dedicated geodetic In SAR mission could obviate coverage problems in the future.

Bulletin of the Seismological Society of America↗

USGS2003

Source model for the Mw 6.7, 23 October 2002, Nenana Mountain earthquake (Alaska) from InSAR

The 23 October 2002 Nenana Mountain Earthquake (M w ∼ 6.7) occurred on the Denali Fault (Alaska), to the west of the M w ∼ 7.9 Denali Earthquake that ruptured the same fault 11 days later. We used 6 interferograms, constructed using radar images from the Canadian Radarsat-1 and European ERS-2 satellites, to determine the coseismic surface deformation and a source model. Data were acquired on ascending and descending satellite passes, with incidence angles between 23 and 45 degrees, and time intervals of 72 days or less. Modeling the event as dislocations in an elastic half space suggests that there was nearly 0.9 m of right-lateral strike-slip motion at depth, on a near-vertical fault, and that the maximum slip in the top 4 km of crust was less than 0.2 m. The Nenana Mountain Earthquake increased the Coulomb stress at the future hypocenter of the 3 November 2002, Denali Earthquake by 30–60 kPa.

Alaska↗