USGS Science⌕ Search

SEARCH · USGS Science

Results for “Investigator Annual 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,405 records · Page 78Linked to original sources

Potential effects of climate change on streamflow, eastern and western slopes of the Sierra Nevada, California and Nevada

Precipitation-runoff models of the East Fork Carson and North Fork American Rivers were developed and calibrated for use in evaluating the sensitivity of streamflow in the north-central Sierra Nevada to climate change. The East Fork Carson River drains part of the rain-shadowed, eastern slope of the Sierra Nevada and is generally higher than the North Fork American River, which drains the wetter, western slope. First, a geographic information system was developed to describe the spatial variability of basin characteristics and to help estimate model parameters. The result was a partitioning of each basin into noncontiguous, but hydrologically uniform, land units. Hydrologic descriptions of these units were developed and the Precipitation- Runoff Modeling System (PRMS) was used to simulate water and energy balances for each unit in response to daily weather conditions. The models were calibrated and verified using historical streamflows over 22-year (Carson River) and 42-year (American River) periods. Simulated annual streamflow errors average plus 10 percent of the observed flow for the East Fork Carson River basin and plus 15 percent for the North Fork American River basin. Interannual variability is well simulated overall, but, at daily scales, wet periods are simulated more accurately than drier periods. The simulated water budgets for the two basins are significantly different in seasonality of streamflow, sublimation, evapotranspiration, and snowmelt. The simulations indicate that differences in snowpack and snowmelt timing can play pervasive roles in determining the sensitivity of water resources to climate change, in terms of both resource availability and amount. The calibrated models were driven by more than 25 hypothetical climate-change scenarios, each 100 years long. The scenarios were synthesized and spatially disaggregated by methods designed to preserve realistic daily, monthly, annual, and spatial statistics. Simulated streamflow timing was not very sensitive to changes in mean precipitation, but was sensitive to changes in mean temperatures. Changes in annual streamflow amounts were amplified reflections of imposed mean precipitation changes, with especially large responses to wetter climates. In contrast, streamflow amount was surprisingly insensitive to mean temperature changes as a result of temporal links between peak snowmelt and the beginning of warm-season evapotranspiration. Comparisons of simulations driven by temporally detailed climate-model changes in which mean temperature changes vary from month to month and simulations in which uniform climate changes were imposed throughout the year indicate that the snowpack accumulates the influences of short-term conditions so that season average climate changes were more important than shorter term changes.

Water-Resources Investigations Report↗

Regional analysis of annual precipitation maxima in Montana

Dimensionless precipitation-frequency curves for estimating precipitation depths having large recurrence intervals were developed for 2-, 6-, and 24-hour storm durations for three homogeneous regions in Montana. Within each homogeneous region, at-site annual precipitation maxima were made dimensionless by dividing by the at-site mean and grouped so that a single frequency curve would be applicable for each duration. L-moment statistics were used to help define the homogeneous regions and to develop the dimensionless precipitation- frequency curves. Data from 459 precipitation stations were used after application of statistical tests to ensure that the data were not serially correlated and were stationary over the general period of data collection (1900-92). The data were found to have a small, but significant, degree of interstation correlation. The GEV distribution was used to construct dimensionless frequency curves of annual precipitation maxima for each duration within each region. Each dimensionless frequency curve was considered to be reliable for recurrence intervals up to the effective record length. Because of significant, though small, interstation correlation in all regions for all durations, and because the selected regions exhibited some heterogeneity, the effective record length was considered to be less than the total number of station-years of data. The effective record length for each duration in each region was estimated using a graphical method and found to range from 500 years for 6-hour duration data in Region 2 to 5,100 years for 24-hour duration data in Region 3.

Water-Resources Investigations Report↗

Depth-Duration Frequency of Precipitation for Oklahoma

A regional frequency analysis was conducted to estimate the depth-duration frequency of precipitation for 12 durations in Oklahoma (15, 30, and 60 minutes; 1, 2, 3, 6, 12, and 24 hours; and 1, 3, and 7 days). Seven selected frequencies, expressed as recurrence intervals, were investigated (2, 5, 10, 25, 50, 100, and 500 years). L-moment statistics were used to summarize depth-duration data and to determine the appropriate statistical distributions. Three different rain-gage networks provided the data (15minute, 1-hour, and 1-day). The 60-minute, and 1-hour; and the 24-hour, and 1-day durations were analyzed separately. Data were used from rain-gage stations with at least 10-years of record and within Oklahoma or about 50 kilometers into bordering states. Precipitation annual maxima (depths) were determined from the data for 110 15-minute, 141 hourly, and 413 daily stations. The L-moment statistics for depths for all durations were calculated for each station using unbiased L-mo-ment estimators for the mean, L-scale, L-coefficient of variation, L-skew, and L-kur-tosis. The relation between L-skew and L-kurtosis (L-moment ratio diagram) and goodness-of-fit measures were used to select the frequency distributions. The three-parameter generalized logistic distribution was selected to model the frequencies of 15-, 30-, and 60-minute annual maxima; and the three-parameter generalized extreme-value distribution was selected to model the frequencies of 1-hour to 7-day annual maxima. The mean for each station and duration was corrected for the bias associated with fixed interval recording of precipitation amounts. The L-scale and spatially averaged L-skew statistics were used to compute the location, scale, and shape parameters of the selected distribution for each station and duration. The three parameters were used to calculate the depth-duration-frequency relations for each station. The precipitation depths for selected frequencies were contoured from weighted depth surfaces to produce maps from which the precipitation depth-duration-frequency curve for selected storm durations can be determined for any site in Oklahoma.

