Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Applied Geophysics”

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 253 records · Page 14Linked to original sources
USGS2022

Synthetic evaluation of infrasonic multipole waveform inversion

Acoustic source inversions estimate the mass flow rate of volcanic explosions or yield of chemical explosions and provide insight into potential source directionality. However, the limitations of applying these methods to complex sources and their ability to resolve a stable solution have not been investigated in detail. We perform synthetic infrasound waveform inversions that use 3-D Green’s functions for a variety of idealized and realistic deployment scenarios using both a flat plane and Yasur volcano, Vanuatu as examples. We investigate the ability of various scenarios to retrieve the input source functions and relative amplitudes for monopole and multipole (monopole and dipole) inversions. Infrasound waveform inversions appear to be a robust method to quantify mass flow rates from simple sources (monopole) using deployments of infrasound sensors placed around a source, but care should be taken when analyzing and interpreting results from more complex acoustic sources (multipole) that have significant directional components. In the examples we consider the solution is stable for monopole inversions with a signal-to-noise ratio greater than five and the dipole component is small. For most scenarios investigated, the vertical dipole component of the multipole explosion source is poorly constrained and can impact the ability to recover the other source term components. Because multipole inversions are ill-posed for many deployments, a low residual does not necessarily mean the proper source vector has been recovered. Synthetic studies can help investigate the limitations and place bounds on information that may be missing using monopole and multipole inversions for potentially directional sources.

Journal of Geophysical Research-Solid Earth↗
USGS2022

A model of the spatiotemporal dynamics of soil carbon following coastal wetland loss applied to a Louisiana salt marsh in the Mississippi River Deltaic Plain

The potential for carbon sequestration in coastal wetlands is high due to protection of carbon (C) in flooded soils. However, excessive flooding can result in the conversion of the vegetated wetland to open water. This transition results in the loss of wetland habitat in addition to the potential loss of soil carbon. Thus, in areas experiencing rapid wetland submergence, such as the Mississippi River Delta, coastal wetlands could become a significant source of carbon emissions if land loss is not mitigated. To accurately assess the capacity of wetlands to store (or emit) carbon in dynamic environments, it is critical to understand the fate of soil carbon following the transition from vegetated wetland to open water. We developed a simple soil carbon model representing soil depths to 1 m using the data collected from a Louisiana coastal salt marsh in the Mississippi River Deltaic Plain to predict soil carbon density and stock following the transition from a vegetated salt marsh to an open water pond. While immediate effects of ponding on the distribution of carbon within the 1-m soil profile were apparent, there were no effects of ponding on the overall, integrated, carbon stocks 14 years, following wetland submergence. Rather, the model predicts that soil carbon losses in the first meter will be realized over long periods of time (∼200 years) due to changes in the source of carbon (biomass vs. mineral sediment) with minimal losses through mineralization.

Louisiana↗
USGS2012

Volcanism, isostatic residual gravity and regional tectonic setting of the Cascade volcanic province

A technique to locate automatically boundaries between crustal blocks of disparate densities was applied to upward continued isostatic residual gravity data. The boundary analysis delineates a narrow gravitational trough that extends the length of the Pliocene and Quaternary volcanic arc from Mount Baker in northern Washington to Lassen Peak in California. Gravitational highs interrupt the trough at two localities: a northwest trending high in southern Washington and a northeast trending high between Mount Shasta and Lassen Peak. The latter anomaly is one of a set of northeast trending anomalies that, within the Quaternary arc, appear related to volcanic segmentation proposed previously on the basis of spatial and compositional distributions of volcanoes. These northeast trending anomalies extend hundreds of kilometers northeast of the arc, are caused by sources in the upper crust, and in some cases are related to exposed pre-Tertiary rocks. Segmentation models invoke geometric characteristics of the subducting plate as the primary factor controlling location and chemistry of volcanism, and these northeast trending gravity sources also may be a product of disturbance of the upper crust by the subduction process. More likely, the gravity sources may reflect upper crustal structures older than the High Cascades, possibly relicts from earlier accretionary events or more recent crustal deformation, that have actively influenced the spatial location of more recent volcanism. Much of the Pliocene and Quaternary volcanism of the Cascade arc has concentrated on or near contacts between crustal blocks of disparate density. These contacts may promote the ascension of magma to the Earth's surface.

Journal of Geophysical Research Solid Earth↗
USGS1997

