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.

964 records · Page 54Linked to original sources

Protocol for Monitoring Fish Assemblages in Pacific Northwest National Parks

Rivers and streams that drain from Olympic, Mount Rainier, and North Cascades National Parks are among the most protected corridors in the lower 48 States, and represent some of the largest tracts of contiguous, undisturbed habitat throughout the range of several key fish species of the Pacific Northwest. These watersheds are of high regional importance as freshwater habitat sanctuaries for native fish, where habitat conditions are characterized as having little to no disturbance from development, channelization, impervious surfaces, roads, diversions, or hydroelectric projects. Fishery resources are of high ecological and cultural importance in Pacific Northwest National Parks, and significantly contribute to economically important recreational, commercial, and tribal fisheries. This protocol describes procedures to monitor trends in fish assemblages, fish abundance, and water temperature in eight rivers and five wadeable streams in Olympic National Park during summer months, and is based on 4 years of field testing. Fish assemblages link freshwater, marine, and terrestrial ecosystems. They also serve as focal resources of national parks and are excellent indicators of ecological conditions of rivers and streams. Despite the vital importance of native anadromous and resident fish populations, there is no existing monitoring program for fish assemblages in the North Coast and Cascades Network. Specific monitoring objectives of this protocol are to determine seasonal and annual trends in: (1) fish species composition, (2) timing of migration of adult fish, (3) relative abundance, (4) age and size structure, (5) extent of non-native and hatchery fish, and (6) water temperature. To detect seasonal and annual trends in fish assemblages in reference sites, we rely on repeated and consistent annual sampling at each monitoring site. The general rationale for the repeated sampling of reference sites is to ensure that we account for the high interannual variability in fish movements and abundances in rivers. One underlying assumption is that the monitoring program is designed in perpetuity, and consequently our capability to detect trends substantially increases with time. The protocol describes sampling designs, methods, training procedures, safety considerations, data management, data analysis, and reporting. The allocation of sampling effort represents a balance between ecological considerations, a sound monitoring approach, and practical limitations caused by logistical constraints and a limited annual budget of $55,000. The widespread declines of native fish species in western North America highlights the importance and urgency of understanding trends in fish assemblages from undisturbed habitats. Seasonal and annual trends in fish assemblages will provide insights at the individual, population, and assemblage level. This protocol will allow managers to detect increases and decreases in abundance of priority management species, and occurrence of non-native, hatchery, and federally listed fish. The detection of trends in fish assemblages will allow for specific management actions that may include: implementation of more appropriate fishing regulations, evaluation of existing hatchery releases, control of non-native fish species, and prioritization of habitat restoration projects. Dissemination and communication of scientific findings on North Coast and Cascades Network fish assemblages will be a core product of this protocol, which will have much relevance to decision makers, park visitors, researchers, and educators.

Techniques and Methods↗

The Farm Process Version 2 (FMP2) for MODFLOW-2005 - Modifications and Upgrades to FMP1

The ability to dynamically simulate the integrated supply-and-demand components of irrigated agricultural is needed to thoroughly understand the interrelation between surface water and groundwater flow in areas where the water-use by vegetation is an important component of the water budget. To meet this need, the computer program Farm Process (FMP1) was updated and refined for use with the U.S. Geological Survey's MODFLOW-2005 groundwater-flow model, and is referred to as MF2005-FMP2. The updated program allows the simulation, analysis, and management of nearly all components of human and natural water use. MF2005-FMP2 represents a complete hydrologic model that fully links the movement and use of groundwater, surface water, and imported water for water consumption of irrigated agriculture, but also of urban use, and of natural vegetation. Supply and demand components of water use are analyzed under demand-driven and supply-constrained conditions. From large- to small-scale settings, the MF2005-FMP2 has the unique set of capabilities to simulate and analyze historical, present, and future conditions. MF2005-FMP2 facilitates the analysis of agricultural water use where little data is available for pumpage, land use, or agricultural information. The features presented in this new version of FMP2 along with the linkages to the Streamflow Routing (SFR), Multi-Node Well (MNW), and Unsaturated Zone Flow (UZF) Packages prevents mass loss to an open system and helps to account for 'all of the water everywhere and all of the time'. The first version, FMP1 for MODFLOW-2000, is limited to (a) transpiration uptake from unsaturated root zones, (b) on-farm efficiency defined solely by farm and not by crop type, (c) a simulation of water use and returnflows related only to irrigated agriculture and not also to non-irrigated vegetation, (d) a definition of consumptive use as potential crop evapotranspiration, (e) percolation being instantly recharged to the uppermost active aquifer, (f) automatic routing of returnflow from runoff either to reaches of tributary stream segments adjacent to a farm or to one reach nearest to the farm's lowest elevation, (g) farm-well pumping from cell locations regardless of whether an irrigation requirement from these cells exists or not, and (h) specified non-routed water transfers from an undefined source outside the model domain. All of these limitations are overcome in MF2005-FMP2. The new features include (a) simulation of transpiration uptake from variably saturated, fully saturated, or ponded root zones (for example, for crops like rice or riparian vegetation), (b) definition of on-farm efficiency not only by farm but also by crop, (c) simulation of water use and returnflow from non-irrigated vegetation (for example, rain-fed agriculture or native vegetation), (d) use of crop coefficients and reference evapotranspiration, (e) simulation of the delay between percolation from farms through the unsaturated zone and recharge into the uppermost active aquifer by linking FMP2 to the UZF Package, (f) an option to manually control the routing of returnflow from farm runoff to streams, (g) an option to limit pumping to wells located only in cells where an irrigation requirement exists, and (h) simulation of water transfers to farms from a series of well fields (for example, recovery well field of an aquifer-storage-and-recovery system, ASR). In addition to the output of an economic budget for each farm between irrigation demand and supply ('Farm Demand and Supply Budget' in FMP1), a new output option called 'Farm Budget' was created for FMP2, which allows the user to track all physical flows into and out of a water accounting unit at all times. Such a unit can represent individual farms, farming districts, natural areas, or urban areas. The example model demonstrates the application of MF2005-FMP2 with delayed recharge through an unsaturated zone, rejected infiltration in a riparian area, changes in de

