USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research, Solid Earth”

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 559 records · Page 31Linked to original sources
USGS2008

Finding concealed active faults: Extending the southern Whidbey Island fault across the Puget Lowland, Washington

The southern Whidbey Island fault zone (SWIF), as previously mapped using borehole data, potential field anomalies, and marine seismic reflection surveys, consists of three subparallel, northwest trending strands extending ∼100 km from near Vancouver Island to the northern Puget Lowland. East of Puget Sound, the SWIF makes landfall between the cities of Seattle and Everett but is concealed beneath a thick mantle of young glacial deposits and vegetation. A ∼20‐km‐wide, northwest trending swath of subparallel, low‐amplitude aeromagnetic anomalies crosses this region of the Puget Lowland and is on strike with the SWIF. The most prominent aeromagnetic anomaly, the Cottage Lake lineament, extends at least 18 km and lies approximately on strike with the SWIF on Whidbey Island. Subtle scarps and topographic lineaments on Pleistocene surfaces, visible on high‐resolution lidar topography at a number of locations along the SWIF, lie on or near these magnetic anomalies. In the field, scarps exhibit northeast‐side‐up and vertical relief of 1 to 5 m. Excavations across several lidar scarps lying on or near magnetic anomalies show evidence for multiple folding and faulting events since deglaciation, most likely above buried reverse/oblique faults. Excavations in areas away from magnetic anomalies do not show evidence of tectonic deformation. In total, paleoseismological evidence suggests that the SWIF produced at least four earthquakes since deglaciation about 16,400 years ago, the most recent less than 2700 years ago.

Washington↗
USGS2008

Variability of passive gas emissions, seismicity, and deformation during crater lake growth at White Island Volcano, New Zealand, 2002-2006

We report on 4 years of airborne measurements of CO2, SO2, and H2S emission rates during a quiescent period at White Island volcano, New Zealand, beginning in 2003. During this time a significant crater lake emerged, allowing scrubbig processes to be investigated. CO2 emissions varied from a baseline of 250 to >2000 t d-1 and demonstrated clear annual cycling that was consistent with numbers of earthquake detections and annual changes in sea level. The annual variability was found to be most likely related to increases in the strain on the volcano during sea level highs, temporarily causing fractures to reduce in size in the upper conduit. SO2 emissions varied from 0 to >400 t d-1 and were clearly affected by scrubbing processes within the first year of take development. Scrubbing caused increases of SO42- and Cl- in lake waters, and the ratio of carbon to total sulphur suggested that elemental sulphur deposition was also significant in the lake during the first year. Careful measurements of the lake level and chemistry allowed estimates of the rate of H2O(g) and HCl(g) input into the lake and suggested that the molar abundances of major gas species (H2O, CO2, SO2, and HCl) during this quiescent phase were similar to fumarolic ratios observed between earlier eruptive periods. The volume of magma estimated from CO2 emissions (0.0 15-0.04 km3) was validated by Cl- increases in the lake, suggesting that the gas and magma are transported from deep to shallow depths as a closed system and likely become open in the upper conduit region. The absence of surface deformation further leads to a necessity of magma convection to supply and remove magma from the degassing depths. Two models of convection configurations are discussed. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Dislocation models of interseismic deformation in the western United States

The GPS-derived crustal velocity field of the western United States is used to construct dislocation models in a viscoelastic medium of interseismic crustal deformation. The interseismic velocity field is constrained by 1052 GPS velocity vectors spanning the ???2500-km-long plate boundary zone adjacent to the San Andreas fault and Cascadia subduction zone and extending ???1000 km into the plate interior. The GPS data set is compiled from U.S. Geological Survey campaign data, Plate Boundary Observatory data, and the Western U.S. Cordillera velocity field of Bennett et al. (1999). In the context of viscoelastic cycle models of postearthquake deformation, the interseismic velocity field is modeled with a combination of earthquake sources on ???100 known faults plus broadly distributed sources. Models that best explain the observed interseismic velocity field include the contributions of viscoelastic relaxation from faulting near the major plate margins, viscoelastic relaxation from distributed faulting in the plate interior, as well as lateral variations in depth-averaged rigidity in the elastic lithosphere. Resulting rigidity variations are consistent with reduced effective elastic plate thickness in a zone a few tens of kilometers wide surrounding the San Andreas fault (SAF) system. Primary deformation characteristics are captured along the entire SAF system, Eastern California Shear Zone, Walker Lane, the Mendocino triple junction, the Cascadia margin, and the plate interior up to ???1000 km from the major plate boundaries.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