Empirical Green's function analysis: Taking the next step

An extension of the empirical Green's function (EGF) method is presented that involves determination of source parameters using standard EGF deconvolution, followed by inversion for a common attenuation parameter for a set of colocated events. Recordings of three or more colocated events can thus be used to constrain a single path attenuation estimate. I apply this method to recordings from the 1995–1996 Ridgecrest, California, earthquake sequence; I analyze four clusters consisting of 13 total events with magnitudes between 2.6 and 4.9. I first obtain corner frequencies, which are used to infer Brune stress drop estimates. I obtain stress drop values of 0.3–53 MPa (with all but one between 0.3 and 11 MPa), with no resolved increase of stress drop with moment. With the corner frequencies constrained, the inferred attenuation parameters are very consistent; they imply an average shear wave quality factor of approximately 20–25 for alluvial sediments within the Indian Wells Valley. Although the resultant spectral fitting (using corner frequency and κ) is good, the residuals are consistent among the clusters analyzed. Their spectral shape is similar to the the theoretical one-dimensional response of a layered low-velocity structure in the valley (an absolute site response cannot be determined by this method, because of an ambiguity between absolute response and source spectral amplitudes). I show that even this subtle site response can significantly bias estimates of corner frequency and κ, if it is ignored in an inversion for only source and path effects. The multiple-EGF method presented in this paper is analogous to a joint inversion for source, path, and site effects; the use of colocated sets of earthquakes appears to offer significant advantages in improving resolution of all three estimates, especially if data are from a single site or sites with similar site response.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Imaging spectroscopy: Earth and planetary remote sensing with the USGS Tetracorder and expert systems

Imaging spectroscopy is a tool that can be used to spectrally identify and spatially map materials based on their specific chemical bonds. Spectroscopic analysis requires significantly more sophistication than has been employed in conventional broadband remote sensing analysis. We describe a new system that is effective at material identification and mapping: a set of algorithms within an expert system decision‐making framework that we call Tetracorder. The expertise in the system has been derived from scientific knowledge of spectral identification. The expert system rules are implemented in a decision tree where multiple algorithms are applied to spectral analysis, additional expert rules and algorithms can be applied based on initial results, and more decisions are made until spectral analysis is complete. Because certain spectral features are indicative of specific chemical bonds in materials, the system can accurately identify and map those materials. In this paper we describe the framework of the decision making process used for spectral identification, describe specific spectral feature analysis algorithms, and give examples of what analyses and types of maps are possible with imaging spectroscopy data. We also present the expert system rules that describe which diagnostic spectral features are used in the decision making process for a set of spectra of minerals and other common materials. We demonstrate the applications of Tetracorder to identify and map surface minerals, to detect sources of acid rock drainage, and to map vegetation species, ice, melting snow, water, and water pollution, all with one set of expert system rules. Mineral mapping can aid in geologic mapping and fault detection and can provide a better understanding of weathering, mineralization, hydrothermal alteration, and other geologic processes. Environmental site assessment, such as mapping source areas of acid mine drainage, has resulted in the acceleration of site cleanup, saving millions of dollars and years in cleanup time. Imaging spectroscopy data and Tetracorder analysis can be used to study both terrestrial and planetary science problems. Imaging spectroscopy can be used to probe planetary systems, including their atmospheres, oceans, and land surfaces.

Journal of Geophysical Research E: Planets↗
USGS2020

Causes of variability in suspended‐sand concentration evaluated using measurements in the Colorado River in Grand Canyon

Rivers commonly exhibit substantial variability in suspended‐sand concentration, even at constant water discharge. Here we derive an approach for evaluating how much of this variability arises from mean bed‐sand grain size. We apply this approach to the Colorado River in Grand Canyon, where discharge‐independent concentration of suspended sand varies by more than a factor of 23 ( N = 1.4 × 10 6 ). Theory predicts that where concentration is controlled by bed‐sand grain size, concentration and grain size in suspension will be inversely correlated (i.e., coarsening of the bed causes suspended sand to become coarser in grain size and lower in concentration). Although the observed correlation is negative, riverbed grain size accounts for only 40% of the variability in concentration. The residuals vary by an order of magnitude; they arise from other processes, such as changes in topography or distribution of sand that cause shear stress to change at constant discharge, changes in the fine tail of bed‐sand grain sizes or changing bedforms. Both bed sand and the other factors influence concentration for durations from less than 1 day to several years. Predictions of concentration based on bed‐sand grain size ( N = 4 × 10 4 ) are less accurate than predictions based on suspended‐sand grain size, probably because suspended sand is a natural integrator of sand‐transporting processes, giving more weight to those areas of the bed that exchange more sand with the flow. Although the causes of variability vary from one river to another, the approach illustrated here is applicable to any river in which concentration varies at constant water discharge.

