USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysics Research”

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 1,117 records · Page 62Linked to original sources

Partly cloudy with a chance of lava flows: Forecasting volcanic eruptions in the 21st century

A primary goal of volcanology is forecasting hazardous eruptive activity. Despite much progress over the last century, however, volcanoes still erupt with no detected precursors, lives and livelihoods are lost to eruptive activity, and forecasting the onsets of eruptions remains fraught with uncertainty. Long‐term forecasts are generally derived from the geological and historical records, from which recurrence intervals and styles of activity can be inferred, while shorter‐term forecasts are derived from patterns in monitoring data. Information from geology and monitoring data can be evaluated and combined using statistical analysis, expert elicitation, and conceptual and or mathematical models. Integrative frameworks, such as event trees, combine this diversity of information to produce probabilistic forecasts that can inform the style and scale of the societal response to a potential future eruption. Several developments show promise to revolutionize the utility and accuracy of these forecasts. These include growth in the quantity and quality of multidisciplinary monitoring data, coupled with increases in computing power; machine learning algorithms, which will allow far better utilization of this growing volume of data; and new physiochemical volcano models and data assimilation algorithms, which take advantage of a wide range of monitoring data and realistic physics to better predict the evolution of a given physical state. Although eruption forecasts may never be as generally reliable as weather forecasts, and great caution must be exercised when attempting to predict highly complex volcanic behavior, these and other innovations—particularly when combined in integrative, fully probabilistic forecasting frameworks—should help volcanologists to better issue warnings of volcanic activity on societally relevant time frames.

Hawaii↗

Depth determination of the 2010 El Mayor‐Cucapah earthquake sequence (M ≥ 4.0)

The 2010 M W 7.2 El Mayor‐Cucapah earthquake ruptured a zone of ~120 km in length in northern Baja California. The geographic distribution of this earthquake sequence was well constrained by waveform relocation. The depth distribution, however, was poorly determined as it is near the edge of, or outside, the Southern California Seismic Network. Here we use two complementary methods to constrain the focal depths of moderate‐sized events ( M ≥ 4.0) in this sequence. We first determine the absolute earthquake depth by modeling the regional depth phases at high frequencies (~1 Hz). We mainly focus on Pn and its depth phases pPn and sPn , which arrive early at regional distance and are less contaminated by crustal multiples. To facilitate depth phase identification and to improve signal‐to‐noise ratio, we take advantage of the dense Southern California Seismic Network and use array analysis to align and stack Pn waveforms. For events without clear depth phases, we further determine their relative depths with respect to those with known depths using differential travel times of the Pn , direct P , and direct S phases recorded for event pairs. Focal depths of 93 out of 122 M ≥ 4.0 events are tightly constrained with absolute uncertainty of about 1 km. Aftershocks are clustered in the depth range of 3–10 km, suggesting a relatively shallow seismogenic zone, consistent with high surface heat flow in this region. Most aftershocks are located outside or near the lower terminus of coseismic high‐slip patches of the main shock, which may be governed by residual strains, local stress concentration, or postseismic slip.

Baja California↗

Assessing seasonal changes in microgravity at Yellowstone caldera

Microgravity time series at active volcanoes can provide an indication of mass change related to subsurface magmatic processes, but uncertainty is often introduced by hydrologic variations and other noise sources that cannot easily be isolated. We empirically assessed seasonality and noise by conducting four surveys over the course of May-October 2017 at Yellowstone caldera, Wyoming. Yellowstone experiences frequent changes in the rates and styles of seismicity and surface deformation, but the mechanisms of these changes are poorly understood because the characteristics of the driving fluids are not clear. Past gravity data from the caldera have yielded ambiguous results, possibly due to hydrologic noise. Given the strong visually observable changes in surface water and snow conditions over the course of our surveys, we expected to see significant variations in gravity. The net change in gravity, however, was less than 20 µGal at most sites, and there was no strong correlation with river and lake levels or snow conditions. Seasonal changes in gravity are therefore small compared to those that would be expected from magmatic activity, although they may be on the same order as those associated with Yellowstone’s hydrothermal system. We did find that noise levels in gravity data were highly dependent on site characteristics, with bedrock sites away from trees yielding the lowest levels of noise, and thin concrete pads in forested areas the highest. These results can be used to plan future surveys at Yellowstone and to reinterpret past data, and they provide guidance in terms of best practices for repeat gravity work on volcanoes worldwide.

