Search USGS⌕ Search

SEARCH · Search USGS

Results for “Environmental Software”

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 271 records · Page 15Linked to original sources
USGS2005

Trends in surface-water quality at selected National Stream Quality Accounting Network (NASQAN) stations, in Michigan

To demonstrate the value of long-term, water-quality monitoring, the Michigan Department of Environmental Quality (MDEQ), in cooperation with the U.S. Geological Survey (USGS), initiated a study to evaluate potential trends in water-quality constituents for selected National Stream Quality Accounting Network (NASQAN) stations in Michigan. The goal of this study is to assist the MDEQ in evaluating the effectiveness of water-pollution control efforts and the identification of water-quality concerns. The study included a total of nine NASQAN stations in Michigan. Approximately 28 constituents were analyzed for trend tests. Station selection was based on data availability, land-use characteristics, and station priority for the MDEQ Water Chemistry Monitoring Project. Trend analyses were completed using the uncensored Seasonal Kendall Test in the computer program Estimate Trend (ESTREND), a software program for the detection of trends in water-quality data. The parameters chosen for the trend test had (1) at least a 5-year period of record (2) about 5 percent of the observations censored at a single reporting limit, and (3) 40 percent of the values within the beginning one-fifth and ending one-fifth of the selected period. In this study, a negative trend indicates a decrease in concentration of a particular constituent, which generally means an improvement in water quality; whereas a positive trend means an increase in concentration and possible degradation of water quality. The results of the study show an overall improvement in water quality at the Clinton River at Mount Clemens, Manistee River at Manistee, and Pigeon River near Caseville. The detected trend for these stations show decreases in concentrations of various constituents such as nitrogen compounds, conductance, sulfate, fecal coliform bacteria, and fecal streptococci bacteria. The negative trend may indicate an overall improvement in agricultural practices, municipal and industrial wastewater-treatment processes, and effective regulations. Phosphorus data for most of the study stations could not be analyzed because of the data limitations for trend tests. The only station with a significant negative trend in total phosphorus concentration is the Clinton River at Mount Clemens. However, scatter-plot analyses of phosphorus data indicate decreasing concentrations with time for most of the study stations. Positive trends in concentration of nitrogen compounds were detected at the Kalamazoo River near Saugatuck and Muskegon River near Bridgeton. Positive trends in both fecal coliform and total fecal coliform were detected at the Tahquamenon River near Paradise. Various different point and nonpoint sources could produce such positive trends, but most commonly the increase in concentrations of nitrogen compounds and fecal coliform bacteria are associated with agricultural practices and sewage-plant discharges. The constituent with the most numerous and geographically widespread significant trend is pH. The pH levels increased at six out of nine stations on all the major rivers in Michigan, with no negative trend at any station. The cause of pH increase is difficult to determine, as it could be related to a combination of anthropogenic activities and natural processes occurring simultaneously in the environment. Trends in concentration of major ions, such as calcium, sodium, magnesium, sulfate, fluoride, chloride, and potassium, were detected at eight out of nine stations. A negative trend was detected only in sulfate and fluoride concentrations; a positive trend was detected only in calcium concentration. The major ions with the most widespread significant trends are sodium and chloride; three positive and two negative trends were detected for sodium, and three negative and two positive trends were detected for chloride. The negative trends in chloride concentrations outnumbered the positive trends. This result indicates a slight improvement in surface-water quality because chloride as a point source in natural water comes from deicing salt, sewage effluents, industrial wastes, and oil fields. For other major ions, such as magnesium and potassium, both positive and negative trends were detected. These changes in trends indicate changes in surface-water quality caused by a variety of point and non-point sources throughout Michigan, as well as natural changes in the environment.

Michigan↗
USGS2000

Monitoring beach morphology changes using small-format aerial photography and digital softcopy photogrammetry

Current methods of monitoring beach morphology changes commonly involve the establishment of Global Positioning System profiles that are surveyed on a regular basis. Although this method produces precise measurements of coastal topography, it is costly in time and effort and may result in large data gaps between profiles. Much of our understanding of coastal dynamics is thus limited by profile spacing and survey frequency. Softcopy photogrammetry is increasingly being used as an alternative to assess shoreline change and study beach morphodynamics. This method of producing three-dimensional topographic models and orthophotographs from digitized aerial photography can aid in filling in the gaps between established profiles. A limiting factor to this technology is the cost of obtaining high-resolution aerial photography. We have developed an aerial mapping system designed to collect data in an efficient and cost-effective way. We use a small-format aerial photography system that can be mounted on a variety of small aircraft on short notice. After a flight, the photographs are scanned, and softcopy photogrammetry software is used to create both DTMs (Digital Terrain Models) and orthophotographs. The DTMs can be compared with existing profiles for accuracy, and volumetric changes can be computed. The orthophotos are used to make precise measurements of the position and morphology of shoreline features. This aerial mapping system is advantageous over previous methods of beach morphology change monitoring because it allows for rapid response to storm events and provides a cost-effective method of establishing a continual monitoring program in erosion hazard areas.

Environmental Geosciences↗
USGS2006

Selected Streamflow Statistics for Streamgaging Stationsin Northeastern Maryland, 2006

