Search USGS⌕ Search

SEARCH · Search USGS

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 217 records · Page 12Linked to original sources

Seasonal patterns of weight, hematology, and serum characteristics of free-ranging female white-tailed deer in Minnesota

Weights, hematology, and serum profIles of white-tailed does in the central Superior National Forest of northeastern Minnesota were examined year-around to determine seasonal patterns of nutritional condition and metabolism. Deer were initially captured by Clover trap or rocket net. Between 15 February 1989 and 23 January 1990, we recaptured 12 adult (> 1.5 years) female deer 1-9 times each (a total of 59 recaptures) using a radio-controlled capture collar. Monthly weights of deer exhibited a cyclic seasonal pattern. Mean weight declined 22 % from February to an annual minimum during May, then steadily increased 45 % to a maximum in October. Seasonal patterns were most evident for hemoglobin concentration, red blood cells, packed cell volume, serum total protein, urea nitrogen, creatinine, the urea N to creatinine ratio, triiodothyronine, cortisol, and potassium. Wide seasonal variations of these characteristics were indicative of shifts in the deer's metabolic physiology. Although seasonal metabolic shifts are partially attributable to an endogenous rhythm, the intensity of, their expression was most likely affected by nutritional changes and concomitant alterations of body condition. Annual changes in seasonal trends of blood characteristics may be useful in investigating nutritional effects of specific environmental and demographic factors. We compare our findings with those reported for deer on ranges farther south.

Minnesota↗

Analysis of the Shallow Groundwater Flow System at Fire Island National Seashore, Suffolk County, New York

Fire Island National Seashore (FIIS) occupies 42 kilometers of the barrier island for which it is named that lies off the southern shore of Suffolk County, N.Y. Freshwater in the highly permeable, sandy aquifer underlying Fire Island is bounded laterally by marine surface waters and at depth by saline groundwater. Interspersed throughout FIIS are 17 pre-existing residential communities that in summer months greatly increase in population through the arrival of summer residents and vacationers; in addition, the National Park Service (NPS) has established several facilities on the island to accommodate visitors to FIIS. The 2.2 million people estimated by the NPS to visit Fire Island annually impact groundwater quality through the release of waste-derived contaminants, such as nutrients, pathogens, and organic compounds, into the environment. Waste-contaminated groundwater can move through the aquifer and threaten the ecological health of the adjacent back-barrier estuaries to which much of the groundwater ultimately discharges. In 2004, the U.S. Geological Survey (USGS), in cooperation with the NPS, began a 3-year investigation to (1) collect groundwater levels and water-quality (nutrient) samples, (2) develop a three-dimensional model of the shallow (water-table) aquifer system and adjacent marine surface waters, and (3) calculate nitrogen loads in simulated groundwater discharges from the aquifer to back-barrier estuaries and the ocean. The hydrogeology of the shallow aquifer system was characterized from the results of exploratory drilling, geophysical surveying, water-level monitoring, and water-quality sampling. The investigation focused on four areas-the communities of Kismet and Robbins Rest, the NPS Visitor Center at Watch Hill, and the undeveloped Otis Pike Fire Island High Dune Wilderness. Thirty-five observation wells were installed within FIIS to characterize subsurface hydrogeology and establish a water-table monitoring network in the four study areas. A variable-density model of the shallow aquifer system and adjacent marine surface waters was developed to simulate groundwater flow patterns and rates. Nitrogen loads from the shallow aquifer system were calculated from representative total nitrogen (TN) concentrations and simulated groundwater discharges to back-barrier estuaries and the ocean. The model simulates groundwater directions, velocities, and discharge rates under 2005 mean annual conditions. Groundwater budgets were developed for recharge areas of similar land use that contribute freshwater to back-barrier estuaries, the ocean, and subsea-discharge zones. Total freshwater discharge from the shallow aquifer system is about 43,500 cubic meters per day (m3/d) (79.8 percent) to back-barrier estuaries and about 10,200 m3/d (18.7 percent) to the ocean; about 836 m3/d (1.5 percent) may exit the system as subsea underflow. The total contribution of fresh groundwater to shoreline discharge zones amounts to about 53,700 m3/d (98.5 percent). The median age of freshwater discharged to back-barrier estuaries and the ocean was 3.4 years, and the 95th-percentile age was 20 years. The TN concentrations and loads under 2005 mean annual conditions for areas that contribute fresh groundwater to back-barrier estuaries and the ocean were calculated for the principal land uses on Fire Island. The overall TN load from the shallow aquifer system to shoreline discharge zones is about 16,200 kilograms per year (kg/yr) (82.2 percent) to back-barrier estuaries and about 3,500 kg/yr (17.8 percent) to the ocean. The overall TN load to marine surface waters amounts to about 19,700 kg/yr-roughly 6 percent of the annual TN load from shallow groundwater entering the South Shore Estuary Reserve (SSER) from the Suffolk County mainland, which is about 345,000 kg/yr. In contrast to the TN load from shallow groundwater for the SSER watershed, which annually yields about 353 kilograms per square kilometer (kg/km2), the overall TN loa

