SEARCH · USGS Science

Results for “Journal of Geophysical Research: Earth Surface”

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

SourceResearch connectionFrom yearThrough year

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

At least 325 records · Page 18Linked to original sources

USGS2009

On near-source earthquake triggering

[1] When one earthquake triggers others nearby, what connects them? Two processes are observed: static stress change from fault offset and dynamic stress changes from passing seismic waves. In the near-source region ( r ≤ 50 km for M ∼ 5 sources) both processes may be operating, and since both mechanisms are expected to raise earthquake rates, it is difficult to isolate them. We thus compare explosions with earthquakes because only earthquakes cause significant static stress changes. We find that large explosions at the Nevada Test Site do not trigger earthquakes at rates comparable to similar magnitude earthquakes. Surface waves are associated with regional and long-range dynamic triggering, but we note that surface waves with low enough frequency to penetrate to depths where most aftershocks of the 1992 M = 5.7 Little Skull Mountain main shock occurred (∼12 km) would not have developed significant amplitude within a 50-km radius. We therefore focus on the best candidate phases to cause local dynamic triggering, direct waves that pass through observed near-source aftershock clusters. We examine these phases, which arrived at the nearest (200–270 km) broadband station before the surface wave train and could thus be isolated for study. Direct comparison of spectral amplitudes of presurface wave arrivals shows that M ∼ 5 explosions and earthquakes deliver the same peak dynamic stresses into the near-source crust. We conclude that a static stress change model can readily explain observed aftershock patterns, whereas it is difficult to attribute near-source triggering to a dynamic process because of the dearth of aftershocks near large explosions.

Journal of Geophysical Research B: Solid Earth↗

USGS2011

Frictional strengths of talc-serpentine and talc-quartz mixtures

Talc is a constituent of faults in a variety of settings, and it may be an effective weakening agent depending on its abundance and distribution within a fault. We conducted frictional strength experiments under hydrothermal conditions to determine the effect of talc on the strengths of synthetic gouges of lizardite and antigorite serpentinites and of quartz. Small amounts of talc weaken serpentinite gouges substantially more than predicted by simple weight averaging. In comparison, mixtures of quartz and talc show a linear trend of strength reduction at talc concentrations 15 wt % and enhanced weakening at higher concentrations. All of the strength data are fit by a modified version of the Reuss mixing law that allows for the dominance of one mineral over the other. The difference in the behavior of serpentinite-talc and quartz-talc mixtures at low talc concentrations is a reflection of their different textures. Lizardite, antigorite, and talc all have platy habits, and displacement within gouges composed of these minerals is localized to narrow shears along which the platy grains have rotated into alignment with the shear surfaces. The shears in the mixed phyllosilicate gouges maximize the proportion of the weaker mineral within them. When mixed with a strong, rounded mineral such as quartz, some minimum concentration of talc is needed to form connected pathways that enhance strength reductions. The typical development of talc by the reaction of Si-rich fluids with serpentinite or dolomite would tend to localize its occurrence in a natural fault and result in enhanced weakening.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Global map of eolian features on Mars

Ten basic categories of eolian features on Mars were identified from a survey of Mariner 9 and Viking orbiter images. The ten features mapped are (1) light streaks (including frost streaks), (2) dark streaks, (3) sand sheets or splotches, (4) barchari dunes, (5) transverse dunes, (6) crescentic dunes, (7) anomalous dunes, (8) yardangs, (9) wind grooves, and (10) deflation pits. The features were mapped in groups, not as individual landforms, and recorded according to their geographic positions and orientations on maps of 1:12.5 million or 1:25 million scale. In the north polar region, light and dark streaks indicate winds from the west and northwest. Frost streaks show northeast and northwest winds. Barchan dunes show southwest and west winds. Transverse dunes show east and west winds. Local erosional features show winds from the northwest or southeast and northeast or southwest. In the middle and low northern latitudes, streaks show northeast winds; yardangs and deflation pits show mostly easterly and northeasterly winds. In the low southern latitudes, light streaks and dunes record northwest and northeast winds; dark streaks record southeast winds. In high southern latitudes, most streaks and dunes record southeast and east winds. In the south polar region, light and frost streaks record southwest and northwest winds, whereas dark streaks and transverse dunes show southeast winds. The patterns recorded for ephemeral features conform to global conditions of strong southern spring and summer (northern fall and winter) wind circulation. Erosional features in bedrock indicate long-term and perhaps ancient wind trends, whereas depositional features may record relatively more recent winds. Wind directions indicated by yardangs, pits, and other erosional features seldom correspond to those shown by streaks. These erosional features indicating winds of different directions than those that formed the streaks may have been carved when the effective regional winds on Mars were different due to relative hemisphere insolation differences, which are related to obliquity variations. Alternatively, many yardangs and pits may be largely controlled by structure and carved along surface structural patterns or features that channeled both ancient and present-day winds. Deflation pits in the mantled terrain may contain the best record of ancient wind trends. These pits seem to have stratigraphically related orientations. Different stratigraphic units composing the mantle contain wind erosion pits that possibly record large-scale wind direction changes because of long-term changes in the obliquity of the Mars.