Streamflow statistics were calculated for 47 U.S. Geological Survey (USGS) streamgaging stations in northeastern Maryland, in cooperation with (1) the University of Maryland, Baltimore County, Center for Urban Environmental Research and Education; (2) the Baltimore City Department of Public Works; and (3) the Baltimore County Department of Environmental Protection and Resource Management. The statistics include the mean, minimum, maximum, and standard deviation of the daily mean discharges for the periods of record at the stations, as well as flow-duration and low-flow frequency statistics. The flow-duration statistics include the 1-, 2-, 5-, 10-, 15-, 20-, 25-, 30-, 40-, 50-, 60-, 70-, 75-, 80-, 85-, 90-, 95-, 98-, and 99-percent duration discharges. The low-flow frequency statistics include the average discharges for 1, 7, 14, and 30 days that recur, on average, once in 1.01, 2, 5, 10, 20, 50, and 100 years. The statistics were computed only for the 25 stations with periods of record of 10 years or more. The statistics were computed from records available through September 30, 2004 using standard methods and computer software developed by the USGS. A comparison between low-flow frequency statistics computed for this study and for a previous study that used data available through September 30, 1989 was done for seven stations. The comparison indicated that, for the 7-day mean low flow, the newer values were 19.8 and 15.3 percent lower for the 20- and 10-year recurrence intervals, respectively, and 2.1 percent higher for the 2-year recurrence interval, than the older values. For the 14-day mean low flow, the newer 20- and 10-year values were 25.2 and 15.5 percent lower, respectively, and the 2-year value was 2.9 percent higher than the older values. For the 30-day mean low flow, the newer 20-, 10-, and 2-year values were 10.8, 7.9, and 0.8 percent lower, respectively, than the older values. The newer values are generally lower than the older ones most likely because two major droughts have occurred since the older study was completed.

Open-File Report↗
USGS1999

Improvements to the DRASTIC ground-water vulnerability mapping method

Ground-water vulnerability maps are designed to show areas of greatest potential for ground-water contamination on the basis of hydrogeologic and anthropogenic (human) factors. The maps are developed by using computer mapping hardware and software called a geographic information system (GIS) to combine data layers such as land use, soils, and depth to water. Usually, ground-water vulnerability is determined by assigning point ratings to the individual data layers and then adding the point ratings together when those layers are combined into a vulnerability map. Probably the most widely used ground-water vulnerability mapping method is DRASTIC, named for the seven factors considered in the method: Depth to water, net Recharge, Aquifer media, Soil media, Topography, Impact of vadose zone media, and hydraulic Conductivity of the aquifer (Aller and others, 1985, p. iv). The DRASTIC method has been used to develop ground-water vulnerability maps in many parts of the Nation; however, the effectiveness of the method has met with mixed success (Koterba and others, 1993, p. 513; U.S. Environmental Protection Agency, 1993; Barbash and Resek, 1996; Rupert, 1997). DRASTIC maps usually are not calibrated to measured contaminant concentrations. The DRASTIC ground-water vulnerability mapping method was improved by calibrating the point rating scheme to measured nitrite plus nitrate as nitrogen (NO2+NO3–N) concentrations in ground water on the basis of statistical correlations between NO2+NO3–N concentrations and land use, soils, and depth to water (Rupert, 1997). This report describes the calibration method developed by Rupert and summarizes the improvements in results of this method over those of the uncalibrated DRASTIC method applied by Rupert and others (1991) in the eastern Snake River Plain, Idaho.

Idaho;Nevada;Utah;Wyoming↗
USGS2002

Testing pop-up satellite tags as a tool for identifying critical habitat for Pacific halibut ( Hippoglossus stenolepis ) in the Gulf of Alaska

To maintain healthy commercial and sport fisheries for Pacific halibut ( Hippoglossus stenolepis ), critical habitat must be defined by determining life history patterns on a daily and seasonal basis. Pop-up satellite archival transmitting (PSAT) tags provide a fisheries-independent method of collecting environmental preference data (depth and ambient water temperature) as well as daily geolocation estimates based on ambient light conditions. In this study, 14 adult halibut (107-165 cm FL) were tagged and released with PSAT tags in and around Resurrection Bay, Alaska. Commercial fishermen recovered two tags, while five tags transmitted data to ARGOS satellites. Horizontal migration was not consistent among fish as three halibut remained in the vicinity of release while four traveled up to 358 km from the release site. Vertical migration was not consistent among fish and over time, but they spent most their time between 150-350 m. The minimum and maximum depths reached by any of the halibut were 2m and 502m, respectively. The fish preferred water temperatures of roughly 6 °C while experiencing ambient temperatures between 4.3 °C and 12.2 °C. Light attenuation with depth prevented existing geolocation software and light sensing hardware from accurately estimating geoposition, however, information from temperature, depth, ocean bathymetry, and pop-off locations provided inference on fish movement in the study area. PSAT tags were a viable tool for determining daily and seasonal behavior and identifying critical halibut habitat, which will aid fisheries managers in future decisions regarding commercial and sport fishing regulations.

Alaska↗
USGS2026

3D semantic mapping of surface geological features

Semantic mapping in 3D is fundamental to a wide range of geoscientific studies and applications, including geomorphology, hazard assessment, and environmental monitoring. However, automatically segmenting geological features from large-scale photogrammetric datasets remains a significant challenge. We present a methodology to address this gap. Using overlapping images collected over environments of interest, Structure-from-Motion (SfM) produces georeferenced point clouds and estimates camera poses. Existing large vision models, such as Segment Anything Model, segment objects in the images, generating pixel-segmentation associations. To produce pixel-point associations, we project the points back onto the camera image planes. As objects are independently segmented across multiple images with different perspectives, we develop a segmentation mosaicking algorithm to build probabilistic point-segmentation associations that combines the pixel-segmentation associations and pixel-point associations. Our methodology is validated using both synthetic data generated by Kubric and real-world UAV-SfM data. The implementation is designed to be compatible with existing SfM software, including Agisoft and OpenDroneMap, for photogrammetry mapping in geoscience studies. As a case study, we apply our method to the semantic mapping of precariously balanced rocks (PBRs), which provide upper-bound constraints on historical ground motion shaking intensity. To support object-level identification of PBRs, we additionally integrated Grounding DINO, enabling text-prompted segmentation of features of interest within UAV imagery. This case study demonstrates the effectiveness of our method in generating a 3D semantic map of PBRs, enabling spatial distribution of PBR fragility for earthquake hazard analysis.