Scientific Investigations Report↗

Analysis of annual dissolved-solids loading from selected natural and irrigated catchments in the Upper Colorado River Basin, 1974-2003

Dissolved-solids loading from 17 natural catchments and 14 irrigated catchments in the Upper Colorado River Basin was examined for the period from 1974 through 2003. In general, dissolved-solids loading increased and decreased concurrently in natural and irrigated catchments but at different magnitudes. Annually, the magnitude of loading in natural catchments changed about 10 percent more, on average, than in irrigated catchments. Measures of variability, or spread, indicate that natural catchments had 35 percent greater annual variability in loading than irrigated catchments. Precipitation and dissolved-solids loads were positively correlated in natural catchments, and a weak positive correlation was determined for irrigated catchments. A weak negative correlation between temperature and dissolved-solids load was determined for both natural and irrigated catchments. In irrigated catchments, the dissolved-solids load response to an above-average precipitation period from 1982 through 1987 generally lagged behind that in the natural catchments. On average, irrigated catchments with reservoir storage had the largest normalized maximum annual loads during the wet period.

Arizona, Colorado, New Mexico, Utah, Wyoming↗

Streamflow statistics for unregulated and regulated conditions for selected locations on the Yellowstone, Tongue, and Powder Rivers, Montana, 1928-2002

Major floods in 1996 and 1997 on the Yellowstone River in Montana intensified public debate over the effects of human activities on the Yellowstone River. In 1999, the Yellowstone River Conservation District Council was formed to address conservation issues on the river. The Yellowstone River Conservation District Council partnered with the U.S. Army Corps of Engineers to conduct a cumulative-effects study on the main stem of the Yellowstone River. The cumulative-effects study is intended to provide a basis for future management decisions in the watershed. Streamflow statistics, such as flow-frequency and flow-duration data calculated for unregulated and regulated streamflow conditions, are a necessary component of the cumulative effects study. The U.S. Geological Survey, in cooperation with the Yellowstone River Conservation District Council and the U.S. Army Corps of Engineers, calculated streamflow statistics for unregulated and regulated conditions for the Yellowstone, Tongue, and Powder Rivers for the 1928–2002 study period. Unregulated streamflow represents flow conditions that might have occurred during the 1928–2002 study period if there had been no water-resources development in the Yellowstone River Basin. Regulated streamflow represents estimates of flow conditions during the 1928–2002 study period if the level of water-resources development existing in 2002 was in place during the entire study period. Peak-flow frequency estimates for regulated and unregulated streamflow were developed using methods described in Bulletin 17B. High-flow frequency and low-flow frequency data were developed for regulated and unregulated streamflows from the annual series of highest and lowest (respectively) mean flows for specified n-day consecutive periods within the calendar year. Flow-duration data, and monthly and annual streamflow characteristics, also were calculated for the unregulated and regulated streamflows.

Montana;North Dakota;Wyoming↗

Special topic—Boreholes

Introduction Installation of instrument packages in deep (several hundred to several thousand meters) boreholes near volcanoes is relatively expensive (a few million to tens of millions of U.S. dollars), but can provide a low-noise, high-quality source of geophysical (seismic, strain, tilt, and pore pressure), physical (temperature and water level), and geochemical data. Observations from instruments at depth have the potential to provide insights into processes associated with magma intrusion, unrest, and eruption that would not otherwise be possible (Lowenstern and others, 2017; Eichelberger, 2020). Examples of instrumented boreholes in volcanic areas include the 3-kilometer (km)-deep Long Valley Exploratory Well (LVEW) in California (for example, Priest and others, 1998; Prejean and Ellsworth, 2001; Fischer and others, 2003; Roeloffs and others, 2003; Sorey and others, 2003), the 1,262 meter-deep NSF Well (commonly referred to as the “Keller Well”) within the summit caldera of Kīlauea, Hawaiʻi (Keller and others, 1979; Myren and others, 2006), and the Caribbean Andesite Lava Island-volcano Precision Seismo-geodetic Observatory (CALIPSO) project at Soufrière Hills, Montserrat, which includes a series of four 200-meter (m)-deep holes (for example, Mattioli and others, 2004; Voight and others, 2006). The Plate Boundary Observatory (PBO) of the National Science Foundation’s Earthscope project placed seismometers, tiltmeters, strainmeters, and pore-pressure sensors at depths of 100 to 250 m in more than 100 boreholes scattered in western North America, including at Mount St. Helens, Washington, and Yellowstone Caldera, Wyoming. The total cost for an instrumented PBO borehole ranged from $250,000 to $270,000 U.S. dollars (USD) and a few thousand USD are required annually for maintenance (David Mencin, UNAVCO, written commun., October 2020).

