Search USGSSearch

SEARCH · Search USGS

Results for “Journal of Geophysical Research - Oceans”

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 217 records · Page 12Linked to original sources

A new strategy for earthquake focal mechanisms using waveform-correlation-derived relative polarities and cluster analysis: Application to the 2014 Long Valley Caldera earthquake swarm

In microseismicity analyses, reliable focal mechanisms can typically be obtained for only a small subset of located events. We address this limitation here, presenting a framework for determining robust focal mechanisms for entire populations of very small events. To achieve this, we resolve relative P and S wave polarities between pairs of waveforms by using their signed correlation coefficients—a by-product of previously performed precise earthquake relocation. We then use cluster analysis to group events with similar patterns of polarities across the network. Finally, we apply a standard mechanism inversion to the grouped data, using either catalog or correlation-derived P wave polarity data sets. This approach has great potential for enhancing analyses of spatially concentrated microseismicity such as earthquake swarms, mainshock-aftershock sequences, and industrial reservoir stimulation or injection-induced seismic sequences. To demonstrate its utility, we apply this technique to the 2014 Long Valley Caldera earthquake swarm. In our analysis, 85% of the events (7212 out of 8494 located by Shelly et al . [2016]) fall within five well-constrained mechanism clusters, more than 12 times the number with network-determined mechanisms. Of the earthquakes we characterize, 3023 (42%) have magnitudes smaller than 0.0. We find that mechanism variations are strongly associated with corresponding hypocentral structure, yet mechanism heterogeneity also occurs where it cannot be resolved by hypocentral patterns, often confined to small-magnitude events. Small (5–20°) rotations between mechanism orientations and earthquake location trends persist when we apply 3-D velocity models and might reflect a geometry of en echelon, interlinked shear, and dilational faulting.

California

Ground-rupturing earthquakes on the northern Big Bend of the San Andreas Fault, California, 800 A.D. to Present

Paleoseismic data on the timing of ground-rupturing earthquakes constrain the recurrence behavior of active faults and can provide insight on the rupture history of a fault if earthquakes dated at neighboring sites overlap in age and are considered correlative. This study presents the evidence and ages for 11 earthquakes that occurred along the Big Bend section of the southern San Andreas Fault at the Frazier Mountain paleoseismic site. The most recent earthquake to rupture the site was the M w 7.7–7.9 Fort Tejon earthquake of 1857. We use over 30 trench excavations to document the structural and sedimentological evolution of a small pull-apart basin that has been repeatedly faulted and folded by ground-rupturing earthquakes. A sedimentation rate of 0.4 cm/yr and abundant organic material for radiocarbon dating contribute to a record that is considered complete since 800 A.D. and includes 10 paleoearthquakes. Earthquakes have ruptured this location on average every ~100 years over the last 1200 years, but individual intervals range from ~22 to 186 years. The coefficient of variation of the length of time between earthquakes (0.7) indicates quasiperiodic behavior, similar to other sites along the southern San Andreas Fault. Comparison with the earthquake chronology at neighboring sites along the fault indicates that only one other 1857-size earthquake could have occurred since 1350 A.D., and since 800 A.D., the Big Bend and Mojave sections have ruptured together at most 50% of the time in M w ≥ 7.3 earthquakes.

California

Frictional strength of wet and dry montmorillonite

Montmorillonite is a common mineral in fault zones, and its low strength relative to other common gouge minerals is important in many models of fault rheology. However, the coefficient of friction, μ , varies with degree of saturation and is not well constrained in the literature due to the difficulty of establishing fully drained or fully dried states in the laboratory. We measured μ of both saturated and oven-dried montmorillonite at normal stresses up to 700 MPa. Care was taken to shear saturated samples slowly enough to avoid pore fluid overpressure. For saturated samples, μ increased from 0.10 to 0.28 with applied effective normal stress, while for dry samples μ decreased from 0.78 to 0.45. The steady state rate dependence of friction, ( a − b ), was positive, promoting stable sliding. The wide disparity in reported frictional strengths can be attributed to experimental procedures that promote differing degrees of partial saturation or overpressured pore fluid conditions.

Journal of Geophysical Research

