SEARCH · Geology Reports

Results for “Journal of Geophysical Research”

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 91 records · Page 5Linked to original sources

USGS2014

Temporal variability of carbon and nutrient burial, sediment accretion, and mass accumulation over the past century in a carbonate platform mangrove forest of the Florida Everglades.

The objective of this research was to measure temporal variability in accretion and mass sedimentation rates (including organic carbon (OC), total nitrogen (TN), and total phosphorous (TP)) from the past century in a mangrove forest on the Shark River in Everglades National Park, USA. The 210 Pb Constant Rate of Supply model was applied to six soil cores to calculate annual rates over the most recent 10, 50, and 100 year time spans. Our results show that rates integrated over longer timeframes are lower than those for shorter, recent periods of observation. Additionally, the substantial spatial variability between cores over the 10 year period is diminished over the 100 year record, raising two important implications. First, a multiple-decade assessment of soil accretion and OC burial provides a more conservative estimate and is likely to be most relevant for forecasting these rates relative to long-term processes of sea level rise and climate change mitigation. Second, a small number of sampling locations are better able to account for spatial variability over the longer periods than for the shorter periods. The site average 100 year OC burial rate, 123 ± 19 (standard deviation) g m -2 yr -1 , is low compared with global mangrove values. High TN and TP burial rates in recent decades may lead to increased soil carbon remineralization, contributing to the low carbon burial rates. Finally, the strong correlation between OC burial and accretion across this site signals the substantial contribution of OC to soil building in addition to the ecosystem service of CO 2 sequestration.

Florida↗

USGS2019

Identifying the molecular signatures of agricultural expansion in Amazonian headwater streams

Agricultural impacts on aquatic ecosystems are well-studied, however, most research has focused on temperate regions, whereas the forefront of agricultural expansion is currently in the tropics. At the vanguard of this growth is the boundary between the Amazon and Cerrado biomes in Brazil, driven primarily by expansion of soybean and corn croplands. Here we examine the impacts of cropland expansion on receiving lowland Amazon Basin headwater streams in terms of dissolved organic carbon (DOC) concentration and dissolved organic matter (DOM) composition via ultra-high resolution mass spectrometry. Streams draining croplands had lower DOC concentrations and DOM molecular signatures enriched in N- and S-containing formula in comparison to forested streams. Cropland streams were also enriched in aliphatic, peptide-like and highly unsaturated and phenolic (low O/C) compound categories in comparison to forest streams (enriched in polyphenolics, condensed aromatics and highly unsaturated and phenolic (high O/C) compound categories) indicative of the shifting of sources from organic-rich surface soils and litter layers to autochthonous and more microbial biomass. Distinct molecular assemblages were strongly correlated with cropland and forest catchments, highlighting headwater streams as sentinels for detecting change. On investigation of unique molecular formulae present in only cropland sites, four cropland markers provided the ability to track agricultural impacts in the region. Overall, these patterns indicate reduced organic matter inputs in croplands and greater microbial degradation at these sites leading to declining DOC concentrations, and DOM of more microbial character in receiving streams that is more biolabile, with clear ramifications for downstream ecology and biogeochemical cycles.

Amazon River basin, Tanguro Ranch, Xinga River↗

USGS2012

Characterizing mercury concentrations and fluxes in a Coastal Plain watershed: Insights from dynamic modeling and data

Mercury (Hg) is one of the leading water quality concerns in surface waters of the United States. Although watershed-scale Hg cycling research has increased in the past two decades, advances in modeling watershed Hg processes in diverse physiographic regions, spatial scales, and land cover types are needed. The goal of this study was to assess Hg cycling in a Coastal Plain system using concentrations and fluxes estimated by multiple watershed-scale models with distinct mathematical frameworks reflecting different system dynamics. We simulated total mercury (Hg T , the sum of filtered and particulate forms) concentrations and fluxes from a Coastal Plain watershed (McTier Creek) using three watershed Hg models and an empirical load model. Model output was compared with observed in-stream Hg T . We found that shallow subsurface flow is a potentially important transport mechanism of particulate Hg T during periods when connectivity between the uplands and surface waters is maximized. Other processes (e.g., stream bank erosion, sediment re-suspension) may increase particulate Hg T in the water column. Simulations and data suggest that variable source area (VSA) flow and lack of rainfall interactions with surface soil horizons result in increased dissolved Hg T concentrations unrelated to DOC mobilization following precipitation events. Although flushing of DOC-Hg T complexes from surface soils can also occur during this period, DOC-complexed Hg T becomes more important during base flow conditions. TOPLOAD simulations highlight saturated subsurface flow as a primary driver of daily Hg T loadings, but shallow subsurface flow is important for Hg T loads during high-flow events. Results suggest limited seasonal trends in Hg T dynamics.

