Search USGSSearch

SEARCH · Search USGS

Results for “Scientific Monthly”

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 379 records · Page 21Linked to original sources

Evapotranspiration over spatially extensive plant communities in the Big Cypress National Preserve, southern Florida, 2007-2010

Evapotranspiration (ET) was quantified over plant communities within the Big Cypress National Preserve (BCNP) using the eddy covariance method for a period of 3 years from October 2007 to September 2010. Plant communities selected for study included Pine Upland, Wet Prairie, Marsh, Cypress Swamp, and Dwarf Cypress. These plant communities are spatially extensive in southern Florida, and thus, the ET measurements described herein can be applied to other humid subtropical locations such as the Everglades. The 3-year mean annual ET was about 1,000, 1,050, 1,100, 930, and 900 mm (millimeters) at the Dwarf Cypress, Wet Prairie, Cypress Swamp, Pine Upland, and Marsh sites, respectively. Spatial differences in annual ET were considerable due to the recovery of the Marsh site from extensive forest fire and drought conditions. Temporal variability in annual ET was relatively small at sites that were well watered (Dwarf Cypress, Wet Prairie, Cypress Swamp, Pine Upland) over the 3-year study. In other words, locations that were well watered appeared to have similar annual ET rates. Diurnal water-level variability was observed in response to ET and was less at flooded sites than at dry sites. For example, surface-water levels declined about 1.5 mm in response to ET at the flooded Cypress Swamp site during July 18-22, 2008 and declined about 10 mm in response to ET at the flooded Dwarf Cypress site from April 18-27, 2008. Specific yield was computed using ET estimates and diurnal water-level variability measured at the dry Pine Upland site as a check on the accuracy of the eddy covariance method. Water levels repeatedly dropped about 15 mm on average in response to ET at the Pine Upland site from April 27 to May 4, 2008. ET was about 3 mm on each of these days, resulting in a reasonable estimate for specific yield of 0.2 for the sandy soils at the Pine Upland site. Monthly ET estimates exhibited seasonal variation. ET was generally greatest between March to October when solar radiation was relatively large, and least from November to February when solar radiation was small. Monthly ET was greatest in the spring and summer at the Cypress Swamp site, reaching rates as large as 140 mm. The large ET rates at this site coincide with the most active period of cypress growth during late spring and early summer. Cypress trees begin to senesce in late summer reducing transpiration. Net radiation and available energy explained most of the variability in ET observed at all five sites. Mean annual and monthly net radiation varied among the sites in response to cloud cover and the albedo of the land surface and plant community. Net radiation was greatest at the Cypress Swamp site, averaging about 130 W/m 2 (watts per square meter) during the 3-year study. Net radiation was generally less at the Dwarf Cypress site, averaging about 115 W/m 2 over 3 years. The Dwarf Cypress site apparently has the largest albedo, which likely is due to the sparse canopy and a highly reflective, calcareous, periphyton-covered land surface. Furthermore, mean annual net radiation was least in the first year of the study, which likely was due to greater cloud cover during a relatively wet year. In contrast, net radiation was greatest in the second year of the study, which likely was due to less cloud cover during a relatively dry year.

Scientific Investigations Report

Transient calibration of a groundwater-flow model of Chimacum Creek Basin and vicinity, Jefferson County, Washington: a supplement to Scientific Investigations Report 2013-5160

A steady-state groundwater-flow model described in Scientific Investigations Report 2013-5160, ”Numerical Simulation of the Groundwater-Flow System in Chimacum Creek Basin and Vicinity, Jefferson County, Washington” was developed to evaluate potential future impacts of growth and of water-management strategies on water resources in the Chimacum Creek Basin. This supplement to that report describes the unsuccessful attempt to perform a calibration to transient conditions on the model. The modeled area is about 64 square miles on the Olympic Peninsula in northeastern Jefferson County, Washington. The geologic setting for the model area is that of unconsolidated deposits of glacial and interglacial origin typical of the Puget Sound Lowlands. The hydrogeologic units representing aquifers are Upper Aquifer (UA, roughly corresponding to recessional outwash) and Lower Aquifer (LA, roughly corresponding to advance outwash). Recharge from precipitation is the dominant source of water to the aquifer system; discharge is primarily to marine waters below sea level and to Chimacum Creek and its tributaries. The model is comprised of a grid of 245 columns and 313 rows; cells are a uniform 200 feet per side. There are six model layers, each representing one hydrogeologic unit: (1) Upper Confining unit (UC); (2) Upper Aquifer unit (UA); (3) Middle Confining unit (MC); (4) Lower Aquifer unit (LA); (5) Lower Confining unit (LC); and (6) Bedrock unit (OE). The transient simulation period (October 1994–September 2009) was divided into 180 monthly stress periods to represent temporal variations in recharge, discharge, and storage. An attempt to calibrate the model to transient conditions was unsuccessful due to instabilities stemming from oscillations in groundwater discharge to and recharge from streamflow in Chimacum Creek. The model as calibrated to transient conditions has mean residuals and standard errors of 0.06 ft ±0.45 feet for groundwater levels and 0.48 ± 0.06 cubic feet per second for flows. Although the expected seasonal trends were observed in model results, the typical observed annual variation of groundwater levels of about 2 feet was not. Streamflow at the most downstream observation point was about three times larger than simulated streamflow. Because the transient version of the model proved inherently unstable, it was not used to simulate forecast conditions for alternate hydrologic or anthropogenic changes. Adaptation of alternate stream simulation packages, such as RIV, or newer versions of MODFLOW, such as MODFLOW-NWT, could possibly assist with achieving calibration to transient conditions.

