USGS Science⌕ Search

SEARCH · USGS Science

Results for “Modeling Report”

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 631 records · Page 35Linked to original sources

Continuous monitoring of nutrient and sediment loads from the Des Plaines River at Route 53 at Joliet, Illinois, water years 2018–20

The Des Plaines River in southern Wisconsin and northern Illinois is the principal conduit for the discharge of wastewater effluent and stormwater runoff from the greater Chicago metropolitan area. In November 2017, the U.S. Geological Survey, in cooperation with the Metropolitan Water Reclamation District of Greater Chicago, installed a continuous monitoring station to measure water quality and streamflow in the Des Plaines River at Joliet, Illinois. Surrogate models encompassing continuous data and discrete water-quality samples were used to estimate loads of nitrate, total phosphorus, and suspended sediment. Comparisons to other major rivers in Illinois show that the Des Plaines River is a substantial contributor to statewide loading estimates for nitrate and total phosphorus but only a minor contributor to suspended sediment. Future loading estimates of total phosphorus could include more research into the effects of combined sewage overflows because these effects likely increased model uncertainty. The results in this report document current loadings and provide a baseline from which to assess future water-quality management decisions.

Illinois↗

Characterization of groundwater resources near the southeastern part of Puget Sound, Washington

More than 1 million people live within the active model area (AMA) in the southeastern part of the lowlands surrounding Puget Sound, or Puget Lowland, Washington, and groundwater is the source for approximately one-half of their public, domestic, and irrigation water demands. The 887-square-mile AMA, located in King and Pierce Counties, represents the area of analysis for the conceptual hydrogeologic framework and numerical groundwater-flow models within the study area and includes the Puyallup River and Chambers-Clover Creek watersheds. To assess the potential hydrologic and anthropogenic impacts to groundwater and the connected surface-water resources, conceptual and numerical groundwater-flow models of groundwater flow were developed by the U.S. Geological Survey Washington Water Science Center in close cooperation with 18 water-resource agencies and stakeholders. This multichapter volume documents the development of the conceptual and numerical groundwater-flow models of groundwater flow. Chapters A, B, and C provide an overall introduction to the multichapter volume (Chapter A), the conceptual hydrogeologic framework (Chapter B), and the groundwater budget (Chapter C). Chapters D and E describe numerical groundwater-flow model construction and calibration (Chapter D) and the numerical groundwater-flow model results (Chapter E). Collectively, these reports present a characterization and simulation tool for groundwater resources near the southeastern part of Puget Sound, Washington.

Scientific Investigations Report↗

ESTIMATION OF URBAN STORM-RUNOFF LOADS.

The United States was divided into three regions, on the basis of mean annual rainfall, to decrease the variability in storm-runoff constituent loads and to improve regression relations with basin and climatic characteristics. Multiple-regression analyses, in progress, are being refined to determine the best regression models for each of the storm-runoff constituent loads in each of the three regions. These techniques, when finalized, can be used to estimate storm-runoff constituent loads for gaged and ungaged urban watersheds. The preliminary standard errors of estimate for five constituents examined to date ranged from 54 to 223 percent, and the coefficients of determination (r**2) ranged from 0. 39 to 0. 94. Total storm rainfall and total contributing drainage area appear to be the most significant independent variables in the regression models. This paper is a progress report on preliminary results of estimating storm-runoff constituent loads for ungaged watersheds.

Conference Paper↗

The influence of precipitation, vegetation and soil properties on the ecohydrology of sagebrush steppe rangelands on the INL site

The INL Site and other landscapes having sagebrush steppe vegetation are experiencing a simultaneous change in climate and floristics that result from increases in exotic species. Determining the separate and combined/interactive effects of climate and vegetation change is important for assessing future changes on the landscape and for hydrologic processes. This research uses the 72 experimental plots established and initially maintained for many years as the “Protective Cap Biobarrier Experiment” by Dr. Jay Anderson and the Stoller ESER program, and the experiment is also now referred to as the “INL Site Ecohydrology Study.” We are evaluating long-term impacts of different plant communities commonly found throughout Idaho subject to different precipitation regimes and to different soil depths. Treatments of amount and timing of precipitation (irrigation), soil depth, and either native/perennial or exotic grass vegetation allow researchers to investigate how vegetation, precipitation and soil interact to influence soil hydrology and ecosystem biogeochemistry. This information will be used to improve a variety of models, as well as provide data for these models.

Report↗

