SEARCH · USGS Science

Results for “Research Square”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18Linked to original sources

USGS2005

Seismicity and tilt associated with the 2003 Anatahan eruption sequence

On May 10, 2003, the first historical eruption of Anatahan volcano in the western Pacific Mariana Islands was fortuitously recorded by a broadband seismograph installed on the island only 4 days prior to the eruption. This station, located 7 km WNW of the active crater, together with another broadband seismograph on Sarigan Island 45 km to the north, continued to operate throughout the 2-month period of major eruptive activity in May and June and throughout the majority of the following year. In June 2003, the Saipan Emergency Management Office and the US Geological Survey installed two telemetered high-gain short-period seismic stations to monitor the activity in real-time. The only earthquakes detected in the 4-day period from the initial seismograph installation until 6 h prior to the eruption occurred approximately 20 km to the northeast of the island on May 8. The first volcano-tectonic (VT) event located near the volcano occurred at 01:53 GMT on May 10. The number of events per hour then increased dramatically and a period of about 80 discrete earthquakes per hour commenced at about 06:20 GMT, immediately prior to the estimated eruption time of 07:30 from the Volcanic Ash Advisory Center. A long-period tilt signal recorded on the horizontal components of the broadband seismograph, indicating upward movement of the crater region, also commenced at about 06:20. Inflation continued until 09:30, when the direction of tilt reversed. Deflation continued until 17:50, coinciding with a reduction in the number of VT events. The larger VT events were located with a linearized least-squares location algorithm. Magnitudes of located VT events on May 10 ranged from 2.0 to 3.2, but a period of larger VT events were recorded on May 11, with the largest M 4.2. After about 36 h of intense earthquake activity, the number of discrete VT events declined and was replaced by nearly continuous volcanic tremor for the next 6 weeks. Differing types of very long-period events may suggest complex non-destructive magmatic source mechanisms, persisting dominantly throughout the first 10 days of the initial eruption. Visual reports indicate that a small craggy dome extruded sometime between May 20 and June 5. From analogy with other volcanic dome extrusions, we believe the dome probably extruded just as the tremor amplitude decreased dramatically about May 24. This dome was then destroyed between June 13 and 14. Reduced displacement of the co-eruption tremor is estimated as approximately 40-80 cm2, suggesting an eruption with a Volcanic Explosivity Index of about 3. 2005 Elsevier B.V. All rights reserved.

Northern Mariana Islands↗

USGS2010

Distribution of the non-native gastropod Melanoides tuberculatus in Biscayne National Park, Florida

