Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Interfaces”

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

Xenon hydrate as an analogue of methane hydrate in geologic systems out of thermodynamic equilibrium

Methane hydrate occurs naturally under pressure and temperature conditions that are not straightforward to replicate experimentally. Xenon has emerged as an attractive laboratory alternative to methane for studying hydrate formation and dissociation in multiphase systems, given that it forms hydrates under milder conditions. However, building reliable analogies between the two hydrates requires systematic comparisons, which are currently lacking. We address this gap by developing a theoretical and computational model of gas hydrates under equilibrium and nonequilibrium conditions. We first compare equilibrium phase behaviors of the Xe·H 2 O and CH 4 ·H 2 O systems by calculating their isobaric phase diagram, and then study the nonequilibrium kinetics of interfacial hydrate growth using a phase field model. Our results show that Xe·H 2 O is a good experimental analog to CH 4 ·H 2 O, but there are key differences to consider. In particular, the aqueous solubility of xenon is altered by the presence of hydrate, similar to what is observed for methane; but xenon is consistently less soluble than methane. Xenon hydrate has a wider nonstoichiometry region, which could lead to a thicker hydrate layer at the gas‐liquid interface when grown under similar kinetic forcing conditions. For both systems, our numerical calculations reveal that hydrate nonstoichiometry coupled with hydrate formation dynamics leads to a compositional gradient across the hydrate layer, where the stoichiometric ratio increases from the gas‐facing side to the liquid‐facing side. Our analysis suggests that accurate composition measurements could be used to infer the kinetic history of hydrate formation in natural settings where gas is abundant.

Geochemistry, Geophysics, Geosystems↗
USGS2020

Semiautomated estimates of directivity and related source properties of small to moderate southern California earthquakes using second seismic moments

We develop a semiautomated method for estimating with second seismic moments the directivity, rupture area, duration, and centroid velocity of earthquakes. The method is applied to 41 southern California earthquakes with magnitude in the range 3.5–5.2 and provides stable results for 28 events. Apparent source time functions (ASTFs) of P and S phases are derived using deconvolution with three stacked empirical Green's functions (seGf). The use of seGf suppresses nongeneric source effects, improves the focal mechanism correspondence to the analyzed earthquakes, and typically allows inclusion of 5 to 15 more ASTFs compared with analysis using a single eGf. Most analyzed earthquakes in the Trifurcation area of the San Jacinto Fault have directivities toward the northwest, while events around Cajon Pass and San Gabriel Mountain tend to propagate toward the southeast. These results are generally consistent with predictions for dynamic rupture on bimaterial interfaces associated with the imaged velocity contrasts in the area. The second moment inversions also provide constraints on the upper and lower bounds of rupture areas in our data set. Stress drops and uncertainties are estimated for elliptical ruptures using the derived characteristic rupture length and width. The semiautomated second moment method with seGfs can be used for routine application to moderate earthquakes in locations with good station coverage.

California↗
USGS2020

A one‐dimensional model for turbulent mixing in the benthic biolayer of stream and coastal sediments

In this paper, we develop and validate a rigorous modeling framework, based on Duhamel's Theorem, for the unsteady one-dimensional vertical transport of a solute across a flat sediment-water interface (SWI) and through the benthic biolayer of a turbulent stream. The modeling framework is novel in capturing the two-way coupling between evolving solute concentrations above and below the SWI and in allowing for a depth-varying diffusivity. Three diffusivity profiles within the sediment (constant, exponentially decaying, and a hybrid model) are evaluated against an extensive set of previously published laboratory measurements of turbulent mass transfer across the SWI. The exponential diffusivity profile best represents experimental observations and its reference diffusivity scales with the permeability Reynolds number, a dimensionless measure of turbulence at the SWI. The depth over which turbulence-enhanced diffusivity decays is of the order of centimeters and comparable to the thickness of the benthic biolayer. Thus, turbulent mixing across the SWI may serve as a universal transport mechanism, supplying the nutrient and energy fluxes needed to sustain microbial growth, and nutrient processing, in the benthic biolayer of stream and coastal sediments.

Water Resources Research↗
USGS2020

Conceptual model for the removal of cold-trapped H2O ice on the Mars northern seasonal springtime polar cap

