Search USGSSearch

SEARCH · Search USGS

Results for “Scientific Investigations 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 1,189 records · Page 66Linked to original sources

Methods to Identify Changes in Background Water-Quality Conditions Using Dissolved-Solids Concentrations and Loads as Indicators, Arkansas River and Fountain Creek, in the Vicinity of Pueblo, Colorado

Effective management of existing water-storage capacity in the Arkansas River Basin is anticipated to help satisfy the need for water in southeastern Colorado. A strategy to meet these needs has been developed, but implementation could affect the water quality of the Arkansas River and Fountain Creek in the vicinity of Pueblo, Colorado. Because no known methods are available to determine what effects future changes in operations will have on water quality, the U.S. Geological Survey, in cooperation with the Southeastern Colorado Water Activity Enterprise, began a study in 2002 to develop methods that could identify if future water-quality conditions have changed significantly from background (preexisting) water-quality conditions. A method was developed to identify when significant departures from background (preexisting) water-quality conditions occur in the lower Arkansas River and Fountain Creek in the vicinity of Pueblo, Colorado. Additionally, the methods described in this report provide information that can be used by various water-resource agencies for an internet-based decision-support tool. Estimated dissolved-solids concentrations at five sites in the study area were evaluated to designate historical background conditions and to calculate tolerance limits used to identify statistical departures from background conditions. This method provided a tool that could be applied with defined statistical probabilities associated with specific tolerance limits. Drought data from 2002 were used to test the method. Dissolved-solids concentrations exceeded the tolerance limits at all four sites on the Arkansas River at some point during 2002. The number of exceedances was particularly evident when streamflow from Pueblo Reservoir was reduced, and return flows and ground-water influences to the river were more prevalent. No exceedances were observed at the site on Fountain Creek. These comparisons illustrated the need to adjust the concentration data to account for varying streamflow. As such, similar comparisons between flow-adjusted data were done. At the site Arkansas River near Avondale, nearly all the 2002 flow-adjusted concentration data were less than the flow-adjusted tolerance limit which illustrated the effects of using flow-adjusted concentrations. Numerous exceedances of the flow-adjusted tolerance limits, however, were observed at the sites Arkansas River above Pueblo and Arkansas River at Pueblo. These results indicated that the method was able to identify a change in the ratio of source waters under drought conditions. Additionally, tolerance limits were calculated for daily dissolved-solids load and evaluated in a similar manner. Several other mass-load approaches were presented to help identify long-term changes in water quality. These included comparisons of cumulative mass load at selected sites and comparisons of mass load contributed at the Arkansas River near Avondale site by measured and unmeasured sources.

Scientific Investigations Report

Regionalized equations for bankfull-discharge and channel characteristics of streams in New York State—Hydrologic Region 6 in the Southern Tier of New York

Equations that relate bankfull discharge and channel characteristics (width, depth, and cross-sectional area) to drainage-area size at gaged sites are needed to define bankfull discharge and channel dimensions at ungaged sites and to provide information for watershed assessments, stream-channel classification, and the design of stream-restoration projects. Such equations are most accurate if derived from streams within an area of uniform hydrologic, climatic, and physiographic conditions and applied only within that region. In New York State, eight hydrologic regions were previously defined on the basis of similar high-flow (flood) characteristics. This report presents drainage areas and associated bankfull characteristics (discharge and channel dimensions) for surveyed streams in southwestern New York (Region 6). Stream-survey data and discharge records from 11 active (currently gaged) sites and 3 inactive (discontinued) sites were used in regression analyses to relate bankfull discharge and bankfull channel width, depth, and cross-sectional area to the size of the drainage area. The resulting equations are: (1) bankfull discharge, in cubic feet per second = 48.0*(drainage area, in square miles) 0.842 ; (2) bankfull channel width, in feet = 16.9*(drainage area, in square miles) 0.419 ; (3) bankfull channel depth, in feet = 1.04*(drainage area, in square miles) 0.244 ; and (4) bankfull channel cross-sectional area, in square feet = 17.6*(drainage area, in square miles) 0.662 . The coefficient of determination (R 2 ) for these four equations were 0.90, 0.79, 0.64, and 0.89, respectively. The high correlation coefficients for bankfull discharge and cross-sectional area indicate that much of the variation in these variables is explained by the size of the drainage area. The smaller correlation coefficients for bankfull channel width and depth indicate that other factors also affect these relations. Recurrence intervals for the estimated bankfull discharge of each stream ranged from 1.01 to 2.35 years; the mean recurrence interval was 1.54 years. The 14 surveyed streams were classified by Rosgen stream type; most were C-type reaches, with occasional B-type reaches. The Region 6 equation (curve) for bankfull discharge was compared with equations previously developed for four other large areas in New York State and southeastern Pennsylvania. The differences among results indicate that, although the equations need to be refined by region before being applied by water-resources managers to local planning and design efforts, similar regions have similar relations between bankfull discharge and channel characteristics.

Scientific Investigations Report

Peak streamflow trends in Minnesota and their relation to changes in climate, water years 1921–2020