Techniques and Methods↗

Ground-water, surface-water, and water-chemistry data, Black Mesa Area, northeastern Arizona — 2000–2001, and performance and sensitivity of the 1988 USGS numerical model of the N aquifer

The N aquifer is the major source of water in the 5,400-square-mile area of Black Mesa in northeastern Arizona. Availability of water is an important issue in this area because of continued industrial and municipal use, a growing population, and precipitation of about 6 to 14 inches per year. The monitoring program in Black Mesa has been operating since 1971 and is designed to determine the long-term effects of ground-water withdrawals from the N aquifer for industrial and municipal uses. The monitoring program includes measurements of (1) ground-water pumping, (2) ground-water levels, (3) spring discharge, (4) surface-water discharge, and (5) ground-water chemistry. In 2000, total ground-water withdrawals were 7,740 acre-feet, industrial use was 4,490 acre-feet, and municipal use was 3,250 acre-feet. From 1999 to 2000, total withdrawals increased by 9 percent, industrial use increased by 7 percent, and municipal use increased by 12 percent. From 1999 to 2001, water levels declined in 10 of 15 wells in the unconfined part of the aquifer, and the median change was -0.4 foot. Water levels declined in 8 of 16 wells in the confined part of the aquifer, and the median change was -0.2 foot. From the prestress period (prior to 1965) to 2001, the median water-level change for 33 wells was -17.2 feet. Median water-level changes were -1.2 feet for 15 wells in the unconfined part of the aquifer and -31.0 feet for 18 wells in the confined part. Discharges were measured once in 1999 and once in 2001 at four springs. Discharges decreased by 5 percent and 33 percent at two springs and increased by 3 percent and 81 percent at two springs. For about the past 10 years, discharges did not significantly change in Burro Spring, the unnamed spring near Dennehotso, and Moenkopi School Spring. The record of discharge from a consistent measuring point for Pasture Canyon Spring is too short for statistical analysis of trends. Continuous records of surface-water discharge have been collected from July 1976 to 2000 at Moenkopi Wash, July 1996 to 2000 at Laguna Creek, June 1993 to 2000 at Dinnebito Wash, and April 1994 to 2000 at Polacca Wash. Median flows for November, December, January, and February of each water year were used as an index of ground-water discharge to those streams. There is no significant trend in the median winter flows for Moenkopi Wash from 1977 to 2000. The records for the other three streams are too short for a statistical analysis of trends. The median winter flows for Dinnebito Wash and Polacca Wash, however, appear to have decreased during the last 6 years. There is no apparent trend in the median winter flows for Laguna Creek since 1997. In 2001, water samples were collected from 12 wells and 4 springs and analyzed for selected chemical constituents. Dissolved-solids concentrations ranged from 102 to 628 milligrams per liter. Water samples from 9 of the wells and from the 4 springs had less than 350 milligrams per liter of dissolved solids. Water-chemistry data with sufficient years of record for a statistical analysis of trends over time are available from 7 wells and 4 springs. From about the mid-1980s or early 1990s to 2001 there are no significant trends in the concentrations of dissolved solids, chloride, and sulfate in water samples from 6 of the 7 wells. The concentration of one tested constituent (dissolved solids) in samples from Rocky Ridge PM3 significantly increased from 1990 to 2001. From the late 1980s to 2001, there are no significant trends in the concentrations of dissolved solids, chloride, and sulfate in water samples from Burro Spring, the unnamed spring near Dennehotso, and Pasture Canyon Spring. From 1987 to 2001, concentrations of chloride and sulfate significantly increased in water samples from Moenkopi School Spring and concentrations of dissolved solids did not significantly change. The performance and sensitivity of the 1988 USGS numerical model of the N aquifer were analyzed. The overall performance of the model in steady-state conditions is reasonable for residuals of heads (difference between observed and simulated steady-state heads); 80 percent of the absolute values of residuals are less than 38 feet. Simulated flows are about 40 percent different than estimated flows at two of three discharge areas; however, this comparison is only a rough approximation of performance because the accuracy of the estimated steady-state flows is uncertain. The overall performance of the model for transient conditions is fair for residuals of changes in head (difference between observed and simulated changes in head from steady state to 1999); 80 percent of the absolute values of residuals are less than 31 feet. The model is biased in two areas. In the Tuba City area, simulated changes in head are more negative than observed changes in head; all six residuals are positive, and three residuals are between 75 and 155 feet. In the confined area of the aquifer, observed changes in head are more negative than simulated changes in head; 12 of the 17 residuals are negative, and 8 residuals are between -57 and -20 feet. Analysis of model sensitivity indicates that recharge, transmissivity, and storage coefficient are the most important parameters for estimating heads, changes in heads, and flows. A strong correlation between recharge and transmissivity and a lack of independent and reliable estimates of recharge, transmissivity, and discharge create a uniqueness problem in model calibration. Several models could be constructed and calibrated with different values of recharge or transmissivity and still have similar fits to the observed data. Information from recent data and studies and more advanced modeling techniques could be used to develop a more representative and less uncertain model. Future data collection and studies should focus on obtaining a better definition of recharge, discharge, transmissivity, and storage coefficient.

Arizona↗

A comparison of runoff quantity and quality from two small basins undergoing implementation of conventional- and low-impact-development (LID) strategies: Cross Plains, Wisconsin, water years 1999–2005

