Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “JGR - 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 55 records · Page 3Linked to original sources
USGS2022

Power-law viscoelastic flow of the lower accretionary prism in the Makran subduction zone following the 2013 Baluchistan Earthquake

Subduction zone accretionary prisms are commonly modeled as elastic structures where permanent deformation is accommodated by faulting and folding of otherwise elastic materials, yet accretionary prisms may exhibit other deformation styles over relatively short time scales. In this study, we use 6.5-year (2014–2021) Sentinel-1 interferometric synthetic aperture radar (InSAR) time-series of post-seismic deformation in the Makran accretionary prism of southeast Pakistan to characterize non-linear viscoelastic deformation within an active accretionary prism on short timescales (months to years). We constructed a series of 3-D finite-element models of the Makran subduction zone, including an accretionary prism, and constrained the elastic thickness of the upper wedge and the flow-law parameters (power-law exponent, activation enthalpy, and pre-exponential constant) of the lower wedge through forward model fits to the InSAR time-series. Our results show that the prism is elastically thin (8–12 km) and the non-linear viscoelastic relaxation of the deep portions of the prism alone can sufficiently explain the post-seismic surface deformation. Our best fitting flow-law parameters ( n = 3.76 ± 0.39, Q = 82.2 ± 37.73 kJ mol −1 , and A = 10 −3.36±4.69 ) are consistent with triggering of low temperature dislocation creep within fluid-saturated siliciclastic rocks. We believe that the fluids necessary for this weakening originate from sedimentary underplating and/or the presence the hydrocarbons. The presence of power-law rheology within the lower wedge impacts the estimated plate coupling and the stress state in the subduction system, with respect to the conventional elastic wedge model, and hence should to be considered in future earthquake cycle models.

Makran accretionary prism↗
USGS2022

Physical properties of the crust influence aftershock locations

Aftershocks do not uniformly surround a mainshock, and instead occur in spatial clusters. Spatially variable physical properties of the crust may influence the spatial distribution of aftershocks. I study four aftershock sequences in Southern California (1992 Landers, 1999 Hector Mine, 2010 El Mayor—Cucapah, and 2019 Ridgecrest) to investigate which physical properties are spatially correlated with aftershock occurrence. I find that aftershocks correlate with several properties, including measures of stress and stress change from the mainshock, fault structure, kinematics, seismic velocity, and heat flow. Aftershock spatial density exhibits an order of magnitude or more variation as a function of these properties. I determine simple empirical relations between each of the properties and the aftershock spatial density, and use these relations to construct new spatial models that describe aftershock locations. The new spatial models are a significant improvement over a simple base model, but do not fully capture the dense spatial clustering of aftershocks. Numerous spatially varying physical properties exhibit no (or poor) correlation with aftershock spatial density, including temperature, rock composition, and rheological properties that might be expected to control aftershock occurrence. These results suggest that while spatially variable physical properties appear to influence aftershock locations, more work is necessary in order to establish the connections between aftershock occurrence and the causative physical properties.

Baja California, California, Nevada↗
USGS2023

On the scale-dependence of fault surface roughness

Defining roughness as the ratio of height to length, the standard approach to characterize amplitudes of single fault, joint and fracture surfaces is to measure average height as a function of profile length. Empirically, this roughness depends strongly on scale. The ratio is approximately 0.01 at a few mm but 10× smaller at a few tens of meters. Surfaces are rougher at small scales. However, these conclusions are metric-dependent. If instead height is averaged over wavelength, roughness is nearly Brown spatial noise, having almost scale-independent apparent surface height to wavelength ratio. The small deviation from scale-independence is of the opposite sense than found using the standard metric; surfaces are slightly rougher at long wavelengths. Some natural surfaces may be Brownian within the measurement uncertainties. These contradictions are curiosities of surfaces that have Hurst exponents between 0.5 and 1, as natural fault surfaces do. The wavelength-based analysis of roughness and how it changes with scale are straight-forward; a normalized Fourier transform approximately preserves amplitude and its scale dependence in the wavelength domain. Among the conclusions from reconsideration of scale dependence are that the scale dependence is weak and much smaller than that of other fault and shear zone properties. Background and aftershock seismicity, jogs and step-overs indicate strong localization (smoothing) with slip and scale. The lack of strong scale dependence to surface roughness suggests it is not the dominant control on brittle shear zone evolution.

