SEARCH · Search USGS

Results for “Journal of Applied Mechanics”

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.

124 records · Page 7Linked to original sources

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↗

USGS1993

The role of acoustic emission in the study of rock fracture

The development of faults and shear fracture systems over a broad range of temperature and pressure and for a variety of rock types involves the growth and interaction of microcracks. Acoustic emission (AE), which is produced by rapid microcrack growth, is a ubiquitous phenomenon associated with brittle fracture and has provided a wealth of information regarding the failure process in rock. This paper reviews the successes and limitations of AE studies as applied to the fracture process in rock with emphasis on our ability to predict rock failure. Application of laboratory AE studies to larger scale problems related to the understanding of earthquake processes is also discussed. In this context, laboratory studies can be divided into the following categories. 1) Simple counting of the number of AE events prior to sample failure shows a correlation between AE rate and inelastic strain rate. Additional sorting of events by amplitude has shown that AE events obey the power law frequency-magnitude relation observed for earthquakes. These cumulative event count techniques are being used in conjunction with damage mechanics models to determine how damage accumulates during loading and to predict failure. 2) A second area of research involves the location of hypocenters of AE source events. This technique requires precise arrival time data of AE signals recorded over an array of sensors that are essentially a miniature seismic net. Analysis of the spatial and temporal variation of event hypocenters has improved our understanding of the progression of microcrack growth and clustering leading to rock failure. Recently, fracture nucleation and growth have been studied under conditions of quasi-static fault propagation by controlling stress to maintain constant AE rate. 3) A third area of study involves the analysis of full waveform data as recorded at receiver sites. One aspect of this research has been to determine fault plane solutions of AE source events from first motion data. These studies show that in addition to pure tensile and double couple events, a significant number of more complex event types occur in the period leading to fault nucleation. 4) P and S wave velocities (including spatial variations) and attenuation have been obtained by artificially generating acoustic pulses which are modified during passage through the sample. ?? 1993.

International Journal of Rock Mechanics and Mining↗

USGS2016

Cross-seasonal effects on waterfowl productivity: Implications under climate change

Previous efforts to relate winter-ground precipitation to subsequent reproductive success as measured by the ratio of juveniles to adults in the autumn failed to account for increased vulnerability of juvenile ducks to hunting and uncertainty in the estimated age ratio. Neglecting increased juvenile vulnerability will positively bias the mean productivity estimate, and neglecting increased vulnerability and estimation uncertainty will positively bias the year-to-year variance in productivity because raw age ratios are the product of sampling variation, the year-specific vulnerability, and year-specific reproductive success. Therefore, we estimated the effects of cumulative winter precipitation in the California Central Valley and the Mississippi Alluvial Valley on pintail ( Anas acuta ) and mallard ( Anas platyrhnchos ) reproduction, respectively, using hierarchical Bayesian methods to correct for sampling bias in productivity estimates and observation error in covariates. We applied the model to a hunter-collected parts survey implemented by the United States Fish and Wildlife Service and band recoveries reported to the United States Geological Survey Bird Banding Laboratory using data from 1961 to 2013. We compared our results to previous estimates that used simple linear regression on uncorrected age ratios from a smaller subset of years in pintail (1961–1985). Like previous analyses, we found large and consistent effects of population size and wetland conditions in prairie Canada on mallard productivity, and large effects of population size and mean latitude of the observed breeding population on pintail productivity. Unlike previous analyses, we report a large amount of uncertainty in the estimated effects of wintering-ground precipitation on pintail and mallard productivity, with considerable uncertainty in the sign of the estimated main effect, although the posterior medians of precipitation effects were consistent with past studies. We found more consistent estimates in the sign of an interaction effect between population size and precipitation, suggesting that wintering-ground precipitation has a larger effect in years of high population size, especially for pintail. When we used the estimated effects in a population model to derive a sustainable harvest and population size projection (i.e., a yield curve), there was considerable uncertainty in the effect of increased or decreased wintering-ground precipitation on sustainable harvest potential and population size. These results suggest that the mechanism of cross-seasonal effects between winter habitat and reproduction in ducks occurs through a reduction in the strength of density dependence in years of above-average wintering-ground precipitation. We suggest additional investigation of the underlying mechanisms and that habitat managers and decision-makers consider the level of uncertainty in these estimates when attempting to integrate habitat management and harvest management decisions. Collection of annual data on the status of wintering-ground habitat in a rigorous sampling framework would likely be the most direct way to improve understanding of mechanisms and inform management.