The transport of H 2 O ice along the retreating north polar seasonal CO 2 ice cap has previously been modeled and observed. Spectral observations show that H 2 O ice forms on the interior of the seasonal cap, while thermal observations show these regions to be consistent with CO 2 ice. Prior to the sublimation of the seasonal CO 2 , the observed H 2 O ice deposits are diminished—and because H 2 O ice sublimation rates are extremely slow while in direct thermal contact with CO 2 ice, an alternate removal process must be operating. We propose a model where the process of removing these H 2 O deposits starts with insolation‐induced basal sublimation of the underlying CO 2 ice. This sublimed gas would “seep” upward and into the interface between the two ices, increasing pressure until the gas pressure fractures the cold‐trapped H 2 O ice. Small fragments would be suspended while larger fragments would be pushed aside, exposing the underlying CO 2 ice.

Geophysical Research Letters↗
USGS2021

A unified model of crustal stress heterogeneity from borehole breakouts and earthquake focal mechanisms

Observations of crustal stress orientation from the regional inversion of earthquake focal mechanisms often conflict with those from borehole breakouts, possibly indicating local stress heterogeneity, either laterally or with depth. To investigate this heterogeneity, we compiled SHmax estimates from previous studies for 57 near‐vertical boreholes with measured breakout azimuths across the Los Angeles region. We identified subsets of earthquake focal mechanisms from established earthquake catalogs centered around each borehole with various criteria for maximum depth and maximum lateral distance from the borehole. Each subset was independently inverted for 3‐D stress orientation and corresponding SHmax probability distributions, then compared with the corresponding borehole breakout‐derived estimate. We find good agreement when both methods sample the basement stress (breakouts are close to the sediment‐basement interface), or when both methods sample the mid‐basin stress (sufficient earthquakes are present within a sedimentary basin). Along sedimentary basin margins, in contrast, we find acceptable agreement only when focal mechanisms are limited to shallow and close earthquakes, implying short‐length‐scale heterogeneity of <20 km. While the region as a whole shows evidence of both lateral and vertical stress orientation heterogeneity, we find a more homogeneous stress state within basement rock, over length scales of 1–35 km. These results reconcile the apparently conflicting observations of short‐length‐scale heterogeneity observed in boreholes, which sample primarily the basins, with the relative homogeneity of stress inferred from focal mechanisms, which sample primarily the basement, and imply distinct regimes for the appropriate use of each type of stress indicator.

JGR Solid Earth↗
USGS2020

Modeling water quality in watersheds: From here to the next generation

In this synthesis, we assess present research and anticipate future development needs in modeling water quality in watersheds. We first discuss areas of potential improvement in the representation of freshwater systems pertaining to water quality, including representation of environmental interfaces, in‐stream water quality and process interactions, soil health and land management, and (peri‐)urban areas. In addition, we provide insights into the contemporary challenges in the practices of watershed water quality modeling, including quality control of monitoring data, model parameterization and calibration, uncertainty management, scale mismatches, and provisioning of modeling tools. Finally, we make three recommendations to provide a path forward for improving watershed water quality modeling science, infrastructure, and practices. These include building stronger collaborations between experimentalists and modelers, bridging gaps between modelers and stakeholders, and cultivating and applying procedural knowledge to better govern and support water quality modeling processes within organizations.

Water Resources Research↗
USGS2020

Unifying advective and diffusive descriptions of bedform pumping in the benthic biolayer of streams

Many water quality and ecosystem functions performed by streams occur in the benthic biolayer, the biologically active upper (~5 cm) layer of the streambed. Solute transport through the benthic biolayer is facilitated by bedform pumping, a physical process in which dynamic and static pressure variations over the surface of stationary bedforms (e.g., ripples and dunes) drive flow across the sediment-water interface. In this paper we derive two predictive modeling frameworks, one advective and the other diffusive, for solute transport through the benthic biolayer by bedform pumping. Both frameworks closely reproduce patterns and rates of bedform pumping previously measured in the laboratory, provided that the diffusion model's dispersion coefficient declines exponentially with depth. They are also functionally equivalent, such that parameter sets inferred from the 2D advective model can be applied to the 1D diffusive model, and vice versa. The functional equivalence and complementary strengths of these two models expand the range of questions that can be answered, for example, by adopting the 2D advective model to study the effects of geomorphic processes (such as bedform adjustments to land use change) on flow-dependent processes and the 1D diffusive model to study problems where multiple transport mechanisms combine (such as bedform pumping and turbulent diffusion). By unifying 2D advective and 1D diffusive descriptions of bedform pumping, our analytical results provide a straightforward and computationally efficient approach for predicting, and better understanding, solute transport in the benthic biolayer of streams and coastal sediments.

