Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Geophysical Research - Oceans”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 739 records · Page 41Linked to original sources
USGS2009

Mapping Curie temperature depth in the western United States with a fractal model for crustal magnetization

We have revisited the problem of mapping depth to the Curie temperature isotherm from magnetic anomalies in an attempt to provide a measure of crustal temperatures in the western United States. Such methods are based on the estimation of the depth to the bottom of magnetic sources, which is assumed to correspond to the temperature at which rocks lose their spontaneous magnetization. In this study, we test and apply a method based on the spectral analysis of magnetic anomalies. Early spectral analysis methods assumed that crustal magnetization is a completely uncorrelated function of position. Our method incorporates a more realistic representation where magnetization has a fractal distribution defined by three independent parameters: the depths to the top and bottom of magnetic sources and a fractal parameter related to the geology. The predictions of this model are compatible with radial power spectra obtained from aeromagnetic data in the western United States. Model parameters are mapped by estimating their value within a sliding window swept over the study area. The method works well on synthetic data sets when one of the three parameters is specified in advance. The application of this method to western United States magnetic compilations, assuming a constant fractal parameter, allowed us to detect robust long-wavelength variations in the depth to the bottom of magnetic sources. Depending on the geologic and geophysical context, these features may result from variations in depth to the Curie temperature isotherm, depth to the mantle, depth to the base of volcanic rocks, or geologic settings that affect the value of the fractal parameter. Depth to the bottom of magnetic sources shows several features correlated with prominent heat flow anomalies. It also shows some features absent in the map of heat flow. Independent geophysical and geologic data sets are examined to determine their origin, thereby providing new insights on the thermal and geologic crustal structure of the western United States.

Journal of Geophysical Research B: Solid Earth↗
USGS2010

Kiholo Bay, Hawaii, earthquake sequence of 2006: Relationship of the main shock slip with locations and source parameters of aftershocks

We study the source process of the Kīholo Bay earthquake ( M W 6.7), which occurred beneath the northwest part of the Island of Hawai‘i on 15 October 2006, and static stress drops of small earthquakes that occurred in 2006 and 2007 around the main shock including aftershocks. We relocate the aftershocks to determine the fault plane from the two nodal planes. The relocated aftershocks define an E‐W trending plane that dips to the south, in good agreement with one of the nodal planes given by the Global Centroid Moment Tensor solution. Waveform inversion is performed with multiple time windows to investigate the rupture speed and the slip distribution of the main shock. Waveforms of an aftershock with M W 5.2 are used to calculate empirical Green's functions. Our results indicate that the rupture propagated unilaterally to the west with a rupture speed greater than 3.0 km/s (63% of the shear wave velocity). This westward rupture is consistent with the fact that aftershocks are distributed predominantly to the west of the main shock epicenter. Most aftershocks are located on the edge of patches with a large slip, or asperities and some also occur inside the patches. We also estimate static stress drops of 39 earthquakes (2.5 < M L < 3.5) that occurred in 2006 and 2007 near the source region of the Kīholo Bay earthquake. Static stress drops range from 0.12 to 8.6 MPa and aftershocks around large slip patches of the main shock likely to have larger stress drops.

Journal of Geophysical Research B: Solid Earth↗
USGS2010

Geodetic evidence for en echelon dike emplacement and concurrent slow slip during the June 2007 intrusion and eruption at Kilauea volcano, Hawaii