Wyoming↗

Surface imaging functions for elastic reverse time migration

Reverse time migration is often used to interpret acoustic or three‐component seismic recordings by creating an image of subsurface seismic reflectors. Here I describe elastic reverse time migration imaging functions that are cast as waveform misfit sensitivity kernels of contrasts in material parameters across hypothetical seismic discontinuities, that is, specular reflectors. The proposed “surface” imaging functions are theoretically applicable to either reflected or converted waves in order to estimate the location and reflectivity of these discontinuities. The surface imaging functions, as well as volumetric sensitivity kernels that target point diffractors, are tested on sets of synthetic surface array recordings that sample the 3‐D seismic wavefield on simple 2‐D structures generated using a 2.5‐D spectral element method. These tests illustrate that in contrast with the volumetric sensitivity kernels, the reflectivity is generally dominated by a high‐amplitude peak that coincides with input locations of discontinuities. Passive recordings of microseismicity, shot gathers, or a combination thereof can be potentially interpreted with the new surface imaging functions to yield useful reflectivity images.

Journal of Geophysical Research Solid Earth↗

Peak ground displacement saturates exactly when expected: Implications for earthquake early warning

The scaling of rupture properties with magnitude is of critical importance to earthquake early warning (EEW) systems that rely on source characterization using limited snapshots of waveform data. ShakeAlert, a prototype EEW system that is being developed for the western United States, provides real-time estimates of earthquake magnitude based on P-wave peak ground displacements measured at stations triggered by the event. The algorithms used in ShakeAlert assume that the displacement measurements at each station are statistically independent and that there exists a linear and time-independent relation between log peak ground displacement and earthquake magnitude. Here we challenge this basic assumption using a comprehensive database of more than 130,000 vertical component waveforms from M4.5-M9 earthquakes occurring near Japan from 1997 through 2017 and recorded by the K-NET and KiK-net strong-motion networks. By analyzing the time-evolution of P-wave peak ground displacements for these earthquakes, we show that there is a break, or saturation, in the magnitude-displacement scaling that depends on the length of the measurement time window. We demonstrate that the magnitude at which this saturation occurs is well-explained by a simple and non-deterministic model of earthquake rupture growth. We then use the predictions of this saturation model to develop a Bayesian framework for estimating posterior uncertainties in real-time magnitude estimates which incorporates the expected time-dependence of the peak displacement measurements.

Journal of Geophysical Research B: Solid Earth↗

Vertical coseismic offsets from differential high-resolution stereogrammetric DSMs: The 2013 Baluchistan, Pakistan earthquake

The recent proliferation of high-resolution (< 3-m spatial resolution) digital topography datasets opens a spectrum of geodetic applications in differential topography, including the quantification of coseismic vertical displacement fields. Most investigations of coseismic vertical displacements to date rely, in part, on pre- or post-event lidar surveys that are intractable or non-existent in many locales. Stereogrammetric digital surface models (DSMs) derived from high-resolution satellite optical imagery provide a new avenue for the retrieval of spatially-dense vertical coseismic displacements on a global scale. In this study, we generated 2-m resolution pre- and post-seismic DSMs from satellite optical imagery spanning the 2013 Mw7.7 Baluchistan strike-slip earthquake that occurred on the Hoshab fault in southern Pakistan. We applied the Iterative Closest Point algorithm to the DSMs to quantify the coseismic vertical displacement field at a spatial resolution of 10-30 m and to generate 3D coseismic strain tensors. We found that across-fault vertical offsets alternated between uplift and subsidence and varied between ~1-3 m in a non-systematic manner along the Hoshab fault. We show that the pre-existing topography and near-fault geomorphology are variably consistent and inconsistent with the displacement kinematics of the 2013 earthquake, and we argue that these relationships highlight varied slip sense history along the Hoshab fault. Notably, topography along the southern extents of the Hoshab fault requires different surface displacement kinematics than occurred in the 2013 earthquake, suggesting that the Hoshab fault accommodates varying senses of slip (bimodal slip) through time.

