USGS Science⌕ Search

SEARCH · USGS Science

Results for “Scientific Data”

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 1,675 records · Page 93Linked to original sources
USGS2015

Simulation of groundwater flow and streamflow depletion in the Branch Brook, Merriland River, and parts of the Mousam River watersheds in southern Maine

Watersheds of three streams, the Mousam River, Branch Brook, and Merriland River in southeastern Maine were investigated from 2010 through 2013 under a cooperative project between the U.S. Geological Survey and the Maine Geological Survey. The Branch Brook watershed previously had been deemed “at risk” by the Maine Geological Survey because of the proportionally large water withdrawals compared to estimates of the in-stream flow requirements for habitat protection. The primary groundwater withdrawals in the study area include a water-supply well in the headwaters of the system and three water-supply wells in the coastal plain near the downstream end of the system. A steady-state groundwater flow model was used to understand the movement of water within the system, to evaluate the water budget and the effect of groundwater withdrawals on streamflows, and to understand streamflow depletion in relation to the State of Maine’s requirements to maintain in-stream flows for habitat protection. Delineation of the simulated groundwater divides compared to the surface-water divides suggests that the groundwater divides in the headwater areas do not exactly correspond to the surface-water divides. Under both pumping and non-pumping conditions, groundwater flows from the headwaters of the Branch Brook watershed into the Mousam River watershed. Pumping in the Mousam River watershed captures a small amount of groundwater from the Branch Brook basin. The cumulative effect of groundwater withdrawals on base flows in two rivers in the study area (Branch Brook and the Merriland River) was evaluated using the groundwater flow model. Streamflow depletion in the headwaters of Branch Brook was 0.12 cubic feet per second (ft 3 /s) for the steady-state simulation, or about 10 percent of the average base flow at that location. Downstream on Branch Brook, the total streamflow depletion from all the wells was 0.59 ft 3 /s, or 3 percent of the average base flow at that location. In the Merriland River downstream from the Merriland River well, the total amount of streamflow depletion was 0.6 ft 3 /s, or about 7 percent of the average base flow. The groundwater model was used to evaluate several different scenarios that could affect streamflow and groundwater discharging to the rivers and streams in the study area. The scenarios were (1) no pumping from the water-supply wells; (2) current pumping from the water-supply wells, but simulated drought conditions (25 percent reduction in recharge); (3) current recharge, but with increased pumping from the large water-supply wells; and (4) drought conditions and increased pumping combined. Simulations of increased pumping in the water-supply wells resulted in streamflow depletion in the headwaters of Branch Brook increasing to 16 percent of the headwater base flow. Simulated increases in the pumping in the coastal plain wells increased the amount of streamflow depletion to 6 percent of the flow in Branch Brook and to 8 percent of the flow in the Merriland River. The additional stress of a drought imposed on the model (25 percent less recharge) had a substantial impact on streamflows, as expected. If the simulated drought occurred simultaneously with an increase in pumping, the base flows would be reduced 48 percent in the headwaters of Branch Brook, compared to the no-pumping scenario. Downstream in Branch Brook, the total reduction in flow would be 29 percent of the simulated base flows in the no-pumping scenario, and in the Merriland River, the reduction would be 33 percent of the base flows in the no-pumping scenario. The study evaluated two different methods of calculating in-stream flow requirements for Branch Brook and the Merriland River—a set of statewide equations used to calculate monthly median flows and the MOVE.1 record-extension technique used on site-specific streamflow measurements. The August median in-stream flow requirement in the Merriland River was calculated as 7.18 ft 3 /s using the statewide equations but was 3.07 ft 3 /s using the MOVE.1 analysis. In Branch Brook, the August median in-stream flow requirements were calculated as 20.3 ft 3 /s using the statewide equations and 11.8 ft 3 /s using the MOVE.1 analysis. In each case, using site-specific data yields an estimate of in-stream flow that is much lower than an estimate the statewide equations provide.

Maine↗
USGS2016

Estimating flood magnitude and frequency at gaged and ungaged sites on streams in Alaska and conterminous basins in Canada, based on data through water year 2012

Estimates of the magnitude and frequency of floods are needed across Alaska for engineering design of transportation and water-conveyance structures, flood-insurance studies, flood-plain management, and other water-resource purposes. This report updates methods for estimating flood magnitude and frequency in Alaska and conterminous basins in Canada. Annual peak-flow data through water year 2012 were compiled from 387 streamgages on unregulated streams with at least 10 years of record. Flood-frequency estimates were computed for each streamgage using the Expected Moments Algorithm to fit a Pearson Type III distribution to the logarithms of annual peak flows. A multiple Grubbs-Beck test was used to identify potentially influential low floods in the time series of peak flows for censoring in the flood frequency analysis. For two new regional skew areas, flood-frequency estimates using station skew were computed for stations with at least 25 years of record for use in a Bayesian least-squares regression analysis to determine a regional skew value. The consideration of basin characteristics as explanatory variables for regional skew resulted in improvements in precision too small to warrant the additional model complexity, and a constant model was adopted. Regional Skew Area 1 in eastern-central Alaska had a regional skew of 0.54 and an average variance of prediction of 0.45, corresponding to an effective record length of 22 years. Regional Skew Area 2, encompassing coastal areas bordering the Gulf of Alaska, had a regional skew of 0.18 and an average variance of prediction of 0.12, corresponding to an effective record length of 59 years. Station flood-frequency estimates for study sites in regional skew areas were then recomputed using a weighted skew incorporating the station skew and regional skew. In a new regional skew exclusion area outside the regional skew areas, the density of long-record streamgages was too sparse for regional analysis and station skew was used for all estimates. Final station flood frequency estimates for all study streamgages are presented for the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities. Regional multiple-regression analysis was used to produce equations for estimating flood frequency statistics from explanatory basin characteristics. Basin characteristics, including physical and climatic variables, were updated for all study streamgages using a geographical information system and geospatial source data. Screening for similar-sized nested basins eliminated hydrologically redundant sites, and screening for eligibility for analysis of explanatory variables eliminated regulated peaks, outburst peaks, and sites with indeterminate basin characteristics. An ordinary least‑squares regression used flood-frequency statistics and basin characteristics for 341 streamgages (284 in Alaska and 57 in Canada) to determine the most suitable combination of basin characteristics for a flood-frequency regression model and to explore regional grouping of streamgages for explaining variability in flood-frequency statistics across the study area. The most suitable model for explaining flood frequency used drainage area and mean annual precipitation as explanatory variables for the entire study area as a region. Final regression equations for estimating the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probability discharge in Alaska and conterminous basins in Canada were developed using a generalized least-squares regression. The average standard error of prediction for the regression equations for the various annual exceedance probabilities ranged from 69 to 82 percent, and the pseudo-coefficient of determination (pseudo-R 2 ) ranged from 85 to 91 percent. The regional regression equations from this study were incorporated into the U.S. Geological Survey StreamStats program for a limited area of the State—the Cook Inlet Basin. StreamStats is a national web-based geographic information system application that facilitates retrieval of streamflow statistics and associated information. StreamStats retrieves published data for gaged sites and, for user-selected ungaged sites, delineates drainage areas from topographic and hydrographic data, computes basin characteristics, and computes flood frequency estimates using the regional regression equations.