Computers & Geosciences↗
USGS2015

User guide to Exploration and Graphics for RivEr Trends (EGRET) and dataRetrieval: R packages for hydrologic data

Evaluating long-term changes in river conditions (water quality and discharge) is an important use of hydrologic data. To carry out such evaluations, the hydrologist needs tools to facilitate several key steps in the process: acquiring the data records from a variety of sources, structuring it in ways that facilitate the analysis, processing the data with routines that extract information about changes that may be happening, and displaying findings with graphical techniques. A pair of tightly linked R packages, called dataRetrieval and EGRET (Exploration and Graphics for RivEr Trends), have been developed for carrying out each of these steps in an integrated manner. They are designed to easily accept data from three sources: U.S. Geological Survey hydrologic data, U.S. Environmental Protection Agency (EPA) STORET data, and user-supplied flat files. The dataRetrieval package not only serves as a “front end” to the EGRET package, it can also be used to easily download many types of hydrologic data and organize it in ways that facilitate many other hydrologic applications. The EGRET package has components oriented towards the description of long-term changes in streamflow statistics (high flow, average flow, and low flow) as well as changes in water quality. For the water-quality analysis, it uses Weighted Regressions on Time, Discharge and Season (WRTDS) to describe long-term trends in both concentration and flux. EGRET also creates a wide range of graphical presentations of the water-quality data and of the WRTDS results. This report serves as a user guide to these two R packages, providing detailed guidance on installation and use of the software, documentation of the analysis methods used, as well as guidance on some of the kinds of questions and approaches that the software can facilitate.

Techniques and Methods↗
USGS2021

Avoidance of cold-, cool-, and warm-water fishes to Zequanox® exposure

Zequanox® is a biopesticide registered by the U.S. Environmental Protection Agency (USEPA) and the Canadian Pest Management Regulatory Agency for controlling dreissenid mussels with demonstrated selective toxicity. However, some research has indicated that Zequanox may impact the body condition and survival of some non-target species. We assessed avoidance behaviors of two species of cold-, cool-, and warm-water fishes to Zequanox at the maximum concentration allowed by the USEPA label (100 mg/L as active ingredient). Naïve, juvenile fish (n = 30 per species) were individually observed in a two-flume choice tank through which Zequanox-treated and untreated water simultaneously flowed in an unobstructed arena. Individual fish were observed during an untreated control period (20 min) and two Zequanox-exposure periods (20 min each). Treatment was alternated between arena sides to account for potential side bias in the test subjects. Positional data were collected and tabulated in real time with EthoVision® XT software. Zequanox concentrations and water quality properties (pH, dissolved oxygen, temperature, and specific conductance) were monitored during each trial. Analysis of treatment response was performed using a contrast within linear mixed-effects models. Our results indicate that Brook Trout, Lake Trout, and Bluegill avoided Zequanox-treated water, Yellow Perch were indifferent to Zequanox-treated water, and Lake Sturgeon and Fathead Minnow were attracted to Zequanox-treated water. These results combined with existing species sensitivity literature may help inform resource managers of potential treatment-related risks.

Management of Biological Invasions↗
USGS2022

Areas contributing recharge to priority wells in valley-fill aquifers in the Neversink River and Rondout Creek drainage basins, New York

In southeastern New York, the villages of Ellenville, Wurtsboro, Woodridge, the hamlet of Mountain Dale, and surrounding communities in the Neversink River and Rondout Creek drainage basins rely on wells that pump groundwater from valley-fill glacial aquifers for public water supply. Glacial aquifers are vulnerable to contamination because they are highly permeable and have a shallow depth to water table. To protect the quality of these water resources, water managers need accurate information about the areas that contribute recharge to production wells that pump from these aquifers. The New York State Department of Environmental Conservation and the New York State Department of Health designated eight priority wells in this region for which water supply protection is of primary concern. The U.S. Geological Survey, in cooperation with the New York State Department of Environmental Conservation and the New York State Department of Health, began an investigation in 2019 with the general objectives of (1) improving understanding of regional groundwater-flow system, (2) delineating areas contributing recharge to eight priority production wells, and (3) quantifying the uncertainty of these contributing areas in a probabilistic way that can be used to inform decision-making related to priority well source-water protection. To complete these objectives, a MODFLOW 6 groundwater model was created encompassing the eight priority wells and the surrounding flow system, which includes parts of the Neversink River and Rondout Creek Basins in Sullivan County and Ulster County, New York. The model was built using Python tools (such as flopy , modflow-setup , and sfrmaker ) that facilitate transparent and repeatable model development using existing datasets. The model parameters were estimated with a stepwise approach using an iterative ensemble smoother implementation of the Parameter ESTimation software PEST++ (version 5.0.0). We evaluated initial “best guess” parameter bounds with a prior Monte Carlo analysis. Results of the first prior Monte Carlo analysis were used to make informed adjustments to model parameter bounds (typically resulting in expanded bounds), and a second prior Monte Carlo analysis was run to identify improved ranges for model parameters during history matching. The history matching effort produced an ensemble of parameter values for the groundwater-flow model that spans the range of values within prior uncertainty bounds. The ensemble is informed by the historical observation data, within a reasonable range of uncertainty on those observations. This history-matched ensemble was used in a particle tracking Monte Carlo analysis to delineate the areas contributing recharge to priority wells. The groundwater-flow and particle tracking (MODPATH7) models were run once for each ensemble member. Deterministic contributing areas computed for each ensemble member were aggregated to produce maps showing the probability that a location contributes recharge to priority wells. Finally, the particle tracking Monte Carlo analysis was repeated for six pumping scenarios, representing a wide range of possible pumping levels, to incorporate uncertainty in future pumping rates related to population growth or other management decisions. Increasing pumping rates generally led to larger contributing recharge areas and larger areas of high probability that a location contributes recharge to priority wells. These maps show the overall uncertainty of the areas contributing recharge to priority wells in the study area and provide a tool for risk-based decision making for protection of well source water.