Washington

Modeling seawater intrusion along the Alabama coastline using physical and machine learning models to evaluate the effects of multiscale natural and anthropogenic stresses

Seawater intrusion threatens groundwater resources in coastal regions, including southern Baldwin County, Alabama, where the freshwater-saltwater interface dynamics remain poorly understood. To address this gap, this study uses combined physics-based and machine-learning models to quantify seawater intrusion caused by natural (storm surges) and anthropogenic (human activities) perturbations. The long short-term memory network and wavelet analysis were used to assess vertical aquifer vulnerabilities, revealing that the shallow part of the Coastal lowlands aquifer system (CL1) in the southern Baldwin County region is more susceptible to sea level rise and groundwater extraction than deeper aquifers. Based on these findings, a cross-sectional numerical model (physics approach) for the CL1 aquifer was developed to evaluate tidal and storm surge effects, using Tropical Storm Claudette (June 2021) as a case study. Results showed that tidal fluctuations had a minimal impact on the saltwater-freshwater interface location, whereas storm surges caused substantial inland movement, with effects lasting for nine months. The steady-state version of the three-dimensional (3D) physical model predicted seawater intrusion across the entire area, and convolutional neural network-based modeling further validated the model results. The 3D physical model was also applied to a smaller area to assess human impact on the saltwater interface due to two groundwater pumping scenarios (± 50% of the baseline pumping rate). Results revealed that a 50% increase in groundwater withdrawals caused seawater to advance ~ 320 m inland, whereas a 50% reduction led to a ~ 270-meter retreat. This study highlights the vulnerability of Alabama’s shallow coastal aquifers to seawater intrusion due to storm surges and human activities, and demonstrates that combining physics-based models with machine learning approaches can improve groundwater predictions, though its accuracy depends on the availability of site-specific data.

Alabama

Estimates of Sediment Load Prior to Dam Removal in the Elwha River, Clallam County, Washington

Years after the removal of the two dams on the Elwha River, the geomorphology and habitat of the lower river will be substantially influenced by the sediment load of the free-flowing river. To estimate the suspended-sediment load prior to removal of the dams, the U.S. Geological Survey collected suspended-sediment samples during water years 2006 and 2007 at streamflow-gaging stations on the Elwha River upstream of Lake Mills and downstream of Glines Canyon Dam at McDonald Bridge. At the gaging station upstream of Lake Mills, discrete samples of suspended sediment were collected over a range of streamflows including a large peak in November 2006 when suspended-sediment concentrations exceeded 7,000 milligrams per liter, the highest concentrations recorded on the river. Based on field measurements in this study and from previous years, regression equations were developed for estimating suspended-sediment and bedload discharge as a function of streamflow. Using a flow duration approach, the average total annual sediment load at the gaging station upstream of Lake Mills was estimated at 327,000 megagrams with a range of uncertainty of +57 to -34 percent (217,000-513,000 megagrams) at the 95 percent confidence level; 77 percent of the total was suspended-sediment load and 23 percent was bedload. At the McDonald Bridge gaging station, daily suspended-sediment samples were obtained using an automated pump sampler, and concentrations were combined with the record of streamflow to calculate daily, monthly, and annual suspended-sediment loads. In water year 2006, an annual suspended-sediment load of 49,300 megagrams was determined at the gaging station at McDonald Bridge, and a load of 186,000 megagrams was determined upstream at the gaging station upstream of Lake Mills. In water year 2007, the suspended-sediment load was 75,200 megagrams at McDonald Bridge and 233,000 megagrams upstream of Lake Mills. The large difference between suspended-sediment loads at both gaging stations shows the extent of sediment trapping by Lake Mills, and a trap efficiency of 0.86 was determined for the reservoir. Pre-dam-removal estimates of suspended-sediment load and sediment-discharge relations will help planners monitor geomorphic and habitat changes in the river as it reaches a dynamic equilibrium following the removal of dams.

Scientific Investigations Report

Metagenomic mapping of cyanobacteria and potential cyanotoxin producing taxa in large rivers of the United States

Cyanobacteria and cyanotoxin producing cyanobacterial blooms are a trending focus of current research. Many studies focus on bloom events in lentic environments such as lakes or ponds. Comparatively few studies have explored lotic environments and fewer still have examined the cyanobacterial communities and potential cyanotoxin producers during ambient, non-bloom conditions. Here we used a metagenomics-based approach to profile non-bloom microbial communities and cyanobacteria in 12 major U.S. rivers at multiple time points during the summer months of 2019. Our data show that U.S. rivers possess microbial communities that are taxonomically rich, yet largely consistent across geographic location and time. Within these communities, cyanobacteria often comprise significant portions and frequently include multiple species with known cyanotoxin producing strains. We further characterized these potential cyanotoxin producing taxa by deep sequencing amplicons of the microcystin E ( mcyE ) gene. We found that rivers containing the highest levels of potential cyanotoxin producing cyanobacteria consistently possess taxa with the genetic potential for cyanotoxin production and that, among these taxa, the predominant genus of origin for the mcyE gene is Microcystis . Combined, these data provide a unique perspective on cyanobacteria and potential cyanotoxin producing taxa that exist in large rivers across the U.S. and can be used to better understand the ambient conditions that may precede bloom events in lotic freshwater ecosystems.

Scientific Reports

Effects of the temporal variability of evapotranspiration on hydrologic simulation in central Florida

