Geology ReportsSearch

Geology topics

Daniel Dzurisin

Publications and source records attributed to Daniel Dzurisin.

At least 19 recordsLinked to original sources

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

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

Volcano deformation: Insights into magmatic systems

Volcano geodesy is the branch of geodetic science that deals with the changing shapes of volcanoes, whether large or small, deep-seated or surficial. Together with seismicity and volcanic gas flux, deformation of the ground surface can be a key indicator of subsurface conditions and processes at volcanoes—information that not only improves scientific understanding of magmatic systems but also is useful for assessing volcano hazards and mitigating their potential consequences. To take full advantage of such information requires detailed characterization of the deformation field in space and time. Currently, no single geodetic technique is capable of providing both the high spatial and temporal resolution required. However, important advances have been made recently by combining information from real-time in situ sensors such as continuous GPS, strainmeters, and tiltmeters with repeated campaign-style GPS, microgravity, interferometric synthetic aperture radar (InSAR), lidar, and photogrammetric observations. Continuous real-time data constrain the timing but not necessarily the spatial extent and pattern of deformation, whereas InSAR provides detailed spatial information but only at intervals of several days to weeks. Repeated microgravity surveys, when combined with independent measurements of surface height, are uniquely sensitive to changes in subsurface mass distribution and therefore can be used to distinguish among processes driven by magma, hydrous fluids, or gas. Simultaneous analysis of multiple geodetic datasets, especially when guided by information from global volcano databases and numerical simulations of volcanic processes and products, can provide a statistical basis for outcome prediction and serve as a guide for additional observations. In the foreseeable future, interactions among magmatic, tectonic, and hydrothermal systems will be monitored and modeled at regional scale in real time, enabling new insights into Earth’s subsurface environment.

Book chapter

Space-based imaging radar studies of U.S. volcanoes

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.

Frontiers in Earth Science

Interferometric synthetic aperture radar study of recent eruptive activity at Shrub mud volcano, Alaska

Shrub mud volcano is one of three large mud volcanoes that comprise the Klawasi Group in the Copper River Basin of southcentral Alaska. Except for minor discharges in the mid-1950s when the group was first described, Shrub was dormant prior to its reactivation in summer 1996. From 1997 to 1999, Shrub vigorously erupted more than 5 x 105 cubic meters of saline mud and carbon dioxide-rich gas at temperatures as high as 54 degrees C. Thereafter, activity waned but continued at least through 2015. We analyzed 192 interferograms derived from 106 synthetic aperture radar (SAR) images acquired by the JERS-1 (L-band), ERS-1/2 (C-band), RADARSAT-1 (C-band), and ALOS PALSAR (L-band) satellites to characterize ground deformation at Shrub before, during, and after its reactivation. Collectively, the interferograms span 1992–2000 and 2006–2011. We fit the observations with two deformation sources: a deflating, steeply dipping, pipe-like body under the summit area and an inflating, shallow-dipping, sill-like body under the southwest flank. Both sources are shallow, with centroids less than 1 km beneath the summit. Prior to reactivation, the flank source inflated ~0.35 x 105 cubic meters per year from July 1992 to May 1996. During eruptive activity, the summit source deflated at higher rates that peaked at ~8.71 x 105 cubic meters per year during May–November 1997 and continued at ~0.95 x 105 cubic meters per year during the 2006–2011 observation window. Cumulative source-volume loss is comparable to the volume of mud erupted. We interpret the summit source as the volcano’s feeder conduit that pressurized prior to the first SAR observation in 1992. Also before 1992, the conduit ruptured to feed a lateral intrusion of mud under the southwest flank, perhaps along a bedding plane in underlying glaciolacustrine deposits. The growing sill caused the southwest flank to inflate while it accommodated the mud supply from depth, which explains why we observed pre-eruptive inflation of the flank but not the summit. The summit began deflating when the conduit ruptured to the surface at the onset of eruptive activity. The flank source did not deflate concurrently because the weight of the thin overburden was insufficient to collapse the sill. There is a suggestion in the modern topography that lateral intrusions under Shrub’s southwest flank are a common feature of activity there.

Alaska

Mount St. Helens retrospective: Lessons learned since 1980 and remaining challenges