This report chapter summarizes the effect of hydroclimatic variability of annual peak streamflow in Minnesota and is part of a larger U.S. Geological Survey multistate study to assess potential nonstationarity in annual peak streamflows across the Midwest. Spatial and temporal patterns were examined for nonstationarity in annual peak streamflow, daily mean streamflow, and modeled climatic data in four analysis periods: (1) a 100-year period, 1921–2020; (2) a 75-year period, 1946–2020; (3) a 50-year period, 1971–2020; and (4) a 30-year period, 1991–2020. Upward trends in annual peak streamflow were detected in northwest to southeast and north to south directions. Downward trends in annual peak streamflow were detected in northeastern and southeastern areas. Trends in peak-flow timing indicated that peak streamflows are being detected later in the water year (the period from October 1 to September 30 designated by the year in which it ends) mainly in the southern areas and earlier in the water year mainly in the northern areas. Changes in climate data point to wetter conditions in southern areas and drier conditions in northern areas. Annual precipitation was determined to be increasing in a northwest to southeast direction and in the east. In contrast, some areas in the north and northwest indicated decreasing annual precipitation. Annual snowfall was determined to be decreasing except in the extreme northeast, where annual snowfall was determined to be increasing. Decreases in annual potential evapotranspiration were detected in the south, and increases were detected in the north. Annual soil moisture increased in southern areas and decreased in northern and eastern areas. The potential spatial and temporal nonstationarity violations detected in the four analysis periods have important implications for flood-frequency analysis and point to the need for guidance on how to incorporate nonstationarities into future flood-frequency analysis in Minnesota.

Minnesota

American white pelicans breeding in the northern plains: productivity, behavior, movements, and migration

Nearly half of American white pelicans (Pelecanus erythrorhynchos; hereafter pelicans) are believed to nest in several large colonies in the northern plains, yet few studies had been conducted on pelicans in this region until research began in 2004 to investigate the impact of West Nile virus (WNV) on their chicks. The work reported here focused on two of the largest colonies in the region, at Bitter Lake, South Dakota, and Chase Lake, North Dakota, during 2005–10. Pelicans usually began arriving at these two breeding colonies in early April. Egg-laying began during mid-April and nest initiations continued through May. The number of nests documented at these colonies reached a high of about 15,400 at Bitter Lake and 17,300 at Chase Lake, both in 2006. During 2006–8, annual variation in hatching success was high (40 to 100 percent) at video-monitored nests, averaging 61 percent for 82 nests at Bitter Lake and 57 percent for 88 nests at Chase Lake. Although most nests contained two eggs, of those where two chicks hatched, both chicks survived to the crèche stage (about 15 days old) at only two nests. Severe weather events, disturbance, and siblicide were documented causes of early season (before mid-July) mortalities. In the late season (mid-July to fledging), WNV was the most important factor contributing to chick mortality. Nests were nearly always attended by one adult during incubation and brooding. Adults typically exchanged places at the nest around mid-day in all 3 years, apparently taking advantage of thermals to gain altitude for travel to and from foraging areas. The mean time of exchange differed by about an hour between Bitter Lake (1328 central standard time, CST) and Chase Lake (1434 CST) colonies. During incubation, nearly 3 days often passed between adult nest exchanges; after chicks hatched, exchanges usually occurred daily. Exchanges were more frequent and chicks were fed more often at successful nests than at failed nests. Adult pelicans with satellite transmitters that incorporated a Global Positioning System (GPS) foraged primarily in shallow areas of lakes and semipermanent wetlands. These areas coincide with typical habitats of crayfish, salamanders, and rough fish, which were also the foods most commonly seen in pelican regurgitates at the colonies. Several satellite-tracked pelicans made frequent round trips between their breeding colony and foraging areas, most likely to provision their chicks. Typical distances travelled to foraging sites ranged from 30 kilometers to over 90 kilometers. Return times to the colonies (about 1300 and 1500 CST at Bitter Lake and Chase Lake, respectively) supported the colony difference documented at video-monitored nests. Of 28 pelicans tagged with GPS satellite transmitters in 2005–6, 26 survived the first summer and migrated south during fall. Nineteen of these returned to the breeding region (defined as north of the latitude of South Dakota’s southern border) in at least 1 year during 2006–9; collectively, they returned to the breeding region 33 times. Very few pelicans returned to the colony where they had been tagged; many did not breed and concentrated their activities at wetland complexes in South Dakota and North Dakota, but few tagged pelicans temporally overlapped at specific sites. During 2005–9, tagged pelicans collectively made 56 migratory trips south in the fall. Most wintered in Mexico, near the gulf coast and elsewhere; others wintered in Texas, Louisiana, Mississippi, and Florida. Individuals typically returned to the same general areas each winter. Individuals rarely followed the same migratory path on their way south and north, but they often roughly repeated southerly or northerly routes among years. Ensuring a sustainable population of American white pelicans requires identification and mitigation of known threats. The work described herein has identified WNV and severe weather as important factors that potentially limit reproductive success and recruitment in the northern plains. Managers in this region could assess the influence of such factors on productivity at key colonies by annually obtaining aerial photographs during peak nesting, and estimating numbers of chicks fledged from aerial photographs or ground counts. Banding a subsample of chicks in late June or early July, followed by a sweep for bands at the end of the season, would allow estimation of mortality rates of older chicks (that normally would fledge) and help track the influence of WNV or other mortality factors over time and varying environmental conditions.

