Search USGS⌕ Search

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.

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

At least 109 records · Page 6Linked to original sources
USGS2003

Source mechanism of long-period events at Kusatsu-Shirane Volcano, Japan, inferred from waveform inversion of the effective excitation functions

We investigate the source mechanism of long-period (LP) events observed at Kusatsu-Shirane Volcano, Japan, based on waveform inversions of their effective excitation functions. The effective excitation function, which represents the apparent excitation observed at individual receivers, is estimated by applying an autoregressive filter to the LP waveform. Assuming a point source, we apply this method to seven LP events the waveforms of which are characterized by simple decaying and nearly monochromatic oscillations with frequency in the range 1-3 Hz. The results of the waveform inversions show dominant volumetric change components accompanied by single force components, common to all the events analyzed, and suggesting a repeated activation of a sub-horizontal crack located 300 m beneath the summit crater lakes. Based on these results, we propose a model of the source process of LP seismicity, in which a gradual buildup of steam pressure in a hydrothermal crack in response to magmatic heat causes repeated discharges of steam from the crack. The rapid discharge of fluid causes the collapse of the fluid-filled crack and excites acoustic oscillations of the crack, which produce the characteristic waveforms observed in the LP events. The presence of a single force synchronous with the collapse of the crack is interpreted as the release of gravitational energy that occurs as the slug of steam ejected from the crack ascends toward the surface and is replaced by cooler water flowing downward in a fluid-filled conduit linking the crack and the base of the crater lake. ?? 2003 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗
USGS2012

Patterns and controlling factors of species diversity in the Arctic Ocean

Aim The Arctic Ocean is one of the last near-pristine regions on Earth, and, although human activities are expected to impact on Arctic ecosystems, we know very little about baseline patterns of Arctic Ocean biodiversity. This paper aims to describe Arctic Ocean-wide patterns of benthic biodiversity and to explore factors related to the large-scale species diversity patterns. Location Arctic Ocean. Methods We used large ostracode and foraminiferal datasets to describe the biodiversity patterns and applied comprehensive ecological modelling to test the degree to which these patterns are potentially governed by environmental factors, such as temperature, productivity, seasonality, ice cover and others. To test environmental control of the observed diversity patterns, subsets of samples for which all environmental parameters were available were analysed with multiple regression and model averaging. Results Well-known negative latitudinal species diversity gradients (LSDGs) were found in metazoan Ostracoda, but the LSDGs were unimodal with an intermediate maximum with respect to latitude in protozoan foraminifera. Depth species diversity gradients were unimodal, with peaks in diversity shallower than those in other oceans. Our modelling results showed that several factors are significant predictors of diversity, but the significant predictors were different among shallow marine ostracodes, deep-sea ostracodes and deep-sea foraminifera. Main conclusions On the basis of these Arctic Ocean-wide comprehensive datasets, we document large-scale diversity patterns with respect to latitude and depth. Our modelling results suggest that the underlying mechanisms causing these species diversity patterns are unexpectedly complex. The environmental parameters of temperature, surface productivity, seasonality of productivity, salinity and ice cover can all play a role in shaping large-scale diversity patterns, but their relative importance may depend on the ecological preferences of taxa and the oceanographic context of regions. These results suggest that a multiplicity of variables appear to be related to community structure in this system.

Journal of Biogeography↗
USGS2009

Hydraulic anisotropy characterization of pneumatic-fractured sediments using azimuthal self potential gradient