Journal of Wildlife Management↗

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↗

USGS2017

Density-dependent vulnerability of forest ecosystems to drought

1. Climate models predict increasing drought intensity and frequency for many regions, which may have negative consequences for tree recruitment, growth and mortality, as well as forest ecosystem services. Furthermore, practical strategies for minimizing vulnerability to drought are limited. Tree population density, a metric of tree abundance in a given area, is a primary driver of competitive intensity among trees, which influences tree growth and mortality. Manipulating tree population density may be a mechanism for moderating drought-induced stress and growth reductions, although the relationship between tree population density and tree drought vulnerability remains poorly quantified, especially across climatic gradients. 2. In this study, we examined three long-term forest ecosystem experiments in two widely distributed North American pine species, ponderosa pine Pinus ponderosa (Lawson & C. Lawson) and red pine Pinus resinosa (Aiton), to better elucidate the relationship between tree population density, growth and drought. These experiments span a broad latitude and aridity range and include tree population density treatments that have been purposefully maintained for several decades. We investigated how tree population density influenced resistance (growth during drought) and resilience (growth after drought compared to pre-drought growth) of stand-level growth during and after documented drought events. 3. Our results show that relative tree population density was negatively related to drought resistance and resilience, indicating that trees growing at lower densities were less vulnerable to drought. This result was apparent in all three forest ecosystems, and was consistent across species, stand age and drought intensity. 4. Synthesis and applications . Our results highlighted that managing pine forest ecosystems at low tree population density represents a promising adaptive strategy for reducing the adverse impacts of drought on forest growth in coming decades. Nonetheless, the broader applicability of our findings to other types of forest ecosystems merits additional investigation.

Journal of Applied Ecology↗

USGS1995

Graphical and PC-software analysis of volcano eruption precursors according to the Materials Failure Forecast Method (FFM)

The Materials Failure Forecasting Method for volcanic eruptions (FFM) analyses the rate of precursory phenomena. Time of eruption onset is derived from the time of “failure” implied by accelerating rate of deformation. The approach attempts to fit data, Ω, to the differential relationship Ω Ω ¨ = A Ω ˙ "> ¨ = AΩ˙ , where the dot superscript represents the time derivative, and the data Ω may be any of several parameters describing the accelerating deformation or energy release of the volcanic system. Rate coefficients, A and α, may be derived from appropriate data sets to provide an estimate of time to “failure”. As the method is still an experimental technique, it should be used with appropriate judgment during times of volcanic crisis. Limitations of the approach are identified and discussed. Several kinds of eruption precursory phenomena, all simulating accelerating creep during the mechanical deformation of the system, can be used with FFM. Among these are tilt data, slope-distance measurements, crater fault movements and seismicity. The use of seismic coda, seismic amplitude-derived energy release and time-integrated amplitudes or coda lengths are examined. Usage of cumulative coda length directly has some practical advantages to more rigorously derived parameters, and RSAM and SSAM technologies appear to be well suited to real-time applications. One graphical and four numerical techniques of applying FFM are discussed. The graphical technique is based on an inverse representation of rate versus time. For α = 2, the inverse rate plot is linear; it is concave upward for α < 2 and concave downward for α > 2. The eruption time is found by simple extrapolation of the data set toward the time axis. Three numerical techniques are based on linear least-squares fits to linearized data sets. The “linearized least-squares technique” is most robust and is expected to be the most practical numerical technique. This technique is based on an iterative linearization of the given rate-time series. The hindsight technique is disadvantaged by a bias favouring a too early eruption time in foresight applications. The “log rate versus log acceleration technique”, utilizing a logarithmic representation of the fundamental differential equation, is disadvantaged by large data scatter after interpolation of accelerations. One further numerical technique, a nonlinear least-squares fit to rate data, requires special and more complex software. PC-oriented computer codes were developed for data manipulation, application of the three linearizing numerical methods, and curve fitting. Separate software is required for graphing purposes. All three linearizing techniques facilitate an eruption window based on a data envelope according to the linear least-squares fit, at a specific level of confidence, and an estimated rate at time of failure.

