USGS Science⌕ Search

SEARCH · USGS Science

Results for “Interfaces”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,495 records · Page 83Linked to original sources

Approaches to highly parameterized inversion-A guide to using PEST for groundwater-model calibration

Highly parameterized groundwater models can create calibration difficulties. Regularized inversion-the combined use of large numbers of parameters with mathematical approaches for stable parameter estimation-is becoming a common approach to address these difficulties and enhance the transfer of information contained in field measurements to parameters used to model that system. Though commonly used in other industries, regularized inversion is somewhat imperfectly understood in the groundwater field. There is concern that this unfamiliarity can lead to underuse, and misuse, of the methodology. This document is constructed to facilitate the appropriate use of regularized inversion for calibrating highly parameterized groundwater models. The presentation is directed at an intermediate- to advanced-level modeler, and it focuses on the PEST software suite-a frequently used tool for highly parameterized model calibration and one that is widely supported by commercial graphical user interfaces. A brief overview of the regularized inversion approach is provided, and techniques for mathematical regularization offered by PEST are outlined, including Tikhonov, subspace, and hybrid schemes. Guidelines for applying regularized inversion techniques are presented after a logical progression of steps for building suitable PEST input. The discussion starts with use of pilot points as a parameterization device and processing/grouping observations to form multicomponent objective functions. A description of potential parameter solution methodologies and resources available through the PEST software and its supporting utility programs follows. Directing the parameter-estimation process through PEST control variables is then discussed, including guidance for monitoring and optimizing the performance of PEST. Comprehensive listings of PEST control variables, and of the roles performed by PEST utility support programs, are presented in the appendixes.

Scientific Investigations Report↗

Simulation of specific conductance and chloride concentration in Abercorn Creek, Georgia, 2000-2009

The City of Savannah operates an industrial and domestic water-supply intake on Abercorn Creek approximately 2 miles from the confluence with the Savannah River upstream from the Interstate 95 bridge. Chloride concentrations are a major concern for the city because industrial customers require water with low chloride concentrations, and elevated chloride concentrations require additional water treatment in order to meet those needs. The proposed deepening of Savannah Harbor could increase chloride concentrations (the major ion in seawater) in the upper reaches of the lower Savannah River estuary, including Abercorn Creek. To address this concern, mechanistic and empirical modeling approaches were used to simulate chloride concentrations at the city's intake to evaluate potential effects from deepening the Savannah Harbor. The first approach modified the mechanistic Environmental Fluid Dynamics Code (EFDC) model developed by Tetra Tech and used for evaluating proposed harbor deepening effects for the Environmental Impact Statement. Chloride concentrations were modeled directly with the EFDC model as a conservative tracer. This effort was done by Tetra Tech under a separate funding agreement with the U.S. Army Corps of Engineers and documented in a separate report. The second approach, described in this report, was to simulate chloride concentrations by developing empirical models from the available data using artificial neural network (ANN) and linear regression models. The empirical models used daily streamflow, specific conductance (field measurement for salinity), water temperature, and water color time series for inputs. Because there are only a few data points that describe the relation between high specific conductance values at the Savannah River at Interstate 95 and the water plant intake, there was a concern that these few data points would determine the extrapolation of the empirical model and potentially underestimate the effect of deepening the harbor on chloride concentrations at the intake. To accommodate these concerns, two ANN chloride models were developed for the intake. The first model (ANN M1e) used all the data. The second model (ANN M2e) only used data when specific conductance at Interstate 95 was less than 175 microsiemens per centimeter at 25 degrees Celsius. Deleting the conductivity data greater than 175 microsiemens per centimeter removed the "plateau" effect observed in the data. The chloride simulations with the ANN M1 model have a low sensitivity to specific conductance (salinity) at Interstate 95, whereas the chloride simulations with the ANN M2 model have a high sensitivity to salinity at Interstate 95. The two modeling approaches (Tetra Tech's EFDC model and the one described in this report) were integrated into a decision support system (DSS) that combines the historical database, output from EFDC, ANN models, ANN model simulation controls, streaming graphics, and model output. The DSS was developed as a Microsoft ExcelTM/Visual Basic for Applications program, which allowed the DSS to be prototyped, easily modified, and distributed in a familiar spreadsheet format. The EFDC and ANN models were used to simulate various harbor deepening scenarios. To accommodate the geometry changes in the harbor, the ANN models used the EFDC model-simulated salinity changes for a historical condition as input. The DSS uses a graphical user interface and allows the user to interrogate the ANN models and EFDC output. Two scenarios were simulated using the Savannah Chloride Model DSS to demonstrate different input options. One scenario decreased winter streamflows to a constant streamflow for 45 days. Streamflows during the period January 1 to February 15 were set to a constant 3,600 cubic feet per second for the simulation period of October 1, 2006, to October 1, 2009. The decreased winter streamflow resulted in predictions of increased specific conductance by as much as 50 microsiemens per centimeter and chlorid

Georgia↗

Phase II modification of the W ater A vailability T ool for E nvironmental R esources (WATER) for Kentucky: The sinkhole-drainage process, point-and-click basin delineation, and results of karst test-basin simulations

This report describes Phase II modifications made to the Water Availability Tool for Environmental Resources (WATER), which applies the process-based TOPMODEL approach to simulate or predict stream discharge in surface basins in the Commonwealth of Kentucky. The previous (Phase I) version of WATER did not provide a means of identifying sinkhole catchments or accounting for the effects of karst (internal) drainage in a TOPMODEL-simulated basin. In the Phase II version of WATER, sinkhole catchments are automatically identified and delineated as internally drained subbasins, and a modified TOPMODEL approach (called the sinkhole drainage process, or SDP-TOPMODEL) is applied that calculates mean daily discharges for the basin based on summed area-weighted contributions from sinkhole drain-age (SD) areas and non-karstic topographically drained (TD) areas. Results obtained using the SDP-TOPMODEL approach were evaluated for 12 karst test basins located in each of the major karst terrains in Kentucky. Visual comparison of simulated hydrographs and flow-duration curves, along with statistical measures applied to the simulated discharge data (bias, correlation, root mean square error, and Nash-Sutcliffe efficiency coefficients), indicate that the SDPOPMODEL approach provides acceptably accurate estimates of discharge for most flow conditions and typically provides more accurate simulation of stream discharge in karstic basins compared to the standard TOPMODEL approach. Additional programming modifications made to the Phase II version of WATER included implementation of a point-and-click graphical user interface (GUI), which fully automates the delineation of simulation-basin boundaries and improves the speed of input-data processing. The Phase II version of WATER enables the user to select a pour point anywhere on a stream reach of interest, and the program will automatically delineate all upstream areas that contribute drainage to that point. This capability enables automatic delineation of a simulation basin of any size (area) and having any level of stream-network complexity. WATER then automatically identifies the presence of sinkholes catchments within the simulation basin boundaries; extracts and compiles the necessary climatic, topographic, and basin characteristics datasets; and runs the SDP-TOPMODEL approach to estimate daily mean discharges (streamflow).

Kentucky↗

Evaluation of effects of changes in canal management and precipitation patterns on salinity in Biscayne Bay, Florida, using an integrated surface-water/groundwater model

