Geology topics

Michael Lisowski

Publications and source records attributed to Michael Lisowski.

At least 37 records · Page 2Linked to original sources

USGS2012

The velocity field along the San Andreas Fault in central and southern California

The velocity field within a 100‐km‐broad zone centered on the San Andreas fault between the Mexican border and San Francisco Bay has been inferred from repeated surveys of trilateration networks in the 1973–1989 interval. The velocity field has the appearance of a shear flow that remains parallel to the local strike of the fault even through such major deflections as the big bend of the San Andreas fault in the Transverse Ranges of southern California. Across‐strike profiles of the fault‐parallel component of velocity exhibit the expected sigmoidal shape, whereas across‐strike profiles of the fault‐normal component of velocity are flat and featureless. No significant convergence upon the fault is observed even along the big bend sector of the fault. Simple dislocation models can explain most of the features of the observed velocity field, but those explanations are not unique. About 35 mm/yr of relative plate motion is accounted for within the span of the trilateration networks. Geologic studies indicate that the secular slip rate on the San Andreas fault is about 35 mm/yr. The agreement between these two estimates implies that most of the strain accumulation is elastic and will be recovered in subsequent earthquakes. The relative motion observed across the San Andreas fault (35 mm/yr) plus that observed across the Eastern California shear zone (8 mm/yr) accounts for most (43 mm/yr) of the observed North America‐Pacific relative plate motion (47 mm/yr).

California↗

USGS2012

Strain accumulation in western Washington

The Juan de Fuca plate is subducted beneath the North American plate off the coast of Washington at a rate of about 40 mm/yr N68°E. The average principal strain rates (extension reckoned positive) measured in northwestern Washington are as follows: Olympic peninsula 25 km south of Port Angeles from 1982 through 1990, and and near Seattle from 1972 through 1985, and . Both strain measurements are consistent with uniaxial contraction in the direction of plate convergence. Uplift rates inferred from tide gage recordings are about 4 mm/yr on the Pacific coast and near 0 mm/yr farther inland near Seattle. These deformation rates are consistent with a model of the Cascadia subduction zone in which the plate interface beneath the continental slope and outer continental shelf is locked but free to slip farther landward. The limited downdip extent of the locked segment of the plate interface is consistent with a shallow depth (∼20 km) of the isotherm (∼450°C) that defines the brittle‐ductile transition. Small thrust events diagnostic of seismic subduction should then occur only offshore and at shallow depths. The principal strain rates measured from 1972 through 1983 in the back arc region near Richland, Washington, are and .

Washington↗

USGS2012

Strain accumulation along the Denali Fault at the Nenana River and Delta River Crossings, Alaska

Surveys of trilateration networks across the Denali fault at the Nenana River in 1982, 1984, and 1988 and at the Delta River in 1975, 1979, 1982, and 1984 indicate a minor (0.10±0.04 μstrain/yr) northeastward uniaxial extension. The component of right‐lateral shear‐strain accumulation across the fault is not significant at the two‐standard‐deviation level. At the Delta River network the strain accumulation rate decreases rapidly with distance from the fault, but evidence for a similar decrease with distance from the fault is lacking at the Nenana River network. The strain accumulation rates inferred from trilateration are consistent with the very long baseline interferometry (VLBI) measurement reported by Ma et al. (1990) and support their contention that significant right‐lateral shear is not accumulating along the Denali fault at the present time. Savage et al. (1981) had earlier concluded erroneously that preliminary geodetic measurements at the Delta River network demonstrated right‐lateral shear strain accumulation. The absence of significant right‐lateral deformation across the Denali fault in the 1975–1988 interval is in marked contrast with the abundant geomorphic evidence for Holocene right‐lateral secular slip at the rate of 10–20 mm/yr on the Denali fault in this sector.

Alaska↗

USGS2012

Strain accumulation across the Wasatch Fault near Ogden, Utah

Deformation of a 70 by 40‐km trilateration network spanning the north trending Wasatch fault near Ogden, Utah, has been monitored from 1972 through 1990. All but nine of the 200 measurements are consistent with deformation that is linear in time. We presume that those nine observations are contaminated by some blunder in making the measurements and that deformation is linear in time. The strain rate over the 1972–1990 interval across a 40‐km‐wide zone west of the Wasatch fault (hanging wall block) is adequately approximated by a 0.04 ± 0.01 μstrain/yr N85°E±5° uniaxial extension. Adequate coverage of the footwall block is available only over the 1981–1990 interval, and the strain rate across a 30‐km‐wide zone there is approximated by a 0.07 ± 0.03 μstrain/yr N20°E±12° uniaxial extension. The observed extension of the hanging wall block in the direction normal to fault strike suggests that the Wasatch fault may be listric; a planar fault would produce a zone of contraction normal to the fault at the surface of the hanging wall block. The component of extension parallel to fault strike observed on the footwall block is unexplained.

