Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Progress Report”

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

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

At least 1,603 records · Page 89Linked to original sources

Monitoring to assess progress toward meeting the total maximum daily load for phosphorus in the Assabet River, Massachusetts: phosphorus loads, 2008 through 2010

Wastewater discharges to the Assabet River contribute substantial amounts of phosphorus, which support accumulations of nuisance aquatic plants that are most evident in the river’s impounded reaches during the growing season. To restore the Assabet River’s water quality and aesthetics, the U.S. Environmental Protection Agency required the major wastewater-treatment plants in the drainage basin to reduce the amount of phosphorus discharged to the river by 2012. From October 2008 to December 2010, the U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection and in support of the requirements of the Total Maximum Daily Load for Phosphorus, collected weekly flow-proportional, composite samples for analysis of concentrations of total phosphorus and orthophosphorus upstream and downstream from each of the Assabet River’s two largest impoundments: Hudson and Ben Smith. The purpose of this monitoring effort was to evaluate conditions in the river before enhanced treatment-plant technologies had effected reductions in phosphorus loads, thereby defining baseline conditions for comparison with conditions following the mandated load reductions. The locations of sampling sites with respect to the impoundments enabled examination of the impoundments’ effects on phosphorus sequestration and on the transformation of phosphorus between particulate and dissolved forms. The study evaluated the differences between loads upstream and downstream from the impoundments throughout the sampling period and compared differences during two seasonal periods of relevance to aquatic plants: April 1 through October 31, the growing season, and November 1 through March 31, the nongrowing season, when existing permit limits allowed average monthly wastewater-treatment-plant-effluent concentrations of 0.75 milligram per liter (growing season) or 1.0 milligram per liter (nongrowing season) for total phosphorus. At the four sampling sites during the growing season, median weekly total phosphorus loads ranged from 110 to 190 kilograms (kg) and median weekly orthophosphorus loads ranged from 17 to 41 kg. During the nongrowing season, median weekly total phosphorus loads ranged from 240 to 280 kg and median weekly orthophosphorus loads ranged from 56 to 66 kg. During periods of low and moderate streamflow, estimated loads of total phosphorus upstream from the Hudson impoundment generally exceeded those downstream during the same sampling periods throughout the study; orthophosphorus loads downstream from the impoundment were typically larger than those upstream. When storm runoff substantially increased the streamflow, loads of total phosphorus and orthophosphorus both tended to be larger downstream than upstream. At the Ben Smith impoundment, both total phosphorus and orthophosphorus loads were generally larger downstream than upstream during low and moderate streamflow, but the differences were not as pronounced as they were at the Hudson impoundment. High flows were also associated with substantially larger total phosphorus and orthophosphorus loads downstream than those entering the impoundment from upstream. In comparing periods of growing- and nongrowing-season loads, the same patterns of loads entering and leaving were observed at both impoundments. That is, at the Hudson impoundment, total phosphorus loads entering the impoundment were greater than those leaving it, and orthophosphorus loads leaving the impoundment were greater than those entering it. At the Ben Smith impoundment, both total phosphorus and orthophosphorus loads leaving the impoundment were greater than those entering it. However, the loads were greater during the nongrowing seasons than during the growing seasons, and the net differences between upstream and downstream loads were about the same. The results indicate that some of the particulate fraction of the total phosphorus loads is sequestered in the Hudson impoundment, where particulate phosphorus probably undergoes some physical and biogeochemical transformations to the dissolved form orthophosphorus. The orthophosphorus may be taken up by aquatic plants or transported out of the impoundments. The results for the Ben Smith impoundment are less clear and suggest net export of total phosphorus and orthophosphorus. Differences between results from the two impoundments may be attributable in part to differences in their sizes, morphology, unmonitored tributaries, riparian land use, and processes within the impoundments that have not been quantified for this study.

Massachusetts↗

Geology and assessment of undiscovered oil and gas resources of the Sverdrup Basin Province, Arctic Canada, 2008

The Sverdrup Basin Province, an area of 515,000 square kilometers on the northern margin of North America, extends 1,300 kilometers across the Canadian Arctic Islands from near the Mackenzie Delta to northern Ellesmere Island. It consists of an intracratonic late Paleozoic to early Cenozoic rift-sag basin and a Mesozoic rift shoulder that bounds it on the north. Basin inception was Mississippian, manifested by deposition of nonmarine strata in rift basins, followed by Pennyslvanian marine transgression, which began with evaporites and progressed to Permian carbonate and clastic deposition at basin fringes and organic-rich marine strata in the basin center. Sediment transport was both northward from North America and southward from a now-subsided or rifted-away landmass to the north. Mesozoic strata indicate continued marine deposition, including both organic-rich, fine-grained rocks deposited during highstands and progradational deltaic sequences. A new episode of rifting began in Middle Jurassic time and culminated in the opening of the Canada Basin by Early Cretaceous seafloor spreading. The Sverdrup Rim formed as the rift shoulder between North America and the thinned, subsided crust to the north. Widespread Upper Cretaceous organic-rich shales were deposited during the major transgression induced by Canada Basin opening, followed by an influx of coarser east-derived detritus. In Paleogene time, incipient North Atlantic seafloor spreading caused deformation in northeasternmost North America, producing uplifts that shed detritus westward across the Sverdrup Basin. Tight folding and thrusting resulting from the Eurekan orogeny took place in the eastern part of the basin during the Eocene, with decreasing intensity of deformation westward. Since deformation ended in late Eocene time, little significant tectonism or deposition has taken place. Two petroleum systems were defined in the Sverdrup Basin Province. Upper Paleozoic marine shale generated petroleum beginning in the Early Triassic, but this petroleum system was not quantitatively assessed because reservoir quality in adjacent strata is poor, the rocks are mostly overmature, and subsequent deformation likely affected trap integrity. The second petroleum system was sourced by Lower Triassic strata rich in oil-prone organic matter. Oil was generated during Paleogene burial synchronous with Eurekan deformation, and the oil migrated into Triassic and Jurassic deltaic, shallow marine and nonmarine strata. However, most of the oil may have escaped during deformation and subsequent uplift and erosion, which probably caused oil to be displaced from traps by gas expansion. The population of undiscovered accumulations was characterized as likely to include stratigraphically trapped and small, structurally trapped accumulations, with a median size of 80 million barrels of oil (MMBO); the number of undiscovered accumulations was estimated to be between 1 and 50, with the most likely number being 10. The resulting estimate of undiscovered, technically recoverable, conventional oil resources is 61 to 1,255 MMBO, with a mean of 427 MMBO. Undiscovered, technically recoverable, conventional gas resources are estimated at 4.95 trillion cubic feet (TCF), with slightly more than half of that in nonassociated gas accumulations. A third petroleum system in the adjacent Amerasia Basin Province to the north was considered somewhat likely to contain accumulations on the Sverdrup Rim. Deeply buried Upper Jurassic, Upper Cretaceous, and Eocene organic-rich strata probably generated oil that may have migrated up the continental slope into Triassic to Paleogene sandstones on the Sverdrup Rim. Based on analogy with the Barrow Arch in Alaska, a median of 20 accumulations was estimated, with accumulation volumes as much as 2,500 MMBO and a median of 100 MMBO. The probability of at least one accumulation of the minimum size assessed (50 MMBO) was estimated at 0.22. The resulting estimate of undiscovered, technically recoverable, conventional oil resources is 0 to 2,679 MMBO, with a mean of 424 MMBO. Mean estimates for associated and nonassociated gas are 1.3 and 2.3 TCF, respectively.

Sverdrup Basin↗

Spatial distribution, temporal variability, and chemistry of the salt wedge in the lower Charles River, Massachusetts, June 1998 to July 1999

The Charles River is of great recreational and ecological value to the Boston metropolitan region and the Commonwealth of Massachusetts. It is also the focus of the U.S. Environmental Protection Agency (USEPA) Region I, Clean Charles 2005 Task Force. The main goal of the Task Force is to make the Charles River 'fishable and swimmable' by the year 2005. Achieving 'fishable and swimmable' conditions will require continued progress in addressing a range of environmental conditions now degrading water quality, including the infiltration of saltwater from Boston Harbor into the freshwater Charles River.To better understand the pattern of saltwater intrusion, the U.S. Geological Survey (USGS), in cooperation with the U.S. Environmental Protection Agency (USEPA), Massachusetts Department of Environmental Management (MADEM), and New England Interstate Water Pollution Control Commission (NEIWPCC), collected data on the spatial distribution, temporal variability, and chemistry of the saltwater that entered the lower Charles River from June 1998 to July 1999. The purpose of this investigation is to extend and complement a regional-scale study of Charles River water quality conducted in 1996 (T. Faber, U.S. Environmental Protection Agency, written commun., 1997), and the ongoing water monitoring activities of the Massachusetts Water Resources Authority (MWRA) and the Charles River Watershed Association (CRWA). The data collected by this investigation supports the Clean Charles 2005 Task Force by providing detailed information concerning a major factor limiting 'fishable and swimmable' conditions in the lower Charles River. Finally, the study will be used to assist current planning efforts of the Metropolitan District Commission (MDC) to restore the historic parklands of the lower Charles River.The 'Basin' is the local term for the reach of the Charles River that begins at the Watertown Dam in Watertown, Mass., and extends about 8 mi through suburban and urban areas to Boston Harbor. Discharge to the harbor is controlled by the 'new' Charles River Dam in Boston (fig. 1). The Basin was created by construction of the 'old' Charles River Dam in 1908 to solve Boston's sanitary problems. Prior to the building of the old Charles River Dam, the lower Charles River was a tidal estuary in which the water levels rose and fell twice daily with the tidal cycle. Low tide would expose untreated sewage that was discharged directly into the river. Exposed sewage created noxious odors and served as a breeding ground for mosquitoes that caused sporadic epidemics of malaria and yellow fever (Jobin, 1998). Damming of the river interrupted the normal tidal cycle and flooded the estuary by creating a freshwater pool (the Basin) that had a constant water elevation of about 0.8 meters (m) above mean sea level. Flooding of the estuary initially improved sanitary conditions and the Basin became a source of enjoyment for the local population and the focus of a large waterfront park in Boston and Cambridge (Jobin, 1998).Although the infiltration of saltwater from the harbor into the Basin was anticipated when the old Charles River Dam was built, neither the magnitude nor the consequences of the infiltration was considered. By 1975, the Metropolitan District Commission (MDC) determined that harbor water covered about 80 percent of the river bottom within the Basin and composed about 50 percent of its depth. The MDC also concluded that fish kills and odors in the spring of 1975 were likely the result of the sulfide-rich saltwater mixing with the overlying freshwater (Metropolitan District Commission, 1975).Saltwater from Boston Harbor that enters the Basin is known as the 'salt wedge' because of the shape it assumes as it moves upstream. Freshwater discharge from upstream pushes against the intruding harbor water until the density differences cause stratification to occur; the freshwater then overrides the denser harbor water (Fischer and others, 1979). The depth from

Water-Resources Investigations Report↗

Investigation of artificial recharge of aquifers in Nebraska