Scientific Investigations Report↗

Changes in high-flow frequency and channel geometry of the Neosho River downstream from John Redmond Dam, southeastern Kansas

The streamflow regimen of the Neosho River downstream from John Redmond Dam in southeastern Kansas has changed significantly since the dam's completion in 1964. The controlled releases from the dam have decreased the magnitudes of peak discharges and increased the magnitudes of low discharges. The trends in river stage for selected discharges also have changed at two of the streamflow-gaging stations--those closest to the dam. There is a significant downward trend in the stages associated with the median annual peak discharges, but no significant trend in the stages associated with the annual mean discharges, which indicates that the river channel is increasing in width but not depth or that the hflow velocity has increased at the streamflow-gaging stations. Because there were not significant trends present in precipitation, mean annual discharge, or annual peak discharge, the changes are attributed to John Redmond Dam.

Kansas↗

Estimating the magnitude and frequency of floods at ungaged locations on urban streams in Tennessee and parts of Alabama, Georgia, Mississippi, North Carolina, and South Carolina, using data through the 2022 water year

In 2024, the U.S. Geological Survey, in cooperation with the Tennessee Department of Transportation, updated the methods for predicting the magnitude and frequency of floods at ungaged locations on streams in urban areas in Tennessee. The study area included 136 streamgages in urban areas in Tennessee, Mississippi, Alabama, Georgia, South Carolina, and North Carolina that had at least 10 percent developed imperviousness in their basins as indicated by data from the 2011 National Land Cover Database. Regression equations were developed to predict streamflows corresponding to the 50-​, 20-​, 10-​, 4-​, 2-​, 1-​, 0.5-​, and 0.2-​percent annual exceedance probabilities (AEPs) and were incorporated into the StreamStats application. In generalized least-​squares regression, the base-​10 logarithm of drainage area, the percentages of the streamgage basins in developed land use, and the percentages of the streamgage basins in the Piedmont and Ridge and Valley Level 3 ecoregions were statistically significant in explaining the variability in annual peak streamflows in the study area. Drainage areas ranged from 0.164 to 93.4 square miles, the percentage of the streamgage basins in developed land use ranged from 26 to 100 percent, and the percentage of the streamgage basins in Piedmont and Ridge and Valley Level 3 ecoregions ranged from 0 to 100 percent. Pseudo R-​squared values for the regression equations ranged from 0.86, or 86 percent, for the 50-​ and 20-​percent AEPs (2-​ and 5-​year floods) to 0.71, or 71 percent, for the 0.2-​percent AEP (500-​year flood). The average variance of prediction (in log base-​10 units) ranged from 0.023 for the 20-​ and 10-​percent AEPs to 0.05 for the 0.2-​percent AEP. The average variance of prediction can be reported as a percentage of the predicted value, known as the standard error of prediction, which ranged from 35.8 percent for the 20-​percent AEP (5-​year flood) to 55.4 percent for the 0.2-​percent AEP (500-​year flood). Methods are presented for estimating annual peak streamflows for gaged locations, ungaged locations on gaged streams, and locations on ungaged streams.

Alabama, Georgia, Mississippi, North Carolina, Sou↗

Abundance and distribution of benthic macroinvertebrates in offshore soft sediments in Western Lake Huron, 2001-2007