The transient response of a hydrologic system can be of concern to water-resource managers, because it is often extreme relatively short-lived events, such as floods or droughts, that profoundly influence the management of the resource. The water available to a hydrologic system for stream flow and aquifer recharge is determined by the difference of precipitation and evapotranspiration (ET). As such, temporal variations in precipitation and ET determine the degree of influence each has on the transient response of the hydrologic system. Meteorological, ET, and hydrologic data collected from 1993 to 2003 and spanning 1- to 3 2/3 -year periods were used to develop a hydrologic model for each of five sites in central Florida. The sensitivities of simulated water levels and flows to simple approximations of ET were quantified and the adequacy of each ET approximation was assessed. ET was approximated by computing potential ET, using the Hargreaves and Priestley-Taylor equations, and applying vegetation coefficients to adjust the potential ET values to actual ET. The Hargreaves and Priestley-Taylor ET approximations were used in the calibrated hydrologic models while leaving all other model characteristics and parameter values unchanged. Two primary factors that influence how the temporal variability of ET affects hydrologic simulation in central Florida were identified: (1) stochastic character of precipitation and ET and (2) the ability of the local hydrologic system to attenuate variability in input stresses. Differences in the stochastic character of precipitation and ET, both the central location and spread of the data, result in substantial influence of precipitation on the quantity and timing of water available to the hydrologic system and a relatively small influence of ET. The temporal variability of ET was considerably less than that of precipitation at each site over a wide range of time scales (from daily to annual). However, when precipitation and ET are of similar magnitude, small errors in ET can produce relatively large errors in available water, and accurate estimates of actual ET are more important. Local hydrologic conditions can also be an important factor influencing the hydrologic response to ET variability. Various points along a flow path in a hydrologic system respond differently to temporal variations in ET. For example, soil moisture contents in the root zone are sensitive to daily variations in ET, whereas spring flow responds to only longer term variations in ET. Both the Hargreaves and Priestley-Taylor equations for potential ET, when applied with an annually invariant monthly vegetation coefficient derived from comparison of measured ET with computed potential ET values, can be used with a hydrologic model to produce reasonable predictions of water levels and flows. Baseline-adjusted modified coefficients of efficiency for simulated water levels ranged from 0.0, indicating that water levels were simulated equally as well with approximated ET as with actual ET values, to -0.6, indicating that water levels were simulated better with actual ET values. Simulations using the Hargreaves approximation consistently yielded larger absolute and relative errors than the Priestley-Taylor approximation. However, the differences between the Hargreaves and Priestley-Taylor simulations generally were much smaller than differences between these simulations and the simulations using actual ET. This suggests that the simpler Hargreaves equation may be an adequate substitute for the more complex Priestley-Taylor equation, depending on the level of accuracy required to satisfy the particular modeling objectives.

Florida

Fena Valley Reservoir watershed and water-balance model updates and expansion of watershed modeling to southern Guam

In 2014, the U.S. Geological Survey, in cooperation with the U.S. Department of Defense’s Strategic Environmental Research and Development Program, initiated a project to evaluate the potential impacts of projected climate-change on Department of Defense installations that rely on Guam’s water resources. A major task of that project was to develop a watershed model of southern Guam and a water-balance model for the Fena Valley Reservoir. The southern Guam watershed model provides a physically based tool to estimate surface-water availability in southern Guam. The U.S. Geological Survey’s Precipitation Runoff Modeling System, PRMS-IV, was used to construct the watershed model. The PRMS-IV code simulates different parts of the hydrologic cycle based on a set of user-defined modules. The southern Guam watershed model was constructed by updating a watershed model for the Fena Valley watersheds, and expanding the modeled area to include all of southern Guam. The Fena Valley watershed model was combined with a previously developed, but recently updated and recalibrated Fena Valley Reservoir water-balance model. Two important surface-water resources for the U.S. Navy and the citizens of Guam were modeled in this study; the extended model now includes the Ugum River watershed and improves upon the previous model of the Fena Valley watersheds. Surface water from the Ugum River watershed is diverted and treated for drinking water, and the Fena Valley watersheds feed the largest surface-water reservoir on Guam. The southern Guam watershed model performed “very good,” according to the criteria of Moriasi and others (2007), in the Ugum River watershed above Talofofo Falls with monthly Nash-Sutcliffe efficiency statistic values of 0.97 for the calibration period and 0.93 for the verification period (a value of 1.0 represents perfect model fit). In the Fena Valley watershed, monthly simulated streamflow volumes from the watershed model compared reasonably well with the measured values for the gaging stations on the Almagosa, Maulap, and Imong Rivers—tributaries to the Fena Valley Reservoir—with Nash-Sutcliffe efficiency values of 0.87 or higher. The southern Guam watershed model simulated the total volume of the critical dry season (January to May) streamflow for the entire simulation period within –0.54 percent at the Almagosa River, within 6.39 percent at the Maulap River, and within 6.06 percent at the Imong River. The recalibrated water-balance model of the Fena Valley Reservoir generally simulated monthly reservoir storage volume with reasonable accuracy. For the calibration and verification periods, errors in end-of-month reservoir-storage volume ranged from 6.04 percent (284.6 acre-feet or 92.7 million gallons) to –5.70 percent (–240.8 acre-feet or –78.5 million gallons). Monthly simulation bias ranged from –0.48 percent for the calibration period to 0.87 percent for the verification period; relative error ranged from –0.60 to 0.88 percent for the calibration and verification periods, respectively. The small bias indicated that the model did not consistently overestimate or underestimate reservoir storage volume. In the entirety of southern Guam, the watershed model has a “satisfactory” to “very good” rating when simulating monthly mean streamflow for all but one of the gaged watersheds during the verification period. The southern Guam watershed model uses a more sophisticated climate-distribution scheme than the older model to make use of the sparse climate data, as well as includes updated land-cover parameters and the capability to simulate closed depression areas. The new Fena Valley Reservoir water-balance model is useful as an updated tool to forecast short-term changes in the surface-water resources of Guam. Furthermore, the now spatially complete southern Guam watershed model can be used to evaluate changes in streamflow and recharge owing to climate or land-cover changes. These are substantial improvements to the previous models of the Fena Valley watershed and Reservoir. Datasets associated with this report are available as a U.S. Geological Survey data release (Rosa and Hay, 2017; DOI:10.5066/F7HH6HV4).