Environmental managers are often faced with the task of designing strategies to accommodate development while minimizing adverse environmental impacts. Low-impact development (LID) is one such strategy that attempts to mitigate environmental degradation commonly associated with impervious surfaces. The U.S. Geological Survey, in cooperation with the Wisconsin Department of Natural Resources, studied two residential basins in Cross Plains, Wis., during water years 1999–2005. A paired-basin study design was used to compare runoff quantity and quality from the two basins, one of which was developed in a conventional way and the other was developed with LID. The conventional-developed basin (herein called “conventional basin”) consisted of curb and gutter, 40-foot street widths, and a fully connected stormwater-conveyance system. The LID basin consisted of grassed swales, reduced impervious area (32-foot street widths), street inlets draining to grass swales, a detention pond, and an infiltration basin. Data collected in the LID basin represented predevelopment through near-complete build-out conditions. Smaller, more frequent precipitation events that produced stormwater discharge from the conventional basin were retained in the LID basin. Only six events with precipitation depths less than or equal to 0.4 inch produced measurable discharge from the LID basin. Of these six events, five occurred during winter months when underlying soils are commonly frozen, and one was likely a result of saturated soil from a preceding storm. In the conventional basin, the number of discharge events, using the same threshold of precipitation depth, was 180, with nearly one-half of those resulting from precipitation depths less than 0.2 inch. Precipitation events capable of producing appreciable discharge in the LID basin were typically those of high intensity or precipitation depth or those that occurred after soils were already saturated. Total annual discharge volume measured from the conventional basin ranged from 1.3 to 9.2 times that from the LID basin. Development of the LID basin did not appreciably alter the hydrologic response to precipitation characterized during predevelopment conditions. Ninety-five percent or more of precipitation in the LID basin was retained during each year of construction from predevelopment through near-complete build-out, surpassing the 90-percent benchmark established for new development by the Wisconsin Department of Natural Resources. The amount of precipitation retained in the conventional basin did not exceed 94 percent and fell below the 90-percent standard 2 of the 6 years monitored. Much of the runoff in the LID basin was retained by an infiltration basin, the largest control structure used to mitigate storm-runoff quantity and quality. The infiltration basin also was the last best-management practice (BMP) used to treat runoff before it left the LID basin as discharge. From May 25, 2002, to September 30, 2005, only 24 of 155 precipitation events exceeded the retention/ infiltrative capacity of the infiltration basin. The overall reduction in runoff volume from these few events was 51 percent. The effectiveness of the infiltration basin decreased as precipitation intensities exceeded 0.5 inch per hour. Annual loads were estimated to characterize the overall effectiveness of low-impact design practices for mitigating delivery of total solids, total suspended solids, and total phosphorus. Annual loads of these three constituents were greater in the LID basin than in the conventional basin in 2000 and 2004. Seventy percent or more of all constituent annual loads were associated with two discharge events in 2000, and a single discharge event produced 50 percent or more of constituent annual loads in 2004. Each of these discharge events was associated with considerable precipitation depths and (or) intensities, ranging from 4.89 to 6.21 inches and from 1.13 to 1.2 inches per hour, respectively. These same storms did not contribute as much of the annual load in the conventional basin. With large storms and saturated soils, the ability of low-impact design techniques to reduce runoff, and thus constituent loads, can be greatly diminished. For both the LID and conventional basins, the temperature of runoff was largely affected by ambient air temperatures. However, the temperature of discharge from the LID basin increased upon runoff cessation. This increase is likely due to solar heating of water that is temporarily stored in the detention pond and infiltration basin.

Wisconsin↗

Results of the application of seismic-reflection and electromagnetic techniques for near-surface hydrogeologic and environmental investigations at Fort Bragg, North Carolina

As part of the U.S. Geological Survey's Resource Conservation and Recovery Act, Facilities Investigations at Fort Bragg, North Carolina, selected geophysical techniques were evaluated for their usefulness as assessment tools for determining subsurface geology, delineating the areal extent of potentially contaminated landfill sites, and locating buried objects and debris of potential environmental concern. Two shallow seismic-reflection techniques (compression and shear wave) and two electromagnetic techniques (ground-penetrating radar and terrain conductivity) were evaluated at several sites at the U.S. Army Base. The electromagnetic techniques also were tested for tolerance to cultural noise, such as nearby fences, vehicles, and power lines. For the terrain conductivity tests, two instruments were used--the EM31 and EM34, which have variable depths of exploration. The shallowest reflection event was 70 feet below land surface observed in common-depth point, stacked compression-wave data from 24- and 12-fold shallow-seismic-reflection surveys. Several reflection events consistent with clay-sand interfaces between 70 and 120 feet below land surface, along with basement-saprolite surfaces, were imaged in the 24-fold, common- depth-point stacked data. 12-fold, common-depth-point stacked data set contained considerably more noise than the 24-fold, common-depth-point data, due to reduced shot-to-receiver redundancy. Coherent stacked reflection events were not observed in the 24-fold, common-depth-point stacked shear-wave data because of the partial decoupling of the shear- wave generator from the ground. At one site, ground-penetrating radar effectively delineated a shallow, 2- to 5-foot thick sand unit bounded by thin (less than 1 foot) clay layers. The radar signal was completely attenuated where the overlying and underlying clay units thickened and the sand unit thinned. The pene- tration depth of the radar signal was less than 10 feet below land surface. A slight increase in electromagnetic conductivity across shallow sampling EM31 and EM34 profiles provided corroborative evidence of the shallow, thickening clay units. Plots of raw EM31 and EM34 data provided no direct interpretable information to delineate sand and clay units in the shallow subsurface. At two sites, the ground-penetrating radar effectively delineated the lateral continuity of surficial sand units 5 to 25 feet in thickness and the tops of their underlying clay units. The effective exploration depth of the ground-penetrating radar was limited by the proximity of clay units to the subsurface and their thickness. The ground-penetrating radar delineated the areal extent and depth of cover at a previously unrecognized extension of a trench-like landfill underlying a vehicle salvage yard. Attenuation of the radar signal beneath the landfill cover and the adjacent subsurface clays made these two mediums indistinguishable by ground-penetrating radar; however, EM31 data indicated that the electrical conductivity of the landfill was higher than the subsurface material adjacent to the landfill. The EM31 and EM34 conductivity surveys defined the areal extent of a landfill whose boundaries were inaccurately mapped, and also identified the locations of an old dumpsite and waste incinerator site at another landfill. A follow-up ground-penetrating radar survey of the abandoned dumpsite showed incongruities in some of the shallow radar reflections interpreted as buried refuse dispersed throughout the landfill. The ground-penetrating radar and EM31 effectively delineated a shallow buried fuel-oil tank. Of the three electromagnetic instruments, the ground-penetrating radar with the shielded 100-megahertz antenna was the least affected by cultural noise followed, in order, by the EM31 and EM34. The combination of terrain- conductivity and ground-penetrating radar for the site assessment of the landfill provided a powerful means to identify the areal extent of the landfill, potenti