Water-Resources Investigations Report↗

Factors affecting post-control reinvasion by seed of an invasive species, Phragmites australis , in the central Platte River, Nebraska.

Invasive plants, such as Phragmites australis , can profoundly affect channel environments of large rivers by stabilizing sediments and altering water flows. Invasive plant removal is considered necessary where restoration of dynamic channels is needed to provide critical habitat for species of conservation concern. However, these programs are widely reported to be inefficient. Post-control reinvasion is frequent, suggesting increased attention is needed to prevent seed regeneration. To develop more effective responses to this invader in the Central Platte River (Nebraska, USA), we investigated several aspects of Phragmites seed ecology potentially linked to post-control reinvasion, in comparison to other common species: extent of viable seed production, importance of water transport, and regeneration responses to hydrology. We observed that although Phragmites seed does not mature until very late in the ice-free season, populations produce significant amounts of viable seed (>50 % of filled seed). Most seed transported via water in the Platte River are invasive perennial species, although Phragmites abundances are much lower than species such as Lythrum salicaria , Cyperus esculentus and Phalaris arundinacea . Seed regeneration of Phragmites varies greatly depending on hydrology, especially timing of water level changes. Flood events coinciding with the beginning of seedling emergence reduced establishment by as much as 59 % compared to flood events that occurred a few weeks later. Results of these investigations suggest that prevention of seed set (i.e., by removal of flowering culms) should be a priority in vegetation stands not being treated annually. After seeds are in the seedbank, preventing reinvasion using prescribed flooding has a low chance of success given that Phragmites can regenerate in a wide variety of hydrologic microsites.

Nebraska↗

Hydrogeology and model-simulated groundwater availability in the Salt Fork Red River aquifer, southwestern Oklahoma, 1980–2015