Scientific Investigations Report

Relations among water levels, specific conductance, and depths of bedrock fractures in four road-salt-contaminated wells in Maine, 2007–9

Data on groundwater-level, specific conductance (a surrogate for chloride), and temperature were collected continuously from 2007 through 2009 at four bedrock wells known to be affected by road salts in an effort to determine the effects of road salting and fractures in bedrock that intersect the well at a depth below the casing on the presence of chloride in groundwater. Dissolved-oxygen data collected periodically also were used to make inferences about the interaction of fractures and groundwater flow. Borehole geophysical tools were used to determine the depths of fractures in each well that were actively contributing flow to the well, under both static and pumped conditions; sample- and measurement-depths were selected to correspond to the depths of these active fractures. Samples of water from the wells, collected at depths corresponding to active bedrock fractures, were analyzed for chloride concentration and specific conductance; from these analyses, a linear relation between chloride concentration and specific conductance was established, and continuous and periodic measurements of specific conductance were assumed to represent chloride concentration of the well water at the depth of measurement. To varying degrees, specific conductance increased in at least two of the wells during winter and spring thaws; the shallowest well, which also was closest to the road receiving salt treatment during the winter, exhibited the largest changes in specific conductance during thaws. Recharge events during summer months, long after application of road salt had ceased for the year, also produced increases in specific conductance in some of the wells, indicating that chloride which had accumulated or sequestered in the overburden was transported to the wells throughout the year. Geophysical data and periodic profiles of water quality along the length of each well’s borehole indicated that the greatest changes in water quality were associated with active fractures; in one case, high concentration of dissolved oxygen at the bottom of the well indicated the presence of a highly transmissive fracture that was in good connection with a surficial feature (stream or atmosphere). Data indicated that fractures have a substantial influence on the transport of chlorides to the subsurface; that elevated specific conductance occurred throughout the year, not just when road salts were applied; and that chloride contamination, as indicated by elevated specific conductance, may persist for years.

Maine

Analysis of regional rainfall-runoff parameters for the Lake Michigan Diversion hydrological modeling

The Lake Michigan Diversion Accounting (LMDA) system has been developed by the U.S. Army Corps of Engineers, Chicago District (USACE-Chicago) and the State of Illinois as a part of the interstate Great Lakes water regulatory program. The diverted Lake Michigan watershed is a 673-square-mile watershed that is comprised of the Chicago River and Calumet River watersheds. They originally drained into Lake Michigan, but now flow to the Mississippi River watershed via three canals constructed in the Chicago area in the early twentieth century. Approximately 393 square miles of the diverted watershed is ungaged, and the runoff from the ungaged portion of the diverted watershed has been estimated by the USACE-Chicago using the Hydrological Simulation Program-FORTRAN (HSPF) program. The accuracy of simulated runoff depends on the accuracy of the parameter set used in the HSPF program. Nine parameter sets comprised of the North Branch, Little Calumet, Des Plaines, Hickory Creek, CSSC, NIPC, 1999, CTE, and 2008 have been developed at different time periods and used by the USACE-Chicago. In this study, the U.S. Geological Survey and the USACE-Chicago collaboratively analyzed the parameter sets using nine gaged watersheds in or adjacent to the diverted watershed to assess the predictive accuracies of selected parameter sets. Six of the parameter sets, comprising North Branch, Hickory Creek, NIPC, 1999, CTE, and 2008, were applied to the nine gaged watersheds for evaluating their simulation accuracy from water years 1996 to 2011. The nine gaged watersheds were modeled by using the three LMDA land-cover types (grass, forest, and hydraulically connected imperviousness) based on the 2006 National Land Cover Database, and the latest meteorological and precipitation data consistent with the current (2014) LMDA modeling framework. Results indicate that the North Branch and Hickory Creek parameter sets, which belong to the original calibration group, attained an overall “satisfactory” rating on monthly runoff volumes based on the three performance statistics selected, but the annual and over-the-period runoff volumes were generally underestimated. Parameter sets CTE and 2008 attained a similar satisfactory rating on monthly runoff volumes but the annual and over-the-period runoff volumes were overestimated in general. Although the percent bias was improved, the CTE and 2008 parameter sets also had increased residuals in monthly runoff volumes and decreased quality of the model fit to the measured streamflows relative to the North Branch and Hickory Creek parameter sets. The NIPC and 1999 parameter sets, on the other hand, had larger percent bias and residuals in monthly runoff volumes, and underestimated the annual and over-the-period runoff volumes. Recalibration of the HSPF parameters to the updated inputs and land covers was completed on two representative watershed models selected from the nine by using a manual method (HSPEXP) and an automatic method (PEST). The objective of the recalibration was to develop a regional parameter set that improves the accuracy in runoff volume prediction for the nine study watersheds. Knowledge about flow and watershed characteristics plays a vital role for validating the calibration in both manual and automatic methods. The best performing parameter set was determined by the automatic calibration method on a two-watershed model. Applying this newly determined parameter set to the nine watersheds for runoff volume simulation resulted in “very good” ratings in five watersheds, an improvement as compared to “very good” ratings achieved for three watersheds by the North Branch parameter set.

