Geology topics

Kurt L. Feigl

Publications and source records attributed to Kurt L. Feigl.

5 recordsLinked to original sources

USGS2024

Seismicity zoning at Coso geothermal field and stress changes from fluid production and migration

The Coso geothermal field is a major geothermal power production site in the western United States. It has been observed that low-magnitude seismic events ( M < 3.71) are unevenly distributed in three distinct zones, namely, nearfield (<3 km), midfield (3–6 km), and farfield (> 6 km) from the Coso geothermal plant. These zones exhibit distinct changes in earthquake location before and during geothermal production episodes that began in 1986. After 1986, the midfield region of the main flank experiences a significantly lower seismicity rate than the surrounding areas before production episodes. During 2014–2019, the farfield earthquakes cluster in the eastern and western parts of the greater Coso area, which is discernably different from how those pre-production earthquake events were distributed along the conjugate NW-SE and SW-NW trending structures across the main flank. Here, we analyze the stage of stress with finite-element-based poroelastic simulations to illustrate how the spatiotemporal evolution of the seismicity is associated with the pattern of stress perturbations caused by fluid migration amid the operations of geothermal power plants. Generally, ∼70% of co-production seismicity is found in zones of increased Coulomb stress between 2014 and 2019 at >99% confidence. Meanwhile, the midfield zone of seismic paucity overlaps with the zone of decreasing pore-fluid pressure. Overall, the results provide a physical explanation of how decadal geothermal operations at Coso have perturbed stress-field changes and contributed to the evolving characteristic seismic pattern, shedding insights into assessing the seismic hazard in other geothermal settings.

California↗

USGS2018

Geomorphic expression of rapid Holocene silicic magma reservoir growth beneath Laguna del Maule, Chile

Large rhyolitic volcanoes pose a hazard, yet the processes and signals foretelling an eruption are obscure. Satellite geodesy has revealed surface inflation signaling unrest within magma reservoirs underlying a few rhyolitic volcanoes. Although seismic, electrical, and potential field methods may illuminate the current configuration and state of these reservoirs, they cannot fully address the processes by which they grow and evolve on geologic time scales. We combine measurement of a deformed paleoshore surface, isotopic dating of volcanism and surface exposure, and modeling to determine the rate of growth of a rhyolite-producing magma reservoir. The numerical approach builds on a magma intrusion model developed to explain the current, decade-long, surface inflation at >20 cm/year. Assuming that the observed 62-m uplift reflects several non-eruptive intrusions of magma, each similar to the unrest over the past decade, we find that ~13 km 3 of magma recharged the reservoir at a depth of ~7 km during the Holocene, accompanied by the eruption of ~9 km 3 of rhyolite. The long-term rate of magma input is consistent with reservoir freezing and pluton formation. Yet, the unique set of observations considered here implies that large reservoirs can be incubated and grow at shallow depth via episodic high-flux magma injections. These replenishment episodes likely drive rapid inflation, destabilize cooling systems, propel rhyolitic eruptions, and thus should be carefully monitored.

Laguna del Maule↗

USGS2016

Volcano deformation source parameters estimated from InSAR: Sensitivities to uncertainties in seismic tomography

The eruption cycle of a volcano is controlled in part by the upward migration of magma. The characteristics of the magma flux produce a deformation signature at the Earth's surface. Inverse analyses use geodetic data to estimate strategic controlling parameters that describe the position and pressurization of a magma chamber at depth. The specific distribution of material properties controls how observed surface deformation translates to source parameter estimates. Seismic tomography models describe the spatial distributions of material properties that are necessary for accurate models of volcano deformation. This study investigates how uncertainties in seismic tomography models propagate into variations in the estimates of volcano deformation source parameters inverted from geodetic data. We conduct finite element model-based nonlinear inverse analyses of interferometric synthetic aperture radar (InSAR) data for Okmok volcano, Alaska, as an example. We then analyze the estimated parameters and their uncertainties to characterize the magma chamber. Analyses are performed separately for models simulating a pressurized chamber embedded in a homogeneous domain as well as for a domain having a heterogeneous distribution of material properties according to seismic tomography. The estimated depth of the source is sensitive to the distribution of material properties. The estimated depths for the homogeneous and heterogeneous domains are 2666 ± 42 and 3527 ± 56 m below mean sea level, respectively (99% confidence). A Monte Carlo analysis indicates that uncertainties of the seismic tomography cannot account for this discrepancy at the 99% confidence level. Accounting for the spatial distribution of elastic properties according to seismic tomography significantly improves the fit of the deformation model predictions and significantly influences estimates for parameters that describe the location of a pressurized magma chamber.

Journal of Geophysical Research B: Solid Earth↗

USGS2014

Dynamics of a large, restless, rhyolitic magma system at Laguna del Maule, southern Andes, Chile

Explosive eruptions of large-volume rhyolitic magma systems are common in the geologic record and pose a major potential threat to society. Unlike other natural hazards, such as earthquakes and tsunamis, a large rhyolitic volcano may provide warning signs long before a caldera-forming eruption occurs. Yet, these signs—and what they imply about magma-crust dynamics—are not well known. This is because we have learned how these systems form, grow, and erupt mainly from the study of ash flow tuffs deposited tens to hundreds of thousands of years ago or more, or from the geophysical imaging of the unerupted portions of the reservoirs beneath the associated calderas. The Laguna del Maule Volcanic Field, Chile, includes an unusually large and recent concentration of silicic eruptions. Since 2007, the crust there has been inflating at an astonishing rate of at least 25 cm/yr. This unique opportunity to investigate the dynamics of a large rhyolitic system while magma migration, reservoir growth, and crustal deformation are actively under way is stimulating a new international collaboration. Findings thus far lead to the hypothesis that the silicic vents have tapped an extensive layer of crystal-poor, rhyolitic melt that began to form atop a magmatic mush zone that was established by ca. 20 ka with a renewed phase of rhyolite eruptions during the Holocene. Modeling of surface deformation, magnetotelluric data, and gravity changes suggest that magma is currently intruding at a depth of ~5 km. The next phase of this investigation seeks to enlarge the sets of geophysical and geochemical data and to use these observations in numerical models of system dynamics.

Laguna del Maule Volcanic Field↗

USGS2006

Geodetic observations of post-seismic transients in the context of the earthquake deformation cycle

Satellite geodetic techniques that can measure displacements with millimeter-level accuracy reveal transient signals in the deformation fields produced by both moderate and large earthquakes. These post-seismic signals exhibit characteristic time scales ranging from weeks to decades and distance scales from hundreds of meters to hundreds of kilometers. By considering them in the context of the earthquake deformation cycle, we can test hypotheses about the processes driving them and constrain the rheology of the lithosphere. We discuss three broad categories of mechanism: afterslip in the plane of the co-seismic rupture (analogous to a rubber eraser), fluid flow in the fault zone (analogous to a water-laden sponge), and ductile flow in a weak substrate (analogous to a pot of honey).

Comptes Rendus - Geoscience↗