Geology topics

Ninfa Bennington

Publications and source records attributed to Ninfa Bennington.

5 recordsLinked to original sources

USGS2020

Linking magma storage and ascent to eruption volume and composition at an arc caldera

Conceptual models of magma storage and transport under calderas favor a connected system of sills and dikes. These features are individually below the resolution of standard seismic tomography, but radial seismic anisotropy can reveal where they exist in aggregate. We model radial anisotropy at Okmok caldera, Alaska, to demonstrate the presence of a caldera-centered stacked sill complex and surrounding dike system. We show that ascending magma, inferred from seismicity, either intersects the sill complex, resulting in a larger volume eruption of evolved magma, or bypasses the overlying sill complex via dikes, resulting in a low-volume mafic eruption. Our results exemplify how the locations of magma storage and paths of transport impact eruption size and composition. As this type of crustal storage is likely common to many calderas, this analysis offers a potential new framework for volcano observatories to forecast the size of impending eruptions.

Alaska↗

USGS2015

Monitoring changes in seismic velocity related to an ongoing rapid inflation event at Okmok volcano, Alaska

Okmok is one of the most active volcanoes in the Aleutian Arc. In an effort to improve our ability to detect precursory activity leading to eruption at Okmok, we monitor a recent, and possibly ongoing, GPS-inferred rapid inflation event at the volcano using ambient noise interferometry (ANI). Applying this method, we identify changes in seismic velocity outside of Okmok’s caldera, which are related to the hydrologic cycle. Within the caldera, we observe decreases in seismic velocity that are associated with the GPS-inferred rapid inflation event. We also determine temporal changes in waveform decorrelation and show a continual increase in decorrelation rate over the time associated with the rapid inflation event. Themagnitude of relative velocity decreases and decorrelation rate increases are comparable to previous studies at Piton de la Fournaise that associate such changes with increased production of volatiles and/ormagmatic intrusion within the magma reservoir and associated opening of fractures and/or fissures. Notably, the largest decrease in relative velocity occurs along the intrastation path passing nearest to the center of the caldera. This observation, along with equal amplitude relative velocity decreases revealed via analysis of intracaldera autocorrelations, suggests that the inflation sourcemay be located approximately within the center of the caldera and represent recharge of shallow magma storage in this location. Importantly, there is a relative absence of seismicity associated with this and previous rapid inflation events at Okmok. Thus, these ANI results are the first seismic evidence of such rapid inflation at the volcano.

Alaska↗

USGS2014

Three-dimensional seismic velocity structure and earthquake relocations at Katmai, Alaska

We invert arrival time data from local earthquakes occurring between September 2004 and May 2009 to determine the three-dimensional (3D) upper crustal seismic structure in the Katmai volcanic region. Waveforms for the study come from the Alaska Volcano Observatory's permanent network of 20 seismic stations in the area (predominantly single-component, short period instruments) plus a densely spaced temporary array of 11 broadband, 3-component stations. The absolute and relative arrival times are used in a double-difference seismic tomography inversion to solve for 3D P- and S-wave velocity models for an area encompassing the main volcanic centers. The relocated hypocenters provide insight into the geometry of seismogenic structures in the area, revealing clustering of events into four distinct zones associated with Martin, Mageik, Trident-Novarupta, and Mount Katmai. The seismic activity extends from about sea level to 2 km depth (all depths referenced to mean sea level) beneath Martin, is concentrated near 2 km depth beneath Mageik, and lies mainly between 2 and 4 km depth below Katmai and Trident-Novarupta. Many new features are apparent within these earthquake clusters. In particular, linear features are visible within all clusters, some associated with swarm activity, including an observation of earthquake migration near Trident in 2008. The final velocity model reveals a possible zone of magma storage beneath Mageik, but there is no clear evidence for magma beneath the Katmai-Novarupta area where the 1912 eruptive activity occurred, suggesting that the storage zone for that eruption may have largely been evacuated, or remnant magma has solidified.

Alaska↗

USGS2012

Earthquake studies reveal the magmatic plumbing system of the Katmai volcanoes

The 1912 eruption of Novarupta was the largest of the 1900s (Fierstein and Hildreth 2001, Hildreth et al. 2003). A century later, fundamental questions remain regarding the source of the magma for that eruption. A previous seismic study of the Katmai area (Jolly et al. 2007) identified a single large area of anomalous structure in the subsurface centered beneath Katmai Pass (Figure 2), but the magma source for the 1912 eruption is thought to have been beneath Mt. Katmai (Hildreth et al. 2003). This mystery was a prime motivation for the research project described here. In summer 2008, scientists and staff from the Alaska Volcano Observatory (AVO) and the University of Wisconsin- Madison installed 11 temporary seismic recording instruments around the Katmai Pass area, complementing the existing AVO seismic network stations (Figure 3). The primary goal of the deployment was to record data from local earthquakes in order to yield an improved model of the three-dimensional structure of the upper crust beneath and surrounding Katmai Pass, using an analysis method known as double-difference seismic tomography (Zhang and Thurber 2003). The method yields a threedimensional image of the velocity of seismic waves in the subsurface, and also produces improved estimates of the locations of the earthquakes beneath the seismic stations. Our main finding is that there is not a single large anomalous zone centered beneath Katmai Pass; rather there are several separate anomalous zones, one each beneath Katmai, Trident-Novarupta, and Martin-Mageik. Furthermore, the earthquakes are tightly clustered beneath the various volcanic centers, and are found to be systematically deeper than previously thought. Linear trends of earthquakes are also revealed, similar to features observed at other volcanoes, possibly outlining previously unidentified fault structures or indicating the path of migrating magma or magmatic fluids and gases.

Alaska↗

USGS2010

Aftershock distribution as a constraint on the geodetic model of coseismic slip for the 2004 Parkfield earthquake

Several studies of the 2004 Parkfield earthquake have linked the spatial distribution of the event’s aftershocks to the mainshock slip distribution on the fault. Using geodetic data, we find a model of coseismic slip for the 2004 Parkfield earthquake with the constraint that the edges of coseismic slip patches align with aftershocks. The constraint is applied by encouraging the curvature of coseismic slip in each model cell to be equal to the negative of the curvature of seismicity density. The large patch of peak slip about 15 km northwest of the 2004 hypocenter found in the curvature-constrained model is in good agreement in location and amplitude with previous geodetic studies and the majority of strong motion studies. The curvature-constrained solution shows slip primarily between aftershock “streaks” with the continuation of moderate levels of slip to the southeast. These observations are in good agreement with strong motion studies, but inconsistent with the majority of published geodetic slip models. Southeast of the 2004 hypocenter, a patch of peak slip observed in strong motion studies is absent from our curvature-constrained model, but the available GPS data do not resolve slip in this region. We conclude that the geodetic slip model constrained by the aftershock distribution fits the geodetic data quite well and that inconsistencies between models derived from seismic and geodetic data can be attributed largely to resolution issues.

California↗