Journal of Volcanology and Geothermal Research↗

USGS1997

A numerical investigation of choked flow dynamics and its application to the triggering mechanism of long-period events at Redoubt Volcano, Alaska

We use numerical simulations of transonic flow through a crack to study the dynamics of the formation of shock waves downstream from a nozzle-like constriction inside the crack. The model solves the full set of Navier-Stokes equations in two dimensions via an explicit multifield finite difference representation. The crack walls are assumed to be perfectly rigid, and elastic coupling to the solid is not considered. The simulations demonstrate how the behavior of unsteady shock waves near the walls can produce recurring step-like pressure transients in the flow, which in turn induce resonance of the fluid-filled crack. The motion of the shock waves is governed primarily by smooth, low-amplitude pressure fluctuations at the outlet of the crack. The force induced on the walls scales with the amplitude of the shock, which is a function of the magnitude of the inlet pressure, aperture of the constriction, and thickness of the boundary layer. The applied force also scales in proportion to the spatial extent of the shock excursion, which depends on the fluctuation rate of outlet pressure. Using the source parameters of long-period (LP) events at Redoubt Volcano, Alaska, as a guide for our simulations, we infer that coupling of the shock to the walls occurs for crack inlet to outlet pressure ratios p i /p o >2.31 and that the position of the shock front becomes most sensitive to outlet pressure fluctuations for flow regimes with p i /p o >2.48. For such regimes, fluctuations of outlet pressure of up to ±0.5 MPa at rates up to 3 MPa/s are sufficient to induce pressure transients with magnitudes up to 12.5 MPa over 0.1–2.5 m of the walls within ∼0.5 s. These flow parameters may be adequate for triggering the LP events in the precursory swarm to the December 14, 1989, eruption of Redoubt. According to the flow model the recurrence rate and amplitudes of L.P events are inferred to be a manifestation of the response of a shallow hydrothermal reservoir to the sustained injection of superheated steam from a magma column roofing below this reservoir.

Journal of Geophysical Research B: Solid Earth↗

USGS2003

Compensatory mechanisms in Great Lakes sea lamprey populations: implications for alternative control strategies

Compensatory mechanisms are demographic processes that tend to increase population growth rates at lower population density. These processes will tend to reduce the effectiveness of actions that use controls on reproductive success to suppress sea lamprey ( Petromyzon marinus ), an economically important pest in the Great Lakes. Historical evidence for compensatory mechanisms in sea lamprey populations was reviewed, and revealed: (1) strong evidence for shifts in sex ratios as sea lamprey abundance was reduced in the early years of the control program; (2) weak and equivocal evidence for increased growth rates of sea lamprey cohorts re-colonizing streams following a lampricide treatment; and (3) suggestions of other compensatory processes, such as earlier ages at metamorphosis, but with little empirical evidence. Larval size distribution data for cohorts in the first and second years following a lampricide treatment (26 pairs of cohorts in 20 streams) was analyzed and did not indicate a consistent pattern of more rapid growth of the first colonizing cohort (only 11 of 33 cases). To test for compensation between spawning and age-1 in sea lamprey populations, data were analyzed for 49 stream-years for which spawning female abundance was known and age-1 abundance was estimated in the following year. A fit of these data to a Ricker stock-recruitment function showed evidence for compensation, measured as reduced survival to age 1 at higher abundance of spawning females. More obvious, however, was a large amount of density-independent variation in survival, which tends to mask evidence for compensatory survival. The results were applied to a simple model that simulates sea lamprey populations and their control in a hypothetical lake. Control strategies that targeted reproductive success performed far less well than comparable strategies that targeted larval populations, because density-independent recruitment variation leads to occasional strong year classes even when spawner abundance is reduced to low levels through alternative control. It is concluded that further study of recruitment variation in lamprey populations is critical to rationalizing alternative controls that target reproductive success, and that recruitment variation needs to be incorporated into models used to evaluate sea lamprey control options.