Water-Resources Investigations Report↗

Predicting sea-level rise vulnerability of terrestrial habitat and wildlife of the Northwestern Hawaiian Islands

If current climate change trends continue, rising sea levels may inundate low-lying islands across the globe, placing island biodiversity at risk. Recent models predict a rise of approximately one meter (1 m) in global sea level by 2100, with larger increases possible in areas of the Pacific Ocean. Pacific Islands are unique ecosystems home to many endangered endemic plant and animal species. The Northwestern Hawaiian Islands (NWHI), which extend 1,930 kilometers (km) beyond the main Hawaiian Islands, are a World Heritage Site and part of the Papahanaumokuakea Marine National Monument. These NWHI support the largest tropical seabird rookery in the world, providing breeding habitat for 21 species of seabirds, 4 endemic land bird species and essential foraging, breeding, or haul-out habitat for other resident and migratory wildlife. In recent years, concern has grown about the increasing vulnerability of the NWHI and their wildlife populations to changing climatic patterns, particularly the uncertainty associated with potential impacts from global sea-level rise (SLR) and storms. In response to the need by managers to adapt future resource protection strategies to climate change variability and dynamic island ecosystems, we have synthesized and down scaled analyses for this important region. This report describes a 2-year study of a remote northwestern Pacific atoll ecosystem and identifies wildlife and habitat vulnerable to rising sea levels and changing climate conditions. A lack of high-resolution topographic data for low-lying islands of the NWHI had previously precluded an extensive quantitative model of the potential impacts of SLR on wildlife habitat. The first chapter (chapter 1) describes the vegetation and topography of 20 islands of Papahanaumokuakea Marine National Monument, the distribution and status of wildlife populations, and the predicted impacts for a range of SLR scenarios. Furthermore, this chapter explores the potential effects of SLR on wildlife breeding habitats for each island. The subsequent chapter (chapter 2) details a study of the Laysan Island ecosystem, describing a quantitative model that incorporates SLR, storm wave, and rising groundwater inundation. Wildlife, storm, and oceanographic data allowed for an assessment of the phenological and spatial vulnerability of Laysan Island's breeding bird species to SLR and storms. Using remote sensing and geospatial techniques, we estimated topography, classified vegetation, modeled SLR, and evaluated a range of climate change scenarios. On the basis of high-resolution airborne data collected during 2010-11 (root-mean-squared error = 0.05-0.18 m), we estimated the maximum elevation of 20 individual islands extending from Kure Atoll to French Frigate Shoals (range: 1.8-39.7 m) and computed the mean elevation (1.7 m, standard deviation 1.1 m) across all low-lying islands. We also analyzed general climate models to describe rainfall and temperature scenarios expected to influence adaptation of some plants and animals for this region. Outcomes for the NWHI predicted an increase in temperature of 1.8-2.6 degrees Celsius (°C) and an annual decrease in precipitation of 24.7-76.3 millimeters (mm) across the NWHI by 2100. Our models of passive SLR (excluding wave-driven effects, erosion, and accretion) showed that approximately 4 percent of the total land area in the NWHI will be lost with scenarios of +1.0 m of SLR and 26 percent will be lost with +2.0 m of SLR. Some atolls are especially vulnerable to SLR. For example, at Pearl and Hermes Atoll our analysis indicated substantial habitat losses with 43 percent of the land area inundated at +1.0 m SLR and 92 percent inundated at +2.0 m SLR. Across the NWHI, seven islands will be completely submerged with +2.0 m SLR. The limited global ranges of some tropical nesting birds make them particularly vulnerable to climate change impacts in the NWHI. Climate change scenarios and potential SLR impacts presented here emphasize the need for early climate change adaptation and mitigation planning, especially for species with limited breeding distributions and/or ranges restricted primarily to the low-lying NWHI: Cyperus pennatiformis var. bryanii , Black-footed Albatross ( Phoebastria nigripes ), Laysan Albatross ( P. immutabilis ), Bonin Petrel ( Pterodroma hypoleuca ), Gray-backed Tern ( Onychoprion lunatus ), Laysan Teal ( Anas laysanensis ), Laysan Finch ( Telespiza cantans ), and Hawaiian monk seal ( Monachus schauinslandi ). Furthermore, SLR scenarios that include the effects of wave dynamics and groundwater rise may indicate amplified vulnerability to climate change driven habitat loss on low-lying islands. In chapter 2, we incorporated the combined effects of SLR, dynamic wave-driven inundation, and rising groundwater in a quantitative study specifically for the Laysan Island ecosystem. This is the first hydrodynamic model to simulate the combined impacts of SLR and wave-driven inundation in the NWHI. We developed a high-resolution digital elevation model (mean vertical accuracy of 0.32 m) for the island. Then using recent satellite imagery, geospatial models, and historical oceanographic, storm, and biological data we estimated potential inundation extent, habitat loss, and wildlife population impacts for a range of potential SLR scenarios (0.00, +0.50, +1.00, +1.50, and +2.00 m) that may occur over the next century. Additionally, we estimated the carrying capacity of Laysan Island for five species based on the available population monitoring data and described how potential losses in nesting habitat could influence population dynamics for Black-footed Albatross, Laysan Albatross, Red-footed Booby (Sula sula), Laysan Teal, and Laysan Finch. For some other seabird populations (Masked Booby, S. dactylatra ; Brown Booby, S. leucogaster ; Great Frigatebird, Fregata minor ; and Sooty Tern, Onychoprion fuscata ), we used recent colony distribution data, land cover maps, and nesting behavior to estimate potential losses of nesting habitat from SLR and wave-driven inundation. We observed far greater potential impacts of SLR to wildlife with the dynamic wave-driven modeling approach than with the passive modeling approach. Depending on SLR scenario and coastal orientation, during storms under a +2.00 m SLR scenario, the wave-driven inundation model predicted three times more inundation than the passive model (17.2 percent of total terrestrial area versus 4.6 percent, respectively). Large-wave events generally added 1 m of water height to passive inundation surfaces, therefore our dynamic models (during storm events) forecasted comparable inundation extents earlier than passive models. Although wave-driven water levels were highest in the northwest quadrant of Laysan Island, the greatest extent of inundation occurred in the southeast where coastal dunes less than 3 m above mean sea level provide little protection from wave-driven inundation. When wave-driven inundation was included in the SLR model for Laysan Island greater nesting habitat loss and potential impacts on wildlife population dynamics were predicted. The consequences of habitat loss due to SLR may be worse for species with colonies in the wave-exposed coastal zones (for example, Black-footed Albatross) and for populations already near the island's carrying capacity (for example, Laysan Teal). Species whose peak incubation and chick-rearing periods coincide with seasonally high wave heights also will be increasingly vulnerable, especially those species nesting on the ground in areas vulnerable to inundation, such as Gray-backed Tern and Black-footed Albatross. Other species that have space for population growth, or are not restricted to a narrow range of habitat types on Laysan (for instance, Sooty Terns), may be less sensitive to habitat loss from SLR over the next century. Our assessments of inundation risk, habitat loss, and wildlife species vulnerability synthesize current knowledge about individual islands and contribute to a broader understanding of the impacts of inundation from SLR and storm-induced waves. Yet, most NWHI and their bird populations lack monitoring data to evaluate adaptations to and impacts of climate change. Exceptions include some data sets from long-term monitoring of wildlife populations, tides, or weather at French Frigate Shoals, Laysan Island, and Midway Atoll. These data sets are potentially valuable baselines, which could be informative for adaptive learning (integrating management and science) to predict, adapt, and mitigate the effects of climate change on NWHI wildlife in the future. This study provides the first quantitative vulnerability assessment for all of the low-lying NWHI, and results identify biological communities, locales, and resident endangered species of Papahanaumokuakea Marine National Monument expected to be at risk from SLR. This report is also intended as a reference for managers and conservation planners, a tool to identify and potentially reduce uncertainty, and a starting place for developing climate change monitoring priorities and future scientific studies.