Large withdrawals of ground water for irrigation are causing progressive declines of ground-water levels in some areas of Nebraska. An investigation was conducted to determine the technical feasibility of artificially recharging aquifers through wells and through surface spreading by means of impoundments, pits, and canals. Information gained from a literature search and from preliminary tests was used to design several artificial-recharge experiments. The experiments showed that large quantities of water can be recharged through wells and by surface spreading if conditions are favorable. In the well experiments, about 0.5 billion gallons of water from an aquifer recharged by the Platte River was transported 3 miles by pipeline and recharged through a well into a Pleistocene sand and gravel aquifer near Aurora where ground-water levels are declining. The recharge rate was about 730 gallons per minute during two tests of 6 and 8 months duration. The rise in ground-water levels due to recharge extended more than a mile from the recharge well in both tests. The pattern of ground-water-level buildup during the 8-month test was similar to that during the 6-month test. Two-thirds of the way through the test the rate of water-level buildup in the recharge well increased greatly because casing failure allowed a large amount of sediment to enter the well. Although the chemical quality of the recharge water was markedly different from that of the native aquifer water, no evidence of clogging due to chemical reaction was detected; also, there was no evidence of clogging due to air entrainment or bacterial growth. Evaluation of water-level changes in the recharge well and in observation wells during the 6-month test indicated some clogging of the aquifer in the immediate vicinity of the recharge well due to a small amount (0.04 milligram per liter) of fine sediment in the recharge water. Analysis of water-level buildup in the recharge well during the 6-month test indicated that recharge could have continued at a rate of about 700 gallons per minute for several years before rehabilitation of the recharge well would have been necessary. In surface-spreading experiments, the maximum sustained infiltration rates from a 24-foot-diameter ring infiltrometer set in loess at the well-recharge site near Aurora were about 0.5 foot per day during a 140day test. The recharge water caused water levels in a perched zone of saturation at depths between 36 and 38 feet to rise 15 feet, to within 21 feet of the surface, indicating that if the test had continued or if the impoundment area had been larger, the water level in the perched zone of saturation might have risen to the surface thereby reducing the infiltration rate. The maximum sustained infiltration rate from a similar experiment in the Sand Hills near Tryon was 11 feet per day. Perching layers also retarded downward infiltration and caused lateral movement of water in the subsurface at this site. Infiltration rates from re-use pits near the Aurora site ranged from 0.01 to 1.60 feet per day, indicating that the permeability of the subsurface material is extremely variable and that perching layers probably are absent in some areas. Flow measurements in an irrigation canal excavated in loess near Farwell indicate an infiltration rate of 0.36 foot per day from a 2.7mile reach.

Nebraska↗

Recent planform changes in the Upper Mississippi River

Geomorphic changes in the Upper Mississippi River (UMR) have long been a concern of river agencies charged with maintaining and restoring river habitat (GREAT 1980; Jackson et al. 1981; USFWS 1992). Large meandering alluvial rivers like the UMR are expected to constantly change and adjust their fluvial landforms within their riparian corridors as a result of the natural interaction of hydrologic processes, sediment movement, and vegetation over time. However, present geomorphic changes in the UMR reflect altered hydrologic, hydraulic, and sediment conditions caused by regulated flows, constructed agricultural levees and navigation dams, altered land use in the watershed, and climate change. Levees reduce lateral hydrologic and sediment connectivity between channels and floodplains on many tributaries and on the Mississippi River downstream of Pool 13. Between each of the dams are a repeating series of landforms associated with tailwater, intermediate, and impounded conditions. The dams maintain a minimum water level, thus creating many off-channel areas that act as sediment traps. Whereas high-head dams cut off sedimentological connectivity longitudinally through the river corridor (Skalak et al., 2013), low head dams on the UMR only slightly altered transport longitudinally. Deltaic-like sedimentation can be common in the impounded sections of dammed rivers. Erosion of relict land surfaces that remained above the raised impounded water levels has been the dominant change in UMR impounded sections due to increased wind fetch leading to increased wave action. Even though upland sources of sediment from tributaries have decreased over the middle to late 20th century, increased annual precipitation, the interplay of increased variability in flood magnitudes from year to year, and more fall and winter flooding have likely changed erosion and sedimentation patterns in the UMR (Belby, et al., 2019). Paradoxically, monitoring and research indicates that the concentration of some water column constituents like total suspended solids and phosphorous has decreased during the 1991 to 2014 time period (Kreiling and Houser, 2016). In areas prone to increased sedimentation, bed elevations rise and thereby water depths are reduced at a given discharge, resulting in loss of fish habitat. Sediment deposition or erosion further influences water exchange rates between main channel and off-channel areas in the river by increasing resistance in connecting channels or enlarging existing connecting channels. Water depth and water exchange rates are the most prominent features describing habitat quality in the UMR (De Jager et al. 2018), and in some cases, the trajectory of planform change from heightened deposition promises to threaten deep backwater habitats particularly important for overwintering fish. Although information on the rate of vertical change in bed elevation is needed for a complete assessment of geomorphic change associated with the loss of deep backwater habitats, mapping planform changes over time (i.e., lateral changes between the land-water boundary) provide needed information on the location, potential cause, and progressive direction of deposition, especially in the mid sections between dams where deltaic processes are the most pronounced. Several types of planform changes have been observed and identified as concerns. For example, island loss in the large impounded areas of the upper part of the UMR was one of the concerns identified by river managers in the 1980s and 90s, and subsequently island construction became a common form of restoration implemented by the Upper Mississippi River Restoration (UMRR) Program (USACE 2012). Other subtler planform changes, such as channel bank erosion and delta formation in backwaters, are perceived to be important, but have largely gone unquantified. A systemwide reconnaissance of the UMR and IWW conducted in 1998 concluded that 14-percent of the river banks were eroding (Nakato and Anderson 1998). However, stabilization of existing river banks has never been widely pursued as a restoration measure, due to the high cost and uncertain benefits. Delta formation reduces the amount of backwater habitat; however, the deltas maintain and create a mix of riparian and aquatic habitats, and that is generally considered to be beneficial for wildlife and fish. If recent hydrologic trends of more frequent and longer duration flood events continue, a better understanding of planform changes can help in describing past changes, and then be used to forecast potential future trajectories of change. If UMR resource managers determine that past and forecasted conditions are undesirable, then UMRR projects could be identified and prioritized to address those concerns. Vegetative cover associations with landform changes have been used to detect and quantify planform changes in many rivers (Johnson 1985; Hiatt 2015; Volte et al. 2015). Freyer and Jefferson (2013) completed such a study in Pool 6 of the UMR using the landcover data from 12 dates over a 115-yr period, including the 1989, 2000, and 2010/2011 landcover/use (LCU) data from the UMRR Program. Planform change detected over the last 20 years represented by the UMRR Program data best reflect present-day geomorphic patterns, rates and processes. Changes occurring prior to dam construction and changes occurring soon after dam construction are likely not the same as those happening now, 50-70 years after dam construction and creation of the impoundments (McHenry et al., 1984; Bhowmik and Adams, 1986; WEST Consultants, 2000). The LCU data from each of the 1989, 2000, and 2010/2011 imagery was developed using similar methods and is available in a Geographical Information System (GIS) for the entire UMR and therefore provides the opportunity for a more comprehensive planform change analysis. This study used GIS overlays of LCU classes to map and quantify changes in planform features over two periods, looking specifically for depositional areas where terrestrial and wetland vegetation expanded at the expense of open water. The land expansion was grouped into four possible process-based types common in large floodplain rivers, some following that used by Lewin et al. (2017). The four types include: crevasse deltas emanating from a breach from a main channel through a natural levee or narrow floodplain into backwaters (crevasse deltas), tributary deltas expanding into backwaters (tributary deltas), deltaic bars at the upstream end of impoundments (impounded deltas), and linear-like bars extending from the downstream ends of narrow levees and remnant floodplains (bar-tail limbs). The methods deployed for change detection addressed possible errors from a variety of sources.

Illinois, Iowa, Minnesota, Missouri, Wisconsin↗

Water-withdrawal and groundwater-level fluctuations, physical and infiltration characteristics of Ohio River sediment, and isotope-based estimates of water sources relative to a riverbank filtration system in an alluvial aquifer, Louisville, Kentucky, 2000–20

The U.S. Geological Survey and the Louisville Water Company, from 2000 to 2020, cooperatively investigated physical and hydrologic conditions in the Ohio River alluvial aquifer and adjacent riverbed sediments in the northeast part of Louisville, Kentucky, to increase knowledge of processes affecting water sources to collector wells during operation of a riverbank filtration system (RBF). The RBF, in 2020, included an initial collector well (CW1) completed in 1999 and a tunnel and collector well system (TCWS) composed of four collector wells (CW2–CW5) completed in late 2010 that were connected by a tunnel that routed water to a pumping station and a treatment plant. Daily mean water-withdrawal rates typically ranged from about 10.4 to 20.7 million gallons per day (Mgal/d) from 2005 to 2007, then were cyclic from 2008 to 2010 and ranged from as low as about 9 to 13 Mgal/d in winter-spring months up to about 19 to 20 Mgal/d in summer-autumn months. From 2011 through 2017, peak daily mean water-withdrawal rates from CW2–CW5 fluctuated between winter-spring lows of about 20–30 Mgal/d to summer-autumn peak rates of about 41.6 Mgal/d in 2013 to about 48.2 Mgal/d in 2012. From 2017 to 2020, TCWS peak daily mean water-withdrawal rates decreased from about 45.7 Mgal/d in 2017 to about 37.5 Mgal/d in 2020 and peak daily mean rates from CW1 ranged from about 20.7 Mgal/d in 2020 to about 21.6 Mgal/d in 2018. Periods in 2019–20 with tunnel pumping station daily mean water-withdrawal rates smaller than about 30 Mgal/d had generally stable water levels in CW2 and indicated the sustainable yield of the TCWS collector wells under the prevailing hydrologic conditions. Drawdown during July–August of 2008 and 2009 was associated with CW1, the sole operating collector well. Groundwater-flow directions outside the CW1 cone of depression during July–August of 2008 and 2009 were from areas near the alluvial valley margin toward the Ohio River and in 2009 toward Harrods Creek. Collector wells CW2 through CW5 derived water during 2011–19 surveys from surface-water infiltration, as indicated by drawdown that extended beneath the Ohio River and Harrods Creek. Increases in drawdown toward the alluvial valley margins in July–August of 2013 and 2015, compared with July 2011, indicated less available recharge from those flow directions and decreased yield of groundwater from the aquifer after progressively longer TCWS water withdrawals. Drawdown at CW2–CW5 and at observation wells between them lessened in October 2019 and indicated that onshore groundwater levels can recover when collector well withdrawals decrease. Overlapping cones of depression in 2011–15 potentiometric surfaces centered around CW2, CW3, CW4, and CW5, and drawdown along divides between those wells indicated well interference between collector wells. Well-interference effects increased during 2013 and 2015 survey periods, despite decreased water-withdrawal rates compared with July 2011, and expanded to include CW1 when its use resumed. Groundwater-level fluctuations in observation wells near the Ohio River were consistent with gradual depletion of groundwater yield during tunnel pumping station operation at river elevations of about 420 feet above the National Geodetic Vertical Datum of 1929 (NGVD 29) and enhanced recharge of the aquifer when river elevations were above flood stage of about 430 feet above NGVD 29. Groundwater-level declines in observation wells near the Ohio River were most apparent in summer and some autumn months between 2012 and 2017 during periods when tunnel pumping station mean daily water-withdrawal rates were about 35 Mgal/d or more. Periods of relatively stable groundwater levels during 2019–20 in observation wells between CW2–CW5, when tunnel pumping station water-withdrawal rates were about 30–33 Mgal/d or less, also indicated that the prevailing sustainable yield was about 30 Mgal/d or less. Wells associated with CW3 and CW4 had the largest increases in slopes of best-fit lines relating water-level difference to water-withdrawal rate between 2012 and 2020, an indication of decreased efficiencies of those wells to supply water or of changes in groundwater yield around the well. Increased slopes of best-fit lines relating water-level difference to water-withdrawal rate from 2019 to 2020, when groundwater levels were otherwise relatively stable, also indicated that the specific capacity of one or more TCWS collector wells may have declined during that period. Riverbed surface elevation, roughness, and hardness changes between hydroacoustic surveys indicated active accretion and erosion of bed sediment offshore from the riverbank filtration system in 2002–4. Riverbed hydraulic conductivity computed from 2013 to 2015 in tests at 11 sites ranged from 0.2 to 230 feet per day. Ten temporary piezometer sites with adjacent operating collector wells had drawdowns that ranged from about 4 to 11.7 feet and hammer blows greater than or equal to 49 counts, which indicated a relatively great density of alluvial sediments and the potential for diminished riverbed infiltration rates. Contributions of Ohio River infiltration to collector wells estimated from deuterium and oxygen-18 isotopic ratios in representative Ohio River and groundwater samples were considerably larger in warmer months than cooler months and ranged from 80–90 percent in CW1–CW5 in May 2020 to 100 percent in May 2019 at CW2, CW4, and CW5. Groundwater contributions to water withdrawals by collector wells increased in cooler autumn months and ranged from 40 percent at CW3 and CW5 in November 2020 to 70 percent at CW2, CW4, and CW5 in November–December 2019. November 2019 and November 2020 isotopic compositions of water from the collector wells had enriched deuterium and oxygen-18 isotopic ratios consistent with mixtures of evaporation-modified water from the Ohio River and meteoric-origin groundwater.