Journal of Great Lakes Research↗

USGS2025

Confocal laser-scanning microscopy (CLSM)-based thermal maturity of Tasmanites and progress in standardization of fluorescence microspectrometry

Evaluation of thermal maturity in vitrinite-free or vitrinite-deficient sediments via fluorescence microspectrometry can provide relevant information related to petroleum exploration and thermal history assessment. However, variation in spectral fluorescence properties of alginite macerals with increasing thermal maturity is largely underexplored. Here, authors of this study have applied confocal laser-scanning microscopy (CLSM) in conjunction with fluorescence microspectrometry to a maturity series of marine Upper Devonian Tasmanites algae from the Ohio Shale (Huron Member) and a single sample from the Marcellus Formation of the Appalachian Basin. Spectral fluorescence properties of Tasmanites were evaluated in relation to orientation, measurement location, and the number of measurements per sample, and were compared to published literature. Emission spectra of Tasmanites from continuous wave laser excitation (405 nm) were acquired from sections perpendicular and parallel to bedding and at different positions within individual Tasmanites bodies. The results showed a progressive red-shift in emission maxima (λ max ) in a large sample sized maturity series ( N = 19), e.g., 493 to 578 nm for the perpendicular section at middle position. Further, blue-shifted apex and mineral-adjacent positions within sections perpendicular to bedding were observed, with the latter being reported here for the first time. While blue-shift at apex positions can be attributed to mechanical deformation-induced reorientation of photoselected fluorophores, the blue-shifted mineral-adjacent positions could result from strain loading and development of a plastic deformation region at the mineral contact zone with Tasmanites . A decrease in standard deviation with increasing number of measured emission maxima is well-observed, and 15 to 20 individual measurements per sample appears sufficient for low standard deviation and coefficient of variance. CLSM-derived thermal maturity parameters indicated that a moderate positive correlation of red/green quotient (Q 650/500 ; R 2 = 0.67) with solid bitumen reflectance (BR o in %) exists. For reproducible results, the determination of λ max and Q 650/500 should be conducted exclusively in the middle position at perpendicular and parallel sections of the polished whole-rock pellets, where the lowest standard deviation in measurement was observed. These results strengthen the suitability and relevance of the CLSM technique in thermal maturity studies of dispersed organic matter (DOM) and contribute to the standardization of fluorescence microspectrometry methods in organic petrology investigation.

Kentucky, Ohio, Virginia, West Virginia↗

USGS2013

Tritium plume dynamics in the shallow unsaturated zone in an arid environment

