Search USGS⌕ Search

SEARCH · Search USGS

Results for “Surface Water Techniques”

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 721 records · Page 40Linked to original sources

Integrating hydrologic and geophysical data to constrain coastal surficial aquifer processes at multiple spatial and temporal scales

Since 1997, repeated, coincident geophysical surveys and extensive hydrologic studies in shallow monitoring wells have been used to study static and dynamic processes associated with surface water-groundwater interaction at a range of spatial scales at the estuarine and ocean boundaries of an undeveloped, permeable barrier island in the Georgia part of the U.S. South Atlantic Bight. Because geophysical and hydrologic data measure different parameters, at different resolution and precision, and over vastly different spatial scales, reconciling the coincident data or even combining complementary inversion, hydrogeochemcial analyses and well-based groundwater monitoring, and, in some cases, limited vegetation mapping to demonstrate the utility of an integrative, multidisciplinary approach for elucidating groundwater processes at spatial scales (tens to thousands of meters) that are often difficult to capture with traditional hydrologic approaches. The case studies highlight regional aquifer characteristics, varying degrees of lateral saltwater intrusion at estuarine boundaries, complex subsurface salinity gradients at the ocean boundary, and imaging of submarsh groundwater discharge and possible free convection in the pore waters of a clastic marsh. This study also documents the use of geophysical techniques for detecting temporal changes in groundwater salinity regimes under natural (not forced) gradients at intratidal to interannual (1998-200 Southeastern U.S.A. drought) time scales.

Book chapter↗

Phreatophytic land-cover map of the northern and central Great Basin Ecoregion: California, Idaho, Nevada, Utah, Oregon, and Wyoming

Increasing water use and changing climate in the Great Basin of the western United States are likely affecting the distribution of phreatophytic vegetation in the region. Phreatophytic plant communities that depend on groundwater are susceptible to natural and anthropogenic changes to hydrologic flow systems. The purpose of this report is to document the methods used to create the accompanying map that delineates areas of the Great Basin that have the greatest potential to support phreatophytic vegetation. Several data sets were used to develop the data displayed on the map, including Shrub Map (a land-cover data set derived from the Regional Gap Analysis Program) and Gap Analysis Program (GAP) data sets for California and Wyoming. In addition, the analysis used the surface landforms from the U.S. Geological Survey (USGS) Global Ecosystems Mapping Project data to delineate regions of the study area based on topographic relief that are most favorable to support phreatophytic vegetation. Using spatial analysis techniques in a GIS, phreatophytic vegetation classes identified within Shrub Map and GAP were selected and compared to the spatial distribution of selected landforms in the study area to delineate areas of phreatophyte vegetation. Results were compared to more detailed studies conducted in selected areas. A general qualitative description of the data and the limitations of the base data determined that these results provide a regional overview but are not intended for localized studies or as a substitute for detailed field analysis. The map is intended as a decision-support aide for land managers to better understand, anticipate, and respond to ecosystem changes in the Great Basin.

Scientific Investigations Map↗

Modeling, mapping, and measuring the risk of freshwater invasive species across Alaska