Invasive species have had major impacts on the Great Lakes. This is especially true of exotic dreissenid mussels which are associated with decreased abundance of native macroinvertebrates and changes in food availability for fish. Beginning in 2001, we added a benthic macroinvertebrate survey to the USGS-Great Lakes Science Center's annual fall prey fish assessment of Lake Huron to monitor abundance of macrobenthos. Mean abundance of Diporeia, the most abundant benthic taxon in Lake Huron reported by previous investigators, declined greatly between 2001 and 2007. Diporeia was virtually absent at 27-m sites by 2001, decreased and was lost completely from 46-m depths by 2006, but remained present at reduced densities at 73-m sites. Dreissenids in our samples were almost entirely quagga mussels Dreissena bugensis. Zebra mussels Dreissena polymorpha were virtually absent from our samples, suggesting that they were confined to nearshore areas shallower than we sampled. Loss of Diporeia at individual sites was associated with arrival of quagga mussels, even when mussel densities were low. Quagga mussel density peaked during 2002, then decreased thereafter. During the study quagga mussels became established at most 46-m sites, but remained rare at 73-m sites. Length frequency distributions suggest that initial widespread recruitment may have occurred during 2001-2002. Like other Great Lakes, Lake Huron quagga mussels were associated with decreased abundance of native taxa, but negative effects occurred even though dreissenid densities were much lower. Dreissenid effects may extend well into deep oligotrophic habitats of Lake Huron.

Journal of Great Lakes Research↗

Low-flow characteristics of streams in Ohio through water year 1997

This report presents selected low-flow and flow-duration characteristics for 386 sites throughout Ohio. These sites include 195 long-term continuous- record stations with streamflow data through water year 1997 (October 1 to September 30) and for 191 low-flow partial-record stations with measurements into water year 1999. The characteristics presented for the long-term continuous-record stations are minimum daily streamflow; average daily streamflow; harmonic mean flow; 1-, 7-, 30-, and 90-day minimum average low flow with 2-, 5-, 10-, 20-, and 50-year recurrence intervals; and 98-, 95-, 90-, 85-, 80-, 75-, 70-, 60-, 50-, 40-, 30-, 20-, and 10-percent daily duration flows. The characteristics presented for the low-flow partial-record stations are minimum observed streamflow; estimated 1-, 7-, 30-, and 90-day minimum average low flow with 2-, 10-, and 20-year recurrence intervals; and estimated 98-, 95-, 90-, 85- and 80-percent daily duration flows. The low-flow frequency and duration analyses were done for three seasonal periods (warm weather, May 1 to November 30; winter, December 1 to February 28/29; and autumn, September 1 to November 30), plus the annual period based on the climatic year (April 1 to March 31).

Ohio↗

Annual compilation and analysis of hydrologic data for Escondido Creek, San Antonio River basin, Texas

Introduction History of Small Watershed Projects in Texas The U.S. Soil Conservation Service is actively engaged in the installation of flood and soil erosion reducing measures in Texas under the authority of the "Flood Control Act of 1936 and 1944" and "Watershed Protection and Flood Prevention Act" (Public Law 566), as amended. The Soil Conservation Service has found a total of approximately 3,500 floodwater-retarding structures to be physically and economically feasible in Texas. As of September 30, 1970, 1,439 of these structures had been built. This watershed-development program will have varying but important effects on the surface and ground-water resources of river basins, especially where a large number of the floodwater-retarding structures are built. Basic hydrologic data under natural and developed conditions are needed to appraise the effects of the structures on the yield and mode of occurrence of runoff. Hydrologic investigations of these small watersheds were begun by the Geological Survey in 1951 and are now being made in 12 study areas (fig. 1). These investigations are being made in cooperation with the Texas Water Development Board, the Soil Conservation Service, the San Antonio River Authority, the city of Dallas, and the Tarrant County Water Control and Improvement District No. 1. The 12 study areas were chosen to sample watershed having different rainfall, topography, geology, and soils. In five of the study areas, (North, Little Elm, Mukewater, little Pond-North Elm, and Pin Oak Creeks), streamflow and rainfall records were collected prior to construction of the floodwater-retarding structures, thus affording the opportunity for analyses of the conditions "before and after" development. A summary of the development of the floodwater-retarding structures in each study areas of September 30, 1970, is shown in table 1. Objectives of the Texas Small Watersheds Project The purpose of these investigations is to collect sufficient data to meeting the following objectives: To determine the net effect of floodwater-retarding structures on the regimen of streamflow at downstream points. To determine the effectiveness of the structures as ground-water recharge facilities. To determine the effect of the structures on the sediment yield at downstream points. To develop relationships between maximum rates and/or volumes of runoff with rainfall in small natural watersheds. To develop a stream-system model for basins with floodwater-retarding structures. To determine the minimum instrumentation necessary for estimating the flood hydrographs below a system of structures, as needed for downstream water-management operation. Purpose and Scope of this Basic-Data Report This report, which is the tenth in a series of basic-data reports published annually for the Escondido Creek study area, contains the rainfall, runoff, and storage data collected during the 1970 water year for the 72.4-square-mile area above the stream-gaging station Escondido Creek at Kenedy, Texas. The location of floodwater-retarding structures and hydrologic-instrument installations in the Escondido Creek study area are shown on figure 2. This investigation is scheduled to continue through a period of both above- and below-normal precipitation to define the various factors used in the analyses of rainfall-runoff relationship. To facilitate the publication and distribution of this report at the earliest feasible time, certain material contained herein does not conform to the formal publication standards of the U.S. Geological Survey.