Annual survival and site fidelity of northern pintails banded on the Yukon-Kuskokwim Delta, Alaska

We banded northern pintails (Anas acuta; n = 13,645) at a single site on the Yukon-Kuskokwim Delta (YKD), Alaska, USA, from 1990 to 2001. We used recaptures from our site in combination with hunter recoveries to model annual survival, recovery rates, and fidelity to our capture location. Most recoveries (>90%) occurred in the Pacific Flyway with 64% reported from California's Central Valley. Our top candidate models allowed survival to vary by sex but not by age or year. Estimated annual survival was 77.6% (95% CI: 73.9-81.0%) for males and 60.2% (95% CI: 53.2-67.0%) for females. Reporting rates varied by age, sex, and year; estimates for adult males exceeded those for adult females by 3.5 times. Within sexes, reporting rates of hatch-year pintails exceeded those of adults. Estimated recovery rates were considerably lower than those estimated during the 1950s-1970s for winter banded pintails (Hestbeck 1993b), but there were no differences in survival rates. This suggests that changes in harvest regulations may not have influenced annual survival in this population. The propensity of banded pintails to return to our capture site (fidelity rate) varied between sexes and was positively correlated with water conditions in prairie Canada. Our estimates of fidelity rates varied from 77.4% to 87.2% for males and 89.8% to 94.3% for females. Our fidelity estimates suggest that some level of subpopulation structuring may exist for northern pintails. Additionally, our estimates of fidelity support previous observations of northern pintails overflying poor wetland habitat conditions on the Canadian prairies.

Journal of Wildlife Management↗

Comparison of groundwater-model construction methods, representations of glacial geology, model designs, and groundwater-model flow simulations within Elkhart County, Indiana

Automated data-processing methods allow hydrologists to efficiently incorporate digital well-record datasets into the construction of hydrostratigraphic frameworks for groundwater-flow models. The method selected to construct the hydrostratigraphic framework can affect the extent of geologic heterogeneity that can be included in the model. The detail generated from a hydrostratigraphic framework can affect groundwater simulation results. The effects of detail on model accuracy, groundwater-flow simulations, and particle-tracking simulations are described in this study. This report compares differences in hydrostratigraphic frameworks and results of groundwater models using (1) a method that incorporates more hydrologic judgment at the expense of using limited lithologic data and (2) a method that is more automated and uses all available lithologic data. The study additionally evaluates the effect of model discretization and inclusion of more (or less) geologic detail on simulation results. Two methods were used to create hydrostratigraphic frameworks of glacial deposits in the St. Joseph River Basin. One method, referred to as the subjective method, manually identifies stratigraphic boundaries using a sample of well logs from State databases and uses two-dimensional kriging to create three model layers of the study area. Indicator kriging is used to define aquifer extent in each layer. The second method, referred to as the objective method, uses three-dimensional kriging to automatically create a detailed heterogeneous model of the study area using all wells logs from the State database. The objective method increases detail in the vertical by greatly increasing the number of computer groundwater model layers from 3 to 30. In Elkhart County, Indiana, a previously published model represents the product of the subjective method, and a newly calibrated model of the same area represents the product of the objective method. An automated calibration procedure was used with the objective model (derived from the objective method) for Elkhart County. The two most-sensitive parameters for the Elkhart County objective model are horizontal hydraulic conductivity of the sand and the combined sand and gravel/gravel deposits. Vertical hydraulic conductivity of the fine-grained and intermediate-sized deposits could not be estimated, possibly indicating major flow paths are along a continuously connected series of sand and gravel deposits and not through a confining layer. The statistics measuring model calibration accuracy for the objective model were slightly better than statistics for the subjective model (model derived from the subjective method) of Elkhart County, but the hydraulic conductivities and flow rates for the two models were different. The mean absolute errors between simulated and measured groundwater levels are 2.04 and 2.16 feet for the objective and subjective models, respectively. Simulated seepage losses from and groundwater discharges to measured stream reaches in the objective model were evenly balanced in terms of over and under simulations of measured values; the subjective model tended to overpredict measured groundwater discharge to streams. The overprediction may be related to the 58 percent greater total inflow and outflow through the subjective model. The greater flow rate through the subjective model results from higher horizontal hydraulic conductivities in the subjective model than in the objective model. Horizontal hydraulic conductivity ranged from 23.9 to 111 feet per day in the objective model and generally ranged from 170 to 370 feet per day in the subjective model. The improvement in calibration statistics for the objective model relative to the subjective model may be from increased detail in how the objective model represents the distribution of fine- and coarse-grained deposits. The improvement also could be associated with the difference in methods used to represent the continuity of the confining unit. The effect of differences in horizontal hydraulic conductivity distributions between the two models for Elkhart County is evident in the groundwater-flow paths simulated by the objective and subjective models. At a withdrawal well location, the flow lines produced by the objective model indicate a wider contributing area than that for the subjective model. The discontinuous confining unit represented in the objective model provided the opportunity for groundwater flow to split into an upper and lower path. The split in flow simulated by the objective model at one location was independently supported by bromide concentrations in groundwater; the subjective model did not duplicate the split in flow.