Alaska↗
USGS2018

Status of groundwater-level altitudes and long-term groundwater-level changes in the Chicot, Evangeline, and Jasper aquifers, Houston-Galveston region, Texas, 2018

Since the early 1900s, most of the groundwater withdrawals in the Houston-Galveston region, Texas, have been from the three primary aquifers that compose the Gulf Coast aquifer system—the Chicot, Evangeline, and Jasper aquifers. Withdrawals from these aquifers are used for municipal supply, industrial, and irrigation purposes. This report, prepared by the U.S. Geological Survey in cooperation with the Harris-Galveston Subsidence District, City of Houston, Fort Bend Subsidence District, Lone Star Groundwater Conservation District, and Brazoria County Groundwater Conservation District, is one in an annual series of reports depicting the status of groundwater-level altitudes and long-term groundwater-level changes in the Chicot, Evangeline, and Jasper aquifers in the Houston-Galveston region. This report contains regional-scale maps depicting approximate 2018 groundwater-level altitudes (represented by measurements made during December 2017 through March 2018) and long-term groundwater-level changes for the Chicot, Evangeline, and Jasper aquifers. In 2018, groundwater-level-altitude contours for the Chicot aquifer ranged from 200 feet (ft) below the North American Vertical Datum of 1988 (hereinafter referred to as “datum”) to 200 ft above datum. The 1977–2018 groundwater-level-change contours for the Chicot aquifer depict a large area of decline in groundwater-level altitudes (120 ft) in northwestern Harris County. The largest rise in groundwater-level altitudes in the Chicot aquifer from 1977 to 2018 (180 ft) was in southeastern Harris County. Groundwater-level-altitude contours for the Evangeline aquifer ranged from 250 ft below datum to 200 ft above datum in 2018. The 1977–2018 groundwater-level-change contours for the Evangeline aquifer depict broad areas where groundwater-level altitudes either declined or rose. The largest decline in groundwater-level altitudes (320 ft) was in southern Montgomery County. The largest rise in groundwater-level altitudes in the Evangeline aquifer from 1977 to 2018 (240 ft) was in southeastern Harris County. In 2018, groundwater-level-altitude contours for the Jasper aquifer ranged from 200 ft below datum to 200 ft above datum. The 2000–18 groundwater-level-change contours for the Jasper aquifer depict groundwater-level declines throughout most of the study area where groundwater-level-altitude data from the Jasper aquifer were collected, with the largest decline (200 ft) in southern Montgomery County.

Texas↗
USGS2018

Modeling hydrodynamics, water temperature, and water quality in Klamath Straits Drain, Oregon and California, 2012–15