JGR Solid Earth↗
USGS2023

Microscale spatial variations in coseismic temperature rise on hematite fault mirrors in the Wasatch fault damage zone

Coseismic temperature rise activates fault dynamic weakening that promotes earthquake rupture propagation. The spatial scales over which peak temperatures vary on slip surfaces are challenging to identify in the rock record. We present microstructural observations and electron backscatter diffraction data from three small-displacement hematite-coated fault mirrors (FMs) in the Wasatch fault damage zone, Utah, to evaluate relations between fault properties, strain localization, temperature rise, and weakening mechanisms during FM development. Millimeter- to cm-thick, matrix-supported, hematite-cemented breccia is cut by ∼25–200 μm-thick, texturally heterogeneous veins that form the hematite FM volume (FMV). Grain morphologies and textures vary with FMV thickness over μm to mm lengthscales. Cataclasite grades to ultracataclasite where FMV thickness is greatest. Thinner FMVs and geometric asperities are characterized by particles with subgrains, serrated grain boundaries, and(or) low-strain polygonal grains that increase in size with proximity to the FM surface. Comparison to prior hematite deformation experiments suggests FM temperatures broadly range from ≥400°C to ≥800–1100°C, compatible with observed coeval brittle and plastic deformation mechanisms, over sub-mm scales on individual slip surfaces during seismic slip. We present a model of FM development by episodic hematite precipitation, fault reactivation, and strain localization, where the thickness of hematite veins controls the width of the deforming zones during subsequent fault slip, facilitating temperature rise and thermally activated weakening. Our data document intrasample coseismic temperatures, resultant deformation and dynamic weakening mechanisms, and the length scales over which these vary on slip surfaces.

Utah↗
USGS2023

A specific earthquake processing workflow for studying long-lived, explosive volcanic eruptions with application to the 2008 Okmok Volcano, Alaska, eruption

By providing unrivaled resolution in both time and space, volcano seismicity helps to chronicle and interpret eruptions. Standard earthquake detection methods are often insufficient as the eruption itself produces continuous seismic waves that obscure earthquake signals. We address this problem by developing an earthquake processing workflow specific to a high-noise volcanic environment and applying it to the explosive 2008 Okmok Volcano eruption. This process includes applying single-channel template matching combined with machine-learning and fingerprint-based techniques to expand the existing earthquake catalog of the eruption. We detected an order of magnitude more earthquakes, then located, relocated, determined locally calibrated magnitudes, and classified the events in the enhanced catalog. This new high-resolution earthquake catalog increases the number of observations by about a factor of 10 and enables the detailed spatiotemporal seismic analysis during a large eruption.

Alaska↗
USGS2022

Wave-driven hydrodynamic processes over fringing reefs with varying slopes, depths, and roughness: Implications for coastal protection

Wave breaking on the steep fore-reef slopes of shallow fringing reefs is effective at dissipating incident sea-swell waves prior to reaching reef shorelines. However, wave setup and free infragravity waves generated during the sea-swell breaking process are often the largest contributors to wave-driven water levels at the shoreline. Laboratory flume experiments and a multi-layer phase-resolving nonhydrostatic wave-flow model, which includes a canopy model to predict drag forces generated by roughness elements, were used to investigate the wave-driven water levels on fringing reefs. Though the model is capable of three dimensional simulations, consistent with the laboratory study, a two-dimensional vertical mode was used. In contrast to many previous studies, both the laboratory experiment and the numerical model account for the effects of large bottom roughness. The numerical model reproduced the observations of the wave transformation and runup over both smooth and rough reef profiles. The numerical model was then extended to quantify the influence of reef geometry and compared to simulations of plane beaches lacking a reef. For a set offshore forcing condition, the fore-reef slope controlled wave runup on reef fronted beaches, whereas the beach slope controlled wave runup on plane beaches. As a result, the coastal protection utility of reefs is dependent on these slopes. For our examples, with a fore-reef slope of 1/5 and a 500 m prototype reef flat length, a beach slope of ~1/30 marked the transition between the reef providing runup reduction for steeper beach slopes and enhancing wave runup for milder slopes. Roughness coverage, spacing, dimensions, and drag coefficient were investigated with results indicating the greatest runup reductions were due to tall roughness elements on the reef flat.