Hawai'i↗

Mineral-deposit model for lithium-cesium-tantalum pegmatites

Lithium-cesium-tantalum (LCT) pegmatites comprise a compositionally defined subset of granitic pegmatites. The major minerals are quartz, potassium feldspar, albite, and muscovite; typical accessory minerals include biotite, garnet, tourmaline, and apatite. The principal lithium ore minerals are spodumene, petalite, and lepidolite; cesium mostly comes from pollucite; and tantalum mostly comes from columbite-tantalite. Tin ore as cassiterite and beryllium ore as beryl also occur in LCT pegmatites, as do a number of gemstones and high-value museum specimens of rare minerals. Individual crystals in LCT pegmatites can be enormous: the largest spodumene was 14 meters long, the largest beryl was 18 meters long, and the largest potassium feldspar was 49 meters long. Lithium-cesium-tantalum pegmatites account for about one-fourth of the world’s lithium production, most of the tantalum production, and all of the cesium production. Giant deposits include Tanco in Canada, Greenbushes in Australia, and Bikita in Zimbabwe. The largest lithium pegmatite in the United States, at King’s Mountain, North Carolina, is no longer being mined although large reserves of lithium remain. Depending on size and attitude of the pegmatite, a variety of mining techniques are used, including artisanal surface mining, open-pit surface mining, small underground workings, and large underground operations using room-and-pillar design. In favorable circumstances, what would otherwise be gangue minerals (quartz, potassium feldspar, albite, and muscovite) can be mined along with lithium and (or) tantalum as coproducts. Most LCT pegmatites are hosted in metamorphosed supracrustal rocks in the upper greenschist to lower amphibolite facies. Lithium-cesium-tantalum pegmatite intrusions generally are emplaced late during orogeny, with emplacement being controlled by pre-existing structures. Typically, they crop out near evolved, peraluminous granites and leucogranites from which they are inferred to be derived by fractional crystallization. In cases where a parental granite pluton is not exposed, one is inferred to lie at depth. Lithium-cesium-tantalum LCT pegmatite melts are enriched in fluxing components including H2O, F, P, and B, which depress the solidus temperature, lower the density, and increase rates of ionic diffusion. This, in turn, enables pegmatites to form thin dikes and massive crystals despite having a felsic composition and temperatures that are significantly lower than ordinary granitic melts. Lithium-cesium-tantalum pegmatites crystallized at remarkably low temperatures (about 350–550 °C) in a remarkably short time (days to years). Lithium-cesium-tantalum pegmatites form in orogenic hinterlands as products of plate convergence. Most formed during collisional orogeny (for example, Kings Mountain district, North Carolina). Specific causes of LCT pegmatite-related magmatism could include: ordinary arc processes; over thickening of continental crust during collision or subduction; slab breakoff during or after collision; slab delamination before, during, or after collision; and late collisional extensional collapse and consequent decompression melting. Lithium-cesium-tantalum pegmatite deposits are present in all continents including Antarctica and in rocks spanning 3 billion years of Earth history. The global age distribution of LCT pegmatites is similar to those of common pegmatites, orogenic granites, and detrital zircons. Peak times of LCT pegmatite genesis at about 2640, 1800, 960, 485, and 310 Ma (million years before present) correspond to times of collisional orogeny and supercontinent assembly. Between these pulses were long intervals when few or no LCT pegmatites formed. These minima overlap with supercontinent tenures at ca. 2450–2225, 1625–1000, 875–725, and 250–200 Ma. Exploration and assessment for LCT pegmatites are guided by a number of observations. In frontier areas where exploration has been minimal at best, the key first-order criteria are an orogenic hinterland setting, appropriate regional metamorphic grades, and the presence of evolved granites and common granitic pegmatites. New LCT pegmatites are most likely to be found near known deposits. Pegmatites tend to show a regional mineralogical and geochemical zoning pattern with respect to the inferred parental granite, with the greatest enrichment in the more distal pegmatites. Mineral-chemical trends in common pegmatites that can point toward an evolved LCT pegmatite include: increasing rubidium in potassium feldspar, increasing lithium in white mica, increasing manganese in garnet, and increasing tantalum and manganese in columbite-tantalite. Most LCT pegmatite bodies show a distinctive internal zonation featuring four zones: border, wall, intermediate (where lithium, cesium, and tantalum are generally concentrated), and core. This zonation is expressed both in cross section and map view; thus, what may appear to be a common pegmatite may instead be the edge of a mineralized body. Neither lithium-cesium-tantalum pegmatites nor their parental granites are likely to cause serious environmental concerns. Soils and country rock surrounding a LCT pegmatite, as well as waste from mining operations, may be enriched in characteristic elements relative to global average soil and bedrock values. These elements may include lithium, cesium, tantalum, beryllium, boron, fluorine, phosphorus, manganese, gallium, rubidium, niobium, tin, and hafnium. Among this suite of elements, however, the only ones that might present a concern for environmental health are beryllium and fluorine, which are included in the U.S. Environmental Protection Agency drinking-water regulations with maximum contaminant levels of 4 micrograms per liter and 4 milligrams per liter, respectively.