New York↗
USGS2020

Statistical methods in water resources

This text began as a collection of class notes for a course on applied statistical methods for hydrologists taught at the U.S. Geological Survey (USGS) National Training Center. Course material was formalized and organized into a textbook, first published in 1992 by Elsevier as part of their Studies in Environmental Science series. In 2002, the work was made available online as a USGS report. The text has now been updated as a USGS Techniques and Methods Report. It is intended to be a text in applied statistics for hydrology, environmental science, environmental engineering, geology, or biology that addresses distinctive features of environmental data. For example, water resources data tend to have many variables with a lower bound of zero, tend to be more skewed than data from many other disciplines, commonly contain censored data (less than values), and assumptions that the data are normally distributed are not appropriate. Computer-intensive methods (bootstrapping and permutation tests) now improve upon and replace the dependence on t-intervals, t-tests, and analysis of variance. A new chapter on sampling design addresses questions such as “How many observations do I need?” The chapter also presents distribution-free methods to help plan sampling efforts. The trends chapter has been updated to include the WRTDS (Weighted Regressions on Time, Discharge, and Season) method for analysis of water-quality data. This new version contains updated graphics and updated guidance on the use of statistical techniques. The text utilizes R, a programming language and open-source software environment, for all exercises and most graphics, and the R code used to generate figures and examples is provided for download.

Techniques and Methods↗
USGS2016

Users and uses of Landsat 8 satellite imagery—2014 survey results

Executive Summary In 2013, Landsat 8 began adding high quality, global, moderate-resolution imagery to the more than 40-year archive of Landsat imagery. To assess the potential effects of the availability of Landsat 8 imagery on users and their work, the U.S. Geological Survey (USGS) Land Remote Sensing Program (LRS) initiated a survey of Landsat users. The objectives of the survey were to 1. Characterize various Landsat user groups, such as United States (U.S.) and international users and Landsat 8 and non-Landsat 8 users; 2. Identify any differences among user groups in uses and preferences; 3. Measure the importance of and satisfaction with Landsat 8 attributes; 4. Assess the importance to users of the frequency of usable imagery; and 5. Determine any challenges in using Landsat 8. The online survey was sent to 51,617 Landsat users registered with USGS in May 2014. Almost 13,000 people responded to the survey for a response rate of 25 percent (n = 12,966). Current Landsat users (users who had used Landsat in their work in the year prior to the survey) composed 89 percent of the sample (n = 11,549) and past Landsat users composed 11 percent (n = 1,417). The results reported here apply to current Landsat users registered with the USGS Earth Resources Observation and Science (EROS) Center. Users from 161 countries responded to the survey. Of those, 19 percent were citizens or permanent residents of the United States and 81 percent resided in other countries. More than 70 percent of current users had used Landsat 8 in the year prior to the survey. The majority of Landsat 8 users (65 percent) were established users who used Landsat imagery regularly both before and after Landsat 8 imagery became available. The average current Landsat user was male, 36 years old, and highly educated, with 9 years of experience using satellite imagery or geographic information system (GIS) software. Landsat 8 users had, on average, two more years of experience than non-Landsat 8 users. Users were employed predominantly by academic institutions (65 percent), followed by private businesses (13 percent), Federal governments (10 percent), State and local governments (6 percent), and nonprofit organizations (6 percent). Of the Landsat imagery obtained in the past year by current users, on average 31 percent came from a Landsat 8 sensor. An equivalent amount came from the Landsat 7 ETM+ sensor (33 percent); slightly less came from Landsats 4 and 5 TM sensors (27 percent). Much less came from Landsats 1 through 5 MSS sensors (5 percent). Overall, more than a third of users’ work used Landsat imagery (38 percent). Of this work, on average, 37 percent of the work was operational. Landsat 8 users considered a greater proportion of their work operational than non-Landsat 8 users (39 percent compared with 29 percent). Environmental sciences and management were the most commonly selected primary applications (selected by 42 percent of users). Land use/land cover (23 percent) was the second most commonly selected primary application, followed by education (12 percent), agriculture (9 percent), and planning and development (6 percent). Landsat 8 users were asked to rank the importance of certain attributes in determining whether to use Landsat 8 imagery in their work. The archive was ranked most important, followed by cost, spatial resolution, extent of coverage, data quality, and frequency of revisit. Users were asked how satisfied they were with these same attributes as they currently apply to Landsat 8 imagery. On average, users were most satisfied with lack of cost, extent of coverage, data quality, and the archive, but they were satisfied with all attributes. Users were asked how often they needed Landsat imagery to meet various requirements for their primary application. The survey question specifically asked how often users needed usable imagery, which differs from how often they would like the Landsat satellites to acquire an image. Users were asked to identify their needed frequency of usable imagery for the following levels: 1. Threshold level—the minimum frequency of usable imagery needed to be of any value to their primary application. 2. Breakthrough level—the frequency of usable imagery that would result in a significant improvement for their primary application of the imagery. 3. Target level—the frequency of usable imagery that would only provide a limited additional increase in the expected performance for their primary application. To meet the threshold level, three-quarters of users needed usable imagery every 17 days or less frequently. At the breakthrough level, two-thirds of users (64 percent) needed a usable image every 5–16 days. The current constellation of two satellites (Landsat 7 and 8) is capable of meeting the threshold and breakthrough needs of most users at least some of the time, but a single satellite would be highly unlikely to do so. Two-fifths of users (40 percent) felt that usable imagery provided every 4 days or more frequently would meet their target level which the current Landsat constellation cannot provide. Landsat 8 users were significantly more likely than non-Landsat 8 users to need usable imagery more frequently to meet their target levels. Additionally, U.S. Landsat 8 users were significantly more likely than other Landsat users to need usable imagery more frequently in order meet both their breakthrough and target levels. To explore the effect of the availability of Landsat 8 imagery on Landsat imagery use in general, established users (those who had consistently used Landsat imagery both before and after Landsat 8 imagery became available) using Landsat 8 imagery were asked about changes in the amount of Landsat imagery they used. The majority of established users using Landsat 8 imagery (60 percent) reported an average increase of 51 percent in the number of scenes obtained after Landsat 8 imagery became available. Landsat 8 users were asked if they had encountered challenges in using Landsat 8 whereas non-Landsat 8 users were asked if such challenges had played a role in why they were not using Landsat 8 imagery. Although many users did not encounter challenges when using or trying to use Landsat 8 data, slightly less than 30 percent did encounter issues with processing the data to a usable point. The most common issue reported was not being able to create or have access to a surface reflectance corrected product. Other challenges were related to the file sizes of images being too large to download, store, or analyze. There were no statistically significant differences between Landsat 8 and non-Landsat 8 users in terms of challenges encountered when using or trying to use the imagery, which indicates that users were not unduly discouraged by the challenges they may have encountered. When asked about potential consequences of not using Landsat 8, more than half of the non-Landsat 8 users did not report detrimental effects on their work from not using the imagery. Of those who did report detrimental effects, decreased quality of work, decreased scope of work, and increased time spent on work were the most common.