The spatiotemporal variability of a tritium plume in the shallow unsaturated zone and the mechanisms controlling its transport were evaluated during a 10-yr study. Plume movement was minimal and its mass declined by 68%. Upward-directed diffusive-vapor tritium fluxes and radioactive decay accounted for most of the observed plume-mass declines. Effective isolation of tritium ( 3 H) and other contaminants at waste-burial facilities requires improved understanding of transport processes and pathways. Previous studies documented an anomalously widespread (i.e., theoretically unexpected) distribution of 3 H (>400 m from burial trenches) in a dry, sub-root-zone gravelly layer (1–2-m depth) adjacent to a low-level radioactive waste (LLRW) burial facility in the Amargosa Desert, Nevada, that closed in 1992. The objectives of this study were to: (i) characterize long-term, spatiotemporal variability of 3 H plumes; and (ii) quantify the processes controlling 3 H behavior in the sub-root-zone gravelly layer beneath native vegetation adjacent to the facility. Geostatistical methods, spatial moment analyses, and mass flux calculations were applied to a spatiotemporally comprehensive, 10-yr data set (2001–2011). Results showed minimal bulk-plume advancement during the study period and limited Fickian spreading of mass. Observed spreading rates were generally consistent with theoretical vapor-phase dispersion. The plume mass diminished more rapidly than would be expected from radioactive decay alone, indicating net efflux from the plume. Estimates of upward 3 H efflux via diffusive-vapor movement were >10× greater than by dispersive-vapor or total-liquid movement. Total vertical fluxes were >20× greater than lateral diffusive-vapor fluxes, highlighting the importance of upward migration toward the land surface. Mass-balance calculations showed that radioactive decay and upward diffusive-vapor fluxes contributed the majority of plume loss. Results indicate that plume losses substantially exceeded any continuing 3 H contribution to the plume from the LLRW facility during 2001 to 2011 and suggest that the widespread 3 H distribution resulted from transport before 2001.

Nevada↗

USGS2014

Seismicity and seismic structure at Okmok Volcano, Alaska

Okmok volcano is an active volcanic caldera located on the northeastern portion of Umnak Island in the Aleutian arc, with recent eruptions in 1997 and 2008. The Okmok area had ~900 locatable earthquakes between 2003 and June 2008, and an additional ~600 earthquakes from the beginning of the 2008 eruption to mid 2009, providing an adequate dataset for seismic tomography. To image the seismic velocity structure of Okmok, we apply waveform cross-correlation using bispectrum verification and double-difference tomography to a subset of these earthquakes. We also perform P-wave attenuation tomography using a spectral decay technique. We examine the spatio-temporal characteristics of seismicity in the opening sequence of the 2008 eruption to investigate the path of magma migration during the establishment of a new eruptive vent. We also incorporate the new earthquake relocations and three-dimensional (3D) velocity model with first-motion polarities to compute focal mechanisms for selected events in the 2008 pre-eruptive and eruptive periods. Through these techniques we obtain precise relocations, a well-constrained 3D P-wave velocity model, and a marginally resolved S-wave velocity model. We image a main low Vp and Vs anomaly directly under the caldera consisting of a shallow zone at 0–2 km depth connected to a larger deeper zone that extends to about 6 km depth. We find that areas of low Qp are concentrated in the central to southwestern portion of the caldera and correspond fairly well with areas of low Vp. We interpret the deeper part of the low velocity anomaly (4–6 km depth) beneath the caldera as a magma body. This is consistent with results from ambient noise tomography and suggests that previous estimates of depth to Okmok's magma chamber based only on geodetic data may be too shallow. The distribution of events preceding the 2008 eruption suggest that a combination of overpressure in the zone surrounding the magma chamber and the introduction of new material from below were jointly responsible for the explosive eruption. Magma escaping from the top of the main magma chamber likely reacted with both a smaller shallow pod of magma and groundwater on its way up below the Cone D area. The earthquakes in the 2008 pre-eruptive and eruptive periods are found to have a mixture of strike-slip, oblique normal, and oblique thrust mechanisms, with a dominant P-axis orientation that is nearly perpendicular to the regional tectonic stress. This may indicate that the stresses related to magmatic activity locally dominated regional tectonic forces during this time period.

Alaska↗

USGS2015

Seismic source dynamics of gas-piston activity at Kı̄lauea Volcano, Hawai‘i