Utah↗

USGS2012

Deformation from 1973 to 1987 in the epicentral area of the 1959 Hebgen Lake, Montana, earthquake (Ms = 7.5)

A 40‐km aperture trilateration network centered on the 1959 Hebgen Lake earthquake epicenter has been surveyed in 1973, 1974, 1976, 1978, 1981, 1984, and 1987. The deformation inferred from those surveys is described roughly by a uniaxial, 0.266 ± 0.014 μstrain/yr, N15°E ± 1°extension that is uniform in both time and space. That extension is orthogonal to the strike (N78°W ± 5°) of the 1959 rupture plane inferred from focal mechanism solutions. The absence of strain accumulation in the N75°W direction suggests that the source of deformation must be greatly elongated in that direction. Looked at in greater detail, the deformation is found to be concentrated in the northern part of the network. The zone of concentrated deformation coincides with a part of a well‐defined trend in seismicity that extends at least 100 km N75°W from the Sour Creek resurgent dome in the Yellowstone caldera. Extension perpendicular to that trend is indicated by focal mechanism solutions. Thus, the zone of rifting identified within the Hebgen Lake network apparently extends east southeast to the Yellowstone caldera. About 8 mm/yr extension is observed across a 20‐km width within the zone.

Montana, Wyoming↗

USGS2012

Deformation from 1973 through 1991 in the epicentral area of the 1992 Landers, California, Earthquake (Ms = 7.5)

Deformation of a 50 × 60 km trilateration network that spans the epicenter of the 1992 Landers earthquake( M s = 7.5) was measured by seven surveys over the 19 years preceding the earthquake. Three moderate earthquakes (1979 Homestead Valley, M s = 5.6; 1986 North Palm Springs, M s = 6.0; and 1992 Joshua Tree, M s = 6.1) occurred within the network during those 19 years. Here we use geodetic and seismic data to construct a dislocation model for each of the three moderate earthquakes. Coseismic changes due to these three moderate earthquakes as predicted by the dislocation models are then removed from the trilateration data. The residual geodetic changes appear to be uniform in time. We take those changes to represent secular strain accumulation. This strain accumulation clearly shows right‐lateral shear across the San Andreas fault but also suggests a northwest‐southeast extension northeast of the eastern end of the “Big Bend” of the San Andreas fault. The Landers earthquake ruptured along a north‐northwest trend across this region of northwest‐southeast extension.

California↗

USGS2012

Postseismic deformation following the 1989 (M = 7.1): Loma Prieta, California, earthquake

Postseismic deformation along a 90‐km profile bisecting the projected surface trace of the coseismic rupture of the 1989 Loma Prieta earthquake has been monitored by frequent GPS surveys for 3.3 years following the earthquake. In addition to the expected deformation associated with secular strain accumulation on the San Andreas and Calaveras faults, deformation associated with postseismic readjustment has been detected. Most of that deformation can be attributed to 1.5 m right‐lateral and 0.9 m reverse postseismic slip on a 5‐km‐wide downdip extension of the Loma Prieta rupture. In addition, there seems to be a 0.1 m postseismic collapse of the Loma Prieta rupture zone in the direction perpendicular to the plane of the rupture. The fault‐normal (N48°E) surface displacements plotted as a function of time exhibit a curvature suggesting a relaxation time of about 1.4 years. Similar plots of the fault‐parallel (N42°W) displacement components do not exhibit significant curvature. Presumably, the deformation shown in those plots is dominated by secular strain accumulation along the San Andreas and Calaveras faults rather than postseismic relaxation.

California↗

USGS2012

Strain accumulation along the Laguna Salada Fault, Baja California, Mexico

