SEARCH · Search USGS
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.
Simulated results of the coupled freshwater-saltwater sharp interface and convective-dispersive numerical models are compared by using steady-state cross-sectional simulations. The results indicate that in some aquifers the calculated sharp interface is located further landward than would be expected.
Radon 222 concentrations in the water and sedimentary columns and radon exchange rates across the sediment-water and air-water interfaces have been measured in a section of south San Francisco Bay. Two independent methods have been used to determine sediment-water exchange rates, and the annual averages of these methods agree within the uncertainty of the determinations, about 20%. The annual average of benthic fluxes from shoal areas is nearly a factor of 2 greater than fluxes from the channel areas. Fluxes from the shoal and channel areas exceed those expected from simple molecular diffusion by factors of 4 and 2, respectively, apparently due to macrofaunal irrigation. Values of the gas transfer coefficient for radon exchange across the air-water interface were determined by constructing a radon mass balance for the water column and by direct measurement using floating chambers. The chamber method appears to yield results which are too high. Transfer coefficients computed using the mass balance method range from 0.4 m/day to 1.8 m/day, with a 6-year average of 1.0 m/day. Gas exchange is linearly dependent upon wind speed over a wind speed range of 3.2–6.4 m/s, but shows no dependence upon current velocity. Gas transfer coefficients predicted from an empirical relationship between gas exchange rates and wind speed observed in lakes and the oceans are within 30% of the coefficients determined from the radon mass balance and are considerably more accurate than coefficients predicted from theoretical gas exchange models.
Along coastal interfaces, components of the Earth system interact to regulate ecosystem functions and Earth’s climate. Between the land and ocean, diverse coastal ecosystem types transform, store, and transport material. A dynamic two-way exchange of energy and matter is driven by hydrological and marine processes such as river and groundwater discharge, tides, waves, and storms. Global models lack representation of coastal processes and related feedbacks on Earth’s climate and ecosystems, impeding their predictions of coastal and global response to change. We recommend leveraging existing monitoring networks and regional models to develop and validate global models that include the coastal interface.
The interface between ground water and surface water in streams is a hotspot for N processing. However, the role of buried organic C in N transformation at this interface is not well understood, and inferences have been based largely on descriptive studies. Our main objective was to determine how buried particulate organic C (POC) affected denitrification and NO 3 − retention in the sediments of an upwelling reach in a sand-plains stream in Wisconsin. We manipulated POC in mesocosms inserted in the sediments. Treatments included low and high quantities of conditioned red maple leaves (buried beneath combusted sand), ambient sediment (sand containing background levels of POC), and a control (combusted sand). We measured denitrification rates in sediments by acetylene-block assays in the laboratory and by changes in N 2 concentrations in the field using membrane inlet mass spectrometry. We measured NO 3 −, NH 4 +, and dissolved organic N (DON) retention as changes in concentrations and fluxes along groundwater flow paths in the mesocosms. POC addition drove oxic ground water to severe hypoxia, led to large increases in dissolved organic C (DOC), and strongly increased denitrification rates and N (NO 3 − and total dissolved N) retention relative to the control. In situ denitrification accounted for 30 to 60% of NO 3 − retention. Our results suggest that buried POC stimulated denitrification and NO 3 − retention by producing DOC and by creating favorable redox conditions for denitrification.
Managing pathogen spillover at the wildlife–livestock interface is a key step towards improving global animal health, food security and wildlife conservation. However, predicting the effectiveness of management actions across host–pathogen systems with different life histories is an on-going challenge since data on intervention effectiveness are expensive to collect and results are system-specific. We developed a simulation model to explore how the efficacies of different management strategies vary according to host movement patterns and epidemic growth rates. The model suggested that fast-growing, fast-moving epidemics like avian influenza were best-managed with actions like biosecurity or containment, which limited and localized overall spillover risk. For fast-growing, slower-moving diseases like foot-and-mouth disease, depopulation or prophylactic vaccination were competitive management options. Many actions performed competitively when epidemics grew slowly and host movements were limited, and how management efficacy related to epidemic growth rate or host movement propensity depended on what objective was used to evaluate management performance. This framework offers one means of classifying and prioritizing responses to novel pathogen spillover threats, and evaluating current management actions for pathogens emerging at the wildlife–livestock interface.
The characterization of pore-space connectivity in porous media at the sediment/water interface is critical to understanding contaminant transport and reactive biogeochemical processes in zones of groundwater and surface-water exchange. Previous in situ studies of dual-domain (i.e., mobile/less-mobile porosity) studies have been limited to solute tracer injections at scales of meters to 100s of meters and subsequent numerical model parameterization using fluid concentration histories. Pairing fine-scale (e.g., sub-meter) geoelectrical measurements with fluid tracer data over time alleviates dependence on flowpath-scale experiments, enabling spatially targeted characterization of shallow sediment/water interface media where biogeochemical reactivity is often high. The Dual-Domain Porosity Apparatus is a field-tested device capable of variable rate-controlled downward flow experiments. The Dual-Domain Porosity Apparatus facilitates meter-scale inference of dual-domain parameters, i.e., mobile/less-mobile exchange rate coefficient and the ratio of less mobile to mobile porosity. The Dual-Domain Porosity Apparatus experimental procedure uses water electrical conductivity as a conservative tracer of differential loading and flushing of pore spaces within the region of measurement. Variable injection rates permit the direct quantification of the flow-dependence of dual-domain parameters, which has been theorized for decades but remains challenging to assess using existing experimental methodologies.
A key step in many earthquake source inversions requires knowledge of the geometry of the fault surface on which the earthquake occurred. Our knowledge of this surface is often uncertain, however, and as a result fault geometry misinterpretation can map into significant error in the final temporal and spatial slip patterns of these inversions. Relying solely on an initial hypocentre and CMT mechanism can be problematic when establishing rupture characteristics needed for rapid tsunami and ground shaking estimates. Here, we attempt to improve the quality of fast finite-fault inversion results by combining several independent and complementary data sets to more accurately constrain the geometry of the seismic rupture plane of subducting slabs. Unlike previous analyses aimed at defining the general form of the plate interface, we require mechanisms and locations of the seismicity considered in our inversions to be consistent with their occurrence on the plate interface, by limiting events to those with well-constrained depths and with CMT solutions indicative of shallow-dip thrust faulting. We construct probability density functions about each location based on formal assumptions of their depth uncertainty and use these constraints to solve for the ‘most-likely’ fault plane. Examples are shown for the trench in the source region of the M w 8.6 Southern Sumatra earthquake of March 2005, and for the Northern Chile Trench in the source region of the November 2007 Antofagasta earthquake. We also show examples using only the historic catalogues in regions without recent great earthquakes, such as the Japan and Kamchatka Trenches. In most cases, this method produces a fault plane that is more consistent with all of the data available than is the plane implied by the initial hypocentre and CMT mechanism. Using the aggregated data sets, we have developed an algorithm to rapidly determine more accurate initial fault plane geometries for source inversions of future earthquakes.
The Mendocino triple junction, the intersection of the Pacific, North American, and Gorda plates, activates a collection of disparate faults that reconcile Cascadia subduction with San Andreas transform motion. The December 20, 2022, Mw 6.4 Ferndale, California earthquake occurred within this complex zone as strike-slip faulting within the subducting Gorda slab. Here, we analyze the seismic and geodetic signatures of the mainshock and aftershock sequence to illuminate its role within complex tectonic surroundings. We find aftershocks on varied fault structures within the uppermost Gorda slab, yet seismicity on the subduction interface itself was notably absent. Nevertheless, we identify small but coherent postseismic deformation that is well modeled by aseismic slip on this interface, likely triggered by stresses generated at the updip limit of coseismic rupture. This sequence demonstrates the potential for interactions between intra-slab earthquakes and slip on the subduction megathrust, highlighting the need to consider this and other subduction zones as coupled systems of interacting faults.
While the recent incursion of highly pathogenic avian influenza into North America has resulted in notable losses to the commercial poultry industry, the mechanism by which virus enters commercial poultry houses is still not understood. One theorized mechanism is that waterfowl shed virus into the environment surrounding poultry farms, such as into retention ponds, and is then transmitted into poultry houses via bridge species. Little is known about if and when wild waterfowl use these retention ponds, leading to uncertainty regarding the potential significance of this interface. To quantify the use of retention ponds on commercial poultry farms by wild waterfowl, we surveyed 12 such ponds across Somerset and Dorchester counties, Maryland, USA. This region was chosen due to the high level of poultry production and its importance for migratory waterfowl. Surveys consisted of recording waterfowl visible on the retention ponds from public roadways at least once per week from 20 September 2022–31 March 2023. Throughout the course of this study, we observed a total of nine species of waterfowl using retention ponds on commercial poultry farms at nine of 12 sites. The number of waterfowl observed at retention ponds varied notably throughout the course of our survey period, with values generally following trends of fall migration within each species indicating that resident birds were not the only individuals to utilize these habitats. Additionally, waterfowl use was highest at sites with little vegetation immediately surrounding the pond, and lowest when ponds were surrounded by trees. Our data suggest that retention ponds on commercial poultry farms present a notable interface for waterfowl to introduce avian influenza viruses to farm sites. However, additional testing and surveys could provide further insight into whether it may be possible to reduce the use of these habitats by wild waterfowl through vegetative management as preliminarily reported here.
Background Characterizing the interface between wild and domestic animal populations is increasingly recognized as essential in the context of emerging infectious diseases (EIDs) that are transmitted by wildlife. More specifically, the spatial and temporal distribution of contact rates between wild and domestic hosts is a key parameter for modeling EIDs transmission dynamics. We integrated satellite telemetry, remote sensing and ground-based surveys to evaluate the spatio-temporal dynamics of indirect contacts between wild and domestic birds to estimate the risk that avian pathogens such as avian influenza and Newcastle viruses will be transmitted between wildlife to poultry. We monitored comb ducks ( Sarkidiornis melanotos melanotos ) with satellite transmitters for seven months in an extensive Afro-tropical wetland (the Inner Niger Delta) in Mali and characterise the spatial distribution of backyard poultry in villages. We modelled the spatial distribution of wild ducks using 250-meter spatial resolution and 8-days temporal resolution remotely-sensed environmental indicators based on a Maxent niche modelling method. Results Our results show a strong seasonal variation in potential contact rate between wild ducks and poultry. We found that the exposure of poultry to wild birds was greatest at the end of the dry season and the beginning of the rainy season, when comb ducks disperse from natural water bodies to irrigated areas near villages. Conclusions Our study provides at a local scale a quantitative evidence of the seasonal variability of contact rate between wild and domestic bird populations. It illustrates a GIS-based methodology for estimating epidemiological contact rates at the wildlife and livestock interface integrating high-resolution satellite telemetry and remote sensing data.
Complex seismic velocity structure near the earthquake source can affect rupture dynamics and strongly modify the seismic waveforms recorded near the fault. Fault‐zone waves are commonly observed in continental crustal settings but are less clear in subduction zones due to the spatial separation between seismic stations and the plate boundary fault. We observed anomalously long duration S waves from earthquake clusters located near the interface of the subducted Gorda plate north of the Mendocino triple junction. In contrast, earthquakes located just a few kilometers below each cluster show impulsive S waves. A nodal array experiment was conducted around the Northern California Seismic Network station KCT for two months to investigate the origin of the complex S waves. Beamforming analysis shows that the S waves contain three arrivals that have different horizontal slownesses, which we term S 1, S 2, and S coda. Similar analysis on P waves also show two arrivals with different horizontal slownesses, which we term P 1 and P 2. P 1 and S 1 have larger horizontal slowness than P 2 and S 2, respectively, indicating that the phase pairs are body waves with different ray paths. Building upon a seismic refraction profile, we construct 1D velocity models and test different thicknesses and V P / V S ratios for the subducted oceanic crust. The arrival times and relative slownesses of P 1/ P 2 and S 1/ S 2 phases indicate that they are the direct and the Moho reflected phases, respectively. Their properties are consistent with a crustal thickness of ∼6 km and a moderate V P / V S ratio (∼1.8). The S coda is more difficult to characterize but has a clear dominant frequency that likely reflects the near‐source velocity and attenuation structure. Our study indicates that waveforms from earthquakes near the interface of the subducted slab can be used to infer detailed structural information about the plate‐boundary zone at seismogenic depths.
A continually increasing number of high‐quality digital strong‐motion records from stations of the National Strong Motion Project (NSMP) of the U.S. Geological Survey, as well as data from regional seismic networks within the United States, calls for automated processing of strong‐motion records with human review limited to selected significant or flagged records. The NSMP has developed the Processing and Review Interface for Strong Motion data (PRISM) software to meet this need. In combination with the Advanced National Seismic System Quake Monitoring System (AQMS), PRISM automates the processing of strong‐motion records. When used without AQMS, PRISM provides batch‐processing capabilities. The PRISM software is platform independent (coded in Java), open source, and does not depend on any closed‐source or proprietary software. The software consists of two major components: a record processing engine composed of modules for each processing step, and a review tool, which is a graphical user interface for manual review, edit, and processing. To facilitate use by non‐NSMP earthquake engineers and scientists, PRISM (both its processing engine and review tool) is easy to install and run as a stand‐alone system on common operating systems such as Linux, OS X, and Windows. PRISM was designed to be flexible and extensible to accommodate implementation of new processing techniques. All the computing features have been thoroughly tested.
Source–station distance is a central input to ground‐motion models (GMMs) for predicting seismic shaking. GMM development uses distance metrics including the Joyner–Boore distance, which is the shortest distance from an observation point to the surface projection of the earthquake rupture, and R rup the shortest distance to the rupture in three dimensions. Thompson and Baltay (2018) proposed the generalized mean rupture distance R p to address observed near‐fault ground‐motion saturation. R p accounts for the contribution to the shaking of all parts of the rupture and provides a simple method for incorporating spatially variable slip. They used R p to develop a GMM for shallow crustal earthquakes, assuming uniform slip. Here, we investigate the improvement offered by an R p ‐based GMM for large subduction interface earthquakes by recalibrating the path term for a published GMM ( Parker et al. , 2022 ) using R p derived from distributed slip models (DSMs). Inspection of within‐event and total residuals indicates that the recalibrated model fits the data at least as well as an alternative in which the path term was recalibrated using R rup and the same dataset. Incorporating slip information in its entirety or trimming the DSM geometry is preferable to assuming uniform moment release on a prescribed model fault that extends beyond the actual ruptured area. R p tuned to minimize model uncertainty is closer to the maximum distance to rupture for peak ground velocity (PGV) than acceleration, possibly reflecting the contribution to high‐frequency shaking of slip on local asperities. However, when the moment is concentrated far from stations, PGV is fit adequately with R p closer to R rup . Our results suggest incorporating slip‐derived moment release through R p could improve GMMs for slab interface events, especially if the implementation of R p ‐based models is refined using a larger dataset of earthquakes with greater geographic diversity to account for regional ground‐motion variations.
The layer-structure minerals serpentine, brucite, and talc are postulated to form in the mantle wedge above a subducting slab as a result of progressive hydration and silica metasomatism. Tectonic mixing at the slab-mantle interface generates serpentinite mélanges that contain blocks of high-pressure (HP) or ultrahigh-pressure (UHP) metamorphic rock derived from the subducting slab. Such serpentinite mélanges may provide a means of exhumation of HP/UHP metamorphic rocks, and may define the lower limit of locked regions on the subduction interface that fail in large earthquakes. We review recently obtained frictional strength data for brucite and talc over the temperature range 25-400°C at 100 MPa effective normal stress and compare them with new data for antigorite. These minerals respond to heating in different ways, causing their frictional strengths to diverge. Water-saturated antigorite strength increases toward the fixed dry value of μ ≈ 0.75-0.80 with heating: μ ≈ 0.50 at 25°C and μ > 0.60 at 400°C. The difference in μ between dry and watersaturated talc gouge also decreases with increasing temperature, but both the dry and watersaturated values of μ are lower at elevated temperatures. For dry talc, μ decreases from 0.35 to 0.25 between 25° and 300°C, whereas for water-saturated talc, μ is approximately 0.20 at 25°C and 0.10-0.15 at elevated temperatures. Weakening of the interlayer bond of talc with heating may be responsible for the overall reduction in its frictional strength. The strength of dry brucite also is fixed at μ = 0.45-0.50, but the water-saturated value of μ decreases from ≈0.30 at 25°C to 0.20-0.25 at 200°-400°C. The water-saturated brucite gouge has extensively recrystallized along the shear surfaces, and its weakening may be attributable to solution-transfer processes.
Researchers examined the fish communities of the Colorado and San Juan river interface zones of Lake Powell during 1999 and 2000. The objectives were to: (1) search for young razorback sucker and Colorado pikeminnow and, if found, RIT tag them; (2) examine the effectiveness of the various collection techniques on juvenile native fish; and (3) describe the fish communities found in these transitory, or actually a??migratorya?? habitats. The San Juan River interface zone was sampled 5 time each year during the spring and summer while the Colorado River site was sampled a total of 3 times over the 2-year period, all in the springa?|
The U.S. Geological Survey Downsizer is a computer application that selects, downloads, verifies, and formats station-based time-series data for environmental-resource models, particularly the Precipitation-Runoff Modeling System. Downsizer implements the client-server software architecture. The client presents a map-based, graphical user interface that is intuitive to modelers; the server provides streamflow and climate time-series data from over 40,000 measurement stations across the United States. This report is the Downsizer user's manual and provides (1) an overview of the software design, (2) installation instructions, (3) a description of the graphical user interface, (4) a description of selected output files, and (5) troubleshooting information.
A continually increasing number of high-quality digital strong-motion records from stations of the National Strong-Motion Project (NSMP) of the U.S. Geological Survey (USGS), as well as data from regional seismic networks within the United States, call for automated processing of strong-motion records with human review limited to selected significant or flagged records. The NSMP has developed the Processing and Review Interface for Strong Motion data (PRISM) software to meet this need. In combination with the Advanced National Seismic System Quake Monitoring System (AQMS), PRISM automates the processing of strong-motion records. When used without AQMS, PRISM provides batch-processing capabilities. The PRISM version 1.0.0 is platform independent (coded in Java), open source, and does not depend on any closed-source or proprietary software. The software consists of two major components: a record processing engine and a review tool that has a graphical user interface (GUI) to manually review, edit, and process records. To facilitate use by non-NSMP earthquake engineers and scientists, PRISM (both its processing engine and review tool) is easy to install and run as a stand-alone system on common operating systems such as Linux, OS X, and Windows. PRISM was designed to be flexible and extensible in order to accommodate new processing techniques. This report provides a thorough description and examples of the record processing features supported by PRISM. All the computing features of PRISM have been thoroughly tested.
This report describes a Graphical-User Interface (GUI) for SUTRA, the U.S. Geological Survey (USGS) model for saturated-unsaturated variable-fluid-density ground-water flow with solute or energy transport,which combines a USGS-developed code that interfaces SUTRA with Argus ONE, a commercial software product developed by Argus Interware. This product, known as Argus Open Numerical Environments (Argus ONE TM ), is a programmable system with geographic-information-system-like (GIS-like) functionality that includes automated gridding and meshing capabilities for linking geospatial information with finite-difference and finite-element numerical model discretizations. The GUI for SUTRA is based on a public-domain Plug-In Extension (PIE) to Argus ONE that automates the use of ArgusONE to: automatically create the appropriate geospatial information coverages (information layers) for SUTRA, provide menus and dialogs for inputting geospatial information and simulation control parameters for SUTRA, and allow visualization of SUTRA simulation results. Following simulation control data and geospatial data input bythe user through the GUI, ArgusONE creates text files in a format required for normal input to SUTRA,and SUTRA can be executed within the Argus ONE environment. Then, hydraulic head, pressure, solute concentration, temperature, saturation and velocity results from the SUTRA simulation may be visualized. Although the GUI for SUTRA discussed in this report provides all of the graphical pre- and post-processor functions required for running SUTRA, it is also possible for advanced users to apply programmable features within Argus ONE to modify the GUI to meet the unique demands of particular ground-water modeling projects.