South Carolina↗

USGS2003

Climate change and Arctic ecosystems: 1. Vegetation changes north of 55°N between the last glacial maximum, mid-Holocene, and present

A unified scheme to assign pollen samples to vegetation types was used to reconstruct vegetation patterns north of 55??N at the last glacial maximum (LGM) and mid-Holocene (6000 years B.P.). The pollen data set assembled for this purpose represents a comprehensive compilation based on the work of many projects and research groups. Five tundra types (cushion forb tundra, graminoid and forb tundra, prostrate dwarf-shrub tundra, erect dwarf-shrub tundra, and low- and high-shrub tundra) were distinguished and mapped on the basis of modern pollen surface samples. The tundra-forest boundary and the distributions of boreal and temperate forest types today were realistically reconstructed. During the mid-Holocene the tundra-forest boundary was north of its present position in some regions, but the pattern of this shift was strongly asymmetrical around the pole, with the largest northward shift in central Siberia (???200 km), little change in Beringia, and a southward shift in Keewatin and Labrador (???200 km). Low- and high-shrub tundra extended farther north than today. At the LGM, forests were absent from high latitudes. Graminoid and forb tundra abutted on temperate steppe in northwestern Eurasia while prostrate dwarf-shrub, erect dwarf-shrub, and graminoid and forb tundra formed a mosaic in Beringia. Graminoid and forb tundra is restricted today and does not form a large continuous biome, but the pollen data show that it was far more extensive at the LGM, while low- and high-shrub tundra were greatly reduced, illustrating the potential for climate change to dramatically alter the relative areas occupied by different vegetation types.

Journal of Geophysical Research D: Atmospheres↗

USGS2002

Draped aeromagnetic survey in Transantarctic Mountains over the area of the Butcher Ridge igneous complex showing extent of underlying mafic intrusion

A draped aeromagnetic survey over the area surrounding the Butcher Ridge igneous complex (BRIC), Transantarctic Mountains, was acquired in 1997–1998 as part of a larger Transantarctic Mountains Aerogeophysical Research Activity survey. The BRIC is a sill-like hypoabyssal intrusion ranging in composition from tholeiitic basalt to rhyolite. An 40 Ar/ 39 Ar age of 174 Ma and the chemical character of the basaltic rocks show the BRIC to be part of the widespread Jurassic Ferrar suite of continental tholeiitic rocks, that extends for 3500 km across Antarctica. The aeromagnetic survey shows a horseshoe-shaped pattern of anomalies reaching amplitudes as great as 1900 nT generally associated with the bedrock topography where it is exposed. It is apparent that the high-amplitude anomaly pattern is more extensive than the 10-km-long exposed outcrop, first crossed by a single 1990 aeromagnetic profile. The highest-amplitude anomalies appear south of the profile acquired in 1990 and extend out of the survey area. The new aeromagnetic data allow determination of the extent of the interpreted Butcher mafic(?) intrusion beneath exposures of Beacon sedimentary rock and ice in the area covered, as well as beneath the small BRIC exposure. The magnetic anomalies show a minimum area of 3000 km 2 , a much greater extent than previously inferred. Magnetic models indicate a minimum thickness of ∼1–2 km for a horizontal intrusion. However, nonunique models with magnetic layers decreasing in apparent susceptibility with depth are consistent with of a 4- to 8-km-thick layered intrusion. These magnetic models indicate progressively deeper erosion of the interpreted mafic-layered body from the south to north. The erosion has removed more magnetic upper layers that mask the magnetic effects of the lower less magnetic layers. The probable minimum volume of the intrusion in the area of the survey is ∼6000 km 3 . An alternate, but less likely, interpretation of a series of dikes can also fit the observed magnetic anomalies.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Shallow structure and deformation along the San Andreas Fault in Cholame Valley, California, based on high-resolution reflection profiling