Texas↗

Estimation of Streamflow Characteristics for Charles M. Russell National Wildlife Refuge, Northeastern Montana

Charles M. Russell National Wildlife Refuge (CMR) encompasses about 1.1 million acres (including Fort Peck Reservoir on the Missouri River) in northeastern Montana. To ensure that sufficient streamflow remains in the tributary streams to maintain the riparian corridors, the U.S. Fish and Wildlife Service is negotiating water-rights issues with the Reserved Water Rights Compact Commission of Montana. The U.S. Geological Survey, in cooperation with the U.S. Fish and Wildlife Service, conducted a study to gage, for a short period, selected streams that cross CMR, and analyze data to estimate long-term streamflow characteristics for CMR. The long-term streamflow characteristics of primary interest include the monthly and annual 90-, 80-, 50-, and 20-percent exceedance streamflows and mean streamflows (Q.90, Q.80, Q.50, Q.20, and QM, respectively), and the 1.5-, 2-, and 2.33- year peak flows (PK1.5, PK2, and PK2.33, respectively). The Regional Adjustment Relationship (RAR) was investigated for estimating the monthly and annual Q.90, Q.80, Q.50, Q.20, and QM, and the PK1.5, PK2, and PK2.33 for the short-term CMR gaging stations (hereinafter referred to as CMR stations). The RAR was determined to provide acceptable results for estimating the long-term Q.90, Q.80, Q.50, Q.20, and QM on a monthly basis for the months of March through June, and also on an annual basis. For the months of September through January, the RAR regression equations did not provide acceptable results for any long-term streamflow characteristic. For the month of February, the RAR regression equations provided acceptable results for the long-term Q.50 and QM, but poor results for the long-term Q.90, Q.80, and Q.20. For the months of July and August, the RAR provided acceptable results for the long-term Q.50, Q.20, and QM, but poor results for the long-term Q.90 and Q.80. Estimation coefficients were developed for estimating the long-term streamflow characteristics for which the RAR did not provide acceptable results. The RAR also was determined to provide acceptable results for estimating the PK1.5., PK2, and PK2.33 for the three CMR stations that lacked suitable peak-flow records. Methods for estimating streamflow characteristics at ungaged sites also were derived. Regression analyses that relate individual streamflow characteristics to various basin and climatic characteristics for gaging stations were performed to develop regression equations to estimate streamflow characteristics at ungaged sites. Final equations for the annual Q.50, Q.20, and QM are reported. Acceptable equations also were developed for estimating QM for the months of February, March, April, June, and July, and Q.50, Q.20, and QM on an annual basis. However, equations for QM for the months of February, March, April, June, and July were determined to be less consistent and reliable than the use of estimation coefficients applied to the regression equation results for the annual QM. Acceptable regression equations also were developed for the PK1.5, PK2, and PK2.33.

Scientific Investigations Report↗

A reassessment of Chao2 estimates for population monitoring of grizzly bears in the Greater Yellowstone Ecosystem