Freshwater ecosystems of the Alaskan Arctic and Subarctic provide resources that are culturally, ecologically, and economically invaluable. Presently, these regions are relatively free of the impacts from invasive species compared to southern latitudes. To date, there have been relatively few verified introductions of aquatic invasive species (AIS) to freshwater ecosystems in Alaska. The expanding list and distribution of AIS has led to significant negative ecological and economic impacts (e.g., waterweed Elodea nuttalli ; E. canadensis and northern pike Esox Lucius introduced outside its native range in Alaska). Escalating human activity across Alaskan lands and waters, coupled with rapidly shifting environmental conditions, increases the potential for new species introductions and subsequent establishment. Creating a proactive framework for well-informed decision-making and action can improve the effectiveness of prevention efforts and bolster decision support tools that help resource managers direct limited resources. Prioritizing AIS that may be introduced and become established, as well as the locations at highest risk of invasion, is foundational to building a proactive invasive species management framework in Alaska. This project sought to identify and prioritize AIS known to be invasive in the contiguous United States, evaluate current and future habitat suitability for AIS in Alaska, and assess potential for AIS to be transported to habitats across Alaska, utilizing similar assessment methods as implemented for Bering Sea marine invasive species and non-native plants in Alaska. To accomplish this goal, the objectives of the project were to: 1) develop a formal ranked list of potential AIS to freshwater systems of Alaska; 2) assess the level of establishment risk for potential AIS by developing habitat suitability models for waterbodies across Alaska; and 3), identify potential pathways and specific vectors for high-risk AIS to invade Alaska and develop a framework for how vector analysis will be completed to understand transport risk. Overall, our goal is horizon scanning which is defined by Roy et al. (2019) as “a systematic examination of potential threats and opportunities, within a given context, and likely future developments, which are at the margin of current thinking and planning.” The scans include pathway analyses and risk screening of species present at pathway origin points, with a focus on identifying species at high risk of being introduced, becoming established, spreading, and causing harm. We refined a list of 28 AIS from a list of hundreds based on characterizations of species’ invasiveness and species’ proximity to Alaska (USGS 2020; GBIF 2022). Next, we evaluated the relative invasiveness of individual species to create an initial AIS ranking. We sought to characterize habitat suitability of AIS by selecting variables that were continental in scale, covering North America to include Alaska as well as the lower 48 states comparing natural discharge, sub-basin average terrain slope (degrees), average silt fraction, average organic carbon, lithological class, and human footprint in sub-basin in 2009. We estimated AIS habitat suitability across the entire state of Alaska using the physiological tolerances of the AIS (Appendix 2). We also evaluated pathways and vectors for the introduction of AIS (Appendix 2). Many pathways and vectors considered did not meet the criteria for Alaska or freshwater systems. Of the 28 ranked species that we categorized as very high, high, and moderate levels of invasiveness; all three risk groups included fish and mollusks (Appendix 2). One commonality of the very high-invasiveness-ranked species was the availability of Ecological Risk Screening Summary documents (USFWS, 2022) produced by U.S. Fish and Wildlife Service (USFWS), except for the goldfish ( Carassius auratus ) and the New Zealand mudsnail ( Potamopyrgus antipodarum ). The Ecological Risk Screening Summary is now available for New Zealand mudsnails. In general, fish species often ranked very high or high in invasiveness and included sportfish and aquarium fish, suggesting the importance of pathways such as aquarium trade, fishing industry, intentional (but illegal) introductions of sportfishes and aquarium fishes for establishment. The technique we used for habitat suitability models necessitated aquatic environmental datasets that were continental in scale, which was often interpolated from very coarse resolution source data layers, particularly in Alaska. Better spatial data representing aquatic environments would likely improve this approach. While the lack of introductions in Alaska and nearby provinces and states is encouraging, the lack of occurrence data for the focal species also created complications for habitat suitability modeling. Despite the challenges, the habitat suitability models indicated limited suitability for warmwater species while some species, such as Brook trout (S alvelinus fontinalis ), have high habitat suitability across Alaska no matter what threshold approach is taken. Some environmental predictors were more important than others. Specifically, the most important predictor variable, ‘frost free days,’ was critical for 15 out of 28 species as expected due to harsh winter conditions in Arctic and Subarctic regions. The second most important predictor was ‘subbasin land surface runoff’, a variable that indicates the amount of discharge and runoff, while the third most important predictor was ‘snow cover’ another indication of winter conditions. Overall, the ability to understand the effect of future climate scenarios on the establishment of AIS was challenging. A detailed dataset of freshwater temperatures and water chemistry (e.g., pH, calcium) would greatly improve the ability to predict invasiveness of freshwater species to Alaska’s ecosystems on a regional basis. Future studies may benefit from a more focused geographic scope examining a group of subbasins or a regional basin rather than the entire state. These drainages could be selected based upon the mostly likely locations of introduction pathways. The two most prevalent pathway risks for AIS are in-state transfer and stowaways/contaminants. Although there are examples of introductions from other pathways, the risk is somewhat mitigated by Alaska’s climate and regulations. However, variable application of protocols for inspection and cleaning of fishing gear, watercraft, and other similar items while traveling into Alaska as well as transferring from waterbody to waterbody within the state creates a substantial risk in introducing invasive species. We plot cumulative invasive vulnerability for all subbasins and for the top 10% of subbasins (Appendix 3).

Alaska↗

Mineralogical, magnetic and geochemical data constrain the pathways and extent of weathering of mineralized sedimentary rocks

The oxidative weathering of sulfidic rock can profoundly impact watersheds through the resulting export of acidity and metals. Weathering leaves a record of mineral transformation, particularly involving minor redox-sensitive phases, that can inform the development of conceptual and quantitative models. In sulfidic sedimentary rocks, however, variations in depositional history, diagenesis and mineralization can change or overprint the distributions of these trace minerals, complicating the interpretation of weathering signatures. Here we show that a combination of bulk mineralogical and geochemical techniques, micrometer-resolution X-ray fluorescence microprobe analysis and rock magnetic measurements, applied to drill core samples and single weathered fractures, can provide data that enable the development of a geochemically consistent weathering model. This work focused on one watershed in the Upper Colorado River Basin sitting within the Mesaverde Formation, a sedimentary sandstone bedrock with disseminated sulfide minerals, including pyrite and sphalerite, that were introduced during diagenesis and subsequent magmatic-hydrothermal mineralization. Combined analytical methods revealed the pathways of iron (Fe), carbonate and silicate mineral weathering and showed how pH controls element retention or release from the actively weathering fractured sandstone. Drill core logging, whole rock X-ray diffraction, and geochemical measurements document the progression from unweathered rock at depth to weathered rock at the surface. X-ray microprobe analyses of a 1-cm size weathering profile along a fracture surface are consistent with the mobilization of Fe(II) and Fe(III) into acidic pore water from the dissolution of primary pyrite, Fe-sphalerite, chlorite, and minor siderite and pyrrhotite. These reactions are followed by the precipitation of secondary minerals such as of goethite and jarosite, a Fe-(oxyhydr)oxide and hydrous Fe(III) sulfate, respectively. Microscale analyses also helped explain the weathering reactions responsible for the mineralogical transformations observed in the top and most weathered section of the drill core. For example, dissolution of feldspar and chlorite neutralizes the acidity generated by Fe and sulfide mineral oxidation, oversaturating the solution in both Fe-oxides. The combination of X-ray spectromicroscopy and magnetic measurements show that the Fe(III) product is goethite, mainly present either as a coatings on fracture surfaces in the actively weathering region of the core or more homogeneously contained within the unconsolidated regolith at the top of the core. Low-temperature magnetic data reveal the presence of ferromagnetic Fe-sulfide pyrrhotite that, although it occurs at trace concentrations, could provide a qualitative proxy for unweathered sulfide minerals because the loss of pyrrhotite is associated with the onset of oxidative weathering. Pyrrhotite loss and goethite formation are detectable through room-temperature magnetic coercivity changes, suggesting that rock magnetic measurements can determine weathering intensity in rock samples at many scales. This work contributes evidence that the weathering of sulfidic sedimentary rocks follows a geochemical pattern in which the abundance of sulfide minerals controls the generation of acidity and dissolved elements, and the pH-dependent mobility of these elements controls their export to the ground- and surface-water.