Since awakening from a 123-year repose in 1980, Mount St. Helens has provided an opportunity to study changes in crustal magma storage at an active arc volcano—a process of fundamental importance to eruption forecasting and hazards mitigation. There has been considerable progress, but important questions remain unanswered. Was the 1980 eruption triggered by an injection of magma into an upper crustal reservoir? If so, when? How did magma rise into the edifice without producing detectable seismicity deeper than ∼2.5 km or measurable surface deformation beyond the volcano’s north flank? Would precursory activity have been recognized earlier if current monitoring techniques had been available? Despite substantial improvements in monitoring capability, similar questions remain after the dome-forming eruption of 2004–2008. Did additional magma accumulate in the reservoir between the end of the 1980–1986 eruption and the start of the 2004–2008 eruption? If so, when? What is the significance of a relative lull in seismicity and surface deformation for several years prior to the 2004–2008 eruption onset? How did magma reach the surface without producing seismicity deeper than ∼2 km or measurable deformation more than a few hundred meters from the vent? Has the reservoir been replenished since the eruption ended, and is it now primed for the next eruption? What additional precursors, if any, should be expected? This paper addresses these questions, explores possible answers, and identifies unresolved issues in need of additional study. The 1980–1986 and 2004–2008 eruptions could have resulted from second boiling during crystallization of magma long-resident in an upper crustal reservoir, rather than from injection of fresh magma from below. If reservoir pressurization and magma ascent were slow enough, resulting strain might have been accommodated by viscoelastic deformation, without appreciable seismicity or surface deformation, until rising magma entered a brittle regime within 2–2.5 km of the surface. Given the remarkably gas-poor nature of the 2004–2008 dome lava, future eruptive activity might require a relatively long period of quiescence and reservoir pressurization or a large injection of fresh magma—an event that arguably has not occurred since the Kalama eruptive period (C.E. 1479–1720).

Washington

Mass addition at Mount St. Helens, Washington, inferred from repeated gravity surveys

Measurements of subtle changes in the Earth’s gravityfield can provideinformation on the addition/loss of mass (e.g., magma or aqueousfluids) beneath a volcano. In this study, wemeasured gravity at Mount St. Helens from 2010 to 2016 to investigate possible mass changes followingthe 2004–2008 dome-forming eruption. The raw gravity measurements were corrected for changes in themass and shape of Crater Glacier and for deformation of the volcanic edifice to obtain residual gravity values.We found positive residual gravity changes that we interpret as evidence for partial recharge of the magmareservoir that fed the 2004–2008 eruption and/or for accumulation of groundwater in one or more shallowaquifers. Most of the gravity signal can be explained by groundwater accumulation; magma recharge is notrequired. However, there is additional evidence for at least partial recharge from seismic, deformation, andgeochemical observations

Washington

Magma supply to Kīlauea Volcano, Hawai‘i, from inception to now: Historical perspective, current state of knowledge, and future challenges

Meticulous field observations are a common underpinning of two landmark studies conducted by Don Swanson dealing with the rate at which magma is supplied to Kīlauea Volcano, Hawai‘i. The first combined effusion rate and ground deformation observations to show that the supply rate to Kīlauea was constant at ~0.11 km 3 /yr during three sustained eruptions from 1952 to 1971, a quiescent period at neighboring Mauna Loa volcano. This rate was also interpreted as the steady supply rate from the mantle to both volcanoes combined throughout historical time. The second breakthrough involved field evidence that activity at Kīlauea alternates between dominantly effusive and explosive styles over time scales of several centuries, and that the magma supply rate during explosive periods is only 1%-2% of the rate during effusive periods. For the historical period, several later studies concluded that the supply rate to Kīlauea has varied by as much as an order of magnitude, contrary to Swanson’s ­earlier suggestion. All such estimates are fraught with uncertainty, given the poorly known amount of magma stored within the volcano’s rift zones as a function of time—an enduring problem and active research topic. Nonetheless, Swanson’s original work remains an important touchstone that spurred many subsequent investigations and refinements. For example, there is strong evidence that Kīlauea experienced a surge in magma supply during 2003–2007 that exceeded the historical average by as much as a factor of two, and that the surge was followed by a comparable lull before the supply rate returned to “normal” by 2016. There is also evidence for supply-rate variations of similar magnitude during the latter part of the twentieth century and possibly earlier, subject to the aforementioned uncertainty in rift-zone storage. The extent to which variations in the magma supply to Kīlauea can be attributed to partitioning between Kīlauea and Mauna Loa, a long-debated topic, remains uncertain. Since Kilauea’s inception, the net magma supply to the volcano (and also to Lō‘ihi Seamount, since it began growing) has increased, while Mauna Loa’s growth rate has slowed, suggesting that the volcanoes compete for the same magma supply. However, geochemical differences between lavas erupted at Kīlauea and Mauna Loa indicate that they do not share a homogeneous mantle source or common lithospheric magma plumbing system. Both ideas might be correct; i.e., Kīlauea and Mauna Loa magmas may be sourced in differing portions of the same melt accumulation zone and ascend through different crustal pathways, but those pathways interact through stress or pressure changes that modulate the supply to each volcano. Currently, magma supply-rate estimates are facilitated by comprehensive imaging of surface deformation and topographic change coupled with measurements of gas emissions. Physics-based models are being developed within a probabilistic framework to provide rigorous estimates of model parameters, including magma supply rate, and their uncertainties. Further refinement will require intensive multiparameter observations of the entire magmatic system—from source to surface and above, and from the volcanoes’ summits to their submerged lower flanks—in order to account fully for a complex magma budget.