The pneumatic fracturing technique is used to enhance the permeability and porosity of tight unconsolidated soils (e.g. clays), thereby improving the effectiveness of remediation treatments. Azimuthal self potential gradient (ASPG) surveys were performed on a compacted, unconsolidated clay block in order to evaluate their potential to delineate contaminant migration pathways in a mechanically-induced fracture network. Azimuthal resistivity (ARS) measurements were also made for comparative purposes. Following similar procedures to those used in the field, compressed kaolinite sediments were pneumatically fractured and the resulting fracture geometry characterized from strike analysis of visible fractures combined with strike data from optical borehole televiewer (BHTV) imaging. We subsequently injected a simulated treatment (electrolyte/dye) into the fractures. Both ASPG and ARS data exhibit anisotropic geoelectric signatures resulting from the fracturing. Self potentials observed during injection of electrolyte are consistent with electrokinetic theory and previous laboratory results on a fracture block model. Visual (polar plot) analysis and linear regression of cross plots show ASPG lobes are correlated with azimuths of high fracture strike density, evidence that the ASPG anisotropy is a proxy measure of hydraulic anisotropy created by the pneumatic fracturing. However, ARS data are uncorrelated with fracture strike maxima and resistivity anisotropy is probably dominated by enhanced surface conduction along azimuths of weak 'starter paths' formed from pulverization of the clay and increases in interfacial surface area. We find the magnitude of electrokinetic SP scales with the applied N2 gas pressure gradient (??PN2) for any particular hydraulically-active fracture set and that the positive lobe of the ASPG anomaly indicates the flow direction within the fracture network. These findings demonstrate the use of ASPG in characterizing the effectiveness of (1) pneumatic fracturing and (2) defining likely flow directions of remedial treatments in unconsolidated sediments and rock. ?? 2008 Elsevier B.V. All rights reserved.

Journal of Contaminant Hydrology↗
USGS2020

Plastic faulting in ice

Plastic faulting is a brittle‐like failure phenomenon exhibited by water ice and several other rock types under confinement. It is suspected to be the mechanism of deep earthquakes and extreme cases of shear localization in shallow rocks. Unlike ordinary Coulombic failure, plastic faulting is characterized by a pressure‐independent failure strength and fault plane oriented 45° to maximum principal stress. To research the question of how the instability initiates, we conducted over 50 constant‐displacement‐rate experiments on polycrystalline ice (phases Ih and II) near the brittle‐to‐ductile (B‐D) transition, at confining pressures P = 0–300 MPa, applied strain rates = 5 × 10 −5 – 7 × 10 −3 s −1 , temperatures T = 105–233 K, and mean grain sizes d = 0.25–1.18 mm. We find that (1) the width of the B‐D transition in variable space is vanishingly narrow, to the point of appearing as a crossover, (2) a plastic fault plane, once formed, is not a zone of subsequent weakness, (3) distributed ice I→II phase transformation in small amounts (<1 vol%) shows no causal relationship to subsequent failure, and (4) plastic faulting also occurs in ice II. We hypothesize that the elusive nucleating “trigger” parallels that of metals and ceramics undergoing severe plastic deformation, wherein transient local structural rearrangement occurs, in turn causing material strength to drop to a level sufficiently low, in a volume sufficiently large, that adiabatic instability is nucleated. Our results do not require and often are inconsistent with phase transformation. Plastic faulting may therefore be available to all solids undergoing severe deformation, and its appearance in so few is simply the result of insufficiently extreme conditions.

Journal of Geophysical Research- Solid Earth↗
USGS2003

Analysis of aquifer mineralization by paleodrainage channels

Mineralization of groundwater resources is a problem in south-central Kansas, due to the penetration of saline water from Permian bedrock formations into the overlying alluvial aquifer. One of the mechanisms involved in the mineralization involves small bedrock features of high permeability located in places occupied by streams and rivers in past geological eras. These geological features are termed 'paleodrainage channels'. The permeability of the overlying aquifer can be significantly smaller than that of the channel fill material. The comparatively fast migration of saline water through these channels of high permeability is associated with the transfer of minerals into the overlying freshwater aquifer. This study applies a set of boundary layer approaches to quantify the process of mineral transfer from the channels into the aquifer. The methods used in the present study provide quick estimation and evaluation of the dilution of the channel flow, as well as mineral concentration profile changes in the mineralized zone created in the overlying aquifer. More generally, the method can also be useful for the analysis and evaluation of various types of groundwater contamination in heterogeneous aquifers. The application of the method is exemplified by a complete set of calculations characterizing the possible mineralization process at a specific channel in south central Kansas. Sensitivity analyses are performed and provide information about the importance of the various parameters that affect the mineralization process. Some possible scenarios for the aquifer mineralization phenomena are described and evaluated. It is shown that the channel mineralization may create either several stream tubes of the aquifer with high mineral concentration, or many stream tubes mineralized to a lesser extent. Characteristics of these two patterns of aquifer mineralization are quantified and discussed. ?? 2003 Published by Elsevier Science B.V.