Executive Summary Located southwest of Klamath Falls, Oregon, Klamath Straits Drain is a 10.1-mile-long canal that conveys water uphill and northward through the use of pumps before discharging to the Klamath River. Klamath Straits Drain traverses an area that historically encompassed Lower Klamath Lake. Currently, the Drain receives water from farmland and from parts of the Lower Klamath Lake National Wildlife Refuge. To support water-quality improvement in Klamath Straits Drain, a hydrodynamic and water-temperature model was constructed and calibrated for calendar years 2012–15 with the two-dimensional model CE-QUAL-W2 (version 4.0). Water quality was calibrated for a subset of that time, from April 1, 2012 to March 31, 2015. Flows in calendar year 2012 were within the normal range, while calendar years 2013–15 were dry years. Significant findings from this study include: In the years studied, only limited flow entered Klamath Straits Drain at the upstream Headworks (KSDH) site. Most flow entered the Drain between KSDH and the E-EE pumps near Township Road through several irrigation channels and ditches. Few data were available to describe the quality of this water for the period of study. The E-EE and F-FF pumps along Klamath Straits Drain mainly operated automatically to keep water levels relatively steady. Ten-minute flow data at streamgage 11509340, downstream of the F-FF pumps, showed high-frequency on/off switching of the F-FF pumps. Combined with daily mean flow data from the F-FF pumps, the downstream 10-minute flow data allowed estimation of 10-minute pumping rates for the F-FF pumps. Paper pump charts showed the existence of short-term variability at the E-EE pumps; however, daily pump data were used at the E-EE pump location in the model. Water temperature in Klamath Straits Drain varied from less than 5 degrees Celsius (°C) (with occasional ice cover in December–January) to greater than 20 °C in May–September. In the years studied, specific conductance was typically 250–850 microsiemens per centimeter, higher than Klamath River specific conductance (typically 100–200 microsiemens per centimeter). Increased chlorophyll a in autumn and winter, along with supersaturated oxygen concentrations, indicated algal blooms in the Drain at that time of year. The blooms were most likely diatoms, based on the timing of blooms sampled elsewhere. Total nitrogen concentration was as much as 5.5 mg/L, with most in dissolved organic and particulate forms, and lower amounts in ammonia and nitrate+nitrite. Total phosphorus concentrations were distributed between orthophophorus (at a median concentration of 0.15 mg/L) and organic and particulate forms (at a median concentration of 0.13 mg/L). Most of the organic carbon in the Klamath Straits Drain was in dissolved rather than particulate form. Newly collected water-quality data for April 1, 2012–March 31, 2015 helped provide the impetus for this modeling study. However, a lack of some data still hindered the construction and calibration of this model. The model would benefit from additional data to describe water-quality boundary conditions, water-quality calibration data upstream of the F-FF pumps, short-term E-EE pump operations, and channel bathymetry in the reach between Highway 97 and the confluence with the Klamath River. Klamath River water mixed upstream into the Klamath Straits Drain, up to the Klamath Straits Drain F-FF pumps at Highway 97, when the F-FF pumps were not operating for periods of hours to days. The F-FF pumps were off for many days during this study, especially during dry years. The boundary between Klamath Straits Drain and the Klamath River was best modeled with an external head condition, which allows exchange of water between the river and the drain in both directions, upstream and downstream. Currently there is a flow gage, water-quality monitor, and a water-quality sampling site located downstream of the F-FF pumps, in the reach where Klamath Straits Drain water can mix with Klamath River water. To sample solely Klamath Straits Drain water, water samples would need to be collected only when the F-FF pumps are actively pumping. Alternately, the sampling location could be moved upstream of the pumps. Interpretation and use of historical water-quality data at the Klamath Straits Drain at Highway 97 site should be done in conjunction with information on pump activity to help inform whether mixing with Klamath River water may have occurred. Total 2014 (a dry year) phosphorus loads from the Drain to the Klamath River were lower and closer to total maximum daily load (TMDL) allocations, as compared to 2013, a year with greater flow and pumping. Modeled travel time through the Klamath Straits Drain, from Headworks to its confluence with the Klamath River, ranged from approximately 24 hours at high flow to 16 days or more, depending on how many days the pumps were turned off. The longer travel times are sufficient for important water-quality transformations, such as algal growth and organic-matter decomposition. This newly constructed model of the Klamath Straits Drain simulates flow, water levels, water temperature, and water quality with acceptable accuracy but with certain data limitations. This model should prove useful in evaluating potential strategies for flow and water-quality management and restoration.

California, Oregon↗
USGS2019

Flood-inundation maps for the Yellow River from River Drive to Centerville Highway, Gwinnett County, Georgia

Digital flood-inundation maps for a 16.4-mile reach of the Yellow River in Gwinnett County, Georgia, from 0.5 mile upstream from River Drive to Centerville Highway (Georgia State Route 124) were developed to depict estimates of the areal extent and depth of flooding corresponding to selected water levels (stages) at two U.S. Geological Survey (USGS) streamgages in the mapped area. The maps for the 9.0-mile reach from 0.5 mile upstream from River Drive to Stone Mountain Highway (U.S. Route 78) are referenced to the streamgage Yellow River near Snellville, Ga. (station 02206500), and the maps for the 7.4-mile reach from Stone Mountain Highway to Centerville Highway are referenced to the streamgage Yellow River at Ga. 124, near Lithonia, Ga. (02207120). Real-time stage information from these streamgages can be used with these maps to estimate near real-time areas of inundation. The forecasted peak-stage information for the USGS streamgages Yellow River near Snellville, Ga. (02206500), and Yellow River at Ga. 124, near Lithonia, Ga. (02207120), can be used in conjunction with the maps developed for this study to show predicted areas of flood inundation. A one-dimensional step-backwater model was developed using the U.S. Army Corps of Engineers Hydrologic Engineering Center's River Analysis System (HEC–RAS) software for the Yellow River and was used to compute flood profiles for a 16.4-mile reach of the Yellow River. The hydraulic model was then used to simulate 16 water-surface profiles at 1.0-foot (ft) intervals at the Yellow River near Snellville streamgage and 17 water-surface profiles at 1.0-ft intervals at the Yellow River near Lithonia streamgage. At the Yellow River near Snellville streamgage, the profiles ranged from a stage of 18.0 ft, which is 819.1 ft above the North American Vertical Datum of 1988 (NAVD 88), to a stage of 33.0 ft, which is 834.1 ft above NAVD 88. At the Yellow River near Lithonia streamgage, the profiles ranged from the National Weather Service action stage of 13.0 ft, which is 732.5 ft above NAVD 88, to a stage of 29.0 ft, which is 748.5 ft above NAVD 88. The simulated water-surface profiles were then combined with a geographic information system digital elevation model—derived from light detection and ranging (lidar) data having a 5.0-ft horizontal resolution—to delineate the area flooded at each 1.0-ft interval of stream stage for both streamgages.

Georgia↗
USGS2019

Status of groundwater-level altitudes and long-term groundwater-level changes in the Chicot, Evangeline, and Jasper aquifers, Houston-Galveston region, Texas, 2019