JGR Oceans↗
USGS2023

Dynamics of the wave-driven circulation in the lee of nearshore reefs

Nearshore rocky reefs with scales of order 10–100 m are common along the world's coastline and often shape wave-driven hydrodynamics and shoreline morphology in their lee. The interaction of waves with these reefs generally results in either two or four-cell mean circulation systems (2CC and 4CC, respectively), with diverging flows behind the reefs and at the shoreline in the 2CC case and flows that diverge in the lee and converge at the shoreline in the 4CC case. By applying a phase-resolving wave-flow model to conduct a detailed analysis of mean momentum balances for waves interacting with nearshore reefs, we develop an understanding of the drivers of 2CC and 4CC flow dynamics and how they vary for different reef geometries and wave and water level conditions. The 2CC or 4CC patterns were primarily driven by alongshore pressure gradients toward the exposed (nonreef fronted) or reef-fronted beach. These alongshore pressure gradients were dependent on the cross-shore setup dynamics governed by the balance between pressure (i.e., related to the setup) and radiation stress gradients, and mean bottom stresses exerted on the water column. If shoreline wave setup in the lee of the reef was less than the exposed beach, a 4CC pattern developed with convergent flow at the shoreline in the lee of the reef; otherwise, a 2CC emerged with divergent flow at the shoreline. Across the parameter space investigated, reef roughness, distance to the shoreline, and beach slope were the three parameters most likely to change the flow patterns between 2CC and 4CC.

JGR-Oceans↗
USGS2022

The Curiosity Rover’s exploration of Glen Torridon, Gale crater, Mars: An overview of the campaign and scientific results

The Mars Science Laboratory rover, Curiosity , explored the clay mineral-bearing Glen Torridon region for one martian year between January 2019 and January 2021, including a short campaign onto the Greenheugh pediment. The Glen Torridon campaign sought to characterize the geology of the area, seek evidence of habitable environments, and document the onset of a potentially global climatic transition during the Hesperian era. Curiosity roved 5 km in total throughout Glen Torridon, from the Vera Rubin ridge to the northern margin of the Greenheugh pediment. Curiosity acquired samples from 11 drill holes during this campaign and conducted the first martian thermochemolytic-based organics detection experiment with the Sample Analysis at Mars instrument suite. The lowest elevations within Glen Torridon represent a continuation of lacustrine Murray formation deposits, but overlying widespread cross bedded sandstones indicate an interval of more energetic fluvial environments and prompted the definition of a new stratigraphic formation in the Mount Sharp group called the Carolyn Shoemaker formation. Glen Torridon hosts abundant phyllosilicates yet remains compositionally and mineralogically comparable to the rest of the Mount Sharp group. Glen Torridon samples have a great diversity and abundance of sulfur-bearing organic molecules, which are consistent with the presence of ancient refractory organic matter. The Glen Torridon region experienced heterogeneous diagenesis, with the most striking alteration occurring just below the Siccar Point unconformity at the Greenheugh pediment. Results from the pediment campaign show that the capping sandstone formed within the Stimson Hesperian aeolian sand sea that experienced seasonal variations in wind direction.

JGR Planets↗
USGS2022

Fundamental science and engineering questions in planetary cave exploration

Nearly half a century ago, two papers postulated the likelihood of lunar lava tube caves using mathematical models. Today, armed with an array of orbiting and fly-by satellites and survey instrumentation, we have now acquired cave data across our solar system—including the identification of potential cave entrances on the Moon, Mars, and at least six other planetary bodies. These discoveries gave rise to the study of planetary caves. To help advance this field, we leveraged the expertise of an interdisciplinary group to identify a strategy to explore caves beyond Earth. Focusing primarily on astrobiology, the cave environment, geology, robotics, instrumentation, and human exploration, our goal was to produce a framework to guide this subdiscipline through at least the next decade. To do this, we first assembled a list of 198 science and engineering questions. Then, through a series of social surveys, 114 scientists and engineers winnowed down the list to the top 53 highest priority questions. This exercise resulted in identifying emerging and crucial research areas that require robust development to ultimately support a robotic mission to a planetary cave—principally the Moon and/or Mars. With the necessary financial investment and institutional support, the research and technological development required to achieve these necessary advancements over the next decade are attainable. Subsequently, we will be positioned to robotically examine lunar caves and search for evidence of life within martian caves; in turn, this will set the stage for human exploration and potential habitation of both the lunar and martian subsurface.