Journal of Hydrology↗
USGS2026

Benchmark dataset of historical annual peak floods classified by causal mechanisms for select US river basins

Considering the causal mechanisms of floods can improve estimates of flood recurrence intervals given that certain flood types can be associated with higher magnitude and more damaging floods. However, few verified datasets of flood types are available to validate the semiautomated and automated classification algorithms needed to apply flood-typing across large hydrologically diverse regions. To address this gap, a benchmark dataset of manually classified flood types was compiled for 1,763 annual maximum flood peaks from 18 stream gauges in six different river basins across the conterminous United States from 1851 to 2022. Within each basin, three representative stream gauges were selected for manual flood typing. A flexible classification framework is introduced that facilitates flood typing across hydrologically diverse regions and accommodates unique combinations of weather and antecedent watershed conditions specific to each region. Floods were manually typed by domain experts using multiple lines of evidence to identify a primary surface water input of each flood (rainfall, snowmelt, or both) and, if relevant, associated storm type and secondary causal mechanisms characterizing antecedent watershed conditions. Across all the study basins, 49% of historical annual maximum flood peaks were attributable to rainfall, 28% to snowmelt, 22% to mixed precipitation, and 1% could not be assigned to a mechanism due to missing or incomplete data. The proposed flood-typing schema supports varying levels of flood typing specificity required for mixed population flood-frequency analysis, flood-type-specific design hydrographs, water quality response studies, and additional applications. This detailed, manually determined benchmark dataset serves as a resource that can be used developing and validating automated or machine learning-based algorithms capable of operationalizing expanded flood peak information.

conterminous United States↗
USGS2021

Developing species-age cohorts from forest inventory and analysis data to parameterize a forest landscape model

Simulating long-term, landscape level changes in forest composition requires estimates of stand age to initialize succession models. Detailed stand ages are rarely available, and even general information on stand history often is lacking. We used data from USDA Forest Service Forest Inventory and Analysis (FIA) database to estimate broad age classes for a forested landscape to simulate changes in landscape composition and structure relative to climate change at Fort Drum, a 43,000 ha U.S. Army installation in northwestern New York. Using simple linear regression, we developed relationships between tree diameter and age for FIA site trees from the host and adjacent ecoregions and applied those relationships to forest stands at Fort Drum. We observed that approximately half of the variation in age was explained by diameter breast height (DBH) across all species studied ( r 2 = 0.42 for sugar maple Acer saccharum to 0.63 for white ash Fraxinus americana ). We then used age-diameter relationships from published research on northern hardwood species to calibrate results from the FIA-based analysis. With predicted stand age, we used tree species life histories and environmental conditions represented by ecological site types to parameterize a stochastic forest landscape model (LANDIS-II) to spatially and temporally model successional changes in forest communities at Fort Drum. Forest stands modeled over 100 years without significant disturbance appeared to reflect expected patterns of increasing dominance by shade-tolerant mesophytic tree species such as sugar maple, red maple ( Acer rubrum ), and eastern hemlock ( Tsuga canadensis ) where soil moisture was sufficient. On drier sandy soils, eastern white pine ( Pinus strobus ), red pine ( P. resinosa ), northern red oak ( Quercus rubra ), and white oak ( Q. alba ) continued to be important components throughout the modeling period with no net loss at the landscape scale. Our results suggest that despite abundant precipitation and relatively low evapotranspiration rates for the region, low soil water holding capacity and fertility may be limiting factors for the spread of mesophytic species on excessively drained soils in the region. Increasing atmospheric temperatures projected for the region could alter moisture regimes for many coarse-textured soils providing a possible mechanism for expansion of xerophytic tree species.

New York↗
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↗
USGS2024

A ‘how-to’ guide for estimating animal diel activity using hierarchical models