A series of complex events at Kīlauea Volcano, Hawaii, 17 June to 19 June 2007, began with an intrusion in the upper east rift zone (ERZ) and culminated with a small eruption (1500 m 3 ). Surface deformation due to the intrusion was recorded in unprecedented detail by Global Positioning System (GPS) and tilt networks as well as interferometric synthetic aperture radar (InSAR) data acquired by the ENVISAT and ALOS satellites. A joint nonlinear inversion of GPS, tilt, and InSAR data yields a deflationary source beneath the summit caldera and an ENE-striking uniform-opening dislocation with ~2 m opening, a dip of ∼80° to the south, and extending from the surface to ~2 km depth. This simple model reasonably fits the overall pattern of deformation but significantly misfits data near the western end of an inferred dike-like source. Three more complex dike models are tested that allow for distributed opening including (1) a dike that follows the surface trace of the active rift zone, (2) a dike that follows the symmetry axis of InSAR deformation, and (3) two en echelon dike segments beneath mapped surface cracks and newly formed steaming areas. The en echelon dike model best fits near-field GPS and tilt data. Maximum opening of 2.4 m occurred on the eastern segment beneath the eruptive vent. Although this model represents the best fit to the ERZ data, it still fails to explain data from a coastal tiltmeter and GPS sites on Kīlauea's southwestern flank. The southwest flank GPS sites and the coastal tiltmeter exhibit deformation consistent with observations of previous slow slip events beneath Kīlauea's south flank, but inconsistent with observations of previous intrusions. Slow slip events at Kīlauea and elsewhere are thought to occur in a transition zone between locked and stably sliding zones of a fault. An inversion including slip on a basal decollement improves fit to these data and suggests a maximum of ~15 cm of seaward fault motion, comparable to previous slow-slip events.

Hawai'i↗
USGS2010

Seismic source mechanism of degassing bursts at Kilauea volcano, Hawaii: Results from waveform inversion in the 10–50 s band

The current (March 2008 to February 2009) summit eruptive activity at Kilauea Volcano is characterized by explosive degassing bursts accompanied by very long period (VLP) seismic signals. We model the source mechanisms of VLP signals in the 10–50 s band using data recorded for 15 bursts with a 10‐station broadband network deployed in the summit caldera. To determine the source centroid location and source mechanism, we minimize the residual error between data and synthetics calculated by the finite difference method for a point source embedded in a homogeneous medium that takes topography into account. The VLP signals associated with the bursts originate in a source region ∼1 km below the eastern perimeter of Halemaumau pit crater. The observed waveforms are well explained by the combination of a volumetric component and a vertical single force component. For the volumetric component, several source geometries are obtained which equally explain the observed waveforms. These geometries include (1) a pipe dipping 64° to the northeast; (2) two intersecting cracks including an east striking crack (dike) dipping 80° to the north, intersecting a north striking crack (another dike) dipping 65° to the east; (3) a pipe dipping 58° to the northeast, intersecting a crack dipping 48° to the west–southwest; and (4) a pipe dipping 57° to the northeast, intersecting a pipe dipping 58° to the west–southwest. Using the dual‐crack model as reference, the largest volume change obtained among the 15 bursts is ∼24,400 m 3 , and the maximum amplitude (peak to peak) of the force is ∼20 GN. Each burst is marked by a similar sequence of deflation and inflation, trailed by decaying oscillations of the volumetric source. The vertical force is initially upward, synchronous with source deflation, then downward, synchronous with source reinflation, followed by oscillations with polarity opposite to the volumetric oscillations. This combination of force and volume change is attributed to pressure and momentum changes induced during a fluid dynamic source mechanism involving the ascent, expansion, and burst of a large slug of gas within the upper ∼150 m of the magma conduit. As the slug expands upon approach to the surface and more liquid becomes wall supported by viscous shear forces, the pressure below the slug decreases, inducing conduit deflation and an upward force on the Earth. The final rapid slug expansion and burst stimulate VLP and LP oscillations of the conduit system, which slowly decay due to viscous dissipation and elastic radiation. Consideration of the fluid dynamic arguments leads us to prefer the dual‐crack VLP source model as it is the only candidate model capable of producing plausible values of length scales and pressure changes. The magnitudes of the vertical forces observed in the 15 bursts appear consistent with slug masses of 10 4 to 10 6 kg.

Hawaii↗
USGS2010

Parametric study of the physical properties of hydrate-bearing sand, silt, and clay sediments: 1. Electromagnetic properties

