USGS Science⌕ Search

SEARCH · USGS Science

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 955 records · Page 53Linked to original sources
USGS2012

A meeting of the waters: interdisciplinary challenges and opportunities in tidal rivers

At the interface of estuarine tides and freshwater rivers lie wetland and aquatic ecosystems, which experience dramatic effects of sea level rise. There, nontidal channels and riparian floodplains are transforming into tidal ecosystems, and tidal freshwater ecosystems are receiving increasing salinity. These river-floodplain systems have both fluvial characteristics, including meandering channels and expansive floodplain forests, and estuarine characteristics, including tides and intertidal wetlands [see Barendregt et al., 2009; Conner et al., 2007, and references therein]. Because tidal rivers lie at the disciplinary divide between fluvial and estuarine science, a knowledge gap has developed in scientists' understanding of the geomorphic and biogeochemical response of these environments to sea level rise, climate change, and anthropogenically driven variations in watershed exports.

Eos, Transactions, American Geophysical Union↗
USGS2012

Carbon dioxide and methane emissions from the Yukon River system

Carbon dioxide (CO 2 ) and methane (CH 4 ) emissions are important, but poorly quantified, components of riverine carbon (C) budgets. This is largely because the data needed for gas flux calculations are sparse and are spatially and temporally variable. Additionally, the importance of C gas emissions relative to lateral C exports is not well known because gaseous and aqueous fluxes are not commonly measured on the same rivers. We couple measurements of aqueous CO 2 and CH 4 partial pressures ( p CO 2 , p CH 4 ) and flux across the water-air interface with gas transfer models to calculate subbasin distributions of gas flux density. We then combine those flux densities with remote and direct observations of stream and river water surface area and ice duration, to calculate C gas emissions from flowing waters throughout the Yukon River basin. CO 2 emissions were 7.68 Tg C yr −1 (95% CI: 5.84 −10.46), averaging 750 g C m −2 yr −1 normalized to water surface area, and 9.0 g C m −2 yr −1 normalized to river basin area. River CH 4 emissions totaled 55 Gg C yr −1 or 0.7% of the total mass of C emitted as CO 2 plus CH 4 and ∼6.4% of their combined radiative forcing. When combined with lateral inorganic plus organic C exports to below head of tide, C gas emissions comprised 50% of total C exported by the Yukon River and its tributaries. River CO 2 and CH 4 derive from multiple sources, including groundwater, surface water runoff, carbonate equilibrium reactions, and benthic and water column microbial processing of organic C. The exact role of each of these processes is not yet quantified in the overall river C budget.

Yukon River Basin↗
USGS2012

Plate interaction in the NE Caribbean subduction zone from continuous GPS observations

Kinematic similarities between the Sumatra and Puerto Rico Trenches highlight the potential for a mega-earthquake along the Puerto Rico Trench and the generation of local and trans-Atlantic tsunamis. We used the horizontal components of continuous GPS (cGPS) measurements from 10 sites on NE Caribbean islands to evaluate strain accumulation along the North American (NA) - Caribbean (CA) plate boundary. These sites move westward and slightly northward relative to CA interior at rates ≤2.5 mm/y. Provided this motion originates in the subduction interface, the northward motion suggests little or no trench-perpendicular thrust accumulation and may in fact indicate divergence north of Puerto Rico, where abnormal subsidence, bathymetry, and gravity are observed. The Puerto Rico Trench, thus, appears unable to generate mega-earthquakes, but damaging smaller earthquakes cannot be discounted. The westward motion, characterized by decreasing rate with distance from the trench, is probably due to eastward motion of CA plate impeded at the plate boundary by the Bahamas platform. Two additional cGPS sites in Mona Passage and SW Puerto Rico move to the SW similar to Hispaniola and unlike the other 10 sites. That motion relative to the rest of Puerto Rico may have given rise to seismicity and normal faults in Mona Rift, Mona Passage, and SW Puerto Rico.

Geophysical Research Letters↗
USGS2012