Strain accumulation observed over the 1978–1991 interval in a 30×100 km aperture trilateration network spanning the Laguna Salada fault is described by the principal strain rates 0.101±0.012 strain/yr N80°E±2° and −0.021±0.012 strain/yr N10°W±2°, extension reckoned positive. These strain accumulation rates have been corrected to remove coseismic effects of the nearby 1979 Imperial Valley ( = 6.5), 1980 Victoria (Baja California) ( = 6.4), 1987 Superstition Hills ( = 6.5), and 1987 Elmore Ranch ( = 5.9) earthquakes. The observed strain rates indicate extension at a rate of about 0.08 strain/yr perpendicular to the trend (N35°W) of the Salton trough as well as a right‐lateral tensor shear strain rate 0.05 strain/yr across it. The extension perpendicular to the trough is observed neither farther north near the Salton Sea nor farther south across the Gulf of California. However, Holocene slip on the Laguna Salada fault, about equal parts right‐lateral and normal slip, is consistent with the observed strain accumulation. A simple dislocation model intended to explain the observed strain accumulation as a product of slip at depth on the Laguna Salada fault would require that the fault be listric.

California↗

USGS2012

Geodetic monitoring of the southern San Andreas Fault, California, 1980-1991

Five geodetic arrays (10 to 40 km aperture) located along the San Andreas fault have been surveyed frequently (several times in most years) over the 1980–1991 interval to detect possible fluctuations in the deformation rate. In each survey of an array the distances between the same four to seven pairs of geodetic monuments were measured. The distances measured (with corresponding standard deviation) range from 8.4 km (3.4 mm) to 38 km (8.2 mm). The data are displayed as plots of measured distance as a function of time. Linear fits in such plots furnish a satisfactory representation of the data. That is, the scatter of the data about the linear fits is within the range expected for the estimated standard deviations in measurement. However, there are coherent low‐amplitude (within the observational error) fluctuations apparent in some of the measured distances. Those fluctuations need not be tectonic effects but rather may be either random patterns or artifacts introduced by systematic drift in instrument calibration or wander of the geodetic monuments. We conclude that the measurements are consistent with steady deformation of the arrays over the 1980–1991 interval, and we find no convincing evidence in the data for fluctuations in the rate of deformation.

California↗

USGS2012

Strain accumulation across the central Nevada seismic zone, 1973–1994

Five trilateration networks extending for 280 km along the central Nevada seismic zone (1915 Pleasant Valley, M = 7.3; 1954 Dixie Valley, M = 6.8; 1954 Stillwater, M = 6.8; 1954 Rainbow Mountain, M = 6.6; 1954 Fairview Peak, M = 7.1; and 1932 Cedar Mountain, M = 7.2) have been surveyed 6 times since 1973 to determine deformation along the zone. Within the precision of measurement the deformation appears uniform along the zone and is described by the principal strain rates 0.036±0.008 μstrain/yr N60°W±3° and −0.031±0.008 μstrain/yr N30°E±3°, extension reckoned positive. The observed strain rates are consistent with simple, right‐lateral, tensor shear at the rate of 0.033 μstrain/yr across a shear zone striking N15°W. This central Nevada shear zone appears to be the northward continuation of the eastern California shear zone. The orientation of the strike‐slip and normal‐slip ruptures within the central Nevada seismic zone are consistent with principal stress axes parallel to the measured principal strain rate axes. Space‐based geodetic measurements (very long baseline interferometry) indicate that the relative motion accommodated across the Basin and Range province west of Ely, Nevada, is about 9.1±1.5 mm/yr N16°W±8° (Dixon et al., 1995.) Notice that the right‐lateral shear zone postulated to explain deformation in the central Nevada seismic zone is properly oriented to accommodate that relative motion. However, a 135‐km effective width of the shear zone would be required to accommodate all of the 9.1 mm/yr relative motion at the strain rates observed in the Nevada seismic zone; only about 3 mm/yr of that relative motion is accommodated within the span of the trilateration networks.

California, Nevada↗

USGS2012

Deformation across the Salton Trough, California, 1973-1977

A trilateration network extending across the San Andreas, San Jacinto, and Elsinore faults in the vicinity of the Salton Sea, California, has been surveyed to very high precision several times in the 5‐year interval 1973–1977. The average strain across the entire network is essentially a uniaxial north‐south contraction at the rate of about 0.3 μstrain/a. There is no substantial strain perpendicular to the Salton trough, indicating no tendency to either open or close that rift. The observed uniaxial north‐south contraction differs from a pure shear parallel to the major faults by a uniform dilatation of about −0.25 μstrain/a that is only partly explained. The shear strain across the network reaches a maximum near the San Jacinto fault and could be explained by right lateral slip at depth on that fault. The displacement pattern also suggests right lateral slip at depth on the San Jacinto fault with total right lateral relative movement near 50 mm/a across the 120‐km breadth of the network. The deformation appears to be uniform in time in the 1973–1977 interval. A dislocation model of the deformation suggests about 50±15 mm/a of relative right lateral slip at depth distributed between the San Andreas, San Jacinto, and Elsinore faults and demonstrates that a significant average dilatation can be generated by pure strike slip on several subparallel faults.