Since 2008, eruptive activity at the summit of Kı̄lauea Volcano, Hawai‘i has been confined to the new Overlook pit crater within the Halema‘uma‘u Crater. Among the broad range of magmatic processes observed in the new pit are recurring episodes of gas pistoning. The gas-piston activity is accompanied by seismic signals that are recorded by a broadband network deployed in the summit caldera. We use raw data recorded with this network to model the source mechanism of representative gas-piston events in a sequence that occurred on 20–25 August 2011 during a gentle inflation of the Kı̄lauea summit. 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–10,000 s. Most of the seismic wavefield produced by gas-pistoning originates in a source region ∼1 km below the eastern perimeter of the Halema‘uma‘u pit crater. The observed waveforms are well explained by a simple volumetric source with geometry composed of two intersecting cracks featuring an east striking crack (dike) dipping 80°to the north, intersecting a north striking crack (another dike) dipping 65° to the east. Each gas-piston event is marked by a similar rapid inflation lasting a few minutes, trailed by a slower deflation ramp extending up to 15 min, attributed to the efficient coupling at the source centroid location of the pressure and momentum changes accompanying the growth and collapse of a layer of foam at 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 : foam thickness (10–50 m), foam cell diameter (0.04–0.10 m), and gas-injection velocity (0.01–0.06 m s −1 ). Gas-piston activity occurs in a narrow pipe with diameter of 6 m connecting the Halema‘uma‘u pit crater to the subjacent dike system. The height of the magma column is estimated at ∼104 m at the start of the sequence based on the period of very long period (VLP) oscillations accompanying the onset of the gas-piston signal. Based on the change in the period of VLP oscillations and tilt evidence, the height of the magma column is inferred to have risen by up to ∼23 m by the end of the 5 day long sequence. A penny-shaped crack model of the dike geometry yields effective diameters of ∼1.2–2.9 km for the east dike and 0.7 km for the north dike. The shallower north dike segment is embedded in a relatively weak medium, compatible with expected mechanical properties in the hydrothermal environment of this dike.

Hawaii↗

USGS2016

RMT focal plane sensitivity to seismic network geometry and faulting style

Modern tectonic studies often use regional moment tensors (RMTs) to interpret the seismotectonic framework of an earthquake or earthquake sequence; however, despite extensive use, little existing work addresses RMT parameter uncertainty. Here, we quantify how network geometry and faulting style affect RMT sensitivity. We examine how data-model fits change with fault plane geometry (strike and dip) for varying station configurations. We calculate the relative data fit for incrementally varying geometries about a best-fitting solution, applying our workflow to real and synthetic seismograms for both real and hypothetical station distributions and earthquakes. Initially, we conduct purely observational tests, computing RMTs from synthetic seismograms for hypothetical earthquakes and a series of well-behaved network geometries. We then incorporate real data and station distributions from the International Maule Aftershock Deployment (IMAD), which recorded aftershocks of the 2010 M W 8.8 Maule earthquake, and a set of regional stations capturing the ongoing earthquake sequence in Oklahoma and southern Kansas. We consider RMTs computed under three scenarios: (1) real seismic records selected for high data quality; (2) synthetic seismic records with noise computed for the observed source-station pairings and (3) synthetic seismic records with noise computed for all possible station-source pairings. To assess RMT sensitivity for each test, we observe the ‘fit falloff’, which portrays how relative fit changes when strike or dip varies incrementally; we then derive the ranges of acceptable strikes and dips by identifying the span of solutions with relative fits larger than 90 per cent of the best fit. For the azimuthally incomplete IMAD network, Scenario 3 best constrains fault geometry, with average ranges of 45° and 31° for strike and dip, respectively. In Oklahoma, Scenario 3 best constrains fault dip with an average range of 46°; however, strike is best constrained by Scenario 1, with a range of 26°. We draw two main conclusions from this study. (1) Station distribution impacts our ability to constrain RMTs using waveform time-series; however, in some tectonic settings, faulting style also plays a significant role and (2) increasing station density and data quantity (both the number of stations and the number of individual channels) does not necessarily improve RMT constraint. These results may be useful when organizing future seismic deployments (e.g. by concentrating stations in alignment with anticipated nodal planes), and in computing RMTs, either by guiding a more rigorous data selection process for input data or informing variable weighting among the selected data (e.g. by eliminating the transverse component when strike-slip mechanisms are expected).