Juan de Fuca slab geometry and its relation to Wadati-Benioff zone seismicity

A new model of the subducted Juan de Fuca plate beneath western North America allows first-order correlations between the occurrence of Wadati-Benioff zone earthquakes and slab geometry, temperature, and hydration state. The geo-referenced 3D model, constructed from weighted control points, integrates depth information from earthquake locations and regional seismic velocity studies. We use the model to separate earthquakes that occur in the Cascadia forearc from those that occur within the underlying Juan de Fuca plate and thereby reveal previously obscured details regarding the spatial distribution of earthquakes. Seismicity within the slab is most prevalent where the slab is warped beneath northwestern California and western Washington suggesting that slab flexure, in addition to expected metamorphic dehydration processes, promotes earthquake occurrence within the subducted oceanic plate. Earthquake patterns beneath western Vancouver Island are consistent with slab dehydration processes. Conversely, the lack of slab earthquakes beneath western Oregon is consistent with an anhydrous slab. Double-differenced relocated seismicity resolves a double seismic zone within the slab beneath northwestern California that strongly constrains the location of the plate interface and delineates a cluster of seismicity 10 km above the surface that includes the 1992 M7.1 Mendocino earthquake. We infer that this earthquake ruptured a surface within the Cascadia accretionary margin above the Juan de Fuca plate. We further speculate that this earthquake is associated with a detached fragment of former Farallon plate. Other subsurface tectonic elements within the forearc may have the potential to generate similar damaging earthquakes.

British Columbia;California;Oregon;Washington↗
USGS2012

Drivers of change in estuarine-coastal ecosystems: Discoveries from four decades of study in San Francisco Bay

Poised at the interface of rivers, ocean, atmosphere and dense human settlement, estuaries are driven by a large array of natural and anthropogenic forces. San Francisco Bay exemplifies the fast-paced change occurring in many of the world's estuaries, bays and inland seas in response to these diverse forces. We use observations from this particularly well-studied estuary to illustrate responses to six drivers that are common agents of change where land and sea meet: water consumption and diversion; human modification of sediment supply; introduction of non-native species; sewage input; environmental policy; and climate shifts. In San Francisco Bay, responses to these drivers include, respectively, shifts in the timing and extent of freshwater inflow and salinity intrusion; decreasing turbidity; restructuring of plankton communities; nutrient enrichment; elimination of hypoxia and reduced metal contamination of biota; and food web changes that decrease resistance of the estuary to nutrient pollution. Detection of these changes and discovery of their causes through environmental monitoring have been essential for establishing and measuring outcomes of environmental policies that aim to maintain high water quality and sustain services provided by estuarine-coastal ecosystems. The wide range of variability time scales and the multiplicity of interacting drivers place heavy demands on estuarine monitoring programs. But the San Francisco Bay case study illustrates why the imperative for monitoring has never been greater.

California↗
USGS2019

The 12 November 2017 Mw 7.3 Ezgeleh–Sarpolzahab (Iran) earthquake and active tectonics of the Lurestan arc

The 12 November 2017 M w 7.3 Ezgeleh‐Sarpolzahab earthquake is the largest instrumentally recorded earthquake in the Zagros Simply Folded Belt by a factor of ∼10 in seismic moment. Exploiting local, regional, and teleseismic data and synthetic aperture radar interferometry imagery, we characterize the rupture, its aftershock sequence, background seismicity, and regional tectonics. The mainshock ruptured slowly (∼2 km/s), unilaterally southward, for ∼40 km along an oblique (dextral‐thrust) fault that dips ∼14°E beneath the northwestern Lurestan arc. Slip is confined to basement depths of ∼12–18 km, resolvably beneath the sedimentary cover which is ∼8 km thick in this area. The gentle dip angle and basement location allow for a broad slip area, explaining the large magnitude relative to earthquakes in the main Fars arc of the Zagros, where shallower, steeper faults are limited in rupture extent by weak sedimentary layers. Early aftershocks concentrate around the southern and western edges of the mainshock slip area and therefore cluster in the direction of rupture propagation, implying a contribution from dynamic triggering. A cluster of events ∼100 km to the south near Mandali (Iraq) reactivated the ∼50° dipping Zagros Foredeep Fault. The basement fault responsible for the Ezgeleh‐Sarpolzahab earthquake probably accounts for the ∼1 km elevation contrast between the Lurestan arc and the Kirkuk embayment but is distinct from sections of the Mountain Front Fault that define frontal escarpments elsewhere in the Zagros. It may be related to a seismic interface underlying the central and southern Lurestan arc, and a key concern is whether or not the more extensive regional structure is also seismogenic.