Baluchistan Province↗

Sea level rise in the Samoan Islands escalated by viscoelastic relaxation after the 2009 Samoa‐Tonga earthquake

The Samoan islands are an archipelago hosting a quarter million people mostly residing in three major islands, Savai'i and Upolu (Samoa), and Tutuila (American Samoa). The islands have experienced sea level rise by 2–3 mm/year during the last half century. The rate, however, has dramatically increased following the M w 8.1 Samoa‐Tonga earthquake doublet (megathrust + normal faulting) in September 2009. Since the earthquake, we found large‐scale gravity increase (0.5 μ Gal/year) around the islands and ongoing subsidence (8–16 mm/year) of the islands from our analysis of Gravity Recovery And Climate Experiment gravity and GPS displacement data. The postseismic horizontal displacement is faster in Samoa, while the postseismic subsidence rate is considerably larger in American Samoa. The analysis of local tide gauge records and satellite altimeter data also identified that the relative sea level rise becomes faster by 7–9 mm/year in American Samoa than Samoa. A simple viscoelastic model with a Maxwell viscosity of 2–3 × 10 18 Pa s for the asthenosphere explained postseismic deformation at nearby GPS sites as well as Gravity Recovery And Climate Experiment gravity change. It is found that the constructive interference of viscoelastic relaxation from both megathrust and normal faulting has intensified the postseismic subsidence at American Samoa, causing ~5 times faster sea level rise than the global average. Our model indicates that this trend is likely to continue for decades and result in sea level rise of 30–40 cm, which is independent of and in addition to anticipated climate‐related sea level rise. It will worsen coastal flooding on the islands leading to regular nuisance flooding.

Journal of Geophysical Research↗

Groundwater inflow toward a preheated volcanic conduit: Application to the 2018 eruption at Kīlauea Volcano, Hawai’i

The many successes in volcano forecasting over the past several decades owe mainly to pattern recognition, both in monitoring data and the geologic record. During the early stages of the 2018 Kīlauea eruption, the conceptual model of Stearns (1925), based on the explosive 1924 Kīlauea eruption, was highly influential. This model postulates that explosions are triggered by liquid-water inflow into a recently vacated magma conduit. Modern quantitative modeling approaches, supplemented by hydrogeologic data unavailable in 1925, yield a more nuanced view. Results demonstrate that liquid-water inflow would likely be delayed by months to years, owing to the inability of liquid water to transit a zone of very hot rock surrounding the conduit. The exercise demonstrates the use of physically based modeling to supplement traditional volcano-forecasting approaches during an ongoing event.

Hawaii↗

Delayed dynamic triggering of disposal-induced earthquakes observed by a dense array in Northern Oklahoma

Recent increases in earthquake occurrence rates in Oklahoma have been linked to the injection of large volumes of saltwater, a byproduct of oil and gas extraction. Here we present a detailed study of remote earthquake triggering in an area of active injection‐induced seismicity in northern Oklahoma using data from the LArge‐n Seismic Survey in Oklahoma (LASSO) temporary array and nearby permanent broadband seismic stations. We estimate changes in earthquake rates and calculate the Coulomb failure stress changes on potential receiver faults due to passing teleseismic surface waves. A statistically significant increase in seismicity is observed ∼8 hr after the 16 April 2016 M w 7.8 Ecuador earthquake. The Coulomb stress changes associated with the Ecuador earthquake are on the order of ∼1 kPa. Physical mechanisms consistent with the observed dynamic stress threshold include failure driven by activation of aseismic slip or hydrological response of the fault system.

Oklahoma↗

Offshore landslide hazard curves from mapped landslide size distributions