California↗

USGS2012

Geodetic measurement of crustal deformation on the San Andreas, Hayward, and Calaveras faults near San Francisco, California

Analysis of a geodetic network of 115 lines crossing the San Andreas, Hayward, and Calaveras faults in the vicinity of San Francisco Bay and measured repeatedly between 1970 and 1980 has revealed details about the accommodation of relative plate motion in this area. The most striking result is that the deformation is not uniformly distributed across the area. In the east bay, along the Hayward and Calaveras faults, all motion appears to take place as slip directly on the fault, with no accumulation of strain in the adjacent crust. On both the Calaveras and the Hayward faults the rate obtained for the 1970 to 1980 period agrees with geologic rates spanning a few million years and with creep rates spanning a few decades. The Hayward fault slip rate is 7 ± 1 mm/yr. The Calaveras fault slip rate is 7 ± 1 mm/yr, with perhaps half of this slip distributed across a zone a few kilometers wide, probably as inelastic deformation of weak near‐surface material. The absence of strain accumulation in the east bay is surprising since the Hayward and Calaveras faults have been the site of large earthquakes in the past. A block located east of the Calaveras fault and south of the Las Positas fault has been rotating clockwise at a rate of 0.3 ± 0.1 μrad/yr with very little internal deformation. Along the San Francisco peninsula no detectable slip occurs (less than 1.5 mm/yr) at the surface, but appreciable strain is accumulating. Near fault shear strain rates are 0.6 ± 0.1 μstrain/yr (engineering) with direction N47°W ± 9. The slip rate near the San Andreas fault is 12.2 ± 3.9 mm/yr distributed across a broad zone. The relative motion across the whole region during the period 1970–1980 is 32.1 ± 7.4 mm/yr.

California↗

USGS2012

Deformation in the White Mountain seismic gap, California-Nevada, 1972-1982

A 100×40 km trilateration network extending from Bishop, California, to near Hawthorne, Nevada, crosses the east end of the Long Valley caldera, site of renewed magma inflation in the 1979–1980 interval, and spans most of the White Mountain seismic gap. The network was surveyed in 1972, 1973, 1976, 1979, 1980, and 1982. The 1980 survey may be contaminated by a scale error. In addition, leveling surveys across the caldera have been run in 1932, 1957, 1975, 1980, 1982, and 1983. Interpretation of the deformation is complicated by the occurrence of the May 1980 Mammoth Lakes earthquake sequence (four earthquakes M L ≥6) at the south edge of the caldera as well as other moderate earthquakes within the White Mountain seismic gap. The vertical deformation is largely accounted for by 0.10‐ to 0.15‐km 3 expansion of a spherical magma chamber 8–10 km beneath the resurgent dome within the Long Valley caldera sometime between July 1979 and September 1980 with an additional expansion of perhaps 0.05 km 3 between September 1980 and July 1982. Some additional sources of deformation within the aftershock zone of the Mammoth Lakes earthquakes seem to be required to explain the horizontal deformation. We show that right‐lateral slip on vertical faults extending WNW from each of the three largest earthquakes in the Mammoth Lakes sequence provides the required additional deformation, but this solution is by no means unique. There are simply too few data to define the rather complex deformation that apparently occurred within the aftershock zone. There is little doubt, however, that inflation of a magma chamber beneath the resurgent dome within the Long Valley caldera was involved in the deformation.

California↗

USGS2010

Cyclic ground tilt associated with the 2004–2008 eruption of Mount St. Helens

The long‐term behavior of the 2004–2008 effusive eruption of Mount St. Helens was characterized by a gradual decline in the rates of seismicity, dome growth, and broad‐scale ground deformation, but shallow near‐periodic “drumbeat” earthquakes over timescales of minutes indicated episodic short‐term behavior. In part to better characterize this behavior and any associated ground deformation, a network of tiltmeters was installed and recorded thousands of cyclic tilt events within the crater. The duration of these events, from minutes to hours, was too long for them to be directly caused by the drumbeat seismicity. Tilt events were characterized by a recoverable, asymmetric pattern of rapid tilt away from the vent followed by a more gradual reversal, were highly correlated between different tiltmeters in the crater, sometimes occurred in association with volcanic seismicity, and ceased at the end of the eruption. Tilt vectors converged on a point just south of the center of the preexisting 1980s lava dome, and the absence of detectable tilt outside the crater suggests a shallow source (<1 km). We examine several models, including cycles of conduit pressurization and plug slip or gas loss and stick‐slip behavior on the interface between the ascending plug and the 1980s lava dome. The small number of stations within the crater prevents a unique determination of source type or geometry, but results are consistent with a mechanism involving extrusion of the semisolid dacite plug and/or cycles of conduit pressurization.