Water Resources Research↗
USGS2021

Labeling poststorm coastal imagery for machine learning: Measurement of interrater agreement

Classifying images using supervised machine learning (ML) relies on labeled training data—classes or text descriptions, for example, associated with each image. Data-driven models are only as good as the data used for training, and this points to the importance of high-quality labeled data for developing a ML model that has predictive skill. Labeling data is typically a time-consuming, manual process. Here, we investigate the process of labeling data, with a specific focus on coastal aerial imagery captured in the wake of hurricanes that affected the Atlantic and Gulf Coasts of the United States. The imagery data set is a rich observational record of storm impacts and coastal change, but the imagery requires labeling to render that information accessible. We created an online interface that served labelers a stream of images and a fixed set of questions. A total of 1,600 images were labeled by at least two or as many as seven coastal scientists. We used the resulting data set to investigate interrater agreement: the extent to which labelers labeled each image similarly. Interrater agreement scores, assessed with percent agreement and Krippendorff's alpha, are higher when the questions posed to labelers are relatively simple, when the labelers are provided with a user manual, and when images are smaller. Experiments in interrater agreement point toward the benefit of multiple labelers for understanding the uncertainty in labeling data for machine learning research.

Earth and Space Science↗
USGS2022

Interaction between climate and tectonics in the northern Lesser Antilles inferred from the last interglacial shoreline on Barbuda island

In the context of increasing evidence of plate interface coupling variability in subduction zones, there is a need to extend the short time window given by instrumental data and to gather data over multiple time and spatial scales. We hence investigated the long-term topography on Barbuda island, located in the northern part of the Lesser Antilles, west of the Caribbean subduction zone. Following pioneering work using a set of marine terraces on the eastern side of the island, we performed the first U-Th dating on 10 corals in growth position from the lowest terrace, for which the highest relative sea-level (RSL) indicator is found at 9 ± 1 m above the mean sea level. We find that this terrace corresponds to the Last Interglacial (LIG) (ages between 122.8 ± 0.3 ka and 128.1 ± 0.3 ka) and we estimate a paleo RSL of 7 ± 2 m above the current mean sea level. The present elevation of the LIG shoreline on Barbuda might imply tectonics as an additional mechanism to eustatic sea level, mantle dynamic topography and glacial isostatic adjustment. East-west morphological asymmetry of Barbuda and difference in LIG shoreline elevation between Barbuda and Antigua suggest a regional tectonic process. As with the proposed westward tilting from the forearc to the volcanic arc of the Guadeloupe archipelago, vertical deformation on Barbuda could be related to plate-scale subduction processes. Long-term uplift of Barbuda might be related to the accumulation of residual coseismic deformation not fully recovered by interseismic subsidence and the corresponding seismogenic segment would extend below the Moho.

Geochemistry, Geophysics, Geosystems↗
USGS2021

Earthquake magnitude distributions on northern Caribbean faults from combinatorial optimization models