Biscayne National Park, located in Biscayne Bay in southeast Florida, is one of the largest marine parks in the country and sustains a large natural marine fishery where numerous threatened and endangered species reproduce. In recent years, the bay has experienced hypersaline conditions (salinity greater than 35 practical salinity units) of increasing magnitude and duration. Hypersalinity events were particularly pronounced during April to August 2004 in nearshore areas along the southern and middle parts of the bay. Prolonged hypersaline conditions can cause degradation of water quality and permanent damage to, or loss of, brackish nursery habitats for multiple species of fish and crustaceans as well as damage to certain types of seagrasses that are not tolerant of extreme changes in salinity. To evaluate the factors that contribute to hypersalinity events and to test the effects of possible changes in precipitation patterns and canal flows into Biscayne Bay on salinity in the bay, the U.S. Geological Survey constructed a coupled surface-water/groundwater numerical flow model. The model is designed to account for freshwater flows into Biscayne Bay through the canal system, leakage of salty bay water into the underlying Biscayne aquifer, discharge of fresh and salty groundwater from the Biscayne aquifer into the bay, direct effects of precipitation on bay salinity, indirect effects of precipitation on recharge to the Biscayne aquifer, direct effects of evapotranspiration (ET) on bay salinity, indirect effects of ET on recharge to the Biscayne aquifer, and maintenance of mass balance of both water and solute. The model was constructed using the Flow and Transport in a Linked Overland/Aquifer Density Dependent System (FTLOADDS) simulator, version 3.3, which couples the two-dimensional, surface-water flow and solute-transport simulator SWIFT2D with the density-dependent, groundwater flow an solute-transport simulator SEAWAT. The model was calibrated by a trial-and-error method to fit observed groundwater heads, estimated base flow, and measured bay salinity and temperatures from 1996 to 2004, as well as the location of the freshwater-saltwater interface in the aquifer, by adjusting ET rate limiters, canal vertical hydraulic conductance, leakage rate coefficients (transition-layer thickness and hydraulic conductivity), Manning's n value, and delineation of rainfall zones. Although flow budget calculations indicate that precipitation, ET, and groundwater flux into the bay represent a small portion of the overall budget, these factors may be important in controlling salinity in some parts of the bay, for example the southern parts of the bay where the canal system is not extensively developed or controlled. The balance of precipitation and ET during the wet season generally results in a reduction of bay salinity, whereas the balance of precipitation and ET during the dry season generally results in an increase in bay salinity. During years when wet season precipitation is lower than average, for example less than 70 percent total precipitation for an average year, ET could outweigh precipitation over the bay for essentially the entire year. Hypersaline conditions are prone to occur near the end of the dry season because precipitation rates are generally lower, canal discharge rates (which are strongly correlated to precipitation rates) are also generally lower, and ET rates are higher than during the rest of the year. The hypersalinity event of 2004 followed several years of relatively low precipitation and correspondingly reduced canal structure releases and was unusually extensive, continuing into July. Thus, hypersalinity is ultimately the result of a cumulative deficit of precipitation. The model was used to test the effects of possible changes in canal flux and precipitation. Simulation results showed that by increasing, reducing, or modifying canal discharge rates, the effects on salinity in the bay were more pronounced in the northern part of the bay where there are more canals and canal-control structures. By doubling and halving precipitation, the effects on bay salinity were more pronounced in the southern part of the bay than in the northern part of the bay where there are fewer canals and canal-control structures. The model is designed to quantify factors that contribute to hypersaline conditions in Biscayne Bay and may be less appropriate for addressing other issues or examining conditions substantially different from those described in this report. Model results must be interpreted in light of model limitations, which include representation of the system and conceptual model, uncertainty in physical properties used to describe the system or processes, the scale and discretization of the system, and representation of the boundary conditions.

Florida↗

The development and application of a decision support system for land management in the Lake Tahoe Basin—The Land Use Simulation Model

This report describes and applies the Land Use Simulation Model (LUSM), the final modeling product for the long-term decision support project funded by the Southern Nevada Public Land Management Act and developed by the U.S. Geological Survey’s Western Geographic Science Center for the Lake Tahoe Basin. Within the context of the natural-resource management and anthropogenic issues of the basin and in an effort to advance land-use and land-cover change science, this report addresses the problem of developing the LUSM as a decision support system. It includes consideration of land-use modeling theory, fire modeling and disturbance in the wildland-urban interface, historical land-use change and its relation to active land management, hydrologic modeling and the impact of urbanization as related to the Lahontan Regional Water Quality Control Board’s recently developed Total Maximum Daily Load report for the basin, and biodiversity in urbanizing areas. The LUSM strives to inform land-management decisions in a complex regulatory environment by simulating parcel-based, land-use transitions with a stochastic, spatially constrained, agent-based model. The tool is intended to be useful for multiple purposes, including the multiagency Pathway 2007 regional planning effort, the Tahoe Regional Planning Agency (TRPA) Regional Plan Update, and complementary research endeavors and natural-resource-management efforts. The LUSM is an Internet-based, scenario-generation decision support tool for allocating retired and developed parcels over the next 20 years. Because USGS staff worked closely with TRPA staff and their “Code of Ordinances” and analyzed datasets of historical management and land-use practices, this report accomplishes the task of providing reasonable default values for a baseline scenario that can be used in the LUSM. One result from the baseline scenario for the model suggests that all vacant parcels could be allocated within 12 years. Results also include: assessment of model functionality, brief descriptions of the 7 basic output tables, assessment of the rate of change in land-use allocation pools over time, locations and amounts of the spatially explicit probabilities of land-use transitions by real estate commodity, and analysis of the state change from today’s existing land cover to potential land uses in the future. Assumptions and limitations of the model are presented. This report concludes with suggested next steps to support the continued utility of the LUSM and additional research avenues.

Nevada↗

Two-dimensional simulation of the June 11, 2010, flood of the Little Missouri River at Albert Pike Recreational Area, Ouachita National Forest, Arkansas

In the early morning hours of June 11, 2010, substantial flooding occurred at Albert Pike Recreation Area in the Ouachita National Forest of west-central Arkansas, killing 20 campers. The U.S. Forest Service needed information concerning the extent and depth of flood inundation, the water velocity, and flow paths throughout Albert Pike Recreation Area for the flood and for streamflows corresponding to annual exceedence probabilities of 1 and 2 percent. The two-dimensional flow model Fst2DH, part of the Federal Highway Administration’s Finite Element Surface-water Modeling System, and the graphical user interface Surface-water Modeling System (SMS) were used to perform a steady-state simulation of the flood in a 1.5-mile reach of the Little Missouri River at Albert Pike Recreation Area. Peak streamflows of the Little Missouri River and tributary Brier Creek served as inputs to the simulation, which was calibrated to the surveyed elevations of high-water marks left by the flood and then used to predict flooding that would result from streamflows corresponding to annual exceedence probabilities of 1 and 2 percent. The simulated extent of the June 11, 2010, flood matched the observed extent of flooding at Albert Pike Recreation Area. The mean depth of inundation in the camp areas was 8.5 feet in Area D, 7.4 feet in Area C, 3.8 feet in Areas A, B, and the Day Use Area, and 12.5 feet in Lowry’s Camp Albert Pike. The mean water velocity was 7.2 feet per second in Area D, 7.6 feet per second in Area C, 7.2 feet per second in Areas A, B, and the Day Use Area, and 7.6 feet per second in Lowry’s Camp Albert Pike. A sensitivity analysis indicated that varying the streamflow of the Little Missouri River had the greatest effect on simulated water-surface elevation, while varying the streamflow of tributary Brier Creek had the least effect. Simulated water-surface elevations were lower than those modeled by the U.S. Forest Service using the standard-step method, but the comparison between the two was favorable with a mean absolute difference of 0.58 feet in Area C and 0.32 feet in Area D. Results of a HEC-RAS model of the Little Missouri River watershed upstream from the U.S. Geological Survey streamflow-gaging station near Langley showed no difference in mean depth in the areas in common between the models, and a difference in mean velocity of only 0.5 foot per second. Predictions of flooding that would result from streamflows corresponding to annual exceedence probabilities of 1 and 2 percent indicated that the extent of inundation of the June 11, 2010, flood exceeded that of the 1 percent flood, and that for both the 1 and 2 percent floods, all of Areas C and D, and parts of Areas A, B, and the Day Use Area were inundated. Predicted water-surface elevations for the 1 and 2 percent floods were approximately 1 foot lower than those predicted by the U.S. Forest Service using a standard-step model.

Arkansas↗

Simulation of salinity intrusion along the Georgia and South Carolina coasts using climate-change scenarios