The mapped active traces of the San Andreas fault are separated by a 1-km-wide right-stepping offset in Cholame Valley. The geometry of this offset, defined in other strike-slip systems as a releasing bend or a dilational jog, has resulted in the formation of a pullapart basin. Various researchers have inferred that this offset served as a rupture terminus for earthquakes on both strands of the San Andreas fault (1966 Parkfield and 1857 Fort Tejon); thus, this en echelon offset may represent a barrier to the propagation of rupture between two segments of the San Andreas fault. We collected 18 km of high-resolution seismic reflection data specifically designed to image the San Andreas fault zone in the shallow crust surrounding this offset. Four short profiles (≤ 3.3 km long) of Mini-Sosie reflection data (1 s of two-way travel time ≈ 1.5 km deep) were collected perpendicular to the San Andreas fault; three of these profiles were tied by a 7-km-long profile that trended northwest through Cholame Valley, subparallel to the San Andreas fault. A zone of incoherent energy, narrow at the surface but widening with depth, underlies the mapped active traces of the San Andreas fault and abruptly terminates shallow reflections on both sides. The reflection profiles and available well data indicate that west of the mapped active traces of the San Andreas fault the shallow subsurface structure of the crust consists of thin (≤400 m thick), offset packages of reflections, laterally coherent on the scale of tens of meters, overlying deformed clastic sedimentary rocks. East of the San Andreas fault, the structure of the shallow crust in southern Cholame Valley is characterized by thick packages of reflections, laterally coherent on the scale of kilometers, overlying the Franciscan complex. All of the strata east of the fault (within Cholame Valley) dip toward the San Andreas fault and the offset, into an approximately 1-km-deep sedimentary basin abutting the south strand of San Andreas fault. The offset in Cholame Valley is characterized by a gentle downwarping of sediments into the offset, the presence of many small faults and discontinuous reflections (consistent with areally distributed shear) between the traces of the main fault, localized subsidence abutting the main strike-slip fault, the formation of a basin, near the offset, that is about as deep as the jog is wide, and the southward propagation of the deformation associated with the offset. Strain field modeling based on simple geometries of the San Andreas and associated faults successfully predicts the general features of the observed topography and subsurface structure of southern Cholame Valley, including subsidence and basin formation near the offset, a discontinuous San Andreas fault plane, and at least one fault in southeastern Cholame Valley.

Journal of Geophysical Research↗

USGS2013

Inferring fault rheology from low-frequency earthquakes on the San Andreas

Families of recurring low-frequency earthquakes (LFEs) within nonvolcanic tremor (NVT) on the San Andreas fault in central California show strong sensitivity to shear stress induced by the daily tidal cycle. LFEs occur at all levels of the tidal shear stress and are in phase with the very small, ~400 Pa, stress amplitude. To quantitatively explain the correlation, we use a model from the existing literature that assumes the LFE sources are small, persistent regions that repeatedly fail during shear of a much larger scale, otherwise aseismically creeping fault zone. The LFE source patches see tectonic loading, creep of the surrounding fault which may be modulated by the tidal stress, and direct tidal loading. If the patches are small relative to the surrounding creeping fault then the stressing is dominated by fault creep, and if patch failure occurs at a threshold stress, then the resulting seismicity rate is proportional to the fault creep rate or fault zone strain rate. Using the seismicity rate as a proxy for strain rate and the tidal shear stress, we fit the data with possible fault rheologies that produce creep in laboratory experiments at temperatures of 400 to 600°C appropriate for the LFE source depth. The rheological properties of rock-forming minerals for dislocation creep and dislocation glide are not consistent with the observed fault creep because strong correlation between small stress perturbations and strain rate requires perturbation on the order of the ambient stress. The observed tidal modulation restricts ambient stress to be at most a few kilopascal, much lower than rock strength. A purely rate dependent friction is consistent with the observations only if the product of the friction rate dependence and effective normal stress is ~ 0.5 kPa. Extrapolating the friction rate strengthening dependence of phyllosilicates (talc) to depth would require the effective normal stress to be ~50 kPa, implying pore pressure is lithostatic. If the LFE source is on the order of tens of meters, as required by the model, rate-weakening friction rate dependence (e.g., olivine) at 400 to 600°C requires that the minimum effective pressure at the LFE source is ~ 2.5 MPa.