Colorado↗

Establishing aquatic restoration priorities using a watershed approach

Since the passage of the Clean Water Act in 1972, the United States has made great strides to reduce the threats to its rivers, lakes, and wetlands from pollution. However, despite our obvious successes, nearly half of the nation’s surface water resources remain incapable of supporting basic aquatic values or maintaining water quality adequate for recreational swimming. The Clean Water Act established a significant federal presence in water quality regulation by controlling point and non-point sources of pollution. Point-sources of pollution were the major emphasis of the Act, but Section 208 specifically addressed non-point sources of pollution and designated silviculture and livestock grazing as sources of non-point pollution. Non-point source pollutants include runoff from agriculture, municipalities, timber harvesting, mining, and livestock grazing. Non-point source pollution now accounts for more than half of the United States water quality impairments. To successfully improve water quality, restoration practitioners must start with an understanding of what ecosystem processes are operating in the watershed and how they have been affected by outside variables. A watershed-based analysis template developed in the Pacific Northwest can be a valuable aid in developing that level of understanding. The watershed analysis technique identifies four ecosystem scales useful to identify stream restoration priorities: region, basin, watershed, and site. The watershed analysis technique is based on a set of technically rigorous and defensible procedures designed to provide information on what processes are active at the watershed scale, how those processes are distributed in time and space. They help describe what the current upland and riparian conditions of the watershed are and how these conditions in turn influence aquatic habitat and other beneficial uses. The analysis is organized as a set of six steps that direct an interdisciplinary team of specialists to examine the biotic and abiotic processes influencing aquatic habitat and species abundance. This process helps develop an understanding of the watershed within the context of the larger ecosystem. The understanding gained can then be used to identify and prioritize aquatic restoration activities at the appropriate temporal and spatial scale. The watershed approach prevents relying solely on site-level information, a common problem with historic restoration efforts. When the watershed analysis process was used in the Whitefish Mountains of northwest Montana, natural resource professionals were able to determine the dominant habitat forming processes important for native fishes and use that information to prioritize, plan, and implement the appropriate restoration activities at the watershed scale. Despite considerable investments of time and resources needed to complete an analysis at the watershed scale, the results can prevent the misdiagnosis of aquatic problems and help ensure that the objectives of aquatic restoration will be met.

Journal of Environmental Management↗

Stratified deposits of the oxides and carbonates of manganese

Compared with the stratified deposits of iron minerals, those of manganese minerals have received little attention until recent years. Before 1930, students of the stratified deposits of manganese minerals have generally concluded that the contained manganese was derived from the decomposition of the rocks that formed the borders of the basins. Only in a few places have geologists recognized that if these rocks were the principal source of the contained manganese the basins should contain also enormous quantities of iron minerals, whereas most of the large deposits of manganese contain little iron. Even though more than 100 years ago some geologists proposed that the iron contained in some deposits of stratified iron minerals was probably derived from hydrothermal waters related to centers of volcanism, not until about 1930 was this source of manganese in similar deposits seriously proposed. This mode of origin was given the name "volcanogene sedimentaire" by French geologists working in Morocco. Since then, other names, such as "exhalative sedimentaire," have been used by geologists working in several European districts. Because of the development of new techniques for the study of the chemical and physical features of the minerals, and because close attention to the lithologic environments of the beds is being given more and more during recent years, it has seemed advisable to review the features of stratified deposits of the manganese minerals in many parts of the world and over a wide range in age. This study indicates that at least three sources may have contributed the manganese in the large deposits of iron-free oxides and carbonates : (1) the rocks that form the borders of the basins, either marine or continental; (2) the nearby underlying sediments, largely of igneous origin, decomposed by warm waters largely derived from depth; and (3) waters of hydrothermal origin derived from great depths during epochs of volcanism from which iron minerals with little manganese are deposited in deep zones, then minerals with much iron and more manganese at intermediate depths, and, finally, manganese minerals with little iron near the surface. As the hot waters of many thermal springs contain more manganese than iron, such waters could yield the pure manganese oxides and carbonates found in stratified deposits . Several kinds of evidence indicate that most of the manganese in the large stratified deposits of the oxides and carbonates of manganese in many parts of the world has been derived from hydrothermal waters from depth related to centers of volcanism. Obviously, manganese derived from decay of the rocks on the lands adjacent to large basins-marine and continental -may have been added to that derived from centers of volcanism to form the sedimentary deposits found in the basins.

Economic Geology↗

An open-water electrical geophysical tool for mapping sub-seafloor heavy placer minerals in 3D and migrating hydrocarbon plumes in 4D