North Dakota;South Dakota

Peak streamflow trends in North Dakota and their relation to changes in climate, water years 1921–2020

Standardized guidelines for completing flood-flow frequency analyses are presented in a U.S. Geological Survey Techniques and Methods report known as Bulletin 17C, https://doi.org/10.3133/tm4B5 . In recent decades (since about 2000), a better understanding of long-term climatic persistence (periods of clustered floods or droughts, or wet or dry periods) and concerns about potential climate change and land-use change have caused a reexamination of the stationarity assumptions underlying methods in Bulletin 17C. Bulletin 17C does not offer guidance on incorporating nonstationarities and further identifies a need for flood-frequency studies that incorporate changing climate or basin characteristics. As part of that reexamination, a study of annual peak streamflow (peak flow) has begun in the Midwest. This chapter of the study summarizes how hydroclimatic variability affects peak flows in North Dakota. In this analysis of peak flow, daily streamflow, and climate metrics, four periods were selected: (1) a 100-year period, 1921–2020; (2) a 75-year period, 1946–2020; (3) a 50-year period, 1971–2020; and (4) a 30-year period, 1991–2020. Output from a monthly water-balance model was used for the climate data. Statistical analysis of peak flow consisted of evaluations of autocorrelation, trends, and change points and was augmented with analyses of seasonality and daily streamflow. The long-term pattern of decreasing peak flow in the west and increasing peak flow in the east is a pattern of opposing signals on either side of the 100th meridian. Analyses indicate that a key factor in changing hydroclimatology is the increase in fall precipitation. The trends in soil moisture closely match the trends in annual precipitation. Nonstationary flood-frequency analysis necessitates detailed exploratory data analysis and additional data and information about climate, land use, and other factors. This study provides extensive exploratory analysis for peak flow, daily streamflow, and climate data for North Dakota, setting the stage for informed nonstationary flood-frequency analysis.

North Dakota

Stochastic model for simulating Souris River Basin regulated streamflow upstream from Minot, North Dakota

The Souris River Basin is a 24,000 square-mile basin in the Provinces of Saskatchewan and Manitoba in Canada, and the State of North Dakota in the United States. Above-average snowpack during the winter of 2010–11, along with record-setting rains in May and June of 2011, led to record flooding that caused extensive damage to Minot, North Dakota, and numerous smaller communities in Saskatchewan, Manitoba, and North Dakota. As a result, the International Souris River Board created the Souris River Flood Task Force to evaluate potential reservoir operation changes and flood control measures to manage future floods and droughts. Part of this evaluation involved identifying a need for a stochastic streamflow model to estimate the likelihood of future flooding or drought. A stochastic natural (unregulated) streamflow simulation model described in a previous report was built upon in this report to include the effects of regulation of four reservoirs (Rafferty, Alameda, and Boundary Reservoirs and Lake Darling) and their operation guidelines. First, a regulated reservoir storage/streamflow routing model was developed and calibrated from when all four reservoirs were in operation until the end of the reconstructed natural streamflow dataset provided by the U.S. Army Corps of Engineers (1992–2011). The regulated reservoir storage/streamflow routing model then was combined with the stochastic natural (unregulated) streamflow model to provide a stochastic regulated streamflow simulation model for the Souris River Basin upstream from Minot, North Dakota. The stochastic regulated streamflow simulation model was used to estimate regulated flood frequency curves, which are useful for feasibility and design of critical structures such as levees or bridges. Three potential future climatic conditions were considered in this analysis: condition A (wet equilibrium), representing wet (similar to 1970–2017) climatic conditions; condition B (transition), representing transition from wet to dry (similar to 1912–69) climatic conditions; and condition C (dry equilibrium), representing dry climatic conditions. Comparison of the estimated flood frequency curves for regulated flow among the three climatic conditions indicated large differences in flood magnitudes for the more extreme (1-percent or less) annual exceedance probabilities. The estimated 0.2-percent annual exceedance probability flood magnitude for the Souris River upstream from Minot, N. Dak., was 29,300 cubic feet per second for condition A (wet equilibrium), compared to 14,800 cubic feet per second for condition C (dry equilibrium). For comparison, the recorded peak streamflow for 2011 for the Souris River upstream from Minot, N. Dak., was 26,900 cubic feet per second. Although it is not possible to predict how long the current (1970–2017) wet climatic conditions may persist, flood risk for at least the next 25 years, or until about 2040, may be represented best by climatic condition A.

North Dakota

Assessing inundation hazards to nuclear powerplant sites using geologically extended histories of riverine floods, tsunamis, and storm surges