Introduction Melanoides tuberculatus (fig. 1), a gastropod that is not native to South Florida, was identified in Biscayne National Park (BNP) while researchers from the U.S. Geological Survey were conducting other studies around the Black Point canals in the summer of 2003. A study to determine the distribution, genetics, and salinity tolerance of this freshwater species began in 2004. For park managers and the recreational users of BNP, the presence of Melanoides tuberculatus is cause for concern because it is the intermediate host for several trematode parasites that affect humans and animals in multiple ways: Human health concerns include parasites (trematode worms) that affect the lung (Paragonimus westermani) and the liver ( Clonorchis sinensis , and Opisthorchis sp. ). Skin irritations and lesions can also be caused by trematodes that normally infect other species. Animal health concerns include parasites that affect the eyes of waterfowl ( Philophthalmus megalurus ), a trematode that burrows into the cartilage of fish and can lead to death ( Centrocestus formosanus ), and a trematode that infects the muscle tissue of fish and causes multiple abnormalities ( Haplorchis sp. ). When parasite-infected fish or crustaceans are eaten by birds or mammals, the next stage of the life cycle is in place (fig. 2). Native snail species concerns include the potential displacement of invertebrate species that make up the natural nearshore benthic communities in BNP. M. tuberculatus has very high reproductive rates, reproduces via parthenogenesis, and is live bearing, factors which increase the survivability of the offspring and increase competition for limited resources. These snails are considered to be freshwater animals in their native habitat of Southeast Asia. However, they have been collected in BNP in both estuarine and marine waters along the western margins of BNP and, as far as 1.7 kilometers (km) from shore at the Black Point canal inflow into Biscayne Bay (383 live per square meter (/m 2 ). In BNP, M. tuberculatus is a benthic inhabitant grazing on micro algal components at the sediment surface. A documented population with as many as 23,000/m 2 was observed at Snapper Creek, near Coral Gables, FL (Roessler and others, 1977), north of BNP.

Florida↗

USGS1997

Daily air temperature interpolated at high spatial resolution over a large mountainous region

Two methods are investigated for interpolating daily minimum and maximum air temperatures ( T min and T max ) at a 1 km spatial resolution over a large mountainous region (830000 km 2 ) in the U.S. Pacific Northwest. The methods were selected because of their ability to (1) account for the effect of elevation on temperature and (2) efficiently handle large volumes of data. The first method, the neutral stability algorithm (NSA), used the hydrostatic and potential temperature equations to convert measured temperatures and elevations to sea-level potential temperatures. The potential temperatures were spatially interpolated using an inverse-squared-distance algorithm and then mapped to the elevation surface of a digital elevation model (DEM). The second method, linear lapse rate adjustment (LLRA), involved the same basic procedure as the NSA, but used a constant linear lapse rate instead of the potential temperature equation. Cross-validation analyses were performed using the NSA and LLRA methods to interpolate T min and T max each day for the 1990 water year, and the methods were evaluated based on mean annual interpolation error (IE). The NSA method showed considerable bias for sites associated with vertical extrapolation. A correction based on climate station/grid cell elevation differences was developed and found to successfully remove the bias. The LLRA method was tested using 3 lapse rates, none of which produced a serious extrapolation bias. The bias-adjusted NSA and the 3 LLRA methods produced almost identical levels of accuracy (mean absolute errors between 1.2 and 1.3°C), and produced very similar temperature surfaces based on image difference statistics. In terms of accuracy, speed, and ease of implementation, LLRA was chosen as the best of the methods tested.

Columbia River basin↗

USGS2024

Self-potential tomography preconditioned by particle swarm optimization— Application to monitoring hyporheic exchange in a bedrock river

A self-potential (SP) data-inversion algorithm was developed and tested on an analytical model of electrical-potential profile data attributed to single and multiple polarized electrical sources. The developed algorithm was then validated by an application to SP-monitoring field data measured on the floodplain of East Fork Poplar Creek, Oak Ridge, Tennessee, to image electrical sources in areas conducive to preferential flow into the flood plain from the bedrock-lined riverbed. The algorithm combined stochastic source-localization by particle-swarm-optimization (PSO) of electrical sources characterized by simplified geometries with source tomography by regularized weighted least-squares minimization of a quadratic objective function. Prior information was incorporated by preconditioning the tomography algorithm by PSO results. Variable percentages of random noise were added to analytical-model data to evaluate the algorithm performance. Results indicated that true parameters of single-source models were inverted and approximated with small residual error, whereas inversion of analytical-model data representing multiple electrical sources accurately approximated the locations of the sources but miscalculated some parameters because of the non-uniqueness of the inverse-model solution. Source tomography applied to analytical model data during testing produced a spatially continuous parameter field that identified the locations of point-scale synthetic dipole sources of electrical current flow with varying degrees of accuracy depending on the prior information incorporated into the tomography. When applied to SP-monitoring field data, the algorithm imaged electrical sources within a known fault that intersects the bedrock riverbed and flood plain of East Fork Poplar Creek and depicted dynamic electrical conditions attributed to hyporheic exchange.

Tennessee↗

USGS2013

Spatial and temporal patterns of dust emissions (2004-2012) in semi-arid landscapes, southeastern Utah, USA

Aeolian dust can influence nutrient availability, soil fertility, plant interactions, and water-holding capacity in both source and downwind environments. A network of 85 passive collectors for aeolian sediment spanning numerous plant communities, soil types, and land-use histories covering approximately 4000 square kilometers across southeastern Utah was used to sample horizontal emissions of aeolian sediment. The sample archive dates to 2004 and is currently the largest known record of field-scale dust emissions for the southwestern United States. Sediment flux peaked during the spring months in all plant communities (mean: 38.1 g m −2 d −1 ), related to higher, sustained wind speeds that begin in the early spring. Dust flux was lowest during the winter period (mean: 5 g m−2 d−1) when surface wind speeds are typically low. Sites dominated by blackbrush and sagebrush shrubs had higher sediment flux (mean: 19.4 g m−2 d−1) compared to grasslands (mean: 11.2 g m−2 d−1), saltbush shrublands (mean: 10.3 g m −2 d −1 ), and woodlands (mean: 8.1 g m −2 d −1 ). Contrary to other studies on dust emissions, antecedent precipitation during one, two, and three seasons prior to sample collection did not significantly influence emission rates. Physical site-scale factors controlling dust emissions were complex and varied from one vegetation type to another.

Utah↗

USGS2002

Identification and synthetic modeling of factors affecting American black duck populations

We reviewed the literature on factors potentially affecting the population status of American black ducks ( Anas rupribes ). Our review suggests that there is some support for the influence of 4 major, continental-scope factors in limiting or regulating black duck populations: 1) loss in the quantity or quality of breeding habitats; 2) loss in the quantity or quality of wintering habitats; 3) harvest, and 4) interactions (competition, hybridization) with mallards ( Anas platyrhychos ) during the breeding and/or wintering periods. These factors were used as the basis of an annual life cycle model in which reproduction rates and survival rates were modeled as functions of the above factors, with parameters of the model describing the strength of these relationships. Variation in the model parameter values allows for consideration of scientific uncertainty as to the degree each of these factors may be contributing to declines in black duck populations, and thus allows for the investigation of the possible effects of management (e.g., habitat improvement, harvest reductions) under different assumptions. We then used available, historical data on black duck populations (abundance, annual reproduction rates, and survival rates) and possible driving factors (trends in breeding and wintering habitats, harvest rates, and abundance of mallards) to estimate model parameters. Our estimated reproduction submodel included parameters describing negative density feedback of black ducks, positive influence of breeding habitat, and negative influence of mallard densities; our survival submodel included terms for positive influence of winter habitat on reproduction rates, and negative influences of black duck density (i.e., compensation to harvest mortality). Individual models within each group (reproduction, survival) involved various combinations of these factors, and each was given an information theoretic weight for use in subsequent prediction. The reproduction model with highest AIC weight (0.70) predicted black duck age ratios increasing as a function of decreasing mallard abundance and increasing acreage of breeding habitat; all models considered involved negative density dependence for black ducks. The survival model with highest AIC weight (0.51) predicted nonharvest survival increasing as a function of increasing acreage of wintering habitat and decreasing harvest rates (additive mortality); models involving compensatory mortality effects received ≈0.12 total weight, vs. 0.88 for additive models. We used the combined model, together with our historical data set, to perform a series of 1-year population forecasts, similar to those that might be performed under adaptive management. Initial model forecasts over-predicted observed breeding populations by ≈25%. Least-squares calibration reduced the bias to ≈0.5% under prediction. After calibration, model-averaged predictions over the 16 alternative models (4 reproduction × 4 survival, weighted by AIC model weights) explained 67% of the variation in annual breeding population abundance for black ducks, suggesting that it might have utility as a predictive tool in adaptive management. We investigated the effects of statistical uncertainty in parameter values on predicted population growth rates for the combined annual model, via sensitivity analyses. Parameter sensitivity varied in relation to the parameter values over the estimated confidence intervals, and in relation to harvest rates and mallard abundance. Forecasts of black duck abundance were extremely sensitive to variation in parameter values for the coefficients for breeding and wintering habitat effects. Model-averaged forecasts of black duck abundance were also sensitive to changes in harvest rate and mallard abundance, with rapid declines in black duck abundance predicted for a range of harvest rates and mallard abundance higher than current levels of either factor, but easily envisaged, particularly given current rates of growth for mallard populations. Because of concerns about sensitivity to habitat coefficients, and particularly in light of deficiencies in the historical data used to estimate these parameters, we developed a simplified model that excludes habitat effects. We also developed alternative models involving a calibration adjustment for reproduction rates, survival rates, or neither. Calibration of survival rates performed best (AIC weight 0.59, % BIAS = -0.280, R 2 =0.679), with reproduction calibration somewhat inferior (AIC weight 0.41, % BIAS = -0.267, R 2 =0.672); models without calibration received virtually no AIC weight and were discarded. We recommend that the simplified model set (4 biological models × 2 alternative calibration factors) be retained as the best working set of alternative models for research and management. Finally, we provide some preliminary guidance for the development of adaptive harvest management for black ducks, using our working set of models.