Indiana↗

Evidence for shallow megathrust slip across the Unalaska seismic gap during the great 1957 Andreanof Islands earthquake, eastern Aleutian Islands, Alaska

We reassess the slip distribution of the 1957 Andreanof Islands earthquake in the eastern part of the aftershock zone where published slip models infer little or no slip. Eyewitness reports, tide gauge data, and geological evidence for 9–23 m tsunami runups imply seafloor deformation offshore Unalaska Island in 1957, in contrast with previous studies that labeled the area a seismic gap. Here, we simulate tsunami dynamics for a suite of deformation models that vary in depth and amount of megathrust slip. Tsunami simulations show that a shallow (5–15 km deep) rupture with ~20 m of slip most closely reproduces the 1957 Dutch Harbor marigram and nearby >18 m runup at Sedanka Island marked by stranded drift logs. Models that place slip >20 km predict waves that arrive too soon. Our results imply that shallow slip on the megathrust in 1957 extended east into an area that presently creeps.

Geophysical Research Letters↗

Evaluation of the groundwater-flow model for the Ohio River alluvial aquifer near Carrollton, Kentucky, updated to conditions in September 2010

The Ohio River alluvial aquifer near Carrollton, Ky., is an important water resource for the cities of Carrollton and Ghent, as well as for several industries in the area. The groundwater of the aquifer is the primary source of drinking water in the region and a highly valued natural resource that attracts various water-dependent industries because of its quantity and quality. This report evaluates the performance of a numerical model of the groundwater-flow system in the Ohio River alluvial aquifer near Carrollton, Ky., published by the U.S. Geological Survey in 1999. The original model simulated conditions in November 1995 and was updated to simulate groundwater conditions estimated for September 2010. The files from the calibrated steady-state model of November 1995 conditions were imported into MODFLOW-2005 to update the model to conditions in September 2010. The model input files modified as part of this update were the well and recharge files. The design of the updated model and other input files are the same as the original model. The ability of the updated model to match hydrologic conditions for September 2010 was evaluated by comparing water levels measured in wells to those computed by the model. Water-level measurements were available for 48 wells in September 2010. Overall, the updated model underestimated the water levels at 36 of the 48 measured wells. The average difference between measured water levels and model-computed water levels was 3.4 feet and the maximum difference was 10.9 feet. The root-mean-square error of the simulation was 4.45 for all 48 measured water levels. The updated steady-state model could be improved by introducing more accurate and site-specific estimates of selected field parameters, refined model geometry, and additional numerical methods. Collection of field data to better estimate hydraulic parameters, together with continued review of available data and information from area well operators, could provide the model with revised estimates of conductance values for the riverbed and valley wall, hydraulic conductivities for the model layer, and target water levels for future simulations. Additional model layers, a redesigned model grid, and revised boundary conditions could provide a better framework for more accurate simulations. Additional numerical methods would identify possible parameter estimates and determine parameter sensitivities.

Kentucky↗

Assessment of groundwater/surface-water interaction and simulation of potential streamflow depletion induced by groundwater withdrawal, Uinta River near Roosevelt, Utah

