USGS ScienceSearch

SEARCH · USGS Science

Results for “Journal of Geophysical Research - Atmospheres”

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

Integrated geophysical characteristics of the 2015 Illapel, Chile, earthquake

On 16 September 2015, a Mw 8.3 earthquake ruptured the subduction zone offshore of Illapel, Chile, generating an aftershock sequence with 14 Mw 6.0–7.0 events. A double source W phase moment tensor inversion consists of a Mw 7.2 subevent and the main Mw 8.2 phase. We determine two slip models for the mainshock, one using teleseismic broadband waveforms and the other using static GPS and InSAR surface displacements, which indicate high slip north of the epicenter and west-northwest of the epicenter near the oceanic trench. These models and slip distributions published in other studies suggest spatial slip uncertainties of ~25 km and have peak slip values that vary by a factor of 2. We relocate aftershock hypocenters using a Bayesian multiple-event relocation algorithm, revealing a cluster of aftershocks under the Chilean coast associated with deep (20–45 km depth) mainshock slip. Less vigorous aftershock activity also occurred near the trench and along strike of the main aftershock region. Most aftershocks are thrust-faulting events, except for normal-faulting events near the trench. Coulomb failure stress change amplitudes and signs are uncertain for aftershocks collocated with deeper mainshock slip; other aftershocks are more clearly associated with loading from the mainshock. These observations reveal a frictionally heterogeneous interface that ruptured in patches at seismogenic depths (associated with many aftershocks) and with homogeneous slip (and few aftershocks) up to the trench. This event likely triggered seismicity separate from the main slip region, including along-strike events on the megathrust and intraplate extensional events.

Illapel

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

Sediment gravity flows triggered by remotely generated earthquake waves

Recent great earthquakes and tsunamis around the world have heightened awareness of the inevitability of similar events occurring within the Cascadia Subduction Zone of the Pacific Northwest. We analyzed seafloor temperature, pressure, and seismic signals, and video stills of sediment-enveloped instruments recorded during the 2011–2015 Cascadia Initiative experiment, and seafloor morphology. Our results led us to suggest that thick accretionary prism sediments amplified and extended seismic wave durations from the 11 April 2012 M w 8.6 Indian Ocean earthquake, located more than 13,500 km away. These waves triggered a sequence of small slope failures on the Cascadia margin that led to sediment gravity flows culminating in turbidity currents. Previous studies have related the triggering of sediment-laden gravity flows and turbidite deposition to local earthquakes, but this is the first study in which the originating seismic event is extremely distant (> 10,000 km). The possibility of remotely triggered slope failures that generate sediment-laden gravity flows should be considered in inferences of recurrence intervals of past great Cascadia earthquakes from turbidite sequences. Future similar studies may provide new understanding of submarine slope failures and turbidity currents and the hazards they pose to seafloor infrastructure and tsunami generation in regions both with and without local earthquakes.

Cascadia subduction zone

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

Behavior of a wave-driven buoyant surface jet on a coral reef

A wave-driven surface buoyant jet exiting a coral reef was studied in order to quantify the amount of water re-entrained over the reef crest. Both moored observations and Lagrangian drifters were used to study the fate of the buoyant jet. To investigate in detail the effects of buoyancy and along-shore flow variations, we developed an idealized numerical model of the system. Consistent with previous work, the ratio of along-shore velocity to jet-velocity and the jet internal Froude number were found to be important determinants of the fate of the jet. In the absence of buoyancy, the entrainment of fluid at the reef crest, creates a significant amount of retention, keeping 60% of water in the reef system. However, when the jet is lighter than the ambient ocean-water, the net effect of buoyancy is to enhance the separation of the jet from shore, leading to a greater export of reef water. Matching observations, our modeling predicts that buoyancy limits retention to 30% of the jet flow for conditions existing on the Moorea reef. Overall, the combination of observations and modeling we present here shows that reef-ocean temperature gradients can play an important role in reef-ocean exchanges.

Journal of Geophysical Research C: Oceans