California↗

USGS2013

Triggered tremor sweet spots in Alaska

To better understand what controls fault slip along plate boundaries, we have exploited the abundance of seismic and geodetic data available from the richly varied tectonic environments composing Alaska. A search for tremor triggered by 11 large earthquakes throughout all of seismically monitored Alaska reveals two tremor “sweet spots”—regions where large-amplitude seismic waves repeatedly triggered tremor between 2006 and 2012. The two sweet spots locate in very different tectonic environments—one just trenchward and between the Aleutian islands of Unalaska and Akutan and the other in central mainland Alaska. The Unalaska/Akutan spot corroborates previous evidence that the region is ripe for tremor, perhaps because it is located where plate-interface frictional properties transition between stick-slip and stably sliding in both the dip direction and laterally. The mainland sweet spot coincides with a region of complex and uncertain plate interactions, and where no slow slip events or major crustal faults have been noted previously. Analyses showed that larger triggering wave amplitudes, and perhaps lower frequencies (<~0.03 Hz), may enhance the probability of triggering tremor. However, neither the maximum amplitude in the time domain or in a particular frequency band, nor the geometric relationship of the wavefield to the tremor source faults alone ensures a high probability of triggering. Triggered tremor at the two sweet spots also does not occur during slow slip events visually detectable in GPS data, although slow slip below the detection threshold may have facilitated tremor triggering.

Alaska↗

USGS2013

Aftershocks of the 2010 M w 7.2 El Mayor-Cucapah earthquake revealcomplex faulting in the Yuha Desert, California

We detect and precisely locate over 9500 aftershocks that occurred in the Yuha Desert region during a 2 month period following the 4 April 2010 M w 7.2 El Mayor-Cucapah (EMC) earthquake. Events are relocated using a series of absolute and relative relocation procedures that include Hypoinverse, Velest, and hypoDD. Location errors are reduced to ~40 m horizontally and ~120 m vertically.Aftershock locations reveal a complex pattern of faulting with en echelon fault segments trending toward the northwest, approximately parallel to the North American-Pacific plate boundary and en echelon, conjugate features trending to the northeast. The relocated seismicity is highly correlated with published surface mapping of faults that experienced triggered surface slip in response to the EMC main shock. Aftershocks occurred between 2 km and 11 km depths, consistent with previous studies of seismogenic thickness in the region. Three-dimensional analysis reveals individual and intersecting fault planes that are limited in their along-strike length. These fault planes remain distinct structures at depth, indicative of conjugate faulting, and do not appear to coalesce onto a throughgoing fault segment. We observe a complex spatiotemporal migration of aftershocks, with seismicity that jumps between individual fault segments that are active for only a few days to weeks. Aftershock rates are roughly consistent with the expected earthquake production rates of Dieterich (1994). The conjugate pattern of faulting and nonuniform aftershock migration patterns suggest that strain in the Yuha Desert is being accommodated in a complex manner.

California↗

USGS2014

Maximum magnitude earthquakes induced by fluid injection

Analysis of numerous case histories of earthquake sequences induced by fluid injection at depth reveals that the maximum magnitude appears to be limited according to the total volume of fluid injected. Similarly, the maximum seismic moment seems to have an upper bound proportional to the total volume of injected fluid. Activities involving fluid injection include (1) hydraulic fracturing of shale formations or coal seams to extract gas and oil, (2) disposal of wastewater from these gas and oil activities by injection into deep aquifers, and (3) the development of enhanced geothermal systems by injecting water into hot, low-permeability rock. Of these three operations, wastewater disposal is observed to be associated with the largest earthquakes, with maximum magnitudes sometimes exceeding 5. To estimate the maximum earthquake that could be induced by a given fluid injection project, the rock mass is assumed to be fully saturated, brittle, to respond to injection with a sequence of earthquakes localized to the region weakened by the pore pressure increase of the injection operation and to have a Gutenberg-Richter magnitude distribution with a b value of 1. If these assumptions correctly describe the circumstances of the largest earthquake, then the maximum seismic moment is limited to the volume of injected liquid times the modulus of rigidity. Observations from the available case histories of earthquakes induced by fluid injection are consistent with this bound on seismic moment. In view of the uncertainties in this analysis, however, this should not be regarded as an absolute physical limit.