The marked decrease in bulk electrical conductivity of sediments in the presence of gas hydrates has been used to interpret borehole electrical resistivity logs and, to a lesser extent, the results of controlled source electromagnetic surveys to constrain the spatial distribution and predicted concentration of gas hydrate in natural settings. Until now, an exhaustive laboratory data set that could be used to assess the impact of gas hydrate on the electromagnetic properties of different soils (sand, silt, and clay) at different effective stress and with different saturations of hydrate has been lacking. The laboratory results reported here are obtained using a standard geotechnical cell and the hydrate-formed tetrahydrofuran (THF), a liquid that is fully miscible in water and able to produce closely controlled saturations of hydrate from dissolved phase. Both permittivity and electrical conductivity are good indicators of the volume fraction of free water in the sediment, which is in turn dependent on hydrate saturation. Permittivity in the microwave frequency range is particularly predictive of free water content since it is barely affected by ionic concentration, pore structure, and surface conduction. Electrical conductivity (or resistivity) is less reliable for constraining water content or hydrate saturation: In addition to fluid-filled porosity, other factors, such as the ionic concentration of the pore fluid and possibly other conduction effects (e.g., surface conduction in high specific surface soils having low conductivity pore fluid), also influence electrical conductivity.

Journal of Geophysical Research↗
USGS2010

Parametric study of the physical properties of hydrate‐bearing sand, silt, and clay sediments: 2. Small‐strain mechanical properties

The small‐strain mechanical properties (e.g., seismic velocities) of hydrate‐bearing sediments measured under laboratory conditions provide reference values for calibration of logging and seismic exploration results acquired in hydrate‐bearing formations. Instrumented cells were designed for measuring the compressional (P) and shear (S) velocities of sand, silts, and clay with and without hydrate and subject to vertical effective stresses of 0.01 to 2 MPa. Tetrahydrofuran (THF), which is fully miscible in water, was used as the hydrate former to permit close control over the hydrate saturation S hyd and to produce hydrate from dissolved phase, as methane hydrate forms in most natural marine settings. The results demonstrate that laboratory hydrate formation technique controls the pattern of P and S velocity changes with increasing S hyd and that the small‐strain properties of hydrate‐bearing sediments are governed by effective stress, σ ′ v and sediment specific surface. The S velocity increases with hydrate saturation owing to an increase in skeletal shear stiffness, particularly when hydrate saturation exceeds S hyd ≈ 0.4. At very high hydrate saturations, the small strain shear stiffness is determined by the presence of hydrates and becomes insensitive to changes in effective stress. The P velocity increases with hydrate saturation due to the increases in both the shear modulus of the skeleton and the bulk modulus of pore‐filling phases during fluid‐to‐hydrate conversion. Small‐strain Poisson's ratio varies from 0.5 in soft sediments lacking hydrates to 0.25 in stiff sediments (i.e., subject to high vertical effective stress or having high S hyd ). At S hyd ≥ 0.5, hydrate hinders expansion and the loss of sediment stiffness during reduction of vertical effective stress, meaning that hydrate‐rich natural sediments obtained through pressure coring should retain their in situ fabric for some time after core retrieval if the cores are maintained within the hydrate stability field.

Journal of Geophysical Research↗
USGS2010

The 2007 M5.4 Alum Rock, California, earthquake: Implications for future earthquakes on the central and southern Calaveras Fault

[1] The similarity of seismograms recorded by two seismic stations demonstrate that the 31 October 2007 moment magnitude M 5.4 Alum Rock earthquake is a repeat of a 1955 M L 5.5 earthquake. Both occurred on Oppenheimer et al.'s (1990) Zone V “stuck patch” on the central Calaveras fault, providing new support for their model of Calaveras fault earthquake activity. We suggest that Zone V fails only in a family of recurring M ∼ 5.4–5.5 earthquakes. The 1955 and 2007 earthquakes are the penultimate and ultimate Zone V events. Earthquakes in 1891 and 1864 are possible earlier Zone V events. The next Zone V event is not expected in the next few decades, assuming a time‐dependent recurrence model: the mean forecast date is 2064 (2035–2104, 95% confidence range). We further suggest that Zones I, II, III, and IV fail in recurring M ∼ 5.1–5.3, M ∼ 5.6–5.8, M ∼ 6.1–6.3, and M ∼ 4.9–5.0 earthquakes, respectively. If our earthquake recurrence model is correct, the next Zone I event is overdue and could occur anytime, and M 5–6 earthquakes should not occur on Zones II, III, and IV before 2014, 2012, and 2026, respectively. We cannot rule out the possibility that Zone VI, which lies at the southern end of the Mission Seismic Trend, where the southern Hayward and central Calaveras faults appear to connect at depth, fails aseismically or in large events on the southern Hayward fault, such as last occurred in 1868, or in large events on the adjoining northern Calaveras fault segment.