Arizona↗
USGS2025

A joint Gaussian process model of geochemistry, geophysics, and temperature for groundwater TDS in the San Ardo Oil Field, California, USA

Decline in availability of fresh groundwater has expanded interest in brackish groundwater resources; however, the distribution of brackish groundwater is poorly understood. Water resources in sedimentary basins across the United States often overlie oil and gas development. Mapping of groundwater total dissolved solids (TDS) using data from oil well geophysical logs has become an important technique for identifying fresh and brackish groundwater. Existing geophysical log analysis methods use porosity and temperature to relate formation resistivity to TDS. Typically, natural geothermal gradients are used to estimate temperature at the location of collected resistivity. However, in thermally enhanced oil fields, steam is injected into the subsurface to mobilize high viscosity oil, creating variable temperature distributions. Furthermore, TDS derived from resistivity also depends on the fractions of dominant ions. Typically, chloride and bicarbonate fractions must be determined. It is also necessary to model TDS across many geologic units with heterogenous porosity distributions. Collectively, each quantity used to estimate TDS (resistivity, porosity, temperature, bicarbonate fraction) varies in space and time, and available data points are rarely collocated. Here, we present a new method of mapping groundwater TDS that continuously models each quantity together with a joint Gaussian process. This method enables mapping fresh and brackish water with practically available data. We apply this method to the San Ardo Oil Field in Monterey County, California, where steam injection occurs. In some areas of the aquifer system overlying the oil zone, the temperature is ∼75 °C, roughly twice the natural background value. Groundwater TDS is typically <1,500 mg/L in the aquifer and increases with depth to ∼9,000 mg/L in the oil-producing zone. A low-permeability clay layer delineates the fresh and brackish water, likely by inhibiting surface recharge from penetrating the deeper zones, allowing higher-TDS connate water to remain in place. Weaker lateral TDS trends may be controlled by recharge patterns associated with the Salinas River. Our model reveals with high certainty that groundwater has freshened in one localized part of the oil-producing zone and suggests with less certainty that more widespread freshening has also occurred. The lowering of TDS was possibly from decades of low-TDS steam injection and the associated fluid production and disposal operations.

California↗
USGS2012

Excitation of a buried magmatic pipe: A seismic source model for volcanic tremor

Recent observations of seismic events at various volcanoes suggest that harmonic tremor results from the sustained occurrence of so-called long-period or low-frequency events. Accordingly, we can view the long-period volcanic event as the elementary process of tremor and interpret it as the impulse response of the tremor-generating system. We present a seismic model in which the source of tremor is the acoustic resonance of a fluid-filled volcanic pipe triggered by excess gas pressure. The model consists of three elements, namely, a triggering mechanism, a resonator, and a radiator. For simplicity, we assume a hemispherical trigger, cylindrial resonator, and circular radiator set in a vertical configuration with the trigger capping the top of the pipe and the disk-shaped radiator shutting off its bottom. Considering the simple case of a source buried in a homogeneous half space, we then apply the discrete wave number method to obtain a complete representation of the ground motion response at near and intermediate distances. The results demonstrate that the displacement attributed to the pipe dominates the near-field motion, while that due to the disk is representative of the intermediate and far fields. The trigger itself has a smaller contribution, mainly limited to the field in the proximity of the source. The characteristics displayed by the free surface response evolve from a strong impulsive signature in the immediate vicinity of the epicenter to a well-developed harmonic wave train dominated by Rayleigh waves at larger distances. No clear shear arrival can be detected in the synthetic seismograms. The displacement spectrum reflects the organ-pipe modes of the conduit, and the bandwidth associated with the dominant spectral peak of motion is controlled by the combined losses due to viscous attenuation in the fluid and elastic radiation into the solid. In the case of the cylindrical magma column considered, the radiation loss is proportional to the square of the pipe radius, while the loss related to viscous damping is inversely proportional to the same factor, indicating that the relative importance of the two loss mechanisms is critically dependent on the geometry of the magma reservoir. The relative importance of the pipe and disk elements, likewise, is a function of the conduit cross section. This suggests the possibility of determining the geometry of the source as well as the radiation loss and in situ magma viscosity from a comparison of near- and far-field observations.