Indiana, Kentucky↗

Tidal-flow, circulation, and flushing changes caused by dredge and fill in Tampa Bay, Florida

Tampa Bay, Florida, underwent extensive physical changes between 1880 and 1972 because of construction of causeways, islands, channels, and shoreline fills. These changes resulted in a progressive reduction in the quantity of tidal water that enters and leaves the bay. Dredging and filling also changed the magnitude and direction of tidal flows in large parts of the bay. A two-dimensional, finite-difference, hydrodynamic model was used to simulate flood, ebb, and residual transport of both water and a dissolved constituent for the physical conditions that existed in Tampa Bay during 1880 and 1972 and for the conditions that are likely to exist in 1985. The calibrated and verified model was used to hindcast water- and constituent-transport changes caused by construction in the bay between 1880 and 1972. The model was used also to forecast changes that can be expected to occur as a result of a major Federal dredging project scheduled for completion in 1985. The model forecasted transport changes caused by the Federal dredging project to be much less areally extensive than the corresponding transport changes caused by construction in Tampa Bay between 1880 and 1972. Dredging-caused changes of more than 50 percent in flood and ebb transport were computed to occur over only 8 or 9 square miles of the bay's 390-square-mile surface area between 1972 and 1985. The model computed that construction between 1880 and 1972 caused changes of similar magnitude over 58 square miles of the bay. Dredging-caused changes of more than 50 percent in residual transport were computed to occur over 58 square miles of the bay between 1972 and 1985. According to the model simulation, construction between 1880 and 1972 caused changes of similar magnitude over 167 square miles. Computations reveal historical tide-induced circulation patterns. The patterns consist of a series of about 20 interconnected circulatory features that range in diameter from 1 to 6 miles. Dredging- and construction-caused changes in size, position, shape, and intensity of the circulatory features increase tide-induced circulation and flushing throughout most of the bay. As a result of past and projected physical changes, the bay can and will more rapidly transfer waterborne constituents that have landward sources to the Gulf of Mexico. Conversely, the bay can and will more rapidly transfer constituents that have their source in the Gulf into the upper parts of the bay. Model results show that the bay can be functionally subdivided into eight circulation zones. The zones near the entrances to Tampa Bay and the major bay subunits, Hillsborough Bay and Old Tampa Bay, have several times greater average circulation than do adjacent and more landward zones. Circulation generally decreases from the Gulf of Mexico to the head of Hillsborough and Old Tampa Bays, with a striking exception in zone 3 in mid-Tampa Bay. This 10-mile section of the bay has significantly lower average circulation than that of adjacent zones, particularly for conditions in 1880. The section is thought to be a circulation constriction that reduces the potential transport of dissolved and suspended constituents. Circulation in the constricted section was computed as having increased 6 percent because of dredge and fill construction between 1880 and 1972. An additional increase of 21 percent, due to the Federal dredging project, is calculated. With these increases, this zone acts as less of a constriction than it did in 1880. This and other circulation increases may have contributed to increased bay salinity and to more rapid flushing of constituents from the bay to the Gulf of Mexico.

Water Supply Paper↗

Geology and ground-water hydrology of the Mokelumne area, California