Journal of Geophysical Research B: Solid Earth↗

USGS2014

Observations of static Coulomb stress triggering of the November 2011 M5.7 Oklahoma earthquake sequence

In November 2011, a M 5.0 earthquake occurred less than a day before a M 5.7 earthquake near Prague, Oklahoma, which may have promoted failure of the mainshock and thousands of aftershocks along the Wilzetta fault, including a M 5.0 aftershock. The M 5.0 foreshock occurred in close proximity to active fluid injection wells; fluid injection can cause a buildup of pore fluid pressure, decrease the fault strength, and may induce earthquakes. Keranen et al . [ 2013 ] links the M 5.0 foreshock with fluid injection, but the relationship between the foreshock and successive events has not been investigated. Here we examine the role of coseismic Coulomb stress transfer on earthquakes that follow the M 5.0 foreshock, including the M 5.7 mainshock. We resolve the static Coulomb stress change onto the focal mechanism nodal plane that is most consistent with the rupture geometry of the three M ≥ 5.0 earthquakes, as well as specified receiver fault planes that reflect the regional stress orientation. We find that Coulomb stress is increased, e.g., fault failure is promoted, on the nodal planes of ~60% of the events that have focal mechanism solutions, and more specifically, that the M 5.0 foreshock promoted failure on the rupture plane of the M 5.7 mainshock. We test our results over a range of effective coefficient of friction values. Hence, we argue that the M 5.0 foreshock, induced by fluid injection, potentially triggered a cascading failure of earthquakes along the complex Wilzetta fault system.

Oklahoma↗

USGS2014

Real-time inversions for finite fault slip models and rupture geometry based on high-rate GPS data

We present an inversion strategy capable of using real-time high-rate GPS data to simultaneously solve for a distributed slip model and fault geometry in real time as a rupture unfolds. We employ Bayesian inference to find the optimal fault geometry and the distribution of possible slip models for that geometry using a simple analytical solution. By adopting an analytical Bayesian approach, we can solve this complex inversion problem (including calculating the uncertainties on our results) in real time. Furthermore, since the joint inversion for distributed slip and fault geometry can be computed in real time, the time required to obtain a source model of the earthquake does not depend on the computational cost. Instead, the time required is controlled by the duration of the rupture and the time required for information to propagate from the source to the receivers. We apply our modeling approach, called Bayesian Evidence-based Fault Orientation and Real-time Earthquake Slip, to the 2011 Tohoku-oki earthquake, 2003 Tokachi-oki earthquake, and a simulated Hayward fault earthquake. In all three cases, the inversion recovers the magnitude, spatial distribution of slip, and fault geometry in real time. Since our inversion relies on static offsets estimated from real-time high-rate GPS data, we also present performance tests of various approaches to estimating quasi-static offsets in real time. We find that the raw high-rate time series are the best data to use for determining the moment magnitude of the event, but slightly smoothing the raw time series helps stabilize the inversion for fault geometry.

Journal of Geophysical Research↗

USGS2014

Crustal earthquake triggering by pre-historic great earthquakes on subduction zone thrusts