JGR Planets↗
USGS2022

Planetary Aeolian landforms: An introduction to the Fifth Planetary Dunes Workshop Special Issue

Aeolian landforms are widespread in our solar system. Understanding the exact nature and processes of formation of these features are challenging tasks necessitating a strong collaboration between scientists with different skills and scientific backgrounds. This paper describes the special issue for the 5 th International Planetary Dunes Workshop, which includes 15 research papers and three commentaries. Among the 18 papers included in this collection, 16 cover Martian aeolian science and two Titan aeolian science. The papers presented focus on bedform morphology and dynamics via remote sensing data, modelling, analogues studies and laboratory experiments. Here we put the main results of the papers in their appropriate scientific context and discuss potential future lines of research.

JGR Planets↗
USGS2022

Planetary caves: A Solar System view of processes and products

We provide the first solar system wide compendium of speleogenic processes and products. An examination of 15 solar system bodies revealed that six cave-forming processes occur beyond Earth including volcanic (cryo and magmatic), fracturing (tectonic and impact melt), dissolution, sublimation, suffusion, and landslides. Although no caves (i.e., confirmed entrances with associated linear passages) have been confirmed, 3,545 SAPs (subsurface access points) have been identified on 11 planetary bodies and the potential for speleogenic processes (and thus SAPs) was observed on an additional four planetary bodies. The bulk of our knowledge on extraterrestrial SAPs is based upon global databases for the Moon and Mars, which are bodies for which high-resolution imagery and other data are available. To further characterize most of the features beyond the Moon and Mars, acquisition (preferably global coverage) and subsequent analysis of high-resolution imagery will be required. The next few decades hold considerable promise for identifying and potentially characterizing caves across the solar system.

JGR Planets↗
USGS2023

Regolith of the crater floor units, Jezero crater, Mars: Textures, composition and implications for provenance

A multi-instrument study of the regolith of Jezero crater floor units by the Perseverance rover has identified three types of regolith: fine-grained, coarse-grained, and mixed-type. Mastcam-Z, WATSON, and SuperCam RMI were used to characterize regolith texture, particle size, and roundedness where possible. Mastcam-Z multispectral and SuperCam LIBS data were used to constrain the composition of the regolith types. Fine-grained regolith is found surrounding bedrock and boulders, comprising bedforms, and accumulating on top of rocks in erosional depressions. Spectral and chemical data show it is compositionally consistent with pyroxene and a ferric-oxide phase. Coarse-grained regolith consists of 1-2 mm well-sorted gray grains that are found concentrated around the base of boulders and bedrock, and armoring bedforms. Its chemistry and spectra indicate it is olivine-bearing, and spatial distribution and roundedness indicate it has been transported, likely by saltation-induced creep. Coarse grains share similarities with the olivine grains observed in the Séítah formation bedrock, making that unit a possible source for these grains. Mixed-type regolith contains fine- and coarse-grained regolith components and larger rock fragments. The rock fragments are texturally and spectrally similar to bedrock within the Máaz and Séítah formations, indicating origins by erosion from those units, although they could also be a lag deposit from erosion of an overlying unit. The fine- and coarse-grained types are compared to their counterparts at other landing sites to inform global, regional, and local inputs to regolith formation within Jezero crater. The regolith characterization presented here informs regolith sampling efforts underway by Perseverance.

JGR-Planets↗
USGS2022

Physical controls on the hydrology of perennially ice-covered lakes, Taylor Valley, Antarctica (1996-2013)