We present a method to calculate landslide hazard curves along offshore margins based on size distributions of submarine landslides. The method analyzes ten different continental margins, that were mapped by high-resolution multibeam sonar with landslide scar areas measured by a consistent GIS procedure. Statistical tests of several different probability distribution models indicate that the lognormal model is most appropriate for these siliciclastic environments, consistent with an earlier study of the U.S. Atlantic margin [Chaytor et al., 2009]. Parameter estimation is performed using the maximum likelihood technique and confidence intervals are determined using likelihood profiles. Pairwise comparison of size distributions for the ten margins indicates that the U.S. Atlantic and Queen Charlotte margins are different than most other margins. These margins represent end members, with the U.S. Atlantic margin having the highest mean scar area and the Queen Charlotte margin, the lowest. We demonstrate that empirical, offshore landslide hazard curves can be developed from the landslide size distributions, if the duration of mapped landslide activity is known. This study indicates that the shape parameter of the size distribution is similar among all ten margins and thus the shape of the hazard curves is also similar. Significant differences in hazard curves among the margins are therefore related to differences in mean sizes and, potentially, differences in the duration of landslide activity.

Journal of Geophysical Research B: Solid Earth↗

Locating shallow seismic sources with waves scattered by surface topography: Validation of the method at the Nevada Test Site

Accurate and robust source location is fundamental in seismology. Previously, we proposed a new full waveform location method using waves scattered by topography near the source, and we demonstrated its potential in obtaining accurate source location using synthetic data. In the work presented here, we validate this method with observed seismic data from the 1993 Non-Proliferation Experiment, a chemical explosion with a precisely known location in a region with moderate topography at the Nevada Test Site. We select the first arrivals (Pn/Pg) and their immediate codas to determine the source location and compare our solution with the known location. We use a collocated-grid finite-difference method to calculate the strain Green’s tensor in a grid-search volume containing the source, and obtain three-component synthetic waveforms at 12 broadband seismic stations at regional distances using source-receiver reciprocity. We assess the solution using a least-squares misfit between the observed and synthetic waveforms. When scattered coda waves are incorporated, the best solution is within a few hundreds of meters of the exact source location, and the estimated uncertainty of the solution is reduced compared to the waveform based solution using only the P waves. The solution is robust to the choice of the frequency content and to the addition of random velocity heterogeneity. We conclude that the full waveform source location method is effective for real seismic data, though more validation tests like this one are needed to further understand its efficacy for source-station geometry, roughness of topography, signal quality and other factors.

Nevada↗

Hydrodynamics and sediment mobility processes over a degraded senile coral reef

Coral reefs can influence hydrodynamics and morphodynamics by dissipating and refracting incident wave energy, modifying circulation patterns, and altering sediment transport pathways. In this study, the sediment and hydrodynamic response of a senile (dead) barrier reef (Crocker Reef, located in the upper portion of the Florida Reef Tract) to storms and quiescent conditions was evaluated using field observations and the Coupled Ocean-Atmosphere-Wave-Sediment Transport model. Waves, circulation, and resultant sediment mobility were modeled across different reef zones. Sediment mobility during quiescent periods and the passage of far-field storms are driven by nonbreaking waves and, to a lesser degree, regional circulation. Spatial variability in these processes produces the present-day distribution of sediment grain size at Crocker Reef, wherein finer-grain material along a shallow central ridge is frequently mobilized (43% to 62% of the time), winnowed away, and deposited along the lower-energy flanks and in the fore reef where sand mobility occurs less frequently (32% to 43% and 1% to 22% of the time, respectively). Analysis of wave conditions for the period of 2006–2014 supports that wave heights rarely exceed the threshold for breaking (0.1% and 0.3% at the reef crest and at the reef flat, respectively), predominantly during the passage of tropical storms. There is a shift to a wave-breaking regime during near-field storms, creating the potential for mobilization of larger material and enhanced reef degradation. Sediment mobility can be enhanced due to wave skewness or the generation of free infragravity waves during periods of depth-induced wave breaking.

Florida↗

Mechanisms of wave‐driven water level variability on reef‐fringed coastlines