Triggering of earthquakes on upper plate faults during and shortly after recent great ( M> 8.0) subduction thrust earthquakes raises concerns about earthquake triggering following Cascadia subduction zone earthquakes. Of particular regard to Cascadia was the previously noted, but only qualitatively identified, clustering of M> ~6.5 crustal earthquakes in the Puget Sound region between about 1200–900 cal yr B.P. and the possibility that this was triggered by a great Cascadia thrust subduction thrust earthquake, and therefore portends future such clusters. We confirm quantitatively the extraordinary nature of the Puget Sound region crustal earthquake clustering between 1200–900 cal yr B.P., at least over the last 16,000. We conclude that this cluster was not triggered by the penultimate, and possibly full-margin, great Cascadia subduction thrust earthquake. However, we also show that the paleoseismic record for Cascadia is consistent with conclusions of our companion study of the global modern record outside Cascadia, that M> 8.6 subduction thrust events have a high probability of triggering at least one or more M> ~6.5 crustal earthquakes.

Cascadia;Puget Sound↗

USGS2014

Volcanic tremor masks its seismogenic source: Results from a study of noneruptive tremor recorded at Mount St. Helens, Washington

On 2 October 2004, a significant noneruptive tremor episode occurred during the buildup to the 2004–2008 eruption of Mount St. Helens (Washington). This episode was remarkable both because no explosion followed, and because seismicity abruptly stopped following the episode. This sequence motivated us to consider a model for volcanic tremor that does not involve energetic gas release from magma but does involve movement of conduit magma through extension on its way toward the surface. We found that the tremor signal was composed entirely of Love and Rayleigh waves and that its spectral bandwidth increased and decreased with signal amplitude, with broader bandwidth signals containing both higher and lower frequencies. Our modeling results demonstrate that the forces giving rise to this tremor were largely normal to conduit walls, generating hybrid head waves along conduit walls that are coupled to internally reflected waves. Together these form a crucial part of conduit resonance, giving tremor wavefields that are largely a function of waveguide geometry and velocity. We find that the mechanism of tremor generation fundamentally masks the nature of the seismogenic source giving rise to resonance. Thus multiple models can be invoked to explain volcanic tremor, requiring that information from other sources (such as visual observations, geodesy, geology, and gas geochemistry) be used to constrain source models. With concurrent GPS and field data supporting rapid rise of magma, we infer that tremor resulted from drag of nearly solid magma along rough conduit walls as magma was forced toward the surface.

Washington↗

USGS2014

A continuous record of intereruption velocity change at Mount St. Helens from coda wave interferometry

In September 2004, Mount St. Helens volcano erupted after nearly 18 years of quiescence. However, it is unclear from the limited geophysical observations when or if the magma chamber replenished following the 1980–1986 eruptions in the years before the 2004–2008 extrusive eruption. We use coda wave interferometry with repeating earthquakes to measure small changes in the velocity structure of Mount St. Helens volcano that might indicate magmatic intrusion. By combining observations of relative velocity changes from many closely located earthquake sources, we solve for a continuous function of velocity changes with time. We find that seasonal effects dominate the relative velocity changes. Seismicity rates and repeating earthquake occurrence also vary seasonally; therefore, velocity changes and seismicity are likely modulated by snow loading, fluid saturation, and/or changes in groundwater level. We estimate hydrologic effects impart stress changes on the order of tens of kilopascals within the upper 4 km, resulting in annual velocity variations of 0.5 to 1%. The largest nonseasonal change is a decrease in velocity at the time of the deep M w  = 6.8 Nisqually earthquake. We find no systematic velocity changes during the most likely times of intrusions, consistent with a lack of observable surface deformation. We conclude that if replenishing intrusions occurred, they did not alter seismic velocities where this technique is sensitive due to either their small size or the finite compressibility of the magma chamber. We interpret the observed velocity changes and shallow seasonal seismicity as a response to small stress changes in a shallow, pressurized system.

Journal of Geophysical Research B: Solid Earth↗

USGS2014

The 17 May 2012 M4.8 earthquake near Timpson, East Texas: An event possibly triggered by fluid injection

