USGS Science⌕ Search

SEARCH · USGS Science

Results for “Nature Scientific Reports”

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,621 records · Page 90Linked to original sources

Hydraulic, geomorphic, and trout habitat conditions of the Lake Fork of the Gunnison River in Hinsdale County, Lake City, Colorado, Water Years 2010-2011

Channel rehabilitation, or reconfiguration, to mitigate a variety of riverine problems has become a common practice in the western United States. However, additional work to monitor and assess the channel response to, and the effectiveness of, these modifications over longer periods of time (decadal or longer) is still needed. The Lake Fork of the Gunnison River has been an area of active channel modification to accommodate the needs of the Lake City community since the 1950s. The Lake Fork Valley Conservancy District began a planning process to assess restoration options for a reach of the Lake Fork in Lake City to enhance hydraulic and ecologic characteristics of the reach. Geomorphic channel form is affected by land-use changes within the basin and geologic controls within the reach. The historic channel was defined as a dynamic, braided channel with an active flood plain. This can result in a natural tendency for the channel to braid. A braided channel can affect channel stability of reconfigured reaches when a single-thread meandering channel is imposed on the stream. The U.S. Geological Survey, in cooperation with the Colorado Water Conservation Board and Colorado River Water Conservation District, began a study in 2010 to quantify existing hydraulic and habitat conditions for a reach of the Lake Fork of the Gunnison River in Lake City, Colorado. The purpose of this report is to quantify existing Lake Fork hydraulic and habitat conditions and establish a baseline against which post-reconfiguration conditions can be compared. This report (1) quantifies the existing hydraulic and geomorphic conditions in a 1.1-kilometer section of the Lake Fork at Lake City that has been proposed as a location for future channel-rehabilitation efforts, (2) characterizes the habitat suitability of the reach for two trout species based on physical conditions within the stream, and (3) characterizes the current riparian canopy density. The FaSTMECH computational flow-model within MD_SWMS was selected to characterize the effects of streamflow on hydraulic and habitat-suitability conditions for a study reach of the Lake Fork. Habitat suitability was evaluated for cutthroat ( Oncorhynchus clarkii ) and brown trout ( Salmo trutta morpha fario ) fry, juveniles, and adults. Microscale (point locations) and mesoscale (reach features) habitats were assessed using the combination of field observations, measurements, and hydraulic simulations within the study reach of the Lake Fork. Microscale trout habitat, presented as weighted usable area, generally increased as streamflow increased for both trout species and all life stages. Areas of suitable microscale habitat occur along the banks for flows of 900 cubic feet per second (ft 3 /s) and less. Out-of-bank areas became more substantial contributors to overall habitat availability for flows of 1,300 ft 3 /s or more when compared to other features. Adult habitat, for both trout species, was the most abundant habitat type for nearly all streamflows. In general, the upper reach provided 2–3 times more available habitat than the lower reach for both trout species. Mesoscale trout habitat of the Lake Fork was assessed based on the conditions present in the 150 ft 3 /s flow simulation as well as field observation. Both the upper and lower reach is primarily characterized as riffle/run habitat. The presence of pool habitat was limited throughout both reaches and occurred along the channel margins. For both reaches, the pool habitat was less than 5 percent of the total wetted area, a percentage that is substantially lower than the recommendations for sustainable populations of 40–70 percent. Areas of cover were adjacent to potential drift feeding areas in the lower reach, and often occurred within the same pool habitat. This may favor energy expenditure ratios of both fish species, wherein little energy is needed to acquire adequate food sources. Sediment mobility is an important process for flushing fine sediments from within the gravel frameworks. Evaluations of channel and flow characteristics at cross-section locations 2–8 show a range of streambed mobility. In general, boundary shear stress and streambed mobility increase with increases in streamflow. Within the cross sections, the greatest boundary shear stress occurs towards the center of the channel. Reach-scale assessment of sediment mobility in the lower reach shows increased streambed mobility. This is due in part to smaller grain sizes in the lower reach, but may also reflect the greater extent of channel alterations, specifically the temporary berms constructed by CDOT in the late 1980s and 1990s, present in this reach.

Colorado↗

Simulations of Ground-Water Flow, Transport, Age, and Particle Tracking near York, Nebraska, for a Study of Transport of Anthropogenic and Natural Contaminants (TANC) to Public-Supply Wells

Contamination of public-supply wells has resulted in public-health threats and negative economic effects for communities that must treat contaminated water or find alternative water supplies. To investigate factors controlling vulnerability of public-supply wells to anthropogenic and natural contaminants using consistent and systematic data collected in a variety of principal aquifer settings in the United States, a study of Transport of Anthropogenic and Natural Contaminants to public-supply wells was begun in 2001 as part of the U.S. Geological Survey National Water-Quality Assessment Program. The area simulated by the ground-water flow model described in this report was selected for a study of processes influencing contaminant distribution and transport along the direction of ground-water flow towards a public-supply well in southeastern York, Nebraska. Ground-water flow is simulated for a 60-year period from September 1, 1944, to August 31, 2004. Steady-state conditions are simulated prior to September 1, 1944, and represent conditions prior to use of ground water for irrigation. Irrigation, municipal, and industrial wells were simulated using the Multi-Node Well package of the modular three-dimensional ground-water flow model code, MODFLOW-2000, which allows simulation of flow and solutes through wells that are simulated in multiple nodes or layers. Ground-water flow, age, and transport of selected tracers were simulated using the Ground-Water Transport process of MODFLOW-2000. Simulated ground-water age was compared to interpreted ground-water age in six monitoring wells in the unconfined aquifer. The tracer chlorofluorocarbon-11 was simulated directly using Ground-Water Transport for comparison with concentrations measured in six monitoring wells and one public supply well screened in the upper confined aquifer. Three alternative model simulations indicate that simulation results are highly sensitive to the distribution of multilayer well bores where leakage can occur and that the calibrated model resulted in smaller differences than the alternative models between simulated and interpreted ages and measured tracer concentrations in most, but not all, wells. Results of the first alternative model indicate that the distribution of young water in the upper confined aquifer is substantially different when well-bore leakage at known abandoned wells and test holes is removed from the model. In the second alternative model, simulated age near the bottom of the unconfined aquifer was younger than interpreted ages and simulated chlorofluorocarbon-11 concentrations in the upper confined aquifer were zero in five out of six wells because the conventional Well Package fails to account for flow between model layers though well bores. The third alternative model produced differences between simulated and interpreted ground-water ages and measured chlorofluorocarbon-11 concentrations that were comparable to the calibrated model. However, simulated hydraulic heads deviated from measured hydraulic heads by a greater amount than for the calibrated model. Even so, because the third alternative model simulates steady-state flow, additional analysis was possible using steady-state particle tracking to assess the contributing recharge area to a public supply well selected for analysis of factors contributing to well vulnerability. Results from particle-tracking software (MODPATH) using the third alternative model indicates that the contributing recharge area of the study public-supply well is a composite of elongated, seemingly isolated areas associated with wells that are screened in multiple aquifers. The simulated age distribution of particles at the study public-supply well indicates that all water younger than 58 years travels through well bores of wells screened in multiple aquifers. The age distribution from the steady-state model using MODPATH estimates the youngest 7 percent of the water to have a flow-weighted mean age

Scientific Investigations Report↗

Physical habitat, water quality, and riverine biological assemblages of selected reaches of the Sheyenne River, North Dakota, 2010