Journal of Geophysical Research Solid Earth↗
USGS2018

Three-dimensional modeling of fine sediment transport by waves and currents in a shallow estuary

A suspended sediment transport model is implemented in the unstructured‐grid SUNTANS model and applied to study fine‐grained sediment transport in South San Francisco Bay. The model enables calculation of suspension of bottom sediment based on combined forcing of tidal currents and wind waves. We show that accurate results can be obtained by employing two‐size classes which are representative of microflocs and macroflocs in the Bay. A key finding of the paper is that the critical calibration parameter is the ratio of the erosion of the microflocs to macroflocs from the bed. Different values of this erosion ratio are needed on the shallow shoals and deeper channels because of the different nature of the sediment dynamics in these regions. Application of a spatially variable erosion ratio and critical shear stress for erosion is shown to accurately reproduce observed suspended sediment concentration at four‐field sites located along a cross‐channel transect. The results reveal a stark contrast between the behavior of the suspended sediment concentration on the shoals and in the deep channel. Waves are shown to resuspend sediments on the shoals, although tidal and wind‐generated currents are needed to mix the thin wave‐driven suspensions into the water column. The contribution to the suspended sediment concentration in the channel by transport from the shoals is similar in magnitude to that due to local resuspension. However, the local contribution is in phase with strong bottom currents which resuspend the sediments, while the contribution from the shoals peaks during low‐water slack tide.

San Francisco Bay↗
USGS1998

Earthquake triggering by transient and static deformations

Observational evidence for both static and transient near-field and far-field triggered seismicity are explained in terms of a frictional instability model, based on a single degree of freedom spring-slider system and rate- and state-dependent frictional constitutive equations. In this study a triggered earthquake is one whose failure time has been advanced by Δ t (clock advance) due to a stress perturbation. Triggering stress perturbations considered include square-wave transients and step functions, analogous to seismic waves and coseismic static stress changes, respectively. Perturbations are superimposed on a constant background stressing rate which represents the tectonic stressing rate. The normal stress is assumed to be constant. Approximate, closed-form solutions of the rate-and-state equations are derived for these triggering and background loads, building on the work of Dieterich [1992, 1994]. These solutions can be used to simulate the effects of static and transient stresses as a function of amplitude, onset time t 0 , and in the case of square waves, duration. The accuracies of the approximate closed-form solutions are also evaluated with respect to the full numerical solution and t 0 . The approximate solutions underpredict the full solutions, although the difference decreases as t 0 approaches the end of the earthquake cycle. The relationship between Δ t and t 0 differs for transient and static loads: a static stress step imposed late in the cycle causes less clock advance than an equal step imposed earlier, whereas a later applied transient causes greater clock advance than an equal one imposed earlier. For equal Δ t , transient amplitudes must be greater than static loads by factors of several tens to hundreds depending on t 0 . We show that the rate-and-state model requires that the total slip at failure is a constant, regardless of the loading history. Thus a static load applied early in the cycle, or a transient applied at any time, reduces the stress at the initiation of failure, whereas static loads that are applied sufficiently late raise it. Rate-and-state friction predictions differ markedly from those based on Coulomb failure stress changes (ΔCFS) in which Δ t equals the amplitude of the static stress change divided by the background stressing rate. The ΔCFS model assumes a stress failure threshold, while the rate-and-state equations require a slip failure threshold. The complete rate-and-state equations predict larger Δ t than the ΔCFS model does for static stress steps at small t 0 , and smaller Δ t than the ΔCFS model for stress steps at large t 0 . The ΔCFS model predicts nonzero Δ t only for transient loads that raise the stress to failure stress levels during the transient. In contrast, the rate-and-state model predicts nonzero Δ t for smaller loads, and triggered failure may occur well after the transient is finished. We consider heuristically the effects of triggering on a population of faults, as these effects might be evident in seismicity data. Triggering is manifest as an initial increase in seismicity rate that may be followed by a quiescence or by a return to the background rate. Available seismicity data are insufficient to discriminate whether triggered earthquakes are “new” or clock advanced. However, if triggering indeed results from advancing the failure time of inevitable earthquakes, then our modeling suggests that a quiescence always follows transient triggering and that the duration of increased seismicity also cannot exceed the duration of a triggering transient load. Quiescence follows static triggering only if the population of available faults is finite.