Book↗

USGS2024

Examining 22 years of ambient seismic wavefield at Mount St. Helens

An increase in seismic activity precedes most volcanic eruptions. Whereas event-based forecasting approaches have been successful, some eruptions remain unanticipated, resulting in casualties and damage. Our study leverages the recent advancements in ambient field seismology. We explore features extracted from continuous ambient fields using traditional methods, for example, peak ground velocity, peak ground acceleration, root mean square, root median square, real-time seismic amplitude measurement, and novel methods (displacement seismic amplitude ratio and spectral width). In addition, we explore unsupervised learning of higher order wavelet features using scattering networks. We find that combining all the methods was necessary to disentangle the effects of seismic sources from structural changes at Mount St. Helens. Although the ambient wavefield-based approach does not yield additional or more significant precursory signals than event-based methods at Mount St. Helens, our study demonstrates that the ambient wavefield provides supplementary information, mainly about structural changes and complements traditional methods. The ambient seismic wavefield offers additional insights into long-lasting processes. We find enhanced wave attenuation correlating with geochemical measurements. We interpret this as ongoing structural changes, such as dome growth or the evolution of the volcanic conduit system. On annual and decadal timescales, we interpret seasonal seismic attenuation in the shallow subsurface as groundwater fluctuations, corroborated by observations at the nearby Spirit Lake level. This multimethod approach at Mount St. Helens sheds light on a volcanic system’s underlying dynamics and structure.