In 2010, data on physical habitat, water quality, and riverine biological assemblages were collected at selected reaches in four locations (Kleven, Sheyenne, Cooperstown, and West Fargo) on the Sheyenne River in east-central North Dakota. Three of the locations (Kleven, Sheyenne, and Cooperstown) are above Baldhill Dam and one location (West Fargo) is below Baldhill Dam on the Sheyenne River. The 2010 data provide information to establish a better understanding of the water-quality and ecological conditions of the Sheyenne River. Concerns were raised about the water-quality and ecological conditions of the Sheyenne River because of the interbasin transfer of water from nearby Devils Lake. The transfer of water from Devils Lake to the Sheyenne River occurs through the Devils Lake State Outlet near Peterson Coulee or, if lake elevations exceed 1,459 feet above National Geodetic Vertical Datum of 1929 (NGVD 29), through a natural outlet, Tolna Coulee. The field measurements of water-quality characteristics and results of chemical analyses generally are comparable to summary statistics calculated for Sheyenne River for 1980 through 2006. Overall, water-quality results show differences between the Kleven, Sheyenne, Cooperstown, and West Fargo reaches. Sulfate concentrations were less than the State of North Dakota criterion of 750 milligrams per liter for the upper Sheyenne River above Baldhill Dam and less than the criterion of 450 milligrams per liter for the lower Sheyenne River below Baldhill Dam. Arsenic concentrations at most reaches exceeded the U.S. Environmental Protection Agency drinking-water standard of 10 micrograms per liter. Nutrient concentrations (nitrogen, phosphorus) were higher in the upper Sheyenne River above Baldhill Dam than below Baldhill Dam where concentrations decreased by about half. In 2010, 35 families and 44 genera of benthic macroinvertebrates were collected and identified. On the basis of the index of biotic intergrity scores for benthic macroinvertebrate communities present in the Sheyenne River, all the reaches were determined to have condition classes of moderately disturbed to most disturbed. The benthic macroinvertebrate communities at the Cooperstown reaches were classed as moderately disturbed, whereas benthic macroinvertebrate communities at the Kleven, Sheyenne, West and Fargo reaches were most disturbed. During data collection, 37 genera and 165 species of periphyton (diatoms and soft-bodied algae) were collected and identified. In periphyton communities, similar taxa species were dominant in the Kleven, Sheyenne, and Cooperstown reaches, and different taxa species were dominant in the West Fargo reaches. For diatoms, the Kleven 3 reach had the lowest species richness value of 33.0, whereas the Cooperstown 8 reach had the highest species richness value of 57.0. For soft-bodied algae, the species richness values ranged from 8.0 at the Sheyenne 4 reach to 20.0 at the West Fargo 10 reach. During the fish collection, 32 species, representing 10 families, were collected in the Sheyenne River. All but two species are native to the Sheyenne River system. Common carp and white crappie are the two introduced species. Of the 32 species, 29 are tolerant to moderately tolerant to changes in water quality and habitat degradation, 16 species are tolerant to moderately tolerant to turbidity, and 16 species are tolerant to moderately tolerant to sensitivity to total dissolved solids, sulfate, and chloride. All fish species were categorized into four trophic groups. The largest group of 19 species was the insectivores (both benthic and general). The predator group consisted of seven species, and the omnivores consisted of six species. More fish were found in the lower Sheyenne River below Baldhill Dam than in the upper Sheyenne River above Baldhill Dam.

North Dakota↗

Numerical simulation of flow in deep open boreholes in a coastal freshwater lens, Pearl Harbor Aquifer, O‘ahu, Hawai‘i

The Pearl Harbor aquifer in southern O‘ahu is one of the most important sources of freshwater in Hawai‘i. A thick freshwater lens overlays brackish and saltwater in this coastal aquifer. Salinity profiles collected from uncased deep monitor wells (DMWs) commonly are used to monitor freshwater-lens thickness. However, vertical flow in DMWs can cause the measured salinity to differ from salinity in the adjacent aquifer or in an aquifer without a DWM. Substantial borehole flow and displacement of salinity in DMWs over several hundred feet have been observed in the Pearl Harbor aquifer. The objective of this study was to evaluate the effects of borehole flow on measured salinity profiles from DMWs. A numerical modeling approach incorporated aquifer hydraulic characteristics and recharge and withdrawal rates representative of the Pearl Harbor aquifer. Borehole flow caused by vertical hydraulic gradients associated with both the natural regional flow system and groundwater withdrawals was simulated. Model results indicate that, with all other factors being equal, greater withdrawal rates, closer withdrawal locations, or higher hydraulic conductivities of the well cause greater borehole flow and displacement of salinity in the well. Borehole flow caused by the natural groundwater-flow system is five orders of magnitude greater than vertical flow in a homogeneous aquifer, and borehole-flow directions are consistent with the regional flow system: downward flow in inland recharge areas and upward flow in coastal discharge areas. Displacement of salinity inside the DMWs associated with the regional groundwater-flow system ranges from less than 1 to 220 ft, depending on the location and assumed hydraulic conductivity of the well. For example, upward displacements of the 2 percent and 50 percent salinity depths in a well in the coastal discharge part of the flow system are 17 and 4.4 ft, respectively, and the average salinity difference between aquifer and borehole is 0.65 percent seawater salinity. Groundwater withdrawals and drawdowns generally occur at shallow depths in the freshwater system with respect to the depth of the DMW and cause upward flow in the DMW. Simulated groundwater withdrawal of 4.3 million gallons per day that is 100 ft from a DMW causes thirty times more borehole flow than borehole flow that is induced by the regional flow field alone. The displacement of the 2 percent borehole salinity depth increases from 17 to 33 ft, and the average salinity difference between aquifer and borehole is 0.85 percent seawater salinity. Peak borehole flow caused by local groundwater withdrawal near DMWs is directly proportional to the pumping rate in the nearby production well. Increasing groundwater withdrawal to 16.7 million gallons per day increases upward displacement of the 50 percent salinity depth (midpoint of the transition zone) from 4.6 to 77 ft, and the average salinity difference between aquifer and borehole is 1.4 percent seawater salinity. Simulated groundwater withdrawal that is 3,000 ft away from DMWs causes less borehole flow and salinity displacements than nearby withdrawal. Simulated effects of groundwater withdrawal from a horizontal shaft and withdrawal from a vertical well in a homogeneous aquifer were similar. Generally, the 50 percent salinity depths are less affected by borehole flow than the 2 percent salinity depths. Hence, measured salinity profiles are useful for calibration of regional numerical models despite borehole-flow effects. Commonly, a 1 percent error in salinity is acceptable in numerical modeling studies. Incorporation of heterogeneity in the model is necessary to simulate long vertical steps observed in salinity profiles in southern O‘ahu. A thick zone of low aquifer hydraulic conductivity limits exchange of water between aquifer and well and creates a long vertical step in the salinity profile. A heterogeneous basalt-aquifer scenario simulates observed vertical salinity steps and borehole flow that is consistent with measured borehole flow from DMWs in southern O‘ahu. However, inclusion of local-scale heterogeneities in regional models generally is not warranted.

Hawai'i↗

Travel Times, Streamflow Velocities, and Dispersion Rates in the Yellowstone River, Montana