Potential changes in climate could alter interactions between environmental and societal systems and adversely affect the availability of water resources in many coastal communities. Changes in streamflow patterns in conjunction with sea-level rise may change the salinity-intrusion dynamics of coastal rivers. Several municipal water-supply intakes are located along the Georgia and South Carolina coast that are proximal to the present day saltwater-freshwater interface of tidal rivers. Increases in the extent of salinity intrusion resulting from climate change could threaten the availability of freshwater supplies in the vicinity of these intakes. To effectively manage these supplies, water-resource managers need estimates of potential changes in the frequency, duration, and magnitude of salinity intrusion near their water-supply intakes that may occur as a result of climate change. This study examines potential effects of climate change, including altered streamflow and sea-level rise, on the dynamics of saltwater intrusion near municipal water-supply intakes in two coastal areas. One area consists of the Atlantic Intracoastal Waterway (AIW) and the Waccamaw River near Myrtle Beach along the Grand Strand of the South Carolina Coast, and the second area is on or near the lower Savannah River near Savannah, Georgia. The study evaluated how future sea-level rise and a reduction in streamflows can potentially affect salinity intrusion and threaten municipal water supplies and the biodiversity of freshwater tidal marshes in these two areas. Salinity intrusion occurs as a result of the interaction between three principal forces—streamflow, mean coastal water levels, and tidal range. To analyze and simulate salinity dynamics at critical coastal gaging stations near four municipal water-supply intakes, various data-mining techniques, including artificial neural network (ANN) models, were used to evaluate hourly streamflow, salinity, and coastal water-level data collected over a period exceeding 10 years. The ANN models were trained (calibrated) to learn the specific interactions that cause salinity intrusions, and resulting models were able to accurately simulate historical salinity dynamics in both study areas. Changes in sea level and streamflow quantity and timing can be simulated by the salinity intrusion models to evaluate various climate-change scenarios. The salinity intrusion models for the study areas are deployed in a decision support system to facilitate the use of the models for management decisions by coastal water-resource managers. The report describes the use of the salinity-intrusion models decision support system to evaluate salinity-intrusion dynamics for various climate-change scenarios, including incremental increases in sea level in combination with incremental decreases in streamflow. Operation of municipal water-treatment plants is problematic when the specific-conductance values for source water are greater than 1,000 to 2,000 microsiemens per centimeter (µS/cm). High specific-conductance values contribute to taste problems that require treatment. Data from a gage downstream from a municipal water intake indicate specific conductance exceeded 1,000 µS/cm about 5.4 percent of the time over the 14-year period from August 1995 to August 2008. Simulations of specific conductance at this gaging station that incorporates sea-level rises resulted in a doubling of the exceedances to 11.0 percent for a 1-foot increase and 17.6 percent for a 2-foot increase. The frequency of intrusion of water with specific conductance values of 1,000 µS/cm was less sensitive to incremental reductions in streamflow than to incremental increases in sea level. Simulations of conditions associated with a 10-percent reduction in streamflow, in combination with a 1-foot rise in sea level, increased the percentage of time specific conductance exceeded 1,000 µS/cm at this site from 11.0 to 13.3 percent, and a 20-percent reduction in streamflow increased the percentage of time to 16.6 percent. Precipitation and temperature data from a global circulation model were used, after scale adjustments, as input to a watershed model of the Yadkin-Pee Dee River basin, which flows into the Waccamaw River and Atlantic Intracoastal Waterway study area in South Carolina. The simulated streamflow for historical conditions and projected climate change in the future was used as input for the ANN model in decision support system. Results of simulations incorporating climate-change projections for alterations in streamflow indicate an increase in the frequency of salinity-intrusion events and a shift in the seasonal occurrence of the intrusion events from the summer to the fall.

South Carolina↗

Plankton communities and summertime declines in algal abundance associated with low dissolved oxygen in the Tualatin River, Oregon

Phytoplankton populations in the Tualatin River in northwestern Oregon are an important component of the dissolved oxygen (DO) budget of the river and are critical for maintaining DO levels in summer. During the low-flow summer period, sufficient nutrients and a long residence time typically combine with ample sunshine and warm water to fuel blooms of cryptophyte algae, diatoms, green and blue-green algae in the low-gradient, slow-moving reservoir reach of the lower river. Algae in the Tualatin River generally drift with the water rather than attach to the river bottom as a result of moderate water depths, slightly elevated turbidity caused by suspended colloidal material, and dominance of silty substrates. Growth of algae occurs as if on a “conveyor belt” of streamflow, a dynamic system that is continually refreshed with inflowing water. Transit through the system can take as long as 2 weeks during the summer low-flow period. Photosynthetic production of DO during algal blooms is important in offsetting oxygen consumption at the sediment-water interface caused by the decomposition of organic matter from primarily terrestrial sources, and the absence of photosynthesis can lead to low DO concentrations that can harm aquatic life. The periods with the lowest DO concentrations in recent years (since 2003) typically occur in August following a decline in algal abundance and activity, when DO concentrations often decrease to less than State standards for extended periods (nearly 80 days). Since 2003, algal populations have tended to be smaller and algal blooms have terminated earlier compared to conditions in the 1990s, leading to more frequent declines in DO to levels that do not meet State standards. This study was developed to document the current abundance and species composition of phytoplankton in the Tualatin River, identify the possible causes of the general decline in algae, and evaluate hypotheses to explain why algal blooms diminish in midsummer. Plankton and water-quality sample data from 2006 to 2008 were combined with parts of a larger discrete-sample and continuous water-quality monitoring dataset and examined to identify patterns in water-quality and algal conditions since 1991, with a particular emphasis on 2003–08. Longitudinal plankton surveys were conducted in 2006–08 at six sites between river miles (RM) 24.5 and 3.4 at 2- to 3-week intervals, or 5–6 per season, and in-situ bioassay experiments were conducted in 2008 to examine the potential effects of wastewater treatment facility (WWTF) effluent and phosphorus additions on phytoplankton biomass and algal photosynthesis. Phytoplankton and zooplankton community composition, streamflow, and water-quality data were analyzed using multivariate statistical techniques to gain insights into plankton dynamics to determine what factors might be most tied to the abundance and characteristics of the phytoplankton assemblages, and identify possible causes of their declines. The connection between low-DO events and algal declines was clearly evident, as bloom crashes were nearly always followed by periods of low DO. Algal blooms occurred each year during 2006–08, producing maximum chlorophyll-a (Chl-a) values in June or July generally in the range of 50–80 micrograms per liter (µg/L). Bloom crashes and absence of sufficient algal photosynthesis in mid- to late-summer contributed to minimum DO concentrations that were less than the State standard of 6.5 milligrams per liter (mg/L) based on the 30-day mean daily concentration, for 62–74 days each year. At times, the absolute minimum State standard (4 mg/L DO) also was not met. To learn more about why low-DO events occurred, specific algal declines during 2003–08 were scrutinized to determine their likely causal factors. From this information, a series of hypotheses were formulated and evaluated in terms of their ability to explain recent declines in algal populations in the river in late summer. Meteorological, streamflow, turbidity, water temperature, and conductance conditions in the Tualatin River during the 2006–08 summer seasons were not atypical. Natural flow comprised the majority (70–80 percent) of flow each year during spring, but then reduced to 38–40 percent during midsummer when WWTF effluent—which contributed as much as 36 percent—and flow augmentation releases comprised a greater fraction of the flow. Summer 2008 was unusual, however, in the prolonged influence from the Wapato Lake agricultural area near Gaston in the upper part of the basin. The previous winter flooding and levee breach at Wapato Lake caused a much greater area of inundation. As a result, drainage from this area continued into July, much later than normal. A subsequent algal bloom in Wapato Lake then seeded the upper Tualatin River, and this drainage had a profound effect on the downstream plankton community. A large blue-green algae bloom developed—the largest in recent memory—prompting a public health advisory for recreational contact for about two weeks. Algal growths and surface blooms are a common feature of the Tualatin River. Most of the dominant algae have growth forms and morphologies that are well suited for planktonic life, employing spines and gas vacuoles to resist settling, forming colonies, and producing mucilage (or toxins) to resist zooplankton grazing. In 2006–08, 143 algal taxa were identified in 117 main-stem samples; diatoms and green algae were more diverse than blue-green, golden, and cryptophyte algae, although these later groups sometimes dominated the overall volumetric abundance (biovolume). The most frequently occurring taxa, occurring in 97–99 percent of samples, were flagellated cryptophytes Cryptomonas erosa and Rhodomonas minuta. Other important algal taxa included centric diatoms Stephanodiscus, Cyclotella, and Melosira species and colonial green algae Scenedesmus and Actinastrum. These taxa comprised the majority of the algal biovolume during much of the growing season. A general seasonal trend in the phytoplankton assemblages was observed, with dominance by filamentous centric diatoms Stephanodiscus and Melosira in spring and early summer, and flagellated cryptophytes and green algae, particularly Chlamydomonas sp., in late-summer; or, in 2008, dominance by blue-green algae Anabaena flos-aquae and Aphanizomenon flos-aquae during the Wapato Lake bloom event. There were 99 zooplankton taxa identified from the Tualatin River in 2006–08, composed primarily of cladocerans, copepods, and rotifers. A seasonal increase in zooplankton abundance was observed in early summer just as or shortly after the phytoplankton population began to increase, with populations growing to 15,000−120,000 organisms per cubic meter in the lower river. Zooplankton abundance showed a predictable and distinct longitudinal downstream increase, particularly downstream of Highway 99W (RM 11.6). Although grazing rates were not measured, the data suggest that, at times, zooplankton grazing may affect algal abundance and species composition in the Tualatin River, with diatoms becoming relatively less abundant and flagellated cryptophytes and green algae relatively more abundant during periods when zooplankton densities were highest. Multivariate statistical analyses identified soluble reactive phosphorus (SRP), natural flow, flow augmentation, and WWTF effluent as important factors influencing Tualatin River phytoplankton populations, with zooplankton density (particularly rotifers and copepods), specific conductance, chloride, and water temperature also having an important influence. Although SRP was highly correlated with the plankton communities, that correlation was likely the result of high or low algal activity (uptake) as SRP concentrations were often reduced to low levels during blooms. While previous studies have already established that phosphorus, among other factors such as flow, places a theoretical cap on the size of the phytoplankton population in the river, sometimes algal declines occur when SRP concentrations are apparently sufficient. To identify alternative causal factors, additional analyses were performed without SRP to focus on other water-quality parameters, zooplankton density, and flow factors. Considering data for all 3 years and including just those samples from the lower Tualatin River not affected by the 2008 Wapato Lake drainage event, three factors (percentage of reservoir flow augmentation, total natural flow, and rotifer density) best explained variations in the phytoplankton assemblages. Analyses focusing on the possible causes of algal declines included the above multivariate analyses, scrutiny of 10 specific instances of declines in algal populations during 2003–08 including several bloom–crash sequences, and analyses of historic routine watershed monitoring data from Clean Water Services. Six factors were hypothesized to be important in causing bloom crashes or impeding blooms from rebounding in August: (1) light limitation from cloudy weather, (2) a reduction in the plankton inocula or “seed” entering the lower river from upstream sources, (3) increased summer streamflows, (4) changes in the dominant sources of flow as the percentage of flow augmentation and WWTF discharges have increased, (5) zooplankton grazing, and (6) low concentrations of bioavailable phosphorus (<0.015 milligram per liter). All of these hypotheses are supported in some fashion by the available data and statistical analyses. Zooplankton grazing, short-term declines in photosynthesis from cloudy weather, total flow as it affects residence time, and the dominant source of flow are primary factors responsible for the low-DO events caused by declines in algae in the lower Tualatin River during late summer. Cloudy weather and increased turbidity are known to inhibit algal growth in the Tualatin River, and slight increases in turbidity in recent years may be a problem. Upstream sources of algae are critical in determining the characteristics and size of downstream populations, as illustrated by the Wapato Lake bloom in 2008, but more data are needed from upstream to fully define the importance of this connection. The sources of flow, through their differential contribution of plankton inocula (quality and amount), were, at times, important factors affecting phytoplankton populations. While SRP concentrations were often most highly correlated with phytoplankton species community, the bioavailability of phosphorus is still somewhat unknown and there are several sources to consider. Preliminary bioassay tests suggested that while treated wastewater effluent may stimulate algae at 30 percent concentrations, negative effects (or decreased stimulation) on Chl-a and DO production may occur at concentrations of 50 percent. Targeted data collection and future research will be needed to further understand the importance of these factors on Tualatin River phytoplankton. While the data and analysis completed for this report provide insights into future research and monitoring that would be useful to continue, additional monitoring of turbidity, Chl-a, and plankton abundance and species composition in the upper part of the basin would enhance our understanding of plankton dynamics and factors affecting phytoplankton abundance in the lower river. Assessment of the key upstream sources of algal inocula via surveys of the major flow sources as well as tributaries and wetlands would provide useful information for the management of river water quality. Other studies that could prove useful for developing management strategies include targeted experiments to evaluate the bioavailability of phosphorus from a variety of sources. New research on phytoplankton–zooplankton interactions, and studies of planktivorous fish, might also provide insight about food web dynamics and potential “top-down” effects of fish predation on the plankton communities. In addition, further development of neural-network or other water-quality models would help to evaluate management strategies and provide forecasts of water-quality conditions. Finally, periodic future reassessments of the available data with the multivariate statistical tools used in this study would be helpful to assess whether and how plankton communities are changing, and to continue to shed light on the importance of factors shaping the plankton. Although certain types and sizes of algal blooms are undesirable, minimum phytoplankton populations are an important part of aquatic food webs and are needed to maintain healthy levels of DO in the river. By understanding the sources, characteristics, causal factors, and responses of the plankton communities, management strategies can be developed to improve DO conditions in the lower Tualatin River during the important summer low-flow period.

