Search USGS⌕ Search

SEARCH · Search USGS

Results for “Statistical Science”

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

Volcano deformation: Insights into magmatic systems

Volcano geodesy is the branch of geodetic science that deals with the changing shapes of volcanoes, whether large or small, deep-seated or surficial. Together with seismicity and volcanic gas flux, deformation of the ground surface can be a key indicator of subsurface conditions and processes at volcanoes—information that not only improves scientific understanding of magmatic systems but also is useful for assessing volcano hazards and mitigating their potential consequences. To take full advantage of such information requires detailed characterization of the deformation field in space and time. Currently, no single geodetic technique is capable of providing both the high spatial and temporal resolution required. However, important advances have been made recently by combining information from real-time in situ sensors such as continuous GPS, strainmeters, and tiltmeters with repeated campaign-style GPS, microgravity, interferometric synthetic aperture radar (InSAR), lidar, and photogrammetric observations. Continuous real-time data constrain the timing but not necessarily the spatial extent and pattern of deformation, whereas InSAR provides detailed spatial information but only at intervals of several days to weeks. Repeated microgravity surveys, when combined with independent measurements of surface height, are uniquely sensitive to changes in subsurface mass distribution and therefore can be used to distinguish among processes driven by magma, hydrous fluids, or gas. Simultaneous analysis of multiple geodetic datasets, especially when guided by information from global volcano databases and numerical simulations of volcanic processes and products, can provide a statistical basis for outcome prediction and serve as a guide for additional observations. In the foreseeable future, interactions among magmatic, tectonic, and hydrothermal systems will be monitored and modeled at regional scale in real time, enabling new insights into Earth’s subsurface environment.

Book chapter↗
USGS2019

State-space analysis of power to detect regional brook trout population trends over time

Threats to aquatic biodiversity are expressed at broad spatial scales, but identifying regional trends in abundance is challenging owing to variable sampling designs, and temporal and spatial variation in abundance. We compiled a regional dataset of brook trout Salvelinus fontinalis counts across their southern range representing 326 sites from eight states between 1982-2014, and conducted a statistical power analysis using Bayesian state-space models to evaluate the ability to detect temporal trends by characterizing posterior distributions with three approaches. A combination of monitoring periods, number of sites and electrofishing passes, decline magnitude and different revisit patterns were tested. Power increased with monitoring periods and decline magnitude. Trends in adults were better detected than young-of-the-year fish, which showed greater inter-annual variation in abundance. The addition of weather covariates to account for the temporal variation increased power only slightly. Single- and three-pass electrofishing methods were similar in power. Finally, power was higher for sampling designs with more frequent revisits over the duration of the monitoring program. Our results provide guidance for broad-scale monitoring designs for temporal trend detection.

Canadian Journal of Fisheries and Aquatic Sciences↗
USGS2021

A framework for assessing the ability to detect macroscale effects on fish growth

Various abiotic and biotic factors affect fish and their habitats at macroscales. For example, changes in global temperatures will likely alter demographic rates, including growth. However, to date, there is no statistical framework for assessing the ability to detect macroscale effects on fish growth under different sampling scenarios. We provide a generalized framework for calculating the frequentist and Bayesian power of detecting macroscale effects on fish growth. We illustrate this framework for a range of sampling scenarios that varied in the number of fish sampled per lake, the number of lakes sampled, and the magnitude of the temperature effect on growth for two case study species. However, the framework can be adapted to investigate other species, sampling scenarios, and environmental drivers. The ability to detect macroscale effects was more affected by the number of lakes sampled rather than the number of fish sampled from each lake. Confidently detecting macroscale effects likely requires sampling hundreds of lakes. This was true for both case study species, despite different life histories and extents of spatial variability in growth.

Canadian Journal of Fisheries and Aquatic Sciences↗
USGS2015

Higher-order statistical moments and a procedure that detects potentially anomalous years as two alternative methods describing alterations in continuous environmental data