The Mokelumne River basin of central California comprises portions of the California Trough and the Sierra Nevada section of the Pacific Mountain system. The California Trough is divisible into four subsections-the Delta tidal plain, the Victor alluvial plain, tlie river flood plains and channels, and the Arroyo Seco dissected pediment. These four subsections comprise the land forms produced by the Mokelumne River and other streams since the Sierra Nevada attained its present height in the Pleistocene epoch. The Victor alluvial plain rises eastward from the Delta plain and abuts on the dissected Arroyo Seco pediment; in the Mokelumne area it is 12 to 16 miles wide and slopes between 5 and 8 feet in a mile. It includes relatively extensive tracts that are intensively cultivated and irrigated with water pumped from wells. The Victor plain has been compounded of overlapping alluvial fans along the western base of the Sierra Nevada. It is prolonged eastward into the pediment by tongues of alluvium along several of the present streams; thus it seems likely that the present stream pattern in the eastern part of the area has been fixed since dissection of the pediment began. Three of the four major streams-the Mokelumne and Cosumnes Rivers and Dry Creek-traverse the Victor plain in trenches which are 15 to 40 feet deep at the heads of their respective alluvial fans but which die out toward the west. The floors of these trenches, the historic flood plains, are from 100 yards to a mile wide. The exceptional major stream, which has not entrenched itself, is the Calaveras River. The Arroyo Seco pediment, which lies east of the Victor plain, was initially at least 8 to 15 miles wide and lay along the western foot of the Sierra Nevada entirely .across the Mokelumne area. Its numerous remnants decline 15 to 35 feet in a mile toward the west. The Sierra Nevada section adjoins and lies east of the California Trough. Its major ridge crests define a volcanic plain whose westward slope is' inferred to have been initially about 90 feet in a mile but is now about 180 feet in a mile, owing to tilting of the Sierra Nevada block in Pleistocene time. In and near the Mokelumne area the Sierra Nevada and California Trough together are roughly coextensive with a single structural unit. The Sierra Nevada constitutes a block that has risen with respect to adjoini;ng valley areas 'by simple rotation or tilting toward the west; it has not been warped or faulted -extensively. It is inferred that this block extends westward beneath the thick alluvial deposits of the trough without material warping or faulting. The oldest rocks of the Mokelumne region are the Carboniferous and Jurassic -rocks that compose the crystalline core of the Sierra Nevada. These are overlain unconformably by sediments of Tertiary age--in upward succession the lone, Valley Springs, Mehrten, and J.Jaguna formations. Of these formations all except the lone are newly discriminated, and type sections are described in the full text. These Tertiary sediments form a great wedge, thinnest along the mountain front to the east, where they have been truncated by erosion. They dip about 2° W. The lone formation (Eocene) consists chiefly of sandstone, clay, and shale; its maximum thickness is 450 feet. The Valley Springs formation (middle? Miocene) overlies the lone formation unconformably. It is composed largely of greenish-gray clay, shale, and sandstone derived from rhyolitic ejectamenta. These rhyolitic deposits are confined to narrow channels in the higher part of the Sierra Nevada, but they spread fanlike over the lower western edge of the mountain block, where they attain a maximum thickness of 525 feet. The Mehrten formation (upper? Miocene and lower Pliocene?) comprises the andesitic rocks that constructed the Sierran volcanic plain. In the Mokelumne area it consists chiefly of sandstone and siltstone but includes, as a minor though conspicuous part of the formation, layers and tongues of resistant breccia or agglomerate, which are presumed to have originated as mud flows. Nonfragmental andesite is not known to occur in the Mokelumne area, although several possible vents occur farther east. In the eastern part of the area the Mehrten formation truncates in turn the Valley Springs and lone formations and the pre-Cretaceous rocks; in the western part the Mehrten formation (andesitic) interfingers with the underlying Valley Springs formation (rhyolitic). Its maximum measured thickness is 400 feet. Few of the irrigation wells are so deep that they can be said with assurance to reach the Mehrten formation. The Laguna formation (Pliocene? and possibly lower Pleistocene) comprises poorly sorted, nonandesitic fluviatile sedimentary that overlie the Mehrten formation. It is inferred to be essentially parallel to and tilted equally with the Mehrten formation and to be about 400 feet thick. The Arroyo Seco gravel (presumably middle Pleistocene) veneers the Arroyo Seco pediment. At its easternmost outcrops the formation is composed of pebbles, cobbles, and boulders in a matrix of brick-red sand and silt; farther west, down the slope of the pediment, it becomes pr9gressively finer. It is inferred that the Arroyo Seco gravel is a coarse fraction of the rock waste that was transported from the Sierra Nevada after the Sierran.block was tilted in Pleistocene time. It is inferred further that the correlative of the Arroyo Seco gravel in the California Trough is a wedge-shaped mass of sediments whose base is the tilted Laguna formation and whose top can be interpolated by projecting a hypothetical surface through the remnants of the pediment. The Victor formation comprises the fluviatile sand, silt, and gravel that built the Victor alluvial plain over the hypothetical equivalent of the Arroyo' Seco gravel along the axis of the California Trough and against the western front of the dissected pediment to the east. The formation is thought to be about 100 feet thick along the western margin of the Mokelumne area, according to an estimate based upon projecting the slope of the Arroyo Seco pediment westward beneath the Victor plain. The Mokelumne area lies on the fertile central plain along the Mokelumne River about the city of Lodi, in northern San Joaquin County, and has been intensively developed for the cultivation of grapes, deciduous fruits, and other crops. Of necessity its great productiveness is maintained by irrigation. Extensive irrigation from wells began about 1907 and has increased steadily until in 1932 about 50,000 acres (80 percent of the area) was watered in that manner. The specific question at issue is the extent to which the supply of ground water and hence the productiveness of the area are dependent upon the water flowing in the Mokelumne River and the extent to which that productiveness may be influenced by regulation of the stream--:in particular, by the substantial regulation of the river that is accomplished by the Pardee Dam of the East Bay Municipal Utility District, which began to function in March 1929. The depth of 1,447 irrigation wells in five townships in the central part of the area (T. 3 N., Rs. 6 and 7 E., and T. 4 N., Rs. 6 to 8 E.) ranges from 20 to 910 feet. About half the wells bottom within a 100-foot zone whose base is 75 feet below the projected Arroyo Seco pediment; essentially that zone constitutes the Victor formation. Only 6 percent of the wells bottom within the next lower 25- foot zone, but the percentage increases sharply for depths still greater; it is inferred that impervious strata are relatively persistent between 75 and 100 feet below the projected pediment and that these are the uppermost part of the Arroyo Seco gravel. Of 580 observation wells known to bottom in the Victor formation, essentially all appear to indicate a regional water-table stage; thus the water is essentially unconfined. On the other hand, nearly all wells so deep that they reach the Arroyo Seco gravel or some underlying formation tap confined water. Near the Mokelumne River the water levels in these deep wells stand below the water table, which is semiperched. In most deep wells remote from the river the water level stands above the water table except during the pumping season. Fluctuations of ground-water levels are ascribed to moving or changing load on the land surface, earthquakes, variation of barometic pressure, ground-water draft by vegetation, infiltration of rain and certain indirect effects of rainfall, infiltration of water applied to the land for irrigation, variation in the discharge of streams, and pumping from wells. In the eastern part of the central district, between Clements and the vicinity of Lockeford, it is inferred that (1) the river and the water in the alluvium of the flood plain are not insulated from the water in the sediments that form the adjacent Victor plain; (2) locally if not generally, however, there are discontinuities in pervious strata along the outer margin of the flood plain, where the water table passes from the alluvium into the enclosing sediments, so that percolation of ground water is impeded materially at that margin; (3) rising river stages set up ground-water waves that store relatively large volumes of water in the alluvium close to the river, whereas falling stages cause much of that stored water to percolate back into the river, weeks and even months lapsing before the ground-water stage becomes steady within the flood plain; and (4) seepage loss from the river into the alluvium tends to be intermittent and to alternate with seepage gain, the rate of loss or gain lagging weeks or months behind the fluctuations of river stage and lagging more for moderate changes at low stage. However, in the succeeding reach downstream as far as Woodbridge, it is inferred that percolation of ground water is not impeded generally along the outer margin of the flood plain and that the river tends to lose almost continuously by seepage rather than intermittently, although the rate of loss fluctuates somewhat in response to changing river stage. The yearly pumpage for irrigation has been as much as 114,600 acre-feet (1928- 29), and there have been as many as 2,500 wells equipped with irrigation pumping plants (1931). Commonly the wells are pumped only in daylight and are idle over week-ends and holidays, also during and after protracted rainstorms in the early part of the season. In a small district near Victor pumping in recent years has begun in January or February, has reached its height in March, and largely has passed by April. In outlying districts general pumping has begun as late as May, reached its height in June or July, and waned by September. Since 1907 the water table appears to have declined steadily in most of the Mokelumne area except along the river. The decline was least in the Woodbridge Irrigation District, where in four typical wells. the average decline from 1907 to 1937 was 3 feet, or 0.15 foot a year. Among 18 shallow wells in the district of most intensive pumping the average recession of the water table from 1907 to 1927 was 11 feet, or 0.55 foot a year; the greatest measured recession was 15 feet, or 0.75 foot a year. From 1927 to 1933 the water table declined 5 feet or more over most of the central pumping district except within 2 miles of the Mokelumne River, and the greatest measured decline was 9 feet. The area of material recession ,extends 4 to 7 miles eastward beyond the central pumping district, whence it is inferred that pumping has drawn gradually on remote ground-water storage. It is inferred that the Mokelumne River ordinarily has been a losing stream between the Mehrten dam site, near Clements, and the Woodbridge Dam, the area that received the percolate having been triangular with its upstream and having included about 5,200 acres of the flood plain and 36,500 acres in outlying districts to the north and to the south. Mean fluctuations of the water table within the area receiving percolate from ihe river are believed to indicate that relatively little water is drawn from outside the area. Accordingly, simple storage methods are competent for a ground-water inventory. It is inferred that the rate of seepage loss from the river depends jointly upon river discharge, stage in the Woodbridge Reservoir, and groundwater pumpage. The foregoing inferences lead to the following conclusions with respect to ground-water replenishment by seepage loss from the river in the intensively cultivated district about Lodi: (1) The annual replenishment has tended to increase for at least two decades, owing to the gradual increase in head between surface water and ground water as ground-water levels have been lowered progressively by pumping; (2) annual replenishment has tended to increase, especially in recent years, owing to gradually prolonged use of the Woodbridge Reservoir, for thereby a relatively large wetted area and great differential head have been maintained for an increasing term; (3) the rate of replenishment tends to be greater under regulation than under the so-called natural regimen, to the extent that regulation has maintained a moderately large wetted area and stage in the river through the later part ·of each pumping season, whi1e the ground-water levels have been lowest. Moreover, for any particular yearly run-off below the Mehrten dam site, the replenishment by seepage would tend to be greater under the regulated regimen to the extent that fluctuations in discharge were suppressed, for the greatest yearly mean stage and mean wetted area would be afforded by constant discharge. -Thus, diverting water out of the Mokelumne River Basin at the Pardee Dam does not necessari1y-entail a diminution in ground-water replenishment by seepage loss along the lower reach of the stream, at least in the replenishment beneath the Victor plain above the gaging station at Woodbridge. Rather. the Pardee Dam affords a means for so regulating the discharge as to effect a maximum ground-water replenishment with-a given run-off in the natural channel. Bodies of ground water perched above the regional water table are common in the Laguna formation, especially in its lower part. Conspicuous bodies occur about 3 miles south of Clay, in a district between 1 mile and 5 miles south of Clements, and along Dry Creek in T. 5 N., Rs. 7 and 8 E. From the relation between the water table and the piezometric surface for water confined in deep aquifers, the area receiving percolate from the Mokelumne River may be divided roughly into (1) a central area, extending not :p1ore than half a mi1e beyond the flood plain, in which the piezometric surface is inferred to have stood below the water table throughout the term of the investigation and hence in which the difference in head has favored the percolation of water from shallow beds into deep beds in all seasons, and (2) an outlying area in which the difference in head likewise favors downward percolation into deep beds during the pumping season but favors upward percolation during the nonpumping season. This outlying area includes about 75 percent of the segment of the Victor plain that receives percolate from the river. From 1927 to 1933 the subartesian head that existed during the nonpumping season in the area remote from the river tended to increase; it is therefore inferred that the relative opportunity for seasonal recharge of the shallow water-bearing beds by underfeeding has likewise tended to increase. On the other hand, the negative differential head in wells near the river also has tended to increase; thus in this central area the opportunity for discharge of water from shallow beds by downward percolation has probably tended to increase. It is believed that ground-water storage within the area near the river is not decreased materially by" discharge westward through deep pervious beds, also that the yearly addition to ground-water storage in the outlying area by deep percolation from a remote easterly source is scant and for all practical purposes is offset by downward percolation along the river.

California↗

Prediction of the probability of elevated nitrate concentrations at groundwater depths used for drinking-water supply in the Puget Sound basin, Washington, 2004–19

The Puget Sound basin encompasses the 13,700-square-mile area that drains to the Puget Sound and the adjacent marine waters of Washington State. Well more than 4 million people live within the basin, with numbers continuing to increase, who rely on the basin’s natural resources including groundwater. The Puget Sound Partnership was created by a Washington State statute to implement a science-based recovery of the Puget Sound to help address impacts to these resources. As part of the recovery, the partnership developed the Puget Sound Vital Signs as measures of ecosystem health that guide the assessment of progress toward Puget Sound recovery goals. The Puget Sound Partnership Leadership Council adopted a Drinking Water Vital Sign associated with human health and quality of life, recognizing certain indicators as integral to the sustainability of Puget Sound recovery efforts. One such Vital Sign indicator was the vulnerability of groundwater throughout the aquifers of the Puget Sound basin to elevated nitrate concentrations as defined by the probability of exceeding 2 milligrams/liter (mg/L) at a specific location and well depth. The U.S. Geological Survey (USGS) led the effort to characterize groundwater vulnerability. For this study, groundwater vulnerability refers to a probability with which a contaminant applied at or near the land surface can migrate to the aquifer of interest for a given set of land-use practices. Nitrate concentration data were selected for evaluation because elevated nitrate concentrations are typically caused by anthropogenic activities and have been associated with deleterious impacts on human health. To identify groundwater vulnerability to elevated nitrate concentrations, logistic regression was used to relate anthropogenic (human associated) and natural variables to the occurrence of elevated nitrate concentrations in untreated groundwater from large public water supply system wells found within the Washington State Department of Health Sentry database. Variables that were analyzed included well depth, soil hydraulic conductivity, precipitation, population density, fertilizer application amounts, and land-use types. Statistically significant models that predicted the probabilities of groundwater nitrate concentrations greater than 2 mg/L based on the predictor variables were created for the time periods 2000–04, 2005–09, 2010–14, and 2015–19. For all time periods, well depth and a measure of the abundance of urban and agricultural land over or near the well consistently helped explain the vulnerability of the well to elevated nitrate concentrations defined as a probability of exceeding 2 mg/L of nitrate. Precipitation and (or) soil hydraulic conductivity were also important predictor variables in the models. The models for each time period were used to create maps of groundwater vulnerability at 150- and 300-foot depths throughout the Puget Sound basin. As expected, the most vulnerable locations were associated with shallower well depths and increased agriculture and urban land cover. Across all four time periods, groundwater vulnerability throughout the Puget Sound was low, with probabilities of exceeding 2 mg/L concentrations of nitrate at depths at 150 and 300 feet typically less than 50 percent. Results also found a slight decrease in probabilities of elevated nitrate concentrations throughout the basin over time. More specifically, additional statistical tests found that groundwater with probabilities of less than about 60 percent declined from 2000 to 2019 and represented more than 75 percent of the modeled Puget Sound basin aquifer. Wells with greater than 60 percent probability increased over the same time period but represented only about 25 percent of the aquifer. The maps and statistical analysis presented in the study provide valuable and informative evaluation of the vulnerability of groundwater in the Puget Sound basin to elevated nitrate concentrations. The probability maps do not represent measured nitrate concentrations in groundwater, but rather they present the probability that nitrate concentrations exceed 2 mg/L. The models and predictions from this study are a viable indicator for the Puget Sound Partnership’s Healthy Human Population—Drinking Water Vital Sign. The logistic regression modeling approach presented here benefits water managers by allowing them to assess temporal trends in a range of probabilities, explore vulnerability changes as new regional land cover and anthropogenic data are generated, and distinguish vulnerabilities at different depths within the aquifer.