In situ Raman spectroscopic investigation of the structure of subduction-zone fluids

In situ Raman spectra of synthetic subduction-zone fluids (KAlSi3O8-H2O system) were measured to 900?? and 2.3 GPa using a hydrothermal diamond-anvil cell. The structures of aqueous fluid and hydrous melt become closer when conditions approach the second critical endpoint. Almost no three-dimensional network was observed in the supercritical fluid above 2 GPa although a large amount of silicate component is dissolved, suggesting that the physical and chemical properties of these phases change drastically at around the second critical endpoint. Our experimental results indicate that the fluids released from a subducting slab change from aqueous fluid to supercritical fluid with increasing depth under the volcanic arcs. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

A model of earthquake triggering probabilities and application to dynamic deformations constrained by ground motion observations

We have used observations from Felzer and Brodsky (2006) of the variation of linear aftershock densities (i.e., aftershocks per unit length) with the magnitude of and distance from the main shock fault to derive constraints on how the probability of a main shock triggering a single aftershock at a point, P(r, D), varies as a function of distance, r, and main shock rupture dimension, D. We find that P(r, D) becomes independent of D as the triggering fault is approached. When r ??? D P(r, D) scales as Dm where m-2 and decays with distance approximately as r-n with n = 2, with a possible change to r-(n-1) at r > h, where h is the closest distance between the fault and the boundaries of the seismogenic zone. These constraints may be used to test hypotheses about the types of deformations and mechanisms that trigger aftershocks. We illustrate this using dynamic deformations (i.e., radiated seismic waves) and a posited proportionality with P(r, D). Deformation characteristics examined include peak displacements, peak accelerations and velocities (proportional to strain rates and strains, respectively), and two measures that account for cumulative deformations. Our model indicates that either peak strains alone or strain rates averaged over the duration of rupture may be responsible for aftershock triggering.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Probabilistic seismic hazard in the San Francisco Bay area based on a simplified viscoelastic cycle model of fault interactions

We construct a viscoelastic cycle model of plate boundary deformation that includes the effect of time-dependent interseismic strain accumulation, coseismic strain release, and viscoelastic relaxation of the substrate beneath the seismogenic crust. For a given fault system, time-averaged stress changes at any point (not on a fault) are constrained to zero; that is, kinematic consistency is enforced for the fault system. The dates of last rupture, mean recurrence times, and the slip distributions of the (assumed) repeating ruptures are key inputs into the viscoelastic cycle model. This simple formulation allows construction of stress evolution at all points in the plate boundary zone for purposes of probabilistic seismic hazard analysis (PSHA). Stress evolution is combined with a Coulomb failure stress threshold at representative points on the fault segments to estimate the times of their respective future ruptures. In our PSHA we consider uncertainties in a four-dimensional parameter space: the rupture peridocities, slip distributions, time of last earthquake (for prehistoric ruptures) and Coulomb failure stress thresholds. We apply this methodology to the San Francisco Bay region using a recently determined fault chronology of area faults. Assuming single-segment rupture scenarios, we find that fature rupture probabilities of area faults in the coming decades are the highest for the southern Hayward, Rodgers Creek, and northern Calaveras faults. This conclusion is qualitatively similar to that of Working Group on California Earthquake Probabilities, but the probabilities derived here are significantly higher. Given that fault rupture probabilities are highly model-dependent, no single model should be used to assess to time-dependent rupture probabilities. We suggest that several models, including the present one, be used in a comprehensive PSHA methodology, as was done by Working Group on California Earthquake Probabilities.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Eruption dynamics at Mount St. Helens imaged from broadband seismic waveforms: Interaction of the shallow magmatic and hydrothermal systems