Most nuclear powerplants in the United States are near rivers, large lakes, or oceans. As evident from the Fukushima Daiichi, Japan, disaster of 2011, these water bodies pose inundation threats. Geologic records can extend knowledge of rare hazards from flooding, storm surges, and tsunamis. This knowledge can aid in assessing the safety of critical structures such as dams and energy plants, for which even remotely possible hazards are pertinent. Quantitative analysis of inundation from geologic records perhaps is most developed for and applied to riverine flood hazards, but because of recent natural disasters, geologic investigations also are now used widely for understanding tsunami hazards and coastal storm surges. Layered sedimentary deposits commonly give the most complete geologic record of large floods, storm surges, and tsunamis. Sedimentary layers may be preserved for hundreds or thousands of years in suitable depositional environments, thereby providing an archive of rare, high-magnitude events. All inundation hazards discussed in this report—riverine floods, tsunamis, and storm surges—have had long records extracted from sedimentary sequences, many specifically for hazard assessment. Geologic records commonly are imprecise, so most hazard assessments benefit from evaluation of many sites and rigorous uncertainty assessment. Despite uncertainties, geologic records commonly can improve knowledge of the types and magnitudes of hazards threatening specific sites or regions. New statistical tools and approaches can efficiently incorporate geologic information into frequency assessments. These tools are most developed for riverine flood hazards, but are to some degree transferable to other episodic natural phenomena such as tsunamis and storm surges. Even with these efficient statistical approaches for examining geologic records, systematic landscape changes may reduce the applicability of retrospective assessments. These non-stationarity issues (such as climate change, sea‑level rise, land-use, dams and flow regulation) may all affect the validity of using past experience—no matter how complete the record—to assess future likelihoods. These issues require site-specific consideration for nearly all hazard assessments drawn from geologic evidence. A screening of the 104 nuclear powerplants in the United States licensed by the Nuclear Regulatory Commission (at 64 sites) indicates several sites for which paleoflood studies likely would provide additional flood-frequency information. Two sites—Duane Arnold, Iowa, on the Cedar River; and David-Besse, Ohio, on the Toussaint River—have geologic conditions suitable for creating and preserving stratigraphic records of flooding and few upstream dams that may complicate flood-frequency analysis. One site—Crystal River, Florida1, on the Withlacoochee River and only 4 kilometers from the coast—has high potential as a candidate for assessing riverine and marine inundation hazards. Several sites on the Mississippi River have high geologic potential, but upstream dams almost certainly now regulate peak flows. Nevertheless, studies on the Mississippi River to evaluate long-term flood frequency may provide results applicable to a wide spectrum of regional hazard issues. Several sites in the southeastern United States have high geologic potential, and studies at these sites also may be helpful in evaluating hazards from outburst floods from landslide dams (river blockages formed by mass movements), which may be a regional hazard. For all these sites, closer investigation and field reconnaissance would be needed to confirm suitable deposits and settings for a complete paleoflood analysis. Similar screenings may help identify high-potential sites for geologic investigations of tsunami and storm-surge hazards.

Scientific Investigations Report

Evaluation of two low-flow releases from Big Tujunga Reservoir, Los Angeles County, California, 2003

Since 1973, the Santa Ana Sucker (Catostomus santaanae) has been listed as a threatened species under the Endangered Species Act. The Lower Big Tujunga Creek, in Los Angeles County, is one of the areas in southern California where the Santa Ana Sucker is still present. This study was designed to assess two flow releases from Big Tujunga Dam that may contribute to favorable habitat conditions for the Santa Ana Sucker. It is important for the Santa Ana Sucker's survival that pools in the lower reach of the study area are replenished periodically. The focus of the study area was on the Lower Big Tujunga Creek within a reach extending approximately 6 miles downstream from the Big Tujunga Reservoir. Six sites were established from the Big Tujunga Dam to Delta Flats day-use area for data collection. This report describes the study design, discharge measurements, and the flow data collected from the two releases. Two scheduled flows (phases 1 and 2) were released from the Big Tujunga Reservoir in August and September 2003. During the first phase, which lasted 50 hours, travel times from the dam to four sites downstream were determined. Arrival times at the four sites were determined on the basis of temperature data. Travel time from the dam to site 6 (the furthest downstream site) was about 51.5 hours. Travel times for subreaches were 3 hours from site 1 to site 2, 6.5 hours from site 2 to site 3, almost 18 hours from site 3 to site 4, and 24 hours from site 4 to site 6. The temperature probe at site 5 was destroyed, and thus the arrival time could not be estimated. A probe that measures stage was placed in one of the many pools downstream from site 4 to evaluate a typical pool response to a low-flow release. Also, discharge measurements were taken at four sites along the study reach. In phase 2, which lasted 5 days (121 hours), flow losses along the 6-mile reach were analyzed. Losses were estimated by measuring difference in flow from the dam to sites 3, 4, 5, and 6, when flow was most stable at each site or when the last measurement made before flow decreased due to flow from dam being shut off. Losses in the plunge pool, directly below the dam were assumed to be negligible for this study. Overall creek loss between the dam and site 6 (the last site) was estimated to be between 4.0 and 4.2 ft3/s (cubic feet per second). Estimated losses between the dam and intermediate sites were about 1.5 ft3/s to site 3; 2.5 ft3/s to site 4; and between 3.7 and 4.1 ft3/s to site 5.

Scientific Investigations Report