Scientific Investigations Report↗

Description of input and examples for PHREEQC version 3: A computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations

PHREEQC version 3 is a computer program written in the C and C++ programming languages that is designed to perform a wide variety of aqueous geochemical calculations. PHREEQC implements several types of aqueous models: two ion-association aqueous models (the Lawrence Livermore National Laboratory model and WATEQ4F), a Pitzer specific-ion-interaction aqueous model, and the SIT (Specific ion Interaction Theory) aqueous model. Using any of these aqueous models, PHREEQC has capabilities for (1) speciation and saturation-index calculations; (2) batch-reaction and one-dimensional (1D) transport calculations with reversible and irreversible reactions, which include aqueous, mineral, gas, solid-solution, surface-complexation, and ion-exchange equilibria, and specified mole transfers of reactants, kinetically controlled reactions, mixing of solutions, and pressure and temperature changes; and (3) inverse modeling, which finds sets of mineral and gas mole transfers that account for differences in composition between waters within specified compositional uncertainty limits. Many new modeling features were added to PHREEQC version 3 relative to version 2. The Pitzer aqueous model ( pitzer.dat database, with keyword PITZER ) can be used for high-salinity waters that are beyond the range of application for the Debye-Hückel theory. The Peng-Robinson equation of state has been implemented for calculating the solubility of gases at high pressure. Specific volumes of aqueous species are calculated as a function of the dielectric properties of water and the ionic strength of the solution, which allows calculation of pressure effects on chemical reactions and the density of a solution. The specific conductance and the density of a solution are calculated and printed in the output file. In addition to Runge-Kutta integration, a stiff ordinary differential equation solver (CVODE) has been included for kinetic calculations with multiple rates that occur at widely different time scales. Surface complexation can be calculated with the CD-MUSIC (Charge Distribution MUltiSIte Complexation) triple-layer model in addition to the diffuse-layer model. The composition of the electrical double layer of a surface can be estimated by using the Donnan approach, which is more robust and faster than the alternative Borkovec-Westall integration. Multicomponent diffusion, diffusion in the electrostatic double layer on a surface, and transport of colloids with simultaneous surface complexation have been added to the transport module. A series of keyword data blocks has been added for isotope calculations— ISOTOPES, CALCULATE_VALUES, ISOTOPE_ALPHAS, ISOTOPE_RATIOS, and NAMED_EXPRESSIONS . Solution isotopic data can be input in conventional units (for example, permil, percent modern carbon, or tritium units) and the numbers are converted to moles of isotope by PHREEQC. The isotopes are treated as individual components (they must be defined as individual master species) so that each isotope has its own set of aqueous species, gases, and solids. The isotope-related keywords allow calculating equilibrium fractionation of isotopes among the species and phases of a system. The calculated isotopic compositions are printed in easily readable conventional units. New keywords and options facilitate the setup of input files and the interpretation of the results. Keyword data blocks can be copied (keyword COPY ) and deleted (keyword DELETE ). Keyword data items can be altered by using the keyword data blocks with the _MODIFY extension and a simulation can be run with all reactants of a given index number (keyword RUN_CELLS ). The definition of the complete chemical state of all reactants of PHREEQC can be saved in a file in a raw data format ( DUMP and _RAW keywords). The file can be read as part of another input file with the INCLUDE$ keyword. These keywords facilitate the use of IPhreeqc, which is a module implementing all PHREEQC version 3 capabilities; the module is designed to be used in other programs that need to implement geochemical calculations; for example, transport codes. Charting capabilities have been added to some versions of PHREEQC. Charting capabilities have been added to Windows distributions of PHREEQC version 3. (Charting on Linux requires installation of Wine.) The keyword data block USER_GRAPH allows selection of data for plotting and manipulation of chart appearance. Almost any results from geochemical simulations (for example, concentrations, activities, or saturation indices) can be retrieved by using Basic language functions and specified as data for plotting in USER_GRAPH . Results of transport simulations can be plotted against distance or time. Data can be added to a chart from tab-separated-values files. All input for PHREEQC version 3 is defined in keyword data blocks, each of which may have a series of identifiers for specific types of data. This report provides a complete description of each keyword data block and its associated identifiers. Input files for 22 examples that demonstrate most of the capabilities of PHREEQC version 3 are described and the results of the example simulations are presented and discussed.