A towed-streamer technology has been developed for mapping placer heavy minerals and dispersed hydrocarbon plumes in the open ocean. The approach uses induced polarization (IP), an electrical measurement that encompasses several different surface-reactive capacitive and electrochemical phenomena, and thus is ideally suited for mapping dispersed or disseminated targets. The application is operated at sea by towing active electrical geophysical streamers behind a ship; a wide area can be covered in three dimensions by folding tow-paths over each other in lawn-mower fashion. This technology has already been proven in laboratory and ocean settings to detect IP-reactive titanium- and rare-earth (REE) minerals such as ilmenite and monazite. By extension, minerals that weather and accumulate/concentrate by a similar mechanism, including gold, platinum, and diamonds, may be rapidly detected and mapped indirectly- even when dispersed and covered with thick, inert sediment. IP is also highly reactive to metal structures such as pipelines and cables. Currently, the only means for mapping an oil-spill plume is to park a large ship in the ocean and drop a sampling string over the side, requiring hours of time per sampling point. The samples must then be chemically analyzed, adding additional time and expense. We believe that an extension of the marine IP technology could also apply to rapidly mapping both seafloor- blanket and disseminated hydrocarbon plumes in the open ocean, as hydrocarbon droplets in conductive seawater are topologically equivalent to a metal-plates-and-dielectric capacitor. Because the effective capacitance would be frequency-dependent on droplet size, the approach we advocate holds the potential to not only map, but also to characterize the evolution and degradation of such a plume over time. In areas where offshore oil field development has been practiced for extended periods, making IP measurements from a towed streamer may be useful for locating buried - nd exposed pipelines, as well as pipeline leaks. We believe this technique will be a more cost-effective method than drop-sampling to map and monitor hydrocarbon plumes in open ocean settings. A marine induced polarization system was used successfully to map a 15 km × 45 km swath of the ocean floor off eastern South Africa with 3-meter sampling along 200-meter-separated profiles. The survey detected titanium-bearing sands up to 15 meters below the seafloor. From preliminary laboratory work it is apparent that we can extend this technology to monitor significant environmental problems including anthropogenic and industrial waste washed into sensitive estuaries and sounds during storm-water runoff episodes, and also to map and characterize dispersed oil plumes in the seawater column in three dimensions, as well as movement and dispersal of both over time.

Conference Paper↗

RIP-ET: A riparian evapotranspiration package for MODFLOW-2005

A new evapotranspiration package for the U.S. Geological Survey's groundwater-flow model, MODFLOW, is documented. The Riparian Evapotranspiration Package (RIP-ET) provides flexibility in simulating riparian and wetland transpiration not provided by the Evapotranspiration (EVT) or Segmented Function Evapotranspiration (ETS1) Packages for MODFLOW 2005. This report describes how the RIP-ET package was conceptualized and provides input instructions, listings and explanations of the source code, and an example. Traditional approaches to modeling evapotranspiration (ET) processes assume a piecewise linear relationship between ET flux and hydraulic head. The RIP-ET replaces this traditional relationship with a segmented, nonlinear dimensionless curve that reflects the eco-physiology of riparian and wetland ecosystems. Evapotranspiration losses from these ecosystems are dependent not only on hydraulic head, but on the plant types present. User-defined plant functional groups (PFGs) are used to elucidate the interaction between plant transpiration and groundwater conditions. Five generalized plant functional groups based on transpiration rates, plant rooting depth, and water tolerance ranges are presented: obligate wetland, shallow-rooted riparian, deep-rooted riparian, transitional riparian and bare ground/open water. Plant functional groups can be further divided into subgroups (PFSGs) based on plant size, density or other characteristics. The RIP-ET allows for partial habitat coverage and mixtures of plant functional subgroups to be present in a single model cell. RIP-ET also distinguishes between plant transpiration and bare-ground evaporation. Habitat areas are designated by polygons; each polygon can contain a mixture of PFSGs and bare ground, and is assigned a surface elevation. This process requires a determination of fractional coverage for each of the plant functional subgroups present in a polygon to account for the mixture of coverage types and resulting transpiration. The fractional cover within a cell has two components: (1) the polygonal fraction of active habitat (excluding area of bare ground, dead trees, or brush) in a cell, and (2) fraction of plant type area or bare ground area in a polygon. RIP-ET determines the transpiration rate for each plant functional group and evaporation from bare ground/open water in a cell, the total ET in the cell, and the total ET rate over the region of simulation.

Techniques and Methods↗

Faecal indicator bacteria enumeration in beach sand: A comparison study of extraction methods in medium to coarse sands

Aims: The absence of standardized methods for quantifying faecal indicator bacteria (FIB) in sand hinders comparison of results across studies. The purpose of the study was to compare methods for extraction of faecal bacteria from sands and recommend a standardized extraction technique. Methods and Results: Twenty-two methods of extracting enterococci and Escherichia coli from sand were evaluated, including multiple permutations of hand shaking, mechanical shaking, blending, sonication, number of rinses, settling time, eluant-to-sand ratio, eluant composition, prefiltration and type of decantation. Tests were performed on sands from California, Florida and Lake Michigan. Most extraction parameters did not significantly affect bacterial enumeration. anova revealed significant effects of eluant composition and blending; with both sodium metaphosphate buffer and blending producing reduced counts. Conclusions: The simplest extraction method that produced the highest FIB recoveries consisted of 2 min of hand shaking in phosphate-buffered saline or deionized water, a 30-s settling time, one-rinse step and a 10 : 1 eluant volume to sand weight ratio. This result was consistent across the sand compositions tested in this study but could vary for other sand types. Significance and Impact of the Study: Method standardization will improve the understanding of how sands affect surface water quality. ?? 2009 The Society for Applied Microbiology.