Journal of Geophysical Research Solid Earth↗

USGS2012

Active high-resolution seismic tomography of compressional wave velocity and attenuation structure at Medicine Lake Volcano, northern California Cascade Range

We determine compressional wave velocity and attenuation structures for the upper crust beneath Medicine Lake volcano in northeast California using a high-resolution active source seismic tomography method. Medicine Lake volcano is a basalt through rhyolite shield volcano of the Cascade Range, lying east of the range axis. The Pg wave from eight explosive sources which has traveled upward through the target volume to a dense array of 140 seismographs provides 1- to 2-km resolution in the upper 5 to 7 km of the crust beneath the volcano. The experiment tests the hypothesis that Cascade Range volcanoes of this type are underlain only by small silicic magma chambers. We image a low-velocity low- Q region not larger than a few tens of cubic kilometers in volume beneath the summit caldera, supporting the hypothesis. A shallower high-velocity high-density feature, previously known to be present, is imaged for the first time in full plan view; it is east-west elongate, paralleling a topographic lineament between Medicine Lake volcano and Mount Shasta. This lineament is interpreted to be the result of an old crustal weakness now affecting the emplacement of magma, both on direct ascent from the lower crust and mantle and in migration from the shallow silicic chamber to summit vents. Differences between this high-velocity feature and the equivalent feature at Newbeny volcano, a volcano in central Oregon resembling Medicine Lake volcano, may partly explain the scarcity of surface hydrothermal features at Medicine Lake volcano. A major low-velocity low- Q feature beneath the southeast flank of the volcano, in an area with no Holocene vents, is interpreted as tephra, flows, and sediments from the volcano deeply ponded on the downthrown side of the Gillem fault, a normal fault mapped at the surface north of the volcano. A high- Q normal-velocity feature beneath the north rim of the summit caldera may be a small, possibly hot, subsolidus intrusion. A high-velocity low- Q region beneath the eastern caldera may be an area of boiling water between the magma chamber and the ponded east flank material. These structural data are useful both for understanding Cascade Range volcanism and for geothermal development in progress on the volcano.

Journal of Geophysical Research Solid Earth↗

USGS1999

Premonitory slip and tidal triggering of earthquakes