Techniques and Methods↗

Estimating peak-flow frequency statistics for selected gaged and ungaged sites in naturally flowing streams and rivers in Idaho

The U.S. Geological Survey, in cooperation with the Idaho Transportation Department, updated regional regression equations to estimate peak-flow statistics at ungaged sites on Idaho streams using recent streamflow (flow) data and new statistical techniques. Peak-flow statistics with 80-, 67-, 50-, 43-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities (1.25-, 1.50-, 2.00-, 2.33-, 5.00-, 10.0-, 25.0-, 50.0-, 100-, 200-, and 500-year recurrence intervals, respectively) were estimated for 192 streamgages in Idaho and bordering States with at least 10 years of annual peak-flow record through water year 2013. The streamgages were selected from drainage basins with little or no flow diversion or regulation. The peak-flow statistics were estimated by fitting a log-Pearson type III distribution to records of annual peak flows and applying two additional statistical methods: (1) the Expected Moments Algorithm to help describe uncertainty in annual peak flows and to better represent missing and historical record; and (2) the generalized Multiple Grubbs Beck Test to screen out potentially influential low outliers and to better fit the upper end of the peak-flow distribution. Additionally, a new regional skew was estimated for the Pacific Northwest and used to weight at-station skew at most streamgages. The streamgages were grouped into six regions (numbered 1_2, 3, 4, 5, 6_8, and 7, to maintain consistency in region numbering with a previous study), and the estimated peak-flow statistics were related to basin and climatic characteristics to develop regional regression equations using a generalized least squares procedure. Four out of 24 evaluated basin and climatic characteristics were selected for use in the final regional peak-flow regression equations. Overall, the standard error of prediction for the regional peak-flow regression equations ranged from 22 to 132 percent. Among all regions, regression model fit was best for region 4 in west-central Idaho (average standard error of prediction=46.4 percent; pseudo-R 2 >92 percent) and region 5 in central Idaho (average standard error of prediction=30.3 percent; pseudo-R 2 >95 percent). Regression model fit was poor for region 7 in southern Idaho (average standard error of prediction=103 percent; pseudo-R 2 <78 percent) compared to other regions because few streamgages in region 7 met the criteria for inclusion in the study, and the region’s semi-arid climate and associated variability in precipitation patterns causes substantial variability in peak flows. A drainage area ratio-adjustment method, using ratio exponents estimated using generalized least-squares regression, was presented as an alternative to the regional regression equations if peak-flow estimates are desired at an ungaged site that is close to a streamgage selected for inclusion in this study. The alternative drainage area ratio-adjustment method is appropriate for use when the drainage area ratio between the ungaged and gaged sites is between 0.5 and 1.5. The updated regional peak-flow regression equations had lower total error (standard error of prediction) than all regression equations presented in a 1982 study and in four of six regions presented in 2002 and 2003 studies in Idaho. A more extensive streamgage screening process used in the current study resulted in fewer streamgages used in the current study than in the 1982, 2002, and 2003 studies. Fewer streamgages used and the selection of different explanatory variables were likely causes of increased error in some regions compared to previous studies, but overall, regional peak‑flow regression model fit was generally improved for Idaho. The revised statistical procedures and increased streamgage screening applied in the current study most likely resulted in a more accurate representation of natural peak-flow conditions. The updated, regional peak-flow regression equations will be integrated in the U.S. Geological Survey StreamStats program to allow users to estimate basin and climatic characteristics and peak-flow statistics at ungaged locations of interest. StreamStats estimates peak-flow statistics with quantifiable certainty only when used at sites with basin and climatic characteristics within the range of input variables used to develop the regional regression equations. Both the regional regression equations and StreamStats should be used to estimate peak-flow statistics only in naturally flowing, relatively unregulated streams without substantial local influences to flow, such as large seeps, springs, or other groundwater-surface water interactions that are not widespread or characteristic of the respective region.

Idaho↗

Borehole geophysical investigation of a formerly used defense site, Machiasport, Maine, 2003-2006