Roosevelt City, Utah, asserts a need for an additional supply of water to meet municipal demands and has identified a potential location for additional groundwater development at the Sprouse well field near the West Channel of the Uinta River. Groundwater is commonly hydraulically linked to surface water and, under some conditions, the pumpage of groundwater can deplete water in streams and other water bodies. In 2008, the U.S. Geological Survey, in cooperation with Roosevelt City, the Utah Department of Natural Resources, and the Ute Indian Tribe, began a study to improve understanding of the local interconnection between groundwater and surface water and to assess the potential for streamflow depletion from future groundwater withdrawals at a potential Roosevelt City development location—the Sprouse well field near the West Channel of the Uinta River. In the study, streamflow gains and losses at the river/aquifer boundary near the well field and changes in those conditions over time were assessed through (1) synoptic measurement of discharge in the stream at multiple sites using tracer-dilution methods, (2) periodic measurement of the vertical hydraulic gradient across the streambed, and (3) continuous measurement of stream and streambed water temperature using heat as a tracer of flow across the streambed. Although some contradictions among the results of the three assessment methods were observed, results of the approaches generally indicated (1) losing streamflow conditions on the West Channel of the Uinta River north of and upstream from the Sprouse well field within the study area, (2) gaining streamflow conditions south of and downstream from the well field, and (3) some seasonal changes in those conditions that correspond with seasonal changes in stream stage and local water-table altitudes. A numerical groundwater flow model was developed on the basis of previously reported observations and observations made during this study, and was used to estimate potential streamflow depletion that might result from future groundwater withdrawals at the Sprouse well field. The model incorporates concepts of transient groundwater flow conditions including fluctuations in groundwater levels and storage, and the distribution of and temporal variations in gains to and losses from streamflow in the West Channel of the Uinta River near the Sprouse well field. Two predictive model simulations incorporated additional future discharge from the Sprouse well field totaling 325 acre-feet annually and biennially during summer months. Results of the predictive model simulations indicate that the water withdrawn by the additional pumping was derived initially from aquifer storage and then, with time, predominantly from streamflow depletion. By the 10th year of the predictive simulation incorporating annual summer pumping from an additional public-supply well in the Sprouse well field, the simulation results indicate that 89 percent of a future annual 325 acre-feet of discharge is derived from depletion of streamflow in the West Channel of the Uinta River. A similar result was observed in a predictive model simulating the same discharge rate but with the new well being pumped every other year.

Utah↗

Simulation of streamflow in the McTier Creek watershed, South Carolina

The McTier Creek watershed is located in the Sand Hills ecoregion of South Carolina and is a small catchment within the Edisto River Basin. Two watershed hydrology models were applied to the McTier Creek watershed as part of a larger scientific investigation to expand the understanding of relations among hydrologic, geochemical, and ecological processes that affect fish-tissue mercury concentrations within the Edisto River Basin. The two models are the topography-based hydrological model (TOPMODEL) and the grid-based mercury model (GBMM). TOPMODEL uses the variable-source area concept for simulating streamflow, and GBMM uses a spatially explicit modified curve-number approach for simulating streamflow. The hydrologic output from TOPMODEL can be used explicitly to simulate the transport of mercury in separate applications, whereas the hydrology output from GBMM is used implicitly in the simulation of mercury fate and transport in GBMM. The modeling efforts were a collaboration between the U.S. Geological Survey and the U.S. Environmental Protection Agency, National Exposure Research Laboratory. Calibrations of TOPMODEL and GBMM were done independently while using the same meteorological data and the same period of record of observed data. Two U.S. Geological Survey streamflow-gaging stations were available for comparison of observed daily mean flow with simulated daily mean flow-station 02172300, McTier Creek near Monetta, South Carolina, and station 02172305, McTier Creek near New Holland, South Carolina. The period of record at the Monetta gage covers a broad range of hydrologic conditions, including a drought and a significant wet period. Calibrating the models under these extreme conditions along with the normal flow conditions included in the record enhances the robustness of the two models. Several quantitative assessments of the goodness of fit between model simulations and the observed daily mean flows were done. These included the Nash-Sutcliffe coefficient of model-fit efficiency index, Pearson's correlation coefficient, the root mean square error, the bias, and the mean absolute error. In addition, a number of graphical tools were used to assess how well the models captured the characteristics of the observed data at the Monetta and New Holland streamflow-gaging stations. The graphical tools included temporal plots of simulated and observed daily mean flows, flow-duration curves, single-mass curves, and various residual plots. The results indicated that TOPMODEL and GBMM generally produced simulations that reasonably capture the quantity, variability, and timing of the observed streamflow. For the periods modeled, the total volume of simulated daily mean flows as compared to the total volume of the observed daily mean flow from TOPMODEL was within 1 to 5 percent, and the total volume from GBMM was within 1 to 10 percent. A noticeable characteristic of the simulated hydrographs from both models is the complexity of balancing groundwater recession and flow at the streamgage when flows peak and recede rapidly. However, GBMM results indicate that groundwater recession, which affects the receding limb of the hydrograph, was more difficult to estimate with the spatially explicit curve number approach. Although the purpose of this report is not to directly compare both models, given the characteristics of the McTier Creek watershed and the fact that GBMM uses the spatially explicit curve number approach as compared to the variable-source-area concept in TOPMODEL, GBMM was able to capture the flow characteristics reasonably well.