The Yellowstone Ecosystem Subcommittee (YES) asked the Interagency Grizzly Bear Study Team (IGBST) to re-assess a technique used in annual population estimation and trend monitoring of grizzly bears in the Greater Yellowstone Ecosystem (GYE). This technique is referred to as the Chao2 approach and estimates the number of females with cubs-of-the-year (hereafter, females with cubs) and, in association with other demographic data, is used by the IGBST to produce annual population estimates. Females with cubs are an easily recognizable population segment, and trends for this reproductive segment of the population are assumed to be representative of trend for the entire population. The overarching objective of the analyses presented in this report was to provide a more accurate representation of the GYE grizzly bear population using the current methodologies in place. Specifically, we addressed two limitations of the current Chao2 approach: 1) underestimation bias associated with a distance criterion used to differentiate annual sightings of females with cubs into unique individuals and 2) limitations of the model-averaging approach to effectively distinguish among potential future population trajectories (decline, stability, and growth). The first issue addressed in this report is the underestimation bias associated with the rule set that Knight et al. (1995) developed to differentiate sightings of females with cubs into unique individuals (i.e., unique family groups). The rule set was originally designed to be conservative by reducing the risk of identifying more females with cubs than actually existed, primarily through use of a distance criterion of 30 km to separate sightings of unique females. This approach resulted in an underestimation bias, and previous research demonstrated that this bias increases with increasing number of females with cubs. Using location data from radio-marked females with cubs, we evaluated alternative distance criteria by simulating scenarios with varying numbers of true females with cubs and sightings. Findings from these analyses demonstrate that bias in estimates of females with cubs can be substantially reduced by changing the 30-km distance criterion in the rule set to 16 km, which produced relatively unbiased estimates. Findings also indicate, however, the importance of adaptability with regard to the distance criteria because of the complex relationships and biases among the various parameters involved in estimation of unique females with cubs. The total number of annual sightings and the true number of females with cubs play particularly important roles. Whereas these analyses remind us that there is no perfect approach to estimating the number of females with cubs from sightings under various scenarios, they provide us with new tools to determine when and how to adapt the monitoring program. The second issue we were tasked to investigate was the potential for improvement of the technique referred to as model-averaging, which serves to smooth relatively high variation in annual estimates. This technique was chosen by YES as the basis for monitoring the Yellowstone grizzly bear population, as described in the 2016 Conservation Strategy. This choice was made in part because the technique has been well documented and population estimates derived from counts of females with cubs are conservative. Using simulations of population trends, we demonstrate why the model-averaging technique currently used cannot distinguish between plausible future trend scenarios. As a suitable alternative to model averaging, we propose the use of generalized additive models (GAMs). Using a suite of simulated trend dynamics relevant to management, we demonstrate GAM performance for tracking trends in females with cubs within the context of the annual monitoring program. We demonstrate the ability to not only document directional changes in population trend but also patterns of stabilization or resiliency after such changes. Furthermore, the proposed monitoring framework provides objective measures useful for early detection of directional changes in trend. The new framework is flexible, allowing retrospective analysis of Chao2-based estimates and future applications to time series of other population metrics, such as vital rates. The aforementioned updates provide us with new tools to determine when and how to adapt the monitoring program. Within the context of current monitoring protocols and effort, and considering the full suite of simulations presented in this report and previous studies, the IGBST plans to incorporate the following changes to the population monitoring protocol: 1) modify the distance criterion, starting with 16 km under current sampling conditions and 2) revise the population monitoring framework using GAMs as the basis for smoothing of annual estimates and detecting trends and changes in trend. Implementation of the 16-km distance criterion combined with use of GAM techniques would affect some of the population metrics (e.g., annual population size and uncertainty, population trend, mortality rates) used to inform management responses. A primary consideration is that the 16-km distance criterion results in total population estimates derived from the Chao2 estimates that are greater than those we have reported in the past. This increase is due to a change in the implementation of the technique and more accurately represents the number of females with cubs in the GYE grizzly bear population. Additionally, interpretation of retrospective trend patterns may change due to the combination of a different distance criterion and enhanced trend monitoring based on the GAM approach we present here. Implementation will require relatively minor changes in the monitoring protocols described in Appendices B and C of the 2016 Conservation Strategy. Finally, we note that the IGBST has ongoing investigations into the merits of an Integrated Population Model (IPM), for which annual Chao2-based estimates are important input data. The IGBST plans to continue those investigations using the 16-km distance criterion to derive Chao2 estimates.

Idaho, Montana, Wyoming↗

Rio Grande transboundary integrated hydrologic model and water-availability analysis, New Mexico and Texas, United States, and Northern Chihuahua, Mexico