Washington↗

Native Prairie Adaptive Management: a multi region adaptive approach to invasive plant management on Fish and Wildlife Service owned native prairies

Much of the native prairie managed by the U.S. Fish and Wildlife Service (FWS) in the Prairie Pothole Region (PPR) of the northern Great Plains is extensively invaded by the introduced cool-season grasses, smooth brome ( Bromus inermis ) and Kentucky bluegrass ( Poa pratensis ). Management to suppress these invasive plants has had poor to inconsistent success. The central challenge to managers is selecting appropriate management actions in the face of biological and environmental uncertainties. In partnership with the FWS, the U.S. Geological Survey (USGS) developed an adaptive decision support framework to assist managers in selecting management actions under uncertainty and maximizing learning from management outcomes. This joint partnership is known as the Native Prairie Adaptive Management (NPAM) initiative. The NPAM decision framework is built around practical constraints faced by FWS refuge managers and includes identification of the management objective and strategies, analysis of uncertainty and construction of competing decision models, monitoring, and mechanisms for model feedback and decision selection. Nineteen FWS field stations, spanning four states of the PPR, have participated in the initiative. These FWS cooperators share a common management objective, available management strategies, and biological uncertainties. Though the scope is broad, the initiative interfaces with individual land managers who provide site-specific information and receive updated decision guidance that incorporates understanding gained from the collective experience of all cooperators. We describe the technical components of this approach, how the components integrate and inform each other, how data feedback from individual cooperators serves to reduce uncertainty across the whole region, and how a successful adaptive management project is coordinated and maintained on a large scale. During an initial scoping workshop, FWS cooperators developed a consensus management objective: increase the composition of native grasses and forbs on native sod while minimizing cost. Cooperators agreed that decision guidance should be provided annually and should account for local, real-time vegetation conditions observed on the ground. Over the course of development, two prototypes of the decision framework were considered. The final framework recognized four alternative actions that managers could take in any given year: (1) Graze—targeted use of grazing ungulates to achieve defoliation, (2) Burn—application of prescribed fire as the single form of defoliation, (3) Burn/Graze—a combination treatment, and (4) Rest—no action. The study area included northern mixed-grass and tallgrass prairie. Native vegetation in mixed–grass prairie has a strong cool-season component and thus the dominant native species have a phenology similar to that of smooth brome and Kentucky bluegrass, making management of those species challenging. In contrast, tallgrass prairie has a strong warm-season native component, leading to an existence of cool-season windows, periods of time in the fall and spring when cool‐season invasive grass species are actively growing and vulnerable to damage via select management actions, but warm‐season grass species are not active and are thus less susceptible to damage via the same actions. This dichotomy between prairie types necessitated the development of separate but parallel decision support systems for mixed-grass and tallgrass biomes. Management units are parcels of native prairie that receive a single management treatment at any one time over their entire extent. At any particular time, the vegetation state of each management unit is characterized by the amount of cover of native grasses and forbs and the type of invasive grass that is dominant. In addition, each unit has a defoliation state which reflects the number of years since the last defoliation event and an index to how intensively the unit was managed during the previous 7 years. State-transition models are used to predict the state of a management unit in year t+1 from its state in year t and a prescribed management action that was applied between the two monitoring events. Alternative models are built around key uncertainties that make choice of a management action difficult. Three uncertainties revolve around whether the effect of management actions depends on (1) type of dominant invader, (2) past defoliation history, and (3) level of invasion. Two additional uncertainties are considered when choosing a management action for tallgrass units: (4) the effectiveness of grazing within the cool-season window as a surrogate for burning when smooth brome is the dominant invader, and (5) the differential effect of active management outside the window as compared to rest. Because data on the probability of transitioning from one state to another under the various models were lacking, expert opinion and elicitation were used to parameterize the models. In addition, cooperators participated in elicitation exercises to extract their beliefs regarding the value of having native prairie compared to the cost of achieving it. Quantifying the subjective expression of utility in this way allowed for mathematical representation of the management objective into an objective function. By maximizing the objective function, cumulative utility is maximized, leading to the identification of a sequence of decisions that will achieve the management objective. The NPAM system adopted a vegetation monitoring protocol that was rapid, inexpensive, and familiar to many of the cooperators. The monitoring protocol served three purposes: (1) determining current vegetation and defoliation states of each unit, (2) evaluating progress toward the management objective, and (3) assessing predictive performance of the alternative models. The management year runs from September 1 to August 31. Management can be applied anytime during that period and monitoring takes places from late June to mid-August. Cooperators enter vegetation data and management information into a centralized database by August 25 of each year. Given the current state of the system (vegetation and defoliation states) and the current understanding of the system (or the belief state), identifying the current best management decision is a matter of looking up the combination (that is, system state and belief state) in the appropriate (mixed-grass or tallgrass) optimal decision table. Given complete uncertainty at the outset of decision-making, initial assignment of equal belief weights to each model was believed reasonable. The decisions in the optimal decision table that correspond to the current belief state constitute the current optimal decision policy. By August 31 of each year, individual cooperators are provided with a recommended management action for each of their management units for the upcoming management year. Upon receiving the management recommendations for their units, managers consider the recommendation, along with other relevant information, and at some point during the year one of the management alternatives is carried out. This iterative cycle of making and implementing a management decision, predicting the response, monitoring the outcome, comparing predicted and observed outcomes, updating model weights, and recommending a management action for the next cycle is expected to result in an accumulation of weight on a representative model of system dynamics, thereby increasing understanding needed to effectively manage native prairies. The NPAM system is now entering its second full year of complete operation, and represents one of only a few fully implemented applications of adaptive management within the U.S. Fish and Wildlife Service. NPAM is truly unique in that it originated from the ground up as a result of the leadership and steadfastness of several refuge biologists and managers confronted with a common problem. These biologists recognized that working together across a large landscape presented perhaps the best opportunity for halting and reversing the invasion of native grasslands by non-native cool-season grasses. Importantly, the NPAM system encapsulates the collective thinking and experience of tens if not hundreds of individuals who have battled this vexing problem for much of their careers. The NPAM initiative is rooted in principles of adaptive management, thereby affording the opportunity for grassland managers to pursue management objectives while acquiring information to reduce uncertainty and improve future management. The project introduced a number of technical innovations that will serve as templates for conservation efforts throughout and beyond the U.S. Fish and Wildlife Service. First, NPAM is an on-the-ground implementation of active adaptive management—possibly the first of its kind in conservation management—in which recommended management actions result from a prospective analysis of future learning (Williams, 1996). Second, by the use of dynamic optimization, NPAM demonstrates how decisions can be made that take into account possible future transitions of the system. Third, NPAM demonstrates how models of partial controllability are an effective means of accommodating unpredictable circumstances that cause a manager to follow a different course than was intended. Finally, the database developed for NPAM is an unparalleled system that enables the rapid integration of data from the field for the generation of ‘just-in-time’ management recommendations. In all, NPAM provides an example of how a science-management partnership can be forged to achieve large-scale conservation objectives.

Minnesota;Montana;North Dakota;South Dakota↗

Monitoring to assess progress toward meeting the Assabet River, Massachusetts, phosphorus total maximum daily load - Aquatic macrophyte biomass and sediment-phosphorus flux

In 2004, the Total Maximum Daily Load (TMDL) for Total Phosphorus in the Assabet River, Massachusetts, was approved by the U.S. Environmental Protection Agency. The goal of the TMDL was to decrease the concentrations of the nutrient phosphorus to mitigate some of the instream ecological effects of eutrophication on the river; these effects were, for the most part, direct consequences of the excessive growth of aquatic macrophytes. The primary instrument effecting lower concentrations of phosphorus was to be strict control of phosphorus releases from four major wastewatertreatment plants in Westborough, Marlborough, Hudson, and Maynard, Massachusetts. The improvements to be achieved from implementing this control were lower concentrations of total and dissolved phosphorus in the river, a 50-percent reduction in aquatic-plant biomass, a 30-percent reduction in episodes of dissolved oxygen supersaturation, no low-flow dissolved oxygen concentrations less than 5.0 milligrams per liter, and a 90-percent reduction in sediment releases of phosphorus to the overlying water. In 2007, the U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection, initiated studies to evaluate conditions in the Assabet River prior to the upgrading of wastewater-treatment plants to remove more phosphorus from their effluents. The studies, completed in 2008, implemented a visual monitoring plan to evaluate the extent and biomass of the floating macrophyte Lemna minor (commonly known as lesser duckweed) in five impoundments and evaluated the potential for phosphorus flux from sediments in impounded and free-flowing reaches of the river. Hydrologically, the two study years 2007 and 2008 were quite different. In 2007, summer streamflows, although low, were higher than average, and in 2008, the flows were generally higher than in 2007. Visually, the effects of these streamflow differences on the distribution of Lemna were obvious. In 2007, large amounts of floating macrophytes accumulated behind bridge constrictions and dams; in 2008, high flows during the early part of the growing season carried floating macrophytes past bridges and over dams, minimizing accumulations. Samples of Lemna were collected and weighed to provide an estimate of Lemna biomass based on areal coverage during the summer growing seasons at eight sites in the five impoundments. Average estimated biomass during 2007 was approximately twice the 2008 biomass in each of the areas monitored. In 2007, in situ hyperspectral and high-resolution, multispectral data from the IKONOS satellite were obtained to evaluate the feasibility of using remote sensing to monitor the extent of aquatic plant growth in Assabet River impoundments. Three vegetation indices based on light reflectance were used to develop metrics with which the hyperspectral and satellite data were compared. The results of the comparisons confirmed that the high-resolution satellite imagery could differentiate among the common aquatic-plant associations found in the impoundments. The use of satellite imagery could counterbalance emphasis on the subjective judgment of a human observer, and airborne hyperspectral data can provide higher resolution imagery than multispectral satellite data. In 2007 and 2008, the potential for sediment flux of phosphorus was examined in free-flowing reaches of the river and in the two largest impoundments-Hudson and Ben Smith. These studies were undertaken to determine in situ flux rates prior to the implementation of the Assabet River Total Maximum Daily Load (TMDL) for phosphorus and to compare these rates with those used in the development and evaluation of the TMDL. Water samples collected from a chamber placed on the river bottom were analyzed for total phosphorus and orthophosphorus. Ambient dissolved oxygen concentrations and seasonal temperature differences appeared to affect the rates of sequestration and sediment release of phosphorus. When dissolved oxygen concentrations remained relatively high in the chambers and when the temperature was relatively low, the tendency was for phosphorus concentrations to decrease in the chambers, indicating sediment sequestration of phosphorus; when dissolved oxygen concentrations dropped to near zero and temperatures were warmest, phosphorus concentrations increased in the chambers, indicating phosphorus flux from the sediment. The rates of release and sequestration in the in situ studies were generally comparable with the rates determined in laboratory studies of Assabet River sediment cores for State and Federal agencies. Sediment-core and chamber studies produced substantial sediment fluxes to the water column only under extremely low-DO or anaerobic conditions rarely found in the Assabet River impoundments; thus, sediment is not likely to be a major phosphorus source, especially when compared to the wastewater effluent, which sustains higher ambient concentrations. The regulatory agencies now (2011) have substantial laboratory and field data with which to determine the required 90-percent reduction in phosphorus flux after the completion of upgrades to the wastewater-treatment plants that discharge to the Assabet River.