Geophysical Journal International↗

USGS2021

The influence of land cover and storm magnitude on hydrologic flowpath activation and runoff generation in steep tropical catchments of central Panama

Despite abundant research documenting that land use/land cover (LULC) have substantial impacts on the hydrology of humid tropical systems, field-based evidence for the physical mechanisms behind these impacts are still lacking. In particular, our understanding of the hydrologic flowpaths that generate runoff in these systems, and how they vary with respect to LULC is insufficient to inform both physically-based hydrologic modeling and land-use decision-making. In this study, we use end-member mixing analysis (EMMA) of stream chemistry, and hydrometric characterizations of hillslope soil moisture to identify hydrologic flowpaths in humid tropical steep-land catchments of varying LULC: mature tropical forest, young secondary tropical forest, cattle pasture. EMMA was applied to data from 14 storm events (six at the mature forest, five at the young secondary forest, and three at the cattle pasture) that were intensively sampled during the 2017 wet season representing a wide range of rainfall magnitudes and intensities. Additionally, volumetric-soil-moisture responses at multiple depths were characterized during and after 74 storm events occurring from 2015 to 2017. EMMA results indicated that lateral preferential flow within the top 30 cm of the soil profile was a dominant source of runoff generation at the two forested catchments, with the contribution of this flow path increasing with rainfall magnitude and intensity. This was corroborated by volumetric-soil-moisture data, that showed that a perched zone of saturation developed at 30 cm at the time of peak storm runoff during the largest events and lasted for the remaining duration of the event. EMMA indicated that runoff was a combination of infiltration-excess overland flow and lateral subsurface flow in the actively grazed pastoral catchment. There, overland flow contributed 62 % of runoff during the highest runoff rate sampled (35.3 mm/hr) and this contribution increased substantially with storm magnitude. This flowpath identification was also supported by volumetric-soil-moisture data at the pasture, with peak saturation at all depths during the largest storm events occurring up to 30 min after peak runoff. These results provide a mechanistic explanation for previously observed hydrologic differences among tropical LULCs. Additionally, the wide range of hydrologic conditions during these storm events provide a basis for understanding how future changes to this, and similar humid tropical regions will impact hydrological processes and water availability.

Journal of Hydrology↗

USGS2025

Demographic mechanisms of snowshoe hare population cycles in Yukon, Canada

One hundred years have elapsed since Charles Elton (1924) described the periodic fluctuations in North American snowshoe hare abundance, yet mechanisms underlying 9–11-year population cycles in snowshoe hares continue to be debated. We applied multistate capture–mark–recapture models to long-term field data (1977–2020) based on >20,000 captures of >7000 unique snowshoe hares ( Lepus americanus ) from Kluane Lake, Yukon, Canada, to estimate and model state-specific demographic parameters. Juveniles had the lowest and reproductive adult females the highest apparent survival. Apparent survival of all sex-age classes was highest during the mid- and late-breeding seasons and was generally better during the increase phase. Conditional probability of females transitioning from non-reproductive to reproductive state, and reproductive females remaining in the reproductive state, increased substantially as the population transitioned from low to increase phase throughout the breeding season. Analysis of stage-structured matrix population models revealed that population-dynamic characteristics were strongly phase-specific, and also varied across seasons, with the increase phases being characterized by high monthly asymptotic population growth rate. Snowshoe hares experienced short stage-specific generation time during the early breeding season across all phases; they experienced relatively long generation time during the increase and low phase of the mid-breeding season, and the increase and peak phase of the late breeding season. Elasticity analyses showed that asymptotic population growth rate was proportionately most sensitive to changes in survival of adult females across all phases and seasons. However, retrospective life table response experiment analysis showed that rapid growth of the snowshoe hare populations during the increase phase was due to improvements in reproductive transitions and pre-weaning survival, whereas population declines are caused primarily by reduced survival (primarily, pre-weaning survival), with reduced reproductive transitions and smaller litter sizes playing a secondary role. Our results suggest that cyclic populations of snowshoe hares are characterized by complex demographic and population-dynamic patterns, depending on phase of the cycle and reproductive season, and that different demographic mechanisms underlie rapid population growth during the increase phase, and swift population declines as the population transitions from the peak to the decline phase. Because our study represents the first comprehensive demographic and population-dynamic study of a cyclic population, similar studies would be needed to test the generalities of our conclusions. Whereas density-dependent predation has been shown to be the primary cause of phase-related changes in survival, future research should focus on identifying mechanisms underlying phase-related changes in reproductive parameters.