The current eruption at Mount St. Helens is characterized by dome building and shallow, repetitive, long-period (LP) earthquakes. Waveform cross-correlation reveals remarkable similarity for a majority of the earthquakes over periods of several weeks. Stacked spectra of these events display multiple peaks between 0.5 and 2 Hz that are common to most stations. Lower-amplitude very-long-period (VLP) events commonly accompany the LP events. We model the source mechanisms of LP and VLP events in the 0.5-4 s and 8-40 s bands, respectively, using data recorded in July 2005 with a 19-station temporary broadband network. The source mechanism of the LP events includes: 1) a volumetric component modeled as resonance of a gently NNW-dipping, steam-filled crack located directly beneath the actively extruding part of the new dome and within 100 m of the crater floor and 2) a vertical single force attributed to movement of the overlying dome. The VLP source, which also includes volumetric and single-force components, is 250 m deeper and NNW of the LP source, at the SW edge of the 1980s lava dome. The volumetric component points to the compression and expansion of a shallow, magma-filled sill, which is subparallel to the hydrothermal crack imaged at the LP source, coupled with a smaller component of expansion and compression of a dike. The single-force components are due to mass advection in the magma conduit. The location, geometry and timing of the sources suggest the VLP and LP events are caused by perturbations of a common crack system.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Physical property changes in hydrate-bearing sediment due to depressurization and subsequent repressurization

Physical property measurements of sediment cores containing natural gas hydrate are typically performed on material exposed, at least briefly, to non-in situ conditions during recovery. To examine the effects of a brief excursion from the gas-hydrate stability field, as can occur when pressure cores are transferred to pressurized storage vessels, we measured physical properties on laboratory-formed sand packs containing methane hydrate and methane pore gas. After depressurizing samples to atmospheric pressure, we repressurized them into the methane-hydrate stability field and remeasured their physical properties. Thermal conductivity, shear strength, acoustic compressional and shear wave amplitudes, and speeds of the original and depressurized/repressurized samples are compared. X– ray computed tomography images track how the gas-hydrate distribution changes in the hydrate-cemented sands owing to the depressurizaton/repressurization process. Because depressurization-induced property changes can be substantial and are not easily predicted, particularly in water-saturated, hydrate-bearing sediment, maintaining pressure and temperature conditions throughout the core recovery and measurement process is critical for using laboratory measurements to estimate in situ properties.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Observations and interpretation of fundamental mode Rayleigh wavefields recorded by the Transportable Array (USArray)

Broadband recordings of the dense Transportable Array (TA) in the western United States provide unparalleled detailed images of long-period seismic surface wavefields. With 400 stations spanning most of the western United States, wavefronts of fundamental mode Rayleigh waves may be visualized coherently across the array at periods ???40 s. In order to constrain the Rayleigh wave phase velocity structure in the western United States, I assemble a data set of vertical component seismograms from 53 teleseismic events recorded by the TA from April 2006 to October 2007. Complex amplitude spectra from these recordings at peni ods 27-100 s are interpreted using the multiplane wave tomographic method of Friederich and Wielandt (1995) and Pollitz (1999). This analysis yields detailed surface wave phase velocity and three-dimensional shear wave velocity patterns across the North American plate boundary zone, elucidating the active processes in the highly heterogeneous western U.S. upper mantle.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Postearthquake relaxation after the 2004 M6 Parkfield, California, earthquake and rate-and-state friction