Open-File Report↗
USGS2018

Revised groundwater-flow model of the glacial aquifer system north of Aberdeen, South Dakota, through water year 2015

The city of Aberdeen, in northeastern South Dakota, requires an expanded and sustainable supply of water to meet current and future demands. Conceptual and numerical models of the glacial aquifer system in the area north of Aberdeen were developed by the U.S. Geological Survey in cooperation with the City of Aberdeen in 2012. The U.S. Geological Survey, in cooperation with the City of Aberdeen, completed a study to revise the original numerical groundwater-flow model using data through water year (WY) 2015 to aid the City of Aberdeen in their development of plans and strategies for a sustainable water supply and to increase understanding of the glacial aquifer system and groundwater-flow system near Aberdeen. The original model was revised to improve the fit between model-simulated values and observed (measured or estimated) data, provide greater insight into surface-water interactions, and improve the usefulness of the model for water-supply planning. The revised groundwater-flow model (hereafter referred to as the “revised model”) presented in this report supersedes the original model. The purpose of this report is to describe a revised groundwater-flow model including data collection, model calibration, and model results for the glacial aquifer system including the Elm, Middle James, and Deep James aquifers north of Aberdeen, South Dakota, using updated hydrologic data through WY 2015. The original numerical model was revised in several ways. The model was modified by adding four new layers, which included a surficial layer, two intervening confining layers, and a shale bedrock layer. The revised model provides an improved understanding of the groundwater-flow system in comparison to the original model. The principal aquifers of the model area include portions of the Elm, Middle James, and Deep James aquifers. The lithologic information used to define and describe the aquifers in the model area was unaltered; however, aquifer properties and boundary conditions were reviewed and updated using geological information reported by the South Dakota Department of Environmental and Natural Resources and information obtained from geophysical investigations for this study. The horizontal extent of the Elm, Middle James, and Deep James aquifers was unaltered from the original model. The thickness of the Deep James aquifer was modified based on interpretations from the geophysical investigations. In general, groundwater in the Elm aquifer flowed from northwest to southeast and locally towards rivers and streams. Similarly, in the Middle James and Deep James aquifers, groundwater also typically flowed southeast. The revisions made to the original model include use of the following MODFLOW stress packages: Recharge, Evapotranspiration, Time-Variant Specified Head, Wells, Drains, and Stream Flow Routing, all of which were updated from the original model except for the Stream Flow Routing Package, which replaced the River Package used in the original model. Model calibration is the process of estimating model parameters to minimize the differences, or residuals, between observed data and simulated values; therefore, Parameter ESTimation (PEST) software was used to optimize model input parameters by matching model-simulated values to observed data. Calibration parameters included horizontal hydraulic conductivity, vertical hydraulic conductivity, specific yield, specific storage, and vertical streambed conductance for stream and drain cells. Multipliers were used to calibrate the recharge and evapotranspiration stresses. Evapotranspiration extinction depth also was adjusted during model calibration. Comparisons to the original model are described to highlight the changes made in the revised model. In general, the revised model adequately simulates the natural system and compares favorably with observed hydrologic data. Simulated water levels were evaluated by comparing them to single water-level observations at selected well locations. The selected wells were the same wells used in the original model. The coefficient of determination value between simulated and observed water levels for the revised model was 0.89 and included simulated and observed values from October 1, 1974 (WY 1975), through September 30, 2015 (WY 2015). The coefficient of determination value for the original model was 0.94 and included simulated and observed values from October 1, 1974, through September 30, 2009. The difference may indicate that the original model could have been overfit to hydraulic head observations because base flow was not simulated. The additional data used in the revised model included some climatically wetter, more extreme periods, such as 2011, in which annual precipitation was 30.9 inches. Average annual precipitation for the original model timeframe, which included data from WYs 1975–2009, was 20.26 inches. Additional precipitation data for WYs 2010–15, included in the revised model timeframe, resulted in an average annual precipitation for WYs 1975–2015 in the model area of 20.6 inches. The larger variability in climate data coupled with the additional water-level data could explain the lower coefficient of determination for water levels in the revised model. The revised model was used to calculate various groundwater-budget components for steady-state and transient conditions for WYs 1975–2015. The time-variant specified-head cells in the revised model had the largest change when compared to the original steady-state model for inflows and outflows. Comparing the transient budget components between the original and the revised models indicated that inflow from recharge and time-variant specified-head cells had the greatest effect on groundwater inflows, and outflow from storage had the greatest effect on groundwater outflows. The simulated potentiometric contours from the revised model were compared with (1) the observed (interpreted) potentiometric surface (layer 2) and the hydraulic head values (layers 4 and 6) and (2) the simulated contours from the original model. The simulated hydraulic gradients and general direction of groundwater flow in the Elm aquifer in the revised model generally matched the observed potentiometric contours, the simulated potentiometric contours from the original model, and general flow directions interpreted to be perpendicular to the contours. Minor discrepancies between simulated potentiometric contours from the revised model and the observed potentiometric contours may be due to the lack of observed data in the model area. The revised model was designed to reduce the limitations of the original model. The revisions were validated by comparing the results of the original model with the revised model. A primary benefit of the revised model is the inclusion of the surficial deposits and the confining units as explicit layers in the model. The addition of the surficial layer was beneficial for three primary reasons: (1) more accurate representation of recharge from precipitation, (2) more accurate representation of groundwater evapotranspiration, and (3) more accurate representation of groundwater and surface-water interactions. The groundwater model is a numeric approximation of a complex physical hydrologic system, and the revised model data were interpolated in regions with sparse data. Additionally, model discretization included averaged and interpolated values for water use, withdrawal rates, and hydraulic conductivity. The revised model provides a useful estimate for hydraulic gradients, groundwater-flow directions, and aquifer response to groundwater withdrawals.