On-fault earthquake magnitude distributions are calculated for northern Caribbean faults using estimates of fault slip and regional seismicity parameters. Integer programming, a combinatorial optimization method, is used to determine the optimal spatial arrangement of earthquakes sampled from a truncated Gutenberg-Richter distribution that minimizes the global misfit in slip rates on a complex fault system. Slip rates and their uncertainty on major faults are derived from a previously published GPS block model for the region, with fault traces determined from offshore geophysical mapping and previously published onshore studies. The optimal spatial arrangement of the sampled earthquakes is compared with the 500-year history of earthquake observations. Rupture segmentation of the subduction interface along the Hispaniola-Puerto Rico Trench (PRT) fault and seismic coupling on the PRT fault appear to exert the primary control over this spatial arrangement. Introducing a rupture barrier for the Hispaniola-PRT fault northwest of Mona Passage, based on geophysical and seismicity observations, and assigning a low slip rate of 2 mm/yr on the PRT fault are most consistent with historical earthquakes in the region. The addition of low slip-rate secondary faults as well as segmentation of the Hispaniola and Septentrional strike-slip fault improves the consistency with historical seismicity. An important observation from the modeling is that varying the slip rate on the PRT fault and different segmentation scenarios result in significant changes to the optimal magnitude distribution on faults farther away. In general, optimal on-fault magnitude distributions are more complex and inter-dependent than is typically assumed in probabilistic seismic hazard analysis and probabilistic tsunami hazard analysis.

JGR Solid Earth↗
USGS2021

Hydrate formation on marine seep bubbles and the implications for water column methane dissolution

Methane released from seafloor seeps contributes to a number of benthic, water column, and atmospheric processes. At seafloor seeps within the methane hydrate stability zone, crystalline gas hydrate shells can form on methane bubbles while the bubbles are still in contact with the seafloor or as the bubbles begin ascending through the water column. These shells reduce methane dissolution rates, allowing hydrate-coated bubbles to deliver methane to shallower depths in the water column than hydrate-free bubbles. Here, we analyze seafloor videos from six deepwater seep sites associated with a diverse range of bubble-release processes involving hydrate formation. Bubbles that grow rapidly are often hydrate-free when released from the seafloor. As bubble growth slows and seafloor residence time increases, a hydrate coating can form on the bubble's gas-water interface, fully coating most bubbles within ∼10 s of the onset of hydrate formation at the seafloor. This finding agrees with water-column observations that most bubbles become hydrate-coated after their initial ∼150 cm of rise, which takes about 10 s. Whether a bubble is coated or not at the seafloor affects how much methane a bubble contains and how quickly that methane dissolves during the bubble's rise through the water column. A simplified model shows that, after rising 150 cm above the seafloor, a bubble that grew a hydrate shell before releasing from the seafloor will have ∼5% more methane than a bubble of initial equal volume that did not grow a hydrate shell after it traveled to the same height.

Journal of Geophysical Research - Oceans↗
USGS2022

Application of recursive estimation to heat tracing for groundwater/surface-water exchange

We present and demonstrate a recursive-estimation framework to infer groundwater/surface-water exchange based on temperature time series collected at different vertical depths below the sediment/water interface. We formulate the heat-transport problem as a state-space model (SSM), in which the spatial derivatives in the convection/conduction equation are approximated using finite differences. The SSM is calibrated to estimate time-varying specific discharge using the Extended Kalman Filter (EKF) and Extended Rauch-Tung-Striebel Smoother (ERTSS). Whereas the EKF is suited to real-time (“online”) applications and uses only the past and current measurements for estimation (filtering), the ERTSS is intended for near-real time or batch-processing (“offline”) applications and uses a window of data for batch estimation (smoothing). The two algorithms are demonstrated with synthetic and field-experimental data and are shown to be efficient and rapid for the estimation of time-varying flux over seasonal periods; further, the recursive approaches are effective in the presence of rapidly changing flux and (or) nonperiodic thermal boundary conditions, both of which are problematic for existing approaches to heat tracing of time-varying groundwater/surface-water exchange.

Water Resources Research↗
USGS2022

Using ensemble data assimilation to estimate transient hydrologic exchange flow under highly dynamic flow conditions