Since the early 1900s, most of the groundwater withdrawals in the Houston-Galveston region, Texas, have been from the three primary aquifers that compose the Gulf Coast aquifer system—the Chicot, Evangeline, and Jasper aquifers. Withdrawals from these aquifers are used for municipal supply, commercial and industrial use, and irrigation. This report, prepared by the U.S. Geological Survey in cooperation with the Harris-Galveston Subsidence District, City of Houston, Fort Bend Subsidence District, Lone Star Groundwater Conservation District, and Brazoria County Groundwater Conservation District, is one in an annual series of reports depicting the status of groundwater-level altitudes and long-term groundwater-level changes in the Chicot, Evangeline, and Jasper aquifers in the Houston-Galveston region. This report contains regional-scale maps depicting approximate 2019 groundwater-level altitudes (represented by measurements made during December 2018 through March 2019) and long-term groundwater-level changes in the Chicot, Evangeline, and Jasper aquifers. In 2019, groundwater-level-altitude contours for the Chicot aquifer ranged from 200 feet (ft) below the North American Vertical Datum of 1988 (hereinafter referred to as “datum”) to 200 ft above datum. The 1977–2019 groundwater-level-change contours for the Chicot aquifer depict a large area of decline in groundwater-level altitudes (100 ft) in northwestern Harris County. The largest rise in groundwater-level altitudes in the Chicot aquifer from 1977 to 2019 (200 ft) was in southeastern Harris County. In 2019, groundwater-level-altitude contours for the Evangeline aquifer ranged from 300 ft below datum to 200 ft above datum. The 1977–2019 groundwater-level-change contours for the Evangeline aquifer depict broad areas where groundwater-level altitudes either declined or rose. The largest decline in groundwater-level altitudes (280 ft) was in southern Montgomery and northern Harris Counties. The largest rise in groundwater-level altitudes in the Evangeline aquifer from 1977 to 2019 (240 ft) was in southeastern Harris County. In 2019, groundwater-level-altitude contours for the Jasper aquifer ranged from 200 ft below datum to 250 ft above datum. The 2000–19 groundwater-level-change contours for the Jasper aquifer depict groundwater-level declines throughout most of the study area where groundwater-level-altitude data from the Jasper aquifer were collected, with the largest decline (200 ft) in southern Montgomery County.

Texas↗
USGS2020

Water-balance techniques for determining available soil-water storage for selected sandy and clay soil study sites in Cass County, North Dakota, 2016–17

The U.S. Geological Survey, in cooperation with the U.S. Department of Agriculture Natural Resources Conservation Service, collected field and remotely sensed data on precipitation, evapotranspiration (ET), and soil-water content to determine available soil-water storage (AWS) at six study sites on sandy and clay soils in Cass County, North Dakota. Data were collected at all the study sites from May 1–October 31, 2016, and from May 1–October 24, 2017. Estimated daily AWS was determined using daily meteorological and potential evapotranspiration (PET) data obtained from various climate stations, and estimated monthly AWS was determined using monthly meteorological and PET data and monthly ET data determined using the Operational Simplified Surface Energy Balance model. AWS during 2016 and 2017 was determined at daily and monthly time steps because of data availability and to assess results using varying time steps. Comparisons of measured and estimated daily values of AWS at the Brewer Lake site indicated poor agreement during May–October 2016 and May–October 2017. Comparisons of measured and estimated daily values of AWS at the Embden East and Embden West sites indicated poor and fair agreement respectively. At the Lynchburg Crop and Lynchburg Grass sites, comparisons of measured and estimated daily values of AWS indicated fair and good relations, respectively, even with the possible effects of soil cracks. Mean estimated values of daily runoff plus soil percolation for the four sandy soil sites indicated that a maximum of about 19 percent of the estimated runoff plus soil percolation could be considered runoff and that the remaining 81 percent could be considered soil percolation, and for the two clay soil sites about 13 percent of the runoff plus soil percolation could have been considered runoff and about 87 percent could have been considered soil percolation. Results indicated little difference between using monthly PET or monthly ET in water-balance equations to estimate monthly AWS for the grouped sandy soil sites, and only slightly better results were obtained using monthly PET than monthly ET to estimate monthly AWS for the grouped clay soil study sites. Overall, the monthly water-balance models did not perform as well as the daily water-balance models for determining AWS at the six study sites. Additional data collection from a longer-period study and adjustments to the models may improve results from the monthly water-balance techniques.

North Dakota↗
USGS2020

Sediment lithology and borehole erosion testing, American and Sacramento Rivers, California