Statistics of central tendency and dispersion may not capture relevant or desired characteristics of the distribution of continuous phenomena and, thus, they may not adequately describe temporal patterns of change. Here, we present two methodological approaches that can help to identify temporal changes in environmental regimes. First, we use higher-order statistical moments (skewness and kurtosis) to examine potential changes of empirical distributions at decadal extents. Second, we adapt a statistical procedure combining a non-metric multidimensional scaling technique and higher density region plots to detect potentially anomalous years. We illustrate the use of these approaches by examining long-term stream temperature data from minimally and highly human-influenced streams. In particular, we contrast predictions about thermal regime responses to changing climates and human-related water uses. Using these methods, we effectively diagnose years with unusual thermal variability and patterns in variability through time, as well as spatial variability linked to regional and local factors that influence stream temperature. Our findings highlight the complexity of responses of thermal regimes of streams and reveal their differential vulnerability to climate warming and human-related water uses. The two approaches presented here can be applied with a variety of other continuous phenomena to address historical changes, extreme events, and their associated ecological responses.

Hydrology and Earth System Sciences↗
USGS2026

Spatio-temporal modeling for assessing geoenergy resources: A workflow applied to gas in place variation in coal beds

The ability to estimate spatio-temporal changes in hydrocarbon reservoir properties and energy resources within pore volumes is essential for optimizing production, reservoir management, geologic energy storage, and safety in underground mining operations. In coal seams, predicting remaining methane gas-in-place (GIP) is critical for quantifying producible gas and improving mine safety and productivity through effective ventilation planning. Although such changes are commonly evaluated using physics-based numerical simulation models, these approaches often require extensive data, calibration effort, and time. This study presents a spatio-temporal geostatistical modeling approach that bridges the gap between purely spatial models and full numerical simulations. The method is applied to a case study of coal seam degasification in the Mary Lee coal group, Black Warrior Basin, Alabama, USA, to estimate GIP evolution over time within a selected mining district. The analysis uses published data from prior natural gas production history-matching of degasification using vertical wells. Empirical spatial and temporal statistics were calculated for reservoir pressure and water saturation, and spatio-temporal variogram models were fitted to experimental variograms. These models provided the structural basis for spatio-temporal kriging, integrated with spatial estimates of time-invariant parameters (porosity, density, and thickness) to estimate GIP. This approach enabled estimation of GIP changes over time, including periods without data. Boxplots of GIP estimates indicated systematic depletion and decreasing spatial variability, reflecting the impacts of degasification. Comparison with cumulative gas production from empirical well records showed approximately 85% agreement based on a relative similarity metric. Spatio-temporal GIP estimates were also used to estimate methane emissions to longwall ventilation systems and compared with reported emissions from the U.S. EPA Greenhouse Gas Reporting Program, showing similar distributions (≈80%) given data limitations. Overall, this integrated modeling approach provides time-dependent GIP estimates with broader implications for resource assessment applications.

Alabama↗
USGS2005

Trends in hydrophobic organic contaminants in urban and reference lake sediments across the United States, 1970-2001

A shift in national policy toward stronger environmental protection began in the United States in about 1970. Conversely, urban land use, population, energy consumption, and vehicle use have increased greatly since then. To assess the effects of these changes on water quality, the U.S. Geological Survey used sediment cores to reconstruct water-quality histories for38 urban and reference lakes across the United States. Cores were age-dated, and concentration profiles of polycyclic aromatic hydrocarbons (PAHs) and chlorinated hydrocarbons were tested statistically. Significant trends in total DDT, p,p???-DDE, and total PCBs were all downward. Trends in chlordane were split evenly between upward and downward, and trends in PAHs were mostly upward. Significant trends did not occur in about one-half of cases tested. Concentrations of p,p???-DDE, p,p???-DDD, and PCBs were about one-half as likely to exceed the probable effect concentration (PEC), a sediment quality guideline, in sediments deposited in the 1990s as in 1965-1975, whereas PAHs were twice as likely to exceed the PEC in the more recently deposited sediments. Concentrations of all contaminants evaluated correlated strongly with urban land use. Upward trends in PAH concentrations, the strong association of PAH with urban settings, and rapid urbanization occurring in the United States suggest that PAHs could surpass chlorinated hydrocarbons in the threat they pose to aquatic biota in urban streams and lakes.