The U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, collected borehole geophysical logs in 18 boreholes and interpreted the data along with logs from 19 additional boreholes as part of an ongoing, collaborative investigation at three environmental restoration sites in Machiasport, Maine. These sites, located on hilltops overlooking the seacoast, formerly were used for military defense. At each of the sites, chlorinated solvents, used as part of defense-site operations, have contaminated the fractured-rock aquifer. Borehole geophysical techniques and hydraulic methods were used to characterize bedrock lithology, fractures, and hydraulic properties. In addition, each geophysical method was evaluated for effectiveness for site characterization and for potential application for further aquifer characterization and (or) evaluation of remediation efforts. Results of borehole geophysical logging indicate the subsurface is highly fractured, metavolcanic, intrusive, metasedimentary bedrock. Selected geophysical logs were cross-plotted to assess correlations between rock properties. These plots included combinations of gamma, acoustic reflectivity, electromagnetic induction conductivity, normal resistivity, and single-point resistance. The combined use of acoustic televiewer (ATV) imaging and natural gamma logs proved to be effective for delineating rock types. Each of the rock units in the study area could be mapped in the boreholes, on the basis of the gamma and ATV reflectivity signatures. The gamma and mean ATV reflectivity data were used along with the other geophysical logs for an integrated interpretation, yielding a determination of quartz monzonite, rhyolite, metasedimentary units, or diabase/gabbro rock types. The interpretation of rock types on the basis of the geophysical logs compared well to drilling logs and geologic mapping. These results may be helpful for refining the geologic framework at depth. A stereoplot of all fractures intersecting the boreholes indicates numerous fractures, a high proportion of steeply dipping fractures, and considerable variation in fracture orientation. Low-dip-angle fractures associated with unloading and exfoliation are also present, especially at a depth of less than 100 feet below the top of casing. These sub-horizontal fractures help to connect the steeply dipping fractures, making this a highly connected fracture network. The high variability in the fracture orientations also increases the connectivity of the fracture network. A preliminary comparison of all fracture data from all the boreholes suggests fracturing decreases with depth. Because all the boreholes were not drilled to the same depth, however, there is a clear sampling bias. Hence, the deepest boreholes are analyzed separately for fracture density. For the deepest boreholes in the study, the intensity of fracturing does not decline significantly with depth. It is possible the fractures observed in these boreholes become progressively tighter or closed with depth, but this is difficult to verify with the borehole methods used in this investigation. The fact that there are more sealed fractures at depth (observed in optical televiewer logs in some of the boreholes) may indicate less opening of the sealed fractures, less water moving through the rock, and less weathering of the fracture infilling minerals. Although the fracture orientation remained fairly constant with depth, differences in the fracture patterns for the three restoration sites indicate the orientation of fractures varies across the study area. The fractures in boreholes on Miller Mountain predominantly strike northwest-southeast, and to a lesser degree they strike northeast. The fractures on or near the summit of Howard Mountain strike predominantly east-west and dip north and south, and the fractures near the Transmitter Site strike northeast-southwest and dip northwest and southeast. The fracture populations for the boreholes on or near the summit of Howard Mountain show more variation than at the other two sites. This variation may be related to the proximity of the fault, which is northeast of the summit of Howard Mountain. In a side-by-side comparison of stereoplots from selected boreholes, there was no clear correspondence between fracture orientation and proximity to the fault. There is, however, a difference in the total populations of fractures for the boreholes on or near the summit of Howard Mountain and the boreholes near the Transmitter Site. Further to the southwest and further away from the fault, the fractures at the Transmitter Site predominantly strike northeast-southwest and northwest-southeast.Heat-pulse flowmeter (HPFM) logging was used to identify transmissive fractures and to estimate the hydraulic properties along the boreholes. Ambient downflow was measured in 13 boreholes and ambient upflow was measured in 9 boreholes. In nine other bedrock boreholes, the HPFM did not detect measurable vertical flow. The observed direction of vertical flow in the boreholes generally was consistent with the conceptual flow model of downward movement in recharge locations and upward flow in discharge locations or at breaks in the slope of land surface. Under low-rate pumping or injection rates [0.25 to 1 gallon per minute (gal/min)], one to three inflow zones were identified in each borehole. Two limitations of HPFM methods are (1) the HPFM can only identify zones within 1.5 to 2 orders of magnitude of the most transmissive zone in each borehole, and (2) the HPFM cannot detect flow rates less than 0.010 + or - 0.005 gal/min, which corresponds to a transmissivity of about 1 foot squared per day (ft2/d). Consequently, the HPFM is considered an effective tool for identifying the most transmissive fractures in a borehole, down to its detection level. Transmissivities below that cut-off must be measured with another method, such as packer testing or fluid-replacement logging. Where sufficient water-level and flowmeter data were available, HPFM results were numerically modeled. For each borehole model, the fracture location and measured flow rates were specified, and the head and transmissivity of each fracture zone were adjusted until a model fit was achieved with the interpreted ambient and stressed flow profiles. The transmissivities calculated by this method are similar to the results of an open-hole slug test; with the added information from the flowmeter, however, the head and transmissivity of discrete zones also can be determined. The discrete-interval transmissivities ranged from 0.16 to 330 ft2/d. The flowmeter-derived open-hole transmissivity, which is the combined total of each of the transmissive zones, ranged from 1 to 511 ft2/d. The whole-well open-hole transmissivity values determined with HPFM methods were compared to the results of open-hole hydraulic tests. Despite the fact that the flowmeter-derived transmissivities consistently were lower than the estimates derived from open-hole hydraulic tests alone, the correlation was very strong (with a coefficient of determination, R2, of 0.9866), indicating the HPFM method provides a reasonable estimate of transmissivities for the most transmissive fractures in the borehole. Geologic framework, fracture characterization, and estimates of hydraulic properties were interpreted together to characterize the fracture network. The data and interpretation presented in this report should provide information useful for site investigators as the conceptual site groundwater flow model is refined. Collectively, the results and the conceptual site model are important for evaluating remediation options and planning or implementing the design of a well field and borehole completions that will be adequate for monitoring flow, remediation efforts, groundwater levels, and (or) water quality. Similar kinds of borehole geophysical logging (specifically the borehole imaging, gamma, fluid logs, and HPFM) should be conducted in any newly installed boreholes and integrated with interpretations of any nearby boreholes. If boreholes are installed close to existing or other new boreholes, cross-hole flowmeter surveys may be appropriate and may help characterize the aquifer properties and connections between the boreholes.

Maine↗