Animal diel activity patterns can aid understanding of (a) how species behaviourally adapt to anthropogenic and natural disturbances, (b) mechanisms of species co-existence through temporal partitioning, and (c) community or ecosystem effects of diel activity shifts. Activity patterns often vary spatially, a feature ignored by the kernel density estimators (KDEs) currently used for estimating diel activity. Ignoring this source of heterogeneity may lead to biased estimates of uncertainty and misleading conclusions regarding the drivers of diel activity. Thus, there is a need for more flexible statistical approaches for estimating activity patterns and testing hypotheses regarding their biotic and abiotic drivers. We illustrate how trigonometric terms and cyclic cubic splines combined with hierarchical models can provide a valuable alternative to KDEs. Like KDEs, these models accommodate circular data, but they can also account for site-to-site and other sources of variability, correlation amongst repeated measures, and variable sampling effort. They can also more readily quantify and test hypotheses related to the effects of covariates on activity patterns. Through empirical case studies, we illustrate how hierarchical models can quantify changes in activity levels due to seasonality and in response to biotic and abiotic factors (e.g. anthropogenic stressors and co-occurrence). We also describe frequentist and Bayesian approaches for quantifying site-specific (conditional) and population-averaged (marginal) activity patterns. We provide guidelines and tutorials with detailed step-by-step instructions for fitting and interpreting hierarchical models applied to time-stamped data, such as those recorded by camera traps and audio recorders. We conclude that this approach offers a viable, flexible, and effective alternative to KDEs when modelling animal activity patterns.

Journal of Animal Ecology↗
USGS2016

Implications of projected climate change for groundwater recharge in the western United States

Existing studies on the impacts of climate change on groundwater recharge are either global or basin/location-specific. The global studies lack the specificity to inform decision making, while the local studies do little to clarify potential changes over large regions (major river basins, states, or groups of states), a scale often important in the development of water policy. An analysis of the potential impact of climate change on groundwater recharge across the western United States (west of 100° longitude) is presented synthesizing existing studies and applying current knowledge of recharge processes and amounts. Eight representative aquifers located across the region were evaluated. For each aquifer published recharge budget components were converted into four standard recharge mechanisms: diffuse, focused, irrigation, and mountain-systems recharge. Future changes in individual recharge mechanisms and total recharge were then estimated for each aquifer. Model-based studies of projected climate-change effects on recharge were available and utilized for half of the aquifers. For the remainder, forecasted changes in temperature and precipitation were logically propagated through each recharge mechanism producing qualitative estimates of direction of changes in recharge only (not magnitude). Several key patterns emerge from the analysis. First, the available estimates indicate average declines of 10–20% in total recharge across the southern aquifers, but with a wide range of uncertainty that includes no change. Second, the northern set of aquifers will likely incur little change to slight increases in total recharge. Third, mountain system recharge is expected to decline across much of the region due to decreased snowpack, with that impact lessening with higher elevation and latitude. Factors contributing the greatest uncertainty in the estimates include: (1) limited studies quantitatively coupling climate projections to recharge estimation methods using detailed, process-based numerical models; (2) a generally poor understanding of hydrologic flowpaths and processes in mountain systems; (3) difficulty predicting the response of focused recharge to potential changes in the frequency and intensity of extreme precipitation events; and (4) unconstrained feedbacks between climate, irrigation practices, and recharge in highly developed aquifer systems.

Journal of Hydrology↗
USGS2014

Formation of fine sediment deposit from a flash flood river in the Mediterranean Sea