Massachusetts↗

Evolving magma storage conditions beneath Mount St. Helens inferred from chemical variations in melt inclusions from the 1980-1986 and current (2004-2006) eruptions

Major element, trace element, and volatile concentrations in 187 glassy melt inclusions and 25 groundmass glasses from the 1980-86 eruption of Mount St. Helens are presented, together with 103 analyses of touching FE-Ti oxide pairs from the same samples. These data are used to evaluate the temporal evolution of the magmatic plumbing system beneath the volcano during 1980-86 and so provide a framework in which to interpret analyses of melt inclusions from the current (2004-2006) eruption. Major and trace element concentrations of all melt inclusions lie at the high SiO 2 end of the data array defined by eruptive products of the late Quaternary age from Mount St. Helens. For several major and trace elements, the glasses define a trend that is oblique to the whole-rock trend, indicating that different mineral assemblages were responsible for the two trends. The whole-rock trend can be ascribed to differentiation of hydrous basaltic parents in a deep-seated magma reservoir, probably at depths great enough to stabilize garnet. In contrast, the glass trends were generated by closed-system crystallization of the phenocryst and microlite mineral assemblages at low pressures. The dissolved H 2 O content of the melt inclusions from 1980-86, as measured by the ion microprobe, ranges from 0 to 6.7 wt. percent, with the highest values obtained from the plinian phase of May 18, 1980. Water contents decrease with increasing SiO 2 , consistent with decompression-driven crystallization. Preliminary data for dissolved CO 2 in melt inclusions from the May 18 plinian phase from August 7, 1980, indicate that X H 2 O in a vapor phase was approximately constant at 0.80, irrespective of H 2 O content, suggestive of closed-system degassing with a high bubble fraction or gas streaming through the subvolcanic system. Temperature and f O 2 estimates for touching Fe-Ti oxides show evidence for heating during crystallization owing to release of latent heat. Consequently, magmas with the highest microlite crystallinities record the highest temperatures. Magmas also become progressively reduced during ascent and degassing, probably as a result of redox equilibria between exsolving S-bearing gases and magmas. The lowest temperature oxides have f O 2 ≈ NNO, similar to high-temperature fumarole gases from the volcano. The temperature and f O 2 of the magma tapped by the plinian phase of May 18, 1980, are 870-875°C and NNO+0.8, respectively. The dissolved volatile contents of the melt inclusions have been used to calculate sealing pressures; that is, the pressure at which chemical exchange between inclusion and matrix melt ceased. These are greatest for the May 18 plinian magma (120 to 320 MPa); lower pressures are recorded by samples of the preplinian cryptodome and by all post-May 18 magmas. Magma crystallinity, calculated from melt-inclusion Rb contents, is negatively correlated with sealing pressure, consistent with decompression crystallization. Elevated contents of Li in melt inclusions from the cryptodome and post-May 18 samples are consistent with transfer of Li in a magmatic vapor phase from deeper parts of the magma system to magma stored at shallower levels. The Li enrichment attains its maximum extent at ~150 MPa, which is ascribed to separation of a single vapor phase into H 2 O-rich gas and dense Li-rich brine at the top of the magma column. There are striking correlations between melt-inclusion chemistry and monitoring data for the 1980-86 eruption. Dissolved SO 2 contents of melt inclusions from any given event, multiplied by the mass of magma erupted during that event, correlate with the measured flux of SO 2 at the surface, suggesting that magma degassing and melt-inclusion sealing are closely related in time and space. Textural and chemical evidence indicates that melt inclusions became effectively sealed (physically or kinetically) shortly before eruption. Thus by converting pressure to depth using a density model and edifice-loading algorithm for the volcano, changing depths of magma extraction with time can be tracked and compared to the seismic record. The plinian eruption of May 18, 1980, involved magma stored 5-11 km below sea level; this is inferred to be the subvolcanic magma chamber. The preceding eruptions, including the May 18, 1980, blast, involved magma withdrawal from the cryptodome and conduit down to 5 km below sea level. Subsequent 1980 eruptions tapped magma down to depths of ≤10 km below sea level. Tapping of magma stored deeper than 2 km below sea level stopped abruptly at the end of 1980, coincident with the onset of extensive shallow seismicity and a change from explosive to effusive eruption style from 1981 to 1986. Overall, the 1980-86 eruption is consistent with the evisceration of a thin, vertically extensive body of magma extending from 5 to at least 11 km below sea level and connected to the surface by a thin conduit. In the absence of sustained high magma-supply rates from depth, decompression crystallization of magma ascending through the system leads eventually to plugging of the conduit. The current eruption of Mount St. Helens shares some similarities with the 1981-86 dome-building phase of the previous eruption, in that there is extensive shallow seismicity and extrusion of highly crystalline material in the form of a sequence of flows and spines. Melt inclusions from the current eruption have low H 2 O contents, consistent with magma extraction from shallow depths. Highly enriched Li in melt inclusions suggests that vapor transport of Li is a characteristic feature of Mount St. Helens. Melt inclusions from the current eruption have subtly different trace-element chemistry from all but one of the 1980-86 melt inclusions, with steeper rareearth-element (REE) patterns and low U, Th, and high-fieldstrength elements (HFSE), indicating addition of a new melt component to the magma system. It is anticipated that increasing involvement of the new melt component will be evident as the current eruption proceeds.

Washington↗

Assessment of the hydraulic connection between ground water and the Peace River, west-central Florida

The hydraulic connection between the Peace River and the underlying aquifers along the length of the Peace River from Bartow to Arcadia was assessed to evaluate flow exchanges between these hydrologic systems. Methods included an evaluation of hydrologic and geologic records and seismic-reflection profiles, seepage investigations, and thermal infrared imagery interpretation. Along the upper Peace River, a progressive long-term decline in streamflow has occurred since 1931 due to a lowering of the potentiometric surface of the Upper Floridan aquifer by as much as 60 feet because of intensive ground-water withdrawals for phosphate mining and agriculture. Another effect from lowering the potentiometric surface has been the cessation of flow at several springs located near and within the Peace River channel, including Kissengen Spring, that once averaged a flow of about 19 million gallons a day. The lowering of ground-water head resulted in flow reversals at locations where streamflow enters sinkholes along the streambed and floodplain. Hydrogeologic conditions along the Peace River vary from Bartow to Arcadia. Three distinctive hydrogeologic areas along the Peace River were delineated: (1) the upper Peace River near Bartow, where ground-water recharge occurs; (2) the middle Peace River near Bowling Green, where reversals of hydraulic gradients occur; and (3) the lower Peace River near Arcadia, where ground-water discharge occurs. Seismic-reflection data were used to identify geologic features that could serve as potential conduits for surface-water and ground-water exchange. Depending on the hydrologic regime, this exchange could be recharge of surface water into the aquifer system or discharge of ground water into the stream channel. Geologic features that would provide pathways for water movement were identified in the seismic record; they varied from buried irregular surfaces to large-scale subsidence flexures and vertical fractures or enlarged solution conduits. Generally, the upper Peace River is characterized by a shallow, buried irregular top of rock, numerous observed sinkholes, and subsidence depressions. The downward head gradient provides potential for the Peace River to lose water to the ground-water system. Along the middle Peace River area, head gradients alternate between downward and upward, creating both recharging and discharging ground-water conditions. Seismic records show that buried, laterally continuous reflectors in the lower Peace River pinch out in the middle Peace River streambed. Small springs have been observed along the streambed where these units pinch out. This area corresponds to the region where highest ground-water seepage volumes were measured during this study. Further south, along the lower Peace River, upward head gradients provide conditions for ground-water discharge into the Peace River. Generally, confinement between the surficial aquifer and the confined ground-water systems in this area is better than to the north. However, localized avenues for surface-water and ground-water interactions may exist along discontinuities observed in seismic reflectors associated with large-scale flexures or subsidence features. Ground-water seepage gains or losses along the Peace River were quantified by making three seepage runs during periods of: (1) low base flow, (2) high base flow, and (3) high flow. Low and high base-flow seepage runs were performed along a 74-mile length of the Peace River, between Bartow and Nocatee. Maximum losses of 17.3 cubic feet per second (11.2 million gallons per day) were measured along a 3.2-mile reach of the upper Peace River. The high-flow seepage run was conducted to quantify losses in the Peace River channel and floodplain between Bartow and Fort Meade. Seepage losses calculated during high-flow along a 7.2-mile reach of the Peace River, from the Clear Springs Mine bridge to the Mobil Mine bridge, were approximately 10 percent of the river flow, or 118 cubic feet per second. Calculated seepages along the Peace River in Hardee and De Soto Counties were inconclusive, because most seepages were within the range of discharge measurement error. Two continuous aerial thermal infrared imagery surveys were conducted to locate sites of ground-water discharge along the Peace River. Although temperature and hydrologic conditions were ideal to observe spring flow using thermal infrared imaging techniques, no sources of ground-water discharge were identified using this method. Diffuse ground-water seepage may, however, provide significant ground-water discharge.

Florida↗

Sediment distribution and hydrologic conditions of the Potomac aquifer in Virginia and parts of Maryland and North Carolina