The McMurdo Dry Valleys, Antarctica, are a polar desert populated with numerous closed-watershed, perennially ice-covered lakes primarily fed by glacial melt. Lake levels have varied by as much as 8 m since 1972 and are currently rising after a decade of decreasing. Precipitation falls as snow, so lake hydrology is dominated by energy available to melt glacier ice and to sublimate lake ice. To understand the energy and hydrologic controls on lake level changes and to explain the variability between neighboring lakes, only a few kilometers apart, we model the hydrology for the three largest lakes in Taylor Valley. We apply a physically based hydrological model that includes a surface energy balance model to estimate glacial melt and lake sublimation to constrain mass fluxes to and from the lakes. Results show that lake levels are very sensitive to small changes in glacier albedo, air temperature, and wind speed. We were able to balance the hydrologic budget in two watersheds using meltwater inflow and sublimation loss from the ice-covered lake alone. A third watershed, closest to the coast, required additional inflow beyond model uncertainties. We hypothesize a shallow groundwater system within the active layer, fed by dispersed snow patches, contributes 23% of the inflow to this watershed. The lakes are out of equilibrium with the current climate. If the climate of our study period (1996-2013) persists into the future, the lakes will reach equilibrium starting in 2300, with levels 2-17 m higher, depending on the lake, relative to the 2020 level.

JGR Earth Surface↗
USGS2023

Sound-side inundation and seaward erosion of a barrier island during hurricane landfall

Barrier islands are especially vulnerable to hurricanes and other large storms, owing to their mobile composition, low elevations, and detachment from the mainland. Conceptual models of barrier-island evolution emphasize ocean-side processes that drive landward migration through overwash, inlet migration, and aeolian transport. In contrast, we found that the impact of Hurricane Dorian (2019) on North Core Banks, a 36-km barrier island on the Outer Banks of North Carolina, was primarily driven by inundation of the island from Pamlico Sound, as evidenced by storm-surge model results and observations of high-water marks and wrack lines. Analysis of photogrammetry products from aerial imagery collected before and after the storm indicate the loss of about 18% of the subaerial volume of the island through the formation of over 80 erosional washout channels extending from the marsh and washover platform, through gaps in the foredunes, to the shoreline. The washout channels were largely co-located with washover fans deposited by earlier events. Net seaward export of sediment resulted in the formation of deltaic bars offshore of the channels, which became part of the post-storm berm recovery by onshore bar migration and partial filling of the washouts with washover deposits within 2 months. This event represents a volumetric setback in the overwash/rollover behavior required for barrier transgression, but the new ponds and lowland habitats may provide beneficial habit for endangered species and will likely persist for years.

North Carolina↗
USGS2023

A model integrating satellite-derived shoreline observations for predicting fine-scale shoreline response to waves and sea-level rise across large coastal regions

Satellite-derived shoreline observations combined with dynamic shoreline models enable fine-scale predictions of coastal change across large spatiotemporal scales. Here, we present a satellite-data-assimilated, “littoral-cell”-based, ensemble Kalman-filter shoreline model to predict coastal change and uncertainty due to waves, sea-level rise (SLR), and other natural and anthropogenic processes. We apply the developed ensemble model to the entire California coastline (approximately 1,760 km), much of which is sparsely monitored with traditional survey methods (e.g., Lidar/GPS). Water-level-corrected, satellite-derived shoreline observations (obtained from the CoastSat toolbox) offer a nearly unbiased representation of in situ surveyed shorelines (e.g., mean sea-level elevation contours) at Ocean Beach, San Francisco. We demonstrate that model calibration with satellite observations during a 20-year hindcast period (1995–2015) provides nearly equivalent model forecast accuracy during a validation period (2015–2020) compared to model calibration with monthly in situ observations at Ocean Beach. When comparing model-predicted shoreline positions to satellite-derived observations, the model achieves an accuracy of <10 m RMSE for nearly half of the entire California coastline for the validation period. The calibrated/validated model is then applied for multi-decadal simulations of shoreline change due to projected wave and sea-level conditions, while holding the model parameters fixed. By 2100, the model estimates that 24%–75% of California's beaches may become completely eroded due to SLR scenarios of 1.0–3.0 m, respectively. The satellite-data-assimilated modeling system presented here is generally applicable to a variety of coastal settings around the world owing to the global coverage of satellite imagery.

JGR Earth Surface↗
USGS2023

Debris-flow process controls on steepland morphology in the San Gabriel Mountains, California