Washington↗

USGS2025

Benthic habitat map of Olowalu Reef, Maui, Hawaii—Geomorphological structure, biological cover, and geologic zonation determined with spectral, lidar, and acoustic data

The fringing coral reef off Olowalu, Maui, Hawaii, has been identified as a local conservation priority site. In 2007, the National Oceanic and Atmospheric Administration (NOAA) produced a benthic habitat map of the Hawaiian Islands that was used as a foundation for this study. To support place-based management of the reef in the future, the U.S. Geological Survey (USGS) mapped the geologic zone, major and dominant geomorphological structure, biological cover type, and percent of biological cover for 11 square kilometers (km 2 ) of Olowalu Reef at a minimum mapping unit (MMU) of 100 square meters (m 2 ) to create a benthic habitat map. Heads-up digitization was employed on 0.50-meter (m) natural color satellite orthoimagery with ancillary 1-m acoustic backscatter imagery from single-scan sonar (sound navigation and ranging). A 1-m, 4-m, and 8-m digital bathymetric model (DBM) was interpolated from bathymetric lidar (light detection and ranging), and various geomorphometric layers derived from the DBMs were used for habitat interpretation. Still-frame imagery of the seafloor extracted from vessel-towed underwater video transects on Olowalu Reef served as ground validation points ( n =870) during active mapping and accuracy assessment points ( n =216) for thematic accuracy assessment. Thematic accuracy was cross-validated by the Hawai‘i Department of Land and Natural Resources Division of Aquatic Resources. Final thematic accuracy was 88.8 percent for major structure, 85.6 percent for dominant structure, 86.0 percent for major biological cover, and 78.6 percent for type and percent of major biological cover. Reef and hardbottom constituted 52 percent of the total mapped habitat, comprising mostly aggregate reef (31 percent) and pavement (11 percent), with large swaths of spur-and-groove (9 percent). Of this hardbottom, 17 percent was covered with moderate (10 to <50 percent) coral and 27 percent with high coral cover (50 to <90 percent). High (50 to <90 percent) macroalgae cover dominated the continuous sand sheets in offshore bank/shelf zones. The map created in this study supplements the NOAA 2007 map and expands on the observations made by USGS sampling of the reef. The NOAA 2007 map and our map differed in total areal extent by a negligible 6 m 2 and were in general thematic agreement. Our map is intended to serve as a baseline for public access, general research, local-level management, and reef change for future studies.

Hawaii↗

USGS2022

Microgravity change during the 2008-2018 Kı̄lauea summit eruption: Nearly a decade of subsurface mass accumulation

Results from nine microgravity campaigns from Kı̄lauea, Hawaiʻi, spanning most of the volcano's 2008–2018 summit eruption, indicate persistent mass accumulation at shallow levels. A weighted least squares approach is used to recover microgravity results from a network of benchmarks around Kı̄lauea's summit, eliminate instrumental drift, and restore suspected data tares. A total mass of 1.9 × 10 11 kg was determined from these microgravity campaigns to have accumulated below Kı̄lauea Caldera during 2009–2015 at an estimated depth of 1.3 km below sea level. Only a fraction of this mass is reflected in surface deformation, and this is consistent with previously reported discrepancies between subsurface mass accumulation and observed surface deformation. The discrepancy, amongst other independent evidence from gas emissions, seismicity, and continuous gravimetry, indicate densification of magma in the reservoirs below the volcano summit. This densification may have been driven by degassing through the summit vent. It is hypothesized that during the final years of the summit eruption, magma densification resulted in a buildup of pressure in the reservoirs that may have contributed to the lower East Rift Zone outbreak of 2018. The observed mass accumulation beneath Kı̄lauea could not have been detected through other techniques and illustrates the importance of microgravity measurements in volcano monitoring.

Hawai'i↗

USGS2010

Modeled and measured glacier change and related glaciological, hydrological, and meteorological conditions at South Cascade Glacier, Washington, balance and water years 2006 and 2007