An unusually complete set of measurements (including rapid rate GPS over the first 10 days) of postseismic deformation is available at 12 continuous GPS stations located close to the epicenter of the 2004 M6.0 Parkfield earthquake. The principal component modes for the relaxation of the ensemble of those 12 GPS stations were determined. The first mode alone furnishes an adequate approximation to the data. Thus, the relaxation at all stations can be represented by the product of a common temporal function and distinct amplitudes for each component (north or east) of relaxation at each station. The distribution in space of the amplitudes indicates that the relaxation is dominantly strike slip. The temporal function, which spans times from about 5 min to 900 days postearthquake, can be fit by a superposition of three creep terms, each of the form ??l loge(1 + t/??l), with characteristic times ??, = 4.06, 0.11, and 0.0001 days. It seems likely that what is actually involved is a broad spectrum of characteristic times, the individual components of which arise from afterslip on different fault patches. Perfettini and Avouac (2004) have shown that an individual creep term can be explained by the spring-slider model with rate-dependent (no state variable) friction. The observed temporal function can also be explained using a single spring-slider model (i.e., single fault patch) that includes rate-and-state-dependent friction, a single-state variable, and either of the two commonly used (aging and slip) state evolution laws. In the latter fits, the rate-and-state friction parameter b is negative.

Journal of Geophysical Research B: Solid Earth↗
USGS2009

Constraining earthquake source inversions with GPS data: 2. A two-step approach to combine seismic and geodetic data sets

We present a new method to combine static and wavefield data to image earthquake ruptures. Our combined inversion is a two-step procedure, following the work of Hernandez et al. (1999), and takes into account the differences between the resolutions of the two data sets. The first step consists of an inversion of the static field, which yields a map of slip amplitude. This inversion exploits a special irregular grid that takes into account the resolution of the static data. The second step is an inversion of the radiated wavefield; it results in the determination of the time evolution of slip on the fault. In the second step, the slip amplitude is constrained to resemble the static slip amplitude map inferred from the GPS inversion. Using this combined inversion, we study the source process of the 2004 M6 Parkfield, California, earthquake. We conclude that slip occurred in two main regions of the fault, each of which displayed distinct rupture behaviors. Slip initiated at the hypocenter with a very strong bilateral burst of energy. Here, slip was localized in a narrow area approximately 10 km long, the rupture velocity was very fast (???3.5 km/s), and slip only lasted a short period of time (<1 s). Then the rupture proceeded to a wider region 12-20 km northwest of the hypocenter. Here, the earthquake developed in a more moderated way: the rupture velocity slowed to ???3.0 km/s and slip lasted longer (1-2 s). The maximum slip amplitude was 0.45 m. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2009

Premonitory acoustic emissions and stick-slip in natural and smooth-faulted Westerly granite

A stick-slip event was induced in a cylindrical sample of Westerly granite containing a preexisting natural fault by loading at constant confining pressure of 150 MPa. Continuously recorded acoustic emission (AE) data and computer tomography (CT)-generated images of the fault plane were combined to provide a detailed examination of microscale processes operating on the fault. The dynamic stick-slip event, considered to be a laboratory analog of an earthquake, generated an ultrasonic signal that was recorded as a large-amplitude AE event. First arrivals of this event were inverted to determine the nucleation site of slip, which is associated with a geometric asperity on the fault surface. CT images and AE locations suggest that a variety of asperities existed in the sample because of the intersection of branch or splay faults with the main fault. This experiment is compared with a stick-slip experiment on a sample prepared with a smooth, artificial saw-cut fault surface. Nearly a thousand times more AE were observed for the natural fault, which has a higher friction coefficient (0.78 compared to 0.53) and larger shear stress drop (140 compared to 68 MPa). However at the measured resolution, the ultrasonic signal emitted during slip initiation does not vary significantly between the two experiments, suggesting a similar dynamic rupture process. We propose that the natural faulted sample under triaxial compression provides a good laboratory analogue for a field-scale fault system in terms of the presence of asperities, fault surface heterogeneity, and interaction of branching faults. ?? 2009.

Journal of Geophysical Research B: Solid Earth↗
USGS2009

Regional three-dimensional seismic velocity model of the crust and uppermost mantle of northern California