South Carolina↗

Effects of water use and land use on streamflow and aquatic habitat in the Sudbury and Assabet River Basins, Massachusetts

Water withdrawals from surface-water reservoirs and groundwater have affected streamflow in the Sudbury and Assabet River Basins. These effects are particularly evident in the upper Sudbury River Basin, which prompted the need to improve the understanding of water resources and aquatic habitat in these basins. In 2004, the U.S. Geological Survey, in cooperation with the Massachusetts Department of Conservation and Recreation, developed a precipitation-runoff model that uses Hydrologic Simulation Program-FORTRAN (HSPF) to evaluate the effects of water use and projected future water-use and land-use change on streamflow. As part of this study, the aquatic habitat in the basins and the effects of streamflow alteration also were evaluated. Chapter 1 of the report covers the development of the HSPF model that focuses on the upper Sudbury River Basin (106 square miles) but covers the entire Sudbury and Assabet River Basins (339 square miles). The model was calibrated to an 11-year period (1993-2003) using observed or estimated streamflow at four streamgages. The model was then used to simulate long-term (1960-2004) streamflows to evaluate the effects of average 1993-2003 water use and projected 2030 water-use and land-use change over long-term climatic conditions. Simulations indicate that the average 1993-2003 withdrawals most altered streamflow relative to no withdrawals in small headwater subbasins where the ratios of mean annual withdrawals to mean annual streamflow are the highest. The effects of withdrawals are also appreciable in other parts of the upper Sudbury River Basin as a result of the perpetuation of the effects of large withdrawals in upstream reaches or in subbasins that also have a high ratio of withdrawal to streamflow. The simulated effects of potential 2030 water-use and land-use change indicate small decreases in flows as a result of increased water demands, but these flow alterations were offset as a result of decreased evapotranspiration associated with the loss of deep-rooted vegetation. Simulations of reactivating production wells near the north end of Lake Cochituate indicate pumping could substantially affect lake levels and flows at the lake outlet or in nearby reaches in the Sudbury River during periods of low flow, but the effects vary depending on the source of the water to the wells, which is largely unknown. Chapter 2 of the report covers the fish-community assessment and comparison of streamflow-setting standards for protecting aquatic habitat. The fish-community assessment indicates the main stems of the Sudbury and Assabet Rivers are dominated by macrohabitat generalists. Water temperatures recorded in seven free-flowing reaches in the upper Sudbury River Basin at three sites unaffected by withdrawals or impoundments are generally suitable for cold-water fish; however, summer temperatures often rose to a level considered critical to long-term survival of brook trout. At four sites downstream from withdrawals or reservoirs, or both, summer water temperatures were often in the upper critical range for brook trout survival. Physically and statistically based methods for determining streamflows for protecting aquatic habitat were applied at 10 selected riffle sites in the Sudbury and Assabet River Basins. Physically based methods, R2Cross and Wetted-Perimeter, use site-specific physical and hydraulic information and a one-dimensional hydraulics model, HEC-RAS, to determine flows that meet the criteria set forth by the method. The median flow that meets 2-of-3 of the R2Cross hydraulic criteria (percentage of bankfull wetted perimeter, average velocity, and mean depth) ranged from about 0.07 to 0.72 cubic feet per second per square mile (ft3/s/mi2) with an overall median of about 0.24 ft3/s/mi2; the median Wetted-Perimeter target flow ranged from about 0.10 to 0.51 ft3/s/mi2 with an overall median of about 0.25 ft3/s/mi2. Statistically based methods?Tennant, New England Aquatic Base Flow (ABF)

Scientific Investigations Report↗

Estimating density and detection of bobcats in fragmented Midwestern landscapes using spatial capture-recapture data from camera traps