Winter snow accumulation and summer snow and ice ablation were measured at South Cascade Glacier, Washington, to estimate glacier mass balance quantities for balance years 2006 and 2007. Mass balances were computed with assistance from a new model that was based on the works of other glacier researchers. The model, which was developed for mass balance practitioners, coupled selected meteorological and glaciological data to systematically estimate daily mass balance at selected glacier sites. The North Cascade Range in the vicinity of South Cascade Glacier accumulated approximately average to above average winter snow packs during 2006 and 2007. Correspondingly, the balance years 2006 and 2007 maximum winter snow mass balances of South Cascade Glacier, 2.61 and 3.41 meters water equivalent, respectively, were approximately equal to or more positive (larger) than the average of such balances since 1959. The 2006 glacier summer balance, -4.20 meters water equivalent, was among the four most negative since 1959. The 2007 glacier summer balance, -3.63 meters water equivalent, was among the 14 most negative since 1959. The glacier continued to lose mass during 2006 and 2007, as it commonly has since 1953, but the loss was much smaller during 2007 than during 2006. The 2006 glacier net balance, -1.59 meters water equivalent, was 1.02 meters water equivalent more negative (smaller) than the average during 1953-2005. The 2007 glacier net balance, -0.22 meters water equivalent, was 0.37 meters water equivalent less negative (larger) than the average during 1953-2006. The 2006 accumulation area ratio was less than 0.10, owing to isolated patches of accumulated snow that endured the 2006 summer season. The 2006 equilibrium line altitude was higher than the glacier. The 2007 accumulation area ratio and equilibrium line altitude were 0.60 and 1,880 meters, respectively. Accompanying the glacier mass losses were retreat of the terminus and reduction of total glacier area. The terminus retreated at a rate of about 13 meters per year during balance year 2006 and at a rate of about 8 meters per year during balance year 2007. Glacier area near the end of balance years 2006 and 2007 was 1.74 and 1.73 square kilometers, respectively. Runoff from the basin containing the glacier and from an adjacent nonglacierized basin was gaged during all or parts of water years 2006 and 2007. Air temperature, wind speed, precipitation, and incoming solar radiation were measured at selected locations on and near the glacier. Air-temperature over the glacier at a height of 2 meters generally was less than at the same altitude in the air mass away from the glacier. Cooling of the air by the glacier increased systematically with increasing ambient air temperature. Empirically based equations were developed to estimate 2-meter-height air temperature over the glacier at five sites from site altitude and temperature at a non-glacier reference site.

Scientific Investigations Report↗

USGS2025

Multiyear crop residue cover mapping using narrow-band vs. broad-band shortwave infrared satellite imagery

Crop residue serves an important role in agricultural systems as high levels of fractional crop residue cover ( f R ) can reduce erosion, preserve soil moisture, and build soil organic carbon. However, the ability to accurately quantify f R at scale has been limited. In this study we produced annual maps of f R for farmland in Maryland, USA using WorldView-3 (WV3) imagery paired with on-farm photographs ( n = 895) classified to f R using SamplePoint software. Univariate linear regressions were used to compare photograph f R to WV3 crop residue indices including: 1) Shortwave Infrared Normalized Difference Residue Index (SINDRI), 2) Shortwave Infrared Difference Residue Index (SIDRI), 3) Normalized Difference Tillage Index (NDTI), and 4) Shortwave Infrared Angle Index (SWIRA). SINDRI and SIDRI are based on narrow bands capable of measuring lignocellulose absorption features. NDTI and SWIRA are based on Landsat-comparable broad bands. Our findings demonstrated that SINDRI outperformed other indices in f R estimation in terms of coefficient of determination ( R 2 = 0.869) and root mean square error (RMSE = 0.111), when R 2 and RMSE were averaged across six individual years. For a univariate analysis combining five years of high-quality WV3 imagery, SINDRI again exhibited the highest f R estimation performance ( R 2 = 0.795; RMSE = 0.141), suggesting that SINDRI can map f R accurately with a singular relationship, potentially reducing the need for labor-intensive ground data collection. For broad-band indices, a multiple linear regression analysis that included a Water Index (WI) and Normalized Difference Vegetation Index (NDVI) as additional predictors increased the accuracy of f R estimation significantly, particularly for SWIRA ( R 2 = 0.767; RMSE = 0.144), but also NDTI ( R 2 = 0.654; RMSE = 0.174). Our findings suggest that while indices computed from narrow-band imagery are most accurate for f R estimation, SWIRA has the potential to improve f R estimation compared to NDTI, especially when used in conjunction with WI and NDVI. An index suite of SWIRA, WI, and NDVI can be computed with Landsat 4–9 imagery, providing a more accurate record of global f R dating back to 1982.

Soil and Tillage Research↗

USGS2023

Train, inform, borrow, or combine? Approaches to process-guided deep learning for groundwater-influenced stream temperature prediction