We present a three-dimensional (3D) tomographic model of the P wave velocity (Vp) structure of northern California. We employed a regional-scale double-difference tomography algorithm that incorporates a finite-difference travel time calculator and spatial smoothing constraints. Arrival times from earthquakes and travel times from controlled-source explosions, recorded at network and/or temporary stations, were inverted for Vp on a 3D grid with horizontal node spacing of 10 to 20 km and vertical node spacing of 3 to 8 km. Our model provides an unprecedented, comprehensive view of the regional-scale structure of northern California, putting many previously identified features into a broader regional context and improving the resolution of a number of them and revealing a number of new features, especially in the middle and lower crust, that have never before been reported. Examples of the former include the complex subducting Gorda slab, a steep, deeply penetrating fault beneath the Sacramento River Delta, crustal low-velocity zones beneath Geysers-Clear Lake and Long Valley, and the high-velocity ophiolite body underlying the Great Valley. Examples of the latter include mid-crustal low-velocity zones beneath Mount Shasta and north of Lake Tahoe. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2009

Inducing in situ, nonlinear soil response applying an active source

[1] It is well known that soil sites have a profound effect on ground motion during large earthquakes. The complex structure of soil deposits and the highly nonlinear constitutive behavior of soils largely control nonlinear site response at soil sites. Measurements of nonlinear soil response under natural conditions are critical to advancing our understanding of soil behavior during earthquakes. Many factors limit the use of earthquake observations to estimate nonlinear site response such that quantitative characterization of nonlinear behavior relies almost exclusively on laboratory experiments and modeling of wave propagation. Here we introduce a new method for in situ characterization of the nonlinear behavior of a natural soil formation using measurements obtained immediately adjacent to a large vibrator source. To our knowledge, we are the first group to propose and test such an approach. Employing a large, surface vibrator as a source, we measure the nonlinear behavior of the soil by incrementally increasing the source amplitude over a range of frequencies and monitoring changes in the output spectra. We apply a homodyne algorithm for measuring spectral amplitudes, which provides robust signal-to-noise ratios at the frequencies of interest. Spectral ratios are computed between the receivers and the source as well as receiver pairs located in an array adjacent to the source, providing the means to separate source and near-source nonlinearity from pervasive nonlinearity in the soil column. We find clear evidence of nonlinearity in significant decreases in the frequency of peak spectral ratios, corresponding to material softening with amplitude, observed across the array as the source amplitude is increased. The observed peak shifts are consistent with laboratory measurements of soil nonlinearity. Our results provide constraints for future numerical modeling studies of strong ground motion during earthquakes.

Journal of Geophysical Research B: Solid Earth↗
USGS2009

Tracking magma volume recovery at Okmok Volcano using GPS and an unscented kalman filter

Changes beneath a volcano can be observed through position changes in a GPS network, but distinguishing the source of site motion is not always straightforward. The records of continuous GPS sites provide a favorable data set for tracking magma migration. Dense campaign observations usually provide a better spatial picture of the overall deformation field, at the expense of an episodic temporal record. Combining these observations provides the best of both worlds. A Kalman filter provides a means for integrating discrete and continuous measurements and for interpreting subtle signals. The unscented Kalman filter (UKF) is a nonlinear method for time‐dependent observations. We demonstrate the application of this technique to deformation data by applying it to GPS data collected at Okmok volcano. Seven years of GPS observations at Okmok are analyzed using a Mogi source model and the UKF. The deformation source at Okmok is relatively stable at 2.5 km depth below sea level, located beneath the center of the caldera, which means the surface deformation is caused by changes in the strength of the source. During the 7 years of GPS observations more than 0.5 m of uplift has occurred, a majority of that during the time period January 2003 to July 2004. The total volume recovery at Okmok since the last eruption in 1997 is ∼60–80%. The UKF allows us to solve simultaneously for the time‐dependence of the source strength and for the location without a priori information about the source. 2009 by the American Geophysical Union.

Alaska↗
USGS2009

Seismic wave triggering of nonvolcanic tremor, episodic tremor and slip, and earthquakes on Vancouver Island