Estimated Depth to Ground Water and Configuration of the Water Table in the Portland, Oregon Area

Reliable information on the configuration of the water table in the Portland metropolitan area is needed to address concerns about various water-resource issues, especially with regard to potential effects from stormwater injection systems such as UIC (underground injection control) systems that are either existing or planned. To help address these concerns, this report presents the estimated depth-to-water and water-table elevation maps for the Portland area, along with estimates of the relative uncertainty of the maps and seasonal water-table fluctuations. The method of analysis used to determine the water-table configuration in the Portland area relied on water-level data from shallow wells and surface-water features that are representative of the water table. However, the largest source of available well data is water-level measurements in reports filed by well constructors at the time of new well installation, but these data frequently were not representative of static water-level conditions. Depth-to-water measurements reported in well-construction records generally were shallower than measurements by the U.S. Geological Survey (USGS) in the same or nearby wells, although many depth-to-water measurements were substantially deeper than USGS measurements. Magnitudes of differences in depth-to-water measurements reported in well records and those measured by the USGS in the same or nearby wells ranged from -119 to 156 feet with a mean of the absolute value of the differences of 36 feet. One possible cause for the differences is that water levels in many wells reported in well records were not at equilibrium at the time of measurement. As a result, the analysis of the water-table configuration relied on water levels measured during the current study or used in previous USGS investigations in the Portland area. Because of the scarcity of well data in some areas, the locations of select surface-water features including major rivers, streams, lakes, wetlands, and springs representative of where the water table is at land surface were used to augment the analysis. Ground-water and surface-water data were combined for use in interpolation of the water-table configuration. Interpolation of the two representations typically used to define water-table position - depth to the water table below land surface and elevation of the water table above a datum - can produce substantially different results and may represent the end members of a spectrum of possible interpolations largely determined by the quantity of recharge and the hydraulic properties of the aquifer. Datasets of depth-to-water and water-table elevation for the current study were interpolated independently based on kriging as the method of interpolation with parameters determined through the use of semivariograms developed individually for each dataset. Resulting interpolations were then combined to create a single, averaged representation of the water-table configuration. Kriging analysis also was used to develop a map of relative uncertainty associated with the values of the water-table position. Accuracy of the depth-to-water and water-table elevation maps is dependent on various factors and assumptions pertaining to the data, the method of interpolation, and the hydrogeologic conditions of the surficial aquifers in the study area. Although the water-table configuration maps generally are representative of the conditions in the study area, the actual position of the water-table may differ from the estimated position at site-specific locations, and short-term, seasonal, and long-term variations in the differences also can be expected. The relative uncertainty map addresses some but not all possible errors associated with the analysis of the water-table configuration and does not depict all sources of uncertainty. Depth to water greater than 300 feet in the Portland area is limited to parts of the Tualatin Mountains, the foothills of the Cascade Range, and muc

Scientific Investigations Report

Use of multidimensional modeling to evaluate a channel restoration design for the Kootenai River, Idaho

River channel construction projects aimed at restoring or improving degraded waterways have become common but have been variously successful. In this report a methodology is proposed to evaluate channel designs before channels are built by using multidimensional modeling and analysis. This approach allows detailed analysis of water-surface profiles, sediment transport, and aquatic habitat that may result if the design is implemented. The method presented here addresses the need to model a range of potential stream-discharge and channel-roughness conditions to best assess the function of the design channel for a suite of possible conditions. This methodology is demonstrated by using a preliminary channel-restoration design proposed for a part of the Kootenai River in northern Idaho designated as critical habitat for the endangered white sturgeon (Acipenser transmontanus) and evaluating the design on the basis of simulations with the Flow and Sediment Transport with Morphologic Evolution of Channels (FaSTMECH) model. This evaluation indicated substantial problems with the preliminary design because boundary conditions used in the design were inconsistent with best estimates of future conditions. As a result, simulated water-surface levels did not meet target levels that corresponded to the designed bankfull surfaces; therefore, the flood plain would not function as intended. Sediment-transport analyses indicated that both the current channel of the Kootenai River and the design channel are largely unable to move the bed material through the reach at bankfull discharge. Therefore, sediment delivered to the design channel would likely be deposited within the reach instead of passing through it as planned. Consequently, the design channel geometry would adjust through time. Despite these issues, the design channel would provide more aquatic habitat suitable for spawning white sturgeon (Acipenser transmontanus) at lower discharges than is currently available in the Kootenai River. The evaluation methodology identified potential problems with the design channel that can be addressed through design modifications to better meet project objectives before channel construction.

Scientific Investigations Report

Water availability and use pilot: Methods development for a regional assessment of groundwater availability, southwest alluvial basins, Arizona