Yukon↗

USGS2012

Resonance of a fluid-driven crack: Radiation properties and implications for the source of long-period events and harmonic tremor

A dynamic source model is presented, in which a three-dimensional crack containing a viscous compressible fluid is excited into resonance by an impulsive pressure transient applied over a small area ΔS of the crack surface. The crack excitation depends critically on two dimensionless parameters called the crack stiffness, C = ( b /μ)( L / d ), and viscous damping loss, F = (12η L )/(ρ ƒ d 2 α), where b is the bulk modulus, η is the viscosity, ρ ƒ is the density of the fluid, μ is the rigidity, α is the compressional velocity of the solid, L is the crack length, and d is the crack thickness. The first parameter characterizes the ability of the crack to vibrate and shapes the spectral signature of the source, and the second quantifies the effect of fluid viscosity on the duration of resonance. Resonance is sustained by a very slow wave trapped in the fluid-filled crack. This guided wave, called the crack wave, is similar to the tube wave propagating in a fluid-filled borehole; it is inversely dispersive, showing a phase velocity that decreases with increasing wavelength, and its wave speed is always lower than the acoustic velocity of the fluid, decreasing rapidly as the crack stiffness increases. The source spectrum shows many sharp peaks characterizing the individual modes of vibration of the crack; the variation of spectral shape, both in the number and width of peaks, is surprisingly complex, reflecting the interference between the lateral and longitudinal modes of resonance, as well as nodes for these modes. The far-field spectrum is marked by narrow-band dominant and subdominant peaks that reflect the interaction of the various source modes. The frequency of the dominant spectral peak radiated by the source is independent of the radiation direction. The frequency, bandwidth, and spacing of the resonant peaks are strongly dependent on the crack stiffness, larger values of the stiffness factor shifting these peaks to lower frequencies and decreasing their bandwidth. The excitation of a particular mode depends on the position of the trigger and on the extent of the crack surface affected by the pressure transient. Fluid viscosity decreases the amplitudes of the main spectral peaks, smears out the finer structure of the spectrum, and greatly reduces the duration of the radiated signal. The energy loss by radiation is stronger for high frequencies, producing a seismic signature that is marked by a high-frequency content near the onset of the signal and dominated by a longer-period component of much longer duration in the signal coda. Such signature is in harmony with those displayed by long-period events observed on active volcanoes and in hydrofracture experiments. The very low velocity which is possible in a crack with high stiffness ( C ≥ 100) also provides an attractive explanation for very long period tremor, such as type 2 tremor at Aso volcano, Japan, without the requirement of an unrealistically large magma container. The standing wave pattern set up on the crack surface by the sustained resonance in the fluid is observable in the near field of the crack, suggesting that the location and extent of the source may be estimated from the mapping of the pattern of nodes and antinodes seen in its vicinity. According to the model, the long-period event and harmonic tremor share the same source but differ in the boundary conditions for fluid flow and in the triggering mechanism setting up the resonance of the source, the former being viewed as the impulse response of the tremor generating system and the latter representing the excitation due to more complex forcing functions.

Journal of Geophysical Research Solid Earth↗