Quantifying dynamic hydrologic exchange flows (HEFs) within river corridors that experience high-frequency flow variations caused by dam regulations is important for understanding the biogeochemical processes at the river water and groundwater interfaces. Heat has been widely used as a tracer to infer steady-state flow velocities through analytical solutions of heat transport defined by the diurnal temperature signals. Under sub-daily dynamic flow conditions, however, such analytical solutions are not applicable due to the violation of their fundamental assumptions. In this study, we developed a data assimilation-based approach to estimate the sub-daily flux under highly dynamic flow conditions using multi-depth temperature observations at a 5-min resolution. If the hydraulic gradient is measured, Darcy's law was used to calculate the flux with permeability estimated from temperature responses below the riverbed. Otherwise, flux was estimated directly by assimilating multi-depth temperature data at 1- or 2-hr time intervals assuming one-dimensional flow and heat transport governing equation. By comparing estimated fluxes with model-generated synthetic truth, we demonstrated that both schemes have robust performance in estimating fluxes under highly dynamic flow conditions. This data assimilation-based flux estimation method was able to capture the vertical sub-daily fluxes using multi-depth high-resolution temperature data alone, even in the presence of multi-dimensional flow. This approach has been successfully applied to real field temperature data collected at the Hanford site, which experiences highly dynamic HEFs. Our study shows the promise of adopting distributed 1-D temperature monitoring to capture spatial and temporal exchange dynamics in river corridors at a watershed scale or beyond.

Washington↗
USGS2023

Long-term changes in concentrations and yield of riverine dissolved silicon from the poles to the tropics

Riverine exports of silicon (Si) influence global carbon cycling through the growth of marine diatoms, which account for ∼25% of global primary production. Climate change will likely alter river Si exports in biome-specific ways due to interacting shifts in chemical weathering rates, hydrologic connectivity, and metabolic processes in aquatic and terrestrial systems. Nonetheless, factors driving long-term changes in Si exports remain unexplored at local, regional, and global scales. We evaluated how concentrations and yields of dissolved Si (DSi) changed over the last several decades of rapid climate warming using long-term data sets from 60 rivers and streams spanning the globe (e.g., Antarctic, tropical, temperate, boreal, alpine, Arctic systems). We show that widespread changes in river DSi concentration and yield have occurred, with the most substantial shifts occurring in alpine and polar regions. The magnitude and direction of trends varied within and among biomes, were most strongly associated with differences in land cover, and were often independent of changes in river discharge. These findings indicate that there are likely diverse mechanisms driving change in river Si biogeochemistry that span the land-water interface, which may include glacial melt, changes in terrestrial vegetation, and river productivity. Finally, trends were often stronger in months outside of the growing season, particularly in temperate and boreal systems, demonstrating a potentially important role of shifting seasonality for the flux of Si from rivers. Our results have implications for the timing and magnitude of silica processing in rivers and its delivery to global oceans.

Global Biogeochemical Cycles↗
USGS2023

Substantial upper plate faulting above a shallow subduction megathrust earthquake: Mechanics and implications of the surface faulting during the 2016 Kaikoura, New Zealand, earthquake

The 2016 moment magnitude 7.8 Kaikoura, New Zealand, earthquake occurred at the southern end of the Hikurangi subduction zone where the upper plate above the shallow megathrust is exposed sub-aerially. As a result, the substantial co-seismic deformation in the upper plate above the megathrust rupture was observed geologically and geodetically. We explore the relationship between this surface faulting and the subduction megathrust rupture and find that the greatest upper plate fault slip occurred coincident (in time and location) with the megathrust rupture. Models of Coulomb stress change demonstrate that these surface faults become positively loaded as the upper plate rebounds during the megathrust event, favoring fault slip. In addition, during the megathrust rupture these faults terminate against an uncoupled subduction plate interface. We simulate the effects of decoupling at the base of these faults and find that very large fault slip is an expected consequence of this decoupling, allowing near-complete strain release. In contrast, typical strike-slip faults, pinned at their base, would have lower amounts of fault slip. These two conditions—increased Coulomb stress and basal decoupling—combine to produce the extreme co-seismic upper plate faulting observed above the shallow Kaikoura megathrust earthquake. Similar conditions occur in other global subduction zones, but in most subduction zones the region above the coupled megathrust is underwater and poorly observed. Our analysis of the Kaikoura earthquake indicates a need to reevaluate patterns of strain accumulation and release in these regions, rather than assuming simple models of elastic rebound.

Tectonics↗
USGS2024

Deep structure of Siletzia in the Puget Lowland: Imaging an obducted plateau and accretionary thrust belt with potential fields