Hawai‘i

Volcano geodesy in the Cascade arc, USA

Experience during historical time throughout the Cascade arc and the lack of deep-seated deformation prior to the two most recent eruptions of Mount St. Helens might lead one to infer that Cascade volcanoes are generally quiescent and, specifically, show no signs of geodetic change until they are about to erupt. Several decades of geodetic data, however, tell a different story. Ground- and space-based deformation studies have identified surface displacements at five of the 13 major Cascade arc volcanoes that lie in the USA (Mount Baker, Mount St. Helens, South Sister, Medicine Lake, and Lassen volcanic center). No deformation has been detected at five volcanoes (Mount Rainier, Mount Hood, Newberry Volcano, Crater Lake, and Mount Shasta), and there are not sufficient data at the remaining three (Glacier Peak, Mount Adams, and Mount Jefferson) for a rigorous assessment. In addition, gravity change has been measured at two of the three locations where surveys have been repeated (Mount St. Helens and Mount Baker show changes, while South Sister does not). Broad deformation patterns associated with heavily forested and ice-clad Cascade volcanoes are generally characterized by low displacement rates, in the range of millimeters to a few centimeters per year, and are overprinted by larger tectonic motions of several centimeters per year. Continuous GPS is therefore the best means of tracking temporal changes in deformation of Cascade volcanoes and also for characterizing tectonic signals so that they may be distinguished from volcanic sources. Better spatial resolution of volcano deformation can be obtained through the use of campaign GPS, semipermanent GPS, and interferometric synthetic aperture radar observations, which leverage the accumulation of displacements over time to improve signal to noise. Deformation source mechanisms in the Cascades are diverse and include magma accumulation and withdrawal, post-emplacement cooling of recent volcanic deposits, magmatic-tectonic interactions, and loss of volatiles plus densification of magma. The Cascade Range thus offers an outstanding opportunity for investigating a wide range of volcanic processes. Indeed, there may be areas of geodetic change that have yet to be discovered, and there is good potential for addressing a number of important questions about how arc volcanoes work before, during, and after eruptions by continuing geodetic research in the Cascade Range.

Cascade Arc

Semipermanent GPS (SPGPS) as a volcano monitoring tool: Rationale, method, and applications

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 ~ 3 month 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 ~ 1 week SGPS measurements. After 5 years of SPGPS observations at Three Sisters, the precision of velocity determinations is estimated to be 0.5 mm/yr in longitude, 0.6 mm/yr in latitude, and 0.8 mm/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.

Journal of Volcanology and Geothermal Research

The 2004–2008 dome-building eruption at Mount St. Helens, Washington: Epilogue

The 2004–2008 dome-building eruption at Mount St. Helens ended during winter 2007–2008 at a time when field observations were hampered by persistent bad weather. As a result, recognizing the end of the eruption was challenging—but important for scientists trying to understand how and why long-lived eruptions end and for public officials and land managers responsible for hazards mitigation and access restrictions. In hindsight, the end of the eruption was presaged by a slight increase in seismicity in December 2007 that culminated on January 12–13, 2008, with a burst of more than 500 events, most of which occurred in association with several tremor-like signals and a spasmodic burst of long-period earthquakes. At about the same time, a series of regular, localized, small-amplitude tilt events—thousands of which had been recorded during earlier phases of the eruption—came to an end. Thereafter, seismicity declined to 10–20 events per day until January 27–28, when a spasmodic burst of about 50 volcano-tectonic earthquakes occurred over a span of 3 h. This was followed by a brief return of repetitive “drumbeat” earthquakes that characterized much of the eruption. By January 31, however, seismicity had declined to 1–2 earthquakes per day, a rate similar to pre-eruption levels. We attribute the tilt and seismic observations to convulsive stagnation of a semisolid magma plug in the upper part of the conduit. The upward movement of the plug ceased when the excess driving pressure, which had gradually decreased throughout the eruption as a result of reservoir deflation and increasing overburden from the growing dome, was overcome by increasing friction as a result of cooling and crystallization of the plug.