Journal of Geophysical Research B: Solid Earth↗
USGS2025

Identifying precursors and tracking pulses of magma ascent in multidisciplinary data during the 2018–2023 phreatomagmatic eruption at Semisopochnoi Island, Alaska

The 2018–2023 phreatomagmatic eruptions at Semisopochnoi Island, Alaska produced abundant long-period (LP) seismicity, harmonic and broadband tremor, and explosion signals over several well-monitored periods of eruption and quiescence. The corresponding dataset provides an excellent opportunity to investigate precursory and syn -eruptive geophysical signals of long-lived phreatomagmatic eruptions using multiparameter observations. We generated explosion and LP event catalogs through novel implementations of the REDPy ( Hotovec-Ellis, 2024 ) repeating event detector in mid-2021 following a network upgrade and the onset of a new phase of the eruption. The hundreds of detected explosions show a high degree of infrasound waveform similarity over more than a year, indicating a repeating source mechanism likely associated with explosive magma-water interaction. The seismic LP catalog shows that events began over a month prior to renewed explosive activity at the beginning of August 2021, and that lower frequency index (FI) LPs were generated in the week prior to the onset of explosions. We applied a recently developed machine learning tool (VOISS-Net, Tan et al., 2024 ) to catalog abundant broadband and harmonic seismic tremor recorded before and during the renewed explosive activity, along with LPs and explosions. The tremor catalogs complement the LP and explosion catalogs by filling out the seismic sequence with the dominant signal types. Together, these catalogs reveal a seismic sequence of renewed unrest that started with several weeks of LP events, followed by LPs with lower FI values and harmonic tremor in the days prior to explosive activity, and finally the onset of discrete explosions and broadband eruption tremor. We interpret this sequence as the ascent of a new pulse of magma that first interacted with the hydrothermal/groundwater system to produce LPs, followed by harmonic tremor, and that ultimately drove explosive magma-water interactions and periods of continuous ash emissions. The 2021 seismic sequence, in combination with long-term records of satellite SO 2 emissions, deformation from interferometric synthetic aperture radar (InSAR) analysis, ash sample analysis, infrasound, and volcano tectonic seismicity, allows us to interpret the entire 9-year period of unrest and eruption that began with an intrusion and earthquake swarm in 2014.

Alaska↗
USGS2012

Interpretation of gravity data in a complex volcano-tectonic setting, southwestern Nevada

This regional gravity study, based on an irregular 2-km data grid, was conducted during the past few years at Yucca Mountain, southern Nye County, Nevada, as part of a program to locate a suitable repository for high-level nuclear waste. About 100 surface rock samples, three borehole gamma-gamma logs, and one borehole gravity study provide excellent density control. A nearly linear increase in density of 0.26 g/cm 3 per kilometer of depth is indicated in the thick tuff sequences that underlie the mountain. Isostatic and 2.0-g/cm 3 Bouguer corrections were applied to the observed gravity values to remove regional gradients and topographic effects, respectively. The Bare Mountain gravity high, with an isostatic anomaly maximum of 48 mGal, is connected with a greater gravity high over the Funeral Mountains, to the southwest; together, these highs result from a continuous block of dense, metamorphosed Precambrian and Paleozoic rocks that stretches across much of the Walker Lane from the east edge of Death Valley to Bare Mountain. The Calico Hills gravity high appears more likely to originate from a northeast trending buried ridge of Paleozoic rocks that extends southwestward beneath Busted Butte, 5 km southeast of the proposed repository, where two- and three-dimensional modeling indicates that the pre-Cenozoic rocks lie less than 1000 m beneath the surface. Tuff, at least 4000 m thick, fills a large steep-sided depression in the pretuff rocks beneath Yucca Mountain and Crater Flat. The gravity low and the thick tuff section lie within a large collapse area that includes the Crater Flat-Timber Mountain-Silent Canyon caldera complexes. Gravity lows in Crater Flat itself are interpreted to coincide with the source areas of the Prow Pass Member, the Bullfrog Member, and the Tram Member of the Crater Flat Tuff; these source areas add nearly 350 km 2 to the previously recognized extent of the local caldera complexes. Southward extension of the broad gravity low associated with Crater Flat into the Amargosa Desert is evidence for sector graben-type collapse segments related to the formation of the Timber Mountain caldera and superimposed on the other volcanic and extensional structures within Crater Flat.