We identify the mechanisms controlling fine deposits on the inner-shelf in front of the Besòs River, in the northwestern Mediterranean Sea. This river is characterized by a flash flood regime discharging large amounts of water (more than 20 times the mean water discharge) and sediment in very short periods lasting from hours to few days. Numerical model output was compared with bottom sediment observations and used to characterize the multiple spatial and temporal scales involved in offshore sediment deposit formation. A high-resolution (50 m grid size) coupled hydrodynamic-wave-sediment transport model was applied to the initial stages of the sediment dispersal after a storm-related flood event. After the flood, sediment accumulation was predominantly confined to an area near the coastline as a result of preferential deposition during the final stage of the storm. Subsequent reworking occurred due to wave-induced bottom shear stress that resuspended fine materials, with seaward flow exporting them toward the midshelf. Wave characteristics, sediment availability, and shelf circulation determined the transport after the reworking and the final sediment deposition location. One year simulations of the regional area revealed a prevalent southwestward average flow with increased intensity downstream. The circulation pattern was consistent with the observed fine deposit depocenter being shifted southward from the river mouth. At the southern edge, bathymetry controlled the fine deposition by inducing near-bottom flow convergence enhancing bottom shear stress. According to the short-term and long-term analyses, a seasonal pattern in the fine deposit formation is expected.

Journal of Geophysical Research: Oceans↗
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↗
USGS2013

Source characterization for an explosion during the 2009 eruption of Redoubt Volcano from very-long-period seismic waves

The 2009 eruption of Redoubt produced several very-long-period (VLP) signals associated with explosions. We invert for the source location and mechanism of an explosion at Redoubt volcano using waveform methods applied to broadband recordings. Such characterization of the source carries information on the geometry of the conduit and the physics of the explosion process. Inversions are carried out assuming the volcanic source can be modeled as a point source, with mechanisms described by a) a set of 3 orthogonal forces, b) a moment tensor consisting of force couples, and c) both forces and moment tensor components. We find that the source of the VLP seismic waves during the explosion is well-described by either a combined moment/force source located northeast of the crater and at an elevation of 1.6 km ASL or a moment source at an elevation of 800 m to the southwest of the crater. The moment tensors for the solutions with moment and force and moment-only share similar characteristics. The source time functions for both moment tensors begin with inflation (pressurization) and execute two cycles of deflation-reinflation (depressurization–repressurization). Although the moment/force source provides a better fit to the data, we find that owing to the limited coverage of the broadband stations at Redoubt the moment-only source is the more robust and reliable solution. Based on the moment-only solution, we estimate a volume change of 19,000 m 3 and a pressure change of 7 MPa in a dominant sill and an out-of-phase volume change of 5000 m 3 and pressure change of 1.8 MPa in a subdominant dike at the source location. These results shed new light on the magmatic plumbing system beneath Redoubt and complement previous studies on Vulcanian explosions at other volcanoes.

Alaska↗
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↗
USGS1981

Sequential extraction techniques applied to a porphyry copper deposit in the basin and range province

Samples of minus-80-mesh (<180 μm) stream sediment, rock containing exposed fracture coatings, and jarosite and chrysocolla were collected from an area surrounding the North Silver Bell porphyry Cu deposit near Tucson, Arizona. The samples were subjected to a series of extractions in a scheme originally designed for use on samples from humid or sub-humid environments, in which the following fractions can effectively be separated: (1) carbonates and exchangeable metals; (2) Mn oxides; (3) organic compounds and sulfides; (4) hydrous Fe oxides; and (5) residual crystalline minerals. Jarosite and chrysocolla, two major minerals of the North Silver Bell area, were found to dissolve over two or more steps of the extraction scheme. The results represent only a limited number of samples from one copper deposit. Nevertheless, they do suggest that in a semiarid to arid environment, where mechanical dispersion of such minerals predominates, uncritical assignment of unique phases, such as Mn oxides or organics to a given extraction would lead to false interpretations of weathering processes. However, the relative proportions of elements dissolved in each step of the jarosite and chrysocolla extractions could be used as a “fingerprint” for recognition of the presence of these two minerals in the stream-sediment and rock samples. The relative abundance of hydrous Fe oxide and jarosite and the alteration zoning could be mapped using data from jarosite and chrysocolla extractions. Manganese oxides were also found to have a greater influence on Zn than on Cu or Pb during supergene alteration. The rapid change in relative importance of the first (1 M -acetic acid) extraction for Cu, Zn, and Pb near the mineralized zone suggested the occurrence of minor hydromorphic processes within the stream sediments. Thus, the acetic acid extraction proved the most effective for pinpointing mineralization in sediments. In contrast, the residual fraction had the longest dispersion train, suggesting that total metal concentrations are most effective in arid environments for reconnaissance surveys.