Washington

InSAR imaging of aleutian volcanoes: Monitoring a volcanic arc from space

Interferometric synthetic aperture radar (InSAR) is a relatively new remote sensing tool that is capable of measuring ground-surface deformation with centimeter-to-subcentimeter precision at a spatial resolution of tens of meters over an area of hundreds to thousands of square kilometers. With its global coverage and all-weather imaging capability, InSAR has become an increasingly important technique for studying volcanoes in remote regions such as the Aleutian Islands. The spatial distribution of surface deformation data derived from InSAR images enables the construction of detailed mechanical models to enhance the study of magmatic processes.

Book

Pre-eruption deformation caused by dike intrusion beneath Kizimen volcano, Kamchatka, Russia, observed by InSAR

Interferometric synthetic aperture radar (InSAR) images reveal a pre-eruption deformation signal at Kizimen volcano, Kamchatka, Russia, where an ongoing eruption began in mid-November, 2010. The previous eruption of this basaltic andesite-to-dacite stratovolcano occurred in 1927–1928. InSAR images from both ascending and descending orbital passes of Envisat and ALOS PALSAR satellites show as much as 6 cm of line-of-sight shortening from September 2008 to September 2010 in a broad area centered at Kizimen. About 20 cm of opening of a nearly vertical dike provides an adequate fit to the surface deformation pattern. The model dike is approximately 14 km long, 10 km high, centered 13 km beneath Kizimen, and strikes NE–SW. Time-series analysis of multi-temporal interferograms indicates that (1) intrusion started sometime between late 2008 and July 2009, (2) continued at a nearly constant rate, and (3) resulted in a volume expansion of 3.2 × 10 7 m 3 by September 2010, i.e., about two months before the onset of the 2010 eruption. Earthquakes located above the tip of the dike accompanied the intrusion. Eventually, magma pressure in the dike exceeded the confining strength of the host rock, triggering the 2010 eruption. Our results provide insight into the intrusion process that preceded an explosive eruption at a Pacific Rim stratovolcano following nearly a century of quiescence, and therefore have implications for monitoring and hazards assessment at similar volcanoes elsewhere.

Kamchatka, Kizimen Volcano

Dynamic deformation of Seguam Island, Alaska, 1992--2008, from multi-interferogram InSAR processing

We generated a time-series of ERS-1/2 and ENVISAT interferometric synthetic aperture radar (InSAR) images to study ground surface deformation at Seguam Island from 1992 to 2008. We used the small baseline subset (SBAS) technique to reduce artifacts associated with baseline uncertainties and atmospheric delay anomalies, and processed images from two adjacent tracks to validate our results. Seguam Island comprises the remnants of two late Quaternary calderas, one in the western caldera of the island and one in the eastern part of the island. The western caldera subsided at a constant rate of ~ 1.6 cm/yr throughout the study period, while the eastern caldera experienced alternating periods of subsidence and uplift: ~ 5 cm/year uplift during January 1993–October 1993 (stage 1), ~ 1.6 cm/year subsidence during October 1993–November 1998 (stage 2), ~ 2.0 cm/year uplift during November 1998–September 2000 (stage 3), ~ 1.4 cm/year subsidence during September 2000–November 2005 (stage 4), and ~ 0.8 cm/year uplift during November 2005– July 2007 (stage 5). Source modeling indicates a deflationary source less than 2 km below sea level (BSL) beneath the western caldera and two sources beneath the eastern caldera: an inflationary source 2.5–6.0 km BSL and a deflationary source less than 2 km BSL. We suggest that uplift of the eastern caldera is driven by episodic intrusions of basaltic magma into a poroelastic reservoir 2.5–6.0 km BSL beneath the caldera. Cooling and degassing of the reservoir between intrusions results in steady subsidence of the overlying surface. Although we found no evidence of magma intrusion beneath the western caldera during the study period, it is the site (Pyre Peak) of all historical eruptions on the island and therefore cooling and degassing of intrusions presumably contributes to subsidence there as well. Another likely subsidence mechanism in the western caldera is thermoelastic contraction of lava flows emplaced near Pyre Peak during several historical eruptions, most recently in 1977 and 1992–93.

Alaska

Mount St. Helens, 1980 to now— What’s going on?