Integrating puffing and explosions in a general scheme for Strombolian-style activity

Strombolian eruptions are among the most common subaerial styles of explosive volcanism worldwide. Distinctive features of each volcano lead to a correspondingly wide range of variations of magnitude and erupted products, but most papers focus on a single type of event at a single volcano. Here, in order to emphasize the common features underlying this diversity of styles, we scrutinize a database from 35 different erupting vents, including 21 thermal infrared videos from Stromboli (Italy), Etna (Italy), Yasur (Vanuatu), and Batu Tara (Indonesia), from puffing, through rapid explosions to normal explosions, with variable ejection parameters and relative abundance of gas, ash, and bombs. Using field observations and high-speed thermal infrared videos processed by a new algorithm, we identify the distinguishing characteristics of each type of activity and how they may relate and interact. In particular, we record that ash-poor normal explosions may be preceded and followed by the onset or the increase of the puffing activity, while ash-rich explosions are emergent, i.e., with inflation of the free surface followed directly by emission of increasingly large gas pockets. Overall, we see that all Strombolian activities form a continuum arising from a common mechanism and are modulated by the combination of two well-established controls: (1) the length of the bursting gas pocket with respect to the vent diameter and (2) the presence and thickness of a high-viscosity layer in the uppermost part of the volcanic conduit.

Journal of Geophysical Research B: Solid Earth

Temporal variation of tectonic tremor activity in southern Taiwan around the 2010 ML6.4 Jiashian Earthquake

Deep tectonic tremor, which is extremely sensitive to small stress variations, could be used to monitor fault-zone processes during large earthquake cycles and aseismic processes before large earthquakes. In this study, we develop an algorithm for the automatic detection and location of tectonic tremor beneath the southern Central Range of Taiwan and examine the spatio-temporal relationship between tremor and the 4 March 2010 ML6.4 Jiashian earthquake, located about 20 km from active tremor sources. We find that tremor in this region has a relatively short duration, short recurrence time, and no consistent correlation with surface GPS data. We find a short-term increase in the tremor rate 19 days before the Jiashian mainshock, and around the time when the tremor rate began to rise, one GPS station recorded a flip in its direction of motion. We hypothesize that tremor is driven by a slow-slip event that preceded the occurrence of the shallower Jiashian mainshock, even though the inferred slip is too small to be observed by all GPS stations. Our study shows that tectonic tremor may reflect stress variation during the pre-nucleation process of a nearby earthquake.

Journal of Geophysical Research B: Solid Earth

Effects of lava heating on volatile-rich slopes on Io

The upper crust of Io may be very rich in volatile sulfur and SO2. The surface is also highly volcanically active, and slopes may be warmed by radiant heat from the lava. This is particularly the case in paterae, which commonly host volcanic eruptions and long-lived lava lakes. Paterae slopes are highly variable, but some are greater than 70°. I model the heating of a volatile slope for two end-member cases: instantaneous emplacement of a large sheet flow, and persistent heating by a long-lived lava lake. In general, single flows can briefly raise sulfur to the melting temperature, or drive a modest amount of sublimation of SO2. Persistently lava-covered surfaces will drive much more significant geomorphic effects, with potentially significant sublimation and slope retreat. In addition to the direct effects, heating is likely to weaken slope materials and may trigger mass wasting. Thus, if the upper crust of Io is rich in these volatile species, future missions with high-resolution imaging are likely to observe actively retreating slopes around lava lakes and other locations of frequent eruptions.

Journal of Geophysical Research E: Planets

Geologic overview of the Mars Science Laboratory rover mission at the Kimberley, Gale crater, Mars