Executive Summary A field investigation intended to measure the potential for erosion of sediments beside the American and Sacramento Rivers near Sacramento, California, is described. The study featured two primary components: (1) drilling and soil sampling to reveal lithology, down to depths matching the local river thalweg, where possible, and (2) borehole erosion tests (BETs) as described by Briaud and others (2017) at many of the same locations. The latter test involved drilling a vertical hole, measuring its diameter profile, inserting a hollow drilling rod to almost the bottom of the hole, and pumping fluid through the drilling rod at a known discharge for a chosen time interval. The hole was then resurveyed to establish an erosion rate (change in borehole radius divided by duration of flow event) as a function of depth, and the test was repeated. This test was performed with water as the erosive fluid at 12 locations, with 1 test repeated with drilling mud. Lithology holes were drilled at these same locations and an additional five locations. Drilling operations took place on river left and river right on the American River and river left (left bank, when looking downstream) on the Sacramento River. The drilling to acquire sediment samples and reveal lithology involved the use of a mobile drilling rig equipped with a 6-inch (in.) auger, a 140-pound pneumatic hammer to drive split spoon and Calmod samplers, and a piston to push Shelby tube samplers to obtain samples of clayey material. Blow count (hammer blows per 6-in. sampler advance) was recorded while sampling, and the process was logged using standard U.S. Army Corps of Engineers (USACE), Sacramento District procedures. Sediment samples were identified and described in the field per ASTM D2488 and then delivered to a USACE laboratory and to Texas A&M University for additional laboratory analysis. The BETs were performed with the same drilling rig that performed the drilling for definition of lithology. In most instances, tests were limited to regions above the water table, to avoid slumping of the borehole and heaving sands pushing into the hole. Most of the tests featured sediments that were primarily silty sand or sandy silt. The testing procedure involved comparing borehole profiles before and after passing an assumed constant discharge through a drilling rod to the bottom of the drilled hole. Discharge and water losses were logged during the testing procedure, and water losses into the walls of the drilled hole were typically less than 5 percent of the introduced volume. For the tests performed with water, the coefficient of variation of the discharge ranged from 4.5 to 28 percent, with a mean of 13 percent, but the mean discharge appeared to be reasonably steady over the typical test duration of 10–30 minutes. It was thus assumed that discharge was constant and water losses during the tests were neglected. Coefficients of variation of the discharge for the three tests performed with drilling mud were much higher (20–50 percent), but erosion rates were much smaller. Resolution of the borehole caliper-reported diameter was 0.1 in. and several of the tests lasted for 10 minutes. With boreholes measured twice, before and after each test, and averaged, these numbers correspond to an apparent erosion rate (radius change divided by test duration) of 0.3 inches per hour (in/hr), which is a theoretical lower bound on what could be measured with this approach and equipment. In practice, 0.5 in/hr appears to be a more realistic lower bound on the detectable erosion rate, based on inspection of computed changes and erosion rates. Three flow speeds (5, 8, and 12 feet per second; ft/s) were targeted for the tests. Because of equipment limitations, it was not possible in the field to reach an average of 12 ft/s throughout any given borehole, although much higher flow speeds were reached locally in some cases. Most tests featured at least two different flow rates, and the borehole was typically surveyed at least twice for each condition, to allow averaging to reduce the influence of random diameter measurement errors. Errors arising from out-of-round boreholes appeared to be uncommon. Briaud and others (2017) recommend stepped increases in the flow rate during a borehole test. This approach was taken during initial testing but proved to be problematic. The drilled hole would be enlarged by the first (smaller) discharge, and then it would be difficult to reach the desired higher flow speed because of the larger annulus between the drilled hole and the drilling rod that supplied the water for testing. This was largely solved by starting with a high discharge and, in many cases, maintaining it for subsequent tests with the average flow speed decreasing as the hole enlarged. Several different measures of erosion rate were computed and investigated by comparison to lithological profiles. The vertically averaged erosion rate for each hole was computed, but this result does not reveal vertical variability of erodibility; and the mean flow speed within the hole is not a good representation of the speed when attempting to determine a relationship between erosion rate and flow speed. Instead, for each 6-inch layer within the hole, vertically averaged erosion rates and local flow speeds were computed and plotted. Where possible, the soil type for each layer was identified. For later laboratory analysis, project protocol dictated collection of Shelby tube samples whenever clay was encountered. Plots of erosion rate versus flow speed displayed scatter that indicate that several other factors influence the erosion potential of the soil. Blow count was not a good predictor variable; it is better correlated with soil type than erodibility. Soils were classified as sand, silt, or clay, depending on which soil type dominated within a sample. In general, those classified as sand and silt did not reveal clear patterns allowing erosion rate to be computed directly from flow speed, but the test results define the range and bounds on the erosion rate. Results for clay were slightly clearer with the erosion rate increasing with flow speed, once a threshold had been reached. In this case, the erosion rate appeared to change near a speed of 7 ft/s; above this threshold, erosion rates jumped from less than 2 in/hr to greater than 3 in/hr. Even for soils with similar classifications, large differences in erodibility were observed between sites and in different layers within an individual hole. One potential means of dealing with this problem would be to perform more tests at each site to allow establishment of relationships between flow speed and erodibility for individual layers within a borehole. The maximum number of tests performed at a site in this study was four, but in some cases, results are available for only one or two flow events. Comparison of data to a set of Erosion Function Apparatus tests that provide better resolution of the vertical variation in the erosion rate versus flow speed relationship would allow further investigation of this idea. It was hypothesized that drilling mud could expand the utility of the test in soft sands by reducing the likelihood of slumping that would be interpreted as erosion. The one test that was performed with drilling mud indicated that it greatly reduced the erosion rate of the soils encountered. It yielded very different results from the test performed at the same site with water. Erosion rate is often expressed as a function of shear stress applied to a soil. In order to compute shear stress on the walls of the drilled hole, one must assume a form for the relationship between flow speed and shear stress and select a friction factor that is often estimated empirically from head loss, observed water-surface profiles, surface roughness, or other data not available in this report. One methodology for computing shear stress from flow speed is discussed in this report, but the test results have been presented in terms of erosion rate versus flow speed to avoid assuming values that are not verifiable via the field data collected in this study. Erosion rate was computed from directly measured values (sequential borehole profiles) and flow speed was computed directly from measured quantities (discharge and borehole geometry). The BET has seen limited application, primarily in clayey soils, whereas most of the soils encountered in this study were primarily sand or silt. The objective of the BET is to determine the erodibility of in situ soil below the ground or riverbed surface. The BET is simple in principle and has the advantage of revealing erodibility of in situ sediments below the ground or riverbed surface; it appears to be very useful in clayey soils, based on previously published work, but is more difficult to apply in sandy soils where slumping and water losses within the hole during testing are more likely to occur. The BET did reveal a large variation in the results both laterally and vertically, even for the same soil-type classification. It is thus recommended that the results be applied considering these spatial variations rather than attempting to universally assign an erosion-rate relationship to a particular soil type. Results have been provided showing the results by site and by sediment classification (sand, silt, and clay), to allow either approach. Where possible, it is important to rely on site-specific results because the erosion-rate relationship for a given soil type varied by site. Data collected during this project have been made publicly available online via the U.S. Geological Survey (USGS) Sciencebase database. The measured borehole profiles, discharge, lithology log sheets, and photos are available in the data release that accompanies this report (see Work and Livsey (2019) in the “Selected References” section for the appropriate link).

California↗
USGS2020

Stressor identification framework of biological impairment in Mississippi streams to support watershed restoration and TMDL development