Wave‐driven water level variability (and runup at the shoreline) is a significant cause of coastal flooding induced by storms. Wave runup is challenging to predict, particularly along tropical coral reef‐fringed coastlines due to the steep bathymetric profiles and large bottom roughness generated by reef organisms, which can violate assumptions in conventional models applied to open sandy coastlines. To investigate the mechanisms of wave‐driven water level variability on a reef‐fringed coastline, we performed a set of laboratory flume experiments on an along‐shore uniform bathymetric profile with and without bottom roughness. Wave setup and waves at frequencies lower than the incident sea‐swell forcing (infragravity waves) were found to be the dominant components of runup. These infragravity waves were positively correlated with offshore wave groups, signifying they were generated in the surf zone by the oscillation of the breakpoint. On the reef flat and at the shoreline, the low‐frequency waves formed a standing wave pattern with energy concentrated at the natural frequencies of the reef flat, indicating resonant amplification. Roughness elements used in the flume to mimic large reef bottom roughness reduced low frequency motions on the reef flat and reduced wave run up by 30% on average, compared to the runs over a smooth bed. These results provide insight into sea‐swell and infragravity wave transformation and wave setup dynamics on steep‐sloped coastlines, and the effect that future losses of reef bottom roughness may have on coastal flooding along reef‐fringed coasts.

Journal of Geophysical Research C: Oceans↗

Hydrodynamics of a tidally‐forced coral reef atoll

The hydrodynamics of a tidally forced semi‐enclosed coral reef atoll (North Scott) at the edge of the continental shelf of northwestern Australia were investigated by combining field observations and numerical modeling. The observations revealed that the spring tidal range outside the atoll reaches 4 m, and as the water level drops below mean sea level, the reef rim surrounding the shallow (~10–15 m) lagoon becomes exposed. During this time, the lagoon can only exchange with the open ocean through two narrow channels, resulting in highly asymmetric water levels and velocities that were most pronounced during spring tide. On average, the ebb tide duration was ~2 hr longer than the flood, with rapid flood velocities in the channel reaching 2 m/s. We applied an unstructured grid model Delft3D‐Flexible Mesh to simulate the atoll hydrodynamics and were able to replicate the asymmetric water levels and complex velocities in the lagoon. The results revealed that at higher tidal stages, a dominant momentum balance exists between the pressure gradient (established by the propagation of the tide on the shelf) and the local flow acceleration of water throughout the interior of the atoll. At lower tidal stages, which coincided with a reversal of the offshore tidal pressure gradient, the lagoon became isolated from offshore dynamics and all momentum terms were negligible. This resulted in a tidally averaged residual westward flow within the lagoon that drove an asymmetric flushing pattern within the atoll, which we propose would be a common flushing mechanism within other tide‐dominated atolls worldwide.

Journal of Geophysical Research C: Oceans↗

Field observations of alongshore runup variability under dissipative conditions in presence of a shoreline sandwave