Illinois, Indiana

Prioritization framework for ranking riverine ecosystem stressors using example sites from the Tualatin River Basin, Oregon

As human populations increase, so does their influence over the environment. Altered terrain, degraded water quality, and threatened or endangered species are all-too-common consequences of a growing anthropogenic influence on the landscape. To help manage these effects, researchers have developed new ways to characterize current environmental conditions and help resource managers seek solutions to bring affected areas back to their best attainable health. Before an ecosystem can be improved, however, its current level of ecological stress must be determined. Characterizing environmental conditions at many sites across a landscape helps managers understand the range of current conditions and prioritize where they might focus restoration and protection efforts. This report details the development of a prioritization framework to score riverine ecosystem stressors in a watershed based on example sites from the Tualatin River Basin in northwestern Oregon. The framework incorporated the most influential site-specific stressors throughout the basin built on a long history of data collection. These stressors were characterized with 13 metrics that were organized into 4 groups: hydrologic, water quality, physical habitat, and biological. Each stressor metric used readily accessible data and was translated to a score between 0 and 10. The higher the score, the healthier the site. This initial application of the framework used field observations and measurements to rank site conditions at two Tualatin River sites and four Tualatin River tributary sites. Given the versatility of this framework, it easily could be expanded to include more sites or new metrics, if necessary. Because stressors varied by season, all metrics for the tributary sites were scored separately during the wet season (November through April) and dry season (May through October). Water-quality data were available over a prolonged period; therefore, water-quality metrics were assessed by season and by decade (1990–99 compared to 2000–12) to evaluate long-term stressor trends. Results for the Tualatin River Basin prioritization framework indicated that the urban tributaries demonstrated the greatest stress throughout the year, especially during the dry summer months. Spatially, the upper Tualatin River was healthier than the lower reaches of the river. Water-quality has improved in the last 10 years, mostly due to improvements in the dry period contaminant scores, but challenges remain with high water temperatures and low dissolved-oxygen conditions. The biggest challenge with this type of research derived from inconsistencies within the available data. Both spatial and temporal data gaps must be addressed to improve the prioritization. Incorporating both discrete and continuous datasets into the prioritization framework remains a challenge because the datasets have slightly different information and criteria and are not always comparable. Regardless, this report provides guidelines for developing a prioritization framework that ranks the ecological health of sites in a watershed and provides guidance on management actions for improving conditions by targeting factors that greatly affect the health of river ecosystems.

Oregon

Development of regression equations to estimate flow durations, low-flow frequencies, and mean flows at ungaged stream sites in Connecticut using data through water year 2022

To aid Federal and State regulatory agencies in the effective management of water resources, the U.S. Geological Survey, in cooperation with the Connecticut Department of Energy and Environmental Protection and the Connecticut Department of Transportation, updated flow statistics for 118 streamgages and developed 47 regression equations to estimate selected flow duration, low flow, and mean flow statistics for the entire State of Connecticut, for the following: 1-, 5-, 10-, 25-, 50-, 75-, 90-, 99-percent flow durations; 7-day, 10-year low-flow frequency and 30-day, 2-year low-flow frequency; and mean flow, spring mean flow, and harmonic mean flow. In addition, regression equations were developed for monthly and seasonal flow durations, ranging from 25 to 99 percent for aquatic biological processes of salmonid spawning (November), overwinter (December–February), clupeid spawning (May), resident spawning (June), and rearing and growth (July–October) periods, and for flow durations ranging from 1 to 99 percent for the habitat forming (March–April) period. Statistics were derived from daily mean streamflow data collected from streamgages with at least 10 years of data through water year 2022 in southern New England and eastern New York. Forty streamgages in Connecticut and adjacent areas of neighboring States were used in the regression analysis. Regression methods of weighted least squares and generalized least squares were used to derive the final coefficients and measures of uncertainty for the regression equations. The equations used to estimate selected streamflow statistics were developed by relating the flow statistics to different basin characteristics (physical, land cover, and climatic) at the 40 streamgages. Nine basin characteristics served as the explanatory variables in the statewide regression equations: drainage area, percentage of area with coarse-grained stratified deposits, stream density, mean basin slope, mean basin elevation, percentage of area with hydrologic soil group A, mean monthly precipitation for November, mean seasonal precipitation in the winter (December, January, and February), and mean annual temperature. The root mean square error of the 47 equations ranged from 7.9 to 121.9 percent, with an average of 27.9 percent. The equations estimate flows most accurately near the mean (50-percent flow duration), become less accurate for low flows, and are the least accurate for extreme low flows. The root mean square error for the 50-percent flow duration is 15.1 percent, with an average of 17.6 percent across the six periods. The extreme low flow statistics of 7-day, 10-year low-flow frequency, 99-percent flow duration, and 99-percent rearing and growth period flow durations have root mean square errors of 121.9, 105.1, and 121.9 percent, respectively. The adjusted coefficient of determination of the 47 equations ranged from 73.4 to 99.5 percent, with an average of 95.1 percent.