Executive Summary: Arizona is located in an arid to semiarid region in the southwestern United States and is one of the fastest growing States in the country. Population in Arizona surpassed 6.5 million people in 2008, an increase of 140 percent since 1980, when the last regional U.S. Geological Survey (USGS) groundwater study was done as part of the Regional Aquifer System Analysis (RASA) program. The alluvial basins of Arizona are part of the Basin and Range Physiographic Province and cover more than 73,000 mi 2 , 65 percent of the State's total land area. More than 85 percent of the State's population resides within this area, accounting for more than 95 percent of the State's groundwater use. Groundwater supplies in the area are expected to undergo further stress as an increasing population vies with the State's important agricultural sector for access to these limited resources. To provide updated information to stakeholders addressing issues surrounding limited groundwater supplies and projected increases in groundwater use, the USGS Groundwater Resources Program instituted the Southwest Alluvial Basins Groundwater Availability and Use Pilot Program to evaluate the availability of groundwater resources in the alluvial basins of Arizona. The principal products of this evaluation of groundwater resources are updated groundwater budget information for the study area and a proof-of-concept groundwater-flow model incorporating several interconnected groundwater basins. This effort builds on previous research on the assessment and mapping of groundwater conditions in the alluvial basins of Arizona, also supported by the USGS Groundwater Resources Program. Regional Groundwater Budget: The Southwest Alluvial Basins-Regional Aquifer System Analysis (SWAB-RASA) study produced semiquantitative groundwater budgets for each of the alluvial basins in the SWAB-RASA study area. The pilot program documented in this report developed new quantitative estimates of groundwater budget components using recent (2000-2007) data and methods of data analysis. Estimates of inflow components, including mountain-front recharge, incidental recharge from irrigation of agriculture, managed recharge from recharge facilities, interbasin underflow from upgradient basins, and streamflow losses, are quantified for recent time periods. Mountain-front recharge is the greatest inflow component to the groundwater system and was estimated using two methods: a basin characteristic model and new precipitation information used in a previously developed regression equation. Annual mountain-front recharge for the study area for 1940-2007 estimated by the two methods is 730,000 acre-ft for the basin characteristic model and 643,000 acre-ft for the regression equation, representing 1.5 percent and 1.3 percent of precipitation, respectively. Outflow components, including groundwater withdrawals, evapotranspiration, and interbasin flow to downgradient basins, are also presented for recent time periods. Groundwater withdrawals accounted for the largest share of the water budget, with nearly 2.4 million acre-ft per year withdrawn from the study area in recent years. Evapotranspiration from groundwater was estimated at nearly 1.3 million acre-ft per year for the study area using a newly developed method incorporating vegetation indices from satellite images and land cover information. For water-budget components with temporal variation that could be assessed from available data, estimates for intervening time periods since before development were also developed. An estimate of aquifer storage change, representing both gains to and losses from the groundwater system since before development, was derived for the most developed basins in the study area using available estimates of groundwater-level changes and storage coefficients. An overall storage loss of 74.5 million acre-ft was estimated for these basins within the study area. Demonstration

Arizona

Using tracers to evaluate streamflow gain-loss characteristics of Terror Creek, in the vicinity of a mine-permit area, Delta County, Colorado, water year 2003

In 2003, the U.S. Geological Survey, in cooperation with Delta County, initiated a study to characterize streamflow gainloss in a reach of Terror Creek, in the vicinity of a mine-permit area planned for future coal mining. This report describes the methods of the study and includes results from a comparison of two sets of streamflow measurements using tracer techniques following the constant-rate injection method. Two measurement sets were used to characterize the streamflow gain-loss associated with reservoir-supplemented streamflow conditions and with natural base-flow conditions. A comparison of the measurement sets indicates that the streamflow gain-loss characteristics of the Terror Creek study reach are consistent between the two hydrologic conditions evaluated. A substantial streamflow gain occurs between measurement locations 4 and 5 in both measurement sets, and streamflow is lost between measurement locations 5 and 7 (measurement set 1, measurement location 6 not visited) and 5 and 6 (measurement set 2). A comparison of the measurement sets above and below the mine-permit area (measurement locations 3 and 7) shows a consistent loss of 0.37 and 0.31 cubic foot per second (representing 5- and 12-percent streamflow losses normalized to measurement location 3) for measurement sets 1 and 2, respectively. This indicates that similar streamflow losses occur both during reservoir-supplemented and natural base-flow conditions, with a mean streamflow loss of 0.34 cubic foot per second for measurement sets 1 and 2. Findings from a previous investigation support the observed streamflow loss between measurement locations 3 and 7 in this study. The findings from the previous investigation indicate a streamflow loss of 0.59 cubic foot per second occurs between these measurement locations. Statistical testing of the differences in streamflow between measurement locations 3 and 7 indicates that there is a discernible streamflow loss. The p-value of 0.0236 for the parametric paired t-test indicates that there is a 2.36-percent probability of observing a sample mean difference of 0.34 cubic foot per second if the population mean is zero. The p-value of 0.125 for the nonparametric exact Wilcoxon signed rank test indicates that there is a 12.5-percent probability of observing a sample mean difference this large if the population mean is zero. The similarity in streamflow gain-loss between measurement sets indicates that the process controlling streamflow may be the same between the two hydrologic conditions evaluated. Gains between measurement locations 4 and 5 may be related to hyporheic flow from tributaries that were dry during the study. No other obvious sources of surface water were identified during the investigation. The cause for the observed streamflow loss between measurement locations 5 and 6 is unknown but may be related to mapped local faulting, 100 years of coal mining in the area, and aquifer recharge.