Although groundwater discharge is a critical stream temperature control process, it is not explicitly represented in many stream temperature models, an omission that may reduce predictive accuracy, hinder management of aquatic habitat, and decrease user confidence. We assessed the performance of a previously-described process-guided deep learning model of stream temperature in the Delaware River Basin (USA). We found lower accuracy (root mean square error [RMSE] of 1.71 versus 1.35°C) and stronger seasonal bias (absolute mean monthly bias of 1.06 vs. 0.68°C) for reaches primarily influenced by deep groundwater as compared to atmospheric conditions. We then tested four approaches for improving groundwater process representation: (a) a custom loss function leveraging the unique patterns of air and water temperature coupling characteristic of different temperature drivers, (b) inclusion of additional groundwater-relevant catchment attributes, (c) incorporation of additional process model outputs, and (d) a composite model. The custom loss function and the additional attributes significantly improved the predictive accuracy in groundwater-dominated reaches (RMSE of 1.37 and 1.26°C) and reduced the seasonal bias (absolute mean monthly bias of 0.44 and 0.48°C), but neither approach could identify holdout groundwater reaches. Variable importance analysis indicates the custom loss function nudges the model to use the existing inputs more efficiently, whereas with the added features the model relies on a broader suite of inputs. This analysis is a substantial step toward more accurately representing groundwater discharge processes in stream temperature models and will improve predictive accuracy and inform habitat management.

Delaware River Basin↗

USGS2002

A National survey of methyl tert -butyl ether and other volatile organic compounds in drinking-water sources: Results of the random source-water survey

Methyl tert -butyl ether (MTBE) was detected in source water used by 8.7 percent of randomly selected community water systems (CWSs) in the United States at concentrations that ranged from 0.2 to 20 micrograms per liter ( μ g/L). The Random Survey conducted by the U.S. Geological Survey, in cooperation with the Metropolitan Water District of Southern California and the Oregon Health & Science University, was designed to provide an assessment of the frequency of detection, concentration, and distribution of MTBE, three other ether gasoline oxygenates, and 62 other volatile organic compounds (VOCs) in ground- and surface-water sources used for drinking-water supplies. The Random Survey was the first of two components of a national assessment of the quality of source water supplying CWSs sponsored by the American Water Works Association Research Foundation. A total of 954 CWSs were selected for VOC sampling from the population of nearly 47,000 active, self-supplied CWSs in all 50 States, Native American Lands, and Puerto Rico based on a statistical design that stratified on CWS size (population served), type of source water (ground and surface water), and geographic distribution (State). At a reporting level of 0.2 μ g/L, VOCs were detected in 27 percent of source-water samples collected from May 3, 1999 through October 23, 2000. Chloroform (in 13 percent of samples) was the most frequently detected of 42 VOCs present in the source-water samples, followed by MTBE. VOC concentrations were generally less than 10 μ g/L 95 percent of the 530 detections and 63 percent were less than 1.0 μ g/L. Concentrations of 1,1-dichloroethene, tetrachloroethene, trichloroethene, vinyl chloride, and total trihalomethanes (TTHMs), however, exceeded drinking-water regulations in eight samples. Detections of most VOCs were more frequent in surface-water sources than in ground-water sources, with gasoline compounds collectively and MTBE individually detected significantly more often in surface water. Use of personal and commercial motorized watercraft on surface-water bodies that are drinking-water sources is probably the reason for the elevated detections of gasoline contaminants relative to ground water. MTBE detections demonstrated a seasonal pattern with more frequent detections in surface water in summer months, which is consistent with seasonal watercraft use. The detection frequency of most VOCs was significantly related to urban land use and population density. Detections of any VOC, non-trihalo-methane compounds, gasoline compounds collectively, the specific gasoline compounds benzene, toluene, ethylbenzene, and xylenes (BTEX), MTBE, solvents, and refrigerants were significantly greater in areas with more than 60 percent urban land use and (or) population density greater than 1,000 people per square mile than in source waters from less urbanized or lower population-density areas. MTBE detections were five times more frequent in source waters from areas with high MTBE use than in source waters from low or no MTBE use, but, unlike other gasoline compounds, MTBE detections were not significantly related to the density of gasoline storage tanks near drinking-water sources.

Water-Resources Investigations Report↗

USGS2016

Remote estimation of surface pCO2 on the West Florida Shelf

Surface p CO 2 data from the West Florida Shelf (WFS) have been collected during 25 cruise surveys between 2003 and 2012. The data were scaled up using remote sensing measurements of surface water properties in order to provide a more nearly synoptic map of p CO 2 spatial distributions and describe their temporal variations. This investigation involved extensive tests of various model forms through parsimony and Principal Component Analysis, which led to the development of a multi-variable empirical surface p CO 2 model based on concurrent MODIS (Moderate Resolution Imaging Spectroradiometer) estimates of surface chlorophyll a concentrations (CHL, mg m −3 ), diffuse light attenuation at 490 nm (Kd_Lee, m −1 ), and sea surface temperature (SST, °C). Validation using an independent dataset showed a p CO 2 Root Mean Square Error (RMSE) of <12 µatm and a 0.88 coefficient of determination (R 2 ) for measured and model-predicted p CO 2 ranging from 300 to 550 µatm. The model was more sensitive to SST than to CHL and Kd_Lee, with a 1 °C change in SST leading to a ~16 µatm change in the predicted p CO 2 . Application of the model to the entire WFS MODIS time series between 2002 and 2014 showed clear seasonality, with maxima (~450 µatm) in summer and minima (~350 µatm) in winter. The seasonality was positively correlated to SST (high in summer and low in winter) and negatively correlated to CHL and Kd_Lee (high in winter and low in summer). Inter-annual variations of p CO 2 were consistent with inter-annual variations of SST, CHL, and Kd_Lee. These results suggest that surface water p CO 2 of the WFS can be estimated, with known uncertainties, from remote sensing. However, while the general approach of empirical regression may work for waters from other areas of the Gulf of Mexico, model coefficients need to be empirically determined in a similar fashion.