Connecticut

Recent (2008-10) water quality in the Barton Springs segment of the Edwards aquifer and its contributing zone, central Texas, with emphasis on factors affecting nutrients and bacteria

The Barton Springs zone, which comprises the Barton Springs segment of the Edwards aquifer and the watersheds to the west that contribute to its recharge, is in south-central Texas, an area with rapid growth in population and increasing amounts of land area affected by development. During November 2008-March 2010, an investigation of factors affecting the fate and transport of nutrients and bacteria in the Barton Springs zone was conducted by the U.S. Geological Survey (USGS), in cooperation with the Texas Commission on Environmental Quality. The primary objectives of the study were to characterize occurrence of nutrients and bacteria in the Barton Springs zone under a range of flow conditions; to improve understanding of the interaction between surface-water quality and groundwater quality; and to evaluate how factors such as streamflow variability and dilution affect the fate and transport of nutrients and bacteria in the Barton Springs zone. The USGS collected and analyzed water samples from five streams (Barton, Williamson, Slaughter, Bear, and Onion Creeks), two groundwater wells (Marbridge and Buda), and the main orifice of Barton Springs in Austin, Texas. During the period of the study, during which the hydrologic conditions transitioned from exceptional drought to wetter than normal, water samples were collected routinely (every 3 to 4 weeks) from the streams, wells, and spring and, in response to storms, from the streams and spring. All samples were analyzed for major ions, nutrients, the bacterium Escherichia coli, and suspended sediment. During the dry period, the geochemistry of groundwater at the two wells and at Barton Springs was dominated by flow from the aquifer matrix and was relatively similar and unchanging at the three sites. At the onset of the wet period, when the streams began to flow, the geochemistry of groundwater samples from the Marbridge well and Barton Springs changed rapidly, and concentrations of most major ions and nutrients and densities of Escherichia coli became more similar to those of samples from the streams relative to concentrations and densities during the dry period. Geochemical modeling indicated that the proportion of Barton Springs discharge composed of stream recharge increased from about 0-8 percent during the dry period to about 80 percent during the wet period. The transition from exceptional drought to wetter-than-normal conditions resulted in a number of marked changes that highlight factors affecting the fate and transport of nutrients and bacteria and the strong influence of stream recharge on water quality in the Barton Springs segment of the Edwards aquifer and had a pronounced effect on the fate of nitrogen species. Organic nitrogen loaded to and stored in soils during the dry period was nitrified to nitrate when the soils were rewetted, resulting in elevated concentrations of nitrate plus nitrite in streams as these constituents were progressively leached during continued wet weather. Estimated mean monthly loads of organic nitrogen and nitrate plus nitrite in stream recharge and Barton Springs discharge, which were relatively low and constant during the dry period, increased during the wet period. Loads of organic nitrogen, on average, were about six times greater in stream recharge than in Barton Springs discharge, indicating that organic nitrogen likely was being converted to nitrate within the aquifer. Loads of total nitrogen (organic nitrogen plus ammonia and nitrate plus nitrite) in stream recharge (162 kilograms per day) and in Barton Springs discharge (157 kilograms per day) for the period of the investigation were not significantly different. Dilution was not an important factor affecting concentrations of nitrate plus nitrite in the streams or in Barton Springs during the period of this investigation: Concentrations of nitrate plus nitrite did not decrease in streams with increasing stream discharge, and nitrate plus nitrite concentrations measured at Barton

Scientific Investigations Report

A comparison of methods to predict historical daily streamflow time series in the southeastern United States

Effective and responsible management of water resources relies on a thorough understanding of the quantity and quality of available water. Streamgages cannot be installed at every location where streamflow information is needed. As part of its National Water Census, the U.S. Geological Survey is planning to provide streamflow predictions for ungaged locations. In order to predict streamflow at a useful spatial and temporal resolution throughout the Nation, efficient methods need to be selected. This report examines several methods used for streamflow prediction in ungaged basins to determine the best methods for regional and national implementation. A pilot area in the southeastern United States was selected to apply 19 different streamflow prediction methods and evaluate each method by a wide set of performance metrics. Through these comparisons, two methods emerged as the most generally accurate streamflow prediction methods: the nearest-neighbor implementations of nonlinear spatial interpolation using flow duration curves (NN-QPPQ) and standardizing logarithms of streamflow by monthly means and standard deviations (NN-SMS12L). It was nearly impossible to distinguish between these two methods in terms of performance. Furthermore, neither of these methods requires significantly more parameterization in order to be applied: NN-SMS12L requires 24 regional regressions—12 for monthly means and 12 for monthly standard deviations. NN-QPPQ, in the application described in this study, required 27 regressions of particular quantiles along the flow duration curve. Despite this finding, the results suggest that an optimal streamflow prediction method depends on the intended application. Some methods are stronger overall, while some methods may be better at predicting particular statistics. The methods of analysis presented here reflect a possible framework for continued analysis and comprehensive multiple comparisons of methods of prediction in ungaged basins (PUB). Additional metrics of comparison can easily be incorporated into this type of analysis. By considering such a multifaceted approach, the top-performing models can easily be identified and considered for further research. The top-performing models can then provide a basis for future applications and explorations by scientists, engineers, managers, and practitioners to suit their own needs.

Scientific Investigations Report

Groundwater response to managed aquifer recharge at the Southeast Houghton Artificial Recharge Project in Tucson, Arizona