California↗
USGS2010

Supraslab earthquake clusters above the subduction plate boundary offshore Sanriku, northeastern Japan: Seismogenesis in a graveyard of detached seamounts?

Thousands of offshore repeating earthquakes with low‐angle thrust focal mechanisms occur along the subduction plate boundary of NE Japan. Double‐difference relocation methods using P ‐ and S ‐wave arrivals reveal clusters of events above these repeating events. To assure good depth control we restrict our study to events that are close to seismic stations. These “supraslab” earthquake clusters are regional features at depths of 25 to 50 km, and most of these clusters are below the depth of the forearc Moho, which we determined from converted waves. Seismicity over this depth range does not occur under the inland area of NE Japan except just below the vicinity of the arc volcanoes. Re‐entrants in the inner trench slope indicate that repeated collisions of seamounts have occurred in the past. Our preliminary interpretation of supraslab clusters is that they represent seismicity in seamounts detached from the Pacific plate during slab descent, driven by the resistance of seamounts to subduction. Detachment during slab descent probably occurs on the sedimented and hydrothermally altered seafloor on which seamounts were originally constructed since these are known as zones of weakness during active island growth. High fluid pressure produced during dehydration of clay minerals and other low‐temperature hydrous minerals could enable detachment at depths. Seamount crust is thus accreted to forearcs, possibly leading to a long‐term component of near‐coastal uplift. Supraslab earthquake clusters may be our most direct evidence of the fates of seamounts and suggest that tectonic underplating is actively occurring in this subduction system.

Northeastern Japan↗
USGS2010

Coherence of Mach fronts during heterogeneous supershear earthquake rupture propagation: Simulations and comparison with observations

We study how heterogeneous rupture propagation affects the coherence of shear and Rayleigh Mach wavefronts radiated by supershear earthquakes. We address this question using numerical simulations of ruptures on a planar, vertical strike-slip fault embedded in a three-dimensional, homogeneous, linear elastic half-space. Ruptures propagate spontaneously in accordance with a linear slip-weakening friction law through both homogeneous and heterogeneous initial shear stress fields. In the 3-D homogeneous case, rupture fronts are curved owing to interactions with the free surface and the finite fault width; however, this curvature does not greatly diminish the coherence of Mach fronts relative to cases in which the rupture front is constrained to be straight, as studied by Dunham and Bhat (2008a). Introducing heterogeneity in the initial shear stress distribution causes ruptures to propagate at speeds that locally fluctuate above and below the shear wave speed. Calculations of the Fourier amplitude spectra (FAS) of ground velocity time histories corroborate the kinematic results of Bizzarri and Spudich (2008a): (1) The ground motion of a supershear rupture is richer in high frequency with respect to a subshear one. (2) When a Mach pulse is present, its high frequency content overwhelms that arising from stress heterogeneity. Present numerical experiments indicate that a Mach pulse causes approximately an ω −1.7 high frequency falloff in the FAS of ground displacement. Moreover, within the context of the employed representation of heterogeneities and over the range of parameter space that is accessible with current computational resources, our simulations suggest that while heterogeneities reduce peak ground velocity and diminish the coherence of the Mach fronts, ground motion at stations experiencing Mach pulses should be richer in high frequencies compared to stations without Mach pulses. In contrast to the foregoing theoretical results, we find no average elevation of 5%-damped absolute response spectral accelerations (SA) in the period band 0.05–0.4 s observed at stations that presumably experienced Mach pulses during the 1979 Imperial Valley, 1999 Kocaeli, and 2002 Denali Fault earthquakes compared to SA observed at non-Mach pulse stations in the same earthquakes. A 20% amplification of short period SA is seen only at a few of the Imperial Valley stations closest to the fault. This lack of elevated SA suggests that either Mach pulses in real earthquakes are even more incoherent that in our simulations or that Mach pulses are vulnerable to attenuation through nonlinear soil response. In any case, this result might imply that current engineering models of high frequency earthquake ground motions do not need to be modified by more than 20% close to the fault to account for Mach pulses, provided that the existing data are adequately representative of ground motions from supershear earthquakes.