The 1973 Oklahoma Water Law (82 OK Stat § 82-1020.5) requires that the Oklahoma Water Resources Board (OWRB) conduct hydrologic investigations of the State’s groundwater basins to support a determination of the maximum annual yield for each groundwater basin (hereinafter referred to as an “aquifer”). The maximum annual yield allocated per acre of land is known as the equal-proportionate-share (EPS) pumping rate. At present (2021), the OWRB has not yet established a maximum annual yield and EPS pumping rate for the Salt Fork Red River aquifer. To provide updated information to the OWRB that could support evaluation and determination of an appropriate maximum annual yield, the U.S. Geological Survey (USGS), in cooperation with the OWRB, conducted a hydrologic investigation and evaluated the effects of potential groundwater withdrawals on groundwater availability in the Salt Fork Red River aquifer. The Salt Fork Red River aquifer in Greer, Harmon, and Jackson Counties of southwestern Oklahoma is composed of about 274.5 square miles of alluvium and terrace deposits associated with the Salt Fork Red River. The mean annual recharge rate to the Salt Fork Red River aquifer for the period 1980–2015 was estimated to be about 2.94 inches per year, or 10.0 percent of the mean annual precipitation for the same period (29.4 inches per year). This 1980–2015 mean annual recharge rate is equivalent to a mean annual recharge rate of about 38,000 acre-feet per year (acre-ft/yr) for the Salt Fork Red River aquifer excluding about 19,764 acres comprising the Mulberry Creek and Horse Creek terraces. The mean annual recharge rates upgradient and downgradient from USGS streamgage 07300500 Salt Fork Red River at Mangum, Okla. (hereinafter referred to as the “Mangum gage”), apportioned by aquifer area (41.5 and 58.5 percent, respectively), were about 16,000 and 22,000 acre-ft/yr, respectively. Mean annual groundwater use for the study period (1980–2015) was 3,532.7 acre-ft/yr; about 77 percent of that groundwater use was for irrigation, and about 23 percent was for public supply. Most groundwater use for irrigation was associated with wells in the Martha terrace. A hydrogeologic framework was developed for the Salt Fork Red River aquifer and included a definition of the aquifer extent and potentiometric surface, as well as a description of the textural and hydraulic properties of aquifer materials. The hydrogeologic framework was used in the construction of the numerical groundwater-flow model of the Salt Fork Red River aquifer described in this report. A conceptual model for the Salt Fork Red River aquifer that reasonably represents the groundwater-flow system was developed to constrain the construction and calibration of the numerical model. The conceptual-model water budget estimated mean annual inflows to, and outflows from, the Salt Fork Red River aquifer for the period 1980–2015 and included a subaccounting of mean annual inflows and outflows for the portions of the aquifer that were upgradient and downgradient from the Mangum gage. The numerical groundwater-flow model of the Salt Fork Red River aquifer was constructed by using MODFLOW-2005 with the Newton formulation solver. The model of the Salt Fork Red River aquifer was spatially discretized into 1,050 rows, 1,125 columns, about 170,000 active cells measuring 200 by 200 feet (ft), and a single convertible layer. The model was temporally discretized into 432 monthly transient stress periods (each with two time steps to improve model stability). An initial steady-state stress period represented mean annual inflows to, and outflows from, the aquifer and produced a solution that was used as the initial condition for subsequent transient stress periods as well as some groundwater-availability scenarios. The model was calibrated to water-table-altitude observations at selected wells and base-flow observations at selected streamgages. The simulated saturated thickness of the Salt Fork Red River aquifer was determined by subtracting the altitude of the aquifer base from the simulated water-table altitude at the end of the numerical-model period (2015). The simulated saturated thickness was more than 75 ft in a paleochannel in the Dodson terrace near the Texas border. The mean aquifer thickness (sum of saturated and unsaturated) was 49.62 ft, and the mean saturated thickness was 28.55 ft. A simulated mean transmissivity of 1,024 feet squared per day was computed from the calibrated hydraulic conductivity and saturated thickness of each cell. The simulated available water in storage at the end of the numerical-model period (2015) was 526,117 acre-feet (acre-ft); about 42 percent of that total was available upgradient from the Mangum gage, and about 58 percent of that total was available downgradient from the Mangum gage (including the Mangum terrace). Three types of groundwater-availability scenarios were run using the calibrated numerical model. These scenarios were used to (1) estimate the EPS pumping rate that ensures a minimum 20-, 40-, and 50-year life of the aquifer, (2) quantify the potential effects of projected well withdrawals on groundwater storage over a 50-year period, and (3) simulate the potential effects of a hypothetical 10-year drought on base flow and groundwater storage. The 20-, 40-, and 50-year EPS pumping rates under normal recharge conditions were about 0.51, 0.48, and 0.48 acre-foot per acre per year, respectively. Given the 155,929-acre modeled aquifer area, these rates correspond to annual yields of about 78,800, 74,900, and 74,700 acre-ft/yr, respectively. For the 20-year EPS scenario, decreasing and increasing recharge by 10 percent resulted in a 6-percent change in the EPS pumping rate in both cases; for the 40- and 50-year EPS scenarios, decreasing and increasing recharge by 10 percent resulted in a 7-percent change in the EPS pumping rate in both cases. Projected 50-year pumping scenarios were used to simulate the effects of selected well withdrawal rates on groundwater storage of the Salt Fork Red River aquifer and base flows in the Salt Fork Red River. The effects of well withdrawals were evaluated by quantifying differences in groundwater storage and base flow in four 50-year scenarios, which applied (1) no groundwater pumping, (2) mean pumping rates for the study period (1980–2015), (3) 2015 pumping rates, and (4) increasing demand pumping rates at simulated wells. The increasing demand pumping rates assumed a cumulative 20.4-percent increase in pumping over 50 years based on 2010–60 demand projections for southwestern Oklahoma. Groundwater storage after 50 years with no pumping was 535,000 acre-ft, or 8,900 acre-ft (1.7 percent) greater than the initial groundwater storage; this groundwater storage increase is equivalent to a mean water-table-altitude increase of 0.48 ft. Groundwater storage after 50 years of pumping at the mean rate for the study period (1980–2015) was 519,900 acre-ft, or 6,200 acre-ft (1.2 percent) less than the initial groundwater storage; this groundwater storage decrease is equivalent to a mean water-table-altitude decline of 0.34 ft. Groundwater storage at the end of the 50-year period with 2015 pumping rates was 513,100 acre-ft, or 13,000 acre-ft (2.5 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean water-table-altitude decline of 0.71 ft. Groundwater storage at the end of the 50-year period with increasing demand pumping rates was 509,700 acre-ft, or 16,500 acre-ft (3.1 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean water-table-altitude decline of 0.89 ft. A hypothetical 10-year drought scenario was used to simulate the effects of a prolonged period of reduced recharge on groundwater storage. The period January 1983–December 1992 was chosen as the simulated drought period. Drought effects were quantified by comparing the results of the drought scenario to those of the calibrated numerical model (no drought) at the end of the simulated drought period (1992). To simulate the hypothetical drought, recharge in the calibrated numerical model was reduced by 50 percent during the simulated drought period (1983–92). Upstream inflows from the Salt Fork Red River, Turkey Creek, and Bitter Creek were reduced by 75 percent. Groundwater storage at the end of the drought period (1992) was 479,200 acre-ft, or 53,200 acre-ft (10.0 percent) less than the groundwater storage of the calibrated numerical model at the end of the drought period. This decrease in groundwater storage is equivalent to a mean water-table-altitude decline of 2.9 ft. At the end of the 10-year hypothetical drought period, simulated base flows at the Mangum gage and USGS streamgage 07301110 Salt Fork Red River near Elmer, Okla., had decreased by about 80 and 70 percent, respectively.

Oklahoma↗

Modified Methodology for Projecting Coastal Louisiana Land Changes over the Next 50 Years

The coastal Louisiana landscape is continually undergoing geomorphologic changes (in particular, land loss); however, after the 2005 hurricane season, the changes were intensified because of Hurricanes Katrina and Rita. The amount of land loss caused by the 2005 hurricane season was 42 percent (562 km2) of the total land loss (1,329 km2) that was projected for the next 50 years in the Louisiana Coastal Area (LCA), Louisiana Ecosystem Restoration Study. The purpose of this study is to provide information on potential changes to coastal Louisiana by using a revised LCA study methodology. In the revised methodology, we used classified Landsat TM satellite imagery from 1990, 2001, 2004, and 2006 to calculate the 'background' or ambient land-water change rates but divided the Louisiana coastal area differently on the basis of (1) geographic regions ('subprovinces') and (2) specific homogeneous habitat types. Defining polygons by subprovinces (1, Pontchartrain Basin; 2, Barataria Basin; 3, Vermilion/Terrebonne Basins; and 4, the Chenier Plain area) allows for a specific erosion rate to be applied to that area. Further subdividing the provinces by habitat type allows for specific erosion rates for a particular vegetation type to be applied. Our modified methodology resulted in 24 polygons rather than the 183 that were used in the LCA study; further, actively managed areas and the CWPPRA areas were not masked out and dealt with separately as in the LCA study. This revised methodology assumes that erosion rates for habitat types by subprovince are under the influence of similar environmental conditions (sediment depletion, subsidence, and saltwater intrusion). Background change rate for three time periods (1990-2001, 1990-2004, and 1990-2006) were calculated by taking the difference in water or land among each time period and dividing it by the time interval. This calculation gives an annual change rate for each polygon per time period. Change rates for each time period were then used to compute the projected change in each subprovince and habitat type over 50 years by using the same compound rate functions used in the LCA study. The resulting maps show projected land changes based on the revised methodology and inclusion of damage by Hurricanes Katrina and Rita. Comparison of projected land change values between the LCA study and this study shows that this revised methodology - that is, using a reduced polygon subset (reduced from 183 to 24) based on habitat type and subprovince - can be used as a quick projection of land loss.

Scientific Investigations Map↗

Hydrology of and Current Monitoring Issues for the Chicago Area Waterway System, Northeastern Illinois

The Chicago Area Waterway System (CAWS) consists of a combination of natural and manmade channels that form an interconnected navigable waterway of approximately 90-plus miles in the metropolitan Chicago area of northeastern Illinois. The CAWS serves the area as the primary drainage feature, a waterway transportation corridor, and recreational waterbody. The CAWS was constructed by the Metropolitan Water Reclamation District of Greater Chicago (MWRDGC). Completion of the Chicago Sanitary and Ship Canal (initial portion of the CAWS) in 1900 breached a low drainage divide and resulted in a diversion of water from the Lake Michigan Basin. A U.S. Supreme Court decree (Consent Decree 388 U.S. 426 [1967] Modified 449 U.S. 48 [1980]) limits the annual diversion from Lake Michigan. While the State of Illinois is responsible for the diversion, the MWRDGC regulates and maintains water level and water quality within the CAWS by using several waterway control structures. The operation and control of water levels in the CAWS results in a very complex hydraulic setting characterized by highly unsteady flows. The complexity leads to unique gaging requirements and monitoring issues. This report provides a general discussion of the complex hydraulic setting within the CAWS and quantifies this information with examples of data collected at a range of flow conditions from U.S. Geological Survey streamflow gaging stations and other locations within the CAWS. Monitoring to address longstanding issues of waterway operation, as well as current (2014) emerging issues such as wastewater disinfection and the threat from aquatic invasive species, is included in the discussion.

Illinois↗

Geochemistry and hydrology of thermal springs in the Idaho Batholith and adjacent areas, central Idaho

The occurrence of nature of thermal springs in the Idaho batholith and adjacent areas suggest a relation between structural controls and deeply circulating hot-water systems. Springs issuing from granitic rocks are associated mostly with major regional fault structures. Springs issuing from other rocks probably are related to local faulting. Individual spring flows and water temperatures are variable and range from less than 1 gallon per minute to 2,710 gallons per minute and from 20.5 degrees to 94.0 degrees Celsius. Annual spring discharge is at least 27,000 acre-feet; heat discharges convectively is estimated to be 5.0 x 10 7 calories per second. Thermal springs discharge relatively dilute water; dissolved solids range from 103 to 839 milligrams per liter. The chemical quality of the water suggests deep circulation of meteoric water. Estimated reservoir temperatures are generally less than 100 degrees Celsius, but temperatures for several springs exceed 150 degrees Celsius. Stable-isotope data suggest that most of the thermal water is not derived from current precipitation. Carbon-14 values indicate that thermal waters are old; apparent residence times range from 9,000 to more than 40,000 years.

Idaho↗

Methods for estimating the magnitude and frequency of peak discharges of rural, unregulated streams in Virginia

Methods are presented for estimating the peak discharges of rural, unregulated streams in Virginia. A Pearson Type III distribution is fitted to the logarithms of the unregulated annual peak-discharge records from 363 stream-gaging stations in Virginia to estimate the peak discharge at these stations for recurrence intervals of 2 to 500 years. Peak-discharge characteristics for 284 unregulated stations are divided into eight regions based on physiographic province, and regressed on basin characteristics, including drainage area, main channel length, main channel slope, mean basin elevation, percentage of forest cover, mean annual precipitation, and maximum rainfall intensity. Regression equations for each region are computed by use of the generalized least-squares method, which accounts for spatial and temporal correlation between nearby gaging stations. This regression technique weights the significance of each station to the regional equation based on the length of records collected at each cation, the correlation between annual peak discharges among the stations, and the standard deviation of the annual peak discharge for each station. Drainage area proved to be the only significant explanatory variable in four regions, while other regions have as many as three significant variables. Standard errors of the regression equations range from 30 to 80 percent. Alternate equations using drainage area only are provided for the five regions with more than one significant explanatory variable. Methods and sample computations are provided to estimate peak discharges at gaged and engaged sites in Virginia for recurrence intervals of 2, 5, 10, 25, 50, 100, 200, and 500 years, and to adjust the regression estimates for sites on gaged streams where nearby gaging-station records are available.

Virginia↗

Monitoring breeding and survival of ring-necked pheasant (Phasianus colchicus) in the Sacramento Valley, Sacramento-San Joaquin River Delta, and Klamath Basin, northern California—Five-year summary, 2013–17

The U.S. Geological Survey Western Ecological Research Center, Pheasants Forever, Mandeville Island Duck Club, and the California Department of Fish and Wildlife collaborated in a reconnaissance study to monitor populations of ring-necked pheasant ( Phasianus colchicus ) using radio-telemetry in the Sacramento Valley, Sacramento-San Joaquin River Delta, and Klamath Basin of northern California. The purpose of this study was to provide agencies and private landowners with a framework of decision-support tools to help manage pheasant populations in California. During winter, spring, and autumn of 2013–17, we radio- or Global Positioning System-marked 227 female pheasant across six study sites. Data collection was focused on investigating nest-site and brood-rearing habitat selection, examining avian predator composition, and estimating population vital rates to improve our understanding of pheasant population dynamics and to identify factors that may contribute to decreases in pheasant populations in California. The cumulative annual adult survival probability across all sites during 2013–17 was 27.6 percent (95-percent confidence interval [CI], 21.9–33.6), and the cumulative nest and brood survival probabilities were 34.5 percent (95-percent CI, 27.0–42.2) and 54.2 percent (95-percent CI, 43.7–63.5), respectively. Evidence from microhabitat surveys completed at nest-sites, brood locations, and random locations suggested that marked female pheasant tended to select increasing vertical cover and residual vegetation cover and tended to avoid areas of increasing bare ground cover regardless of life-history stage. However, females at nest-sites selected increasing grass cover and height, whereas brood-rearing females tended to select increasing forb cover and height. Only perennial grass cover and perennial grass height were shown to have a positive influence on nest survival, which suggests that increasing perennial grass cover in areas occupied by pheasant may increase nest survival. Analysis of environmental factors linked to vital rate information are ongoing and will continue with investigations at increased spatial scales (that is, macro-habitat) to develop integrated population models that can incorporate abundance estimates from crow count data with vital rates from telemetry data. This report includes results from 5 years of data collection and should be interpreted with caution, as these findings are preliminary.

California↗

Summit crater lake observations, and the location, chemistry, and pH of water samples near Mount Chiginagak volcano, Alaska: 2004-2012

Mount Chiginagak is a hydrothermally active volcano on the Alaska Peninsula, approximately 170 km south–southwest of King Salmon, Alaska (fig. 1). This small stratovolcano, approximately 8 km in diameter, has erupted through Tertiary to Permian sedimentary and igneous rocks (Detterman and others, 1987). The highest peak is at an elevation of 2,135 m, and the upper ~1,000 m of the volcano are covered with snow and ice. Holocene activity consists of debris avalanches, lahars, and lava flows. Pleistocene pyroclastic flows and block-and-ash flows, interlayered with andesitic lava flows, dominate the edifice rocks on the northern and western flanks. Historical reports of activity are limited and generally describe “steaming” and “smoking” (Coats, 1950; Powers, 1958). Proximal tephra collected during recent fieldwork suggests there may have been limited Holocene explosive activity that resulted in localized ash fall. A cluster of fumaroles on the north flank, at an elevation of ~1,750 m, commonly referred to as the “north flank fumarole” have been emitting gas throughout historical time (location shown in fig. 2). The only other thermal feature at the volcano is the Mother Goose hot springs located at the base of the edifice on the northwestern flank in upper Volcano Creek, at an elevation of ~160 m (fig. 2, near sites H1, H3, and H4). Sometime between November 2004 and May 2005, a ~400-m-wide, 100-m-deep lake developed in the snow- and ice-filled summit crater of the volcano (Schaefer and others, 2008). In early May 2005, an estimated 3 million cubic meters (3×106 m3) of sulfurous, clay-rich debris and acidic water exited the crater through tunnels at the base of a glacier that breaches the south crater rim. More than 27 km downstream, these acidic flood waters reached approximately 1.3 m above normal water levels and inundated a fertile, salmon-spawning drainage, acidifying the entire water column of Mother Goose Lake from its surface waters to its maximum depth of 45 m (resulting pH ~2.9), and preventing the annual salmon run in the King Salmon River. A simultaneous release of gas and acidic aerosols from the crater caused widespread vegetation damage along the flow path. Since 2005, we have been monitoring the crater lake water that continues to flow into Mother Goose Lake by collecting surface water samples for major cation and anion analysis, measuring surface-water pH of affected drainages, and photo-documenting the condition of the summit crater lake. This report describes water sampling locations, provides a table of chemistry and pH measurements, and documents the condition of the summit crater between 2004 and 2011. In September 2013, the report was updated with results of water-chemistry samples collected in 2011 and 2012, which were added as an addendum.

Alaska↗

Geochemistry of the Mattole River in Northern California

The chemical composition of streams can vary greatly with changing discharge during storm runoff. These chemical changes are related to the pathways of various water parcels from the time they fall as rain until they enter the stream, and to the interactions between water and sediment during transport downstream. In order to understand better the chemical variations during storms, an extensive investigation was made of the Mattole River, a chemically clean coastal stream in Mendocino County, California. The Mattole drains a topographically mature basin of 620 sw km which has relief of about 1200 m, a long summer dry season, and mean annual rainfall of about 2300 mm. The stream flow is composed of seasonally varying proportions of four flow components, namely, surface runoff, quick-return flow (rainfall having brief and intimate contact with the soil before entering the surface drainage), delayed-return flow, and base runoff. Each component is identified by its characteristic chemistry and by the time delay between rainfall and entrance into the stream. Information is also presented on rain chemistry, adsorption reactions of suspended sediments in the fresh and brackish environments, and compositional variation of river sediments with particle size. (Woodard-USGS)

Open-File Report↗

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

Frequency analysis on annual peak streamflow (hereinafter, peak flow) is essential to water-resources management applications, including critical structure design (for example, bridges and culverts) and floodplain mapping. Nonstationarity is a statistical property of a peak-flow series such that the distributional properties (the mean, variance, or skew) change either gradually (monotonic trend) or abruptly (shift, step change or change point) through time. Not incorporating or accounting for observed nonstationarity into peak-flow frequency analysis might result in a poor representation of the true probability of large floods and thus misrepresent the actual flood risks to life and property. This report summarizes how hydroclimatic variability might affect the temporal and spatial distributions of peak-flow data in the State of Montana (and northern Wyoming) and is part of a larger study to document peak-flow nonstationarity and hydroclimatic changes across a nine-State region consisting of Illinois, Iowa, Michigan, Minnesota, Missouri, Montana, North Dakota, South Dakota, and Wisconsin. A wide range of analyses and statistical approaches are applied to document the primary mechanisms controlling floods and characterize temporal changes in hydroclimatic variables and peak flows. This study was completed in cooperation with the Montana Department of Natural Resources and Conservation. The purpose of this report is to characterize temporal and spatial patterns of nonstationarity in peak flows and hydroclimatology in Montana and northern Wyoming. In this evaluation, peak-flow, daily streamflow, and model-simulated gridded climatic data were examined for monotonic trends, change points, and other statistical properties that might indicate changing climatic and environmental conditions. This report includes background information on the study area, the history of U.S. Geological Survey peak-flow data collection and frequency analysis in Montana, and the review of research relating to hydroclimatic variability and change in Montana. This study might help provide a framework for addressing potential nonstationarity issues in peak-flow frequency updates that commonly are completed by the U.S. Geological Survey in cooperation with other agencies throughout the Nation. The analytical structure of this study includes analyses of monotonic trends and change points in numerous hydroclimatic variables in assigned 30-, 50-, 75-, and 100-year analysis periods. For Montana and part of Wyoming, the 30-, 50-, 75, and 100-year analyses included 157, 70, 48, and 12 streamgages, respectively. For those streamgages, nonstationarities were analyzed in the following variables: (1) climatic variables, including annual and seasonal (winter, spring, summer, and fall) temperature and precipitation; (2) daily streamflow variables, including the annual center of volume duration, annual center of volume median, and peaks over threshold with a mean of four events per year; and (3) annual peak-flow variables, including peak-flow timing and magnitude. A likelihood approach was used to express statistical confidence and assign the nonstationarity results as likely upward or downward (highest statistical confidence), somewhat likely upward or downward (less statistical confidence), or about as likely as not (little statistical confidence; hereinafter, neutral). For the nonstationarity analyses of the climatic, daily streamflow, and peak-flow variables, the results are presented in detail and discussed with respect to statewide patterns and geographic variability. For each of the 30-, 50-, and 75-year analyses, peak-flow change-point and monotonic trend analyses were compiled for streamgages classified with likely downward or likely upward trends. For those streamgages, the associated basin characteristics and nonstationarity results for peak-flow timing, daily streamflow, and climatic variables were investigated and statistically compared to discern associations among other variables that might contribute to the peak-flow nonstationarity results. The 50- and 75-year peak-flow nonstationarities identified in this study are mostly downward, in association with mostly upward temperature and potential evapotranspiration:precipitation monotonic trends. For the 50-, 75-, and 100-year analyses, the peak-flow change points are predominantly downward and are concentrated in the 1970s and 1980s, which indicates general consistency among the longer trend periods. These findings are in association with substantial research documenting globally rising temperature and atmospheric greenhouse gas concentrations that might be largely attributed to anthropogenic activities. Anthropogenic effects might represent long-term (on the order of several decades to more than a century) climate changes that might happen within highly variable natural climate fluctuations. Several paleo studies in the north-central United States have indicated that hydroclimatic extremes (that is, low- and high-streamflow conditions) before European settlement have been outside of extremes since the 1900s. Depending on the interactions of anthropogenic effects and natural climate variability, extreme high-streamflow conditions might occur in the future, even in the presence of long-term downward peak-flow trends.

Montana, Wyoming↗

Rates of evapotranspiration, recharge from precipitation beneath selected areas of native vegetation, and streamflow gain and loss in Carson Valley, Douglas County, Nevada, and Alpine County, California

Rapid growth and development in Carson Valley is causing concern over the continued availability of water resources to sustain such growth into the future. A study to address concerns over water resources and to update estimates of water-budget components in Carson Valley was begun in 2003 by the U.S. Geological Survey, in cooperation with Douglas County, Nevada. This report summarizes micrometeorologic, soil-chloride, and streambed-temperature data collected in Carson Valley from April 2003 through November 2004. Using these data, estimates of rates of discharge by evapotranspiration (ET), rates of recharge from precipitation in areas of native vegetation on the eastern and northern sides of the valley, and rates of recharge and discharge from streamflow infiltration and seepage on the valley floor were calculated. These rates can be used to develop updated water budgets for Carson Valley and to evaluate potential effects of land- and water-use changes on the valley's water budget. Data from eight ET stations provided estimates of annual ET during water year 2004, the sixth consecutive year of a drought with average or below average precipitation since 1999. Estimated annual ET from flood-irrigated alfalfa where the water table was from 3 to 6 feet below land surface was 3.1 feet. A similar amount of ET, 3.0 feet, was estimated from flood-irrigated alfalfa where the water table was about 40 feet below land surface. Estimated annual ET from flood-irrigated pasture ranged from 2.8 to 3.2 feet where the water table ranged from 2 to 5 feet below land surface, and was 4.4 feet where the water table was within 2 feet from land surface. Annual ET estimated from nonirrigated pasture was 1.7 feet. Annual ET estimated from native vegetation was 1.9 feet from stands of rabbitbrush and greasewood near the northern end of the valley, and 1.5 feet from stands of native bitterbrush and sagebrush covering alluvial fans along the western side of the valley. Uncertainty in most ET estimates is about 12 percent, but ranged from +30 and +50 percent to -20 and -40 percent for nonirrigated pasture and native bitterbrush and sagebrush. Estimated rates for water year 2004 likely are less than those during years of average, or above average precipitation when the water table would be closer to land surface. Test holes drilled in areas of native vegetation on the northern and eastern sides of Carson Valley had high concentrations of soil chloride at depths ranging from 4 to 18 feet below land surface at six locations on the eastern side of the valley. The high chloride concentrations indicate that modern-day precipitation at the six locations does not percolate deeper than the root zone of native vegetation. Estimates of the time required to accumulate the measured amount of chloride to depths of about 30 feet below land surface at the six test holes ranged from about 3,000 to 12,000 years. Low concentrations of soil chloride in two test holes on the northern end of Carson Valley and in a test hole on the eastern side of Fish Spring Flat indicate that a small amount of recharge from modern-day precipitation is taking place. Estimated annual recharge from precipitation at the two locations was 0.03 and 0.04 foot on the northern end of the valley and 0.02 foot on the eastern side of Fish Spring Flat. Uncertainty in the estimated recharge rates was about ?0.01 foot. Estimates of the time required to accumulate the measured amount of chloride to depths of about 30 feet below land surface at the three test holes ranged from about 100 to 700 years. The two test holes near the northern end of the valley are in gravel and eolian sand deposits and recharge from precipitation may be taking place at similar rates in other areas with gravel and eolian sand deposits. Based on results from other test holes, recharge at the rate estimated for the test hole on the eastern side of Fish Spring Flat is not likely applicable to a large area. Data from 37 site

Scientific Investigations Report↗

Floods of September 16, 1975 in the Tallaboa Valley, Puerto Rico

This report provides a record of the flood of September 16, 1975, and associated hydrologic data. These data can be used in making rational decisions in formulating effective flood-plain regulations that would minimize flood problems in the Tallaboa Valley. The Tallaboa Valley (fig. 1) lies on the southwestern coast of Puerto Rico about 12 km (7.5 mi) west of the city of Ponce. The climate is semiarid with a mean annual precipitation of about 1,020 mm (40 in). It is a relatively small valley which until the late fifties was dedicated mainly to the product ion of sugarcane. Since then, the economy started changing from agriculture to industry. At present the lower part of the valley is completely industrialized. The upper part of the valley is used for agriculture--mostly as pastureland. The reach from Highway 132 downstream to the new Highway 2 is mainly agricultural with a few scattered communities. The lower part of the valley that includes the area from the Highway 2 bridge to the mouth is covered with heavy industry concerned with the refinement of oil, and associated industries. The continuous growth of industry and consequently of the urban area of Peñuelas has encouraged growth onto the flood plain. The study area is subject to flooding by the Rio Tallaboa and the Rio Guayanes, its main tributary. The Rio Guayanes is the- cause of flooding in the urban area of Penuelas, encompassing the reach from the northern edge of town downstream to Highway 132 where it joins the Rio Tallaboa. Data are generally referred to in SI (International System) units followed by corresponding inch-pound units in parentheses. The SI units.may be converted to inch-pound units by multiplying the units given by the factors shown.

Tallaboa Valley↗

Evaluation of thermograph data for California streams

Statistical analysis of water-temperature data from California streams indicates that, for most purposes, long-term operation of thermographs (automatic water-temperature recording instruments) does not provide a more useful record than either short-term operation of such instruments or periodic measurements. Harmonic analyses were made of thermograph records 5 to 14 years in length from 82 stations. More than 80 percent of the annual variation in water temperature is explained by the harmonic function for 77 of the 82 stations. Harmonic coefficients based on 8 years of thermograph record at 12 stations varied only slightly from coefficients computed using two equally split 4-year records. At five stations where both thermograph and periodic (10 to 23 measurements per year) data were collected concurrently, harmonic coefficients for periodic data were defined nearly as well as those for thermograph data. Results of this analysis indicate that, except where detailed surveillance of water temperatures is required or where there is a chance of temporal change, thermograph operations can be reduced substantially without affecting the usefulness of temperature records.

California↗

Cost-effectiveness of the U.S. Geological Survey stream-gaging program in Indiana

Analysis of the stream gaging program in Indiana was divided into three phases. The first phase involved collecting information concerning the data need and the funding source for each of the 173 surface water stations in Indiana. The second phase used alternate methods to produce streamflow records at selected sites. Statistical models were used to generate stream flow data for three gaging stations. In addition, flow routing models were used at two of the sites. Daily discharges produced from models did not meet the established accuracy criteria and, therefore, these methods should not replace stream gaging procedures at those gaging stations. The third phase of the study determined the uncertainty of the rating and the error at individual gaging stations, and optimized travel routes and frequency of visits to gaging stations. The annual budget, in 1983 dollars, for operating the stream gaging program in Indiana is $823,000. The average standard error of instantaneous discharge for all continuous record gaging stations is 25.3%. A budget of $800,000 could maintain this level of accuracy if stream gaging stations were visited according to phase III results. A minimum budget of $790,000 is required to operate the gaging network. At this budget, the average standard error of instantaneous discharge would be 27.7%. A maximum budget of $1 ,000,000 was simulated in the analysis and the average standard error of instantaneous discharge was reduced to 16.8%. (Author 's abstract)

Water-Resources Investigations Report↗

Hydrogeology of, and simulation of ground-water flow in, a mantled carbonate-rock system, Cumberland Valley, Pennsylvania

The U.S. Geological Survey conducted a study in a highly productive and complex regolith-mantled carbonate valley in the northeastern part of the Cumberland Valley, Pa., as part of its Appalachian Valleys and Piedmont Regional Aquifer-system Analysis program. The study was designed to quantify the hydrogeologic characteristics and understand the ground-water flow system of a highly productive and complex thickly mantled carbonate valley. The Cumberland Valley is characterized by complexly folded and faulted carbonate bedrock in the valley bottom, by shale and graywacke to the north, and by red-sedimentary and diabase rocks in the east-southeast. Near the southern valley hillslope, the carbonate rock is overlain by wedge-shaped deposit of regolith, up to 450 feet thick, that is composed of residual material, alluvium, and colluvium. Locally, saturated regolith is greater than 200 feet thick. Seepage-run data indicate that stream reaches, near valley walls, are losing water from the stream, through the regolith, to the ground-water system. Results of hydrograph-separation analyses indicate that base flow in stream basins dominated by regolith-mantled carbonate rock, carbonate rock, and carbonate rock and shale are 81.6, 93.0, and 67.7 percent of total streamflow, respectively. The relative high percentage for the regolith-mantled carbonate-rock basin indicates that the regolith stores precipitation and slowly, steadily releases this water to the carbonate-rock aquifer and to streams as base flow. Anomalies in water-table gradients and configuration are a result of topography and differences in the character and distribution of overburden material, permeability, rock type, and geologic structure. Most ground-water flow is local, and ground water discharges to nearby springs and streams. Regional flow is northeastward to the Susquehanna River. Average-annual water budgets were calculated for the period of record from two continuous streamflow-gaging stations. Average-annual precipitation range from 39.0 to 40.5 inches, and averages about 40 inches for the model area. Average-annual recharge, which was assumed equal to the average-annual base flow, ranged from 12 inches for the Conodoguinet Creek, and 15 inches for the Yellow Breeches Creek. The thickly-mantled carbonate system was modeled as a three- dimensional water-table aquifer. Recharge to, ground-water flow through, and discharge from the Cumberland Valley were simulated. The model was calibrated for steady-state conditions using average recharge and discharge data. Aquifer horizontal hydraulic conductivity was calculated from specific-capacity data for each geologic unit in the area. Particle-tracking analyses indicate that interbasin and intrabasin flows of groundwater occur within the Yellow Breeches Creek Basin and between the Yellow Breeches and Conodoguinet Creek Basins.

Pennsylvania↗