We have conducted a series of laboratory simulations of earthquakes using granite cylinders containing precut bare fault surfaces at 50 MPa confining pressure. Axial shortening rates between 10 −4 and 10 −6 mm/s were imposed to simulate tectonic loading. Average loading rate was then modulated by the addition of a small-amplitude sine wave to simulate periodic loading due to Earth tides or other sources. The period of the modulating signal ranged from 10 to 10,000 s. For each combination of amplitude and period of the modulating signal, multiple stick-slip events were recorded to determine the degree of correlation between the timing of simulated earthquakes and the imposed periodic loading function. Over the range of parameters studied, the degree of correlation of earthquakes was most sensitive to the amplitude of the periodic loading, with weaker dependence on the period of oscillations and the average loading rate. Accelerating premonitory slip was observed in these experiments and is a controlling factor in determining the conditions under which correlated events occur. In fact, some form of delayed failure is necessary to produce the observed correlations between simulated earthquake timing and characteristics of the periodic loading function. The transition from strongly correlated to weakly correlated model earthquake populations occurred when the amplitude of the periodic loading was approximately 0.05 to 0.1 MPa shear stress (0.03 to 0.06 MPa Coulomb failure function). Lower-amplitude oscillations produced progressively lower correlation levels. Correlations between static stress increases and earthquake aftershocks are found to degrade at similar stress levels. Typical stress variations due to Earth tides are only 0.001 to 0.004 MPa, so that the lack of correlation between Earth tides and earthquakes is also consistent with our findings. A simple extrapolation of our results suggests that approximately 1% of midcrustal earthquakes should be correlated with Earth tides. Triggered seismicity has been reported resulting from the passage of surface waves excited by the Landers earthquake. These transient waves had measured amplitudes in excess of 0.1 MPa at frequencies of 0.05 to 0.2 Hz in regions of notable seismicity increase. Similar stress oscillations in our laboratory experiments produced strongly correlated stick-slip events. We suggest that seemingly inconsistent natural observations of triggered seismicity and absence of tidal triggering indicate that failure is amplitude and frequency dependent. This is the expected result if, as in our laboratory experiments, the rheology of the Earth's crust permits delayed failure.

Journal of Geophysical Research B: Solid Earth↗

USGS1997

Three-dimensional upper crustal velocity structure beneath San Francisco Peninsula, California

This paper presents new seismic data from, and crustal models of the San Francisco Peninsula. In much of central California the San Andreas fault juxtaposes the Cretaceous granitic Salinian terrane on its west and the Late Mesozoic/Early Tertiary Franciscan Complex on its east. On San Francisco Peninsula, however, the present-day San Andreas fault is completely within a Franciscan terrane, and the Pilarcitos fault, located southwest of the San Andreas, marks the Salinian-Franciscan boundary. This circumstance has evoked two different explanations: either the Pilarcitos is a thrust fault that has pushed Franciscan rocks over Salinian rocks or the Pilarcitos is a transform fault that has accommodated significant right-lateral slip. In an effort to better resolve the subsurface structure of the peninsula faults, we established a temporary network of 31 seismographs arrayed across the San Andreas fault and the subparallel Pilarcitos fault at ∼1–2 km spacings. These instruments were deployed during the first 6 months of 1995 and recorded local earthquakes, air gun sources set off in San Francisco Bay, and explosive sources. Travel times from these sources were used to augment earthquake arrival times recorded by the Northern California Seismic Network and were inverted for three-dimensional velocity structure. Results show lateral velocity changes at depth (∼0.5–7 km) that correlate with downward vertical projections of the surface traces of the San Andreas and Pilarcitos faults. We thus interpret the faults as high-angle to vertical features (constrained to a 70°–110° dip range). From this we conclude that the Pilarcitos fault is probably an important strike-slip fault that accommodated much of the right-lateral plate boundary strain on the peninsula prior to the initiation of the modern-day San Andreas fault in this region sometime after about 3.0 m.y. ago.

Journal of Geophysical Research B: Solid Earth↗

USGS1997

The physical basis of glacier volume-area scaling

Ice volumes are known for only a few of the roughly 160,000 glaciers worldwide but are important components of many climate and sea level studies which require water flux estimates. A scaling analysis of the mass and momentum conservation equations shows that glacier volumes can be related by a power law to more easily observed glacier surface areas. The relationship requires four closure choices for the scaling behavior of glacier widths, slopes, side drag and mass balance. Reasonable closures predict a volume-area scaling exponent which is consistent with observations, giving a physical and practical basis for estimating ice volumes. Glacier volume is insensitive to perturbations in the mass balance scaling, but changes in average accumulation area ratios reflect significant changes in the scaling of both mass balance and ice volume.

Journal of Geophysical Research B: Solid Earth↗

USGS1999

Assessment of Mars Pathfinder landing site predictions