Errata **September 28, 2018: The purpose of a USGS Open-file report (OFR) is dissemination of information that must be released immediately to fill a public need or information that is not sufficiently refined to warrant publication in one of the other USGS series. As part of that refinement process, an error was discovered in one of the input data sets of the Rio Grande Transboundary Integrated Hydrologic Model (RGTIHM) that this OFR was based upon. The error involved the assignment of storage properties to “phantom cells.” Phantom cells are required for most variants of MODFLOW that use a structured finite-difference grid when individual stratigraphic layers are represented as separate layers. Using phantom cells is a common practice that allows separate model layers to be maintained without having to combine stratigraphic layers into equivalent model layers or to use an unstructured grid. Typically, phantom cell horizontal hydraulic conductivities and storage properties are set to a small number and vertical hydraulic conductivities are set to a number large enough to allow vertical flow between the vertically adjacent layers. In the RGTIHM, the specific storage properties of the phantom cells for the upper (RGTIHM layers 3 and 4), middle (RGTIHM layers 5 and 6), and lower (RGTIHM layers 7 and 8) members of the Santa Fe Group were inadvertently assigned a value of 1 feet -1 . The revision of these specific storage values to a small number (1.0 x 10 -09 feet -1 ) required additional trial-and-error model calibration and a new sensitivity analysis. After calibration, the overall model fit remained similar to the fit described in the OFR, but the fit for many individual features such as project water available for diversions at the American Canal and Acequia Madre improved due to the reduction in flow coming from lower layers. Overall, there is still an average net depletion of groundwater flow, and the conclusions of the report are not changed. The revised average annual groundwater flow depletion simulated for the period 1953-2014 is -1,480 acre-feet/year for the entire model region, and -3,660 acre-feet/year for the portion of the model in the United States. The final version of the model will be the basis of the USGS Scientific Investigations Report that will supersede this OFR. An updated Model Archive of RGTIHM is available upon request to the USGS California Water Science Center. The corrected version of the model WAS the basis for the USGS Scientific Investigations Report that SUPERSEDED this Open-File Report.** Abstract Changes in population, agricultural development and practices (including shifts to more water-intensive crops), and climate variability are increasing demands on available water resources, particularly groundwater, in one of the most productive agricultural regions in the Southwest—the Rincon and Mesilla Valley parts of Rio Grande Valley, Doña Ana and Sierra Counties, New Mexico, and El Paso County, Texas. The goal of this study was to produce an integrated hydrological simulation model to help evaluate water-management strategies, including conjunctive use of surface water and groundwater for historical conditions, and to support long-term planning for the Rio Grande Project. This report describes model construction and applications by the U.S. Geological Survey, working in cooperation and collaboration with the Bureau of Reclamation. This model, the Rio Grande Transboundary Integrated Hydrologic Model, simulates the most important natural and human components of the hydrologic system, including selected components related to variations in climate, thereby providing a reliable assessment of surface-water and groundwater conditions and processes that can inform water users and help improve planning for future conditions and sustained operations of the Rio Grande Project (RGP) by the Bureau of Reclamation. Model development included a revision of the conceptual model of the flow system, construction of a Transboundary Rio Grande Watershed Model (TRGWM) water-balance model using the Basin Characterization Model (BCM), and construction of an integrated hydrologic flow model with MODFLOW-One-Water Hydrologic Flow Model (referred to as One Water). The hydrologic models were developed for and calibrated to historical conditions of water and land use, and parameters were adjusted so that simulated values closely matched available measurements (calibration). The calibrated model was then used to assess the use and movement of water in the Rincon Valley, Mesilla Basin, and northern part of the Conejos-Médanos Basin, with the entire region referred to as the “Transboundary Rio Grande” or TRG. These tools provide a means to understand hydrologic system response to the evolution of water use in the region, its availability, and potential operational constraints of the RGP. The conceptual model identified surface-water and groundwater inflows and outflows that included the movement and use of water both in natural and in anthropogenic systems. The groundwater-flow system is characterized by a layered geologic sedimentary sequence combined with the effects of groundwater pumping, operation of the RGP, natural runoff and recharge, and the application of irrigation water at the land surface that is captured and reused in an extensive network of canals and drains as part of the conjunctive use of water in the region. Historical groundwater-level fluctuations followed a cyclic pattern that were aligned with climate cycles, which collectively resulted in alternating periods of wet or dry years. Periods of drought that persisted for one or more years are associated with low surface-water availability that resulted in higher rates of groundwater-level decline. Rates of groundwater-level decline also increased during periods of agricultural intensification, which necessitated increasing use of groundwater as a source of irrigation water. Agriculture in the area was initially dominated by alfalfa and cotton, but since 1970 more water-intensive pecan orchards and vegetable production have become more common. Groundwater levels substantially declined in subregions where drier climate combined with increased demand, resulting in periods of reduced streamflows. Most of the groundwater was recharged in the Rio Grande Valley floor, and most of the pumpage and aquifer storage depletion was in Mesilla Basin agricultural subregions. A cyclic imbalance between inflows and outflows resulted in the modeled cyclic depletion (groundwater withdrawals in excess of natural recharge) of the groundwater basin during the 75-year simulation period of 1940–2014. Changes in groundwater storage can vary considerably from year to year, depending on land use, pumpage, and climate conditions. Climatic drivers of wet and dry years can greatly affect all inflows, outflows, and water use. Although streamflow and, to a minor extent, precipitation during inter-decadal wet-year periods replenished the groundwater historically, contemporary water use and storage depletion could have reduced the effects of these major recharge events. The average net groundwater flow-rate deficit for 1953–2014 was estimated to be about 8,990 acre-feet per year.