Journal of Geophysical Research Solid Earth↗
USGS2010

Viscoelastic-cycle model of interseismic deformation in the northwestern United States

We apply a viscoelastic cycle model to a compilation of GPS velocity fields in order to address the kinematics of deformation in the northwestern United States. A viscoelastic cycle model accounts for time-dependent deformation following large crustal earthquakes and is an alternative to block models for explaining the interseismic crustal velocity field. Building on the approach taken in Pollitz et al., we construct a deformation model for the entire western United States-based on combined fault slip and distributed deformation-and focus on the implications for the Mendocino triple junction (MTJ), Cascadia megathrust, and western Washington. We find significant partitioning between strike-slip and dip-slip motion near the MTJ as the tectonic environment shifts from northwest-directed shear along the San Andreas fault system to east-west convergence along the Juan de Fuca Plate. By better accounting for the budget of aseismic and seismic slip along the Cascadia subduction interface in conjunction with an assumed rheology, we revise a previous model of slip for the M~ 9 1700 Cascadia earthquake. In western Washington, we infer slip rates on a number of strike-slip and dip-slip faults that accommodate northward convergence of the Oregon Coast block and northwestward convergence of the Juan de Fuca Plate. Lateral variations in first order mechanical properties (e.g. mantle viscosity, vertically averaged rigidity) explain, to a large extent, crustal strain that cannot be rationalized with cyclic deformation on a laterally homogeneous viscoelastic structure. Our analysis also shows that present crustal deformation measurements, particularly with the addition of the Plate Boundary Observatory, can constrain such lateral variations.

Geophysical Journal International↗
USGS2013

Very long period conduit oscillations induced by rockfalls at Kilauea Volcano, Hawaii

Eruptive activity at the summit of Kilauea Volcano, Hawaii, beginning in 2010 and continuing to the present time is characterized by transient outgassing bursts accompanied by very long period (VLP) seismic signals triggered by rockfalls from the vent walls impacting a lava lake in a pit within the Halemaumau pit crater. We use raw data recorded with an 11-station broadband network to model the source mechanism of signals accompanying two large rockfalls on 29 August 2012 and two smaller average rockfalls obtained by stacking over all events with similar waveforms to improve the signal-to-noise ratio. To determine the source centroid location and source mechanism, we minimize the residual error between data and synthetics calculated by the finite difference method for a point source embedded in a homogeneous medium that takes topography into account. We apply a new waveform inversion method that accounts for the contributions from both translation and tilt in horizontal seismograms through the use of Green's functions representing the seismometer response to translation and tilt ground motions. This method enables a robust description of the source mechanism over the period range 1–1000 s. The VLP signals associated with the rockfalls originate in a source region ∼1 km below the eastern perimeter of the Halemaumau pit crater. The observed waveforms are well explained by a simple volumetric source with geometry composed of two intersecting cracks including an east striking crack (dike) dipping 80° to the north, intersecting a north striking crack (another dike) dipping 65° to the east. Each rockfall is marked by a similar step-like inflation trailed by decaying oscillations of the volumetric source, attributed to the efficient coupling at the source centroid location of the pressure and momentum changes induced by the rock mass impacting the top of the lava column. Assuming a simple lumped parameter representation of the shallow magmatic system, the observed pressure and volume variations can be modeled with the following attributes: rockfall volume (200–4500 m 3 ), length of magma column (120–210 m), diameter of pipe connecting the Halemaumau pit crater to the subjacent dike system (6 m), average thickness of the two underlying dikes (3–6 m), and effective magma viscosity (30–210 Pa s). Most rockfalls occur during episodes of sustained deflation of the Kilauea summit. The mass loss rate in the shallow magmatic system is estimated to be 1400–15,000 kg s −1 based on measurements of the temporal variation of VLP period in the two large rockfalls that occurred on 29 August 2012.

Hawaii↗
USGS2018

Findings and lessons learned from the assessment of the Mexico-United States transboundary San Pedro and Santa Cruz aquifers: The utility of social science in applied hydrologic research