Journal of Geochemical Exploration↗
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↗
USGS2012

The mechanics of ground deformation precursory to dome-building extrusions at Mount St. Helens 1981-1982

Detailed monitoring at Mount St. Helens since 1980 has enabled prediction of the intermittent eruptive activity (mostly dome growth) with unprecedented success. During 1981 and 1982, accelerating deformation of the crater floor around the vent (including radial cracks, thrust faults, and ground tilt) was the earliest indicator of impending activity. Numerical experiments using the finite element method to model the mechanics of the crater floor show that all deformation features can be explained by a uniform shear-stress boundary condition along the conduit wall. The magnitude of the shear stress required to match observed displacements (1–7 MPa) is inversely proportional to the conduit diameter (estimated to be 25–100 m). The most probable source of this shear stress is the flow of viscous magma up the conduit and into the lava dome. We propose a model in which the accelerating deformation, beginning as much as 4 weeks before extrusions, is caused by the increasing velocity of ascending magma in the conduit. This model is examined by using deformation data of the dome before four extrusions in 1981 and 1982 to estimate the volumetric flow rate through the conduit. This flow rate and an estimate of the effective viscosity of the magma enable calculation of an ascent velocity and an applied shear stress that, again, depend on the conduit diameter. The results of these calculations are consistent with the finite element experiments and show that the proposed model is feasible. In light of this model, events observed just before or near the time extrusions began, such as reversals of ground tilt direction from outward to inward and the sudden decrease in the number of shallow earthquakes, may indicate an abrupt decrease of shear stress in the conduit. This could be explained by a decrease in either the ascent velocity, or the effective viscosity of the magma ascending through the shallow conduit, or both, near the time of extrusion. Precursory deformation like that measured at Mount St. Helens should be observable at similar volcanoes elsewhere because it is caused by the fundamental process of magma ascent.

Journal of Geophysical Research Solid Earth↗
USGS2012

Spatial patterns of aftershocks of shallow focus earthquakes in California and implications for deep focus earthquakes

Previous workers have pioneered statistical techniques to study the spatial distribution of aftershocks with respect to the focal mechanism of the main shock. Application of these techniques to deep focus earthquakes failed to show clustering of aftershocks near the nodal planes of the main shocks. To better understand the behavior of these statistics, this study applies them to the aftershocks of six large shallow focus earthquakes in California (August 6, 1979, Coyote Lake; May 2, 1983, Coalinga; April 24, 1984, Morgan Hill; August 4, 1985, Kettleman Hills; July 8, 1986, North Palm Springs; and October 1, 1987, Whittier Narrows). The large number of aftershocks accurately located by dense local networks allows us to treat these aftershock sequences individually instead of combining them, as was done for the deep earthquakes. The results for individual sequences show significant clustering about the closest nodal plane and the strike direction for five of the sequences and about the presumed fault plane for all six sequences. This implies that the previously developed method does work properly. Nonrandom behavior was also found about the slip directions, the P axis, the T axis, and the B axis, but this is probably caused by the lack of independence between these axes and the previously mentioned features of the focal mechanisms. Given that the method does work and that deep aftershocks were not shown to cluster about the main shock nodal planes, the shallow focus data were used to simulate the deep focus study. The goal is to determine if there are artificial factors that make clustering in the deep focus data unobservable. To more closely mimic the work on deep earthquakes, the largest aftershocks from each of the six sequences were combined and studied with respect to their respective main shock focal mechanisms. This reduced the significance of the clustering about the focal mechanism parameters, but not below 95% confidence. Gaussian noise was then added to the aftershock hypocenters in order to determine if the larger hypocentral and focal mechanism errors in the deep focus data could account for the previous negative result. The conclusion is that the following reasons are sufficient to explain the lack of clustering about the main shock nodal planes for the deep focus aftershocks: the need to combine aftershocks from several sequences, the size of the hypocentral location and focal mechanism errors, and the alignment of distant aftershocks with the Wadati-Benioff zone.

Journal of Geophysical Research Solid Earth↗