The Clean Water Act (CWA) requires States to identify waters that are impaired for designated uses. These waters are published through a State’s §303(d) list. The CWA also requires that a total maximum daily load (TMDL) be completed for each water body to calculate the maximum amount of contaminants that can be present in that water body and still meet water-quality standards. The Mississippi Department of Environmental Quality (MDEQ) uses a statewide monitoring and assessment strategy to collect benthic macroinvertebrate community data to assess the health of streams and rivers and to identify impaired waters. Waters that are found to be impaired based on the macroinvertebrate community data are listed on the Mississippi §303(d) list, and the cause of impairment is listed as “biological impairment.” Although the CWA requires TMDLs to be developed for applicable contaminants identified in the §303(d) list, TMDLs cannot be computed for stream reaches in Mississippi listed for biological impairment because the actual stressors causing the impairment have not yet been determined. The MDEQ and other water-resource managers in Mississippi require a framework for stressor identification in biologically impaired streams and rivers. This report is organized to (1) provide a general overview of biological impairment and stressor identification in stream ecosystems and (2) provide a detailed framework for stressor identification of Mississippi streams that are biologically impaired. The intent is for the framework to reduce subjectivity, provide consistency, and allow for adaptation as the science evolves. The stressor identification framework for Mississippi involves six key steps: Define the impairment, List the candidate causes of impairment and develop a conceptual model, Compile all relevant data, Evaluate the data, Identify probable causes of impairment by using a weight-of-evidence approach, and Generate a report of results.

Mississippi↗
USGS2024

Simulation of hydrodynamics and water temperature in a 21-mile reach of the upper Illinois River, Illinois, 2020–22

This report describes the development of a CE-QUAL-W2 river hydrodynamics and temperature model of a 21-mile reach of the Illinois River including a 3-mile reach of a major tributary, the Fox River. Model outputs consist of streamflow, water velocity, water-surface elevation, and water-temperature time series that can be used to simulate summer conditions in years with and without extensive development of harmful algal blooms (HABs). These analyses may provide a better understanding of some complex factors contributing to HAB development along the Illinois River. Such an understanding may provide more accurate HAB timing and location predictions and may help determine potential mitigating activities to prevent or limit the size and duration of HABs. Using the observed and simulated hydrodynamic conditions in the Illinois River study reach, it was possible to compare and contrast streamflow, velocity, and temperature conditions in years with varying HAB distributions. Occurrences of extensive HABs were documented in the study reach in June 2020 and June 2021, but only a small HAB restricted to the Marseilles Lock and Dam pool occurred in the summer of 2022. The objective then was to find similarities in site conditions between 2020 and 2021 that may contrast with the conditions in 2022. Among the 3 years included in the study, the variability in simulated water temperature exceeded variability in observed streamflow and simulated velocities. The longest period of water temperatures greater than 27 degrees Celsius (°C) in the selected locations in June of the three analysis years was in the second half of June 2022, yet no study-area wide HAB was documented in 2022. Simulations indicated that after warm water temperatures were established in the reach in June 2022, a cooling period broke up the warming period. This period of cooling was greater in magnitude and duration downstream from the location of a localized HAB perhaps limiting the spread of the bloom. Residence times differed substantially in segments representing different channel features; values ranged from 0.28 to 17.3 (days per 500 meters of channel) between the main stem and backwater areas, respectively. Variation in average June residence times was also greater among different channel features than among different years in the study period. The HABs in 2020 and 2021 at Starved Rock Dam were documented when water temperatures were about 26 °C. River backwater areas at some locations did attain these temperatures 2 to 3 days before the conditions in the main stem. Residence times in the backwater areas, however, generally exceeded 9 days, thus limiting the exchange of water carrying algal biomass into the main channel. Hydrodynamic model calibration involved adjusting model parameters until observed and simulated daily water-surface elevations, daily streamflows, discrete velocities, and channel areas were similar. Temperature calibration was done with near-surface continuous time-series data and discrete vertical profile temperatures. Observed and simulated water temperatures generally were within 1 °C at all monitoring locations.

Illinois↗
USGS2024

Correlation analysis of groundwater and hydrologic data, Kaloko-Honokōhau National Historical Park, Hawai‘i

Designated in 1978, Kaloko-Honokōhau National Historical Park is located on the west coast of the Island of Hawaiʻi. The Kaloko-Honokōhau National Historical Park encompasses about 1,200 acres of coastal land and nearshore ecosystems, which include wetlands, anchialine pools (landlocked bodies of brackish water with hydrologic connections to the ocean), fishponds, a fishtrap, and coral reefs. These nearshore ecosystems are dependent on groundwater discharge with a freshwater component and provide habitat for threatened and endangered, endemic species, such as the orangeblack Hawaiian damselfly ( Megalagrion xanthomelas ) and the Hawaiian coot (ʻAlae keʻokeʻo, Fulica alai ). The populations of these native species, however, are threatened because of habitat loss related to urban development and environmental changes. Kaloko-Honokōhau National Historical Park is within the Keauhou aquifer system and the North Kona District, which experienced a 52 percent resident-population increase between 2000 and 2020 and a 41 percent visitor increase between 2008 and 2019. To support the current water demand associated with this growing population, groundwater is the primary source of freshwater used in the North Kona District, with about 15 million gallons of groundwater withdrawn from the Keauhou aquifer system per day since 2009. With anticipated development, future (2015–35) groundwater withdrawal from the Keauhou aquifer system is projected to be about 55 percent greater than recent (2012–14) withdrawal. Because Kaloko-Honokōhau National Historical Park is located within a coastal aquifer, natural and human-induced changes can affect the quality and quantity of groundwater, which can threaten groundwater-dependent ecosystems. To improve understanding of recent groundwater conditions, the U.S. Geological Survey, in cooperation with the National Park Service, undertook this study to document correlations between hydrologic time-series datasets from sites in and near Kaloko-Honokōhau National Historical Park using the nonparametric (distribution-free) Kendall’s tau statistical test. For the statistical analyses, dependent variables representing the groundwater system include groundwater level, the groundwater-level difference between pairs of sites, and specific conductance, and independent variables include datasets of sea level, rainfall, and groundwater withdrawal. About 34 percent of the 140 non-time-lagged Kendall’s tau statistical tests evaluated in this report are statistically significant ( p -value ≤ 0.050) with generally weak (0.1 ≤ tau ≤ 0.2) to moderate (0.2 ≤ tau ≤ 0.3) correlations. Groundwater levels measured at monitoring sites have the strongest correlation with the multivariate El Niño–Southern Oscillation index and withdrawal from production wells at the nearby Kohanaiki Private Club Community. Specific conductance is not consistently and significantly correlated with the independent hydrologic variables investigated in this report. Because the relations between hydrologic variables are commonly not instantaneous, a second set of correlations was evaluated after applying a range of time lags to the independent variable datasets. Relative to the non-time-lagged case (the set of correlations that did not use time-lagged independent variables), some of the time-lagged independent variables improved correlations with some of the dependent variables. For a particular independent variable, similar time lags were expected between the independent variable and dependent variable at all four monitoring sites. However, different time lags among the four sites sometimes produced the strongest correlations. This study identified several correlations that are statistically significant and hydrologically plausible, but the correlations could indicate that multiple concurrent factors are controlling the observed groundwater-system response, which might be better addressed using multivariate analyses. This study only investigates bivariate correlations, which may not explain all the variance in the data. The correlations analyzed in this report are limited by the quantity of available hydrologic data in the area near Kaloko-Honokōhau National Historical Park and are based on 14 years of time-series data, which were aggregated to a relatively coarse monthly temporal resolution that represents the minimum resolution common to all datasets.