Oregon↗

Hydrogeology, groundwater seepage, nitrate distribution, and flux at the Raleigh hydrologic research station, Wake County, North Carolina, 2005-2007

rom 2005 to 2007, the U.S. Geological Survey and the North Carolina Department of Environment and Natural Resources, Division of Water Quality, conducted a study to describe the geologic framework, measure groundwater quality, characterize the groundwater-flow system, and describe the groundwater/surface-water interaction at the 60-acre Raleigh hydrogeologic research station (RHRS) located at the Neuse River Waste Water Treatment Plant in eastern Wake County, North Carolina. Previous studies have shown that the local groundwater quality of the surficial and bedrock aquifers at the RHRS had been affected by high levels of nutrients. Geologic, hydrologic, and water-quality data were collected from 3 coreholes, 12 wells, and 4 piezometers at 3 well clusters, as well as from 2 surface-water sites, 2 multiport piezometers, and 80 discrete locations in the streambed of the Neuse River. Data collected were used to evaluate the three primary zones of the Piedmont aquifer (regolith, transition zone, and fractured bedrock) and characterize the interaction of groundwater and surface water as a mechanism of nutrient transport to the Neuse River. A conceptual hydrogeologic cross section across the RHRS was constructed using new and existing data. Two previously unmapped north striking, nearly vertical diabase dikes intrude the granite beneath the site. Groundwater within the diabase dike appeared to be hydraulically isolated from the surrounding granite bedrock and regolith. A correlation exists between foliation and fracture orientation, with most fractures striking parallel to foliation. Flowmeter logging in two of the bedrock wells indicated that not all of the water-bearing fractures labeled as water bearing were hydraulically active, even when stressed by pumping. Groundwater levels measured in wells at the RHRS displayed climatic and seasonal trends, with elevated groundwater levels occurring during the late spring and declining to a low in the late fall. Vertical gradients in the groundwater discharge area near the Neuse River were complex and were affected by fluctuations in river stage, with the exception of a well completed in a diabase dike. Water-quality data from the wells and surface-water sites at the RHRS were collected continuously as well as during periodic sampling events. Surface-water samples collected from a tributary were most similar in chemical composition to groundwater found in the regolith and transition zone. Nitrate (measured as nitrite plus nitrate, as nitrogen) concentrations in the sampled wells and tributary ranged from about 5 to more than 120 milligrams per liter as nitrogen. Waterborne continuous resistivity profiling conducted on the Neuse River in the area of the RHRS measured areas of low apparent resistivity that likely represent groundwater contaminated by high concentrations of nitrate. These areas were located on either side of a diabase dike and at the outfall of two unnamed tributaries. The diabase dike preferentially directed the discharge of groundwater to the Neuse River and may isolate groundwater movement laterally. Discrete temperature measurements made within the pore water beneath the Neuse River revealed seeps of colder groundwater discharging into warmer surface water near a diabase dike. Water-quality samples collected from the pore water beneath the Neuse River indicated that nitrate was present at concentrations as high as 80 milligrams per liter as nitrogen on the RHRS side of the river. The highest concentrations of nitrate were located within pore water collected from an area near a diabase dike that was identified as a suspected seepage area. Hydraulic head was measured and pore water samples were collected from two 140-centimeter-deep (55.1-inch-deep) multiport piezometers that were installed in bed sediments on opposite sides of a diabase dike. The concentration of nitrate in pore water at a suspected seepage area ranged from 42 to 82 milligrams per liter as nitrogen with a median concentration of 79 milligrams per liter as nitrogen. On the opposite side of the dike, concentrations of nitrate in pore water samples ranged from 3 to 91 milligrams per liter as nitrogen with a median concentration of 52 milligrams per liter. At one of the multiport piezometers the vertical gradient of hydraulic head between the Neuse River and the groundwater was too small to measure. At the multiport piezometer located in the suspected seepage area, an upward gradient of about 0.1 was present and explains the occurrence of higher concentrations of nitrate near the sediment/water interface. Horizontal seepage flux from the surficial aquifer to the edge of the Neuse River was estimated for 2006. Along a 130-foot flow path, the estimated seepage flux ranged from –0.52 to 0.2 foot per day with a median of 0.09 foot per day. The estimated advective horizontal mass flux of nitrate along a 300-foot reach of the Neuse River ranged from –10.9 to 5 pounds per day with a median of 2.2 pounds per day. The total horizontal mass flux of nitrate from the surficial aquifer to the Neuse River along the 130-foot flow path was estimated to be about 750 pounds for all of 2006. Seepage meters were deployed on the bed of the Neuse River in the areas of the multiport piezometers on either side of the diabase dike to estimate rates of vertical groundwater discharge and flux of nitrate. The average estimated daily seepage flux differed by two orders of magnitude between seepage areas. The potential vertical flux of nitrate from groundwater to the Neuse River was estimated at an average of 2.5 grams per day near one of the multiport piezometers and an average of 784 grams per day at the other. These approximations suggest that under some hydrologic conditions there is the potential for substantial quantities of nitrate to discharge from the groundwater to the Neuse River.

