USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research - Planets”

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 919 records · Page 51Linked to original sources
USGS2019

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↗
USGS2019

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↗
USGS2019

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↗
USGS2019

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↗
USGS2019

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↗
USGS2019

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↗
USGS2019

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↗
USGS2018

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↗
USGS2018

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↗
USGS2018

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↗
USGS2018

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↗
USGS2018

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↗
USGS2019

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↗
USGS2019

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↗
USGS2019

Laboratory experiments of volcanic ash resuspension by wind

Fresh volcanic eruption deposits tend to be loose, bare, and readily resuspended by wind. Major resuspension events in Patagonia, Iceland, and Alaska have lofted ash clouds with potential to impact aircraft, infrastructure, and downwind communities. However, poor constraints on this resuspension process limit our ability to model this phenomenon. Here, we present laboratory experiments measuring threshold shear velocities and emission rates of resuspended ash under different environmental conditions, including relative humidity of 25–75% and simulated rainfall with subsequent drying. Eruption deposits were replicated using ash collected from two major eruptions: the 18 May 1980 eruption of Mount St. Helens and the 1912 eruption of Novarupta, in Alaska's Valley of Ten Thousand Smokes. Samples were conditioned in a laboratory chamber and prepared with bulk deposit densities of 1,300–1,500 kg/m 3 . A control sample of dune sand was included for comparison. The deposits were subjected to different wind speeds using a modified PI‐SWERL® instrument. Under a constant relative humidity of 50% and shear velocities 0.4–0.8 m/s, PM 10 emission by resuspension ranged from 10 to >100 mg·m −2 ·s −1 . Addition of liquid water equivalent to 5 mm of rainfall had little lasting effect on Mount St. Helens wind erosion potential, while the Valley of Ten Thousand Smokes deposits exhibited lower emissions for at least 12 days. The results indicate that particle resuspension due to wind erosion from ash deposits potentially exceeds that of most desert surfaces and approaches some of the highest emissions ever measured.

Journal of Geophysical Research - Atmospheres↗
USGS2018

Floor-fractured craters on Ceres and implications for interior processes

Several of the impact craters on Ceres have sets of fractures on their floors. These fractures appear similar to those found within a class of lunar craters referred to as floor-fractured craters (FFCs). We have cataloged the Ceres FFCs according to the classification scheme designed for the Moon. An analysis of the depth to diameter ratio for Ceres craters shows that, like lunar FFCs, the Ceres FFCs are anomalously shallow. Large (>50 km) Ceres FFCs are most consistent with Class 1 lunar FFCs, while smaller craters on Ceres are more consistent with Class 4 lunar FFCs. This suggests that Ceres FFCs may similarly be undergoing fracturing due to the intrusion of a low-density material below the craters. While on the Moon (and Mars) the intrusive material is hypothesized to be silicate magma, cryomagmatic intrusions are more likely responsible for the formation of the Ceres FFCs. However, new models suggest that at least some of the FFC fractures may have formed due to the solid-state flow of a low-viscosity and low-density material into the crater wall.

Journal of Geophysical Research↗
USGS2018

Submarine deposition of a subaerial landslide in Taan Fiord, Alaska

A large subaerial landslide entered Taan Fiord, Alaska, on 17 October 2015 producing a tsunami with runup to 193 m. We use LiDAR data to show the slide volume to be 76 + 3/−4 million cubic meters and that 51,000,000 m 3 entered Taan Fiord. In 2016, we mapped the fjord with multibeam bathymetry and high‐resolution seismic data. Landslide and postlandslide deposits extend 6 km downfjord, are up to 70 ± 11 m thick, and have a total volume of ~147,000,000 m 3 . Seismic data image a blocky landslide unit and two units deposited immediately after the landslide. The blocky landslide unit is ~65,000,000 m 3 . We infer it consists dominantly of subaerially derived material and secondarily of fjord floor sediment. The overlying units are likely megaturbidites presumably deposited within minutes to days after the landslide. We infer that these deposits dominantly consist of fjord floor material mobilized and suspended as the slide entered and traveled downfjord. The lower postlandslide unit is up to 35 ± 6 m thick, and the upper unit is up to 12 ± 3 m thick. These deposits are distinctive and will leave a lasting record of the event. This subaerial‐to‐submarine landslide deposit is distinct from other submarine landslide deposits studied in Alaskan fjords because it has a much greater thickness, larger and more angular blocks, distinctive postlandslide megaturbidites, and a higher‐amplitude acoustic signature of the blocky deposit. The tight constraints on the landslide source and deposit volumes, topography, bathymetry, and tsunami runup heights and flow directions should make this a benchmark site for landslide‐tsunami models.

Alaska↗
USGS2018

Salt marsh loss affects tides and sediment budget in shallow bays

The current paradigm is that salt marshes and their important ecosystem services are threatened by global climate change; indeed, large marsh losses have been documented worldwide. Morphological changes associated with salt marsh erosion are expected to influence the hydrodynamics and sediment dynamics of coastal systems. Here the influence of salt marsh erosion on the tidal hydrodynamics and sediment storage capability of shallow bays is investigated. Hydrodynamics, sediment transport, and vegetation dynamics are simulated using the numerical framework Coupled Ocean‐Atmosphere‐Wave‐Sediment Transport in the Barnegat Bay‐Little Egg Harbor system, USA. We show that salt marsh erosion influences the propagation of tides into back‐barrier basins, reducing the periodic inundation and sediment delivery to marsh platforms. As salt marshes erode, the sediment trapping potential of marsh platforms decreases exponentially. In this test case, up to 50% of the sediment mass trapped by vegetation is lost once a quarter of the marsh area is eroded. Similarly, without salt marshes the sediment budget of the entire bay significantly declines. Therefore, a positive feedback might be triggered such that as the salt marsh retreats the sediment storage capacity of the system declines, which could in turn further exacerbate marsh degradation.

New Jersey↗