The Yellowstone River is a vital natural resource to the residents of southeastern Montana and is a primary source of water for irrigation and recreation and the primary source of municipal water for several cities. The Yellowstone River valley is the primary east-west transportation corridor through southern Montana. This complex of infrastructure makes the Yellowstone River especially vulnerable to accidental spills from various sources such as tanker cars and trucks. In 2008, the U.S. Geological Survey (USGS), in cooperation with the Montana Department of Environmental Quality, initiated a dye-tracer study to determine instream travel times, streamflow velocities, and dispersion rates for the Yellowstone River from Lockwood to Glendive, Montana. The purpose of this report is to describe the results of this study and summarize data collected at each of the measurement sites between Lockwood and Glendive. This report also compares the results of this study to estimated travel times from a transport model developed by the USGS for a previous study. For this study, Rhodamine WT dye was injected at four locations in late September and early October 2008 during reasonably steady streamflow conditions. Streamflows ranged from 3,490 to 3,770 cubic feet per second upstream from the confluence of the Bighorn River and ranged from 6,520 to 7,570 cubic feet per second downstream from the confluence of the Bighorn River. Mean velocities were calculated for each subreach between measurement sites for the leading edge, peak concentration, centroid, and trailing edge at 10 percent of the peak concentration. Calculated velocities for the centroid of the dye plume for subreaches that were completely laterally mixed ranged from 1.83 to 3.18 ft/s within the study reach from Lockwood Bridge to Glendive Bridge. The mean of the completely mixed centroid velocity for the entire study reach, excluding the subreach between Forsyth Bridge and Cartersville Dam, was 2.80 ft/s. Longitudinal dispersion rates of the dye plume for this study ranged from 0.06 ft/s for the subreach upstream from Forsyth Bridge to 2.25 ft/s for the subreach upstream from Calyspo Bridge for subreaches where the dye was completely laterally mixed. A relation was determined between travel time of the peak concentration and time for the dye plume to pass a site (duration). This relation can be used to estimate when the receding concentration of a potential contaminant reaches 10 percent of its peak concentration for accidental spills into the Yellowstone River. Data from this dye-tracer study were used to evaluate velocity and concentration estimates from a transport model developed as part of an earlier USGS study. Comparison of the estimated and calculated velocities for the study reach indicate that the transport model estimates the velocities of the Yellowstone River between Huntley Bridge and Glendive Bridge with reasonable accuracy. Velocities of the peak concentration of the dye plume calculated for this study averaged 10 percent faster than the most probable velocities and averaged 12 percent slower than the maximum probable velocities estimated from the transport model. Peak Rhodamine WT dye concentrations were consistently lower than the transport model estimates except for the most upstream subreach of each dye injection. The most upstream subreach of each dye injection is expected to have a higher concentration because of incomplete lateral mixing. Lower measured peak concentrations for all other sites were expected because Rhodamine WT dye deteriorates when exposed to sunlight and will sorb onto the streambanks and stream bottom. Velocity-streamflow relations developed by using routine streamflow measurements at USGS gaging stations and the transport model can be used to estimate mean streamflow velocities throughout a range of streamflows. The variation in these velocity-streamflow relations emphasizes the uncertainty in estimating the mean streamflow veloc

Scientific Investigations Report↗

Occurrence of trihalomethanes in the nation's ground water and drinking-water supply wells, 1985-2002

This report describes the occurrence of trihalomethanes (THMs) in the Nation's ground water and drinking-water supply wells based on analysis of 5,642 samples of untreated ground water and source water collected or compiled during 1985-2002 by the U.S. Geological Survey National Water-Quality Assessment (NAWQA) Program. THMs are a group of volatile organic compounds (VOCs) with natural and anthropogenic sources that are of interest because they are associated with acute and chronic health problems in humans. THMs occur in water primarily from chlorination and are classified as disinfection by-products. In this report, the four THMs are discussed in the order of chloroform, bromodichloromethane, dibromochloromethane, and then bromoform; this sequence corresponds to largest to smallest chlorine content and smallest to largest bromine content. Four trihalomethanes were detected in less than 20 percent of samples from studies of (1) aquifers, (2) shallow ground water in agricultural areas, (3) shallow ground water in urban areas, (4) domestic wells, and (5) public wells. Detection frequencies for individual THMs in the five studies ranged from zero for shallow ground water in agricultural areas to 19.5 percent for shallow ground water in urban areas. None of the samples from aquifer studies, domestic wells, or public wells had total THM concentrations (the sum of the concentrations of chloroform, bromodichloromethane, dibromochloromethane, and bromoform) greater than or equal to the U.S. Environmental Protection Agency Maximum Contaminant Level of 80 micrograms per liter (?g/L). Comparisons of results among studies of aquifers, shallow ground water in agricultural areas, and shallow ground water in urban areas were used to describe the occurrence of the four THMs in ground water for three different land-use settings-mixed, agricultural, and urban, respectively. At the 0.2-?g/L assessment level, one or more of the four THMs were detected in 7.9 percent of the samples from aquifer studies, 2.2 percent of the samples from shallow ground water in agricultural areas, and 19.5 percent of the samples from shallow ground water in urban areas. In general, detection frequencies and concentrations of the four THMs were greater in shallow ground water in urban areas compared to aquifer studies and to shallow ground water in agricultural areas. For all three of these studies, the most common two-THM mixture at the 0.2-?g/L assessment level was chloroform-bromodichloromethane, and this was the only two-THM mixture found in samples of shallow ground water in agricultural areas. Comparisons of results between studies of domestic wells and public wells were used to describe the occurrence of the four THMs in two different supplies of ground water used for drinking water. At the 0.2-?g/L assessment level, one or more of the four THMs were detected in 5.2 percent of the domestic well samples and in 14.7 percent of the public well samples. In general, detection frequencies and THM concentrations were greater in samples from public wells than from domestic wells. At the 0.2-?g/L assessment level, the six possible two-THM mixtures occurred about six times more frequently in samples from public wells than from domestic wells. One of the most common two-THM mixtures in samples from domestic and public wells was bromodichloromethane-dibromochloromethane. Detection frequency is associated with the chlorine content of the THM compound. In general, for each of the five studies, as the chlorine content of the THM compound decreased, the detection frequency at the 0.2-?g/L assessment level also decreased. The exception was the study of public wells in which the detection frequency of the THMs decreased in the following order: chloroform, bromoform, dibromochloromethane, and bromodichloromethane. At the 0.2-?g/L assessment level, the median concentration for one or more of the four THMs ranged from 0.3 ?g/L (shallow ground water in agricultural a

Scientific Investigations Report↗

Modeling hydrodynamics, water temperature, and water quality in the Klamath River upstream of Keno Dam, Oregon, 2006-09

A hydrodynamic, water temperature, and water-quality model was constructed for a 20-mile reach of the Klamath River downstream of Upper Klamath Lake, from Link River to Keno Dam, for calendar years 2006-09. The two-dimensional, laterally averaged model CE-QUAL-W2 was used to simulate water velocity, ice cover, water temperature, specific conductance, dissolved and suspended solids, dissolved oxygen, total nitrogen, ammonia, nitrate, total phosphorus, orthophosphate, dissolved and particulate organic matter, and three algal groups. The Link-Keno model successfully simulated the most important spatial and temporal patterns in the measured data for this 4-year time period. The model calibration process provided critical insights into water-quality processes and the nature of those inputs and processes that drive water quality in this reach. The model was used not only to reproduce and better understand water-quality conditions that occurred in 2006-09, but also to test several load-reduction scenarios that have implications for future water-resources management in the river basin. The model construction and calibration process provided results concerning water quality and transport in the Link-Keno reach of the Klamath River, ranging from interesting circulation patterns in the Lake Ewauna area to the nature and importance of organic matter and algae. These insights and results include: * Modeled segment-average water velocities ranged from near 0.0 to 3.0 ft/s in 2006 through 2009. Travel time through the model reach was about 4 days at 2,000 ft 3 /s and 12 days at 700 ft3/s flow. Flow direction was aligned with the upstream-downstream channel axis for most of the Link-Keno reach, except for Lake Ewauna. Wind effects were pronounced at Lake Ewauna during low-flow conditions, often with circulation in the form of a gyre that rotated in a clockwise direction when winds were towards the southeast and in a counterclockwise direction when winds were towards the northwest. * Water temperatures ranged from near freezing in winter to near 30 degrees C at some locations and periods in summer; seasonal water temperature patterns were similar at the inflow and outflow. Although vertical temperature stratification was not present at most times and locations, weak stratification could persist for periods up to 1-2 weeks, especially in the downstream parts of the reach. Thermal stratification was important in controlling vertical variations in water quality. * The specific conductance, and thus density, of tributaries within the reach usually was higher than that of the river itself, so that inflows tended to sink below the river surface. This was especially notable for inflows from the Klamath Straits Drain, which tended to sink to the bottom of the Klamath River at its confluence and not mix vertically for several miles downstream. * The model was able to capture most of the seasonal changes in the algal population by modeling that population with three algal groups: blue-green algae, diatoms, and other algae. The blooms of blue-green algae, consisting mostly of Aphanizomenon flos aquae that entered from Upper Klamath Lake, were dominant, dwarfing the populations of the other two algae groups in summer. A large part of the blue-green algae population that entered this reach from upstream tended to settle out, die, and decompose, especially in the upper part of the Link-Keno reach. Diatoms reached a maximum in spring and other algae in midsummer. * Organic matter, occurring in both dissolved and particulate forms, was critical to the water quality of this reach of the Klamath River, and was strongly tied to nutrient and dissolved-oxygen dynamics. Dissolved and particulate organic matter were subdivided into labile (quickly decaying) and refractory (slowing decaying) groups for modeling purposes. The particulate matter in summer, consisting largely of dead blue-green algae, decayed quickly. Consequently, this particulate matt