New Mexico, Northern Chihuahua, Texas↗

Magnitude, frequency, and trends of floods at gaged and ungaged sites in Washington, based on data through water year 2014

An investigation into the magnitude and frequency of floods in Washington State computed the annual exceedance probability (AEP) statistics for 648 U.S. Geological Survey unregulated streamgages in and near the borders of Washington using the recorded annual peak flows through water year 2014. This is an updated report from a previous report published in 1998 that used annual peak flows through the water year 1996. New in this report, a regional skew coefficient was developed for the Pacific Northwest region that includes areas in Oregon, Washington, Idaho and western Montana within the Columbia River drainage basin south of the United States-Canada border, the coastal areas of Oregon and western Washington, and watersheds draining into Puget Sound, Washington. The skew coefficient is an important term in the Log Pearson Type III equation used to define the distribution of the log-transformed annual peaks. The Expected Moments Algorithm was used to fit historical and censored peak-flow data to the log Pearson Type III distribution. A Multiple Grubb-Beck test was employed to censor low outliers of annual peak flows to improve on the frequency distribution. This investigation also includes a section on observed trends in annual peak flows that showed significant trends ( p -value < 0.05) in 21 of 83 long-term sites, but with small magnitude Kendall tau values suggesting a limited monotonic trend in the time series of annual peaks. Most of the sites with a significant trend in western Washington were positive and all the sites with significant trends (three sites) in eastern Washington were negative. Multivariate regression analysis with measured basin characteristics and the AEP statistics at long-term, unregulated, and un-urbanized (defined as drainage basins with less than 5 percent impervious land cover for this investigation) streamgages within Washington and some in Idaho and Oregon that are near the Washington border was used to develop equations to estimate AEP statistics at ungaged basins. Washington was divided into four regions to improve the accuracy of the regression equations; a set of equations for eight selected AEPs and for each region were constructed. Selected AEP statistics included the annual peak flows that equaled or exceeded 50, 20, 10, 4, 2, 1, 0.5 and 0.2 percent of the time equivalent to peak flows for peaks with a 2-, 5-, 10-, 25-, 50-, 100-, 200-, and 500-year recurrence intervals, respectively. Annual precipitation and drainage area were the significant basin characteristics in the regression equations for all four regression regions in Washington and forest cover was significant for the two regression regions in eastern Washington. Average standard error of prediction for the regional regression equations ranged from 70.19 to 125.72 percent for Regression Regions 1 and 2 on the eastern side of the Cascade Mountains and from 43.22 to 58.04 percent for Regression Regions 3 and 4 on the western side of the Cascade Mountains. The pseudo coefficient of determination (where a value of 100 signifies a perfect regression model) ranged from 68.39 to 90.68 for Regression Regions 1 and 2, and 92.35 to 95.44 for Regions 3 and 4. The calculated AEP statistics for the streamgages and the regional regression equations are expected to be incorporated into StreamStats after the publication of this report. StreamStats is the interactive Web-based map tool created by the U.S. Geological Survey to allow the user to choose a streamgage and obtain published statistics or choose ungaged locations where the program automatically applies the regional regression equations and computes the estimates of the AEP statistics.

Washington↗

Precipitation data for burial grounds 5 and 6, Oak Ridge National Laboratory, Tennessee, 1976-1980

As part of a hydrogeologic investigation, precipitation data were collected at two stations, one each in Burial Grounds 5 and 6 at the Oak Ridge National Laboratory, Tennessee. Daily, monthly, and annual values are reported herein for the period from January 1976 through December 1980. During this period, annual values ranged from about 25 percent above to about 25 percent below the calculated mean of 51.96 inches at Burial Ground 5 and 49.60 inches at Burial Ground 6.

Tennessee↗