The Mars Science Laboratory (MSL) Curiosity rover completed a detailed investigation at the Kimberley waypoint within Gale crater from sols 571-634 using its full science instrument payload. From orbital images examined early in the Curiosity mission, the Kimberley region had been identified as a high-priority science target based on its clear stratigraphic relationships in a layered sedimentary sequence that had been exposed by differential erosion. Observations of the stratigraphic sequence at the Kimberley made by Curiosity are consistent with deposition in a prograding, fluvio-deltaic system during the late Noachian to early Hesperian, prior to the existence of most of Mt. Sharp. Geochemical and mineralogic analyses suggest that sediment deposition likely took place under cold conditions with relatively low water-to-rock ratios. Based on elevated K2O abundances throughout the Kimberley formation, an alkali feldspar protolith is likely one of several igneous sources from which the sediments were derived. After deposition, the rocks underwent multiple episodes of diagenetic alteration with different aqueous chemistries and redox conditions, as evidenced by the presence of Ca-sulfate veins, Mn-oxide fracture-fills, and erosion-resistant nodules. More recently, the Kimberley has been subject to significant aeolian abrasion and removal of sediments to create modern topography that slopes away from Mt. Sharp, a process that has continued to the present day.

Journal of Geophysical Research E: Planets

A probabilistic approach to remote compositional analysis of planetary surfaces

Reflected light from planetary surfaces provides information, including mineral/ice compositions and grain sizes, by study of albedo and absorption features as a function of wavelength. However, deconvolving the compositional signal in spectra is complicated by the nonuniqueness of the inverse problem. Trade-offs between mineral abundances and grain sizes in setting reflectance, instrument noise, and systematic errors in the forward model are potential sources of uncertainty, which are often unquantified. Here we adopt a Bayesian implementation of the Hapke model to determine sets of acceptable-fit mineral assemblages, as opposed to single best fit solutions. We quantify errors and uncertainties in mineral abundances and grain sizes that arise from instrument noise, compositional end members, optical constants, and systematic forward model errors for two suites of ternary mixtures (olivine-enstatite-anorthite and olivine-nontronite-basaltic glass) in a series of six experiments in the visible-shortwave infrared (VSWIR) wavelength range. We show that grain sizes are generally poorly constrained from VSWIR spectroscopy. Abundance and grain size trade-offs lead to typical abundance errors of ≤1 wt % (occasionally up to ~5 wt %), while ~3% noise in the data increases errors by up to ~2 wt %. Systematic errors further increase inaccuracies by a factor of 4. Finally, phases with low spectral contrast or inaccurate optical constants can further increase errors. Overall, typical errors in abundance are <10%, but sometimes significantly increase for specific mixtures, prone to abundance/grain-size trade-offs that lead to high unmixing uncertainties. These results highlight the need for probabilistic approaches to remote determination of planetary surface composition.

Journal of Geophysical Research E: Planets

Probability distributions of bed load particle velocities, accelerations, hop distances, and travel times informed by Jaynes's principle of maximum entropy

We describe the most likely forms of the probability distributions of bed load particle velocities, accelerations, hop distances, and travel times, in a manner that formally appeals to inferential statistics while honoring mechanical and kinematic constraints imposed by equilibrium transport conditions. The analysis is based on E. Jaynes's elaboration of the implications of the similarity between the Gibbs entropy in statistical mechanics and the Shannon entropy in information theory. By maximizing the information entropy of a distribution subject to known constraints on its moments, our choice of the form of the distribution is unbiased. The analysis suggests that particle velocities and travel times are exponentially distributed and that particle accelerations follow a Laplace distribution with zero mean. Particle hop distances, viewed alone, ought to be distributed exponentially. However, the covariance between hop distances and travel times precludes this result. Instead, the covariance structure suggests that hop distances follow a Weibull distribution. These distributions are consistent with high-resolution measurements obtained from high-speed imaging of bed load particle motions. The analysis brings us closer to choosing distributions based on our mechanical insight.

Journal of Geophysical Research F: Earth Surface

On extracting sediment transport information from measurements of luminescence in river sediment