Washington↗

USGS2010

Geodetic constraints on magma movement and withdrawal during the 2006 eruption of Augustine Volcano: Chapter 17 in The 2006 eruption of Augustine Volcano, Alaska

For the first time in the United States, a modern geodetic network of continuously recording Global Positioning System (GPS) receivers has measured a complete eruption cycle at a stratovolcano, Augustine Volcano in Alaska, from the earliest precursory unrest through the return to background quiescence. The on-island network consisted of five continuously recording, telemetered GPS stations, four continuously recording, nontelemetered stations, and about 10 campaign bench marks. The continuous network recorded several distinct and conspicuous signals over the course of the unrest and eruption, starting with a months-long precursory inflation centered beneath the volcano at around sea level. Nearly coincident with the highest volumetric eruption rates, this inflation gave way to a more deep seated deflation that we interpret as a major withdrawal (approx. 25 million m3 of compressed magma) from a nearly cylindrical magma reservoir centered about 5 km below sea level. Detailed analysis of the geodetic time series reveals additional nuance, including the probable upward propagation of a small dike into the edifice in the 60 days or so before the onset of large-scale explosive activity. Comparisons of the geodetic data and their resulting interpretations with other data, such as earthquake hypocenters and petrologically inferred magma-pressure histories, reveal a kinematic, if not mechanical, account of the 2006 eruption that details the shape and location of the magma source region, the means and velocity of magma transport, and the establishment of a short lived volume- (or pressure-) buffering capability held within the magma reservoir. The cumulative deformation over the course of the eruption shows a large signal close in and high on the edifice that decays rapidly with distance. This pattern indicates a small permanent increase in the edifice volume (beyond that added by the surficial lava dome) but also shows that without close-in (<2.5 km from the summit) stations, the eruption might have been invisible to campaign GPS stations alone.

Alaska↗

USGS2009

Continuing inflation at Three Sisters volcanic center, central Oregon Cascade Range, USA, from GPS, leveling, and InSAR observations

Uplift of a broad area centered ~6 km west of the summit of South Sister volcano started in September 1997 (onset estimated from model discussed in this paper) and was continuing when surveyed in August 2006. Surface displacements were measured whenever possible since August 1992 with satellite radar interferometry (InSAR), annually since August 2001 with GPS and leveling surveys, and with continuous GPS since May 2001. The average maximum displacement rate from InSAR decreased from 3–5 cm/yr during 1998–2001 to ~1.4 cm/yr during 2004–2006. The other datasets show a similar pattern, i.e., surface uplift and extension rates decreased over time but deformation continued through August 2006. Our best-fit model to the deformation data is a vertical, prolate, spheroidal point-pressure source located 4.9–5.4 km below the surface. The source inflation rate decreased exponentially during 2001–2006 with a 1/ e decay time of 5.3 ± 1.1 years. The net increase in source volume from September 1997 to August 2006 was 36.5–41.9 x 10 6 m 3 . A swarm of ~300 small ( M max = 1.9) earthquakes occurred beneath the deforming area in March 2004; no other unusual seismicity has been noted. Similar deformation episodes in the past probably would have gone unnoticed if, as we suspect, most are small intrusions that do not culminate in eruptions.

Oregon↗

USGS2008

Constraints and conundrums resulting from ground-deformation measurements made during the 2004-2005 dome-building eruption of Mount St. Helens, Washington