Remote sensing data at scales of kilometers and an Earth analog were used to accurately predict the characteristics of the Mars Pathfinder landing site at a scale of meters. The surface surrounding the Mars Pathfinder lander in Ares Vallis appears consistent with orbital interpretations, namely, that it would be a rocky plain composed of materials deposited by catastrophic floods. The surface and observed maximum clast size appears similar to predictions based on an analogous surface of the Ephrata Fan in the Channeled Scabland of Washington state. The elevation of the site measured by relatively small footprint delay-Doppler radar is within 100 m of that determined by two-way ranging and Doppler tracking of the spacecraft. The nearly equal elevations of the Mars Pathfinder and Viking Lander 1 sites allowed a prediction of the atmospheric conditions with altitude (pressure, temperature, and winds) that were well within the entry, descent, and landing design margins. High-resolution (~38 m/pixel) Viking Orbiter 1 images showed a sparsely cratered surface with small knobs with relatively low slopes, consistent with observations of these features from the lander. Measured rock abundance is within 10% of that expected from Viking orbiter thermal observations and models. The fractional area covered by large, potentially hazardous rocks observed is similar to that estimated from model rock distributions based on data from the Viking landing sites, Earth analog sites, and total rock abundance. The bulk and fine-component thermal inertias measured from orbit are similar to those calculated from the observed rock size-frequency distribution. A simple radar echo model based on the reflectivity of the soil (estimated from its bulk density), and the measured fraction of area covered by rocks was used to approximate the quasi-specular and diffuse components of the Earth-based radar echos. Color and albedo orbiter data were used to predict the relatively dust free or unweathered surface around the Pathfinder lander compared to the Viking landing sites. Comparisons with the experiences of selecting the Viking landing sites demonstrate the enormous benefit the Viking data and its analyses and models had on the successful predictions of the Pathfinder site. The Pathfinder experience demonstrates that, in certain locations, geologic processes observed in orbiter data can be used to infer surface characteristics where those processes dominate over other processes affecting the Martian surface layer. Copyright 1999 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗

USGS1997

The roles of time and displacement in velocity-dependent volumetric strain of fault zones

The relationship between measured friction μ A and volumetric strain during frictional sliding was determined using a rate and state variable dependent friction constitutive equation, a common work balance relating friction and volume change, and two types of experimental faults: initially bare surfaces of Westerly granite and rock surfaces separated by a 1 mm layer of less than 90 μm Westerly granite gouge. The constitutive equation is the sum of a constant term representing the nominal resistance to sliding and two smaller terms: a rate dependent term representing the shear viscosity of the fault surface (direct effect), and a term which represents variations in the area of contact (evolution effect). The work balance relationship requires that μ A differs from the frictional resistance that leads to shear heating by the derivative of fault normal displacement with respect shear displacement, d δ n / d δ s . An implication of this relationship is that the rate dependence of d δ n / d δ s contributes to the rate dependence of μ A . Experiments show changes in sliding velocity lead to changes in both fault strength and volume. Analysis of data with the rate and state equations combined with the work balance relationship preclude the conventional interpretation of the direct effect in the rate and state variable constitutive equations. Consideration of a model bare surface fault consisting of an undeformable indentor sliding on a deformable surface reveals a serious flaw in the work balance relationship if volume change is time-dependent. For the model, at zero slip rate indentation creep under the normal load leads to time-dependent strengthening of the fault surface but, according to the work balance relationship, no work is done because compaction or dilatancy can only be induced by shearing. Additional tests on initially bare surfaces and gouges show that fault normal strain in experiments is time-dependent, consistent with the model. This time-dependent fault normal strain, which is not accounted for in the work balance relationship, explains the inconsistency between the constitutive equations and the work balance. For initially bare surface faults, all rate dependence of volume change is due to time dependence. Similar results are found for gouge. We conclude that μ A reflects the frictional resistance that results in shear heating, and no correction needs to be made for the volume changes. The result that time-dependent volume changes do not contribute to μ A is a general result and extends beyond these experiments, the simple indentor model and particular constitutive equations used to illustrate the principle.

Journal of Geophysical Research B: Solid Earth↗

USGS1998

Ductile shear zones beneath strike-slip faults: Implications for the thermomechanics of the San Andreas fault zone