North Carolina↗

An evaluation of seepage gains and losses in Indian Creek Reservoir, Ada County, Idaho, April 2010–November 2011

The U.S. Geological Survey, in cooperation with the Idaho Department of Water Resources, conducted an investigation on Indian Creek Reservoir, a small impoundment in east Ada County, Idaho, to quantify groundwater seepage into and out of the reservoir. Data from the study will assist the Idaho Water Resources Department’s Comprehensive Aquifer Management Planning effort to estimate available water resources in Ada County. Three independent methods were utilized to estimate groundwater seepage: (1) the water-budget method; (2) the seepage-meter method; and (3) the segmented Darcy method. Reservoir seepage was quantified during the periods of April through August 2010 and February through November 2011. With the water-budget method, all measureable sources of inflow to and outflow from the reservoir were quantified, with the exception of groundwater; the water-budget equation was solved for groundwater inflow to or outflow from the reservoir. The seepage-meter method relies on the placement of seepage meters into the bottom sediments of the reservoir for the direct measurement of water flux across the sediment-water interface. The segmented-Darcy method utilizes a combination of water-level measurements in the reservoir and in adjacent near-shore wells to calculate water-table gradients between the wells and the reservoir within defined segments of the reservoir shoreline. The Darcy equation was used to calculate groundwater inflow to and outflow from the reservoir. Water-budget results provided continuous, daily estimates of seepage over the full period of data collection, while the seepage-meter and segmented Darcy methods provided instantaneous estimates of seepage. As a result of these and other difference in methodologies, comparisons of seepage estimates provided by the three methods are considered semi-quantitative. The results of the water-budget derived estimates of seepage indicate seepage to be seasonally variable in terms of the direction and magnitude of flow. The reservoir tended to gain water from seepage of groundwater in the early spring months (March–May), while seepage losses to groundwater from the reservoir occurred in the drier months (June–October). Net monthly seepage rates, as computed by the water-budget method, varied greatly. Reservoir gains from seepage ranged from 0.2 to 59.4 acre-feet per month, while reservoir losses to seepage ranged from 1.6 and 26.8 acre-feet per month. An analysis of seepage meter estimates and segmented-Darcy estimates qualitatively supports the seasonal patterns in seepage provided by the water-budget calculations, except that they tended to be much smaller in magnitude. This suggests that actual seepage might be smaller than those estimates made by the water-budget method. Although the results of all three methods indicate that there is some water loss from the reservoir to groundwater, the seepage losses may be due to rewetting of unsaturated near-shore soils, possible replenishment of a perched aquifer, or both, rather than through percolation to the local aquifer that lies 130 feet below the reservoir. A lithologic log from an adjacent well indicates the existence of a clay lithology that is well correlated to the original reservoir’s base elevation. If the clay lithologic unit extends beneath the reservoir basin underlying the fine-grain reservoir bed sediments, the clay layer should act as an effective barrier to reservoir seepage to the local aquifer, which would explain the low seepage loss estimates calculated in this study.

Idaho↗

Settlement of the USS Arizona, Pearl Harbor, Hawaii

The U.S. Geological Survey, in collaboration with the National Park Service Submerged Resources Center, undertook investigations at the USS Arizona Memorial at Pearl Harbor, Hawaii, in 2002, 2003, and 2005 to characterize geological factors affecting the deterioration and movement of the hull of the USS Arizona. Since sinking on the morning of December 7, 1941, the hull of the USS Arizona has been slowly but steadily disappearing below the surface of Pearl Harbor. Continuous sediment coring at three of four locations around the hull of the Arizona was only partially successful, but it was sufficient to identify a varied sedimentary substrate beneath the hull. A boring near the stern reveals a thick, continuous sequence of soft, gray clay to the bottom of the boring. In contrast, borings near the bow and starboard side, below about 5 meters subbottom depth, indicate the presence of very stiff, brown clay and coral debris and an absence of soft clay. Multisensor core logger scanning of the recovered cores distinguishes the lower density of the soft, gray clay at the stern from the higher density of the stiff, brown clays and coral debris at the bow and starboard side. Uniaxial consolidation testing of the soft gray clay indicates a normally consolidated sequence, whereas the stiff, brown clay and coral debris are overconsolidated. Profiles of shear wave velocity vs. depth obtained through spectral analysis of interface wave testing around the perimeter of the hull in 2005 identified areas of higher velocity, stiffer sediment at the bow and starboard side, which correspond to the dense, stiff clay recovered near the bow and starboard borings. Low shear-wave velocities at the port midship and quarter of the hull correlate with the lower density, softer sediment recovered from the boring at the stern. Cross sections of the subbottom of the Memorial combine results from the sediment borings and geophysical surveys and depict a wedge of soft clay unconformably overlying the stiff clays and coral debris beneath the aft half of the USS Arizona and thickening toward the stern. The 2008 position of the hull has been documented using both tide-based and differential Global Positioning System (GPS) measuring systems. Analysis of historical and recent photographs was done to create a record of settlement from the time of sinking in 1941 to the present. By examining shadows in suitable photos, the sun azimuth, local time of day, and tide levels were determined to derive tide-adjusted and sea-level-rise-corrected elevations for structures on the hull and from these elevations to obtain settlement and tilt trends. The settlement trends, most complete for barbette 3, have two components. An early, nonlinear component ends on December 9, 1941, and represents the initial penetration and displacement of the bottom sediment by the hull. A linear, long-term trend of normal consolidation continues to the present day. Long-term settlement rates are greatest at the stern and decrease linearly to the midship, showing that the aft half of the hull is moving as an intact, rigid body. The recent rate of settlement at the stern is about 3.5 mm/year; rates at the starboard midship and forward part of the hull are less than one-third of the stern rate. The aft half of the USS Arizona hull presently tilts about 2 degrees to port, an increase of at least 1.5 degrees since the initial sinking of the ship. The results of this study identify differential settlement of the Arizona hull, due to the wedge of soft clay underlying the aft half of the hull, as the cause of the movement of the hull beneath the surface of Pearl Harbor. Calculation of sediment consolidation using lab-determined properties of the soft clay demonstrates that the observed settlements can be reproduced by projecting appropriate clay thicknesses beneath the hull. Several of the high-quality photographs analyzed for the historical settlement analysis highlight some of the limitations of this retrospective technique for determining tide-based elevations. In these cases, calculated structure elevations do not conform to the settlement trend, indicating that there can be complicating factors affecting the interpretation of the photos. Conflicting dates for events during the salvage operations were also encountered.

Hawai'i↗

Geohydrology, geochemistry, and groundwater simulation (1992-2011) and analysis of potential water-supply management options, 2010-60, of the Langford Basin, California