A prolonged period of dome growth at Mount St. Helens starting in September-October 2004 provides an opportunity to study how the volcano deforms before, during, and after an eruption by using modern instruments and techniques, such as global positioning system (GPS) receivers and interferometric synthetic aperture radar (InSAR), together with more traditional ones, including tiltmeters, triangulation, photogrammetry, and time-lapse photography. No precursory ground deformation was detected by campaign GPS measurements made in 2000 and 2003, nor by a continuous GPS station (JRO1) operating ~9 km to the north-northwest of the vent area since May 1997. However, JRO1 abruptly began moving downward and southward, toward a source centered about 8 km beneath the volcano, concurrently with the start of a shallow earthquake swarm on September 23, 2004. The JRO1 velocity slowed from ~0.5 millimeters per day (mm/d) in late September–early October 2004 until spring 2005. Thereafter, it was essentially constant at ~0.04 mm/d through December 2005. In similar fashion, the growth rate of the welt on the south crater floor slowed from 8.9 m 3 /s during October 4–11 to 6.4 m 3 /s during October 11-13, 2004; this trend continued after emergence of the first lava spine on October 11. The volumetric extrusion rate decreased from 5.9 m 3 /s during October 13-November 4, 2004, to 2.5 m 3 /s during December 11, 2004-January 3, 2005, and for the remainder of 2005, it was in the range 2.0-0.7 m 3 /s. Fifteen continuous GPS stations, installed soon after the eruption began, showed radially inward and downward ground motions through December 2005. Likewise, InSAR observations spanning the first year of the eruption indicate broad subsidence centered near the vent. Model-derived estimates of source-volume decrease from September 23, 2004, to October 31, 2006, are 16-24×10 6 m 3 , substantially less than the volume erupted during the same period (87×10 6 m 3 through October 21, 2006). The discrepancy can be explained by a combination of magma expansion and recharge in the source region. Lack of precursory deformation at JRO1 suggests that the conduit is poorly coupled to the rest of the edifice, so the rising magma column was able to push ahead older conduit material rather than intruding it. Constraints on conduit length and radius require that reservoir magma (as opposed to conduit-filling magma) reached the surface early during the eruption, probably soon after CO 2 emission rates peaked in early October 2004. If rapid emergence of spine 3 (the first whaleback-shaped extrusion) in late October 2004 marked the arrival of reservoir magma, then the volume of conduit material flushed from the system was about 20×10 6 m 3 --the volume of surface deformation plus spines on November 4, 2004. The corresponding radius for a cylinder extending from the surface to depth d = 5 km is 35.7 m, or 28.2 m for d = 8 km. The average ascent rate through the conduit, assuming reservoir magma began its rise on September 23, 2004, was 120 m/d for d = 5 km, or 190 m/d for d = 8 km. Observed lineal extrusion rates were 2-10 m/d, so the conduit must widen considerably near the surface. Equating magma flux through the conduit to that at the surface, we obtain a vent radius of 125 m and an extrusion rate of 5.7 m 3 /s--both values representative of the early part of the eruption. Lack of precursory inflation suggests that the volcano was poised to erupt magma already stored in a crustal reservoir when JRO1 was installed in 1997. Trilateration and campaign GPS data indicate surface dilatation, presumably caused by reservoir expansion between 1982 and 1991, but no measurable deformation between 1991 and 2003. We conclude that all three of the traditionally reliable eruption precursors (seismicity, ground deformation, and volcanic gas emission) failed to provide warning that an eruption was imminent until a few days before a visible welt appeared at the surface--a situation reminiscent of the 1980 north-flank bulge at Mount St. Helens.

Washington↗

USGS2008

Analysis of GPS-measured deformation associated with the 2004-2006 dome-building eruption of Mount St. Helens, Washington

Detecting far-field deformation at Mount St. Helens since the crater-forming landslide and blast in 1980 has been difficult despite frequent volcanic activity and improved monitoring techniques. Between 1982 and 1991, the systematic extension of line lengths in a regional GPS trilateration network is consistent with recharge of a deep magma chamber during that interval. The rate of extension, however, averages only 3 mm/yr, and some of this apparent deformation may result from systematic scale error in the electronic distance measurements. Subsequent GPS surveys and data from a continuous GPS station, located 9 km north of Mount St. Helens and operating since 1997, show no significant volcanic deformation until the start of unrest on September 23, 2004. The current eruption has been accompanied by subtle but widespread inward and downward movement of GPS monitoring stations, exponentially decreasing with time and totaling as much as 30 mm. The observed deformation is consistent with the predictions of an elastic half-space model of a vertically elongate magma chamber with its center at a depth of around 7 to 8 km and with a total cavity-volume loss of about 16–24×10 6 m 3 . The discrepancy between the estimated cavity-volume loss and the >83×10 6 -m 3 volume of the erupted dome can be explained, for the most part, by exsolution of gas in the stored magma and by minor input of new magma during the eruption.

Washington↗