Managed aquifer recharge is a widespread practice for storing water in the subsurface as groundwater. At a managed aquifer recharge facility in southern Arizona, groundwater-level and repeat microgravity data were collected to monitor aquifer response. These data were used to inform parameter identification for an unsaturated-zone flow model used to simulate the recharge process. The facility, the Southeast Houghton Artificial Recharge Project (SHARP), consists of 3 surface basins (about 27,600 square meters [6.8 acres] total surface area) where recycled water is distributed in recharge cycles lasting several months, with dry periods in between. During the study period, December 2020–December 2022, Tucson Water (the City of Tucson’s water utility) reported 6.56×10 6 cubic meters of water (5,320 acre-feet) recharged. Monitoring included groundwater-level observations at 3 monitoring wells and repeat microgravity measurements at as many as 22 locations (some stations were destroyed between surveys). Six gravity surveys were carried out using absolute- and relative-gravity meters. Large gravity increases, more than 250 microgals, were observed during the first repeat survey, 3.5 months after the start of recharge, but only in the immediate vicinity of the recharge basins. Data show that water moved downward to the water table, and storage changes in the unsaturated zone away from the facility were likely minimal. Gravity decreased at stations more than 1 kilometer from the facility, consistent with regional groundwater-level changes. Groundwater-level increases in wells adjacent to the recharge basins began 2 months after the second repeat gravity survey, and 5.5 months after recharge began. Unsaturated-zone flow modeling was carried out using software that simulates water movement and parameter estimation. Model calibration was carried out by minimizing an objective function calculated from the differences between simulated and observed groundwater levels, and between simulated and observed repeat microgravity data. Including repeat microgravity data in the objective function reduced the uncertainty in estimated parameter values for saturated hydraulic conductivity and saturated water content. Modeling indicated that the unsaturated zone between the recharge basins and the water table does not become saturated even after 685 days of simulated infiltration. This gradual wetting may account for increasing infiltration rates over time, as hydraulic conductivity increases with increasing water content. Unsaturated-zone water content decreased rapidly between recharge cycles. Model-simulated groundwater mounding extended about 1 kilometer from the center of SHARP after the 685-day period following the onset of recharge.

Arizona

Effects of Canals and Roads on Hydrologic Conditions and Health of Atlantic White Cedar at Emily and Richardson Preyer Buckridge Coastal Reserve, North Carolina, 2003-2006

The effects of canals and roads on hydrologic conditions and on the health of Atlantic white cedar at the Emily and Richardson Preyer Buckridge Coastal Reserve in North Carolina were evaluated by using data collected from the 1980s to 2006. Water levels were monitored along two transects established perpendicular to roads and canals in areas of healthy and unhealthy Atlantic white cedar as part of a study conducted from February 2003 through March 2006. Because of the low hydraulic gradient at the Reserve, the rate and direction of water movement are sensitive to disturbance. Canals increased drainage and contributed to lower water levels in some parts of the Reserve, whereas roads, depending on orientation, impeded drainage. Canals also appeared to facilitate movement of brackish water from the Alligator River into the interior of the Reserve during storms and wind tides. Data indicate that an influx of brackish water occurred in mid-September 2005 several days after the passage of Hurricane Ophelia. Although precipitation amounts and wind speeds associated with Hurricane Ophelia were not large, substantial changes in specific conductance occurred at the canal site on the unhealthy Atlantic white cedar transect. No corresponding increase in specific conductance was observed at the canal site on the healthy Atlantic white cedar transect. The specific conductance of water samples from canals and piezometers was highly correlated with concentrations of chloride and sodium. Ion ratios of some of the water samples, particularly samples with high specific conductance, were similar to those of seawater. Thermal and chemical stratification of water in the canals occurred during summer and winter months, and turnover and mixing occurred in the spring and fall. Upwelling of ground water as a result of excavation for roads did not appear to have a significant effect on the water quality of samples from the canals or piezometers. The specific conductance of water samples from piezometers installed in the root zone of healthy stands of Atlantic white cedar generally was lower than in water samples from unhealthy stands. This pattern also was observed in samples from piezometers installed on the transects and in other areas of the Reserve. Roads appear to have isolated some areas of the Reserve from the high-conductivity water in nearby canals. The paths by which brackish water entered the Reserve cannot be determined from the data obtained during this investigation. It appears that water can enter the Reserve from various directions, depending on wind patterns and water levels in the Alligator River.

North Carolina

Suspended-sediment characteristics for the Johnson Creek basin, Oregon, water years 2007-10

Significant Findings An analysis of suspended-sediment transport in the Johnson Creek basin, Oregon, during water years 2007–10 indicated that: Streamflow characteristics for the 4 years of study were not extremely dry or wet, and represented near-average conditions. Computed average annual suspended-sediment loads were 1,890 and 4,640 tons at the Gresham and Milwaukie stations, respectively. More than 70 percent of suspended-sediment transport in the watershed occurred during the high-flow months of November, December, and January. Less than 10 percent of suspended-sediment transport in the watershed occurred during April–October. About 50 percent of all suspended-sediment load is transported during the highest 1 percent of streamflows. The January 2009 streamflow peak was the third highest in the 70-year record for Johnson Creek. About 50 percent of suspended-sediment transport in water year 2009 occurred in January. The drainage area upstream of the Gresham streamflow-gaging station constitutes about 30 percent of the drainage area at the Milwaukie station, but accounted for about 40 percent of the suspended sediment and 45 percent of the streamflow at the Milwaukie station. On an annual basis, most of the higher sediment yield at the Gresham station, relative to the Milwaukie station, can be explained by the higher streamflow yield at the Gresham station rather than by higher suspended-sediment concentration.