Oregon↗

Historical changes to channel planform and bed elevations downstream from dams along Fall Creek and Middle Fork Willamette River, Oregon, 1926–2016

Operation of large, multipurpose dams within the Middle Fork Willamette River Basin, Oregon, including the Fall Creek sub-basin, have disrupted natural streamflow and sediment transport regimes and fish passage along the river corridors. Documenting channel morphology, including channel planform, landforms, vegetation cover, and river channel elevations at multiple points in time spanning the 20th and early 21st centuries, is useful for characterizing net changes occurring in response to construction and operation of these dams. The U.S. Geological Survey assessed historical channel changes that occurred within the past century in response to the construction and operation of flood-control dams by evaluating planimetric datasets (from 1926 plan and profile surveys and 1936 and 2016 aerial photographs) and elevation datasets (from 1926 plan and profile surveys and 2015 light detection and ranging [lidar]). This study specifically focuses on the lower 27.3 kilometers (km) of the Middle Fork Willamette River and the lower 11.5 km of Fall Creek, or the reaches downstream from the U.S. Army Corps of Engineers Dexter Dam and Fall Creek Dam, to the confluence with Coast Fork Willamette River. Altogether, compilation and evaluation of datasets for Fall Creek and the Middle Fork Willamette River downstream from the dams provide a foundation for understanding: channel morphology and patterns of geomorphic stability prior to dam construction in 1926 and 1936; channel morphology and patterns of lateral and vertical stability of the early 21st century that reflect present-day (post-dam) streamflow and sediment regimes as of 2015–16; and geomorphic transformations of the river corridors in the decades following dam construction, including changes in planform and bed elevation (determined from water-surface elevations). Findings from this study can be used to provide historical and geomorphic context for geomorphic responses to deep reservoir drawdowns on Fall Creek Lake that mobilize reservoir sediment downstream and informing other restoration and river-management activities; this report summarizes one component of a larger research effort to document the magnitude and spatial distribution of geomorphic responses to sediment releases from draining Fall Creek Lake. As of 2016, the modern Fall Creek flows through a narrow, semi-alluvial channel that efficiently conveys water and sediment at typical streamflows downstream from Fall Creek Dam. This channel planform, including the positions and distributions of bars and secondary water features (side channels, alcoves, and ponds), generally reflects pre-dam conditions in 1936, suggesting relatively modest morphological adjustments resulted from reductions in sediment supply and alterations to peak streamflow after dam construction. The most substantial morphologic change detected over this period was a reduction in unvegetated gravel bars. As of 2016, the modern Middle Fork Willamette River is a large, gravel-bed river that, despite substantial transformations in channel morphology and reduction in lateral dynamism following the construction of multiple upstream dams, remains a dominantly alluvial river. Prior to dam construction in 1926 and 1936, the reaches of the Middle Fork Willamette River downstream from Dexter Dam were laterally active with multi-thread and single-thread channels flanked by large, shifting gravel bars. Since streamflow regulation and other channel modifications in the mid-20th century, these reaches have become less laterally active and encompass a narrower floodplain corridor as abundant former gravel bars were converted to low-elevation floodplains colonized by young, dense forests. The Middle Fork Willamette River downstream from Dexter Dam has remained mostly vertically stable between 1926 and 2015, although localized segments possibly decreased in elevation as much as 2.3 meters.

Oregon↗

Effects of a drawdown on plant communities in a freshwater impoundment at Lacassine National Wildlife Refuge, Louisiana

Disturbance is an important natural process in the creation and maintenance of wetlands. Water depth manipulation and prescribed fire are two types of disturbance commonly used by humans to influence vegetation succession and composition in wetlands with the intention of improving wildlife habitat value. A 6,475-hectare (ha) impoundment was constructed in 1943 on Lacassine National Wildlife Refuge in southwest Louisiana to create freshwater wetlands as wintering waterfowl habitat. Ten years after construction of the impoundment, called Lacassine pool, was completed, refuge staff began expressing concerns about increasing emergent vegetation cover, organic matter accumulation, and decreasing area of open water within the pool. Because the presence of permanent standing water impedes actions that can address these concerns, a small impoundment within the pool where it was possible to manipulate water depth was created. The 283-ha subimpoundment called Unit D was constructed in 1989. Water was pumped from Unit D in 1990, and the unit was permanently reflooded about 3 years later. Four prescribed fires were applied during the drawdown. A study was initiated in 1990 to investigate the effect of the experimental drawdown on vegetation and soils in Unit D. Four plant community types were described, and cores were collected to measure the depth of the soil organic layer. A second study of Unit D was conducted in 1997, 4 years after the unit was reflooded, by using the same plots and similar sampling methods. This report presents an analysis and synthesis of the data from the two studies and provides an evaluation of the impact of the management techniques applied. We found that plant community characteristics often differed among the four communities and varied with time. Species richness increased in two of the communities, and total aboveground biomass increased in all four during the drawdown. These changes, however, did not persist when Unit D was reflooded; by 1997, species richness and aboveground biomass were equivalent to values before the drawdown. The change in waterfowl food value of the plant communities during the drawdown varied; it did not change in two communities, increased in one, and decreased in one. A consistent pattern noted was that waterfowl food value was higher in communities that contained open water than in those dominated by emergent plants, both soon after the drawdown was initiated in Unit D and 4 years after reflooding. A reduction in depth of the soil organic layer became apparent 20 months after drawdown was initiated, and this reduction persisted in 1997, 4 years after reflooding. A separate 2003 study on soil characteristics in Lacassine pool found that the depth to the clay layer was lower in Unit D than in the rest of the pool. We were not able to establish a cause-and-effect relation between any changes noted and the fact water levels in the unit were drawn down because the initial study in 1990 did not include control plots. Changes in vegetation and soil organic layer depth identified in Unit D may have occurred in the surrounding Lacassine pool habitat as well. Similarly, we were unable to form any conclusions about the effect of the prescribed fire treatments because there was no information on which plots were burned. Because of the known relation between anaerobic soil conditions and reduced decomposition of organic matter, however, it is likely that the drawdown in Unit D resulted in an increased decomposition rate and a reduction in the depth of the soil organic layer.

Louisiana↗

Geologic map of the Agnesi quadrangle (V-45), Venus