This study summarizes our investigation of the 17 May 2012 MW-RMT4.8 earthquake near Timpson, Texas, the largest earthquake recorded historically in eastern Texas. To identify preshocks and aftershocks of the 17 May event we examined the arrivals recorded at Nacogdoches (NATX) 30 km from the 17 May epicenter, at nearby USArray Transportable Array stations, and at eight temporary stations deployed between 26 May 2012 and mid-2013. At NATX we identified seven preshocks, the earliest occurring in April 2008. Reliably located aftershocks recorded by the temporary stations lie along a 6 km long NW-SE linear trend corresponding to a previously mapped basement fault that extends across the highest-intensity (MMI VII) region of the 17 May main shock. Earthquakes in this sequence are relatively shallow—with focal depths ranging from 1.6 to 4.6 km. Evidence supporting these depths include: hypocentral locations of exceptionally well-recorded aftershocks, S-P intervals at the nearest stations, and comparisons of synthetics and observed seismograms. Within 3 km of the linear trend of aftershock activity there are two Class II injection disposal wells injecting at 1.9 km depth beginning in August 2006 and February 2007, with injection rates averaging 42,750 m3/mo and 15,600 m3/mo, respectively. Several observations support the hypothesis that fluid injection triggered the Timpson sequence: well-located epicenters are situated near a mapped basement fault and near high-volume injection wells, focal depths are at or below the depths of injection, and the earliest preshock (April 2008) occurred after the onset of injection in 2006.

Texas↗

USGS2014

Triggering and modulation of geyser eruptions in Yellowstone National Park by earthquakes, earth tides, and weather

We analyze intervals between eruptions (IBEs) data acquired between 2001 and 2011 at Daisy and Old Faithful geysers in Yellowstone National Park. We focus our statistical analysis on the response of these geysers to stress perturbations from within the solid earth (earthquakes and earth tides) and from weather (air pressure and temperature, precipitation, and wind). We conclude that (1) the IBEs of these geysers are insensitive to periodic stresses induced by solid earth tides and barometric pressure variations; (2) Daisy (pool geyser) IBEs lengthen by evaporation and heat loss in response to large wind storms and cold air; and (3) Old Faithful (cone geyser) IBEs are not modulated by air temperature and pressure variations, wind, and precipitation, suggesting that the subsurface water column is decoupled from the atmosphere. Dynamic stress changes of 0.1−0.2 MPa resulting from the 2002 M-7.9 Denali, Alaska, earthquake surface waves caused a statistically significant shortening of Daisy geyser's IBEs. Stresses induced by other large global earthquakes during the study period were at least an order of magnitude smaller. In contrast, dynamic stresses of >0.5 MPa from three large regional earthquakes in 1959, 1975, and 1983 caused lengthening of Old Faithful's IBEs. We infer that most subannual geyser IBE variability is dominated by internal processes and interaction with other geysers. The results of this study provide quantitative bounds on the sensitivity of hydrothermal systems to external stress perturbations and have implications for studying the triggering and modulation of volcanic eruptions by external forces.

Wyoming↗

USGS2014

Three-dimensional seismic velocity structure of Mauna Loa and Kilauea volcanoes in Hawaii from local seismic tomography

We present a new three-dimensional seismic velocity model of the crustal and upper mantle structure for Mauna Loa and Kilauea volcanoes in Hawaii. Our model is derived from the first-arrival times of the compressional and shear waves from about 53,000 events on and near the Island of Hawaii between 1992 and 2009 recorded by the Hawaiian Volcano Observatory stations. The V p model generally agrees with previous studies, showing high-velocity anomalies near the calderas and rift zones and low-velocity anomalies in the fault systems. The most significant difference from previous models is in V p / V s structure. The high- V p and high- V p / V s anomalies below Mauna Loa caldera are interpreted as mafic magmatic cumulates. The observed low- V p and high- V p / V s bodies in the Kaoiki seismic zone between 5 and 15 km depth are attributed to the underlying volcaniclastic sediments. The high- V p and moderate- to low- V p / V s anomalies beneath Kilauea caldera can be explained by a combination of different mafic compositions, likely to be olivine-rich gabbro and dunite. The systematically low- V p and low- V p / V s bodies in the southeast flank of Kilauea may be caused by the presence of volatiles. Another difference between this study and previous ones is the improved V p model resolution in deeper layers, owing to the inclusion of events with large epicentral distances. The new velocity model is used to relocate the seismicity of Mauna Loa and Kilauea for improved absolute locations and ultimately to develop a high-precision earthquake catalog using waveform cross-correlation data.

Hawaii↗