Environmental Science & Technology↗
USGS2011

Social network models predict movement and connectivity in ecological landscapes

Network analysis is on the rise across scientific disciplines because of its ability to reveal complex, and often emergent, patterns and dynamics. Nonetheless, a growing concern in network analysis is the use of limited data for constructing networks. This concern is strikingly relevant to ecology and conservation biology, where network analysis is used to infer connectivity across landscapes. In this context, movement among patches is the crucial parameter for interpreting connectivity but because of the difficulty of collecting reliable movement data, most network analysis proceeds with only indirect information on movement across landscapes rather than using observed movement to construct networks. Statistical models developed for social networks provide promising alternatives for landscape network construction because they can leverage limited movement information to predict linkages. Using two mark-recapture datasets on individual movement and connectivity across landscapes, we test whether commonly used network constructions for interpreting connectivity can predict actual linkages and network structure, and we contrast these approaches to social network models. We find that currently applied network constructions for assessing connectivity consistently, and substantially, overpredict actual connectivity, resulting in considerable overestimation of metapopulation lifetime. Furthermore, social network models provide accurate predictions of network structure, and can do so with remarkably limited data on movement. Social network models offer a flexible and powerful way for not only understanding the factors influencing connectivity but also for providing more reliable estimates of connectivity and metapopulation persistence in the face of limited data.

Proceedings of the National Academy of Sciences of↗
USGS2019

A stratigraphic approach to inferring depositional ages from detrital geochronology data