Scientific Investigations Report

Lava lake activity at the summit of Kīlauea Volcano in 2016

The ongoing summit eruption at Kīlauea Volcano, Hawai‘i, began in March 2008 with the formation of the Overlook crater, within Halema‘uma‘u Crater. As of late 2016, the Overlook crater contained a large, persistently active lava lake (250 × 190 meters). The accessibility of the lake allows frequent direct observations, and a robust geophysical monitoring network closely tracks subtle changes at the summit. These conditions present one of the best opportunities worldwide for understanding persistent lava lake behavior and the geophysical signals associated with open-vent basaltic eruptions. In this report, we provide a descriptive and visual summary of lava lake activity during 2016, a year consisting of continuous lava lake activity. The lake surface was composed of large black crustal plates separated by narrow incandescent spreading zones. The dominant motion of the surface was normally from north to south, but spattering produced transient disruptions to this steady motion. Spattering in the lake was common, consisting of one or more sites on the lake margin. The Overlook crater was continuously modified by the deposition of spatter (often as a thin veneer) on the crater walls, with frequent collapses of this adhered lava into the lake. Larger collapses, involving lithic material from the crater walls, triggered several small explosive events that deposited bombs and lapilli around the Halema‘uma‘u Crater rim, but these did not threaten public areas. The lava lake level varied over several tens of meters, controlled primarily by changes in summit magma reservoir pressure (in part driven by magma supply rates) and secondarily by fluctuations in spattering and gas release from the lake (commonly involving gas pistoning). The lake emitted a persistent gas plume, normally averaging 1,000–8,000 metric tons per day (t/d) of sulfur dioxide (SO 2 ), as well as a constant fallout of small juvenile and lithic particles, including Pele’s hair and tears. The gas emissions created volcanic air pollution (vog) that affected large areas of the Island of Hawai‘i. The summit eruption has been a major attraction for visitors in Hawai‘i Volcanoes National Park. During 2016, the rising lake levels allowed the lake and its spattering to be more consistently visible from public viewing areas, enhancing the visitor experience. The U.S. Geological Survey’s Hawaiian Volcano Observatory (HVO) closely monitors the summit eruption and keeps emergency managers and the public informed of activity.

Hawai'i

Flow characteristics and salinity patterns in tidal rivers within the northern Ten Thousand Islands, southwest Florida, water years 2007–19

Freshwater flow to the Ten Thousand Islands (TTI) estuary has been altered by the construction of the Tamiami Trail and construction of features in the now defunct Southern Golden Gate Estates development. This development included four associated canals that combine into the Faka Union Canal, which discharges into the TTI estuary. The Picayune Strand Restoration Project (PSRP) was initiated in 2007 to improve freshwater delivery to the TTI estuary by removing hundreds of miles of roads, emplacing hundreds of canal plugs, removing exotic vegetation, and constructing three pump stations. Quantifying the tributary flows and salinity patterns prior to, during, and after the restoration is essential to assessing the effectiveness of upstream restoration efforts. The U.S. Geological Survey, in cooperation with U.S. Army Corps of Engineers, initiated an ongoing study in 2006 to assess flow and salinity patterns in the TTI estuary. This is the second report by the U.S. Geological Survey describing flow characteristics and salinity patterns in the TTI area as part of the PSRP. This report describes flow characteristics and salinity patterns for the monitoring stations at Faka Union River, Pumpkin River, and East River and includes an assessment of salinity data from the Faka Union Boundary and Blackwater River water-quality stations for water years 2007–19. A water year is defined as the 12-month period from October 1 for any given year to September 30 of the following year. Annual and monthly variations in flow and salinity are often related to variations in rainfall with high and low annual flows (and below average and above average salinities) typically occurring during years with high and low annual rainfall, respectively. Monthly flows typically begin increasing in June and peak in September. Over the study period, positive trends in rainfall-adjusted monthly flow were detected at Faka Union River and East River, whereas no significant trend in flow was detected at Pumpkin River. Faka Union River is the largest contributor of freshwater to the TTI estuary, providing over 80 percent of the annual freshwater inflow to the estuary. The Faka Union Canal is expected to be the largest contributor of freshwater because until the PSRP is completed, the Faka Union Canal receives substantial drainage from multiple canals, which is not the case for Pumpkin and East Rivers. East River was the second largest contributor, followed by Pumpkin River. East River is downstream of the Fakahatchee Stand, which is a larger contributing area than the current contributing area for Pumpkin River. Monthly mean salinities were lowest at Faka Union River and East River, indicating that they received a greater amount of freshwater than the stations to the west. Negative trends in rainfall-adjusted salinity monthly means were observed at all monitoring stations during the study period. Increased trends in flow and decreased trends in salinity are attributed to increases in flow from upstream canals.

Florida

Assessment of streamflow trends in the eastern Dakotas, water years 1960–2019