South Dakota↗
USGS2022

A Central Asia hydrologic monitoring dataset for food and water security applications in Afghanistan

From the Hindu Kush mountains to the Registan Desert, Afghanistan is a diverse landscape where droughts, floods, conflict, and economic market accessibility pose challenges for agricultural livelihoods and food security. The ability to remotely monitor environmental conditions is critical to support decision making for humanitarian assistance. The Famine Early Warning Systems Network (FEWS NET) Land Data Assimilation System (FLDAS) global and Central Asia data streams provide information on hydrologic states for routine integrated food security analysis. While developed for a specific project, these data are publicly available and useful for other applications that require hydrologic estimates of the water and energy balance. These two data streams are unique because of their suitability for routine monitoring, as well as for being a historical record for computing relative indicators of water availability. The global stream is available at ∼ 1-month latency, and monthly average outputs are on a 10 km grid from 1982–present. The second data stream, Central Asia (21–56 ∘ N, 30–100 ∘ E), at ∼ 1 d latency, provides daily average outputs on a 1 km grid from 2000–present. This paper describes the configuration of the two FLDAS data streams, background on the software modeling framework, selected meteorological inputs and parameters, and results from previous evaluation studies. We also provide additional analysis of precipitation and snow cover over Afghanistan. We conclude with an example of how these data are used in integrated food security analysis. For use in new and innovative studies that will improve understanding of this region, these data are hosted by U.S. Geological Survey data portals and the National Aeronautics and Space Administration (NASA). The Central Asia data described in this paper can be accessed via the NASA repository at https://doi.org/10.5067/VQ4CD3Y9YC0R (Jacob and Slinski, 2021), and the global data described in this paper can be accessed via the NASA repository at https://doi.org/10.5067/5NHC22T9375G (McNally, 2018).

Earth System Science Data↗
USGS2018

Methods used to estimate daily streamflow and water availability in the Massachusetts Sustainable-Yield Estimator version 2.0

The Massachusetts Sustainable-Yield Estimator is a decision support tool that provides estimates of daily unaltered streamflow, water-use-adjusted streamflow, and water availability for ungaged, user-defined basins in Massachusetts. Daily streamflow at the ungaged site is estimated for unaltered (no water use) and water-use scenarios. The procedure for estimating streamflow was developed previously and has been implemented with minor changes and updated water-use data in version 2.0 of the Massachusetts Sustainable-Yield Estimator. Unaltered streamflow at selected exceedance probabilities is estimated by previously published regression equations. Streamflow is interpolated between the regressed quantiles to produce a continuous flow duration curve. A daily streamflow time series is produced for the ungaged site by relating the estimated flow duration curve at the ungaged site to a flow duration curve at a gaged reference site and then transferring the dates from the reference site to the ungaged site. Minor refinements were made to the previously published methods to estimate unaltered and water-use-adjusted streamflow, including a procedure to enforce the monotonic structure of the regression-based unaltered flow quantiles, improvements to the interpolation method used for computing the estimated flow duration curve, and updates to the methods used to compute time-lagged stream alterations from groundwater pumping or discharges. Additionally, a procedure was developed to estimate prediction intervals for daily and monthly unregulated streamflow time series at an ungaged site. The Massachusetts Sustainable-Yield Estimator computes water-use-adjusted streamflow using water-use data provided by the Massachusetts Department of Environmental Protection. Available water-use data included monthly withdrawal and wastewater discharge volumes from 2010 to 2014 for surface-water and groundwater sources. Water-use-adjusted streamflow represents the potential effect of current water use on natural streamflow in the basin over the range of historical hydrologic conditions. Georeferenced water withdrawal and discharge volumes were incorporated into the Massachusetts StreamStats web application for use in version 2.0 of the Massachusetts Sustainable-Yield Estimator. To compute water-use-adjusted streamflow, mean daily withdrawals and discharges within a user-defined basin are subtracted and added to the unaltered time series, respectively. Surface-water volumes are applied directly to the equation. Time-lagged streamflow alterations from groundwater withdrawal or wastewater discharge sources are estimated by using a response-coefficient method developed from results of previously published, calibrated groundwater models in Massachusetts. The Massachusetts Sustainable-Yield Estimator was updated to version 2.0 to improve software stability and usability. The version 2.0 software application was developed in Microsoft Access with a graphical user interface. All geoprocessing steps, including basin delineation and compilation of basin characteristics and water use within the basin, were completed in the Massachusetts StreamStats web application and exported for use by the Massachusetts Sustainable-Yield Estimator version 2.0.