Camera-trapping data analyzed with spatially explicit capture–recapture (SCR) models can provide a rigorous method for estimating density of small populations of elusive carnivore species. We sought to develop and evaluate the efficacy of SCR models for estimating density of a presumed low-density bobcat ( Lynx rufus ) population in fragmented landscapes of west-central Illinois, USA. We analyzed camera-trapping data from 49 camera stations in a 1,458-km 2 area deployed over a 77-day period from 1 February to 18 April 2017. Mean operational time of cameras was 52 days (range = 32–67 days). We captured 23 uniquely identifiable bobcats 113 times and recaptured these same individuals 90 times; 15 of 23 (65.2%) individuals were recaptured at ≥2 camera traps. Total number of bobcat capture events was 139, of which 26 (18.7%) were discarded from analyses because of poor image quality or capture of only a part of an animal in photographs. Of 113 capture events used in analyses, 106 (93.8%) and 7 (6.2%) were classified as positive and tentative identifications, respectively; agreement on tentative identifications of bobcats was high (71.4%) among 3 observers. We photographed bobcats at 36 of 49 (73.5%) camera stations, of which 34 stations were used in analyses. We estimated bobcat density at 1.40 individuals (range = 1.00–2.02)/100 km 2 . Our modeled bobcat density estimates are considerably below previously reported densities (30.5 individuals/100 km 2 ) within the state, and among the lowest yet recorded for the species. Nevertheless, use of remote cameras and SCR models was a viable technique for reliably estimating bobcat density across west-central Illinois. Our research establishes ecological benchmarks for understanding potential effects of colonization, habitat fragmentation, and exploitation on future assessments of bobcat density using standardized methodologies that can be compared directly over time. Further application of SCR models that quantify specific costs of animal movements (i.e., least-cost path models) while accounting for landscape connectivity has great utility and relevance for conservation and management of bobcat populations across fragmented Midwestern landscapes.

Illinois↗

Comparison of Estimated Areas Contributing Recharge to Selected Springs in North-Central Florida by Using Multiple Ground-Water Flow Models

Areas contributing recharge to springs are defined in this report as the land-surface area wherein water entering the ground-water system at the water table eventually discharges to a spring. These areas were delineated for Blue Spring, Silver Springs, Alexander Springs, and Silver Glen Springs in north-central Florida using four regional ground-water flow models and particle tracking. As expected, different models predicted different areas contributing recharge. In general, the differences were due to different hydrologic stresses, subsurface permeability properties, and boundary conditions that were used to calibrate each model, all of which are considered to be equally feasible because each model matched its respective calibration data reasonably well. To evaluate the agreement of the models and to summarize results, areas contributing recharge to springs from each model were combined into composite areas. During 1993-98, the composite areas contributing recharge to Blue Spring, Silver Springs, Alexander Springs, and Silver Glen Springs were about 130, 730, 110, and 120 square miles, respectively. The composite areas for all springs remained about the same when using projected 2020 ground-water withdrawals.

Open-File Report↗

Observing systems, modeling, and forecasting

Predicting harmful algal blooms (HABs) requires integrating physical, chemical, and biological data collected from observing networks and then assimilating these data into models, which are used to generate forecasts. In 2005, the Harmful Algal Research and Response: A National Environmental Science Strategy 2005-2015 (HARRNESS, 2005) made recommendations on how to improve HAB modeling and forecasting over the next decade. Key HARRNESS recommendations related to sensing, networking, and modeling HABs included: ● Support the development and validation of new and improved technologies for remote cell and toxin detection, and for modeling and forecasting, ● Improve coordination of monitoring/ and modeling efforts, both at national and regional levels, ● Improve the use of networking technologies for monitoring and modeling efforts, ● Conduct sustained time series measurements of the biotic, chemical, and physical environments impacted by HABs, ● Develop food web models on the ecosystem fate and effects of toxins, ● Develop and improve species-specific models that link to physical-biological models. Here we review HAB observing, modeling, and forecasting advances and technologies and recommend research and management priorities for the next decade and beyond. Our report encompasses sensing technologies, sensor networking and data management, models and forecasts, and the paths to operationalize forecasts. Continued improvements of deployable sensors are foundational to improving early warning indicators, models, and forecasts, which are only as good as the underlying data. Sensing technology has advanced considerably in the last decade; for example, more capable fluorometric pigment sensors can track changes in bloom biomass in real-time. Additionally, automated imaging/classifying systems to identify and quantify key harmful algal (HA) taxa are being routinely deployed. However, deployable toxin sensors are available for only some HAB toxins and continue to be identified as a critical need by researchers and managers. As more and improved sensors and technologies become available, the data quality associated with each sensor needs to be assessed. Data quality encompasses the reliability, accuracy, and uncertainty associated with sensor-generated data. These data need to be of known quality so that researchers, managers, and end-users can reliably determine if the information is appropriate for their intended applications. Many of the data quality recommendations from HARRNESS (2005) are still relevant and have been reiterated within the management community. Understanding and documenting data quality, and when applicable, standardizing best practices for sensor use, continue to be recommended.