Journal of Applied Microbiology↗

Observation of the geology and geomorphology of the 1999 Marsokhod test site

The Marsokhod rover returned data from six stations that were used to decipher the geomorphology and geology of a region not previously visited by members of the geomorphology field team. Satellite images and simulated descent images provided information about the regional setting. The landing zone was on an alluvial apron flanking a mountain block to the west and a playa surface to the east. Rover color images, infrared spectra analysis of the mountains, and the apron surface provided insight into the rock composition of the nearby mountains. From the return data the geomorphology team interpreted the region to consist of compressionally deformed, ancient marine sediments and igneous rocks exposed by more recent extensional tectonics. Unconsolidated alluvial materials blanket the lower flanks of the mountains. An ancient shoreline cut into alluvial material marks a high stand of water during a past, wetter climate period. Playa sediments floor a present-day, seasonally, dry lake. Observations made by the rover using panoramic and close-up (hand specimens-scale) image data and color scene data confirmed the presence of boulders, cobbles, and fines of various provinces. Rover traverses to sites identified as geologically distinct, such as fan, channel, shoreline, and playa, provided useful clues to the geologic interpretations. Analysis of local rocks was given context only through comparison with distant geologic features. These results demonstrated the importance of a multifaceted approach to site interpretation through comparison of interpretations derived by differing geologic techniques.

California↗

Development and evaluation of a boat-mounted RFID antenna for monitoring freshwater mussels

Development of radio frequency identification (RFID) technology and passive integrated transponder (PIT) tags has substantially increased the ability of researchers and managers to monitor populations of aquatic organisms. However, use of transportable RFID antenna systems (i.e., backpack-mounted) is currently limited to wadeable aquatic environments (<1.4 m water depth). We describe the design, construction, and evaluation of a boat-mounted RFID antenna to detect individually PIT-tagged benthic aquatic organisms (mussels). We evaluated the effects of tag orientation on detection distances in water with a 32-mm half-duplex PIT tag. Detection distances up to 50 cm from the antenna coils were obtained, but detection distance was dependent on tag orientation. We also evaluated detection distance of PIT tags beneath the sediment to simulate detection of burrowing mussels with 23- and 32-mm tags. In sand substrate, the maximum detection distance varied from 3.5 cm and 4.5 cm (vertical tag orientation) to 24.7 cm and 39.4 cm (45&deg; tag orientation) for the 23- and 32-mm PIT tags, respectively. Our results suggest a 1.4-m total detection width for tagged mussels on the substrate surface by the boat-mounted antenna system regardless of tag orientation. However, burrowed mussels may require multiple passes to increase detection that would be influenced by depth, tag orientation, and tag size. Construction of the boat-mounted antenna was relatively low in cost (<500 USD) and had several advantages (less labor and time intensive, increased safety) over traditional mussel sampling techniques (diving, snorkeling) in nonwadeable habitats.

Freshwater Science↗

Hydrogeology, hydrologic effects of development, and simulation of groundwater flow in the Borrego Valley, San Diego County, California

Executive Summary The Borrego Valley is a small valley (110 square miles) in the northeastern part of San Diego County, California. Although the valley is about 60 miles northeast of city of San Diego, it is separated from the Pacific Ocean coast by the mountains to the west and is mostly within the boundaries of Anza-Borrego Desert State Park. From the time the basin was first settled, groundwater has been the only source of water to the valley. Groundwater is used for agricultural, recreational, and municipal purposes. Over time, groundwater withdrawal through pumping has exceeded the amount of water that has been replenished, causing groundwater-level declines of more than 100 feet in some parts of the basin. Continued pumping has resulted in an increase in pumping lifts, reduced well efficiency, dry wells, changes in water quality, and loss of natural groundwater discharge. As a result, the U.S. Geological Survey began a cooperative study of the Borrego Valley with the Borrego Water District (BWD) in 2009. The purpose of the study was to develop a greater understanding of the hydrogeology of the Borrego Valley Groundwater Basin (BVGB) and to provide tools to help evaluate the potential hydrologic effects of future development. The objectives of the study were to (1) improve the understanding of groundwater conditions and land subsidence, (2) incorporate this improved understanding into a model that would assist in the management of the groundwater resources in the Borrego Valley, and (3) use this model to test several management scenarios. This model provides the capability for the BWD and regional stakeholders to quantify the relative benefits of various options for increasing groundwater storage. The study focuses on the period 1945&ndash;2010, with scenarios 50 years into the future. This report documents and presents (1) an analysis of the conceptual model, (2) a description of the hydrologic features, (3) a compilation and analysis of water-quality data, (4) the measurement and analysis of land subsidence by using geophysical and remote sensing techniques, (5) the development and calibration of a two-dimensional borehole-groundwater-flow model to estimate aquifer hydraulic conductivities, (6) the development and calibration of a three-dimensional (3-D) integrated hydrologic flow model, (7) a water-availability analysis with respect to current climate variability and land use, and (8) potential future management scenarios. The integrated hydrologic model, referred to here as the &ldquo;Borrego Valley Hydrologic Model&rdquo; (BVHM), is a tool that can provide results with the accuracy needed for making water-management decisions, although potential future refinements and enhancements could further improve the level of spatial and temporal resolution and model accuracy. Because the model incorporates time-varying inflows and outflows, this tool can be used to evaluate the effects of temporal changes in recharge and pumping and to compare the relative effects of different water-management scenarios on the aquifer system. Overall, the development of the hydrogeologic and hydrologic models, data networks, and hydrologic analysis provides a basis for assessing surface and groundwater availability and potential water-resource management guidelines.