Journal of Geophysical Research B: Solid Earth↗
USGS2010

Glimpses of East Antarctica: Aeromagnetic and satellite magnetic view from the central Transantarctic Mountains of East Antarctica

Aeromagnetic and satellite magnetic data provide glimpses of the crustal architecture within the Ross Sea sector of the enigmatic, ice-covered East Antarctic shield critical for understanding both global tectonic and climate history. In the central Transantarctic Mountains (CTAM), exposures of Precambrian basement, coupled with new high-resolution magnetic data, other recent aeromagnetic transects, and satellite magnetic and seismic tomography data, show that the shield in this region comprises an Archean craton modified both by Proterozoic magmatism and early Paleozoic orogenic basement reactivation. CTAM basement structures linked to the Ross Orogeny are imaged 50–100 km farther west than previously mapped, bounded by inboard upper crustal Proterozoic granites of the Nimrod igneous province. Magnetic contrasts between craton and rift margin sediments define the Neoproterozoic rift margin, likely reactivated during Ross orogenesis and Jurassic extension. Interpretation of satellite magnetic and aeromagnetic patterns suggests that the Neoproterozoic rift margin of East Antarctica is offset by transfer zones to form a stepwise series of salients tracing from the CTAM northward through the western margin of the Wilkes Subglacial Basin to the coast at Terre Adélie. Thinned Precambrian crust inferred to lie east of the rift margin cannot be imaged magnetically because of modification by Neoproterozoic and younger tectonic events.

Journal of Geophysical Research B: Solid Earth↗
USGS2010

Postseismic relaxation following the 1989 MS7.1 Loma Prieta earthquake, central California

The postseismic relaxation (postseismic displacement less displacement that would have occurred at the preseismic rate) measured by GPS and leveling following the 1989 M S 7.1 Loma Prieta earthquake is reexamined. The temporal dependence of the relaxation over the first 1200 days postseismic is well described by 1 − e − t / τ , where τ = 414 ± 92 days. (Quoted uncertainties are standard deviations.) That temporal dependence appears to be a linear function of the cumulative number of M > 2.5 aftershocks that have occurred. The relaxation is attributed to afterslip (1.56 ± 0.20 m dextral strike slip and 0.60 ± 0.04 m reverse slip) on the downdip extension (depth, 16–21 km) of the coseismic rupture plus a collapse (0.11 ± 0.02 m fault-normal displacement) of the rupture zone (depth, 5–16 km). Because the postseismic uplift was determined by leveling over a route with relief in excess of 1000 m, an allowance (18 ± 2 ppm of height above the base elevation) for excess (beyond corrections already applied) height-dependent error in the measured uplift was estimated simultaneously with the afterslip and collapse parameters. This new solution for afterslip and collapse on the plane of the rupture provides an alternative explanation to the suggestion by Bürgmann et al. (1997) that the postseismic deformation was due to afterslip on the coseismic rupture and on a nearby, shallow, 30° dipping thrust.

California↗
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

Ground surface deformation patterns, magma supply, and magma storage at Okmok volcano, Alaska, from InSAR analysis: 2. Coeruptive deflation, July-August 2008

A hydrovolcanic eruption near Cone D on the floor of Okmok caldera, Alaska, began on 12 July 2008 and continued until late August 2008. The eruption was preceded by inflation of a magma reservoir located beneath the center of the caldera and ∼3 km below sea level (bsl), which began immediately after Okmok's previous eruption in 1997. In this paper we use data from several radar satellites and advanced interferometric synthetic aperture radar (InSAR) techniques to produce a suite of 2008 coeruption deformation maps. Most of the surface deformation that occurred during the eruption is explained by deflation of a Mogi-type source located beneath the center of the caldera and 2–3 km bsl, i.e., essentially the same source that inflated prior to the eruption. During the eruption the reservoir deflated at a rate that decreased exponentially with time with a 1/ e time constant of ∼13 days. We envision a sponge-like network of interconnected fractures and melt bodies that in aggregate constitute a complex magma storage zone beneath Okmok caldera. The rate at which the reservoir deflates during an eruption may be controlled by the diminishing pressure difference between the reservoir and surface. A similar mechanism might explain the tendency for reservoir inflation to slow as an eruption approaches until the pressure difference between a deep magma production zone and the reservoir is great enough to drive an intrusion or eruption along the caldera ring-fracture system.