The Agnesi quadrangle (V–45), named for centrally located Agnesi crater, encompasses approximately 6,500,000 km 2 extending from lat 25° to 50° S. and from long 30° to 60° E. The V–45 quadrangle lies within Venus’ lowland broadly between highlands Ovda Regio to the northeast and Alpha Regio to the west. The region ranges in altitude from 6,051 to 6,054 km, with an average of ~6,052 km, which is essentially mean planetary radius. The quadrangle displays a wide range of features including large to small arcuate exposures of ribbon-tessera terrain (Hansen and Willis, 1998), ten lowland coronae, two montes, 13 pristine impact craters, and long but localized volcanic flows sourced to the west in V–44. Shield terrain (Hansen, 2005) occurs across much of the V–45 quadrangle. Although V–45 lies topographically within the lowland, it includes only one planitia (Fonueha Planitia), perhaps because the features mentioned decorate it. Geologic mapping of the Agnesi quadrangle (V–45) provides an opportunity (1) to examine the nature of lowland ribbon-tessera terrain, as compared to highland ribbon-tessera terrain; (2) to examine the nature and history of lowland coronae and montes to evaluate hypotheses for the evolution of these features; and (3) to evaluate global catastrophic/episodic resurfacing hypotheses. Ribbon-tessera terrain (Hansen and Willis, 1998), a structurally distinctive terrain marked by roughly orthogonal fold axes and ribbon troughs, characterizes Venusian crustal plateaus (Hansen and others, 1999), although this distinctive fabric is also variably preserved across some lowland regions (Hansen and Willis, 1996, 1998), including the Agnesi quadrangle. Isolated kipukas of ribbon-tessera terrain protrude through the shield-terrain veneer and preserve evidence of local surface processes that predated shield-terrain emplacement, including those leading to regional development of tessera terrain fabrics. The ribbon-tessera fabrics are similar to the fabrics that characterize high-standing crustal plateaus, such as nearby Ovda and Alpha Regiones. Ribbon-bearing kipukas may preserve evidence of ancient collapsed crustal plateaus (Phillips and Hansen, 1994; Ivanov and Head, 1996; Hansen and Willis, 1998; Ghent and Tibuleac, 2002), or they may record different, but rheologically similar, processes. Geologic mapping of the V–45 quadrangle constrains the geologic history of 12 lowland coronae and montes. Lowland coronae are relatively rare (10% of coronae); they do not form chains, typical of mesoland coronae (Phillips and Hansen, 1994), nor do they cluster, typical of highland coronae associated with volcanic rises (Stofan and others, 2001). The evolution of all coronae, whether by a single mechanism or a range of mechanisms, is a topic of debate (for example, Stofan and others, 1992, 2001; Vita-Finzi and others, 2005; Hamilton, 2005). The V–45 quadrangle hosts about 20 percent of the lowland coronae on Venus; the geologic history of these features should provide critical insight toward understanding coronae evolution. Two general classes of hypotheses have emerged to address the near random spatial distribution of ~970 apparently pristine impact craters across the surface of Venus: (1) catastrophic/episodic resurfacing and (2) equilibrium/evolutionary resurfacing. Catastrophic/episodic hypotheses propose that a global-scale, temporally punctuated event or events dominated Venus’ evolution and that the generally uniform impact crater distribution (Schaber and others, 1992; Phillips and others, 1992; Herrick and others, 1997) reflects craters that accumulated during relative global quiescence since that event (for example, Strom and others, 1994; Herrick, 1994; Turcotte and others, 1999). Equilibrium/evolutionary hypotheses suggest instead that the near random crater distribution results from relatively continuous, but spatially localized, resurfacing in which volcanic and (or) tectonic processes occur across the planet through time, although the style of operative processes may have varied temporally and spatially (for example, Phillips and others, 1992; Guest and Stofan, 1999; Hansen and Young, 2007). Geologic relations within the map area allow us to test the catastrophic/episodic versus equilibrium/evolutionary resurfacing hypotheses.

Scientific Investigations Map↗

Apparent Resistivity and Estimated Interaction Potential of Surface Water and Groundwater along Selected Canals and Streams in the Elkhorn-Loup Model Study Area, North-Central Nebraska, 2006-07

In 2005, the State of Nebraska adopted new legislation that in part requires local Natural Resources Districts to include the effect of groundwater use on surface-water systems in their groundwater management plan. In response the U.S. Geological Survey, in cooperation with the Upper Elkhorn, Lower Elkhorn, Upper Loup, Lower Loup, Middle Niobrara, Lower Niobrara, Lewis and Clark, and Lower Platte North Natural Resources Districts, did a study during 2006-07 to investigate the surface-water and groundwater interaction within a 79,800-square-kilometer area in north-central Nebraska. To determine how streambed materials affect surface-water and groundwater interaction, surface geophysical and lithologic data were integrated at four sites to characterize the hydrogeologic conditions within the study area. Frequency-domain electromagnetic and waterborne direct- current resistivity profiles were collected to map the near-surface hydrogeologic conditions along sections of Ainsworth Canal near Ainsworth, Nebraska; Mirdan and Geranium Canals near Ord, Nebraska; North Loup River near Ord, Nebraska; and Middle Loup River near Thedford, Nebraska. Lithologic data were collected from test holes at each site to aid interpretation of the geophysical data. Geostatistical analysis incorporating the spatial variability of resistivity was used to account for the effect of lithologic heterogeneity on effective hydraulic permeability. The geostatistical analysis and lithologic data descriptions were used to make an interpretation of the hydrogeologic system and derive estimates of surface-water/groundwater interaction potential within the canals and streambeds. The estimated interaction potential at the Ainsworth Canal site and the Mirdan and Geranium Canal site is generally low to moderately low. The sediment textures at nearby test holes typically were silt and clay and fine-to-medium sand. The apparent resistivity values for these sites ranged from 2 to 120 ohm-meters. The vertical and horizontal variability of the apparent resistivity data were consistently low. Low resistive variability indicates little lithologic heterogeneity for either canal site. The surface-water/groundwater interaction-potential estimates are in agreement with the narrow frequency distribution of resistivity, low apparent resistivities, low spatial heterogeneity, and test-hole grain-size ranges. The estimated surface-water/groundwater interaction potential at the North Loup and Middle Loup River sites is moderate to moderately high. The sediment textures at nearby test holes were predominantly fine, medium, and coarse sand with some silt and silty to sandy clay. The apparent resistivity values for these sites ranged from 34 to 1,338 ohm-meters. The vertical variability of the resistivity data was moderately high. The horizontal variability at these sites is low to moderately low. The higher resistive variability at these sites indicates generally greater lithologic heterogeneity than at either the Ainsworth Canal site or the Mirdan and Geranium Canal site. The surface-water/groundwater interaction-potential estimates are in agreement with the generally moderate to high apparent resistivity, the greater spatial heterogeneity, and the variable lithologic texture. A higher interaction potential as compared to the canal sites is expected because of the higher subsurface resistivity and greater lithologic heterogeneity.

Scientific Investigations Report↗

Determination of dilution factors for discharge of aluminum-containing wastes by public water-supply treatment facilities into lakes and reservoirs in Massachusetts