[1] We explore the physical conditions that enable triggering of nonvolcanic tremor and earthquakes by considering local seismic activity on Vancouver Island, British Columbia during and immediately after the arrival of large-amplitude seismic waves from 30 teleseismic and 17 regional or local earthquakes. We identify tremor triggered by four of the teleseismic earthquakes. The close temporal and spatial proximity of triggered tremor to ambient tremor and aseismic slip indicates that when a fault is close to or undergoing failure, it is particularly susceptible to triggering of further events. The amplitude of the triggering waves also influences the likelihood of triggering both tremor and earthquakes such that large amplitude waves triggered tremor in the absence of detectable aseismic slip or ambient tremor. Tremor and energy radiated from regional/local earthquakes share the same frequency passband so that tremor cannot be identified during these smaller, more frequent events. We confidently identify triggered local earthquakes following only one teleseism, that with the largest amplitude, and four regional or local events that generated vigorous aftershock sequences in their immediate vicinity. Earthquakes tend to be triggered in regions different from tremor and with high ambient seismicity rates. We also note an interesting possible correlation between large teleseismic events and episodic tremor and slip (ETS) episodes, whereby ETS events that are “late” and have built up more stress than normal are susceptible to triggering by the slight nudge of the shaking from a large, distant event, while ETS events that are “early” or “on time” are not.

Vancouver Island↗
USGS2009

Volcano-tectonic implications of 3-D velocity structures derived from joint active and passive source tomography of the island of Hawaii

We present a velocity model of the onshore and offshore regions around the southern part of the island of Hawaii, including southern Mauna Kea, southeastern Hualalai, and the active volcanoes of Mauna Loa, and Kilauea, and Loihi seamount. The velocity model was inverted from about 200,000 first-arrival traveltime picks of earthquakes and air gun shots recorded at the Hawaiian Volcano Observatory (HVO). Reconstructed volcanic structures of the island provide us with an improved understanding of the volcano-tectonic evolution of Hawaiian volcanoes and their interactions. The summits and upper rift zones of the active volcanoes are characterized by high-velocity materials, correlated with intrusive magma cumulates. These high-velocity materials often do not extend the full lengths of the rift zones, suggesting that rift zone intrusions may be spatially limited. Seismicity tends to be localized seaward of the most active intrusive bodies. Low-velocity materials beneath parts of the active rift zones of Kilauea and Mauna Loa suggest discontinuous rift zone intrusives, possibly due to the presence of a preexisting volcanic edifice, e.g., along Mauna Loa beneath Kilauea's southwest rift zone, or alternatively, removal of high-velocity materials by large-scale landsliding, e.g., along Mauna Loa's western flank. Both locations also show increased seismicity that may result from edifice interactions or reactivation of buried faults. New high-velocity regions are recognized and suggest the presence of buried, and in some cases, previously unknown rift zones, within the northwest flank of Mauna Loa, and the south flanks of Mauna Loa, Hualalai, and Mauna Kea. Copyright 2009 by the American Geophysical Union.

Hawaii↗
USGS2008

Seismicity and deformation induced by magma accumulation at three basaltic volcanoes

We analyzed the evolution of volcano‐tectonic (VT) seismicity and deformation at three basaltic volcanoes (Kilauea, Mauna Loa, Piton de la Fournaise) during phases of magma accumulation. We observed that the VT earthquake activity displays an accelerating evolution at the three studied volcanoes during the time of magma accumulation. At the same times, deformation rates recorded at the summit of Kilauea and Mauna Loa volcanoes were not accelerating but rather tend to decay. To interpret these observations, we propose a physical model describing the evolution of pressure produced by the accumulation of magma into a reservoir. This variation of pressure is then used to force a simple model of damage, where damage episodes are equivalent to earthquakes. This model leads to an exponential increase of the VT activity and to an exponential decay of the deformation rate during accumulation phases. Seismicity and deformation data are well fitted by such an exponential model. The time constant, deduced from the exponential increase of the seismicity, is in agreement with the time constant predicted by the model of magma accumulation. This VT activity can thus be a direct indication of the accumulation of magma at depth, and therefore can be seen as a long‐term precursory phenomenon, at least for the three studied basaltic volcanoes. Unfortunately, it does not allow the prediction of the onset of future eruptions, as no diverging point (i.e., critical time) is present in the model.

Journal of Geophysical Research B: Solid Earth↗