Continental Shelf Research↗

USGS2014

Pharmacokinetics of buprenorphine hydrochloride following intramuscular and intravenous administration to American kestrels ( Falco sparverius )

Objective —To determine the pharmacokinetics of buprenorphine hydrochloride after IM and IV administration to American kestrels ( Falco sparverius ). Animals —13 healthy 3-year-old captive-bred American kestrels. Procedures —Buprenorphine hydrochloride (0.6 mg/kg) was administered IM to all birds. Blood samples were collected at 9 times, ranging from 5 minutes to 9 hours after drug administration. Plasma buprenorphine concentrations were measured by use of tandem liquid chromatography–mass spectrometry. Pharmacokinetic parameters were determined by use of least squares linear regression and noncompartmental analysis of naïve pooled data. After a washout period of 2 weeks, the same dose of buprenorphine was administered IV to all birds and blood samples were collected at the same times after drug administration. Results —Maximum plasma buprenorphine concentration was achieved within 5 minutes after IM administration. For IM administration, bioavailability was 94.8% and elimination half-life was 92.1 minutes. For IV administration, steady-state volume of distribution was 4,023.8 mL/kg, plasma clearance was 49.2 mL/min/kg, and elimination half-life was 105.5 minutes. Conclusions and Clinical Relevance —Buprenorphine was rapidly absorbed, and bioavailability was good after IM administration to American kestrels. Plasma buprenorphine concentrations were > 1 ng/mL for 9 hours after both IM and IV administration. These results, in combination with those of a pharmacodynamic study, suggested that the analgesic effects of buprenorphine could last at least 6 to 9 hours in this species. Further investigations of the duration of analgesic effects, multiple-dose protocols, and potential adverse effects of buprenorphine are warranted in American kestrels and other raptors.

American Journal of Veterinary Research↗

USGS2013

Detection of environmental DNA of bigheaded carps in samples collected from selected locations in the St. Croix River and in the Mississippi River

The use of molecular methods, such as the detection of environmental deoxyribonucleic acid (eDNA), have become an increasingly popular tool in surveillance programs that monitor for the presence of invasive species in aquatic systems. One early application of these methods in aquatic systems was surveillance for DNA of Asian carps (specifically bighead carp Hypophthalmichthys nobilis and silver carp H. molitrix) in water samples taken from the Chicago Area Waterway System. The ability to identify DNA of a species in an environmental sample presents a potentially powerful tool because these sensitive analyses can presumably detect the presence of DNA in water even when the species is not abundant or are difficult to catch or monitor with traditional gear. Prior to research presented in this report, an initial eDNA surveillance effort was completed in selected locations in the Upper Mississippi and St. Croix Rivers in 2011 after the capture of a bighead carp in the St. Croix River near Prescott, WI. Data presented in this report were developed to duplicate the 2011 monitoring results from the Upper Mississippi and St. Croix Rivers and to provide critical insight into the technique to inform future work in these locations. We specifically sought to understand the potential confounding effects of other pathways of eDNA movement (e.g., fish-eating birds, watercraft) on the variation in background DNA by collecting water samples from (1) sites within the St. Croix River and the upper Mississippi River where the DNA of silver carp was previously detected, (2) sites considered to be free of Asian carp, and (3) a site known to have a large population of Asian carp. We also sought to establish a baseline Asian carp eDNA signature to which future eDNA sampling efforts could be compared. All samples taken as part of this effort were processed using conventional polymerase chain reaction (PCR) according to procedures outlined in the U.S. Army Corps of Engineers Quality Assurance Project Plan with minor deviations designed to enhance the rigor of our data. Presence of DNA in PCR-positive samples was confirmed by Sanger sequencing (forward and reverse) and sequences were considered positive only if sequences (forward and reverse) of ≥150 base pairs had a match of ≥95% to those of published sequences for bighead carp or silver carp. The DNA of bighead carp and silver carp was not detected in environmental samples collected above and below St. Croix Falls Dam on the St. Croix River, above and below the Coon Rapids Dam and below Lock and Dam 1 on the Upper Mississippi River, and from two negative control lakes, Square Lake and Lake Riley. The DNA of silver carp was detected in environmental samples collected below Lock and Dam 19 at Keokuk, Iowa, a reach of the river with high silver carp abundance. The portion (68%) of environmental samples taken below Lock and Dam 19 that were determined to contain the DNA of silver carp was similar to that reported in the scientific literature for other abundant species. The DNA of bighead carp, however, was not detected in environmental samples collected below Lock and Dam 19, a reach of the river known to have bighead carp. Previous reported detections of the DNA of silver carp in samples collected in 2011 were not replicated in this study. Additional analyses are planned for the DNA extracted from the samples collected in 2012. Those analyses may provide additional information regarding the lack of amplification of bighead carp DNA and the lengths of the sequences of silver carp DNA present in samples taken below Lock and Dam 19. These additional analyses may help inform the use of eDNA monitoring in large, complex systems like the Mississippi River.