Hawaii↗
USGS2025

Methods for peak-flow frequency analysis for streamgages in or near Montana, North Dakota, South Dakota, and Wyoming

The U.S. Geological Survey, in cooperation with the Montana Department of Natural Resources and Conservation, North Dakota Department of Water Resources, South Dakota Department of Transportation, and the Wyoming Water Development Office, has developed standard methods of peak-flow frequency analysis for studies in Montana, North Dakota, South Dakota, and Wyoming. These methods describe the implementation of national flood frequency guidelines described in Bulletin 17C ( https://doi.org/10.3133/tm4B5 ) for the four States and deviations from Bulletin 17C standard procedures to accommodate unusual hydrologic conditions. A U.S. Geological Survey data release accompanying this report ( https://doi.org/10.5066/P1WHRK8H ) provides example peak-flow frequency analyses for selected streamgages in the study area. The methods described in this report can be used to publish similar data releases for other streamgages in the study area.

Montana, North Dakota, South Dakota, Wyoming↗
USGS2017

Reconstruction of spatio-temporal temperature from sparse historical records using robust probabilistic principal component regression

Scientific records of temperature and precipitation have been kept for several hundred years, but for many areas, only a shorter record exists. To understand climate change, there is a need for rigorous statistical reconstructions of the paleoclimate using proxy data. Paleoclimate proxy data are often sparse, noisy, indirect measurements of the climate process of interest, making each proxy uniquely challenging to model statistically. We reconstruct spatially explicit temperature surfaces from sparse and noisy measurements recorded at historical United States military forts and other observer stations from 1820 to 1894. One common method for reconstructing the paleoclimate from proxy data is principal component regression (PCR). With PCR, one learns a statistical relationship between the paleoclimate proxy data and a set of climate observations that are used as patterns for potential reconstruction scenarios. We explore PCR in a Bayesian hierarchical framework, extending classical PCR in a variety of ways. First, we model the latent principal components probabilistically, accounting for measurement error in the observational data. Next, we extend our method to better accommodate outliers that occur in the proxy data. Finally, we explore alternatives to the truncation of lower-order principal components using different regularization techniques. One fundamental challenge in paleoclimate reconstruction efforts is the lack of out-of-sample data for predictive validation. Cross-validation is of potential value, but is computationally expensive and potentially sensitive to outliers in sparse data scenarios. To overcome the limitations that a lack of out-of-sample records presents, we test our methods using a simulation study, applying proper scoring rules including a computationally efficient approximation to leave-one-out cross-validation using the log score to validate model performance. The result of our analysis is a spatially explicit reconstruction of spatio-temporal temperature from a very sparse historical record.

Advances in Statistical Climatology, Meteorology a↗
USGS1995

USGS supports ecosystem management in the San Francisco bay and delta

In the past several years, the Department of the Interior has ,l?laced particular emphasis on "ecosystem management" - the integration of scientific knowledge of ecological relationships with resource management practices to sustain ecological, cultural, and economic systems in broad habitat areas; eg., forest, desert, and aquatic habitats. The goal of ecosystem management is to understand the habitat requirements of many species and thereby assess the impact of a variety of human activities on regional biodiversity rather than to manage for individual species. In support of the management of the natural resources of San Francisco Bay and the delta, the USGS is augmenting its ongoing research program to provide additional research, monitoring, and assessment information and to improve the availability of that information. Following a workshop with resource and regulatory agency representatives, hosted by the USGS in early 1994, USGS researchers identified specific areas where existing and new scientific information could contribute to resolving some of the critical bay and delta land and aquatic management issues. The "ecosystems initiative" effort now underway addresses questions in three areas: Freshwater flow, including sediment and contaminant transport. Wetlands. Information and data access. Following is a brief description of activities being supported by this initiative.

Interagency Ecological Program Newsletter↗
USGS2008

Scientific results from Gulf of Mexico Gas Hydrates Joint Industry Project Leg 1 drilling: Introduction and overview

The Gulf of Mexico Gas Hydrates Joint Industry Project (JIP) is a consortium of production and service companies and some government agencies formed to address the challenges that gas hydrates pose for deepwater exploration and production. In partnership with the U.S. Department of Energy and with scientific assistance from the U.S. Geological Survey and academic partners, the JIP has focused on studies to assess hazards associated with drilling the fine-grained, hydrate-bearing sediments that dominate much of the shallow subseafloor in the deepwater (>500 m) Gulf of Mexico. In preparation for an initial drilling, logging, and coring program, the JIP sponsored a multi-year research effort that included: (a) the development of borehole stability models for hydrate-bearing sediments; (b) exhaustive laboratory measurements of the physical properties of hydrate-bearing sediments; (c) refinement of new techniques for processing industry-standard 3-D seismic data to constrain gas hydrate saturations; and (d) construction of instrumentation to measure the physical properties of sediment cores that had never been removed from in situ hydrostatic pressure conditions. Following review of potential drilling sites, the JIP launched a 35-day expedition in Spring 2005 to acquire well logs and sediment cores at sites in Atwater Valley lease blocks 13/14 and Keathley Canyon lease block 151 in the northern Gulf of Mexico minibasin province. The Keathley Canyon site has a bottom simulating reflection at ???392 m below the seafloor, while the Atwater Valley location is characterized by seafloor mounds with an underlying upwarped seismic reflection consistent with upward fluid migration and possible shoaling of the base of the gas hydrate stability (BGHS). No gas hydrate was recovered at the drill sites, but logging data, and to some extent cores, suggest the occurrence of gas hydrate in inferred coarser-grained beds and fractures, particularly between 220 and 330 m below the seafloor at the Keathley Canyon site. This paper provides an overview of the results of the initial phases of the JIP work and introduces the 15 papers that make up this special volume on the scientific results related to the 2005 logging and drilling expedition.

Marine and Petroleum Geology↗
USGS2026

Decision-centered systems for the future of wildlife conservation

The future of wildlife conservation will be determined by how well we decide, not simply by how much we know. Advances in quantitative ecology, sensor networks, and artificial intelligence have enhanced the field's ability to estimate population states and predict ecological dynamics. Yet, these technological and analytical advancements have not translated into proportional improvements in wildlife conservation decision making. The prevailing science-centered paradigm, which prioritizes research production as the primary route to better management, was appropriate when limited ecological understanding was the primary constraint on decision making. Today, wildlife conservation challenges are more often restricted by inadequate decision-making processes rather than by insufficient information. This paper defines a decision-centered system as an institutional and organizational arrangement in which the generation, integration, and application of scientific knowledge are explicitly structured around the decisions required for wildlife conservation. Decision systems evaluate success not only through scientific rigor, but through the quality of decision-making processes and outcomes they support. Artificial intelligence (AI) is likely to deepen this mismatch by rapidly expanding analytical capacity while leaving the underlying challenges of decision making unresolved. As AI data synthesis, modeling, and estimation become increasingly automated, the comparative advantage of conservation professionals is likely to shift toward tasks requiring human judgment (e.g., eliciting stakeholder values, weighting competing management objectives). The conceptual basis for this shift already exists in structured decision making, adaptive management, and similar frameworks. The remaining barriers to implementation are both institutional and organizational. The rules, norms, and incentive structures governing what counts as good conservation work, including what gets funded, published, and rewarded, will benefit from shifting alongside the agencies, universities, and professional societies that train and employ wildlife professionals. Conditions for this transition include redesigning incentives, training, and governance so that organizations are prepared not only to generate knowledge, but to consistently translate it into effective action.

The Wildlife Society Bulletin↗
USGS2023

Field observations and long short-term memory modeling of spectral wave evolution at living shorelines in Chesapeake Bay, USA

Living shorelines as a nature-based solution for climate change adaptation were constructed in many places around the world. The success of this type of projects requires long-term monitoring for adaptive management. The paper presents a novel framework leveraging scientific machine learning methods for accurate and rapid prediction of long-term hydrodynamic forcing impacting living shorelines using short-term measurements of water levels and wind waves in the largest estuary in the U.S. Different from existing data-driven wave prediction models focusing on significant wave heights, this study is focused on the prediction of wave energy spectra in shallow water using winds and tides as the input feature and short-term measurements of wave spectra and water depths as the label. Long Short-Term Memory (LSTM) models were developed using four-month wave measurements in the stormy seasons to predict integral wave parameters and energy spectra for multiple years. The developed models accurately predicted wave heights, peak periods, and energy spectra around the living shorelines, capturing complex wave dynamics, such as wave generation by wind, nonlinear wave-wave interactions, and depth-limited wave breaking in the shallow water of a large estuary. The validated models were then used to determine the long-term wave forcing impacting the living shorelines based on the modeled wave characteristics and spectra. Model results show that the surrogate models utilizing LSTM to predict wave spectra in the frequency domain enable long-term predictions of spectral wave evolution with a minimal computational cost. Our findings provide valuable insights into the efficacy of living shorelines in attenuating wave energy and demonstrate the utility of this approach in assessing the effectiveness of such living shoreline structures.

Applied Ocean Research↗
USGS2017

Accounting for sampling patterns reverses the relative importance of trade and climate for the global sharing of exotic plants

Aim The distributions of exotic species reflect patterns of human-mediated dispersal, species climatic tolerances and a suite of other biotic and abiotic factors. The relative importance of each of these factors will shape how the spread of exotic species is affected by ongoing economic globalization and climate change. However, patterns of trade may be correlated with variation in scientific sampling effort globally, potentially confounding studies that do not account for sampling patterns. Location Global. Time period Museum records, generally from the 1800s up to 2015. Major taxa studied Plant species exotic to the United States. Methods We used data from the Global Biodiversity Information Facility (GBIF) to summarize the number of plant species with exotic occurrences in the United States that also occur in each other country world-wide. We assessed the relative importance of trade and climatic similarity for explaining variation in the number of shared species while evaluating several methods to account for variation in sampling effort among countries. Results Accounting for variation in sampling effort reversed the relative importance of trade and climate for explaining numbers of shared species. Trade was strongly correlated with numbers of shared U.S. exotic plants between the United States and other countries before, but not after, accounting for sampling variation among countries. Conversely, accounting for sampling effort strengthened the relationship between climatic similarity and species sharing. Using the number of records as a measure of sampling effort provided a straightforward approach for the analysis of occurrence data, whereas species richness estimators and rarefaction were less effective at removing sampling bias. Main conclusions Our work provides support for broad-scale climatic limitation on the distributions of exotic species, illustrates the need to account for variation in sampling effort in large biodiversity databases, and highlights the difficulty in inferring causal links between the economic drivers of invasion and global patterns of exotic species occurrence.

Global Ecology and Biogeography↗