With the increasing use of detrital geochronology data for provenance analyses, we have also developed new constraints on the age of otherwise undateable sedimentary deposits. Because a deposit can be no older than its youngest mineral constituent, the youngest defensible detrital mineral age defines the maximum depositional age of the sampled bed. Defining the youngest `defensible' age in the face of uncertainty (e.g., analytical and geological uncertainty, or sample contamination) is challenging. The current standard practice of finding multiple detrital minerals with indistinguishable ages provides confidence that a given age is not an artifact; however, we show how requiring this overlap reduces the probability of identifying the true youngest component age. Barring unusual complications, the principle of superposition dictates that sedimentary deposits must get younger upsection. This fundamental constraint can be incoporated into the analysis of depositional ages in sedimentary sections through the use of Bayesian statistics, allowing for the inference of bounded estimates of true depositional ages and uncertainties from detrital geochronology so long as some minimum age constraints are present. We present two approaches for constructing a Bayesian model of deposit ages, first solving directly for the ages of deposits with the prior constraint that the ages of units must obey stratigraphic ordering, and second describing the evolution of ages with a curve that represents the sediment accumulation rate. Using synthetic examples we highlight how this method preforms in less-than-ideal circumstances. In an example from the Magallanes Basin of Patagonia, we demonstrate how introducing other age information from the stratigraphic section (e.g., fossil assemblages or radiometric dates) and formalizing the stratigraphic context of samples provides additional constraints on and information regarding depositional ages or derived quantities (e.g., sediment accumulation rates) compared to isolated analysis of individual samples.

Frontiers in Earth Science↗
USGS2007

Surface-Water Conditions in Georgia, Water Year 2005

INTRODUCTION The U.S. Geological Survey (USGS) Georgia Water Science Center-in cooperation with Federal, State, and local agencies-collected surface-water streamflow, water-quality, and ecological data during the 2005 Water Year (October 1, 2004-September 30, 2005). These data were compiled into layers of an interactive ArcReaderTM published map document (pmf). ArcReaderTM is a product of Environmental Systems Research Institute, Inc (ESRI?). Datasets represented on the interactive map are * continuous daily mean streamflow * continuous daily mean water levels * continuous daily total precipitation * continuous daily water quality (water temperature, specific conductance dissolved oxygen, pH, and turbidity) * noncontinuous peak streamflow * miscellaneous streamflow measurements * lake or reservoir elevation * periodic surface-water quality * periodic ecological data * historical continuous daily mean streamflow discontinued prior to the 2005 water year The map interface provides the ability to identify a station in spatial reference to the political boundaries of the State of Georgia and other features-such as major streams, major roads, and other collection stations. Each station is hyperlinked to a station summary showing seasonal and annual stream characteristics for the current year and for the period of record. For continuous discharge stations, the station summary includes a one page graphical summary page containing five graphs, a station map, and a photograph of the station. The graphs provide a quick overview of the current and period-of-record hydrologic conditions of the station by providing a daily mean discharge graph for the water year, monthly statistics graph for the water year and period of record, an annual mean streamflow graph for the period of record, an annual minimum 7-day average streamflow graph for the period of record, and an annual peak streamflow graph for the period of record. Additionally, data can be accessed through the layer's link to the National Water Inventory System Web (NWISWeb) Interface.

Georgia↗
USGS2002

Spatio-temporal dynamics of species richness in coastal fish communities

Determining patterns of change in species richness and the processes underlying the dynamics of biodiversity are of key interest within the field of ecology, but few studies have investigated the dynamics of vertebrate communities at a decadal temporal scale. Here, we report findings on the spado-temporal variability in the richness and composition of fish communities along the Norwegian Skagerrak coast having been surveyed for more than half a century. Using statistical models incorporating non-detection and associated sampling variance, we estimate local species richness and changes in species composition allowing us to compute temporal variability in species richness. We tested whether temporal variation could be related to distance to the open sea and to local levels of pollution. Clear differences in mean species richness and temporal variability are observed between fjords that were and were not exposed to the effects of pollution. Altogether this indicates that the fjord is an appropriate scale for studying changes in coastal fish communities in space and time. The year-to-year rates of local extinction and turnover were found to be smaller than spatial differences in community composition. At the regional level, exposure to the open sea plays a homogenizing role, possibly due to coastal currents and advection.

Proceedings of the Royal Society B: Biological Sci↗
USGS2014

Sensor data as a measure of native freshwater mussel impact on nitrate formation and food digestion in continuous-flow mesocosms

Native freshwater mussels can influence the aquatic N cycle, but the mechanisms and magnitude of this effect are not fully understood. We assessed the effects of Amblema plicata and Lampsilis cardium on N transformations over 72 d in 4 continuous-flow mesocosms, with 2 replicates of 2 treatments (mesocosms with and without mussels), equipped with electronic water-chemistry sensors. We compared sensor data to discrete sample data to assess the effect of additional sensor measurements on the ability to detect mussel-related effects on NO 3 – formation. Analysis of 624 sensor-based data points detected a nearly 6% increase in NO 3 – concentration in overlying water of mesocosms with mussels relative to mesocosms without mussels (p < 0.05), whereas analysis of 36 discrete samples showed no statistical difference in NO 3 – between treatments. Mussels also significantly increased NO 2 – concentrations in the overlying water, but no significant difference in total N was observed. We used the sensor data for phytoplankton-N and NH 4 + to infer that digestion times in mussels were 13 ± 6 h. The results suggest that rapid increases in phytoplankton-N levels in the overlying water can lead to decreased lag times between phytoplankton-N and NH 4 + maxima. This result indicates that mussels may adjust their digestion rates in response to increased levels of food. The adjustment in digestion time suggests that mussels have a strong response to food availability that can disrupt typical circadian rhythms. Use of sensor data to measure directly and to infer mussel effects on aquatic N transformations at the mesocosm scale could be useful at larger scales in the future.

Freshwater Science↗
USGS2012

Archive of digital chirp subbottom profile data collected during USGS cruise 10BIM04 offshore Cat Island, Mississippi, September 2010

In September of 2010, the U.S. Geological Survey (USGS), in cooperation with the U.S. Army Corps of Engineers (USACE), conducted a geophysical survey to investigate the geologic controls on barrier island framework of Cat Island, Miss., as part of a broader USGS study on Barrier Island Mapping (BIM). These surveys were funded through the Mississippi Coastal Improvements Program (MsCIP) and the Northern Gulf of Mexico (NGOM) Ecosystem Change and Hazard Susceptibility Project as part of the Holocene Coastal Evolution of the Mississippi-Alabama Region Subtask. This report serves as an archive of unprocessed digital chirp subbottom data, trackline maps, navigation files, GIS files, Field Activity Collection System (FACS) logs, and formal FGDC metadata. Gained (showing a relative increase in signal amplitude) digital images of the seismic profiles are also provided. Refer to the Acronyms page for expansions of acronyms and abbreviations used in this report. The USGS Saint Petersburg Coastal and Marine Science Center (SPCMSC) assigns a unique identifier to each cruise or field activity. For example, 10BIM04 tells us the data were collected in 2010 during the fourth field activity for that project in that calendar year. Refer to http://walrus.wr.usgs.gov/infobank/programs/html/definition/activity.html for a detailed description of the method used to assign the field activity identification (ID). All chirp systems use a signal of continuously varying frequency; the EdgeTech SB-512i system used during this survey produces high-resolution, shallow-penetration (typically less than 50 milliseconds (ms)) profile images of sub-seafloor stratigraphy. The towfish contains a transducer that transmits and receives acoustic energy; it was housed within a float system (built at the SPCMSC), which allows the towfish to be towed at a constant depth of 1.07 meters (m) below the sea surface. As transmitted acoustic energy intersects density boundaries, such as the seafloor or sub-surface sediment layers, some energy is reflected back toward the transducer, received, and recorded by a Personal Computer (PC)-based seismic acquisition system. This process is repeated at regular time intervals (for example, 0.125 seconds (s)), and returned energy is recorded for a specific duration (for example, 50 ms). In this way, a two-dimensional (2-D) vertical image of the shallow geologic structure beneath the ship track is produced. Figure 1 displays the acquisition geometry. Refer to table 1 for a summary of acquisition parameters and table 2 for trackline statistics. The archived trace data are in standard Society of Exploration Geophysicists (SEG) SEG Y rev. 0 format (Barry and others, 1975); the first 3,200 bytes of the card image header are in American Standard Code for Information Interchange (ASCII) format instead of Extended Binary Coded Decimal Interchange Code (EBCDIC) format. The SEG Y files may be downloaded and processed with commercial or public domain software such as Seismic Unix (SU) (Cohen and Stockwell, 2010). See the How To Download SEG Y Data page for download instructions. The printable profiles provided here are GIF images that were processed and gained using SU software, and they can be viewed from the Profiles page or from links located on the trackline maps; refer to the Software page for links to example SU processing scripts. The SEG Y files are available on the DVD version of this report or on the Web, downloadable via the USGS Coastal and Marine Geoscience Data System (http://cmgds.marine.usgs.gov). The data are also available for viewing using GeoMapApp (http://www.geomapapp.org) and Virtual Ocean (http://www.virtualocean.org) multi-platform open source software.

Mississippi↗
USGS2018

Wetlands inform how climate extremes influence surface water expansion and contraction

Effective monitoring and prediction of flood and drought events requires an improved understanding of how and why surface water expansion and contraction in response to climate varies across space. This paper sought to (1) quantify how interannual patterns of surface water expansion and contraction vary spatially across the Prairie Pothole Region (PPR) and adjacent Northern Prairie (NP) in the United States, and (2) explore how landscape characteristics influence the relationship between climate inputs and surface water dynamics. Due to differences in glacial history, the PPR and NP show distinct patterns in regards to drainage development and wetland density, together providing a diversity of conditions to examine surface water dynamics. We used Landsat imagery to characterize variability in surface water extent across 11 Landsat path/rows representing the PPR and NP (images spanned 1985–2015). The PPR not only experienced a 2.6-fold greater surface water extent under median conditions relative to the NP, but also showed a 3.4-fold greater change in surface water extent between drought and deluge conditions. The relationship between surface water extent and accumulated water availability (precipitation minus potential evapotranspiration) was quantified per watershed and statistically related to variables representing hydrology-related landscape characteristics (e.g., infiltration capacity, surface storage capacity, stream density). To investigate the influence stream connectivity has on the rate at which surface water leaves a given location, we modeled stream-connected and stream-disconnected surface water separately. Stream-connected surface water showed a greater expansion with wetter climatic conditions in landscapes with greater total wetland area, but lower total wetland density. Disconnected surface water showed a greater expansion with wetter climatic conditions in landscapes with higher wetland density, lower infiltration and less anthropogenic drainage. From these findings, we can expect that shifts in precipitation and evaporative demand will have uneven effects on surface water quantity. Accurate predictions regarding the effect of climate change on surface water quantity will require consideration of hydrology-related landscape characteristics including wetland storage and arrangement.

Prairie Pothole Regino↗
USGS2018

Ecosystem classification and mapping of the Laurentian Great Lakes

Owing to the enormity and complexity of the Laurentian Great Lakes, an ecosystem classification is needed to better understand, protect, and manage this largest freshwater ecosystem in the world. Using a combination of statistical analyses, published knowledge, and expert opinion, we identified key driving variables and their ecologically relevant thresholds and delineated and mapped aquatic systems for the entire Great Lakes. We identified and mapped 77 aquatic ecological units (AEUs) that depict unique combinations of depth, thermal regime, hydraulic, and landscape classifiers. Those 77 AEU types were distributed across 1997 polygons (patches) ranging from 1 to >48 000 km 2 in area and were most diverse in the nearshore (35 types), followed by the coastal margin (26), and then the offshore (16). Our classification and mapping of ecological units captures gradients that characterize types of aquatic systems in the Great Lakes and provides a geospatial accounting framework for resource inventory, status and trend assessment; research for ecosystem questions; and management and policy-making.

Laurentian Great Lakes↗
USGS2017

Uptake and disposition of select pharmaceuticals by bluegill exposed at constant concentrations in a flow-through aquatic exposure system

The increasing use of pharmaceuticals has led to their subsequent input into and release from wastewater treatment plants, with corresponding discharge into surface waters that may subsequently exert adverse effects upon aquatic organisms. Although the distribution of pharmaceuticals in surface water has been extensively studied, the details of uptake, internal distribution, and kinetic processing of pharmaceuticals in exposed fish have received less attention. For this research, we investigated the uptake, disposition, and toxicokinetics of five pharmaceuticals (diclofenac, methocarbamol, rosuvastatin, sulfamethoxazole, and temazepam) in bluegill sunfish ( Lepomis macrochirus ) exposed to environmentally relevant concentrations (1000–4000 ng L –1 ) in a flow-through exposure system. Temazepam and methocarbamol were consistently detected in bluegill biological samples with the highest concentrations in bile of 4, 940, and 180 ng g –1 , respectively, while sulfamethoxazole, diclofenac, and rosuvastatin were only infrequently detected. Over 30-day exposures, the relative magnitude of mean concentrations of temazepam and methocarbamol in biological samples generally followed the order: bile ≫ gut > liver and brain > muscle, plasma, and gill. Ranges of bioconcentration factors (BCFs) in different biological samples were 0.71–3960 and 0.13–48.6 for temazepam and methocarbamol, respectively. Log BCFs were statistically positively correlated to pH adjusted log K ow (that is, log D ow ), with the strongest relations for liver and brain ( r 2 = 0.92 and 0.99, respectively), implying that bioconcentration patterns of ionizable pharmaceuticals depend on molecular status, that is, whether a pharmaceutical is un-ionized or ionized at ambient tissue pH. Methocarbamol and temazepam underwent rapid uptake and elimination in bluegill biological compartments with uptake rate constants ( K u ) and elimination rate constants ( K e ) at 0.0066–0.0330 h –1 and 0.0075–0.0384 h –1 , respectively, and half-lives at 18.1–92.4 h. Exposure to mixtures of diclofenac, methocarbamol, sulfamethoxazole, and temazepam had little or no influence on the uptake and elimination rates, suggesting independent multiple uptake and disposition behaviors of pharmaceuticals by fish would occur when exposed to effluent-influenced surface waters.

Environmental Science & Technology↗
USGS2020

A multi-model approach toward understanding iron fouling at rock-fill drainage sites along roadways in New Hampshire, USA

Factors affecting iron fouling in wet areas adjacent to roadways were investigated by collecting field rock cut and aqueous physicochemical data; developing exploratory predictive models; and developing geochemical models. Basic data included the identification of iron fouling from aerial imagery and field visits at 374 New Hampshire rock cut locations, and their associated rock-fill sites. Based on field water quality measurements from wet areas at 36 of the rock-fill sites, the occurrence of iron fouling was associated with higher values of specific conductance, lower concentrations of dissolved oxygen and lower pH compared to areas without iron fouling. A statistical model, using boosted regression trees, was developed to predict the occurrence of iron fouling in wet areas adjacent to roadways where rock-fill from nearby rock cuts was used in roadway construction. The model was used to develop a continuous iron fouling probability map for the state of New Hampshire that can be used to better understand the occurrence of iron fouling. Geochemical models illustrate how iron fouling of waters increases along roadways built with fill from sulfidic rock cuts as a result of acid generation from pyrite dissolution and ferrous iron (Fe 2+ ) oxidation and increases in areas with greater specific conductance from deicing runoff caused by cation exchange. More iron is precipitated as goethite in simulations that include pyrite, and in simulations with deicing salts added, indicating that rock-fill sites with rocks that contain pyrite and water with greater salt content could have enhanced iron fouling.

New Hampshire↗
USGS2020

Fluoride occurrence in United States groundwater

Data from 38,105 wells were used to characterize fluoride (F) occurrence in untreated United States (U.S.) groundwater. For domestic wells (n = 11,032), water from which is generally not purposely fluoridated or monitored for quality, 10.9% of the samples have F concentrations >0.7 mg/L (U.S. Public Health Service recommended optimal F concentration in drinking water for preventing tooth decay) (87% are <0.7 mg/L); 2.6% have F > 2 mg/L (EPA Secondary Maximum Contaminant Level, SMCL); and 0.6% have F > 4 mg/L (EPA MCL). The data indicate the biggest concern with F in domestic wells at the national scale could be one of under consumption of F with respect to the oral-health benchmark (0.7 mg/L). Elevated F concentrations relative to the SMCL and MCL are regionally important, particularly in the western U.S. Statistical comparisons of potentially important controlling factors in four F-concentration categories (<0.1–0.7 mg/L; >0.7–2 mg/L; >2–4 mg/L; >4 mg/L) at the national scale indicate the highest F-concentration category is associated with groundwater that has significantly greater pH values, TDS and alkalinity concentrations, and well depths, and lower Ca/Na ratios and mean annual precipitation, than the lowest F-concentration category. The relative importance of the controlling factors appears to be regionally variable. Three case studies illustrate the spatial variability in controlling factors using groundwater-age (groundwater residence time), water-isotope (evaporative concentration), and water-temperature (geothermal processes) data. Populations potentially served by domestic wells with F concentrations <0.7, >0.7, >2, and >4 mg/L are estimated to be ~28,200,000, ~3,110,000; ~522,000; and ~172,000 people, respectively, in 40 principal aquifers with at least 25 F analyses per aquifer.

Science of the Total Environment↗