Report↗

Comparison of hydrologic data and water budgets between 2003–08 and 2018–23 for the eastern part of the Arbuckle-Simpson aquifer, south-central Oklahoma

The Arbuckle-Simpson aquifer is divided spatially into three parts (eastern, central, and western). The largest groundwater withdrawals are from the eastern part of the Arbuckle-Simpson aquifer, which provides water to approximately 39,000 people in Ada and Sulphur, Oklahoma, and surrounding areas. The Arbuckle-Simpson aquifer, including the eastern part, is designated a sole source aquifer for its service area. Based primarily on data collected between 2003 and 2008, a series of comprehensive hydrologic studies of the Arbuckle-Simpson aquifer was published to provide the information necessary to perform groundwater-flow model simulations so that the Oklahoma Water Resources Board could determine how much water could be withdrawn from the aquifer while maintaining flow to springs and streams. As part of the Phase 1 studies, an aquifer water budget was developed from a numerical model for the period 2003–08. For this report, Phase 1 refers to the 2003–08 data collection period, although for some of the analyses, data collected prior to 2003 were used to inform model development work. Allocation of water from this aquifer was then established by the Oklahoma Water Resources Board in 2013. Additional well-spacing rules were also established by the Oklahoma Water Resources Board for sensitive sole source groundwater basins. To determine how the water budget for the eastern part of the Arbuckle-Simpson aquifer has changed over time, recently collected hydrologic data (2018–23) were compared to data collected during 2003–08. The analysis of changes in the aquifer water budget from 2003–08 to 2018–23 could help resource managers better understand changes in the overall balance of water in storage and the potential effects on streamflow, changes in groundwater levels, and the effects of different water uses in the aquifer area on available water in the eastern part of the Arbuckle-Simpson aquifer and streams overlying the eastern part of the Arbuckle-Simpson aquifer.

Oklahoma↗

Evaluating coastal landscape response to sea-level rise in the northeastern United States: approach and methods

The U.S. Geological Survey is examining effects of future sea-level rise on the coastal landscape from Maine to Virginia by producing spatially explicit, probabilistic predictions using sea-level projections, vertical land movement rates (due to isostacy), elevation data, and land-cover data. Sea-level-rise scenarios used as model inputs are generated by using multiple sources of information, including Coupled Model Intercomparison Project Phase 5 models following representative concentration pathways 4.5 and 8.5 in the Intergovernmental Panel on Climate Change Fifth Assessment Report. A Bayesian network is used to develop a predictive coastal response model that integrates the sea-level, elevation, and land-cover data with assigned probabilities that account for interactions with coastal geomorphology as well as the corresponding ecological and societal systems it supports. The effects of sea-level rise are presented as (1) level of landscape submergence and (2) coastal response type characterized as either static (that is, inundation) or dynamic (that is, landform or landscape change). Results are produced at a spatial scale of 30 meters for four decades (the 2020s, 2030s, 2050s, and 2080s). The probabilistic predictions can be applied to landscape management decisions based on sea-level-rise effects as well as on assessments of the prediction uncertainty and need for improved data or fundamental understanding. This report describes the methods used to produce predictions, including information on input datasets; the modeling approach; model outputs; data-quality-control procedures; and information on how to access the data and metadata online.

Open-File Report↗

Groundwater/surface-water interactions in the Bad River Watershed, Wisconsin