Groundwater withdrawals began in 1992 from the Langford Basin within the Fort Irwin National Training Center (NTC), California. From April 1992 to December 2010, approximately 12,300 acre-feet of water (averaging about 650 acre-feet per year) has been withdrawn from the basin and transported to the adjacent Irwin Basin. Since withdrawals began, water levels in the basin have declined by as much as 40 feet, and the quality of the groundwater withdrawn from the basin has deteriorated. The U.S. Geological Survey collected geohydrologic data from Langford Basin during 1992–2011 to determine the quantity and quality of groundwater available in the basin. Geophysical surveys, including gravity, seismic refraction, and time-domain electromagnetic induction surveys, were conducted to determine the depth and shape of the basin, to delineate depths to the Quaternary-Tertiary interface, and to map the depth to the water table and changes in water quality. Data were collected from existing wells and test holes, as well as 11 monitor wells that were installed at 5 sites as part of this study. Water-quality samples collected from wells in the basin were used to determine the groundwater chemistry within the basin and to delineate potential sources of poor-quality groundwater. Analysis of stable isotopes of oxygen and hydrogen in groundwater indicates that present-day precipitation is not a major source of recharge to the basin. Tritium and carbon-14 data indicate that most of the basin was recharged prior to 1952, and the groundwater in the basin has an apparent age of 12,500 to 30,000 years. Recharge to the basin, estimated to be less than 50 acre-feet per year, has not been sufficient to replenish the water that is being withdrawn from the basin. A numerical groundwater-flow model was developed for the Langford Basin to better understand the aquifer system used by the Fort Irwin NTC as part of its water supply, and to provide a tool to help manage groundwater resources at the NTC. Measured groundwater-level declines since the initiation of withdrawals (1992–2011) were used to calibrate the groundwater-flow model. The simulated recharge was about 46 acre-feet per year, including approximately 6 acre-feet per year of natural recharge derived from precipitation runoff and as much as 40 acre-feet per year of underflow from the Irwin Basin. Between April 1992 and December 2010, an average of about 650 acre-feet per year of water was withdrawn from the Langford Basin. Groundwater withdrawals in excess of natural recharge resulted in a net loss of 11,670 acre-feet of groundwater storage within the basin for the simulation period. The Fort Irwin NTC is considering various groundwater-management options to address the limited water resources in the Langford Basin. The calibrated Langford Basin groundwater-flow model was used to evaluate the hydrologic effects of four groundwater-withdrawal scenarios being considered by the Fort Irwin NTC over the next 50 years (January 2011 through December 2060). Continuation of the 2010 withdrawal rate in the three existing production wells will result in 70 feet of additional drawdown in the central part of the basin. Redistributing the 2010 withdrawal rate equally to the three existing wells and two proposed new wells in the northern and southern parts of the basin would result in about 10 feet less drawdown in the central part of the basin but about 100 feet of additional drawdown in the new well in the northern part of the basin and about 50 feet of additional drawdown in the new well in the southern part of the basin. Reducing the withdrawals from the three existing production wells in the central part of the basin from about 45,000 acre-feet to about 32,720 acre-feet would result in about 40 feet of additional drawdown in the central basin near the pumping wells, about 25 feet less than if withdrawals were not reduced. The combination of reducing and redistributing the cumulative withdrawals to the three existing and two proposed new wells results in about 40 feet of additional drawdown in the central and southern parts of the basin and about 70 feet in the northern part of the basin. These results show that reducing and redistributing the groundwater withdrawals would maintain the upper aquifer at greater than 50 percent of its predevelopment saturated thickness throughout the groundwater basin. The scenarios simulated for this study demonstrate how the calibrated model can be utilized to evaluate the hydrologic effects of different water-management strategies.

California↗

Simulation of groundwater flow in the "1,500-foot" sand and "2,000-foot" sand, with scenarios to mitigate saltwater migration in the "2,000-foot" sand of the Baton Rouge area, Louisiana

Groundwater withdrawals have caused saltwater to encroach into freshwater-bearing aquifers beneath Baton Rouge, Louisiana. Groundwater investigations in the 1960s identified a freshwater-saltwater interface located at the Baton Rouge Fault, across which abrupt changes in water levels occur. Aquifers south of the fault generally contain saltwater, and aquifers north of the fault contain freshwater, though limited saltwater encroachment has been detected within 7 of the 10 aquifers north of the fault. The 10 aquifers beneath the Baton Rouge area, which includes East and West Baton Rouge Parishes, Pointe Coupee Parish, and East and West Feliciana Parishes, provided about 167 million gallons per day (Mgal/d) for public supply and industrial use in 2010. Groundwater withdrawals from the “2,000-foot” sand in East Baton Rouge Parish have caused water-level drawdown as great as 356 feet (ft) and induced saltwater movement northward across the fault. Saltwater encroachment threatens industrial wells that are located about 3 miles north of the fault. Constant and variable-density groundwater models were developed with the MODFLOW and SEAWAT groundwater modeling codes to evaluate strategies to control saltwater migration, including changes in the distribution of groundwater withdrawals and installation of “scavenger” wells to intercept saltwater before it reaches existing production wells. Six hypothetical scenarios simulated the effects of different groundwater withdrawal options on groundwater levels within the “1,500-foot” sand and the “2,000-foot” sand and the transport of saltwater within the “2,000-foot” sand during 2008–47. Scenario 1 is considered a base case for comparison to the other five scenarios and simulates continuation of 2007 reported groundwater withdrawals. Scenario 2 simulates discontinuation of withdrawals from seven selected industrial wells located in the northwest corner of East Baton Rouge Parish and predicts that water levels within the “1,500-foot” sand will be about 10 to 12 ft higher with this withdrawal reduction than under scenario 1. Scenario 3 simulates the effects of a scavenger well on water levels and chloride concentrations within the “2,000-foot” sand. The scavenger well, which withdraws water from the base of the “2,000-foot” sand at a rate of 2.0 Mgal/d, is simulated at two possible locations. In comparison to the concentrations simulated in scenario 1, operation of the scavenger well at the locations specified in scenario 3 reduces the chloride concentrations at all existing chloride-observation well locations. Scenario 4 simulates a 3.6 Mgal/d reduction in total groundwater withdrawals from selected wells screened in the “2,000-foot” sand that are located in the Baton Rouge industrial district. Under scenario 4, chloride concentrations decrease in the leading portion of the plume south of the industrial district but increase in areas farther east. Scenario 5 simulates the effects of total cessation of withdrawals from the “2,000-foot” sand in the industrial district, which causes a change in the groundwater-flow direction toward municipal supply wells and increased chloride concentrations in the area where municipal supply wells are located. Scenario 6 simulates the combined effect of withdrawal reductions from the “2,000-foot” sand and operation of a scavenger well and was most effective at decreasing the size of the plume area and median and mean chloride concentrations within the “2000-foot” sand in the Baton Rouge area.

Louisiana↗

Origins and delineation of saltwater intrusion in the Biscayne aquifer and changes in the distribution of saltwater in Miami-Dade County, Florida

