Geology Reports⌕ Search

SEARCH · Geology Reports

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 901 records · Page 50Linked to original sources
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

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

Degradation of 100‐m‐scale rocky ejecta craters at the InSight Landing Site on Mars and implications for surface processes and erosion rates in the hesperian and amazonian

Rocky ejecta craters (RECs) at the Interior Exploration using Seismic Investigations, Geodesy and Heat Transport (InSight) landing site on Elysium Planitia, Mars, provide constraints on crater modification and rates for the Hesperian and Amazonian. The RECs are between 10 m and 1.2 km in diameter and exhibit five classes of preservation. Class 1 represents pristine craters with sharp rims and abundant ejected rocks. From Classes 2 to 5, rims become more subdued, craters are infilled, and the ejecta become discontinuously distributed. High‐Resolution Imaging Science Experiment digital elevation models indicate a maximum depth to diameter ratio of ~0.15, which is lower than pristine models for craters of similar size. The low ratio is related to the presence of a loosely consolidated regolith and early‐stage eolian infill. Rim heights have an average height to diameter ratio of ~0.03 for the most pristine class. The size‐frequency distribution of RECs, plotted using cumulative and differential methods, indicates that crater classes within the diameter range of 200 m to 1.2 km are separated by ~100 to 200 Myr. Smaller craters degrade faster, with classes separated by <100 Myr. Rim erosion can be entirely modeled by nonlinear diffusional processes using the calculated timescales and a constant diffusivity of 8 × 10 −7 m 2 /year for craters 200 to 500 m in diameter. Diffusion models only partly capture depth‐related degradation, which requires eolian infill. Depth degradation and rim erosion rates are 10 −2 to 10 −3 m/Myr, respectively. The rates are consistent with relatively slow modification that is typical of the last two epochs of Martian history.

Journal of Geophysical Research - Planets↗
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↗
USGS2019

Overview of spirit microscopic imager results

This paper provides an overview of Mars Exploration Rover Spirit Microscopic Imager (MI) operations and the calibration, processing, and analysis of MI data. The focus of this overview is on the last five Earth years (2005–2010) of Spirit 's mission in Gusev crater, supplementing the previous overview of the first 450 sols of the Spirit MI investigation. Updates to radiometric calibration using in‐flight data and improvements in high‐level processing are summarized. Released data products are described, and a table of MI observations, including target/feature names and associated data sets, is appended. The MI observed natural and disturbed exposures of rocks and soils as well as magnets and other rover hardware. These hand‐lens‐scale observations have provided key constraints on interpretations of the formation and geologic history of features, rocks, and soils examined by Spirit . MI images complement observations by other Spirit instruments, and together show that impact and volcanic processes have dominated the origin and evolution of the rocks in Gusev crater, with aqueous activity indicated by the presence of silica‐rich rocks and sulfate‐rich soils. The textures of some of the silica‐rich rocks are similar to terrestrial hot spring deposits, and observations of subsurface cemented layers indicate recent aqueous mobilization of sulfates in places. Wind action has recently modified soils and abraded many of the rocks imaged by the MI, as observed at other Mars landing sites.

Journal of Geophysical Research E: Planets↗
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

Geomorphic evolution of a gravel‐bed river under sediment‐starved vs. sediment‐rich conditions: River response to the world's largest dam removal

Understanding river response to sediment pulses is a fundamental problem in geomorphic process studies, with myriad implications for river management. However, because large sediment pulses are rare and usually unanticipated, they are seldom studied at field scale. We examine fluvial response to a massive (~20 Mt) sediment pulse released by the largest dam removal globally, on the Elwha River, Washington, United States, in an 11‐year before‐after/control‐impact study of channel morphology and grain size. We test the hypothesis that for a given flow magnitude, greater geomorphic change occurs under sediment‐rich conditions than under sediment‐starved conditions. Channel response to flow forcing was significantly different during the sediment‐pulse peak, 1–2 years after dam removal began, than earlier or later. During peak sediment supply our hypothesis was supported; major geomorphic change occurred under low flows and unit stream power ≤60 W/m 2 . However, by 4–6 years after dam removal began, rates of geomorphic change and sensitivity to stream power had decreased substantially such that our hypothesis was no longer unequivocally supported. These findings are consistent with a two‐phase conceptual model of dam‐removal response, involving a transport‐limited state followed by a more supply‐limited state. From comparisons with other dam removals and natural sediment pulses, we infer that the longevity of sediment‐pulse signals in gravel‐bed rivers depends upon gradient, river discharge, valley morphology, and sediment grain size. Stream power associated with substantial geomorphic change varies with sediment supply, such that assigning a general threshold stream power to gravel‐bed rivers may be untenable.

Washington↗
USGS2019

Downstream‐propagating channel responses to decadal‐scale climate variability in a glaciated river basin

Regional climate is an important control on the rate of coarse sediment mobilization and transport in alpine river systems. Changes in climate are then expected to cause a cascade of geomorphic responses, including adjustments in downstream channel morphology. However, the mechanics and sensitivity of channel response to short‐term climate variability remain poorly documented. In the Nooksack River, which drains a glaciated stratovolcano in Washington State, bed elevation changes were inferred from shifting stage–discharge relations at seven USGS stream gages. Decadal‐scale elevation trends at most sites can be explained as a downstream‐propagating channel response to regional climate variability, where periods of persistent warm, dry [cool, wet] conditions corresponded to periods of aggradation [incision]. The channel elevation response propagated downstream at a rate of one to four kilometers per year; propagation rate scaled closely with channel slope. Historical trends in glacier extent and flood intensity both show some potential to explain climate–sediment linkages, though assessing causation is complicated by the shared climate signal in both records. Results show the influence of the Pacific Decadal Oscillation, with relatively high coarse sediment yields prior to 1950 and since 1980, and notably lower sediment yields from 1950 to 1980. Measured sediment yields from nearby glaciated basins corroborate this history, suggesting a regional coherence to these climate–sediment linkages. These results document consistent relations between climate, sediment supply and downstream channel response at the basin‐scale, with channel responses propagating downstream over periods of decades with little apparent attenuation.

Washington↗
USGS2019

Prediction of ice‐free conditions for a perennially ice‐covered Antarctic lake

Although perennially ice‐covered Antarctic lakes have experienced variable ice thicknesses over the past several decades, future ice thickness trends and associated aquatic biological responses under projected global warming remain unknown. Heat stored in the water column in chemically stratified Antarctic lakes that have middepth temperature maxima can significantly influence the ice thickness trends via upward heat flux to the ice/water interface. We modeled the ice thickness of the west lobe of Lake Bonney, Antarctica, based on possible future climate scenarios utilizing a 1D thermodynamic model that accounts for surface radiative fluxes as well as the heat flux associated with the temperature evolution of the water column. Model results predict that the ice cover of Lake Bonney will shift from perennial to seasonal within one to four decades, a change that will drastically influence ecosystem processes within the lake.

Journal of Geophysical Research F: Earth Surface↗