We have carried out two-dimensional (2-D) numerical experiments on the bulk flow of a layer of fluid that is driven in a strike-slip sense by constant velocities applied at its boundaries. The fluid has the (linearized) conventional rheology assumed to apply to lower crust/upper mantle rocks. The temperature dependence of the effective viscosity of the fluid and the shear heating that accompanies deformation have been incorporated into the calculations, as has thermal conduction in an overlying crustal layer. Two end-member boundary conditions have been considered, corresponding to a strong upper crust driving a weaker ductile substrate and a strong ductile layer driving a passive, weak crust. In many cases of practical interest, shear heating is concentrated close to the axial plane of the shear zone for either boundary condition. For these cases, the resulting steady state temperature field is well approximated by a cylindrical heat source embedded in a conductive half-space at a depth corresponding to the top of the fluid layer. This approximation, along with the application of a theoretical result for one-dimensional shear zones, permits us to obtain simple analytical approximations to the thermal effects of 2-D ductile shear zones for a range of assumed rheologies and crustal geotherms, making complex numerical calculations unnecessary. Results are compared with observable effects on heat flux near the San Andreas fault using constraints on the slip distribution across the entire fault system. Ductile shearing in the lower crust or upper mantle can explain the observed increase in surface heat flux southeast of the Mendocino triple junction and match the amplitude of the regional heat flux anomaly in the California Coast Ranges. Because ductile dissipation depends only weakly on slip rate, faults moving only a few millimeters per year can be important heat sources, and the superposition of effects of localized ductile shearing on both currently active and now inactive strands of the San Andreas system can explain the breadth of the heat flux anomaly across central California.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Geology of the Mount St. Helens area: Record of discontinuous volcanic and plutonic activity in the Cascade Arc of southern Washington

The Quaternary edifice of Mount St. Helens volcano was built upon a deeply eroded terrane of gently folded and altered volcanic and plutonic rocks that represent the core of the Tertiary Cascade magmatic arc. These rocks constitute an east dipping homoclinal sequence, several kilometers thick, of subaerially erupted mafic to silicic flows and volcaniclastic strata; K‐Ar ages from this section range from about 28 to 23 Ma (late Oligocene and earliest Miocene), which corresponds to an apparent lull in Cascade volcanism to the north of Mount Rainier. Volcanism was essentially continuous during this period of time, and neither a well‐defined base nor top of the section is exposed within the mapped area. Basalt and basaltic andesite dominate the lower part of the mid‐Tertiary section, whereas andesitic and dacitic rocks comprise most of the upper part. This section was intruded by numerous mafic to silicic dikes, sills, and irregular plutonic bodies, most no more than a few million years younger than their host rocks, and subjected to pervasive burial metamorphism and widespread hydrothermal alteration. Large areas of hornfelsed rock surrounding even relatively small intrusions indicate that the proportion of plutonic rock becomes significantly greater at shallow depth beneath the existing erosion surface. A large granitic pluton intruded the mid‐Tertiary section north of Spirit Lake at about 21 Ma. The Earl porphyry copper deposit occurs within the pluton but appears too young (17 Ma) to be genetically related to it. In contrast to the rather continuous and voluminous Oligocene to early Miocene activity, volcanism since then in the Mount St. Helens area has been localized and Volumetrically minor. Products of three younger eruptive periods have been recognized: a sequence of 15 m.y. old pyroxene andesite flows resting unconformably on mid‐Tertiary strata south of Mount St. Helens, widespread shallow dikes and sills of pyroxene andesite between 10 and 8 m.y. old, and compositionally diverse rocks erupted during the past 3 m.y. The Quaternary lavas are more potassic than the Tertiary lavas and typically contain phenocrysts of hornblende and biotite, which are absent from the older rocks. A number of Tertiary structures define a broad NNE trending zone that may reflect a deep‐seated lithospheric flaw that has controlled the locus of Cascade magmatism in southern Washington for the past 25 m.y. Mount St. Helens lies within this zone at the intersection of the NNW striking St. Helens seismic zone (SHZ) and an ENE trending alignment of Pleistocene silicic plug‐domes. No surface breakage has been detected along the SHZ, which is apparently very young. The linear zone of silicic vents is probably controlled by a fault that has been interpreted from seismic records to occur directly beneath the volcano. This zone parallels the directions of regional maximum horizontal compressive stress and North America/Juan de Fuca plate convergence. Mount St. Helens is an example of a low‐volume tectonically controlled magmatic system in an early stage of development.