Sediments of the heavily used Potomac aquifer broadly contrast across major structural features of the Atlantic Coastal Plain Physiographic Province in eastern Virginia and adjacent parts of Maryland and North Carolina. Thicknesses and relative dominance of the highly interbedded fluvial sediments vary regionally. Vertical intervals in boreholes of coarse-grained sediment commonly targeted for completion of water-supply wells are thickest and most widespread across the central and southern parts of the Virginia Coastal Plain. Designated as the Norfolk arch depositional subarea, the entire sediment thickness here functions hydraulically as a single interconnected aquifer. By contrast, coarse-grained sediment intervals are thinner and less widespread across the northern part of the Virginia Coastal Plain and into southern Maryland, designated as the Salisbury embayment depositional subarea. Fine-grained intervals that are generally avoided for completion of water-supply wells are increasingly thick and widespread northward. Fine-grained intervals collectively as thick as several hundred feet comprise two continuous confining units that hydraulically separate three vertically spaced subaquifers. The subaquifers are continuous northward but merge southward into the single undivided Potomac aquifer. Lastly, far southeastern Virginia and northeastern North Carolina are designated as the Albemarle embayment depositional subarea, where both coarse- and fine-grained intervals are of only moderate thickness. The entire sediment thickness functions hydraulically as a single interconnected aquifer. A substantial hydrologic separation from overlying aquifers is imposed by the upper Cenomanian confining unit. Potomac aquifer sediments were deposited by a fluvial depositional complex spanning the Virginia Coastal Plain approximately 100 to 145 million years ago. Westward, persistently uplifted granite and gneiss source rocks sustained a supply of coarse-grained sand and gravel. Immature, high-gradient braided streams deposited longitudinal bars and channel fills across the Norfolk arch subarea. By contrast, across the Salisbury and Albemarle embayment subareas, mature, medium- to low-gradient meandering streams deposited medium- to coarse-grained channel fills and point bars segregated from fine-grained overbank deposits. The Virginia depositional complex merged northward across the Salisbury embayment subarea with another complex in Maryland. Here, additional sediments were received from schist source rocks that underwent three cycles of initial uplift and rapid erosion followed by crustal stability and erosional leveling. Because of the predominance of coarse-grained sediments, transmissivity, hydraulic conductivity, and regional velocities of lateral flow through the Potomac aquifer are greatest across the Norfolk arch depositional subarea, but decrease progressively northward with increasingly fine-grained sediments. Confining units hydraulically separate the Potomac aquifer from overlying aquifers, as indicated by large vertical hydraulic gradients. By contrast, most of the Potomac aquifer internally functions hydraulically as a single interconnected aquifer, as indicated by uniformly small vertical gradients. Most fine-grained sediments within the aquifer do not hydraulically separate overlying and underlying coarse-grained sediments. Across the Salisbury embayment depositional subarea, however, hydraulic separation among the vertically spaced subaquifers is imposed by the intervening confining units. The Potomac aquifer is the largest and most heavily used source of groundwater in the Virginia Coastal Plain. Water-level declines as great as 200 feet create the potential for saltwater intrusion. Conventional stratigraphic correlation has been generally ineffective at accurately characterizing complexly distributed fluvial sediments that compose the Potomac aquifer. Consequently, the aquifer’s internal hydraulic connectivity and overall hydrologic function have not been well understood. Water-supply planning and development efforts have been hampered, and interpretations of regulatory criteria for allowable water-level declines have been ambiguous. An investigation undertaken during 2010–11 by the U.S. Geological Survey, in cooperation with the Virginia Department of Environmental Quality, provides a comprehensive regional description of the spatial distribution of Potomac aquifer sediments and their relation to hydrologic conditions. Altitudes and thicknesses of 2,725 vertical sediment intervals represent the spatial distribution of Potomac aquifer sediments in the Virginia Coastal Plain and adjacent parts of Maryland and North Carolina. Sediment intervals are designated as either dominantly coarse or fine grained and were determined by interpretation of geophysical logs and ancillary information from 456 boreholes. Sediment-interval and borehole summary statistical data indicate regional trends in sediment lithology and stratigraphic continuity, upon which three structurally based and hydrologically distinct sediment depositional subareas are designated. Broad patterns of sediment deposition over time are inferred from published sediment pollen-age data. Discrepancies in previously drawn hydrostratigraphic relations between southeastern Virginia and northeastern North Carolina are partly resolved based on borehole geophysical logs and a recently documented geologic map and corehole. A conceptual model theorizes the depositional history of the sediments and geologically accounts for their distribution. Documented pumping tests of the Potomac aquifer at 197 locations produced 336 values of transmissivity and 127 values of storativity. Based on effective aquifer thicknesses, 296 values of sediment hydraulic conductivity and 113 values of sediment specific storage are calculated. Vertical hydraulic gradients are calculated from 9,479 pairs of water levels measured between November 17, 1953, and October 4, 2011, in 129 closely spaced pairs of wells. Borehole sediment-interval and related data provide a means to achieve high yielding production wells in the Potomac aquifer by site-specific targeting of drilling operations toward water-bearing coarse-grained sand and gravel. Advance knowledge of the potential of different parts of the aquifer also aids in planning optimal groundwater-development areas. Depositional subareas further provide a possible context for resource management. Current (2013) regulatory limits on water-level declines are relative to top surfaces of subdivided upper, middle, and lower Potomac aquifers across the entire Virginia Coastal Plain, but have the potential to exceed the same limit relative to a single undivided Potomac aquifer. By contrast, designation of the sediments as a single aquifer in the Norfolk arch and Albemarle embayment subareas—and as a series of vertically spaced subaquifers and intervening confining units in the Salisbury embayment subarea—best reflects understanding of the Potomac aquifer and can avoid the potential for excessive water-level declines. Simulation modeling to evaluate effects of groundwater withdrawals could be designed similarly, including vertical discretization and (or) zonation of the Potomac aquifer based on depositional subareas and a geostatistical distribution of aquifer properties derived from borehole sediment-interval data. Further resource-management information needs extend beyond the developed part of the Potomac aquifer, particularly across the Northern Neck and Middle Peninsula where only the shallowest part of the aquifer is known, and include structural aspects such as faults, basement bedrock, and the Chesapeake Bay impact crater.

Maryland, North Carolina, Virginia↗

The Yellowstone hotspot, Greater Yellowstone ecosystem, and human geography

Active geologic processes associated with the Yellowstone hotspot are fundamental in shaping the landscapes of the greater Yellowstone ecosystem (GYE), a high volcanic plateau flanked by a crescent of still higher mountainous terrain. The processes associated with the Yellowstone hotspot are volcanism, faulting, and uplift and are observed in the geology at the surface. We attribute the driving forces responsible for the northeastward progression of these processes to a thermal plume rising through the Earth’s mantle into the base of the southwest-moving North American plate. This progression began 16 million years ago (Ma) near the Nevada-Oregon border and arrived at Yellowstone about 2 Ma. Before arrival of the hotspot, an older landscape existed, particularly mountains created during the Laramide orogeny about 70–50 Ma and volcanic terrain formed by Absaroka andesitic volcanism mostly between 50–45 Ma. These landscapes were more muted than the present, hotspot-modified landscape because the Laramide-age mountains had worn down and an erosion surface of low relief had developed on the Absaroka volcanic terrain. The Yellowstone Plateau was built by hotspot volcanism of rhyolitic lavas and caldera-forming rhyolite tuffs (ignimbrites). Streams eroding back into the edges of this plateau have created scenic waterfalls and canyons such as the Grand Canyon of the Yellowstone and Lewis Canyon. Rhyolite is poor in plant nutrients and forms sandy, well-drained soils that support the monotonous, fire-adapted lodgepole pine forests of the Yellowstone Plateau. Non-rhyolitic rocks surround this plateau and sustain more varied vegetation, including spruce, fir, and whitebark pine forests broken by grassy meadows. Heat from the hotspot rises upward and drives Yellowstone’s famed geysers, hot springs, and mudpots. These thermal waters are home to specialized, primitive ecosystems, rich in algae and bacteria. The rock alteration associated with hydrothermal systems creates the bright colors of Yellowstone’s Grand Canyon. Basin-and-range-style faulting has accompanied migration of the hotspot to Yellowstone and formed the linear mountains and valleys that occur north and south of the hotspot track, which is the present-day eastern Snake River Plain. High rates of basin-and-range faulting occurred adjacent to the migrating Yellowstone hotspot, creating distinctive landscapes within the GYE such as the Teton Range/Jackson Hole, with characteristic rugged, forested ranges and adjacent flat-floored grassy valleys. The difference in altitude between the mountains and valleys provides a topographic gradient in which vegetation maturation advances with altitude; animal-migration patterns also follow this trend. The valleys provide natural meadows, agricultural land, town sites, and corridors for roads. Uplift of the GYE by as much as 1 km (3,000 ft) during the last 5 million years has resulted in ongoing erosion of deep, steep-walled valleys. Many prominent ecological characteristics of Yellowstone derive from this hotspot-induced uplift, including the moderate- to high- altitude terrain and associated cool temperatures and deep snowfall. Modern and Pleistocene climate and associated vegetation patterns strongly relate to the topography created by the hotspot and its track along the eastern Snake River Plain. Winter air masses from the moist northern Pacific Ocean traverse the topographic low of the Snake River Plain to where orographic rise onto the Yellowstone Plateau and adjacent mountains produces deep snow. A winter precipitation shadow forms on the lee (eastern) sides of the GYE. During Pleistocene glacial times, this moisture conduit provided by the hotspot-track-produced ice-age glaciers that covered the core of the present GYE. These glaciers sculpted bedrock and produced glacial moraines that are both forested and unforested, sand and gravel of ice-marginal streams and outwash gravels that are commonly covered with sagebrush-grassland, and silty lake sediments that are commonly covered by lush grassland such as Hayden Valley. The effects of the Yellowstone hotspot also profoundly shaped the human history in the GYE. Uplift associated with the hotspot elevates the GYE to form the Continental Divide, and streams drain radially outward like spokes from a hub. Inhabitants of the GYE 12,000–10,000 years ago, as well as more recent inhabitants, followed the seasonal green-up of plants and migrating animals up into the mountain areas. During European immigration, people settled around Yellowstone in the lower parts of the drainages and established roads, irrigation systems, and cultural associations. The core Yellowstone highland is too harsh for agriculture and inhospitable to people in the winter. Beyond this core, urban and rural communities exist in valleys and are separated by upland areas. The partitioning inhibits any physical connection of communities, which in turn complicates pursuit of common interests across the whole GYE. Settlements thus geographically isolated evolved as diverse, independent communities

Idaho, Montana, Utah, Wyoming↗

Groundwater discharge by evapotranspiration, Dixie Valley, west-central Nevada, March 2009-September 2011