Accurately quantifying sediment transport rates in rivers remains an important goal for geomorphologists, hydraulic engineers, and environmental scientists. However, current techniques for measuring long-time scale (102–106 years) transport rates are laborious, and formulae to predict transport are notoriously inaccurate. Here we attempt to estimate sediment transport rates by using luminescence, a property of common sedimentary minerals that is used by the geoscience community for geochronology. This method is advantageous because of the ease of measurement on ubiquitous quartz and feldspar sand. We develop a model from first principles by using conservation of energy and sediment mass to explain the downstream pattern of luminescence in river channel sediment. We show that the model can accurately reproduce the luminescence observed in previously published field measurements from two rivers with very different sediment transport styles. The model demonstrates that the downstream pattern of river sand luminescence should show exponential-like decay in the headwaters which asymptotes to a constant value with further downstream distance. The parameters from the model can then be used to estimate the time-averaged virtual velocity, characteristic transport lengthscale, storage time scale, and floodplain exchange rate of fine sand-sized sediment in a fluvial system. The sediment transport values predicted from the luminescence method show a broader range than those reported in the literature, but the results are nonetheless encouraging and suggest that luminescence demonstrates potential as a sediment transport indicator. However, caution is warranted when applying the model as the complex nature of sediment transport can sometimes invalidate underlying simplifications.

Journal of Geophysical Research F: Earth Surface

Constraining the relative importance of raindrop- and flow-driven sediment transport mechanisms in postwildfire environments and implications for recovery time scales

Mountain watersheds recently burned by wildfire often experience greater amounts of runoff and increased rates of sediment transport relative to similar unburned areas. Given the sedimentation and debris flow threats caused by increases in erosion, more work is needed to better understand the physical mechanisms responsible for the observed increase in sediment transport in burned environments and the time scale over which a heightened geomorphic response can be expected. In this study, we quantified the relative importance of different hillslope erosion mechanisms during two postwildfire rainstorms at a drainage basin in Southern California by combining terrestrial laser scanner-derived maps of topographic change, field measurements, and numerical modeling of overland flow and sediment transport. Numerous debris flows were initiated by runoff at our study area during a long-duration storm of relatively modest intensity. Despite the presence of a well-developed rill network, numerical model results suggest that the majority of eroded hillslope sediment during this long-duration rainstorm was transported by raindrop-induced sediment transport processes, highlighting the importance of raindrop-driven processes in supplying channels with potential debris flow material. We also used the numerical model to explore relationships between postwildfire storm characteristics, vegetation cover, soil infiltration capacity, and the total volume of eroded sediment from a synthetic hillslope for different end-member erosion regimes. This study adds to our understanding of sediment transport in steep, postwildfire landscapes and shows how data from field monitoring can be combined with numerical modeling of sediment transport to isolate the processes leading to increased erosion in burned areas.

California

Debris flow runup on vertical barriers and adverse slopes

Runup of debris flows against obstacles in their paths is a complex process that involves profound flow deceleration and redirection. We investigate the dynamics and predictability of runup by comparing results from large-scale laboratory experiments, four simple analytical models, and a depth-integrated numerical model (D-Claw). The experiments and numerical simulations reveal the important influence of unsteady, multidimensional flow on runup, and the analytical models highlight key aspects of the underlying physics. Runup against a vertical barrier normal to the flow path is dominated by rapid development of a shock, or jump in flow height, associated with abrupt deceleration of the flow front. By contrast, runup on sloping obstacles is initially dominated by a smooth flux of mass and momentum from the flow body to the flow front, which precedes shock development and commonly increases the runup height. D-Claw simulations that account for the emergence of shocks show that predicted runup heights vary systematically with the adverse slope angle and also with the Froude number and degree of liquefaction (or effective basal friction) of incoming flows. They additionally clarify the strengths and limitations of simplified analytical models. Numerical simulations based on a priori knowledge of the evolving dynamics of incoming flows yield quite accurate runup predictions. Less predictive accuracy is attained in ab initio simulations that compute runup based solely on knowledge of static debris properties in a distant debris flow source area. Nevertheless, the paucity of inputs required in ab initio simulations enhances their prospective value in runup forecasting.

Journal of Geophysical Research F: Earth Surface

Numerical simulation of large-scale bed load particle tracer advection-dispersion in rivers with free bars