Study Region This study region encompasses the Transboundary San Pedro and Santa Cruz aquifers which are shared between the states of Sonora (Mexico) and Arizona (US). Special regional considerations include a semi-arid climate, basin-fill aquifers with predominantly montane recharge areas, economic drivers in the mining, trade, and military sectors, groundwater-dependent cities with expanding cones of depression, interbasin groundwater transfers, ground- and surface-water contamination, and protected aquatic and riparian habitats that act as significant migration corridors for hundreds of species, including some that are threatened and endangered. Study Focus We focus on lessons learned from the hydrologic assessment of the Transboundary San Pedro and Santa Cruz aquifers. We conducted the work, in two phases: (1) laying the groundwork and (2) implementation. The “laying the groundwork” phase consisted of binational meetings with stakeholders and key actors (agencies and individuals), and the development of an understanding of the physical, institutional, historical, and socio-political context. This led to signing of the binational Transboundary Aquifer Assessment Program (TAAP) agreement in 2009 and detailed the process for cooperation and coordination in the assessment of shared aquifers. The implementation phase began with an agreement to proceed with the study of four “focus” aquifers (Santa Cruz, San Pedro, Mesilla (Conejos-Médanos in Mexico), and Hueco Bolson (Bolsón del Hueco in Mexico)) and development of associated technical teams. Though we do include a brief discussion of the lessons learned from the physical science portion of the study, the results have been described and published elsewhere. The bulk of the paper instead focuses on the findings and lessons learned from the integration of social-science perspectives into a largely physical-science based program, since there is a growing recognition of the need for this type of approach especially in the management and assessment of transboundary aquifers. New Hydrological Insights for the Region The Sonora-Arizona effort succeeded because both countries were adequately represented, and because of flexibility of skills and ability of teams comprising both university and government scientists. Teams included social and earth scientists. Including the social sciences was critical to research design and implementation, and to addressing the cultural, institutional, and socio-political contexts of transboundary aquifer assessment. Significant components of the continuing implementation phase include strategic planning, data compilation and analysis, cross-border integration of datasets, geophysical and geochemical surveys, and internal, peer, and stakeholder engagement.

Arizona, Sonora↗
USGS2015

Pore-pressure sensitivities to dynamic strains: observations in active tectonic regions

Triggered seismicity arising from dynamic stresses is often explained by the Mohr-Coulomb failure criterion, where elevated pore pressures reduce the effective strength of faults in fluid-saturated rock. The seismic response of a fluid-rock system naturally depends on its hydro-mechanical properties, but accurately assessing how pore-fluid pressure responds to applied stress over large scales in situ remains a challenging task; hence, spatial variations in response are not well understood, especially around active faults. Here I analyze previously unutilized records of dynamic strain and pore-pressure from regional and teleseismic earthquakes at Plate Boundary Observatory (PBO) stations from 2006 through 2012 to investigate variations in response along the Pacific/North American tectonic plate boundary. I find robust scaling-response coefficients between excess pore pressure and dynamic strain at each station that are spatially correlated: around the San Andreas and San Jacinto fault systems, the response is lowest in regions of the crust undergoing the highest rates of secular shear strain. PBO stations in the Parkfield instrument cluster are at comparable distances to the San Andreas fault (SAF), and spatial variations there follow patterns in dextral creep rates along the fault, with the highest response in the actively creeping section, which is consistent with a narrowing zone of strain accumulation seen in geodetic velocity profiles. At stations in the San Juan Bautista (SJB) and Anza instrument clusters, the response depends non-linearly on the inverse fault-perpendicular distance, with the response decreasing towards the fault; the SJB cluster is at the northern transition from creeping-to-locked behavior along the SAF, where creep rates are at moderate to low levels, and the Anza cluster is around the San Jacinto fault, where to date there have been no statistically significant creep rates observed at the surface. These results suggest that the strength of the pore pressure response in fluid-saturated rock near active faults is controlled by shear strain accumulation associated with tectonic loading, which implies a strong feedback between fault strength and permeability: dynamic triggering susceptibilities may vary in space and also in time.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Waveform inversion of very long period impulsive signals associated with magmatic injection beneath Kilauea Volcano, Hawaii