Massachusetts↗
USGS2021

Capacity assessment for Earth Monitoring, Analysis, and Prediction (EarthMAP) and future integrated monitoring and predictive science at the U.S. Geological Survey

Executive Summary Managers of our Nation’s resources face unprecedented challenges driven by the convergence of increasing, competing societal demands and a changing climate that affects the stability, vulnerability, and predictability of those resources. To help meet these challenges, the scientific community must take advantage of all available technologies, data, and integrative Earth systems modeling capacity to better inform resource and risk management decisions. This is the overarching goal of the U.S. Geological Survey (USGS) Earth Monitoring, Analysis, and Prediction (EarthMAP) vision: “By 2030, the USGS will deliver well integrated observations and predictions of the future state of natural systems—water, ecosystems, energy, minerals, hazards—at regional and national scales, working primarily with federal, state, and academic partners to develop and operate the capability” (U.S. Geological Survey, 2021). Providing more integrated Earth systems science and actionable information to decision makers, stakeholders, and the public requires a better understanding of the depth and distribution of existing capacity (capabilities, tools, and techniques) across the Bureau. Identifying existing capacity is also a critical first step toward gap analysis and targeted investments to increase capacity over time. The USGS formed a Capacity Assessment Team (CAT) and charged it with (1) conducting a Request for Information (RFI) to identify existing USGS expertise and activities supportive of integrated and predictive science to inform decision making, (2) developing a strategy and proof-of-concept for a continuously updated capacity assessment capability, and (3) identifying lessons learned to inform development of best practices for future capacity assessment efforts. The RFI took the form of a survey, with content guided by the science and technology needs identified in a USGS report titled “Grand Challenges for Integrated U.S. Geological Survey Science—A Workshop Report” (Jenni and others, 2017). The 44-question survey provided respondents the ability to rate their level of experience with a suite of priority disciplines, analysis and modeling approaches, technologies, and stakeholder engagement strategies and to enter optional narrative text for supporting context. An introductory portion focused on general science capacity assessment, followed by three sections targeting capabilities related to the foundational components of EarthMAP: (1) data and information integration, (2) integrated predictive science, and (3) actionable information. The survey results provided a high-level snapshot of USGS capacity in the targeted areas. Respondents (1,035 individuals) represented approximately 13 percent of the USGS across all mission areas and regions. Seventy-four percent of the respondents held a science-focused position title and the remainder had position titles in information technology, computer science, management, administrative, or other (contractors, volunteers, emeritus, and unknown). To provide greater insight into respondent capabilities and activities, information from the U.S. Department of the Interior and USGS enterprise information systems were used to further characterize topical expertise and organizational associations of survey respondents. To address the ongoing need to assess the Bureau’s capacity to address integrated predictive science priorities, the CAT developed a software-based proof-of-concept called the Integrated Science Assessment Information Database (iSAID) for assembling various information sources together toward making the full extent of USGS capabilities and scientific assets available for routine capacity assessment. This proof-of-concept is intended to serve as a catalyst for further development. The process of implementing the EarthMAP capacity assessment survey, analyzing survey responses, and developing the proof-of-concept resulted in lessons learned, findings, and recommendations. Example scenarios throughout the report demonstrate how capacity assessment data can inform science planning. Three overarching findings and recommendations are: (1) Finding: Capacity is limited in some critical disciplines, skills, and technology applications, but “sufficient” depends on the question and the need relative to availability at a given point in time. Recommendation: Develop an on-demand capacity assessment framework that enables rapid identification and evaluation of existing and available expertise to support decision needs as they arise. (2) Finding: Institutional barriers and lack of awareness constrain the ability of USGS staff to adopt new technologies, collaborate across administrative boundaries, and deliver actionable information to stakeholders in a timely manner. However, these barriers are not universally experienced. Recommendation: Pursue more targeted inquiries to clarify which institutional barriers are obstructing the adoption of new technologies and approaches or the sharing of expertise and equipment across organizational and regional boundaries. These inquiries should inform USGS leadership, mission areas, and regions whether policies can be revised or whether a lack of understanding is creating perceived obstacles. Highlight cases when staff have successfully adopted new technologies and approaches to advance EarthMAP priorities and provide actionable information in a timely manner to spread awareness of how perceived obstacles can be navigated and overcome when appropriate. (3) Finding: Examples of people and projects integrating across disciplines and scales and applying advanced approaches to meet complex stakeholder needs exist. Such examples provide transfer value across the spectrum from approach to decision making. Many projects, already underway, appear to meet elements of the EarthMAP vision, and the USGS has people who can provide leadership in multiple types of specific integrated science efforts. Recommendation: Use these findings as a starting point for near-term strategic planning for integrated science. Highlight, incentivize, and build on existing interdisciplinary predictive science and information delivery activities across the USGS to advance toward further realization of an EarthMAP capacity. The CAT efforts to develop and assess existing USGS capacity to advance the EarthMAP vision revealed a fundamental challenge for not only this effort but any effort to assess existing capacity: A considerable amount of thought, time, and effort is required to survey and assess capabilities and tools available to support a given need, yet best results are still likely to provide an incomplete assessment. To better meet the frequent need to assess capabilities, tools, products, and projects that address an expressed strategic priority, the CAT proposes the concept of an on-demand capacity assessment framework supported by a software package that dynamically pulls and integrates information from existing USGS information systems and public domain registries. Although existing USGS enterprise information systems currently lack the structure, cross-system consistency, interoperability, and stability to support a continuously updated capacity assessment capability, we identify reasonable near-term steps to improve the utility of information gathered on expertise and project capacity and to improve the consistency and completeness of information and the ability of USGS systems to share that information. The ability to search and characterize this information will make future assessments of capacity faster, more complete, more efficient, and more targeted. This approach would grow the Bureau’s capacity knowledge over time, iteratively improving the ability to access, leverage, and synthesize existing capabilities and assets as well as identify and fill critical gaps. The greatest promise for developing integrated science could lie in linking across existing projects and expertise to create a multi-project capacity for addressing large, complex environmental issues.