Hydrologic extremes, whether periods of drought or flooding, are occurring more frequently with greater severity and can have substantial economic impacts. Along with flooding, the timing and volume of streamflow also is changing across the United States. The focus of this report is to characterize a unique trend in mean annual streamflow occurring in eastern North and South Dakota, hereafter referred to as the eastern Dakotas, that is not being observed anywhere else in the conterminous United States. Streamflow records for 1,853 U.S. Geological Survey streamgages obtained from the U.S. Geological Survey National Water Information System database with a continuous record of mean annual streamflow during water years 1960–2019 were included in this study. Using a Kendall tau statistical test ( p -value less than or equal to 0.10), 573 streamgages had a statistically significant upward trend in mean annual streamflow and are primarily located in the Midwest and northeastern United States. Of the streamgages, 182 had a statistically significant downward trend and are located primarily in the western and southeastern States. Several sites had increases in streamflow between 100 and 500 percent. Most of the streamgages with the highest increases in mean annual streamflow are along the same rivers in the eastern Dakotas, regardless of basin size. A comparison of mean annual streamflow of the last decade (2010–19) to the first decade (1960–69) of the study period shows that the largest increases in annual streamflow volumes in the United States also are in the eastern Dakotas. Among all 1,853 streamgages in the United States, the Sheyenne River near Warwick, North Dakota (U.S. Geological Survey station 05056000), has the greatest percent change, with an increase of 486 percent. Several factors may be contributing to increasing trends in streamflow in the eastern Dakotas and may include, in part, precipitation changes owing to climatic variation within the region, geologic makeup of the subsurface, and land-use changes. A better understanding of these research areas will help producers, resource managers, and infrastructure engineers to make more informed environmental and economic decisions.

North Dakota, South Dakota

Assessment of long-term changes in surface-water extent within Klamath Marsh, south-central Oregon, 1985–2021

The annual maximum extent of surface water in Klamath Marsh has naturally fluctuated in response to periods of wet and dry conditions in the surrounding basin. Field observations during the 2010s indicate that the annual maximum extent of surface water has been declining and the marsh is not responding to hydrologic inputs as it had historically. This report describes the results of a hydrologic evaluation of Klamath Marsh to characterize and understand multi-year declines in the surface-water extent and increased intermittency of streamflow exiting the marsh. Landsat imagery collected during 1985–2021 was processed to create a time series of annual maximum surface-water extent to assess changes in surface-water inundation within the marsh. A 50-percent decrease in the mean surface area of annual total open-water extent (TOWE) during the latter half of the study period (2003–21) compared to the first half (1985–2003) was observed in this 37-year time-series dataset. The change in open-water extent was offset by a corresponding increase in dry land in the marsh. Time series of streamflow, groundwater level, total annual precipitation, annual mean temperature, and anthropogenic water use and water management were compiled and evaluated to improve understanding of the factors affecting TOWE. Statistically significant downward trends in the regional groundwater table and streamflow into and out of the marsh were identified as well as statistically significant upward trends in annual mean temperature. Statistically significant correlations among TOWE, streamflow, and groundwater level also were identified. The decreasing trends could not be attributed to changes in total annual precipitation or changing anthropogenic groundwater use within the study area. Declines in the open-water extent of Klamath Marsh since 2000 principally are due to a decoupling of the groundwater and surface-water system beneath the marsh because of regional declines in groundwater level. Regional increases in air temperature and the reestablishment of more than 55,000 acres of forested land within the study area have likely contributed to increasing evapotranspiration, leaving less water available for groundwater recharge and stream base flow and resulting in basin-wide declines in streamflow and groundwater levels.

Oregon

Estimation of shallow ground-water recharge in the Great Lakes basin

This report presents the results of the first known integrated study of long-term average ground-water recharge to shallow aquifers (generally less than 100 feet deep) in the United States and Canada for the Great Lakes, upper St. Lawrence, and Ottawa River Basins. The approach used was consistent throughout the study area and allows direct comparison of recharge rates in disparate parts of the study area. Estimates of recharge are based on base-flow estimates for streams throughout the Great Lakes Basin and the assumption that base flow in a given stream is equal to the amount of shallow ground-water recharge to the surrounding watershed, minus losses to evapotranspiration. Base-flow estimates were developed throughout the study area using a single model based on an empirical relation between measured base-flow characteristics at streamflow-gaging stations and the surficial-geologic materials, which consist of bedrock, coarse-textured deposits, fine-textured deposits, till, and organic matter, in the surrounding surface-water watershed. Model calibration was performed using base-flow index (BFI) estimates for 959 stations in the U.S. and Canada using a combined 28,784 years of daily streamflow record determined using the hydrograph-separation software program PART. Results are presented for watersheds represented by 8-digit hydrologic unit code (HUC, U.S.) and tertiary (Canada) watersheds. Recharge values were lowest (1.6-4.0 inches/year) in the eastern Lower Peninsula of Michigan; southwest of Green Bay, Wisconsin; in northwestern Ohio; and immediately south of the St. Lawrence River northeast of Lake Ontario. Recharge values were highest (12-16.8 inches/year) in snow shadow areas east and southeast of each Great Lake. Further studies of deep aquifer recharge and the temporal variability of recharge would be needed to gain a more complete understanding of ground-water recharge in the Great Lakes Basin.

Great Lakes Basin

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