Oregon

The 2006 eruption of Augustine Volcano, Alaska

Augustine Volcano, the most historically active volcano in Alaska’s Cook Inlet region, again showed signs of life in April 2005. Escalating seismic unrest, ground deformation, and gas emissions culminated in an eruption from January 11 to mid-March of 2006, the fifth major eruption in 75 years. The eruption began with a series of 13 short-lived blasts over 20 days that sent pyroclastic flows; snow, rock, and ice avalanches; and lahars down the volcano’s snow clad flanks; ash clouds drifted hundreds of kilometers downwind. Punctuated explosive activity gave way to effusion of lava and emplacement of thick block-and-ash flows on the volcano’s north flank that continued through mid-February. In mid-March renewed extrusion resulted in the building of a new, higher summit lava dome and two blocky lava flows on the north and northeast flanks of the cone. The eruption resulted in ash fall on many south-central Alaskan communities and disrupted air traffic in the region. Augustine’s frequent eruptions and relatively easy access have long drawn volcanologists to study the accumulation, ascent, and eruption of andesitic to dacitic magma. Studies of the most recent activity before 2006, in 1976 and 1986, revealed that the volcano lately produces explosive eruptions that are preceded by months of unrest and injection of new magma into a storage region in the upper several kilometers of the crust. Each of these eruptions then followed a similar progression from explosive to effusive behavior over several months. Petrologic and geophysical observations suggest that these three eruptions were triggered by similar magma mixing events and that the subsequent ascent and eruption of magma was governed by processes that were roughly constant from one eruption to the next. Geologic studies of the island show that in the more distant past parts of Augustine’s edifice have failed repeatedly, resulting in debris avalanches that entered the sea and, at least once, in 1883, caused a tsunami that hit surrounding Cook Inlet coastlines. Such edifice failures and resultant local tsunamis should be expected in the future. Recognition of Augustine’s frequent activity and hazardous nature led to the installation of a network of telemetered seismometers beginning in 1971, the establishment of a geodetic network in 1988, and the installation of other new instrumentation such as pressure sensors, broadband seismometers, and cameras by the Alaska Volcano Observatory (AVO), and the selection of Augustine for geodetic instrumentation through the EarthScope/Plate Boundary Observatory program in 2004. In addition, remote sensing techniques, such as airborne thermal imaging and the advanced spaceborne thermal emission and reflection radiometer (ASTER), provided novel and often critical information as the 2006 eruption progressed. The combination of a long-term seismic network and an array of new monitoring techniques has provided a breadth and depth of understanding of Augustine’s most recent activity that was not possible in the past. This volume contains 28 chapters reporting on a diverse suite of new scientific observations and investigations that were motivated by the 2006 eruption. Understanding the magmatic processes that drive eruptions, identifying eruptive events, tracking the movement of ash clouds, and communicating the resultant hazards to other government agencies and the public are all critical tasks for AVO, and chapters touch upon all of these topics. One goal in this compilation is to synthesize the diverse information into as complete an understanding of the magmatic and eruptive processes as possible. An equally important goal is to provide a framework for diagnosing periods of unrest and formulating forecasts of eruptions that will certainly take place at Augustine in the future. This latter goal is especially important, as Augustine’s frequent eruptive activity suggests that another eruption can be expected within the next several decades. Consequently, the investigations in this volume are intended to provide both a means to better forecast future eruptive episodes and also an opportunity to formulate and test future hypotheses for magmatic and eruptive processes. Future eruptions may follow a course similar to those observed in 1976, 1986, and 2006. However, a major perturbation that upsets conditions within the magmatic system could occur, owing perhaps to the rise of a much larger or different parental magma or to a large edifice failure similar to the 1883 sector collapse. In such events, the comprehensive study of past eruptions will provide data critical to assessing the current state of the magmatic system. In assembling this volume we have sought as consistent and accurate a portrayal of the 2006 eruption as possible. We have asked all authors to refer to the same basic eruption chronology, unless their observations and data require alternative explanations. Naturally, not all techniques or methodologies produce a completely consistent set of observations, nor do the precise conclusions in every paper support one another. We have grouped chapters on the basis of discipline. Papers that focus on specific techniques, methodology, or instrumentation are placed throughout the volume where they best fit with others that rely on their results.

Alaska

Invasive brown treesnakes (Boiga irregularis) move short distances and have small activity areas in a high prey environment

Animal movements reflect temporal and spatial availability of resources as well as when, where, and how individuals access such resources. To test these relationships for a predatory reptile, we quantified the effects of prey abundance on the spatial ecology of invasive brown treesnakes ( Boiga irregularis ) on Guam. Five months after toxicant-mediated suppression of a brown treesnake population, we simultaneously used visual encounter surveys to generate relative rodent abundance and radiotelemetry of snakes to document movements of surviving snakes. After snake suppression, encounter rates for small mammals increased 22-fold and brown treesnakes had smaller mean daily movement distances (24 ± 13 m/day, x ¯ "> x ± SD) and activity areas (5.47 ± 5 ha) than all previous observations. Additionally, snakes frequenting forest edges, where our small mammal encounters were the highest, had smaller mean daily movement distances and three-dimensional activity volumes compared to those within the forest interior. Collectively, these results suggest that reduced movements by snakes were in part a response to increased prey availability. The impact of prey availability on snake movement may be a management consideration when attempting to control cryptic invasive species using tools that rely on movement of the target species to be effective,

Guam