We use data from broadband seismometers deployed around the summit of Kilauea Volcano to quantify the mechanism associated with a transient in the flow of magma feeding the east rift eruption of the volcano. The transient is marked by rapid inflation of the Kilauea summit peaking at 22 μrad 4.5 hours after the event onset, followed by slow deflation over a period of 3 days. Superimposed on the summit inflation is a series of sawtooth displacement pulses, each characterized by a sudden drop in amplitude lasting 5–10 s followed by an exponential recovery lasting 1–3 min. The sawtooth waveforms display almost identical shapes, suggesting a process involving the repeated activation of a fixed source. The particle motion associated with each sawtooth is almost linear, and its major swing shows compressional motion at all stations. Analyses of semblance and particle motion are consistent with a point source located 1 km beneath the northeast edge of the Halemaumau pit crater. To estimate the source mechanism, we apply a moment tensor inversion to the waveform data, assuming a point source embedded in a homogeneous half-space with compressional and shear wave velocities representative of the average medium properties at shallow depth under Kilauea. Synthetic waveforms are constructed by a superposition of impulse responses for six moment tensor components and three single force components. The origin times of individual impulses are distributed along the time axis at appropriately small, equal intervals, and their amplitudes are determined by least squares. In this inversion, the source time functions of the six tensor and three force components are determined simultaneously. We confirm the accuracy of the inversion method through a series of numerical tests. The results from the inversion show that the waveform data are well explained by a pulsating transport mechanism operating on a subhorizontal crack linking the summit reservoir to the east rift of Kilauea. The crack acts like a buffer in which a batch of fluid (magma and/or gas) accumulates over a period of 1–3 min before being rapidly injected into a larger reservoir (possibly the east rift) over a timescale of 5–10 s. The seismic moment and volume change associated with a typical batch of fluid are approximately 10 14 N m and 3000 m 3 , respectively. Our results also point to the existence of a single force component with amplitude of 10 9 N, which may be explained as the drag force generated by the flow of viscous magma through a narrow constriction in the flow path. The total volume of magma associated with the 4.5-hour-long activation of the pulsating source is roughly 500,000 m 3 in good agreement with the integrated volume flow rate of magma estimated near the eruptive site.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

The global short-period wavefield modelled with a Monte Carlo seismic phonon method

At high frequencies (∼1 Hz), much of the seismic energy arriving at teleseismic distances is not found in the main phases (e.g. P , PP , S , etc.) but is contained in the extended coda that follows these arrivals. This coda results from scattering off small-scale velocity and density perturbations within the crust and mantle and contains valuable information regarding the depth dependence and strength of this heterogeneity as well as the relative importance of intrinsic versus scattering attenuation. Most analyses of seismic coda to date have concentrated on S -wave coda generated from lithospheric scattering for events recorded at local and regional distances. Here, we examine the globally averaged vertical-component, 1-Hz wavefield (>10° range) for earthquakes recorded in the IRIS FARM archive from 1990 to 1999. We apply an envelope-function stacking technique to image the average time–distance behavior of the wavefield for both shallow (≤50 km) and deep (≥500 km) earthquakes. Unlike regional records, our images are dominated by P and P coda owing to the large effect of attenuation on PP and S at high frequencies. Modelling our results is complicated by the need to include a variety of ray paths, the likely contributions of multiple scattering and the possible importance of P -to- S and S -to- P scattering. We adopt a stochastic, particle-based approach in which millions of seismic phonons are randomly sprayed from the source and tracked through the Earth. Each phonon represents an energy packet that travels along the appropriate ray path until it is affected by a discontinuity or a scatterer. Discontinuities are modelled by treating the energy normalized reflection and transmission coefficients as probabilities. Scattering probabilities and scattering angles are computed in a similar fashion, assuming random velocity and density perturbations characterized by an exponential autocorrelation function. Intrinsic attenuation is included by reducing the energy contained in each particle as an appropriate function of traveltime. We find that most scattering occurs in the lithosphere and upper mantle, as previous results have indicated, but that some lower-mantle scattering is likely also required. A model with 3 to 4 per cent rms velocity heterogeneity at 4-km scale length in the upper mantle and 0.5 per cent rms velocity heterogeneity at 8-km scale length in the lower mantle (with intrinsic attenuation of Q α = 450 above 200 km depth and Q α = 2500 below 200 km) provides a reasonable fit to both the shallow- and deep-earthquake observations, although many trade-offs exist between the scale length, depth extent and strength of the heterogeneity.

Geophysical Journal International↗