Sediment transport in the presence of large reef bottom roughness

The presence of large bottom roughness, such as that formed by benthic organisms on coral reef flats, has important implications for the size, concentration, and transport of suspended sediment in coastal environments. A 3 week field study was conducted in approximately 1.5 m water depth on the reef flat at Ningaloo Reef, Western Australia, to quantify the cross-reef hydrodynamics and suspended sediment dynamics over the large bottom roughness (∼20–40 cm) at the site. A logarithmic mean current profile consistently developed above the height of the roughness; however, the flow was substantially reduced below the height of the roughness (canopy region). Shear velocities inferred from the logarithmic profile and Reynolds stresses measured at the top of the roughness, which are traditionally used in predictive sediment transport formulations, were similar but much larger than that required to suspend the relatively coarse sediment present at the bed. Importantly, these stresses did not represent the stresses imparted on the sediment measured in suspension and are therefore not relevant to the description of suspended sediment transport in systems with large bottom roughness. Estimates of the bed shear stresses that accounted for the reduced near-bed flow in the presence of large roughness vastly improved the relationship between the predicted and observed grain sizes that were in suspension. Thus, the impact of roughness, not only on the overlying flow but also on bed stresses, must be accounted for to accurately estimate suspended sediment transport in regions with large bottom roughness, a common feature of many shallow coastal ecosystems.

Journal of Geophysical Research C: Oceans

Rip currents and alongshore flows in single channels dredged in the surf zone

To investigate the dynamics of flows near nonuniform bathymetry, single channels (on average 30 m wide and 1.5 m deep) were dredged across the surf zone at five different times, and the subsequent evolution of currents and morphology was observed for a range of wave and tidal conditions. In addition, circulation was simulated with the numerical modeling system COAWST, initialized with the observed incident waves and channel bathymetry, and with an extended set of wave conditions and channel geometries. The simulated flows are consistent with alongshore flows and rip-current circulation patterns observed in the surf zone. Near the offshore-directed flows that develop in the channel, the dominant terms in modeled momentum balances are wave-breaking accelerations, pressure gradients, advection, and the vortex force. The balances vary spatially, and are sensitive to wave conditions and the channel geometry. The observed and modeled maximum offshore-directed flow speeds are correlated with a parameter based on the alongshore gradient in breaking-wave-driven-setup across the nonuniform bathymetry (a function of wave height and angle, water depths in the channel and on the sandbar, and a breaking threshold) and the breaking-wave-driven alongshore flow speed. The offshore-directed flow speed increases with dissipation on the bar and reaches a maximum (when the surf zone is saturated) set by the vertical scale of the bathymetric variability.

Journal of Geophysical Research: Oceans

Observations and a linear model of water level in an interconnected inlet-bay system

A system of barrier islands and back-barrier bays occurs along southern Long Island, New York, and in many coastal areas worldwide. Characterizing the bay physical response to water level fluctuations is needed to understand flooding during extreme events and evaluate their relation to geomorphological changes. Offshore sea level is one of the main drivers of water level fluctuations in semienclosed back-barrier bays. We analyzed observed water levels (October 2007 to November 2015) and developed analytical models to better understand bay water level along southern Long Island. An increase (∼0.02 m change in 0.17 m amplitude) in the dominant M 2 tidal amplitude (containing the largest fraction of the variability) was observed in Great South Bay during mid-2014. The observed changes in both tidal amplitude and bay water level transfer from offshore were related to the dredging of nearby inlets and possibly the changing size of a breach across Fire Island caused by Hurricane Sandy (after December 2012). The bay response was independent of the magnitude of the fluctuations (e.g., storms) at a specific frequency. An analytical model that incorporates bay and inlet dimensions reproduced the observed transfer function in Great South Bay and surrounding areas. The model predicts the transfer function in Moriches and Shinnecock bays where long-term observations were not available. The model is a simplified tool to investigate changes in bay water level and enables the evaluation of future conditions and alternative geomorphological settings.

Journal of Geophysical Research C: Oceans

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