Journal of Geophysical Research: Solid Earth↗
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↗
USGS2019

Seasonal and spatial variation in the location and reactivity of a nitrate-contaminated groundwater discharge zone in a lakebed

Groundwater discharge delivering anthropogenic N from surrounding watersheds can impact lake nutrient budgets. However, upgradient groundwater processes and changing dynamics in N biogeochemistry at the groundwater-lake interface are complex and difficult to resolve. In this study, hydrograph variations in a groundwater flow-through lake altered discharge patterns of a wastewater-derived, groundwater contaminant plume, thereby affecting biogeochemical processes controlling N transport. Groundwater geochemistry 15 cm under the lakebed along transects perpendicular to shore varied from oxic to anoxic with increasing nitrate concentrations (10-75 M) and corresponding gradients in nitrite and nitrous oxide. Porewater depth profiles of nitrate concentrations and stable isotope compositions largely reflected upgradient groundwater N sources and N-cycle processes, with minor additional nitrate reduction in the shallowest lakebed sediments. Potential denitrification rates determined in laboratory microcosms were 10-100 fold higher in near-surface sediments (0-5 cm) than in deeper sediments (5-30 cm) and were correlated with sediment carbon content and abundance of denitrification genes (nirS, nosZI, and nosZII). Potential anammox-driven N2 production was highest in deeper anoxic sediments. Injection of bromide and nitrite in the lake sediments indicated a vertical porewater velocity of 4-5 cm hr-1, with highest nitrite consumption rates above 10 cm. However, short residence times in the shallow sediments allowed only a small fraction of the contaminant nitrate to be removed before discharging into the lake. Results demonstrate the importance of resolving local versus upgradient biogeochemical processes affecting contaminant distribution in discharge areas, and transient migration of local gradients and processes in response to changing lake levels and groundwater flow paths.

Journal of Geophysical Research: Biogeosciences↗
USGS2019

Residence time controls on the fate of nitrogen in flow‐through lakebed sediments

For many glacial lakes with highly permeable sediments, water exchange rates control hydrologic residence times within the sediment‐water interface (SWI) and the removal of reactive compounds such as nitrate, a common pollutant in lakes and groundwater. Here we conducted a series of focused tracer injection experiments in the upper 20 cm of the naturally downwelling SWI in a flow‐through lake on Cape Cod, MA. We systematically varied residence time and reactant controls on nitrate processing, using isotopically labeled 15 N nitrate to monitor the effect of these changes on nitrate removal via denitrification. The addition of acetate, a labile carbon compound, triggered the lake SWI to switch from net production to net removal of nitrate. When acetate was combined with increased residence time created by controlled reductions in water flux, we observed a fivefold increase in nitrate removal, a 26‐fold increase in N 2 production, and a 42‐fold increase in N 2 O production. We demonstrate that water residence time is an important control on the fate of nitrate in these lake SWIs and illustrate that seasonal conditions that alter lake exchange rates and variability in lake carbon may predict dynamic nitrate removal across the SWI. Additionally, observed N 2 O production during the oxic pore water experiments paired with geophysical characterization of the sediment porosity revealed that the lake SWI has less mobile pores occupying upward of 50% of the total porosity volume, which function as reactive microzones for nitrate processing.

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