Dilution of aluminum discharged to reservoirs in filter-backwash effluents at water-treatment facilities in Massachusetts was investigated by a field study and computer simulation. Determination of dilution is needed so that permits for discharge ensure compliance with water-quality standards for aquatic life. The U.S. Environmental Protection Agency chronic standard for aluminum, 87 micrograms per liter (μg/L), rather than the acute standard, 750 μg/L, was used in this investigation because the time scales of chronic exposure (days) more nearly match rates of change in reservoir concentrations than do the time scales of acute exposure (hours). Whereas dilution factors are routinely computed for effluents discharged to streams solely on the basis of flow of the effluent and flow of the receiving stream, dilution determination for effluents discharged to reservoirs is more complex because (1), compared to streams, additional water is available for dilution in reservoirs during low flows as a result of reservoir flushing and storage during higher flows, and (2) aluminum removal in reservoirs occurs by aluminum sedimentation during the residence time of water in the reservoir. Possible resuspension of settled aluminum was not considered in this investigation. An additional concern for setting discharge standards is the substantial concentration of aluminum that can be naturally present in ambient surface waters, usually in association with dissolved organic carbon (DOC), which can bind aluminum and keep it in solution. A method for dilution determination was developed using a mass-balance equation for aluminum and considering sources of aluminum from groundwater, surface water, and filter-backwash effluents and losses caused by sedimentation, water withdrawal, and spill discharge from the reservoir. The method was applied to 13 reservoirs. Data on aluminum and DOC concentrations in reservoirs and influent water were collected during the fall of 2009. Complete reservoir volume was determined to be available for mixing on the basis of vertical and horizontal aluminum-concentration profiling. Losses caused by settling of aluminum were assumed to be proportional to aluminum concentration and reservoir area. The constant of proportionality, as a function of DOC concentration, was established by simulations in each of five reservoirs that differed in DOC concentration. In addition to computing dilution factors, the project determined dilution factors that would be protective with the same statistical basis (frequency of exceedance of the chronic standard) as dilutions computed for streams at the 7-day-average 10-year-recurrence annual low flow (the 7Q10). Low-flow dilutions are used for permitting so that receiving waters are protected even at the worst-case flow levels. The low-flow dilution factors that give the same statistical protection are the lowest annual 7-day-average dilution factors with a recurrence of 10 years, termed 7DF10s. Determination of 7DF10 values for reservoirs required that long periods of record be simulated so that dilution statistics could be determined. Dilution statistics were simulated for 13 reservoirs from 1960 to 2004 using U.S. Geological Survey Firm-Yield Estimator software to model reservoir inputs and outputs and present-day values of filter-effluent discharge and aluminum concentration. Computed settling velocities ranged from 0 centimeters per day (cm/d) at DOC concentrations of 15.5 milligrams per liter (mg/L) to 21.5 cm/d at DOC concentrations of 2.7 mg/L. The 7DF10 values were a function of aluminum effluent discharged. At current (2009) effluent discharge rates, the 7DF10 values varied from 1.8 to 115 among the 13 reservoirs. In most cases, the present-day (2009) discharge resulted in receiving water concentrations that did not exceed the standard at the 7DF10. Exceptions were one reservoir with a very small area and three reservoirs with high concentrations of DOC. Maximum permissible discharges were determined for water-treatment plants by adjusting discharges upward in simulations until the 7DF10 resulted in reservoir concentrations that just met the standard. In terms of aluminum flux, these discharges ranged from 0 to 28 kilograms of aluminum per day.

Massachusetts↗

Salinity trends in the upper Colorado River basin upstream from the Grand Valley Salinity Control Unit, Colorado, 1986-2003

In 1974, the Colorado River Basin Salinity Control Act was passed into law. This law was enacted to address concerns regarding the salinity content of the Colorado River. The law authorized various construction projects in selected areas or "units" of the Colorado River Basin intended to reduce the salinity load in the Colorado River. One such area was the Grand Valley Salinity Control Unit in western Colorado. The U. S. Geological Survey has done extensive studies and research in the Grand Valley Salinity Control Unit that provide information to aid the U.S. Bureau of Reclamation and the Natural Resources Conservation Service in determining where salinity-control work may provide the best results, and to what extent salinity-control work was effective in reducing salinity concentrations and loads in the Colorado River. Previous studies have indicated that salinity concentrations and loads have been decreasing downstream from the Grand Valley Salinity Control Unit, and that the decreases are likely the result of salinity control work in these areas. Several of these reports; however, also document decreasing salinity loads upstream from the Grand Valley Salinity Control Unit. This finding was important because only a small amount of salinity-control work was being done in areas upstream from the Grand Valley Salinity Control Unit at the time the findings were reported (late 1990’s). As a result of those previous findings, the U.S. Bureau of Reclamation entered into a cooperative agreement with the U.S. Geological Survey to investigate salinity trends in selected areas bracketing the Grand Valley Salinity Control Unit and regions upstream from the Grand Valley Salinity Control Unit. The results of the study indicate that salinity loads were decreasing upstream from the Grand Valley Salinity Control Unit from 1986 through 2003, but the rates of decrease have slowed during the last 10 years. The average rate of decrease in salinity load upstream from the Grand Valley Salinity Control Unit was 10,700 tons/year. This accounts for approximately 27 percent of the decrease observed downstream from the Grand Valley Salinity Control Unit. Salinity loads were decreasing at the fastest rate (6,950 tons/year) in Region 4, which drains an area between the Colorado River at Cameo, Colorado (station CAMEO) and Colorado River above Glenwood Springs, Colorado (station GLEN) streamflow-gaging stations. Trends in salinity concentration and streamflow were tested at station CAMEO to determine if salinity concentration, streamflow, or both are controlling salinity loads upstream from the Grand Valley Salinity Control Unit. Trend tests of individual ion concentrations were included as potential indicators of what sources (based on mineral composition) may be controlling trends in the upper Colorado. No significant trend was detected for streamflow from 1986 to 2003 at station CAMEO; however, a significant downward trend was detected for salinity concentration. The trend slope indicates that salinity concentration is decreasing at a median rate of about 3.54 milligrams per liter per year. Five major ions (calcium, magnesium, sodium, sulfate, and chloride) were tested for trends. The results indicate that processes within source areas with rock and soil types (or other unidentified sources) bearing calcium, sodium, and sulfate had the largest effect on the downward trend in salinity load upstream from station CAMEO. Downward trends in salinity load resulting from ground-water sources and/or land-use change were thought to be possible reasons for the observed decreases in salinity loads; however, the cause or causes of the decreasing salinity loads are not fully understood. A reduction in the amount of ground-water percolation from Region 4 (resulting from work done through Federal irrigation system improvement programs as well as privately funded irrigation system improvements) has helped reduce annual salinity load from Region 4 by approximately 7,400 tons. This amount is equal to about 5.9 percent of the total decrease (125,000 tons, or about 6,950 tons/year) estimated to have occurred in Region 4 during water years 1986 through 2003. A geographic information system was used to quantify the change in the amount of irrigated land upstream from the Grand Valley Salinity Control Unit from 1993 through 2000. These data indicated that the amount of irrigated land did not change substantially, thus indicating that the downward trends in salinity load did not result from land-use change.

Colorado↗

Assessing spatial variability of nutrients, phytoplankton, and related water-quality constituents in the California Sacramento–San Joaquin Delta at the landscape scale—2018 high resolution mapping surveys