Minnesota, Wisconsin↗

USGS1997

Spatial earthquake hazard assessment of Evansville, Indiana

The earthquake hazard has been evaluated for a 150-square-kilometer area around Evansville, Indiana. GIS-QUAKE, a system that combines liquefaction and ground motion analysis routines with site-specific geological, geotechnical, and seismological information, was used for the analysis. The hazard potential was determined by using 586 SPT borings, 27 CPT sounding, 39 shear-wave velocity profiles and synthesized acceleration records for body-wave magnitude 6.5 and 7.3 mid-continental earthquakes, occurring at distances of 50 km and 250 km, respectively. The results of the GIS-QUAKE hazard analyses for Evansville identify areas with a high hazard potential that had not previously been identified in earthquake zonation studies. The Pigeon Creek area specifically is identified as having significant potential for liquefaction-induced damage. Damage as a result of ground motion amplification is determined to be a moderate concern throughout the area. Differences in the findings of this zonation study and previous work are attributed to the size and range of the database, the hazard evaluation methodologies, and the geostatistical interpolation techniques used to estimate the hazard potential. Further, assumptions regarding the groundwater elevations made in previous studies are also considered to have had a significant effect on the results.

Indiana↗

USGS2016

Groundwater science relevant to the Great Lakes Water Quality Agreement: A status report

When the Great Lakes Water Quality Agreement (GLWQA) was signed in 1972 by the Governments of Canada and the United States (the “Parties”) (Environment Canada, 2013a), groundwater was not recognized as important to the water quality of the Lakes. At that time, groundwater and surface water were still considered as two separate systems, with almost no appreciation for their interaction. When the GLWQA was revised in 1978 (US Environmental Protection Agency (USEPA), 2012), groundwater contamination, such as that reported at legacy industrial sites such as those at Love Canal near the Niagara River, was squarely in the news. Consequently, the potential impacts of contaminated groundwater from such sites on Great Lakes water quality became a concern (Beck, 1979), and Annex 16 was added to the agreement, to address “pollution from contaminated groundwater” (Francis, 1989). However, no formal process for reporting under this annex was provided. The GLWQA Protocol in 1987 modified Annex 16 and called for progress reports beginning in 1988 (USEPA, 1988). The Protocol in 2012 provided a new Annex 8 to address groundwater more holistically (Environment 2 Canada, 2013b). Annex 8 (Environment Canada, 2013b) commits the Parties to coordinate groundwater science and management actions; as a first step, to “publish a report on the relevant and available groundwater science” by February 2015 (this report); and to “identify priorities for science activities and actions for groundwater management, protection, and remediation…” The broader mandate of Annex 8 is to (1) “identify groundwater impacts on the chemical, physical and biological integrity of the Waters of the Great Lakes;” (2) “analyze contaminants, including nutrients in groundwater, derived from both point and non-point sources impacting the Waters of the Great Lakes;” (3) “assess information gaps and science needs related to groundwater to protect the quality of the Waters of the Great Lakes;” and (4) “analyze other factors, such as climate change, that individually or cumulatively affect groundwater’s impact on the quality of the Waters of the Great Lakes.” A binational Annex 8 Subcommittee was formed to lead efforts to fulfill the mandate of this annex (members listed on p. i of this report). In turn, this subcommittee has recruited a task team to prepare this report (listed as authors of each chapter). This report addresses all of the above four objectives, based on a compilation of the “relevant and available groundwater science.” Specifically, the second objective (to “analyze contaminants”) is addressed by incorporating information obtained in ongoing monitoring and research activities conducted by the Parties, and by various other members of the Great Lakes Executive Committee.

Illinois, Indiana, Michigan, Minnesota, New York, ↗