With increasing population growth and land-use change, urban communities in the desert Southwest are progressively looking toward remote basins to supplement existing water supplies. Pending applications by Churchill County for groundwater appropriations from Dixie Valley, Nevada, a primarily undeveloped basin east of the Carson Desert, have prompted a reevaluation of the quantity of naturally discharging groundwater. The objective of this study was to develop a revised, independent estimate of groundwater discharge by evapotranspiration (ET g ) from Dixie Valley using a combination of eddy-covariance evapotranspiration (ET) measurements and multispectral satellite imagery. Mean annual ET g was estimated during water years 2010 and 2011 at four eddy-covariance sites. Two sites were in phreatophytic shrubland dominated by greasewood, and two sites were on a playa. Estimates of total ET and ET g were supported with vegetation cover mapping, soil physics considerations, water‑level measurements from wells, and isotopic water sourcing analyses to allow partitioning of ET g into evaporation and transpiration components. Site-based ET g estimates were scaled to the basin level by combining remotely sensed imagery with field reconnaissance. Enhanced vegetation index and brightness temperature data were compared with mapped vegetation cover to partition Dixie Valley into five discharging ET units and compute basin-scale ET g . Evapotranspiration units were defined within a delineated groundwater discharge area and were partitioned as (1) playa lake, (2) playa, (3) sparse shrubland, (4) moderate-to-dense shrubland, and (5) grassland. Groundwater ET is influenced primarily by phreatophytic vegetative cover, salinity of soil and groundwater within the playa, depth to groundwater, solar radiation, and air temperature. The annual groundwater contribution to site‑scale ET ranged from 24 to 61 percent of total ET at vegetated sites and 4 to 15 percent of total ET at playa sites. Mean annual ET g from vegetated sites ranged from 53 millimeters (mm) (0.17 foot [ft], 7.3 percent vegetative cover) to 225 mm (0.74 ft, 24.8 percent vegetative cover). Cumulative liquid‑water fluxes in the unsaturated zone indicate that ET g at vegetated sites was influenced primarily by plant transpiration. Binary mixing analyses of oxygen-18 isotopes in groundwater and shallow soil water indicate that plants predominantly use groundwater throughout the year. Groundwater fractions in greasewood stem water varied seasonally and ranged from 0.63 to 1.0. Mean annual playa ET g ranged from about 11 mm (0.04 ft) at the inner playa site (near-surface volumetric water content of 37–53 percent) to about 20 mm (0.07 ft) at the outer playa site located within 2 kilometers of the playa edge (near-surface volumetric water content of 25–38 percent), but playa ET g estimates were within the probable error (plus or minus [±] 20–23 mm; 0.06–0.08 ft). Varying playa ET g was influenced predominantly by salinity rather than depth to groundwater. Osmotic resistance and physical impediments to ET (such as surface salt crusts and salt precipitate in the soil pore space) increased with increasing salinity toward the playa center, whereas vapor pressure decreased. Mean annual basin-scale ET g totaled about 28 million cubic meters (Mm 3 ) (23,000 acre-feet [acre-ft]), and represents the sum of ET g from all ET units. Annual groundwater ET from vegetated areas totaled about 26 Mm 3 (21,000 acre-ft), and was dominated by the moderate-to-dense shrubland ET unit (54 percent), followed by sparse shrubland (37 percent) and grassland (9 percent) ET units. Senesced grasses observed in the northern most areas of the moderate-to-dense ET unit likely confounded the vegetation index and led to an overestimate of ET g for this ET unit. Therefore, mean annual ET g for moderate-to-dense shrubland presented here is likely an upper bound. Annual groundwater ET from the playa ET unit was 2.2 Mm 3 (1,800 acre-ft), whereas groundwater ET from the playa lake ET unit was 0–0.1 Mm 3 (0–100 acre-ft). Oxygen-18 and deuterium data indicate discharge from the playa center predominantly represents removal of local precipitation-derived recharge. The playa lake estimate, therefore, is considered an upper bound. Mean annual ET g estimates for Dixie Valley are assumed to represent the pre‑development, long-term ET g rates within the study area.

Nevada↗

Insights and strategic opportunities from the USGS 2024 Per- and Polyfluoroalkyl Substances (PFAS) Interagency Workshop

Introduction In 2021, the U.S. Geological Survey (USGS) published Circular 1490 titled, “Integrated Science for the Study of Perfluoroalkyl and Polyfluoroalkyl Substances (PFAS) in the Environment: A Strategic Science Vision for the U.S. Geological Survey” (Tokranov and others, 2021). Circular 1490 was created to be a resource for USGS scientists prioritizing and planning research related to per- and polyfluoroalkyl substances (PFAS) and to be a guide for developing partnerships with other scientists, State and Federal agencies, and stakeholders engaged in PFAS research and management and mitigation of the environmental and human-health effects of PFAS. This USGS PFAS Strategic Science Vision document was intended to be the foundation for a “living strategic vision,” periodically providing updates on the state of USGS PFAS research, emerging PFAS data gaps and needs, and progress on interagency and stakeholder PFAS partnerships and priorities. To meet this objective, the USGS planned to host an Interagency and Stakeholder PFAS Workshop every 2–3 years. During September 10–12, 2024, the USGS hosted the first Interagency and Stakeholder PFAS Workshop in Reston, Virginia. The Workshop brought together experts from other Federal agencies (U.S. Environmental Protection Agency, National Institute of Environmental Health Sciences, Food and Drug Administration, Department of Defense [Air Force, Army]), State agencies (Washington Fish and Wildlife, Virginia Department of Transportation), and academia (Harvard University, University of Maryland) to address key challenges relating to the measurement and modeling of PFAS and the implications for environmental health. Participants engaged in in-depth discussions centered around six pivotal topics related to PFAS: (1) sampling protocols, methods and interpretation; (2) environmental sources, source apportionment, and occurrence; (3) environmental fate and transport; (4) human and wildlife exposure routes and risk; (5) bioconcentration, bioaccumulation, and biomagnification; and (6) ecotoxicology and effects. Each topic had three breakout sessions. A recurrent theme of workshop discussions was how data on a nationwide scale for PFAS occurrence in various environmental matrices, including air, water, food crops, biota, soil, and streambed sediment could help to advance scientific understanding. Participants noted significant geospatial data gaps, particularly in the midwestern and southern United States and the Pacific Northwest. PFAS data collection tends to be more robust along the eastern seaboard and in California. Participants stressed how enhancing the integration of large and small datasets across various agencies could help to support national scale understanding of PFAS. To address these gaps, attendees suggested leveraging datasets from Federal entities like the USGS and the U.S. Department of Defense, State agencies, and municipal utility services to develop predictive contaminant detection and transport models. Improved coordination between water quality programs and USGS research could help to facilitate access to valuable data, leading to comprehensive databases that inform PFAS point (wastewater treatment plants and landfills) and nonpoint (runoff from land, atmospheric deposition, food packaging) sources, environmental transport mechanisms, environmental detection and concentrations, potential exposure routes, and health effects on different biota, including humans. A specific request was made to develop a map demarking the depth of modern (1953 or later) groundwater, which is susceptible to surface-derived anthropogenic (that is, human-made) contamination, based on tritium-age dating. Emphasis was placed on incorporation of hydrology, groundwater flow paths, groundwater–surface water interactions, and landscape factors in predictive statistical models as a step to improve contaminant source identification and tracking. Molecular fingerprinting approaches garnered attention as techniques to link specific PFAS mixtures detected in a sample to environmental sources and levels in biota (Dávila-Santiago and others, 2022). Integrating data from abiotic (that is, water, soil, and air) and biotic (that is, living organisms) systems identified as a research opportunity. For example, understanding the composition of soils and sediments, which include a mixture of mineral, plant, and animal components, could advance understanding of exposure pathways. The discussions highlighted opportunities to explore and understand the potential redistribution and biotic exposures of PFAS from biosolid and wastewater treatment plant effluent land application practices, in addition to atmospheric releases and discharges from landfill and wastewater treatment plants. Participants identified research gaps surrounding how these sources may contribute to contamination and may affect surrounding ecosystems, including a better definition of anthropogenic background concentrations. Moving forward, the collection of co-occurrence data was noted as a means to improve understanding of complex mixtures and to leverage companion modeling efforts focused on areas with high and low contamination levels to identify areas of concern and unaffected resources. Participants emphasized how centralized USGS databases and the establishment of sample-metadata archives can help to ensure that samples are preserved and accessible for future research. In conclusion, the workshop participants identified opportunities to bridge data gaps and improve measurement techniques, modeling frameworks, databases, and communication, to enhance the understanding of PFAS and their effects on environmental and human health. Upon completion of the workshop, participants indicated an interest in developing strategic data collection, modeling, and analytical approaches to address these challenges.

Open-File Report↗

Water-level monitoring along San Andreas and San Jacinto faults, southern California, during fiscal year 1980

Beginning in October 1976, a program of water-level monitoring of abandoned water wells was initiated in the Palmdale area with the purpose of identifying possible water-level changes premonitory to a major earthquake on the San Andreas fault. In October 1977, the program was extended southeastward along the rift zone to the Valyermo area. In November 1977, the monitoring of water wells along the San Jacinto fault was initiated with the expectation of experiencing a moderate size earthquake while monitoring was in progress. Currently over thirty wells are being monitored. Eleven wells are monitored continuously with Stevens Type F recorders, two of which have been modified to operate with the Caltech Remote Observatory Support Systems (TIMS). The remaining wells are probed weekly, or in some cases semi-weekly or daily, by volunteers. We are endeavoring to improve the volunteer program by increasing the frequency of measurements and simplifying the procedure to minimize measurement errors. Weekly water-level data are displayed on computer-generated hydrographs for each well. Rainfall and earthquakes are plotted on the graphs for direct comparison with water levels. The hydrographs are updated and reviewed weekly. Weekly hydrographs are also prepared from recorder charts on two wells maintained by W. R. Moyle, Jr., of the U.S. Geological Survey, Water Resources Division Office, Laguna Niguel, California. A M5.5 earthquake occurred on 25 February 1980, between the Buck Ridge and San Jacinto strands of the San Jacinto fault zone, about 8 miles east-southeast of Anza. This is the largest earthquake which has occurred on the faults within our monitoring network. The Stevens recorder chart for well number 11S/6E-3N4 located in Borrego Valley about 20 miles southeast of the epicenter indicates that, during a period of about four hours on 21 February, the water level rose 1.5 to 1.6 feet and returned to its prior level. This is one of the most remarkable short-term water-level fluctuations observed during our monitoring program; it occurred about 88 hours prior to the earthquake. This well has been monitored since October 1977 and has had a Stevens recorder since October 1978. Long-term water levels have been remarkably steady compared to those in other wells we are monitoring. This well shows a strong response to earth tides. The Stevens record on a second well in Borrego Valley (11S/6E-1C1) showed a much smaller (0.1 foot) spike in water level at the same time as the spike in well number 11S/6E-3N4. Other continuously recording wells in Anza and Ocotillo Wells showed no identifiable water-level anomalies. The two wells in Borrego Valley which showed possible strain-induced water-level spikes appear to be more sensitive to earth tides than wells in the same area which did not show spikes in water level. We cannot conceive of any nontectonic cause for the spike on the Stevens record of well 11S/6E3N4. A creep event on the Coyote Creek fault about four miles from the well is one possibility. Because the spikes in the two wells in Borrego Valley are unique for the long-term record of these wells, they may represent precursors to the 25 February earthquake. Local seismic activity was also recorded in water levels of the two wells in Borrego Springs with continuous water-level recorders. One well showed peculiar high frequency water-level fluctuations from 22 July to 21 August 1980; since then the record has appeared normal. Several wells in the Palmdale-Valyermo area have shown peculiar changes in water level within the past year or so. The long-term hydrograph of well number 5N/12W-4H1 shows the most unusual behavior. The well showed no response to the 1977 and very heavy 1978 rainstorms, yet began to rise in early 1979 and continued to rise through the dry season of 1979. By comparison, most wells show a more normal seasonal response to rainfall. Eight wells in the Palmdale-Valyermo area have shown water-level changes which are different than would have been predicted from the previous history of water-level changes and seasonal rainfall. A longer period of observation is required to determine whether the changes are anomalous. Changes in the strain pattern and other geophysical phenomena have also been observed in southern California during about the same period as the water-level changes. If the water-level changes are the result of tectonic strain, the mechanism and significance are unknown. However, it is interesting that five of the eight wells with possibly anomalous water-level changes have been identified as good strain meters based on their response to earth tides. Six wells showing an unexpected rise in water level are located west of the earthquake swarm which occurred in 1976-1977 (McNally et al, 1978), whereas the two which show water levels lower than would have been predicted are located east of the earthquake swarm.

California↗