Alaska↗
USGS2010

Anisotropy, repeating earthquakes, and seismicity associated with the 2008 eruption of Okmok Volcano, Alaska

We use shear wave splitting (SWS) analysis and double-difference relocation to examine temporal variations in seismic properties prior to and accompanying magmatic activity associated with the 2008 eruption of Okmok volcano, Alaska. Using bispectrum cross-correlation, a multiplet of 25 earthquakes is identified spanning five years leading up to the eruption, each event having first motions compatible with a normal fault striking NE–SW. Cross-correlation differential times are used to relocate earthquakes occurring between January 2003 and February 2009. The bulk of the seismicity prior to the onset of the eruption on 12 July 2008 occurred southwest of the caldera beneath a geothermal field. Earthquakes associated with the onset of the eruption occurred beneath the northern portion of the caldera and started as deep as 13 km. Subsequent earthquakes occurred predominantly at 3 km depth, coinciding with the depth at which the magma body has been modeled using geodetic data. Automated SWS analysis of the Okmok catalog reveals radial polarization outside the caldera and a northwest-southeast polarization within. We interpret these polarizations in terms of a magma reservoir near the center of the caldera, which we model with a Mogi point source. SWS analysis using the same input processing parameters for each event in the multiplet reveals no temporal changes in anisotropy over the duration of the multiplet, suggesting either a short-term or small increase in stress just before the eruption that was not detected by GPS, or eruption triggering by a mechanism other than a change of stress in the system.

Alaska↗
USGS2010

Lessons from (triggered) tremor

I test a “clock-advance” model that implies triggered tremor is ambient tremor that occurs at a sped-up rate as a result of loading from passing seismic waves. This proposed model predicts that triggering probability is proportional to the product of the ambient tremor rate and a function describing the efficacy of the triggering wave to initiate a tremor event. Using data mostly from Cascadia, I have compared qualitatively a suite of teleseismic waves that did and did not trigger tremor with ambient tremor rates. Many of the observations are consistent with the model if the efficacy of the triggering wave depends on wave amplitude. One triggered tremor observation clearly violates the clock-advance model. The model prediction that larger triggering waves result in larger triggered tremor signals also appears inconsistent with the measurements. I conclude that the tremor source process is a more complex system than that described by the clock-advance model predictions tested. Results of this and previous studies also demonstrate that (1) conditions suitable for tremor generation exist in many tectonic environments, but, within each, only occur at particular spots whose locations change with time; (2) any fluid flow must be restricted to less than a meter; (3) the degree to which delayed failure and secondary triggering occurs is likely insignificant; and 4) both shear and dilatational deformations may trigger tremor. Triggered and ambient tremor rates correlate more strongly with stress than stressing rate, suggesting tremor sources result from time-dependent weakening processes rather than simple Coulomb failure.

Cascadia↗
USGS2010

Volcano monitoring using GPS: Developing data analysis strategies based on the June 2007 Kīlauea Volcano intrusion and eruption

The global positioning system (GPS) is one of the most common techniques, and the current state of the art, used to monitor volcano deformation. In addition to slow (several centimeters per year) displacement rates, GPS can be used to study eruptions and intrusions that result in much larger (tens of centimeters over hours-days) displacements. It is challenging to resolve precise positions using GPS at subdaily time intervals because of error sources such as multipath and atmospheric refraction. In this paper, the impact of errors due to multipath and atmospheric refraction at subdaily periods is examined using data from the GPS network on Kīlauea Volcano, Hawai'i. Methods for filtering position estimates to enhance precision are both simulated and tested on data collected during the June 2007 intrusion and eruption. Comparisons with tiltmeter records show that GPS instruments can precisely recover the timing of the activity.

Hawai'i↗
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↗