A groundwater-flow model was developed for the Bad River Watershed and surrounding area by using the U.S. Geological Survey (USGS) finite-difference code MODFLOW-NWT. The model simulates steady-state groundwater-flow and base flow in streams by using the streamflow routing (SFR) package. The objectives of this study were to: (1) develop an improved understanding of the groundwater-flow system in the Bad River Watershed at the regional scale, including the sources of water to the Bad River Band of Lake Superior Chippewa Reservation (Reservation) and groundwater/surface-water interactions; (2) provide a quantitative platform for evaluating future impacts to the watershed, which can be used as a starting point for more detailed investigations at the local scale; and (3) identify areas where more data are needed. This report describes the construction and calibration of the groundwater-flow model that was subsequently used for analyzing potential locations for the collection of additional field data, including new observations of water-table elevation for refining the conceptualization and corresponding numerical model of the hydrogeologic system. The study area can be conceptually divided into three primary hydrogeologic environments. The first encompasses the southern uplands with relatively low topographic relief, where groundwater-flow is unconfined and occurs primarily in sandy till and glacial outwash overlying Archean-aged crystalline bedrock. The second includes a transitional area of higher topographic relief and shallow depth to bedrock, in the vicinity of ridges formed by steeply dipping, early-Proterozoic aged metasedimentary units of the Marquette Range Supergroup (including the Ironwood Formation), and late-Proterozoic igneous units associated with the Midcontinent Rift System (MRS). Groundwater-flow in this area likely occurs primarily through connected networks of bedrock fractures that are not well characterized, and also in isolated pockets of Quaternary deposits. The third and last hydrogeologic environment includes lowlands along Lake Superior where a deep sandstone aquifer is confined by thick deposits of clay-rich till. Model input was compiled by using both published and unpublished data. Constant flux boundary conditions for the model perimeter were developed from a regional analytic element model described in appendix 1 of this report. Pumping from 26 high-capacity wells within the model area was included. The SFR stream network was developed from the National Hydrography Dataset (NHDPlus Version 2) and hydrography from the Wisconsin Department of Natural Resources (WDNR). Hydraulic conductivity values were determined for each model cell by interpolation from a network of pilot points, within zones representing major hydrogeologic units. Recharge to the groundwater system was estimated on a cell-by-cell basis by using the Soil Water Balance code (SWB), with gridded daily temperature and precipitation data for the period 1980–2011, and GIS coverages of soil and land-surface conditions. Estimated recharge varies considerably, following spatial patterns in the precipitation and soil hydrologic group inputs. The lowest recharge values occur in the Superior lowlands, whereas the highest values occur in the upland areas, especially those underlain by sandy soils, and in the vicinity of bedrock hills. The model was calibrated to groundwater-levels and base flows obtained from the USGS National Water Information System (NWIS) database, and groundwater-levels obtained from the WDNR and Band River Band well-construction databases. Calibration was performed via nonlinear regression by using the parameter-estimation software suite PEST. Groundwater levels and base-flow observations in the calibration dataset were well simulated by the calibrated model, with reasonable values of hydraulic conductivity. The pilot-point parameters that were most constrained by observations during model calibration coincided with the locations containing the most wells (head observations)—especially the population centers of Ashland, Mellen, and other communities along the major highway corridors. Results from the calibrated model illustrate differences in the nature of groundwater-flow within the watershed. In the southern part of the watershed, where bedrock is shallow, groundwater flow paths are relatively short, extending from local recharge areas to adjacent first and second-order streams. In contrast, laterally continuous deposits of clay-rich till covering the Superior Lowlands isolate most smaller streams from the sandstone aquifer, allowing for longer flow paths toward larger streams such as the Bad, Marengo, and White Rivers. Approximately three-quarters of all first-order stream cells were dry in the Superior Lowlands, compared to only half of first-order stream cells in the southern bedrock uplands. The model was used to delineate the groundwatershed for the Bad and Kakagon Rivers. “Groundwatershed” is defined as the area contributing groundwater discharge to one of these streams and their tributaries. The groundwatershed was found to align closely with the surface-watershed, with the most notable exception occurring along the southwestern half of Birch Hill, where surface water drains southwest towards the Potato River, and groundwater flows north and east towards Lake Superior. Similarly, the contributing area of groundwater-flow to the Reservation was delineated. Results indicate the off-Reservation groundwater contributing area to be limited in comparison to the extent of the watershed, extending southward into the highlands underlain by MRS igneous rock units, but not further into the area underlain by the Marquette Range Supergroup. Stable isotope samples were collected from 54 wells within the watershed, to investigate sources of groundwater. Oxygen-18 (δ 18 O) values lower than -13.0 per mil were documented in the sampling, and likely indicate the presence of recharge water from the last glacial period (>9,500 years old) beneath the northern portion of the Reservation, in the vicinity of Odanah, Wisconsin. Finally, a new data-worth analysis of potential new monitoring-well locations was performed by using the model. The relative worth of new measurements was evaluated based on their ability to increase confidence in model predictions of groundwater levels and base flows at 35 locations, under the condition of a proposed open-pit iron mine. Results of the new data-worth analysis, and other inputs and outputs from the Bad River model, are available through an online dynamic web mapping service at ( http://wim.usgs.gov/badriver/ ).

Wisconsin↗