Mount St. Helens seized the world’s attention in 1980 when the largest historical landslide on Earth and a powerful explosive eruption reshaped the volcano, created its distinctive crater, and dramatically modified the surrounding landscape. An enormous lava dome grew episodically in the crater until 1986, when the volcano became relatively quiet. A new glacier grew in the crater, wrapping around and partly burying the lava dome. From 1987 to 2003, sporadic earthquake swarms and small steam explosions indicated that magma (molten rock) was being replenished deep underground. In 2004, steam-and-ash explosions heralded the start of another eruption. A quieter phase of continuous lava extrusion followed and lasted until 2008, building a new dome and doubling the volume of lava on the crater floor. Scientists with the U.S. Geological Survey and University of Washington’s Pacific Northwest Seismograph Network maintain constant watch for signs of renewed activity at Mount St. Helens and other Cascade volcanoes. Now is an ideal time for both actual and virtual visitors to Mount St. Helens to learn more about dramatic changes taking place on and beneath this active volcano.

Washington

Crustal subsidence and extension and Medicine Lake volcano, northern California

The pattern of historical ground deformation, seismicity, and crustal structure near Medicine Lake volcano illustrates a close relation between magmatism and tectonism near the margin of the Cascade volcanic chain and the Basin and Range tectonic province. Between leveling surveys in 1954 and 1989 the summit of Medicine Lake volcano subsided 389±43 mm with respect to a reference bench mark 40 km to the southwest (average rate = 11.1±1.2 mm/yr). A smaller survey across the summit caldera in 1988 suggests that the subsidence rate was 15–28 mm/yr during 1988–1989. Swarms of shallow earthquakes ( M ≤ 4.6) occurred in the region during August 1978, January–February 1981, and September 1988. Except for the 1988 swarm, which occurred beneath Medicine Lake caldera, most historical earthquakes were located at least 25 km from the summit. The spatial relation between subsidence and seismicity indicates (1) radially symmetric downwarping of the volcano's summit and flanks centered near the caldera and (2) downfaulting of the entire edifice along regional faults located 25–30 km from the summit. We propose that contemporary subsidence, seismicity, and faulting are caused by (1) loading of the crust by more than 600 km 3 of erupted products plus a large volume of mafic intrusives; (2) east‐west extension in the western Basin and Range province; and, to a lesser extent, (3) crystallization or withdrawal of magma beneath the volcano. Thermal weakening of the subvolcanic crust by mafic intrusions facilitates subsidence and influences the distribution of earthquakes. Subsidence occurs mainly by aseismic creep within 25 km of the summit, where the crust has been heated and weakened by intrusions, and by normal faulting during episodic earthquake swarms in surrounding, cooler terrain.

California

Variations in magma supply rate at Kilauea Volcano, Hawaii

When an eruption of Kilauea lasts more than 4 months, so that a well-defined conduit has time to develop, magma moves freely through the volcano from a deep source to the eruptive site at a constant rate of 0.09 km 3 /yr. At other times, the magma supply rate to Kilauea, estimated from geodetic measurements of surface displacements, may be different. For example, after a large withdrawal of magma from the summit reservoir, such as during a rift zone eruption, the magma supply rate is high initially but then lessens and exponentially decays as the reservoir refills. Different episodes of refilling may have different average rates of magma supply. During four year-long episodes in the 1960s, the annual rate of refilling varied from 0.02 to 0.18 km 3 /yr, bracketing the sustained eruptive rate of 0.09 km 3 /yr. For decade-long or longer periods, our estimate of magma supply rate is based on long-term changes in eruptive rate. We use eruptive rate because after a few dozen eruptions the volume of magma that passes through the summit reservoir is much larger than the net change of volume of magma stored within Kilauea. The low eruptive rate of 0.009 km 3 /yr between 1840 and 1950, compared to an average eruptive rate of 0.05 km 3 /yr since 1950, suggests that the magma supply rate was lower between 1840 and 1950 than it has been since 1950. An obvious difference in activity before and since 1950 was the frequency of rift zone eruptions: eight rift zone eruptions occurred between 1840 and 1950, but more than 20 rift zone eruptions have occurred since 1950. The frequency of rift zone eruptions influences magma supply rate by suddenly lowering pressure of the summit magma reservoir, which feeds magma to rift zone eruptions. A temporary drop of reservoir pressure means a larger-than-normal pressure difference between the reservoir and a deeper source, so magma is forced to move upward into Kilauea at a faster rate.

Journal of Geophysical Research B: Solid Earth