California↗

Hydrogeology of, and ground-water flow in, a valley-fill and carbonate-rock aquifer system near Long Valley in the New Jersey Highlands

The hydrogeology of and ground-water flow in a valley-fill and carbonate-rock aquifer system were evaluated by using numerical-modeling techniques and geochemical interpretations to address concerns about the adequacy of the aquifer system to meet increasing demand for water. The study was conducted during 1987-90 by the U.S. Geological Survey, in cooperation with the New Jersey Department of Environmental Protection and Energy. The effects of recent and anticipated ground-water withdrawals on water levels, stream base flows, and water budgets were estimated. Simulation results indicate that recent withdrawals of 4.7 million gallons per day have resulted in water-level declines of up to 35 feet. Under conditions of increases in withdrawals of 121 percent, water levels would decline up to an additional 28 feet. The magnitude of predicted average base-flow depletion, when compared with historic low flows, indicates that projected increases in withdrawals may substantially deplete seasonal low flow of Drakes Brook and South Branch Raritan River. Results of a water-budget analysis indicate that the sources of water to additional supply wells would include leakage from the overlying valley-fill aquifer and induced leakage of surface water into the aquifer system. Results of water-quality analyses indicate that human activities are affecting the quality of the ground water. With the exception of an elevated iron concentration in water from one well, concentrations of inorganic constituents in water from 75 wells did not exceed New Jersey primary or secondary drinking-water regulations. Volatile organic compounds were detected in water from several wells; in two samples, concentrations of specific compounds exceeded drinking-water regulations.

Water-Resources Investigations Report↗

The effect of discharge and water quality of the Alafia River, Hillsborough River, and the Tampa Bypass Canal on nutrient loading to Hillsborough Bay, Florida

Techniques to measure discharge and nutrient loads in the tidally affected portions of two major rivers tributary to Tampa Bay, the Alafia River and the Hillsborough River, were developed and tested. Discharge, water quality, and total phosphorus and total nitrogen loads for the period April 1, 1991, through March 31, 1992, were evaluated and compared with discharge,water quality, and loads at long-term, nontidal gages in the basins. Long-term discharge and water-quality characteristics at selected sites in the Alafia river and Hillsborough River basins were evaluated. A long-term, decreasing trend in annual-mean discharge was observed for discharges at the Alafia River, Sulphur Springs, and Hillsborough River. Low-flow and high-flow characteristics in the Alafia River and Hillsborough River have changed as well. The decreasing trend in the Alafia River discharges is not due to deficient rainfall but probably is due to decreased ground-water inflow to the river because of long-term declines in the potentiometric surface of the Upper Floridan aquifer. Daily-mean discharges at the mouth of the Alafia River were more variable than discharges at the long-term gage upstream. Daily-mean discharge near the mouth of the river was negative at times, indicating a net loss of water from the river. Daily-mean discharge from the Hillsborough River was minimal from Apil to May 1991, and from late September 1991 to March 1992. During these periods, discharge from Sulphur Springs was a major source of freshwater to the tidally affected reach of the river. Concentrations of total phosphorus and orthophosphorus in the Alafia River above Lithia Springs were the greatest in the 1960's and have generally declined since then. Total nitrogen concentrations have been declining since about 1981. However, increases in nitrate plus nitrite nitrogen concentrations are occurring in Lithia Springs, a second-magnitude spring that flows into the Alafia River. Specific conductance of water discharging from Sulphur Springs to the Hillsborough River has increased from about 124 to more than 2,000 microsiemens per centimeter since 1945. Water quality at the mouth of the Alafia River and Hillsborough River is the result of mixing of freshwater and estuarine water from Hillsborough Bay. Large daily variations in water quality occur at these site because of tidal currents, and vertical stratification of specific conductance is a common feature. Concentrations of phosphorus, nitrate plus nitrite nitrogen, organic carbon, and silica are inversely related to specific conductance at the mouth of the Alafia River. Constituent concentration and discharge data were used to compute loads during the study period. Average daily phosphorus loads were 2.4 tons per day at the mouth of the Alafia River; 0.35 ton per day at the mouth of the Hillsborough River; and 0.06 ton per day at the Tampa Bypass Canal. Average daily nitrogen loads were 1.7 tons per day at the mouth of the Alafia River; 0.86 ton per day at the mouth of the Hillsborough River; and 0.26 ton per day at the Tampa Bypass Canal. The greatest annual loads of phosphorus and nitrogen from the major tributaries to Hillsborough Bay are from the Alafia River, with the greatest loads at the river mouth. Total phosphorus load from the Alafia River was about 894 tons during April 1991 through March 1992, more than six times greater than phosphorus loads form the Hillsborough River. Annual total nitrogen load at the mouth of the Alafia River was about 630 tons, two times greater than at the mouth of the Hillsborough River and more than six times greater than loads from the Tampa Bypass Canal. Basinwide yields of total phosphorus during April 1991 through March 1992 were about 2 tons per square mile at the mouth of the Alafia river basin and were about 0.2 ton per square mile at the mouth of the Hillsborough River. Total nitrogen yield was about 1.5 tons per square mile at the mouth of the Alafia River and about 0.5 ton per s