Washington↗

USGS2018

Quantifying uncertainty in cumulative surface slip along the Cucamonga Fault, a crustal thrust fault in southern California

Studies of historic earthquake ground surface ruptures show that displacements along strike are spatially variable. As a result, latest Quaternary slip rates developed from a spatially restricted set of cumulative displacement measurements may not accurately represent fault velocity. Here we examine the uncertainties associated with slip on the Cucamonga Fault, which is part of a network of faults that have generated damaging historical earthquakes in and around Los Angeles, California. Numerous scarps along its ~25‐km length are well expressed on alluvial fans. We make 310 measurements of vertical separation across the scarps using lidar data. We show that the dispersion of the vertical separations cannot be explained by our best estimates of analytical uncertainties alone. Additional epistemic uncertainties are required. We find that the magnitude of the required epistemic uncertainty is typically larger than analytical uncertainty by a factor of 3 and typically about 22% of the maximum vertical separation. These relationships appear to hold at several spatial scales. We examine three potential sources of epistemic uncertainty and find that none among surface age uncertainty, fault dip, and anthropogenic landscape alteration is likely sufficient to explain the overdispersion of the data, which suggests differences in cumulative strain along the strike of the fault. We calculate a range of dip‐slip rates between 0.4 and 2.6 mm/year. In light of our results, we suggest that future thrust‐fault slip‐rate studies adopt an epistemic uncertainty of 22% of the maximum value for vertical separation measurements, unless there are sufficient data to demonstrate otherwise.

California↗

USGS2012

Experimental deformation of topaz crystals: Possible embrittlement by intracrystalline water

Crystallographically oriented single‐crystal prisms of gem quality topaz (composition AlSiO (OH− F) where x = 0.04 ± 0.01) were deformed at a confining pressure of 1.50 GPa, a temperature of 800°C, and a strain rate of 2×10 s. Under nearly identical conditions, all crystals of anhydrous rock‐forming minerals that have been tested to date, such as olivine, quartz, feldspars, pyroxenes, and refractory oxides, deform plastically; in contrast, our topaz crystals failed by brittle fracture regardless of the orientation of the compression direction. No optical evidence for plastic deformation was detected. Another suite of experiments with compression perpendicular to the (001) cleavage at = 100°–950°C and a strain rate of 2×10 s displayed two regimes of behavior: (1) at >400°C, fracture strength was independent of temperature, and fracture occurred on one or two surfaces parallel to {103}; (2) at <400°C, the fracture strength increased rapidly with decreasing temperature, no macroscopic stress drop was observed, and many closely spaced conjugate fractures formed on (103) and (103). The anomalous brittleness of topaz compared to anhydrous silicate and oxide crystals indicates that intracrystalline “water” plays a role in the embrittlement. We suggest that water within the topaz crystals promotes fracture in ways similar to the mechanisms of slow crack growth aided by environmental moisture.

Journal of Geophysical Research B: Solid Earth↗

USGS2009

Observations of periglacial landforms in Utopia Planitia with the High Resolution Imaging Science Experiment (HiRISE)

The region of western Utopia Planitia (80-105 degreesE, 40-55 degrees N) displays several types of landforms similar to Earth periglacial features, including scallop-shaped depressions and networks of polygonal terrains. The scalloped depressions have been proposed to originate from thermokarstic processes such as sublimation and/or melting of near-surface ground ice. Using HiRISE imagery, we characterize these depressions and several associated, distinct polygon networks in unprecedented morphologic and topographic detail and investigate support for an ice-based degradation process. The scalloped depressions and interior polygons and ridges are found to evolve together, mainly influenced by sublimation, local proximity of ground ice to the surface, and obliquity variations.

Journal of Geophysical Research↗

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↗

USGS2015