Video measurements of runup were collected at low tide along several profiles covering an alongshore distance of 500 m. The morphology displayed a complex shape with a shoreline sandwave in the lower beach face of about 250 m long mirrored in the inner sandbar. Wave conditions were stationary and moderate (offshore height of 2 m and peak period of nearly 13 s) but yet dissipative. Runup energy was dominated by infragravity frequencies. Alongshore variations in runup (by a factor up to 3) observed both in the incident and infragravity bands were much higher than reported previously (e.g., Guedes et al., 2012, https://doi.org/10.1016/j.csr.2012.08.022 ; Ruggiero et al., 2004, https://doi.org/10.1029/2003JC002160 ) while the alongshore variations in other environmental parameters (e.g., foreshore beach slope) appear to be much lower. Our data suggest that the beach morphology in the inner surf zone plays a crucial role by inducing rapid and significant modification in the incident wave pattern and the alongshore coherence length scales were consistent with the typical alongshore length scale of the morphology.

Journal of Geophysical Research C: Oceans↗

Spatial and temporal variability of pCO2, carbon fluxes and saturation state on the West Florida Shelf

The West Florida Shelf (WFS) is a source of uncertainty for the Gulf of Mexico carbon budget. Data from the synthesis of approximately 135,000 pCO2 values from over 96 cruises from the WFS show that the shelf waters fluctuate between being a weak source to a weak sink of carbon with the atmosphere. Overall, the shelf acts as a weak source of CO2 at 0.32 ± 1.5 mol m-2 yr-1. Subregions, however, reveal slightly different trends, where surface waters associated with 40 m – 200 m isobaths in the northern and southern WFS are generally weak sinks all year, except for summer when they act as sources of CO2. Conversely, nearshore waters (< 40 m) are a source of CO2 are a source all year round, particularly the southern shallow waters. The pCO2 of seawater has been increasing at a rate of approximately 5.26 µatm yr-1 as compared to atmospheric pCO2 which has increased at a rate of about 1.7 µatm yr-1 from 1996 to 2016. The pCO2 and CO2 flux on the shelf from 1996 - 2016 have increased about 49 µatm, and 1.08 mol m-2, respectively. The WFS is emitting 9.23 Tg C yr-1, with the southern nearshore region emitting the most at 9.01 Tg C yr-1 and the northern region acting as a sink of -1.96 Tg C yr-1.

West Florida shelf↗

Submarine permafrost map in the arctic modelled using 1D transient heat flux (SuPerMAP)

Offshore permafrost plays a role in the global climate system, but observations of permafrost thickness, state, and composition are limited to specific regions. The current global permafrost map shows potential offshore permafrost distribution based on bathymetry and global sea level rise. As a first‐order estimate, we employ a heat transfer model to calculate the subsurface temperature field. Our model uses dynamic upper boundary conditions that synthesize Earth System Model air temperature, ice mass distribution and thickness, and global sea level reconstruction and applies globally distributed geothermal heat flux as a lower boundary condition. Sea level reconstruction accounts for differences between marine and terrestrial sedimentation history. Sediment composition and pore water salinity are integrated in the model. Model runs for 450 ka for cross‐shelf transects were used to initialize the model for circumarctic modeling for the past 50 ka. Preindustrial submarine permafrost (i.e., cryotic sediment), modeled at 12.5‐km spatial resolution, lies beneath almost 2.5 ×10 6 km 2 of the Arctic shelf. Our simple modeling approach results in estimates of distribution of cryotic sediment that are similar to the current global map and recent seismically delineated permafrost distributions for the Beaufort and Kara seas, suggesting that sea level is a first‐order determinant for submarine permafrost distribution. Ice content and sediment thermal conductivity are also important for determining rates of permafrost thickness change. The model provides a consistent circumarctic approach to map submarine permafrost and to estimate the dynamics of permafrost in the past.

Journal of Geophysical Research C: Oceans↗

Seasonal, spring-neap, and tidal variation in cohesive sediment transport parameters in estuarine shallows

Numerical models for predicting sediment concentrations and transport rely on parameters such as settling velocity and bed erodibility that describe sediment characteristics, yet these parameters are rarely probed directly. We investigated temporal and spatial variation in sediment parameters in the shallows of San Pablo Bay, CA. Flow, turbulence, and suspended sediment data were measured at sites located at 1 and 2 m below mean lower low water (MLLW) from November 2013 through April 2015, supplemented by monthlong periods in 2011, 2012, and 2016. Maximum current velocities were 0.40–0.47 m s -1 at these depths; the strongest currents decreased to 0.27–0.34 m s -1 during neap periods. Winters 2013–2014 and 2014–2015 experienced strong drought conditions, limiting the potential for seasonal impact on sediment conditions during this experiment. Despite this, the more storm‐influenced site showed clear changes during the winter: the roughness parameter decreased from 10 −4 to 10 −5 m, from hydrodynamically rough to smooth conditions, and bed erodibility increased by an order of magnitude. Median settling velocity was 2.05·10 −4 m s -1 ; it varied twofold within a tidal cycle, decreasing as current velocity grew during flood and ebb. This tidal control on floc size affected settling velocity on the spring‐neap timescale, possibly driving a spring‐neap oscillation in erodibility. Our findings highlight variation in sediment dynamics that is commonly ignored in numerical models and the need for field observations to ground truth ongoing modeling efforts.

California↗