Asymptotic characteristics of the transport of bed load tracer particles in rivers have been described by advection-dispersion equations. Here we perform numerical simulations designed to study the role of free bars, and more specifically single-row alternate bars, on streamwise tracer particle dispersion. In treating the conservation of tracer particle mass, we use two alternative formulations for the Exner equation of sediment mass conservation: the flux-based formulation, in which bed elevation varies with the divergence of the bed load transport rate, and the entrainment-based formulation, in which bed elevation changes with the net deposition rate. Under the condition of no net bed aggradation/degradation, a 1-D flux-based deterministic model that does not describe free bars yields no streamwise dispersion. The entrainment-based 1-D formulation, on the other hand, models stochasticity via the probability density function (PDF) of particle step length, and as a result does show tracer dispersion. When the formulation is generalized to 2-D to include free alternate bars, however, both models yield almost identical asymptotic advection-dispersion characteristics, in which streamwise dispersion is dominated by randomness inherent in free bar morphodynamics. This randomness can result in a heavy-tailed PDF of waiting time. In addition, migrating bars may constrain the travel distance through temporary burial, causing a thin-tailed PDF of travel distance. The superdiffusive character of streamwise particle dispersion predicted by the model is attributable to the interaction of these two effects.

Journal of Geophysical Research F: Earth Surface

Comment on “The reduction of friction in long-runout landslides as an emergent phenomenon” by Brandon C. Johnson et al.

Results from a highly idealized, 2-D computational model indicate that dynamic normal-stress rarefactions might cause friction reduction in long-runout landslides, but the physical relevance of the idealized dynamics has not been confirmed by experimental tests. More importantly, the model results provide no evidence that refutes alternative hypotheses about friction reduction mechanisms. One alternative hypothesis, which is strongly supported by field evidence, experimental data, and the predictions of a well-constrained computational model, involves development of high pore fluid pressures in deforming landslide material or overridden bed material. However, no scientific basis exists for concluding that a universal mechanism is responsible for friction reduction in all long-runout landslides.

Journal of Geophysical Research F: Earth Surface

A model integrating longshore and cross-shore processes for predicting long-term shoreline response to climate change

We present a shoreline change model for coastal hazard assessment and management planning. The model, CoSMoS-COAST (Coastal One-line Assimilated Simulation Tool), is a transect-based, one-line model that predicts short-term and long-term shoreline response to climate change in the 21 st century. The proposed model represents a novel, modular synthesis of process-based models of coastline evolution due to longshore and cross-shore transport by waves and sea-level rise. Additionally, the model uses an extended Kalman filter for data assimilation of historical shoreline positions to improve estimates of model parameters and thereby improve confidence in long-term predictions. We apply CoSMoS-COAST to simulate sandy shoreline evolution along 500 km of coastline in Southern California, which hosts complex mixtures of beach settings variably backed by dunes, bluffs, cliffs, estuaries, river mouths, and urban infrastructure, providing applicability of the model to virtually any coastal setting. Aided by data assimilation, the model is able to reproduce the observed signal of seasonal shoreline change for the hindcast period of 1995-2010, showing excellent agreement between modeled and observed beach states. The skill of the model during the hindcast period improves confidence in the model's predictive capability when applied to the forecast period (2010-2100) driven by GCM-projected wave and sea-level conditions. Predictions of shoreline change with limited human intervention indicate that 31% to 67% of Southern California beaches may become completely eroded by 2100 under sea-level rise scenarios of 0.93 to 2.0 m.

California

Terrestrial ecosystem model performance in simulating productivity and its vulnerability to climate change in the northern permafrost region

Realistic projection of future climate-carbon (C) cycle feedbacks requires better understanding and an improved representation of the C cycle in permafrost regions in the current generation of Earth system models. Here we evaluated 10 terrestrial ecosystem models for their estimates of net primary productivity (NPP) and responses to historical climate change in permafrost regions in the Northern Hemisphere. In comparison with the satellite estimate from the Moderate Resolution Imaging Spectroradiometer (MODIS; 246 ± 6 g C m −2 yr −1 ), most models produced higher NPP (309 ± 12 g C m −2 yr −1 ) over the permafrost region during 2000–2009. By comparing the simulated gross primary productivity (GPP) with a flux tower-based database, we found that although mean GPP among the models was only overestimated by 10% over 1982–2009, there was a twofold discrepancy among models (380 to 800 g C m −2 yr −1 ), which mainly resulted from differences in simulated maximum monthly GPP (GPP max ). Most models overestimated C use efficiency (CUE) as compared to observations at both regional and site levels. Further analysis shows that model variability of GPP and CUE are nonlinearly correlated to variability in specific leaf area and the maximum rate of carboxylation by the enzyme Rubisco at 25°C ( V c max_25 ), respectively. The models also varied in their sensitivities of NPP, GPP, and CUE to historical changes in climate and atmospheric CO 2 concentration. These results indicate that model predictive ability of the C cycle in permafrost regions can be improved by better representation of the processes controlling CUE and GPP max as well as their sensitivity to climate change.