Water-Resources Investigations Report↗

Use of advanced borehole geophysical techniques to delineate fractured-rock ground-water flow and fractures along water-tunnel facilities in northern Queens County, New York

Advanced borehole geophysical methods were used to assess the geohydrology of crystalline bedrock along the course of a new water tunnel for New York City. The logging methods include natural gamma, spontaneous potential, single-point resistance, mechanical and acoustic caliper, focused electromagnetic induction, electromagnetic resistivity, magnetic susceptibility, borehole-fluid temperature and conductance, differential temperature, heat-pulse flowmeter, acoustic televiewer, borehole deviation, optical televiewer, and borehole radar. Integrated interpretation of the geophysical logs from an 825-foot borehole (1) provided information on the extent, orientation, and structure (foliation and fractures) within the entire borehole, including intensely fractured intervals from which core recovery may be poor; (2) delineated transmissive fracture zones intersected by the borehole and provided estimates of their transmissivity and hydraulic head; and (3) enabled mapping of the location and orientation of structures at distances as much as 100 ft from the borehole. Analyses of the borehole-wall image and the geophysical logs from the borehole on Crescent Street, in northern Queens County, are presented here to illustrate the application of the methods. The borehole penetrates gneiss and other crystalline bedrock that has predominantly southeastward dipping foliation and nearly horizontal and southeastward-dipping fractures. The heat-pulse flowmeter logs obtained under pumping and nonpumping conditions, together with the other geophysical logs, indicate five transmissive fracture zones. More than 90 percent of the open-hole transmissivity is associated with a fracture zone 272 feet BLS (below land surface). A transmissive zone at 787 feet BLS that consists of nearly parallel fractures lies within the projected tunnel path; here the hydraulic head is 12 to 15 feet lower than that of transmissive zones above the 315-foot depth. The 60-megahertz directional borehole radar logs indicate the location and orientation of two closely spaced radar reflectors that would intersect the projection of the borehole below its drilled depth. Subsequent excavation of the tunnel past the borehole allowed comparison of the log analysis with conditions observed in the tunnel. The tunnel was found to intersect gneiss with southeastward dipping foliation; many nearly horizontal fractures; and a southeastward dipping fracture zone whose location, character, and orientation was consistent with that of the mapped radar reflectors. The fracture zone produced inflow to the tunnel at a rate of 50 to 100 gallons per minute. All conditions indicated by the logging methods were consistent with those observed within the tunnel.

New York↗

Equations for estimating selected streamflow statistics in Rhode Island

Regional regression equations were developed for estimating selected natural&mdash;unaffected by alteration&mdash;streamflows of specific flow durations and low-flow frequency statistics for ungaged stream sites in Rhode Island. Selected at-site streamflow statistics are provided for 41 long-term streamgages, 21 short-term streamgages, and 135 partial-record stations in Rhode Island, eastern Connecticut, and southeastern and south-central Massachusetts. The regression equations for estimating selected streamflow statistics and the at-site statistics estimated for each of the 197 sites may be used by Federal, State, and local water managers in addressing water issues in and near Rhode Island. Multiple and simple linear regression equations were developed to estimate the 99-, 98-, 95-, 90-, 85-, 80-, 75-, 70-, 60-, 50-, 40-, 30-, 25-, 20-, 15-, 10-, 5-, 2-, and 1-percent flow durations and the 7Q2 (7-day, 2-year) and 7Q10 (7-day, 10-year) low-flow-frequency statistics. An additional 49 selected statistics, for which regression equations were not developed, also were estimated for the long- and short-term streamgages and partial-record stations for flow durations between the 99.99 and 0.01 percent and for the mean annual, mean monthly, and median monthly streamflows. A total of 70 selected streamflow statistics were estimated for 41 long-term streamgages, 21 short-term streamgages, and 135 partial-record stations in and near Rhode Island. Estimates of the long-term streamflow statistics for the 21 short-term streamgages and 135 partial-record stations were developed by the Maintenance of Variance Extension, type 1 (MOVE.1), record-extension technique. The equations used to estimate selected streamflow statistics were developed by relating the 19 flow-duration and 2 low-flow-frequency statistics to 31 different basin characteristics (physical, land-cover, and climatic) at the 41 long-term and 19 of 21 short-term streamgages (a total of 60 streamgages) in and near Rhode Island. The 135 partial-record stations were not used in the regression analyses. The regression analyses were done by using a user-weighted least-squares technique in the weighted-multiple-linear regression program for the 90- to 1-percent flow-duration statistics. For the 99-, 98-, and 95-percent flow durations and the 7Q2 and 7Q10 statistics, left-censored regression analyses were used to account for zero flows at a few streamgages. The regression analyses determined that two basin characteristics&mdash;drainage area and stream density&mdash;were the only significant explanatory variables for 16 of the 19 flow-duration and the 2 low-flow regression equations. For the 10-, 15-, and 20-percent flow-duration regression equations, drainage area was the only significant explanatory variable. The standard error of the estimate for the 21 regression equations ranged from 17.58 to 141.83 percent. The 99- to 85-percent flow durations and the low-flow statistics 7Q2 and 7Q10 had the highest standard errors of the estimate, ranging from 48.68 to 141.83 percent. The standard error of the estimate for the medium- to high-flow statistics&mdash;the 80- to 1-percent flow durations&mdash;ranged from 17.58 to 37.65 percent, with the standard errors for the 60- to 1-percent flow durations all being less than about 21 percent. Data also are provided to allow the user to calculate the 90-percent prediction intervals for the 21 streamflow statistics. The equations, which are based on data from streams with little to no flow alterations, will provide an estimate of the natural flows for a selected site. They will not estimate flows for altered sites with dams, surface-water withdrawals, groundwater withdrawals (pumping wells), diversions, and wastewater discharges. If the equations are used to estimate streamflow statistics for altered sites, the user should adjust the flow estimates for the alterations. The regression equations should be used only for ungaged sites with drainage areas between 0.52 and 294 square miles and stream densities between 0.94 and 3.49 miles per square mile; these are the ranges of the explanatory variables in the equations.