Changes in seismic velocity during the first 14 months of the 2004–2008 eruption of Mount St. Helens, Washington

Mount St. Helens began erupting in late 2004 following an 18 year quiescence. Swarms of repeating earthquakes accompanied the extrusion of a mostly solid dacite dome over the next 4 years. In some cases the waveforms from these earthquakes evolved slowly, likely reflecting changes in the properties of the volcano that affect seismic wave propagation. We use coda-wave interferometry to quantify small changes in seismic velocity structure (usually <1%) between two similar earthquakes and employed waveforms from several hundred families of repeating earthquakes together to create a continuous function of velocity change observed at permanent stations operated within 20 km of the volcano. The high rate of earthquakes allowed tracking of velocity changes on an hourly time scale. Changes in velocity were largest near the newly extruding dome and likely related to shallow deformation as magma first worked its way to the surface. We found strong correlation between velocity changes and the inverse of real-time seismic amplitude measurements during the first 3 weeks of activity, suggesting that fluctuations of pressure in the shallow subsurface may have driven both seismicity and velocity changes. Velocity changes during the remainder of the eruption likely result from a complex interplay of multiple effects and are not well explained by any single factor alone, highlighting the need for complementary geophysical data when interpreting velocity changes.

Journal of Geophysical Research B: Solid Earth↗

USGS2005

The U-tube: A novel system for acquiring borehole fluid samples from a deep geologic CO2 sequestration experiment

A novel system has been deployed to obtain geochemical samples of water and gas, at in situ pressure, during a geologic CO 2 sequestration experiment conducted in the Frio brine aquifer in Liberty County, Texas. Project goals required high-frequency recovery of representative and uncontaminated aliquots of a rapidly changing two-phase fluid (supercritical CO 2 and brine) fluid from 1.5 km depth. The data sets collected, using both the liquid and gas portions of the downhole samples, provide insights into the coupled hydrogeochemical issues affecting CO 2 sequestration in brine-filled formations. While the basic premise underlying the U-tube sampler is not new, the system is unique because careful consideration was given to the processing of the recovered two-phase fluids. In particular, strain gauges mounted beneath the high-pressure surface sample cylinders measured the ratio of recovered brine to supercritical CO 2 . A quadrupole mass spectrometer provided real-time gas analysis for perfluorocarbon and noble gas tracers that were injected along with the CO 2 . The U-tube successfully acquired frequent samples, facilitating accurate delineation of the arrival of the CO 2 plume, and on-site analysis revealed rapid changes in geochemical conditions.

Journal of Geophysical Research Solid Earth↗

USGS1996

Analysis of fractures intersecting Kahi Puka Well 1 and its relation to the growth of the island of Hawaii

As part of the Hawaii Scientific Drilling Project, Kahi Puka Well 1 penetrated about 275 m of Mauna Loa basalts overlying a sequence of Mauna Kea flow units as it was drilled and cored to a total depth of 1053 m below land surface. A borehole televiewer (BHTV) was run in most of the well in successive stages prior to casing in order to obtain magnetically oriented acoustic images of the borehole wall. A total of 283 individual fractures were identified from this log and characterized in terms of strike and dip. These data are divided into three vertical sections based upon age and volcanic source, and lower hemisphere stereographic plots identify two predominant, subparallel fracture subsets common to each section. Assuming that most of the steeply dipping fractures observed in the BHTV log are tensile features generated within basalt flows during deposition and cooling, this fracture information can be combined with models of the evolution of the island of Hawaii to investigate the depositional history of these Mauna Loa and Mauna Kea basalts over the past 400 kyr. The directions of high-angle fractures appear to be generally parallel to topography or to the coastline at the time of deposition, as is supported by surface mapping of modern flows. Consequently, an overall counterclockwise rotation of about 75?? in the strike of these fractures from the bottom to the top of the well represents a systematic change in depositional slope direction over time. We attribute the observed rotation in the orientations of the two predominant fracture subsets over the past 400 kyr to changes in the configurations of volcanic sources during shield building and to the structural interference of adjacent volcanoes that produces shifts in topographic patterns.

Journal of Geophysical Research B: Solid Earth↗