Journal of Geophysical Research: Biogeosciences

DOM composition and transformation in boreal forest soils: The effects of temperature and organic-horizon decomposition state

The boreal region stores large amounts of organic carbon (C) in organic-soil horizons, which are vulnerable to destabilization via warming and disturbance. Decomposition of soil organic matter (SOM) contributes to the production and turnover of dissolved organic matter (DOM). While temperature is a primary control on rates of SOM and DOM cycling, little is known about temperature effects on DOM composition in soil leachate. Here we conducted a 30 day incubation to examine the effects of temperature (20 versus 5°C) and SOM decomposition state (moss versus fibric versus amorphous horizons) on DOM composition in organic soils of interior Alaska. We characterized DOM using bulk dissolved organic C (DOC) concentration, chemical fractionation, optical properties, and ultrahigh-resolution mass spectrometry. We observed an increase in DOC concentration and DOM aromaticity in the 20°C treatment compared to the 5°C treatment. Leachate from fibric horizons had higher DOC concentration than shallow moss or deep amorphous horizons. We also observed chemical shifts in DOM leachate over time, including increases in hydrophobic organic acids, polyphenols, and condensed aromatics and decreases in low-molecular weight hydrophilic compounds and aliphatics. We compared ultrahigh-resolution mass spectrometry and optical data and observed strong correlations between polyphenols, condensed aromatics, SUVA 254 , and humic-like fluorescence intensities. These findings suggest that biolabile DOM was preferentially mineralized, and the magnitude of this transformation was determined by kinetics (i.e., temperature) and substrate quality (i.e., soil horizon). With future warming, our findings indicate that organic soils may release higher concentrations of aromatic DOM to aquatic ecosystems.

Journal of Geophysical Research: Biogeosciences

Nitrate reduction mechanisms and rates in an unconfined eogenetic karst aquifer in two sites with different redox potential

This study integrates push-pull tracer tests (PPTT) with microbial characterization of extracted water via quantitative polymerase chain reaction (qPCR) and reverse transcriptase qPCR (RT-qPCR) of selected functional N transformation genes to quantify nitrate reduction mechanisms and rates in sites with different redox potential in a karst aquifer. PPTT treatments with nitrate (AN) and nitrate-fumarate (ANC) were executed in two wells representing anoxic and oxic geochemical end-members. Oxic aquifer zero-order nitrate loss rates (mmol L −1 h −1 ) were similar for AN and ANC treatment, ranging from 0.03 ± 0.01 to 0.05 ± 0.01. Anoxic aquifer zero-order nitrate loss rates ranged from 0.03 ± 0.02 (AN) to 0.13 ± 0.02 (ANC). Microbial characterization indicates mechanisms influencing these rates were dissimilatory nitrate reduction to ammonium (DNRA) at the anoxic site with AN treatment, assimilatory reduction of nitrate to ammonium (ANRA) with ANC treatment in the water column at both sites, and additional documented nitrate reduction that occurred in unsampled biofilms. With carbon treatment, total numbers of microbes (16S rRNA genes) significantly increased (fourteenfold to thirtyfold), supporting stimulated growth with resulting ANRA. Decreased DNRA gene concentrations ( nrfA DNA) and increased DNRA activity ratio (nrfA-cDNA/DNA) supported the assertion that DNRA occurred in the anoxic zone with AN and ANC treatment. Furthermore, decreased DNRA gene copy numbers at the anoxic site with ANC treatment suggests that DNRA microbes in the anoxic site are chemolithoautotrophic. Increased RT-qPCR denitrification gene expression ( nirK and nirS ) was not observed in water samples, supporting that any observed NO 3 -N loss due to denitrification may be occurring in unsampled microbial biofilms.

Florida