Intrusion of saltwater into parts of the shallow karst Biscayne aquifer is a major concern for the 2.5 million residents of Miami-Dade County that rely on this aquifer as their primary drinking water supply. Saltwater intrusion of this aquifer began when the Everglades were drained to provide dry land for urban development and agriculture. The reduction in water levels caused by this drainage, combined with periodic droughts, allowed saltwater to flow inland along the base of the aquifer and to seep directly into the aquifer from the canals. The approximate inland extent of saltwater was last mapped in 1995. An examination of the inland extent of saltwater and the sources of saltwater in the aquifer was completed during 2008–2011 by using (1) all available salinity information, (2) time-series electromagnetic induction log datasets from 35 wells, (3) time-domain electromagnetic soundings collected at 79 locations, (4) a helicopter electromagnetic survey done during 2001 that was processed, calibrated, and published during the study, (5) cores and geophysical logs collected from 8 sites for stratigraphic analysis, (6) 8 new water-quality monitoring wells, and (7) analyses of 69 geochemical samples. The results of the study indicate that as of 2011 approximately 1,200 square kilometers (km 2 ) of the mainland part of the Biscayne aquifer were intruded by saltwater. The saltwater front was mapped farther inland than it was in 1995 in eight areas totaling about 24.1 km 2 . In many of these areas, analyses indicated that saltwater had encroached along the base of the aquifer. The saltwater front was mapped closer to the coast than it was in 1995 in four areas totaling approximately 6.2 km 2 . The changes in the mapped extent of saltwater resulted from improved spatial information, actual movement of the saltwater front, or a combination of both. Salinity monitoring in some of the canals in Miami-Dade County between 1988 and 2010 indicated influxes of saltwater, with maximum salinities ranging from 1.4 to 32 practical salinity units (PSU) upstream of the salinity control structures. Time-series electromagnetic induction log data from monitoring wells G–3601, G–3608, and G–3701, located adjacent to the Biscayne, Snapper Creek, and Black Creek Canals, respectively, and upstream of the salinity control structures, indicated shallow influxes of conductive water in the aquifer that likely resulted from leakage of brackish water or saltwater from these canals. The determination that saltwater influxes were recent is supported by the similarity in the oxygen and hydrogen stable isotope composition in samples from the Snapper Creek Canal, 1.6 kilometers (km) inland of a salinity control structure, and in samples from well G–3608, which is adjacent to the canal, as well as by the relative ages of the water sampled from well G–3608 and other wells open to the aquifer below the saltwater interface. Historical and recent salinity information from the Card Sound Road Canal, monitoring well FKS8 located adjacent to the canal, and the 2001 helicopter electromagnetic survey indicated that saltwater may occasionally leak from this canal as far inland as 15 km. This leakage may be prevented or reduced by a salinity control structure that was installed in May 2010. Saltwater also may have leaked from the Princeton Canal. Results of geochemical sampling and analysis indicate a close correspondence between droughts and saltwater intrusion. Tritium/helium-3 apparent (piston-flow) ages determined from samples of saltwater with chloride concentrations of about 1,000 milligrams per liter (mg/L) or greater generally corresponded to a period during which droughts were frequent. Comparison of average daily air temperatures in Miami, Florida, with estimates of recharge temperatures determined from the dissolved gas composition in water samples indicated that saltwater likely entered the aquifer in April or early May when water levels are typically at their lowest during the year. Conversely, most of the samples of freshwater with chloride concentrations less than about 1,000 mg/L indicate recharge temperatures corresponding to air temperatures in mid to late May when rainfall and water levels in the aquifer increase, and the piston-flow ages of these samples correspond to wet years. The piston-flow ages of freshwater samples generally were younger than ages of samples of saltwater. Saltwater samples that were depleted in boron, magnesium, potassium, sodium, and sulfate, and enriched in calcium relative to the concentrations theoretically produced by freshwater/seawater mixing, generally were found to be associated with areas where saltwater had recently intruded. The calcium to (bicarbonate + sulfate) molar ratios (Ca/(HCO 3 +SO 4 )) of these samples generally were greater than 1. Saltwater samples from some of the monitoring wells, however, indicated little or no enrichment or depletion of these ions relative to the theoretical freshwater/seawater mixing line, and the Ca/(HCO 3 +SO 4 ) molar ratios of these samples generally were less than 1. Results indicated that aquifer materials are approaching equilibrium with seawater at these well locations.

Florida↗

Simulation of the effects of rainfall and groundwater use on historical lake water levels, groundwater levels, and spring flows in central Florida

The urbanization of central Florida has progressed substantially in recent decades, and the total population in Lake, Orange, Osceola, Polk, and Seminole Counties more than quadrupled from 1960 to 2010. The Floridan aquifer system is the primary source of water for potable, industrial, and agricultural purposes in central Florida. Despite increases in groundwater withdrawals to meet the demand of population growth, recharge derived by infiltration of rainfall in the well-drained karst terrain of central Florida is the largest component of the long-term water balance of the Floridan aquifer system. To complement existing physics-based groundwater flow models, artificial neural networks and other data-mining techniques were used to simulate historical lake water level, groundwater level, and spring flow at sites throughout the area. Historical data were examined using descriptive statistics, cluster analysis, and other exploratory analysis techniques to assess their suitability for more intensive data-mining analysis. Linear trend analyses of meteorological data collected by the National Oceanic and Atmospheric Administration at 21 sites indicate 67 percent of sites exhibited upward trends in air temperature over at least a 45-year period of record, whereas 76 percent exhibited downward trends in rainfall over at least a 95-year period of record. Likewise, linear trend analyses of hydrologic response data, which have varied periods of record ranging in length from 10 to 79 years, indicate that water levels in lakes (307 sites) were about evenly split between upward and downward trends, whereas water levels in 69 percent of wells (out of 455 sites) and flows in 68 percent of springs (out of 19 sites) exhibited downward trends. Total groundwater use in the study area increased from about 250 million gallons per day (Mgal/d) in 1958 to about 590 Mgal/d in 1980 and remained relatively stable from 1981 to 2008, with a minimum of 559 Mgal/d in 1994 and a maximum of 773 Mgal/d in 2000. The change in groundwater-use trend in the early 1980s and the following period of relatively slight trend is attributable to the concomitant effects of increasing public-supply withdrawals and decreasing use of water by the phosphate industry and agriculture. On the basis of available historical data and exploratory analyses, empirical lake water-level, groundwater-level, and spring-flow models were developed for 22 lakes, 23 wells, and 6 springs. Input time series consisting of various frequencies and frequency-band components of daily rainfall (1942 to 2008) and monthly total groundwater use (1957 to 2008) resulted in hybrid signal-decomposition artificial neural network models. The final models explained much of the variability in observed hydrologic data, with 43 of the 51 sites having coefficients of determination exceeding 0.6, and the models matched the magnitude of the observed data reasonably well, such that models for 32 of the 51 sites had root-mean-square errors less than 10 percent of the measured range of the data. The Central Florida Artificial Neural Network Decision Support System was developed to integrate historical databases and the 102 site-specific artificial neural network models, model controls, and model output into a spreadsheet application with a graphical user interface that allows the user to simulate scenarios of interest. Overall, the data-mining analyses indicate that the Floridan aquifer system in central Florida is a highly conductive, dynamic, open system that is strongly influenced by external forcing. The most important external forcing appears to be rainfall, which explains much of the multiyear cyclic variability and long-term downward trends observed in lake water levels, groundwater levels, and spring flows. For most sites, groundwater use explains less of the observed variability in water levels and flows than rainfall. Relative groundwater-use impacts are greater during droughts, however, and long-term trends in water levels and flows were identified that are consistent with historical groundwater-use patterns. The sensitivity of the hydrologic system to rainfall is expected, owing to the well-drained karst terrain and relatively thin confinement of the Floridan aquifer system in much of central Florida. These characteristics facilitate the relatively rapid transmission of infiltrating water from rainfall to the water table and contribute to downward leakage of water to the Floridan aquifer system. The areally distributed nature of rainfall, as opposed to the site-specific nature of groundwater use, and the generally high transmissivity and low storativity properties of the semiconfined Floridan aquifer system contribute to the prevalence of water-level and flow patterns that mimic rainfall patterns. In general, the data-mining analyses demonstrate that the hydrologic system in central Florida is affected by groundwater use differently during wet periods, when little or no system storage is available (high water levels), compared to dry periods, when there is excess system storage (low water levels). Thus, by driving the overall behavior of the system, rainfall indirectly influences the degree to which groundwater use will effect persistent trends in water levels and flows, with groundwater-use impacts more prevalent during periods of low water levels and spring flows caused by low rainfall and less prevalent during periods of high water levels and spring flows caused by high rainfall. Differences in the magnitudes of rainfall and groundwater use during wet and dry periods also are important determinants of hydrologic response. An important implication of the data-mining analyses is that rainfall variability at subannual to multidecadal timescales must be considered in combination with groundwater use to provide robust system-response predictions that enhance sustainable resource management in an open karst aquifer system. The data-driven approach was limited, however, by the confounding effects of correlation between rainfall and groundwater use, the quality and completeness of the historical databases, and the spatial variations in groundwater use. The data-mining analyses indicate that available historical data when used alone do not contain sufficient information to definitively quantify the related individual effects of rainfall and groundwater use on hydrologic response. The knowledge gained from data-driven modeling and the results from physics-based modeling, when compared and used in combination, can yield a more comprehensive assessment and a more robust understanding of the hydrologic system than either of the approaches used separately.

Florida↗

Estimation of unaltered daily mean streamflow at ungaged streams of New York, excluding Long Island, water years 1961-2010