Open-File Report↗
USGS2023

SaTSeaD: Satellite Triangulated Sea Depth open-source bathymetry module for NASA Ames Stereo Pipeline

We developed the first-ever bathymetric module for the NASA Ames Stereo Pipeline (ASP) open-source topographic software called Satellite Triangulated Sea Depth, or SaTSeaD, to derive nearshore bathymetry from stereo imagery. Correct bathymetry measurements depend on water surface elevation, and whereas previous methods considered the water surface horizontal, our bathymetric module accounts for the curvature of the Earth in the imagery. The process is semiautomatic, reliable, and repeatable, independent of any external bathymetry data eliminating user bias in selecting bathymetry calibration points, and it can generate a fully integrated and seamless topo-bathymetry digital elevation model (TBDEM) in the same coordinate system, comparable with the band-ratio method irrespective of the regression method used for the band-ratio algorithm. The ASP output can be improved by applying a camera bundle adjustment to minimize reprojection errors and by alignment to a more accurate topographic (above water) surface without any bathymetric input since the derived TBDEM is a rigid surface. These procedures can decrease bathymetry root mean square errors from 30 to 80 percent, depending on environmental conditions, the quality of satellite imagery, and the spectral band used (e.g., blue, green, or panchromatic).

Remote Sensing↗
USGS2018

An experimental comparison of composite and grab sampling of stream water for metagenetic analysis of environmental DNA

Use of environmental DNA (eDNA) to assess distributions of aquatic and semi-aquatic macroorganisms is promising, but sampling schemes may need to be tailored to specific objectives. Given the potentially high variance in aquatic eDNA among replicate grab samples, compositing smaller water volumes collected over a period of time may be more effective for some applications. In this study, we compared eDNA profiles from composite water samples aggregated over three hours with grab water samples. Both sampling patterns were performed with identical autosamplers paired at two different sites in a headwater stream environment, augmented with exogenous fish eDNA from an upstream rearing facility. Samples were filtered through 0.8 μm cellulose nitrate filters and DNA was extracted with a cetyl trimethylammonium bromide procedure. Eukaryotic and bacterial community profiles were derived by amplicon sequencing of 12S ribosomal, 16S ribosomal, and cytochrome oxidase I loci. Operational taxa were assigned to genus with a lowest common ancestor approach for eukaryotes and to family with the RDP Classifier software for prokaryotes. Eukaryotic community profiles were more consistent with composite sampling than grab sampling. Downstream, rarefaction curves suggested faster taxon accumulation for composite samples, and estimated richness was higher for composite samples as a set than for grab samples. Upstream, composite sampling produced lower estimated richness than grab samples, but with overlapping standard errors. Furthermore, a bimodal pattern of richness as a function of sequence counts suggested the impact of clumped particles on upstream samples. Bacterial profiles were insensitive to sample method, consistent with the more even dispersion expected for bacteria compared with eukaryotic eDNA. Overall, samples composited over 3 h performed equal to or better than triplicate grab sampling for quantitative community metrics, despite the higher total sequencing effort provided to grab replicates. On the other hand, taxon-specific detection rates did not differ appreciably and the two methods gave similar estimates of the ratio of the common fish genera Salmo and Coregonus at each site. Unexpectedly, Salmo eDNA dropped out substantially faster than Coregonus eDNA between the two sites regardless of sampling method, suggesting that differential settling affects the estimation of relative abundance. We identified bacterial patterns that were associated with eukaryotic diversity, suggesting potential roles as biomarkers of sample representativeness.

PeerJ↗
USGS2011

Moderate-resolution sea surface temperature data for the nearshore North Pacific

Coastal sea surface temperature (SST) is an important environmental characteristic in determining the suitability of habitat for nearshore marine and estuarine organisms. This publication describes and provides access to an easy-to-use coastal SST dataset for ecologists, biogeographers, oceanographers, and other scientists conducting research on nearshore marine habitats or processes. The data cover the Temperate Northern Pacific Ocean as defined by the 'Marine Ecosystems of the World' (MEOW) biogeographic schema developed by The Nature Conservancy. The spatial resolution of the SST data is 4-km grid cells within 20 km of the shore. The data span a 29-year period - from September 1981 to December 2009. These SST data were derived from Advanced Very High Resolution Radiometer (AVHRR) instrument measurements compiled into monthly means as part of the Pathfinder versions 5.0 and 5.1 (PFSST V50 and V51) Project. The processing methods used to transform the data from their native Hierarchical Data Format Scientific Data Set (HDF SDS) to georeferenced, spatial datasets capable of being read into geographic information systems (GIS) software are explained. In addition, links are provided to examples of scripts involved in the data processing steps. The scripts were written in the Python programming language, which is supported by ESRI's ArcGIS version 9 or later. The processed data files are also provided in text (.csv) and Access 2003 Database (.mdb) formats. All data except the raster files include attributes identifying realm, province, and ecoregion as defined by the MEOW classification schema.

Open-File Report↗