Steep landscapes evolve largely by debris flows, in addition to fluvial and hillslope processes. Abundant field observations document that debris flows incise valley bottoms and transport substantial sediment volumes, yet their contributions to steepland morphology remain uncertain. This has, in turn, limited the development of debris-flow incision rate formulations that produce morphology consistent with natural landscapes. In many landscapes, including the San Gabriel Mountains (SGM), California, steady-state fluvial channel longitudinal profiles are concave-up and exhibit a power-law relationship between channel slope and drainage area. At low drainage areas, however, valley slopes become nearly constant. These topographic forms result in a characteristically curved slope-area signature in log-log space. Here, we use a one-dimensional landform evolution model that incorporates debris-flow erosion to reproduce the relationship between this curved slope-area signature and erosion rate in the SGM. Topographic analysis indicates that the drainage area at which steepland valleys transition to fluvial channels correlates with measured erosion rates in the SGM, and our model results reproduce these relationships. Further, the model only produces realistic valley profiles when parameters that dictate the relationship between debris-flow erosion, valley-bottom slope, and debris-flow depth are within a narrow range. This result helps place constraints on the mathematical form of a debris-flow incision law. Finally, modeled fluvial incision outpaces debris-flow erosion at drainage areas less than those at which valleys morphologically transition from near-invariant slopes to concave profiles. This result emphasizes the critical role of debris-flow incision for setting steepland form, even as fluvial incision becomes the dominant incisional process.

California↗
USGS2023

Characterizing performance of freshwater wetland methane models across time scales at FLUXNET-CH4 sites using wavelet analyses

Process-based land surface models are important tools for estimating global wetland methane (CH 4 ) emissions and projecting their behavior across space and time. So far there are no performance assessments of model responses to drivers at multiple time scales. In this study, we apply wavelet analysis to identify the dominant time scales contributing to model uncertainty in the frequency domain. We evaluate seven wetland models at 23 eddy covariance tower sites. Our study first characterizes site-level patterns of freshwater wetland CH 4 fluxes (FCH 4 ) at different time scales. A Monte Carlo approach was developed to incorporate flux observation error to avoid misidentification of the time scales that dominate model error. Our results suggest that (a) significant model-observation disagreements are mainly at multi-day time scales (<15 days); (b) most of the models can capture the CH 4 variability at monthly and seasonal time scales (>32 days) for the boreal and Arctic tundra wetland sites but have significant bias in variability at seasonal time scales for temperate and tropical/subtropical sites; (c) model errors exhibit increasing power spectrum as time scale increases, indicating that biases at time scales <5 days could contribute to persistent systematic biases on longer time scales; and (d) differences in error pattern are related to model structure (e.g., proxy of CH 4 production). Our evaluation suggests the need to accurately replicate FCH 4 variability, especially at short time scales, in future wetland CH 4 model developments.

JGR Biogeosciences↗
USGS2023

Scaling microseismic cloud shape during hydraulic stimulation using in-situ stress and permeability

Forecasting microseismic cloud shape as a proxy of stimulated rock volume may improve the design of an energy extraction system. The microseismic cloud created during hydraulic stimulation of geothermal reservoirs is known empirically to extend in the general direction of the maximum principal stress. However, this empirical relationship is often inconsistent with reported results, and the cloud growth process remains poorly understood. This study investigates microseismic cloud growth using data obtained from a hydraulic stimulation project in Basel, Switzerland, and explores its correlation with measured in situ stress. We applied principal component analysis to a time series of microseismicity for macroscopic characterization of microseismic cloud growth in two- and three-dimensional space. The microseismic cloud, in addition to extending in the general direction of maximum principal stress, expanded in the direction of intermediate principal stress. The orientation of the least microseismic cloud growth was stable and almost identical to the minimum principal stress direction. Further, microseismic cloud shape ratios showed good agreement when compared with in situ stress magnitude ratios. The permeability tensor estimated from microseismicity also provided a good correlation in terms of direction and magnitude with the microseismic cloud growth. We show that in situ stress plays a dominant role by controlling the permeability of each existing fracture in the reservoir fracture system. Consequently, microseismic cloud growth can be scaled by in situ stress as a first-order approximation if there is sufficient variation in the orientation of existing faults.

Basel↗