Rhode Island↗

User's guide to the Variably Saturated Flow (VSF) process to MODFLOW

A new process for simulating three-dimensional (3-D) variably saturated flow (VSF) using Richards' equation has been added to the 3-D modular finite-difference ground-water model MODFLOW. Five new packages are presented here as part of the VSF Process--the Richards' Equation Flow (REF1) Package, the Seepage Face (SPF1) Package, the Surface Ponding (PND1) Package, the Surface Evaporation (SEV1) Package, and the Root Zone Evapotranspiration (RZE1) Package. Additionally, a new Adaptive Time-Stepping (ATS1) Package is presented for use by both the Ground-Water Flow (GWF) Process and VSF. The VSF Process allows simulation of flow in unsaturated media above the ground-water zone and facilitates modeling of ground-water/surface-water interactions. Model performance is evaluated by comparison to an analytical solution for one-dimensional (1-D) constant-head infiltration (Dirichlet boundary condition), field experimental data for a 1-D constant-head infiltration, laboratory experimental data for two-dimensional (2-D) constant-flux infiltration (Neumann boundary condition), laboratory experimental data for 2-D transient drainage through a seepage face, and numerical model results (VS2DT) of a 2-D flow-path simulation using realistic surface boundary conditions. A hypothetical 3-D example case also is presented to demonstrate the new capability using periodic boundary conditions (for example, daily precipitation) and varied surface topography over a larger spatial scale (0.133 square kilometer). The new model capabilities retain the modular structure of the MODFLOW code and preserve MODFLOW's existing capabilities as well as compatibility with commercial pre-/post-processors. The overall success of the VSF Process in simulating mixed boundary conditions and variable soil types demonstrates its utility for future hydrologic investigations. This report presents a new flow package implementing the governing equations for variably saturated ground-water flow, four new boundary condition packages unique to unsaturated flow, the Adaptive Time-Stepping Package for use with both the GWF Process and the new VSF Process, detailed descriptions of the input and output files for each package, and six simulation examples verifying model performance.

Techniques and Methods↗

Land Disturbance Associated with Oil and Gas Development and Effects of Development-Related Land Disturbance on Dissolved-Solids Loads in Streams in the Upper Colorado River Basin, 1991, 2007, and 2025

Oil and gas resource development in the Upper Colorado River Basin (UCRB) has increased substantially since the year 2000. The UCRB encompasses several significant oil and gas producing areas that have the potential for continued oil and gas resource development. Land disturbance associated with oil and gas resource development is caused by activities related to constructing drill pads to contain drilling and well maintenance equipment and roads to access the drill pad. Land disturbed by oil and gas development has the potential to cause increased erosion, stream degradation, habitat fragmentation and alteration, and increase public use of areas that may be environmentally sensitive. Land disturbance resulting from oil and gas resource development has not been monitored and mapped on a regional scale in the UCRB. However, information on the location and age of oil and gas wells in the UCRB is available. These data combined with geographic data analysis and modeling techniques were used to estimate the total area of disturbed land associated with oil and gas resource development in 1991 and in 2007 in the UCRB. Additional information about anticipated oil and gas development in the UCRB was used to project land disturbance to the year 2025. Results of the analysis indicate that approximately 117,500 acres (183 mi2) of total land disturbance was associated with drill pads and related roads in the UCRB in 1991. The estimated area of disturbed land associated with oil and gas development increased 53 percent to 179,400 acres (280 mi2) in 2007. Projecting oil and gas development through 2025 results in a potential near doubling of the land surface disturbance to approximately 319,300 acres (500 mi2). Estimated land disturbance for 1991 and 2007 were input to a contaminant transport model developed for the UCRB to assess the statistical significance of energy-related land disturbance to contributing dissolved solids to basin streams. The statistical assessment was an observational study based on an existing model and available water-quality monitoring data for the basin. No new data were collected for the analysis. The source coefficient calibrated for the disturbed lands associated with oil and gas development in 2007 was zero, which indicated that estimated land disturbance from oil and gas development is not statistically significant in explaining dissolved solids in UCRB streams. The lack of significance in the contaminant transport modeling framework may be due to the amount of available monitoring data, the spatial distribution of monitoring sites with respect to land disturbance, or the overall quantity of land disturbance associated with oil and gas development basin wide. Finally, dissolved-solids loads derived from natural landscapes may be similar to loads derived from lands disturbed by oil and gas resource development. The model recalibration done for this study confirms calibration results from Kenney and others (2009): the most significant contributor to dissolved solids in the UCRB is irrigated agricultural land, which covers an area substantially larger than the estimated area disturbed by oil and gas development and is subjected to artificially applied water.

Scientific Investigations Report↗