Executive Summary This study examined the abundance and distribution of nutrients and phytoplankton in the tidal aquatic environments of the Sacramento–San Joaquin Delta (Delta) and Suisun Bay, comprising three spatial surveys conducted in May, July, and October of 2018 that used continuous underway high frequency sampling and measurements onboard a high-speed boat to characterize spatial variation across the extent of the Delta. The method used involves simultaneously collecting information about the concentration and spatial distribution of all major nutrient forms with analogous information about the major classes of phytoplankton and associated water-quality conditions. The results showed substantial variation across space and time, providing an unprecedented snapshot of the dynamic environmental processes that shape the ways nutrients interact with and affect aquatic habitats in the Delta. The purposes of this study were to improve our understanding of how hydrodynamics, landscape features, and aquatic primary productivity interact to drive nutrient cycling and transport in the Delta and to provide insights into the underlying processes most directly responsible for the conditions at the time of this study, and thus into the range of conditions that may be expected following the wide array of prospective future changes to the Delta. One major anticipated change at the time of this study was the planned upgrade to the Sacramento Regional Wastewater Treatment Plant, but the study also informs our understanding of potential effects from other changes to the Delta, such as those caused by other nutrient-management actions, flow actions, large-scale wetland restoration, drought, flood, levee failure, and changes to water management. Nutrient loading is the primary driver of nutrient concentrations in the Delta, but several other major drivers interact to shape their distribution and effects: geomorphology, hydrodynamics, landscape features, and aquatic productivity. Hydrodynamics affect timescales of transport and dilution of nutrient loads in the Delta. During transit through the system, channel geometry, tidal mixing, and water exports affect hydrodynamics in diverse ways that influence water-residence and transport times, thereby markedly affecting the range of times during which natural internal cycling can alter nutrient concentrations and forms. Channel geometry and location shape tidal energy and river currents into these observed dynamics. Interactions with Delta aquatic landscapes such as herbaceous tidal marsh, submerged aquatic vegetation, and large expanses of intertidal or subtidal sediments (all highly productive landscapes) exert demand on available nutrient supplies but can also simultaneously transform and generate nutrients. Finally, while phytoplankton require nutrients to sustain production and thus are a potential nutrient sink, the amount and form of nutrients also can influence the occurrence of harmful algal blooms (HABs) that adversely affect aquatic organisms as well as affect the occurrence of beneficial algal blooms that result in production of algae that are favorable for imperiled Delta pelagic aquatic food webs. The surveys revealed a complex mosaic of spatial variation, with nutrient concentrations varying from near zero to well above concentrations considered eutrophic; nutrient concentrations were more often related to the extent of hydrologic transport and mixing than to specific geographic locations or to specific landscape features. Similarly, the surveys identified phytoplankton abundance ranging from near detection to the level of large phytoplankton blooms, with large variation in phytoplankton community composition. Although the study occurred during a period of low bloom activity, phytoplankton productivity appeared to be the strongest potential sink for inorganic nutrients in the Delta, indicating that it is a larger control on nutrient concentrations and distribution than previously understood. Cycling and transformation within the water column only appeared to substantially lower total nutrient concentrations at the longest estimated transport timescales. Contrary to expectations, we did not observe substantial nutrient depletion near landscape-scale features such as open-water habitats, submerged aquatic vegetation beds, extensive wetlands, or exposed sediments, indicating that these habitat types did not act as major sinks for nutrients in the Delta during these surveys. These results indicated that nutrient reduction efforts may have the greatest effect on pelagic phytoplankton productivity in the more productive reaches of the Delta and estuary, but these effects are unlikely to be magnified by changes to nutrient loss within the Delta over conceivable changes in flow conditions, Delta water management actions, or large-scale wetland restoration activities. Nevertheless, local processes were shown to cause substantial loss, and thus integrating of nutrient effects with other indicators of aquatic habitat conditions will help inform planning future actions at specific sites. Finally, we note that the primary contribution of this study was intended to be the survey data themselves. Aside from the results highlighted in this report, the surveys are a benchmark against which future environmental change may be evaluated, including changes to nutrient management or water exports, drought, large-scale wetland restoration, and climate change. Further, although we highlight some of the main findings from the surveys in this report, the necessarily limited scope precludes examination of many topics for which these surveys may be highly informative. To facilitate the utility of these data to stakeholders, managers, and researchers, we have released the data online (Bergamaschi and others, 2020) and created an online data exploration portal ( https:​//ca.water​.usgs.gov/​bay-​delta/​2018-​delta-​wide-​mapping-​surveys.html ) where users may query the surveys in a variety of ways to test hypotheses, examine relationships, assess spatial trends, and download data. The data exploration portal is intended to be an immersive experience that allows users to gain greater understanding of the complex interactions that shape Delta aquatic environments. This report is intended as a companion to the portal, allowing the reader to challenge and further explore the highlighted findings. This study was a collaboration between the U.S. Geological Survey and the Delta Regional Monitoring Program, with additional funding provided from U.S. Geological Survey Cooperative Matching Funds Program.

California↗

Potential effects of alterations to the hydrologic system on the distribution of salinity in the Biscayne aquifer in Broward County, Florida

To address concerns about the effects of water-resource management practices and rising sea level on saltwater intrusion, the U.S. Geological Survey in cooperation with the Broward County Environmental Planning and Community Resilience Division, initiated a study to examine causes of saltwater intrusion and predict the effects of future alterations to the hydrologic system on salinity distribution in eastern Broward County, Florida. A three-dimensional, variable-density solute-transport model was calibrated to conditions from 1970 to 2012, the period for which data are most complete and reliable, and was used to simulate historical conditions from 1950 to 2012. These types of models are typically difficult to calibrate by matching to observed groundwater salinities because of spatial variability in aquifer properties that are unknown, and natural and anthropogenic processes that are complex and unknown; therefore, the primary goal was to reproduce major trends and locally generalized distributions of salinity in the Biscayne aquifer. The methods used in this study are relatively new, and results will provide transferable techniques for protecting groundwater resources and maximizing groundwater availability in coastal areas. The model was used to (1) evaluate the sensitivity of the salinity distribution in groundwater to sea-level rise and groundwater pumping, and (2) simulate the potential effects of increases in pumping, variable rates of sea-level rise, movement of a salinity control structure, and use of drainage recharge wells on the future distribution of salinity in the aquifer. Results from the simulation of historical conditions indicate that the model generally represents the observed greater westward extent of elevated salinity in the central part of the intruded area relative to the northern and southernmost parts of the intruded area. Results of sensitivity testing indicate that the extent of elevated salinity is most sensitive to pumping in areas where the source of saltwater is largely offshore, from the Atlantic Ocean, and is most sensitive to sea-level rise in areas where the source of salinity is downward leakage of brackish water from canals. Simulations of future scenarios indicate that increases in pumping near the existing interface may cause the interface to advance and decreases in pumping may cause it to retreat. Climatic effects, such as periods of prolonged drought or high precipitation, may augment or counteract long-term effects of changes in pumping on aquifer salinity at well fields. With increasing rates of sea-level rise, the freshwater-saltwater interface advances progressively inland, and flow-averaged salinities at well fields near the existing interface increase commensurately. Hypothetical southeastward (downstream) re-positioning of the existing G–54 salinity-control structure may prevent the interface from moving northwestward along and near the North New River canal, but beneficial effects are localized. Implementation of freshwater recharge wells in the city of Hallandale Beach may also have only a localized freshening effect in the aquifer and little appreciable effect on the freshwater-saltwater interface or on concentrations of salinity at well fields. Model accuracy and use are limited by uncertainty in the physical properties and boundary conditions of the system, uncertainty in historical and future conditions, and generalizations made in the mathematical relationships used to describe the physical processes of groundwater flow and transport. Because of these limitations, model results should be considered in relative rather than absolute terms. Nonetheless, model results do provide useful information on the relative scale of response of the system to changes in pumping distribution, sea-level rise, and mitigation activities.

Florida↗

Ground-Water Quality and Potential Effects of Individual Sewage Disposal System Effluent on Ground-Water Quality in Park County, Colorado, 2001-2004

In 2000, the U.S. Geological Survey, in cooperation with Park County, Colorado, began a study to evaluate ground-water quality in the various aquifers in Park County that supply water to domestic wells. The focus of this study was to identify and describe the principal natural and human factors that affect ground-water quality. In addition, the potential effects of individual sewage disposal system (ISDS) effluent on ground-water quality were evaluated. Ground-water samples were collected from domestic water-supply wells from July 2001 through October 2004 in the alluvial, crystalline-rock, sedimentary-rock, and volcanic-rock aquifers to assess general ground-water quality and effects of ISDS's on ground-water quality throughout Park County. Samples were analyzed for physical properties, major ions, nutrients, bacteria, and boron; and selected samples also were analyzed for dissolved organic carbon, human-related (wastewater) compounds, trace elements, radionuclides, and age-dating constituents (tritium and chlorofluorocarbons). Drinking-water quality is adequate for domestic use throughout Park County with a few exceptions. Only about 3 percent of wells had concentrations of fluoride, nitrate, and (or) uranium that exceeded U.S. Environmental Protection Agency national, primary drinking-water standards. These primary drinking-water standards were exceeded only in wells completed in the crystalline-rock aquifers in eastern Park County. Escherichia coli bacteria were detected in one well near Guffey, and total coliform bacteria were detected in about 11 percent of wells sampled throughout the county. The highest total coliform concentrations were measured southeast of the city of Jefferson and west of Tarryall Reservoir. Secondary drinking-water standards were exceeded more frequently. About 19 percent of wells had concentrations of one or more constituents (pH, chloride, fluoride, sulfate, and dissolved solids) that exceeded secondary drinking-water standards. Currently (2004), there is no federally enforced drinking-water standard for radon in public water-supply systems, but proposed regulations suggest a maximum contaminant level of 300 picocuries per liter (pCi/L) and an alternative maximum contaminant level of 4,000 pCi/L contingent on other mitigating remedial activities to reduce radon levels in indoor air. Radon concentrations in about 91 percent of ground-water samples were greater than or equal to 300 pCi/L, and about 25 percent had radon concentrations greater than or equal to 4,000 pCi/L. Generally, the highest radon concentrations were measured in samples collected from wells completed in the crystalline-rock aquifers. Analyses of ground-water-quality data indicate that recharge from ISDS effluent has affected some local ground-water systems in Park County. Because roughly 90 percent of domestic water used is assumed to be recharged by ISDS's, detections of human-related (wastewater) compounds in ground water in Park County are not surprising; however, concentrations of constituents associated with ISDS effluent generally are low (concentrations near the laboratory reporting levels). Thirty-eight different organic wastewater compounds were detected in 46 percent of ground-water samples, and the number of compounds detected per sample ranged from 1 to 17 compounds. Samples collected from wells with detections of wastewater compounds also had significantly higher (p-value < 0.05) chloride and boron concentrations than samples from wells with no detections of wastewater compounds. ISDS density (average subdivision lot size used to estimate ISDS density) was related to ground-water quality in Park County. Chloride and boron concentrations were significantly higher in ground-water samples collected from wells located in areas that had average subdivision lot sizes of less than 1 acre than in areas that had average subdivision lot sizes greater than or equal to 1 acre. For wells completed in the crystalline-

Scientific Investigations Report↗

Lunar remote sensing and measurements

Remote sensing and measurements of the Moon from Apollo orbiting spacecraft and Earth form a basis for extrapolation of Apollo surface data to regions of the Moon where manned and unmanned spacecraft have not been and may be used to discover target regions for future lunar exploration which will produce the highest scientific yields. Orbital remote sensing and measurements discussed include (1) relative ages and inferred absolute ages, (2) gravity, (3) magnetism, (4) chemical composition, and (5) reflection of radar waves (bistatic). Earth-based remote sensing and measurements discussed include (1) reflection of sunlight, (2) reflection and scattering of radar waves, and (3) infrared eclipse temperatures. Photographs from the Apollo missions, Lunar Orbiters, and other sources provide a fundamental source of data on the geology and topography of the Moon and a basis for comparing, correlating, and testing the remote sensing and measurements. Relative ages obtained from crater statistics and then empirically correlated with absolute ages indicate that significant lunar volcanism continued to 2.5 b.y. (billion years) ago-some 600 m.y. (million years) after the youngest volcanic rocks sampled by Apollo-and that intensive bombardment of the Moon occurred in the interval of 3.84 to 3.9 b.y. ago. Estimated fluxes of crater-producing objects during the last 50 m.y. agree fairly well with fluxes measured by the Apollo passive seismic stations. Gravity measurements obtained by observing orbiting spacecraft reveal that mare basins have mass concentrations and that the volume of material ejected from the Orientale basin is near 2 to 5 million km 3 depending on whether there has or has not been isostatic compensation, little or none of which has occurred since 3.84 b.y. ago. Isostatic compensation may have occurred in some of the old large lunar basins, but more data are needed to prove it. Steady fields of remanent magnetism were detected by the Apollo 15 and 16 subsatellites, and the lunar dipole field was revised to no more than 6x 10 19 gauss. High-resolution mapping of fields of weak remanent magnetism (to 0.1 gamma) was made possible by the Apollo plasma and energetic-particle experiment. Although the causes of remanent magnetism are poorly understood, correlations with geologic units suggest the results may ultimately have farreaching significance to lunar history. Maria are much less structured by strong surface magnetic anomalies than the highlands. The strongest anomalies are associated with ejecta of farside basins, plains materials filling pre-Imbrian craters, and other old Imbrian to pre-Imbrian units. The high remanent fields could be due to cooling of ejecta units in an ancient magnetic field, lunar regolith maturity, extensive reworking and disruption of a magnetized layer, or simply surface roughness. Orbital geochemical experiments have shown that lunar high lands have larger Al: Si ratios and smaller Mg: Si ratios than maria. These two ratios are inversely related on a regional basis. With the exception of fresh craters, albedo and Al : Si ratios vary directly, showing that compositional differences as well as exposure of fresh materials are responsible for high albedos. Statistically treated data show that geologic contacts and compositional boundaries are concentric and can be roughly matched. Some craters on mare material have penetrated the mare fill, bringing highland-type materials to the surface. Natural radioactivity from thorium, potassium, and uranium is inversely correlated with elevation. Mare regions are enriched in iron, titanium, and magnesium relative to the highlands. Orbital bistatic-radar results provide estimates of surface roughness at two scale lengths (about 30 m and 250 m), which agree with visual estimates of roughness. The dielectric constant of the lunar surface, where sampled, is uniform to 13-cm radar and near 3. Slope frequency distributions measured by the radar vary and

Professional Paper↗

Water and Streambed-Sediment Quality in the Upper Elk River Basin, Missouri and Arkansas, 2004-06

The U.S. Geological Survey, in cooperation with the Missouri Department of Natural Resources, collected water and streambedsediment samples in the Upper Elk River Basin in southwestern Missouri and northwestern Arkansas from October 2004 through December 2006. The samples were collected to determine the stream-water quality and streambed-sediment quality. In 1998, the Missouri Department of Natural Resources included a 21.5-mile river reach of the Elk River on the 303(d) list of impaired waters in Missouri as required by Section 303(d) of the Federal Clean Water Act. The Elk River is on the 303(d) list for excess nutrient loading. The total phosphorus distribution by decade indicates that the concentrations since 2000 have increased significantly from those in the 1960s, 1980s, and 1990s. The nitrate as nitrogen (nitrate) concentrations also have increased significantly in post-1985 from pre-1985 samples collected at the Elk River near Tiff City. Concentrations have increased significantly since the 1960s. Concentrations in the 1970s and 1980s, though similar, have increased from those in the 1960s, and the concentrations from the 1990s and 2000s increased still more. Nitrate concentrations significantly increased in samples that were collected during large discharges (greater than 355 cubic feet per second) from the Elk River near Tiff City. Nitrate concentrations were largest in Indian Creek. Several sources of nitrate are present in the basin, including poultry facilities in the upper part of the basin, effluent inflow from communities of Anderson and Lanagan, land-applied animal waste, chemical fertilizer, and possible leaking septic systems. Total phosphorus concentrations were largest in Little Sugar Creek. The median concentration of total phosphorus from samples from Little Sugar Creek near Pineville was almost four times the median concentration in samples from the Elk River near Tiff City. Median concentrations of nutrient species were greater in the stormwater samples than the median concentrations in the ambient samples. Nitrate concentrations in stormwater samples ranged from 133 to 179 percent of the concentration in the ambient samples. The total phosphorus concentrations in the stormwater samples ranged from about 200 to more than 600 percent of the concentration in the ambient samples. Base-flow conditions as reflected by the seepage run of the summer of 2006 indicate that 52 percent of the discharge at the Elk River near Tiff City is contributed by Indian Creek. Little Sugar Creek contributes 32 percent and Big Sugar Creek 9 percent of the discharge in the Elk River near Tiff City. Only about 7 percent of the discharge at Tiff City comes from the mainstem of the Elk River. Concentrations of dissolved ammonia plus organic nitrogen as nitrogen, dissolved ammonia as nitrogen, dissolved phosphorus, and dissolved orthophosphorus were detected in all streambed-sediment leachate samples. Concentrations of leachable nutrients in streambed-sediment samples generally tended to be slightly larger along the major forks of the Elk River as compared to tributary sites, with sites in the upper reaches of the major forks having among the largest concentrations. Concentrations of leachable nutrients in the major forks generally decreased with increasing distance downstream.

Scientific Investigations Report↗