Detailed understanding of crustal components and tectonic history of forearcs is important due to their geological complexity and high seismic hazard. The principal component of the Cascadia forearc is Siletzia, a composite basaltic terrane of oceanic origin. Much is known about the lithology and age of the province. However, glacial sediments blanketing the Puget Lowland obscure its lateral extent and internal structure, hindering our ability to fully understand its tectonic history and its influence on modern deformation. In this study, we apply map-view interpretation and two-dimensional modeling of aeromagnetic and gravity data to the magnetically stratified Siletzia terrane revealing its internal structure and characterizing its eastern boundary. These analyses suggest the contact between Siletzia (Crescent Formation) and the Eocene accretionary prism trends northward under Lake Washington. North of Seattle, this boundary dips east where it crosses the Kingston arch, whereas south of Seattle the contact dips west where it crosses the Seattle uplift (SU). This westward dip is opposite the dip of the Eocene subduction interface, implying obduction of Siletzia upper crust at this southern location. Elongate pairs of high and low magnetic anomalies over the SU suggest imbrication of steeply-dipping, deeply rooted slices of Crescent Formation within Siletzia. We hypothesize these features result from duplication of Crescent Formation in an accretionary fold-thrust belt during the Eocene. The active Seattle fault divides this Eocene fold-thrust belt into two zones with different structural trends and opposite frontal ramp dips, suggesting the Seattle fault may have originated as a tear fault during accretion.

Washington↗
USGS2024

Subsurface redox interactions regulate ebullitive methane flux in heterogeneous Mississippi River deltaic wetland

As interfaces connecting terrestrial and ocean ecosystems, coastal wetlands develop temporally and spatially complex redox conditions, which drive uncertainties in greenhouse gas emission as well as the total carbon budget of the coastal ecosystem. To evaluate the role of complex redox reactions in methane emission from coastal wetlands, a coupled reactive-transport model was configured to represent subsurface biogeochemical cycles of carbon, nitrogen, and sulfur, along with production and transport of multiple gas species through diffusion and ebullition. This model study was conducted at multiple sites along a salinity gradient in the Barataria Basin at the Mississippi River Deltaic Plain. Over a freshwater to saline gradient, simulated total flux of methane was primarily controlled by its subsurface production and consumption, which were determined by redox reactions directly (e.g., methanogenesis, methanotrophy) and indirectly (e.g., competition with sulfate reduction) under aerobic and/or anaerobic conditions. At fine spatiotemporal scales, surface methane fluxes were also strongly dependent on transport processes, with episodic ebullitive fluxes leading to higher spatial and temporal variability compared to the gradient-driven diffusion flux. Ebullitive methane fluxes were determined by methane fraction in total ebullitive gas and the frequency of ebullitive events, both of which varied with subsurface methane concentrations and other gas species. Although ebullition thresholds are constrained by local physical factors, this study indicates that redox interactions not only determine gas composition in ebullitive fluxes but can also regulate ebullition frequency through gas production.

Journal of Advances in Modeling Earth Systems↗
USGS2024

Using stochastic point pattern analysis to track regional orientations of magmatism during the transition to cenozoic extension and Rio Grande rifting, Southern Rocky Mountains

The southern Rocky Mountains in Colorado and northern New Mexico hosted intracontinental magmatism that developed during a tectonic transition from shortening (Laramide orogeny, ca. 75 to 40 Ma) through extension and rifting. We present a novel approach that uses stochastic weighted bootstrap simulations of a large set of new and historical geochronology data to better understand how regional anisotropies responsible for focusing magma emplacement evolved through time. This technique can detect subtle trends in directional distributions, including multi-modal orientations, and can be filtered from regional to local scales. Our results indicate that magmatism followed first the northeast trend of the Colorado mineral belt between 75 and 40 Ma and deviated afterward. These deviations vary depending on the scale of the analysis. At the smallest scale we evaluated (<75 km), the orientation of magmatism from 45 to 30 Ma rotated counter-clockwise before aligning with the north-south trend of the modern Rio Grande rift. Larger, regional-scale analyses indicate magma centers between 40 to 35 Ma and 25 to 20 Ma were dominantly oriented southwest-northeast, whereas magmatism between 35 and 25 Ma had north-south orientation. The large areal footprint of magmatism and shifting regional patterns suggest that ancient zones of weakness in the North American lithosphere accommodated magma flow at different moments in time, rather than controlled by a retreating interface of the Farallon and North American plates.

Tectonics↗