The lakes, rivers, and streams of New York State provide an essential water resource for the State. The information provided by time series hydrologic data is essential to understanding ways to promote healthy instream ecology and to strengthen the scientific basis for sound water management decision making in New York. The U.S. Geological Survey, in cooperation with The Nature Conservancy and the New York State Energy Research and Development Authority, has developed the New York Streamflow Estimation Tool to estimate a daily mean hydrograph for the period from October 1, 1960, to September 30, 2010, at ungaged locations across the State. The New York Streamflow Estimation Tool produces a complete estimated daily mean time series from which daily flow statistics can be estimated. In addition, the New York Streamflow Estimation Tool provides a means for quantitative flow assessments at ungaged locations that can be used to address the objectives of the Clean Water Act&mdash;to restore and maintain the chemical, physical, and biological integrity of the Nation&rsquo;s waters. The New York Streamflow Estimation Tool uses data from the U.S. Geological Survey streamflow network for selected streamgages in New York (excluding Long Island) and surrounding States with shared hydrologic boundaries, and physical and climate basin characteristics to estimate the natural unaltered streamflow at ungaged stream locations. The unaltered streamflow is representative of flows that are minimally altered by regulation, diversion, or mining, and other anthropogenic activities. With the streamflow network data, flow-duration exceedance probability equations were developed to estimate unaltered streamflow exceedance probabilities at an ungaged location using a methodology that equates streamflow as a percentile from a flow-duration curve for a particular day at a hydrologically similar reference streamgage with streamflow as a percentile from the flow-duration curve for the same day at an ungaged location. The reference streamgage is selected using map correlation, a geostatistical method in which variogram models are developed that correlate streamflow at one streamgage with streamflows at all other locations in the study area. Regression equations used to predict 17 flow-duration exceedance probabilities were developed to estimate the flow-duration curves at ungaged locations for New York using geographic information system-derived basin characteristics. A graphical user interface, with an integrated spreadsheet summary report, has been developed to estimate and display the daily mean streamflows and statistics and to evaluate different water management or water withdrawal scenarios with the estimated monthly data. This package of regression equations, U.S. Geological Survey streamgage data, and spreadsheet application produces an interactive tool to estimate an unaltered daily streamflow hydrograph and streamflow statistics at ungaged sites in New York. Among other uses, the New York Streamflow Estimation Tool can assist water managers with permitting water withdrawals, implementing habitat protection, estimating contaminant loads, or determining the potential affect from chemical spills.

New York↗

Regional chloride distribution in the Northern Atlantic Coastal Plain aquifer system from Long Island, New York, to North Carolina

The aquifers of the Northern Atlantic Coastal Plain are the principal source of water supply for the region&rsquo;s nearly 20 million residents. Water quality and water levels in the aquifers, and maintenance of streamflow, are of concern because of the use of this natural resource for water supply and because of the possible effects of climate change and changes in land use on groundwater. The long-term sustainability of this natural resource is a concern at the local community scale, as well as at a regional scale, across state boundaries. In 2010, the U.S. Geological Survey (USGS) began a regional assessment of the Northern Atlantic Coastal Plain aquifers. An important part of this assessment is a regional interpretation of the extent of saltwater and the proximity of saltwater to fresh-groundwater resources and includes samples and published interpretations of chloride concentrations newly available since the last regional chloride assessment in 1989. This updated assessment also includes consideration of chloride samples and refined interpretations that stem from the 1994 discovery of the buried 35 million year old Chesapeake Bay impact structure that has substantially altered the understanding of the hydrogeologic framework and saltwater distribution in eastern Virginia. In this study, the regional area of concern for the chloride samples and interpretations extends from the Fall Line in the west to the outer edge of the Continental Shelf in the east and from the eastern tip of Long Island in the north to about halfway down the North Carolina coast in the south. Discussions of chloride distribution are presented for each of the 10 regional aquifer layers of the Northern Atlantic Coastal Plain, including the offshore extents. Maps of interpreted lines of equal concentration or isochlors were manually prepared for nine of the regional aquifers; a map was not prepared for the surficial regional aquifer. The isochlor interpretations include the offshore extent of the nine regional aquifers and are presented on a 1:2,000,000 scale base map. Vertically, the chloride samples and interpretations range from deepest (oldest) to shallowest (youngest)&mdash;Potomac-Patuxent, Potomac-Patapsco, Magothy, Matawan, Monmouth-Mount Laurel, Aquia, Piney Point, Lower Chesapeake, and Upper Chesapeake regional aquifers. The approach of this study maximizes the overall density of chloride information and data by assessing relevant published interpretations, all USGS chloride samples, and all relevant offshore samples in one comprehensive interpretation. Published isochlors, where they were interpreted by regional aquifer, were used as much as possible for this regional isochlor assessment. Publication dates for the isochlors used range from 1982 to 2015, and the scales for the isochlors range from local (county or municipality) to state (sub-regional) to regional. The USGS National Water Information System database provided well sample data for the parts of aquifers that are mainly beneath the land areas and yielded 37,517 water-quality records for 1903 through 2011. Published data reports from four phases of research-related offshore coring (1976, 1993, 1997, 2009) were the main source of water-quality data for the parts of aquifers from the shoreline to the outer edge of the Continental Shelf and yielded samples from multiple depths of each of 13 cores. This study also used interpretations and offshore core data from the last regional chloride assessment (1989) which, in addition to 7 offshore cores, included water-quality data from about 500 wells, and borehole geophysics interpretations from a subset of 11 wells. All published information and data that were used in this study were considered time independent and did not assess the published interpretations or data for temporal trends. The approach used here examined only published interpretations and available chloride data, and did not directly use supplemental techniques that can provide insight into the distribution of saltwater, such as geochemical characterization, borehole geophysical information, and geochronology. Isochlor maps for this study are limited to manual interpretations of the 250-milligram per liter (mg/L) and 10,000-mg/L boundaries developed for 9 of the 10 regional aquifers that constitute the regional hydrogeologic framework of the Northern Atlantic Coastal Plain. For a given aquifer, the approach was to initially consider published isochlor interpretations, where available, then to modify the published interpretations, if necessary, to the extent indicated by the well and core samples. The final step was to interpolate isochlors to the full extent of each aquifer layer in areas with sufficient samples or cited interpretations, or to extrapolate isochlors in areas with no samples or where samples were sparse. The principal limitation of this study is that, because of its regional extent, data and information density can vary greatly, and thus confidence in interpretations can vary widely for onshore and offshore areas across the study area. In areas of sparse data, some samples of elevated chloride could be misinterpreted as being part of a regional elevated chloride trend, and in other cases, an elevated concentration could be misinterpreted as being of only local importance. The interpretive work of this study was applied to a 1:2,000,000 scale base map. Locations of isochlors, wells, cores, political boundaries, and shorelines are meant to be considered approximate. The isochlors presented in this study were manually interpreted for each aquifer unit as a conceptual representation of an equal concentration line approximately in the middle of an aquifer&rsquo;s thickness. Differences in chloride concentration lines between the top and bottom of an aquifer could be substantial, especially for the thick parts of aquifers, but that information is not presented in this regional assessment. Although additional offshore chloride data are available compared to 27 years ago (1989), the offshore information remains sparse, resulting in less confidence in the offshore interpretations than in the onshore interpretations. Regionally, the 250- and 10,000-mg/L isochlors tend to map progressively eastward from the deepest to the shallowest aquifers across the Northern Atlantic Coastal Plain aquifer system but with some exceptions. The additional data, conceptual understanding, and interpretations in the vicinity of the buried Chesapeake Bay impact structure in eastern Virginia resulted in substantial refinement of isochlors in that area. Overall, the interpretations in this study are updates of the previous regional study from 1989 but do not comprise major differences in interpretation and do not indicate regional movement of the freshwater-saltwater interface since then.

Delaware, Maryland, New Jersey, New York, North Ca↗

UFINCH: A method for simulating unit and daily flows in networks of channels described by NHDPlus using continuous flow data at U.S. Geological Survey streamgages

The UFINCH (Unit Flows In Networks of Channels) computer application can be used to simulate daily and unit flows in networks of streams based on geospatial data in the National Hydrography Dataset NHDPlus (with value added attributes), and U.S. Geoogical Survey daily streamflow data from a downstream (or base) streamgage. Among streamflow augmentation methods, UFINCH has the unique capability to estimate time series of flows from a single base (downstream) streamgage to many upstream reaches, while conserving flows within the basin. UFINCH also provides a simple statistical model to adjust simulated flows to better match continuous flows from data at an upstream streamgage. Parameters of the statistical model are estimated using overlapping periods of record at the two streamgages, but the adjustment can be applied to all years of record available at the base streamgage. This report describes the main features of UFINCH and presents results from a sample application. Interactive graphical user interfaces and automated geographical information processing facilitate flow-data retrievals provide an intuitive environment for efficient and effective generation of flow information in a network. UFINCH is coded in the Matlab programming language and can be run in the Matlab programming environment, with supporting statistical, optimization, and mapping toolboxes, or from compiled code on a Microsoft Windows computer.

Scientific Investigations Report↗