USGS Science⌕ Search

SEARCH · USGS Science

Results for “Agricultural Systems”

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

Geology and ground-water resources of Rock County, Wisconsin

Rock County is in south-central Wisconsin adjacent to the Illinois State line. The county has an area of about 723 square miles and had a population of about 113,000 in 1957 ; it is one of the leading agricultural and industrial counties in the State. The total annual precipitation averages about 32 inches, and the mean annual temperature is about 48 ? F. Land-surface altitudes are generally between 800 and 00 feet, but range from 731 feet, where the Rock River flows into Illinois, to above 1,080 feet, at several places in the northwestern part of the county. The northern part of Rock County consists of the hills and kettles of a terminal moraine which slopes southward to a flat, undissected outwash plain. The southeastern part of the county is an area of gentle slopes, whereas the southwestern part consists of steep-sided valleys and ridges. Rock County is within the drainage basin of the Rock River, which flows southward through the center of the county. The western and southwestern parts of ,the county are drained by the Sugar River und Coon Creek, both of which flow into the Pecatonica River in Illinois and thence into the Rock River. The southeastern part of the county is drained by Turtle Creek, which also flows into Illinois before joining the Rock River. Nearly all the lakes and ponds are in the northern one-third of the county, the area of most recent glaciation. The aquifers in Rock County are of sedimentary origin and include deeply buried sandstones, shales, and dolomites of the Upper Cambrian series. This series overlies crystalline rocks of Precambrian age and supplies water to all the cities and villages in the county. The St. Peter sandstone of Ordovician age underlies all Rock County except where the formation has been removed by erosion in the Rock and Sugar River valleys, and perhaps in Coon Creek valley. The St. Peter sandstone is the principal source of water for domestic, stock, and small industrial wells in the western half of the county. This sandstone also yields some water to uncased wells that tap the deeper rocks of the Upper Cambrian series. East of the Rock River the Platteville, Decorah, and Galena formations undifferentiated, or Platteville-Galena unit, is the principal source of water for domestic and stock wells. Unconsolidated deposits of glacial origin cover most of Rock County and supply water to many small wells. In the outwash deposits along the Rock River, wells of extremely high capacity have been developed for industrial and municipal use. The most significant feature of the bedrock surface in Rock County is the ancestral Rock River valley, which has been filled with glacial outwash to a depth of at least 396 feet below the present land surface. East of the buried valley the bedrock has a fiat, relatively undissected surface. West of the valley the bedrock surface is rugged and greatly dissected. Ground water in Rock County occurs under both water-table and artesian conditions; however, because of the interconnection and close relation of all ground water in the county, the entire system is considered to be a single groundwater body whose surface may be represented by one piezometric map. Recharge occurs locally, throughout the county. Nearly all recharge is derived directly from precipitation that percolates downward to become a part of the groundwater body. Natural movement of water in the consolidated water-bearing units is generally toward the buried Rock and Sugar River valleys. Movement of water in the sandstones of Cambrian age was calculated to be about 44 million gallons a day toward the Rock River. Discharge from wells in Rock County in 1957 was about 23 million gallons a day. Nearly 90 percent of this water was drawn from the area along the Rock River. Drilled wells, most of which were drilled by the cable-tool method, range in diameter from 3 to 26 inches, and in depth from 46 to 1,225 feet. Driven wells in alluvium and glacial drift are usually 1? to 2? in

Wisconsin↗

Integrated science strategy for assessing and monitoring water availability and migratory birds for terminal lakes across the Great Basin, United States

Executive Summary In 2022, the U.S. Geological Survey (USGS) established the Saline Lake Ecosystems Integrated Water Availability Assessment (IWAAs) to monitor and assess the hydrology of terminal lakes in the Great Basin and the migratory birds and other wildlife dependent on those habitats. Scientists from across the USGS (with specialties in water quantity, water quality, limnology, avian biology, data science, landscape ecology, and science communication) formed the Saline Lake Ecosystems IWAAs Team. The team has developed this regional strategic science plan to guide data collection and assessment activities at terminal lakes in the Great Basin. The U.S. Congress requested the USGS to establish the Saline Lake Ecosystems IWAAs in response to historically low water levels at terminal lakes and associated wetlands across the Great Basin. Not all Great Basin terminal lakes have high salinity; however, all terminal lakes occur in endorheic, closed, basins with no surface-water outflow. Low lake levels across the Great Basin are the result of increased water use for agriculture and municipalities, drought conditions, and a warming climate. Great Basin terminal lake water extents have decreased by as much as 90 percent over the last 150 years, and terminal lake wetlands have decreased in area by as much as 47 percent since 1984. Lake elevations and wetland areas are primarily supported by freshwater inputs from snowmelt feeding upgradient rivers, streams, and springs. These freshwater inputs have been severely reduced because of continued and increased surface-water diversions and surface-water capture through groundwater pumping for agriculture, mining, and public supply as well as unprecedented drought conditions and warming temperatures related to climate change. Water quality, specifically salinity, is highly variable for terminal lakes of the Great Basin, and this variability is a result of the balance between freshwater inflow and evaporation. Variability of salinity at each of the terminal lakes can be affected by lake morphology, hydrogeologic features of the basin, annual variability in weather patterns, and changes in upgradient water use. Hypersaline terminal lakes provide abundant food resources such as brine shrimp and brine flies that support nesting and migrating birds. The density and composition of invertebrates are closely tied to lake salinity. Increased salinity can exceed the tolerance of invertebrates, severely limiting their biomass. In contrast, decreased salinity can lead to altered invertebrate community composition, reducing the abundance of optimal avian prey resources. Great Basin terminal lake ecosystems, including open-water and adjacent aquatic and terrestrial environments, provide resources necessary to sustain many animal populations throughout the year. Although a variety of taxa use terminal lakes, these ecosystems are of acute importance for the millions of migratory waterbirds (for example, shorebirds, wading birds, and waterfowl) dependent on the network of terminal lakes and their associated wetlands. Migratory birds transiting the Pacific and Central Flyways use Great Basin terminal lake ecosystems throughout the year to feed, nest, and transit between wintering and breeding ranges. As such, successful conservation of birds and their habitats requires coordinated management of water and habitats across the Great Basin network of terminal lakes and wetlands. The linkages between water availability and ecosystem vulnerability of terminal lakes in the Great Basin are not well understood. The vulnerability of terminal lakes is related to the factors driving change and adaptive capacity of the lake ecosystem. Saline lake ecosystems are vulnerable when changes in water quantity affect ecosystem function. Water quantity affects salinity, which affects food webs and habitat; these linkages can be investigated with water-quality and food web monitoring. Water quantity also affects inundated habitat, which can be quantified through remote sensing. It is necessary to quantify hydroclimatic and water use controls on water availability to terminal lakes to assess the response of the ecosystems. Remotely sensed data can provide a broad-scale and long-term synoptic view of terminal lake hydrologic characteristics, but ground observations are required to interpret changes in water quality and ecological functions. Some terminal lake basins have ongoing monitoring and modeling efforts within the Great Basin (for example, Great Salt Lake, Carson River Basin), yet most monitoring locations are hydrologically upgradient and too far away from lake inflows to provide an accurate assessment of hydrological trends for the lake ecosystems. Other terminal lakes have no long-term hydrological monitoring in their respective watersheds (for example, Lake Abert). Ecological data collection in the Great Basin is also insufficient to understand how many birds exist on the landscape, how birds use the mosaic of terminal-lake habitats as an interconnected system, and how Great Basin terminal lakes are linked to the larger continental system of the Pacific and Central Flyways. Across agencies and organizations, tracking bird movement, abundance, and diversity is inconsistent, with some lakes having once- or twice-a-year bird survey efforts and a few locations having more intensive ecological data-gathering efforts (for example, Great Salt Lake, Lake Abert). Bridging hydrological and ecological information gaps will improve understanding of the trends in water supply and water quality, habitat availability and usage, and impacts on vulnerable waterbird species, all of which would be used by managers in coordinated conservation of this unique network of terminal-lake habitats. The terminal lakes of the Great Basin are part of the Basin and Range physiographic province that extends from the Colorado Plateau on the east to the Sierra Nevada on the west, and from the Snake River Plain on the north to the Garlock fault and the Mojave block on the south. The Great Basin is larger than 650,000 square kilometers and encompasses most of the State of Nevada but also extends to western Utah, eastern California, southeastern Idaho, southwestern Wyoming, and southeastern Oregon. The climate is arid to semiarid with a hydrologic regime that is snowmelt dominated, providing as much as 75 percent of total annual runoff for the region. Terminal lakes of the Great Basin occupy the lowest areas of closed (endorheic) drainage basins, such that lake levels and water quality respond rapidly to surface-water inflow. Terminal lakes provide local and regional economic value to the States in the Great Basin, including mineral extraction, aquaculture, public works, and recreational uses. As an example, assessments of Great Salt Lake’s ecological health and economic impact find hemispheric importance for the former and regional importance for the latter. Great Salt Lake creates about 7,000 jobs and $2 billion of economic output per year, most of which would be lost with further declines in lake level. The objectives of this Science Strategy are threefold: (1) to identify how changing water availability affects the quality, diversity, and abundance of habitats supporting continental waterbird populations; (2) to highlight the scientific monitoring and assessment needs of Great Basin terminal lakes; and (3) to support coordinated management and conservation actions to benefit those ecosystems, migratory birds, and other wildlife. There are long-term hydrological, ecological, and societal challenges associated with terminal lakes ecosystems in the Great Basin. This Science Strategy benefits partners by providing a conceptual model, nested at different spatial extents, that identifies key scientific information needs to inform coordinated implementation of management and conservation plans within and among hydrologic basins to address these complex challenges.

California, Idaho, Nevada, Oregon, Utah, Wyoming↗

Geologic framework and hydrogeology of the middle Carson River Basin, Eagle, Dayton, and Churchill Valleys, West-Central Nevada

Changes in land use and water use and increasing development of water resources in the middle Carson River basin may affect flow of the river and, in turn, affect downstream water users dependent on sustained river flows to Lahontan Reservoir. The U.S. Geological Survey, in cooperation with the Bureau of Reclamation, began a study in 2008 of the middle Carson River basin, extending from Eagle Valley to Churchill Valley. Various types of geologic and hydrologic data were compiled from previous studies, collected for this study, and compiled and analyzed to provide a framework for development of a numerical model of the groundwater and surface-water flow systems of the basin. Geologic units that are assumed to have similar hydrologic characteristics were grouped into hydrogeologic units comprised of consolidated rocks of pre-Cenozoic age that underlie a unit of consolidated volcanic rock and semi-consolidated sediments of Tertiary age. The principal aquifer in the study area is comprised of unconsolidated sediments of Quaternary age. The Quaternary sediments include alluvial fan, fluvial, and lake sediments, and were grouped into a basin-fill hydrogeologic unit that overlies the pre-Cenozoic and Tertiary hydrologic units. The thickness of the combined section of Tertiary volcanic and sedimentary rocks and Quaternary basin-fill deposits previously was estimated to range from zero where pre-Cenozoic rocks are exposed to greater than 10,000 feet in the Bull Canyon subbasin, and greater than 6,000 feet on the western side of Churchill Butte and beneath the Desert Mountains. The thickness of Quaternary basin-fill sediments was estimated using gravity data and lithologic descriptions from driller’s logs. The most permeable parts of basin-fill sediments are greater than 1,000 feet thick in the Carson Plains subbasin, greater than 800 feet and 600 feet thick in the western and northeastern parts of the Stagecoach subbasin, and greater than 1,000 feet and 800 feet thick in the northern and southern parts of Churchill Valley, respectively. The distribution of aquifer properties was estimated for basin-fill sediments using slug-test and aquifer test data, and the lithologic descriptions of previously mapped geologic units. Slug-test data show hydraulic conductivity is greater than 10 to greater than 100 feet per day for fluvial sediments near the flood plain, less than 10 feet per day for basin-fill sediments outside the flood plain, and less than 1 foot per day for consolidated rocks. Estimates of transmissivity exceed 20,000 feet squared per day near the Carson River in Dayton, Churchill, and western Lahontan Valleys and in the northern part of the Stagecoach subbasin, and exceed 10,000 feet squared per day in the western part of Churchill Valley. A transmissivity of 90,000 feet squared per day was estimated from results of an aquifer test in the Carson Plains subbasin, indicating that permeable gravel and cobble zones at depths greater than 400 feet supplied water to the pumping well. Estimates of specific yield ranged from less than 1 to 2 percent for most consolidated rocks, from 1 to 15 percent for semi-consolidated Tertiary sediments, and from 10 to 40 percent for unconsolidated basin-fill sediments. Water-level altitude maps based on measurements at about 300 wells in 2009 show water levels have declined as much as 70 feet since 1964 on the northwestern side of Eagle Valley, about 10 feet since 1995 near Dayton in the Carson Plains subbasin, and from 5 to 10 feet since 1982 in the western and northeastern parts of the Stagecoach subbasin and the northwestern part of Churchill Valley. The declines are likely the result of municipal and agricultural pumping. The maps show a groundwater divide between the Carson Plains and Stagecoach subbasins, and a continuous hydraulic gradient between the Stagecoach subbasin and Churchill Valley. Groundwater flow directions are uncertain beneath parts of the boundary of Churchill Valley. The altitude of the top of pre-Cenozoic rocks shows thick sections of saturated Tertiary rocks and sediments south of the Dead Camel Mountains and beneath the eastern part of the Desert Mountains through which groundwater flow between Churchill Valley, Mason Valley, and Lahontan Valley may take place. North of Lahontan reservoir, beneath the Dead Camel Mountains, and beneath the southern part of Adrian Valley, the altitude of pre-Cenozoic rocks indicates groundwater flow between the three valleys is minimal. Streamflow measurements, supported by data on the deuterium content and specific conductance of surface-water samples, indicate a loss of Carson River streamflow in the Riverview subbasin, streamflow gains in the Moundhouse subbasin and the eastern part of the Carson Plains subbasin, and streamflow losses in the Bull Canyon subbasin. Comparisons of fluctuations in groundwater levels to those in stream stage in the Carson Plains subbasin indicate that streamflow lost to infiltration from the Carson River, from irrigation ditches, and from irrigated fields is an important source of groundwater recharge. Fluctuations in groundwater levels compared with the stage of Lahontan Reservoir in Churchill Valley indicate losses to infiltration from the reservoir during high stage and groundwater seepage to the reservoir during low stage.

Nevada↗

Relations between total phosphorus and orthophosphorus concentrations and rainfall, surface-water discharge, and groundwater levels in Big Cypress Seminole Indian Reservation, Florida, 2014–16

The Seminole Tribe of Florida (the Tribe) is partnering with the U.S. Environmental Protection Agency to develop a numeric phosphorus criterion for the 52,000-acre Big Cypress Seminole Indian Reservation (BCSIR), which is located downgradient of the Everglades Agricultural Area, and of other public and private lands, in southeastern Hendry County and northwestern Broward County in southern Florida. The U.S. Geological Survey (USGS), in cooperation with the Tribe, used water-quality data collected between October 2014 and September 2016 by the Tribe and the South Florida Water Management District (SFWMD), along with data from rainfall gages, surface-water stage and discharge gages, and groundwater monitoring wells, to (1) examine the relations between local hydrology and measured total phosphorus (TP) and orthophosphorus (OP) concentrations and (2) identify explanatory variables for TP concentrations. Of particular concern were conditions when TP exceeded 10 parts per billion (ppb) (0.01 milligram per liter [mg/L]) given that the State of Florida and the Miccosukee Tribe of Indians Alligator Alley Reservation (located downstream of the BCSIR) have adopted a 10-ppb maximum TP criterion for surface waters. From October 2014 to September 2016, the Tribe collected 47–52 samples at each of nine water-quality sites for analysis of TP and OP, except at one site where 28 samples were collected. For all sites sampled, concentrations of TP (as phosphorus [P]) ranged from less than 0.002 mg/L (2 ppb) to a maximum of nearly 0.50 mg/L (500 ppb), whereas concentrations of OP (as P), the reactive form of inorganic phosphorus readily absorbed by plants and (or) abiotically absorbed, ranged from less than 0.003 mg/L (3 ppb) to a maximum of 0.24 mg/L (240 ppb). The median and interquartile ranges of concentrations of TP and OP in the samples collected in 2014–16 by the Tribe were similar to the median and interquartile ranges of concentrations in samples collected by the SFWMD at nearby sites during the same period. Differences in concentrations can likely be explained by differences in sample collection methods, sampling locations, sample collection time, and the hydrology during sampling or by the number of samples collected. A major limitation of this study was the short duration of sample collection, which covers a limited range of hydrologic conditions within the BCSIR. The effect of surface-water and groundwater hydrologic conditions on TP and OP concentrations was assessed by using rainfall data and surface-water stage and discharge records. The highest TP and OP concentrations occurred during peak surface-water flows in the canals following long dry periods. Concentrations of TP and OP increased internal to the BCSIR in the western half of the BCSIR during wet periods, but increased concentrations tended to lag behind rainfall events, likely because control structures upstream of sampling sites do not release flows until the water levels in the canals reach predetermined levels. This pattern may indicate that bed sediments in the canals contain high concentrations of phosphorus that becomes resuspended during high flows or that phosphorus salts that had accumulated on dry land during dry periods are carried into the canals by runoff. The largest TP spikes usually occurred at the beginning of high-flow events, but then quickly tapered off even when flows remained high. Groundwater flows were assessed in the BCSIR by using groundwater level observations from two preexisting USGS monitoring well clusters, each characterized by a shallow well installed in the surficial aquifer system and a deeper well installed in the intermediate aquifer system. Groundwater levels were evaluated with respect to surface-water levels and discharge in the BCSIR during the period of surface-water sampling. During dry conditions water levels in canals were often higher than groundwater levels in the surficial aquifer, indicating the potential for surface water to recharge the surficial aquifer. During wetter conditions, this trend reversed, and there was potential for shallow groundwater discharge into the canals. From October 2014 to September 2016, concentrations of TP tended to decrease as surface-water inflows moved across the BCSIR from north to south. In both the western and eastern halves of the reservation, the mean concentration of TP was lower in the surface-water outflows from the BCSIR than in the inflows. The mean concentration of TP in the inflows to the western reservation was 0.04 mg/L (40 ppb), whereas the mean concentration of TP in the outflows was 0.03 mg/L (30 ppb). In the eastern reservation, the mean concentration of TP in the inflows was 0.07 mg/L (70 ppb), whereas the mean concentration of TP in the outflows was 0.04 mg/L (40 ppb). TP and OP concentrations were evaluated relative to other water-quality parameters, including turbidity, suspended solids, nitrate plus nitrite, dissolved oxygen, pH, and specific conductance, to determine if any relations existed between TP and other variables. Weak relations were indicated for turbidity and suspended solids at two sites, which indicates that there may be a relation of increased TP to mobilization of sediment.

Florida↗

Sediment load and distribution in the lower Skagit River, Skagit County, Washington

The Skagit River delivers about 40 percent of all fluvial sediment that enters Puget Sound, influencing flood hazards in the Skagit lowlands, critically important estuarine habitat in the delta, and some of the most diverse and productive agriculture in western Washington. A total of 175 measurements of suspended-sediment load, made routinely from 1974 to 1993, and sporadically from 2006 to 2009, were used to develop and evaluate regression models of sediment transport (also known as “sediment-rating curves”) for estimating suspended-sediment load as a function of river discharge. Using a flow-range model and 75 years of daily discharge record (acquired from 1941 to 2015), the mean annual suspended-sediment load for the Skagit River near Mount Vernon, Washington, was estimated to be 2.5 teragrams (Tg, where 1 Tg = 1 million metric tons). The seasonal model indicates that 74 percent of the total annual suspended‑sediment load is delivered to Puget Sound during the winter storm season (from October through March), but also indicates that discharge is a poor surrogate for suspended‑sediment concentration (SSC) during the summer low-flow season. Sediment-rating curves developed for different time periods revealed that the regression model slope of the SSC-discharge relation increased 66 percent between the periods of 1974–76 and 2006–09 when suspended-sediment samples were collected, implying that changes in sediment supply, channel hydraulics, and (or) basin hydrology occurred between the two time intervals. In the relatively wet water year 2007 (October 1, 2006, through September 30, 2007), an automated sampler was used to collect daily samples of suspended sediment from which an annual load of 4.5 Tg was calculated, dominated by a single large flood event that contributed 1.8 Tg, or 40 percent of the total. In comparison, the annual load calculated for water year 2007 using the preferred flow-range model was 4.8 Tg (+6.7 percent), in close agreement with the measured value. Particle size affects sediment transport, fate and distribution across watersheds, and therefore is important for predicting how coastal environments, particularly deltas and beaches, will respond to changes in climate and sea-level. Particle-size analysis of winter storm samples indicated that about one-half of the suspended-sediment load consisted of fines (that is, silt- and clay-sized particles smaller than 0.0625 mm in diameter), and the remainder consisted of mostly fine- to medium-sized sand (0.0625–0.5 mm), whereas bedload during winter storm flows (about 1–3 percent of total sediment load) was predominantly composed of medium to coarse sand (0.25–1 mm). A continuous turbidity record from the Anacortes Water Treatment Plant (water years 1999–2013), used as a surrogate for the concentration of fines ( R 2 = 0.93, p = 4.2E-10, n = 17), confirms that about one-half of the mean annual suspended-sediment load is composed of fines. The distribution of flow through the delta distributaries (that is, the channels into which the main stem splits as it approaches the delta) is dynamic, with twice as much flow through the North Fork of the Skagit River relative to the South Fork during low-flow conditions, and close to equal flows in the two channels during high-flow conditions. Turbidity, monitored at several locations in the lower river in spring 2009, was essentially uniform among sites, indicating that fines are well mixed in the lower Skagit River system (defined as the Skagit River and all its distributaries downstream of the Mount Vernon streamgage). A strong relation ( R 2 = 0.95, p = 3.2E-14, n = 21; linear regression) between the concentration of fines and turbidity measured at various locations in summer 2009 indicates that turbidity is an effective surrogate for the concentration of fines, independent of location in the river, under naturally well-mixed fluvial conditions. This relation is especially useful for monitoring suspended sediment in western Washington rivers that are seasonally dominated by glacier meltwater because glacial melting typically produces suspended-sediment concentrations that are not well correlated with discharge. These results provide a comprehensive set of tools to estimate sediment delivery and delta responses of interest to scientists and resource managers including decision-makers examining options for flood hazard mitigation, estuary restoration, and climate change adaptation.

Washington↗

Ground water in the Eola-Amity Hills area, northern Willamette Valley, Oregon

The Eola-Amity Hills area ,comprises about 230 square miles on the west side of the Willamette Valley between Salem and McMinnville, Oreg. The area is largely rural, and agriculture is the principal occupation. Rocks ranging in age from Eocene to Recent underlie the area. The oldest rocks are a sequence more than 5,000 feet thick of marine-deposited shale and siltstone strata, with thin interbeds of sandstone that range in age from Eocene to middle Oligocene. They are widely exposed in and west of the Eola-Amity Hills and underlie younger sedimentary and volcanic rocks throughout the study area. In the Eola-Amity Hills and Red Hills of Dundee, the Columbia River Group, a series of eastward-dipping basaltic lava flows locally of Miocene age, and conformably overlies the marine sedimentary rocks. The Columbia River Group ranges in thickness from less than 1 foot to about 900 feet and has an average thickness of about 200 feet. The formation is exposed in the Eola-Amity Hills and Red Hills of Dundee and, at places, extends to the east beneath younger rocks. Overlying the Columbia River Group and marine sedimentary rocks are nonmarine sedimentary deposits that range in thickness from less than 1 foot, where they lap up (to an altitude of about 200 ft) on the flanks of the higher hills, to several hundred feet along the east margin of the study area. These deposits include the Troutdale Formation of Pliocene age, the Willamette Silt of late Pleistocene age, and alluvium of the Willamette River and its tributaries. The Troutdale Formation and the alluvium of the Willamette River contain the most productive aquifers in the Eola-Amity Hills area. These aquifers, which consist mainly of sand and gravel, generally yield moderate to large quantities of water to properly constructed wells. Basalt of the Columbia River Group yields small to moderate quantities of water to wells, and the marine sedimentary rocks and Willamette Silt generally yield small but adequate quantities of water for domestic and stock supplies. Ground water from the Columbia River Group and nonmarine sedimentary rocks is chemically suitable for irrigation and other uses, as is the water from shallow depths in the marine sedimentary rocks. However, water from depths of more than several hundred feet in the marine sedimentary rocks contains large amounts of chloride and other dissolved mineral constitutents that make it unsuitable for most uses. Samples from three fairly closely spaced wells obtaining water from depth zones of 50 to 77, 191 to 201, and about 2,000 feet contained 172, 1,160, and 26,000 ppm (parts per million) of chloride, respectively. About 6,100 acre-feet of ground water was pumped from wells and withdrawn from springs for various uses during 1964; of this amount about 4,800 acre-feet was used for irrigation. The total volume of ground water withdrawn and put to beneficial use each year is small compared with the amount that discharges naturally by evapotranspiration and through undeveloped seeps and springs. Much of the natural discharge could be intercepted and put to beneficial use by pumping from wells. Major problems affecting the development of ground water in the area include a) uneven areal distribution of permeable rocks, (b) undesirable chemical quality of the ground water locally in the marine sedimentary rocks, and (c) fine sand entering wells that tap the Troutdale Formation and thereby causing loss of well efficiency and costly wear on pumps and water-supply systems.

Water Supply Paper↗

Comparative analyses of hydrological responses of two adjacent watersheds to climate variability and change using the SWAT model

Water quality problems in the Chesapeake Bay Watershed (CBW) are expected to be exacerbated by climate variability and change. However, climate impacts on agricultural lands and resultant nutrient loads into surface water resources are largely unknown. This study evaluated the impacts of climate variability and change on two adjacent watersheds in the Coastal Plain of the CBW, using the Soil and Water Assessment Tool (SWAT) model. We prepared six climate sensitivity scenarios to assess the individual impacts of variations in CO 2 concentration (590 and 850 ppm), precipitation increase (11 and 21 %), and temperature increase (2.9 and 5.0 °C), based on regional general circulation model (GCM) projections. Further, we considered the ensemble of five GCM projections (2085–2098) under the Representative Concentration Pathway (RCP) 8.5 scenario to evaluate simultaneous changes in CO 2 , precipitation, and temperature. Using SWAT model simulations from 2001 to 2014 as a baseline scenario, predicted hydrologic outputs (water and nitrate budgets) and crop growth were analyzed. Compared to the baseline scenario, a precipitation increase of 21 % and elevated CO 2 concentration of 850 ppm significantly increased streamflow and nitrate loads by 50 and 52 %, respectively, while a temperature increase of 5.0 °C reduced streamflow and nitrate loads by 12 and 13 %, respectively. Crop biomass increased with elevated CO 2 concentrations due to enhanced radiation- and water-use efficiency, while it decreased with precipitation and temperature increases. Over the GCM ensemble mean, annual streamflow and nitrate loads showed an increase of ∼ 70 % relative to the baseline scenario, due to elevated CO 2 concentrations and precipitation increase. Different hydrological responses to climate change were observed from the two watersheds, due to contrasting land use and soil characteristics. The watershed with a larger percent of croplands demonstrated a greater increased rate of 5.2 kg N ha −1 in nitrate yield relative to the watershed with a lower percent of croplands as a result of increased export of nitrate derived from fertilizer. The watershed dominated by poorly drained soils showed increased nitrate removal due do enhanced denitrification compared to the watershed dominated by well-drained soils. Our findings suggest that increased implementation of conservation practices would be necessary for this region to mitigate increased nitrate loads associated with predicted changes in future climate.

Greensboro Watershed, Tuckahoe Creek Watershed↗

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

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

Oregon↗

Sources, instream transport, and trends of nitrogen, phosphorus, and sediment in the lower Tennessee River basin, 1980-96

In 1997, the U.S. Geological Survey (USGS) began an assessment of the lower Tennessee River Basin as part of the National Water-Quality Assessment Program. Existing nutrient and sediment data from 1980 to 1996 were compiled, screened, and interpreted to estimate watershed inputs from nutrient sources, provide a general description of the distribution and transport of nutrients and sediments in surface water, and evaluate trends in nutrient and sediment concentrations in the lower Tennessee (LTEN) River Basin. Nitrogen inputs from major sources varied widely among tributary basins in the LTEN River Basin. Point source wastewater discharges contributed between 0 and 0.61 tons per square mile per year [(tons/mi 2 )/yr]. Of the nonpoint sources of nitrogen for which inputs were estimated (atmospheric deposition, nitrogen fixation, fertilizer application, and livestock waste) livestock waste contributed the largest input in about two-thirds (7 out of 11) of the tributary basins, and fertilizer application contributed the largest input in the remaining 4 basins. Nitrogen input from fertilizer application was the most variable spatially among the nonpoint sources of nitrogen, ranging from 1.5 to 23 (tons/mi 2 )/yr. Atmospheric deposition estimates varied the least from basin to basin, ranging from 1.6 to 2.0 (tons/mi 2 )/yr. Estimates of nitrogen input from livestock waste ranged between 2.0 to 13 (tons/mi 2 )/yr. The percentage of the input from each of these nonpoint sources that entered the surface-water system is not known. Wastewater discharge contributed between 0 and 0.14 (ton/mi 2 )/yr of phosphorus to tributary basins. Livestock waste contributed most of the input in 8 out of the 11 basins, and fertilizer application contributed the most in the remaining 3 basins. Estimates of phosphorus input for fertilizer application ranged from 0.35 to 5.1 (tons/mi 2 )/yr and from 0.62 to 4.3 (tons/mi 2 )/yr from livestock waste. Reservoirs on the main stem of the Tennessee River and on the Duck and Elk Rivers affect nutrient transport because hydrodynamic conditions in the reservoirs promote assimilation by aquatic plants and deposition of particulate matter. Observed decreases in total nitrite plus nitrate and dissolved-orthophosphorus concentrations in reservoirs or at sites downstream of reservoirs during summer months were probably related to seasonality of plant growth. Nutrient and sediment data used to estimate annual instream loads and yields were compiled from various water-quality monitoring programs and represent the best available data in the LTEN River Basin, but these data have several characteristics that limit accuracy of load estimates. Many of the monitoring programs were not designed with the objective of annual load estimation, and data representing storm transport are, therefore, sparse; sampling and analytical methods varied through time and among the monitoring programs, hampering spatial and temporal comparisons. The load estimates computed from these data are useful for evaluating broad spatial patterns of instream load, and comparisons of instream load to inputs, but may not be sufficiently accurate for local-scale evaluations of water quality. Estimates of the mean annual instream load of total nitrogen entering (Chattanooga, Tenn.) and leaving (Paducah, Ky.) the LTEN River Basin were 29,000 and 60,000 tons per year (tons/yr), respectively. These estimates represent a gain of 31,000 tons/yr, on average, across the area (18,930 mi 2 ) between these inlet and outlet sites. The sum of the mean annual instream load from gaged tributaries to the main stem within the study unit was 14,000 tons/yr; however, this number cannot be directly compared with the gain between the inlet and outlet sites because (1) the gaged area represents only 30 percent of the total area and (2) the period of record at many tributary sites did not correspond with the period of record at the inlet or outlet sites. Estimates of mean annual instream load of total phosphorus at the inlet and outlet sites of the LTEN River Basin were 1,300 and 5,000 tons/yr, respectively, representing a gain of 3,700 tons/yr, on average, across the study unit. The sum of the gaged tributary load, representing only 28 percent of the area contributing to the main stem, was 4,300 tons/yr. Although this number cannot be closely compared with the gain throughout the study unit, for the same reasons given for total nitrogen, a general comparison suggests that the main stem of the Tennessee River and the tributary embayments along the main stem function as a sink for total phosphorus, removing a substantial amount from the water column through deposition or assimilation. The estimates of inputs can be compared and correlated with yields (area-normalized instream loads); significant correlations between estimates of inputs and yields might be useful as predictive tools for instream water quality where monitoring data are not available. Yields of nitrogen correlated moderately well with inputs from nonpoint sources, based on 1992 estimates. Nitrogen yield was highest [3.5 (tons/mi 2 )/yr] for Town Creek, for which the balance of nonpoint-source inputs to agricultural lands (fertilizer application plus nitrogen fixation plus livestock waste minus harvest) was also the highest [15 (tons/mi 2 )/yr]. Nitrogen yield was low [1.0 (tons/mi 2 )/yr] for the Buffalo River, for which the balance of agricultural nonpoint-source input was correspondingly low [3.2 (tons/mi 2 )/yr, the second lowest]. Correlation of wastewater discharge with yield was poor, and contrasted with the significant correlation between wastewater discharge and median nitrogen concentration during low streamflow. The poor correlation between wastewater discharge and annual yield was expected, however, as wastewater discharge is a small fraction compared with annual yield. In contrast with nitrogen, phosphorus yield did not correlate well with any estimated inputs or land-use types for the tributary basins. Phosphorus yield was highest [1.1 and 0.93 (tons/mi 2 )/yr] at two sites along the Duck River and at Elk River near Prospect [0.89 (ton/mi 2 )/yr]; however, estimates of inputs at these sites were in the middle of their respective ranges. The influence of the outcrop of phosphatic limestone formations of the brown-phosphate districts in the lower Duck and lower Elk River Basins might be responsible for the poor correlation between estimated inputs and yields of phosphorus. The outcrop pattern of these phosphatic limestones are an important factor to consider as regional boundaries are established for attainable, region-specific water-quality criteria for total phosphorus. Estimates of sediment input from cropland soil erosion in 1992 ranged from 51 to 540 (tons/mi 2 )/yr among the major hydrologic units in the LTEN River Basin. Information was not available to estimate this input for individual tributaries. Sediment yield estimates ranged from 65 to 263 (tons/mi 2 )/yr for the three tributary monitoring basins for which instream data were available, and from 17 to 26 (tons/mi 2 )/yr for the Tennessee River at South Pittsburg and at Pickwick Landing Dam, respectively. Lower sediment yields for the main stem sites compared with the tributary sites is probably due to sediment deposition in the main stem of the Tennessee River and tributary embayments along the main stem. Most of the significant trends in nutrient concentrations from about 1985 to about 1995 were decreasing trends, except for total nitrite plus nitrate, which increased at one site on the Elk River. The spatial distribution of decreasing trends of total nitrogen and total ammonia corresponds with the spatial variation among basins in wastewater loading rate. The time period of observed trends corresponds to the period of improvements in municipal treatment, thus decreases in wastewater effluent concentrations of nitrogen might be responsible for the decreasing trend in instream concentrations at these sites. Concentrations of total phosphorus did not decrease during this period at these sites, as might have been expected considering the reductions in wastewater input of phosphorus during this period.

Alabama, Georgia, Kentucky, Mississippi, Tennessee↗

Simulation of groundwater flow, effects of artificial recharge, and storage volume changes in the Equus Beds aquifer near the city of Wichita, Kansas well field, 1935–2008

The Equus Beds aquifer is a primary water-supply source for Wichita, Kansas and the surrounding area because of shallow depth to water, large saturated thickness, and generally good water quality. Substantial water-level declines in the Equus Beds aquifer have resulted from pumping groundwater for agricultural and municipal needs, as well as periodic drought conditions. In March 2006, the city of Wichita began construction of the Equus Beds Aquifer Storage and Recovery project to store and later recover groundwater, and to form a hydraulic barrier to the known chloride-brine plume near Burrton, Kansas. In October 2009, the U.S. Geological Survey, in cooperation with the city of Wichita, began a study to determine groundwater flow in the area of the Wichita well field, and chloride transport from the Arkansas River and Burrton oilfield to the Wichita well field. Groundwater flow was simulated for the Equus Beds aquifer using the three-dimensional finite-difference groundwater-flow model MODFLOW-2000. The model simulates steady-state and transient conditions. The groundwater-flow model was calibrated by adjusting model input data and model geometry until model results matched field observations within an acceptable level of accuracy. The root mean square (RMS) error for water-level observations for the steady-state calibration simulation is 9.82 feet. The ratio of the RMS error to the total head loss in the model area is 0.049 and the mean error for water-level observations is 3.86 feet. The difference between flow into the model and flow out of the model across all model boundaries is -0.08 percent of total flow for the steady-state calibration. The RMS error for water-level observations for the transient calibration simulation is 2.48 feet, the ratio of the RMS error to the total head loss in the model area is 0.0124, and the mean error for water-level observations is 0.03 feet. The RMS error calculated for observed and simulated base flow gains or losses for the Arkansas River for the transient simulation is 7,916,564 cubic feet per day (91.6 cubic feet per second) and the RMS error divided by (/) the total range in streamflow (7,916,564/37,461,669 cubic feet per day) is 22 percent. The RMS error calculated for observed and simulated streamflow gains or losses for the Little Arkansas River for the transient simulation is 5,610,089 cubic feet per day(64.9 cubic feet per second) and the RMS error divided by the total range in streamflow (5,612,918/41,791,091 cubic feet per day) is 13 percent. The mean error between observed and simulated base flow gains or losses was 29,999 cubic feet per day (0.34 cubic feet per second) for the Arkansas River and -1,369,250 cubic feet per day (-15.8 cubic feet per second) for the Little Arkansas River. Cumulative streamflow gain and loss observations are similar to the cumulative simulated equivalents. Average percent mass balance difference for individual stress periods ranged from -0.46 to 0.51 percent. The cumulative mass balance for the transient calibration was 0.01 percent. Composite scaled sensitivities indicate the simulations are most sensitive to parameters with a large areal distribution. For the steady-state calibration, these parameters include recharge, hydraulic conductivity, and vertical conductance. For the transient simulation, these parameters include evapotranspiration, recharge, and hydraulic conductivity. The ability of the calibrated model to account for the additional groundwater recharged to the Equus Beds aquifer as part of the Aquifer Storage and Recovery project was assessed by using the U.S. Geological Survey subregional water budget program ZONEBUDGET and comparing those results to metered recharge for 2007 and 2008 and previous estimates of artificial recharge. The change in storage between simulations is the volume of water that estimates the recharge credit for the aquifer storage and recovery system. The estimated increase in storage of 1,607 acre-ft in the basin storage area compared to metered recharge of 1,796 acre-ft indicates some loss of metered recharge. Increased storage outside of the basin storage area of 183 acre-ft accounts for all but 6 acre-ft or 0.33 percent of the total. Previously estimated recharge credits for 2007 and 2008 are 1,018 and 600 acre-ft, respectively, and a total estimated recharge credit of 1,618 acre-ft. Storage changes calculated for this study are 4.42 percent less for 2007 and 5.67 percent more for 2008 than previous estimates. Total storage change for 2007 and 2008 is 0.68 percent less than previous estimates. The small difference between the increase in storage from artificial recharge estimated with the groundwater-flow model and metered recharge indicates the groundwater model correctly accounts for the additional water recharged to the Equus Beds aquifer as part of the Aquifer Storage and Recovery project. Small percent differences between inflows and outflows for all stress periods and all index cells in the basin storage area, improved calibration compared to the previous model, and a reasonable match between simulated and measured long-term base flow indicates the groundwater model accurately simulates groundwater flow in the study area. The change in groundwater level through recent years compared to the August 1940 groundwater level map has been documented and used to assess the change of storage volume of the Equus Beds aquifer in and near the Wichita well field for three different areas. Two methods were used to estimate changes in storage from simulation results using simulated change in groundwater levels in layer 1 between stress periods, and using ZONEBUDGET to calculate the change in storage in the same way the effects of artificial recharge were estimated within the basin storage area. The three methods indicate similar trends although the magnitude of storage changes differ. Information about the change in storage in response to hydrologic stresses is important for managing groundwater resources in the study area. The comparison between the three methods indicates similar storage change trends are estimated and each could be used to determine relative increases or decreases in storage. Use of groundwater level changes that do not include storage changes that occur in confined or semi-confined parts of the aquifer will slightly underestimate storage changes; however, use of specific yield and groundwater level changes to estimate storage change in confined or semi-confined parts of the aquifer will overestimate storage changes. Using only changes in shallow groundwater levels would provide more accurate storage change estimates for the measured groundwater levels method. The value used for specific yield is also an important consideration when estimating storage. For the Equus Beds aquifer the reported specific yield ranges between 0.08 and 0.35 and the storage coefficient (for confined conditions) ranges between 0.0004 and 0.16. Considering the importance of the value of specific yield and storage coefficient to estimates of storage change over time, and the wide range and substantial overlap for the reported values for specific yield and storage coefficient in the study area, further information on the distribution of specific yield and storage coefficient within the Equus Beds aquifer in the study area would greatly enhance the accuracy of estimated storage changes using both simulated groundwater level, simulated groundwater budget, or measured groundwater level methods.

Kansas↗

Identification of metrics to monitor salt marsh integrity on National Wildlife Refuges in relation to conservation and management objectives

Executive Summary Most salt marshes in the US have been degraded by human activities, and threats from physical alterations, surrounding land-use, species invasions, and global climate change persist. Salt marshes are unique and highly productive ecosystems with high intrinsic value to wildlife, and many National Wildlife Refuges (NWRs) have been established in coastal areas to protect large tracts of salt marsh and wetland-dependent species. Various management practices are employed routinely on coastal NWRs to restore and enhance marsh integrity and ensure ecosystem sustainability. Prioritizing NWR salt marshes for application of management actions and choosing among multiple management options requires scientifically-based methods for assessing marsh condition. Monitoring is integral to structured decision-making (SDM), a formal process for decomposing a decision into its essential elements. Within a natural resource context, SDM involves identifying management objectives, alternative management actions, and expected management outcomes. The core of SDM is a set of criteria for measuring system performance and evaluating management responses. Therefore, use of SDM to frame natural resource decisions leads to logical selection of monitoring attributes that are linked explicitly to management needs. We used SDM to guide selection of variables for monitoring the ecological integrity of salt marshes within the National Wildlife Refuge System (NWRS). Our objectives were to identify indicators of salt marsh integrity that are effective across large geographic regions, responsive to a wide range of threats, and feasible to implement within funding and staffing constraints of the NWRS. In April, 2008, we engaged interdisciplinary experts in a week-long rapid prototyping SDM workshop to define the essential elements of salt marsh management decisions on refuges throughout the northeastern, southwestern, and northwestern US, corresponding to respective Regions 5, 2, and 1 of the US Fish and Wildlife Service (FWS). Through this process we identified measurable attributes for monitoring salt marsh ecosystems that are integrated into conservation practice and target management objectives. The following salt marsh attributes were identified through the SDM process either for describing state condition to determine management needs or for evaluating the achievement of management objectives: historical condition and geomorphic setting; ditch density; surrounding land use; ratio of open water area to vegetation area; rate of pesticide application; environmental contaminant concentration; change in marsh surface elevation relative to sea level rise; tidal range and groundwater level; surface topography; salinity; and species composition and abundance of vegetation, invasive species, invertebrates, nekton, and breeding and wintering birds. The identified attributes were too broadly defined to serve as operational monitoring variables. Therefore, we tested specific metrics for quantifying most of these attributes in summers of 2008 and 2009. The first four attributes in the above list can be characterized by office-based analysis of existing GIS data layers. The remaining attributes require field-based methods for assessment. We were forced to exclude a small number of attributes from field tests due to inconsistent data (pesticide application rate, environmental contaminant concentrations) or requirements that exceeded the scope of this project (change in marsh surface elevation; surface topography; benthic invertebrates; wintering birds). We evaluated potential metrics for evaluating all remaining field attributes. In partnership with NWRS biologists, we tested rapid versus intensive metrics for monitoring field attributes (tidal range and groundwater level; marsh surface elevation; salinity; and species composition and abundance of vegetation, invasive species, nekton, and breeding birds) at coastal refuges throughout FWS Region 5. Seven refuges participated in metric testing in 2008: Rachel Carson (ME), Parker River (MA), Wertheim (NY), E. B. Forsythe (NJ), Bombay Hook (DE), Prime Hook (DE), and Eastern Shore of Virginia Complex (VA). These seven and two additional refuges participated in metric testing in 2009: Rhode Island Complex (RI) and Stewart B. McKinney (CT). We based all field metrics on existing protocols for salt marsh assessment. Sampling locations were determined randomly within delineated marsh study units (MSUs) at each refuge. Detailed field methods are provided in appendices to this report. Measurements for individual metrics were averaged across samples within MSUs during each year of sampling. Each year, correlation or regression analysis was conducted on average measurements across MSUs within each attribute set to identify redundant metrics. Statistical redundancy between a pair of metrics within an attribute set (i.e., correlation or regression slopes with p-values < 0.05) was considered justification for eliminating one of the pair from the regional set of monitoring metrics. Decisions regarding metric elimination versus retention were based on feasibility of monitoring, considering such factors as sampling time, resources required, and potential for regional standardization in implementation. The result of these tests is a reduced suite of monitoring metrics that targets NWRS management decisions and is practicable for implementing on a regional scale. Based on these tests, we recommend the following list of metrics for monitoring integrity of NWRS salt marshes (marsh attribute category is in parentheses): (historical condition and geomorphic setting) position of marsh in the landscape, marsh shape, degree of fill and/or fragmentation, degree of tidal flushing, amount of aquatic edge; (ditch density) ranking of ditch density from none to severe; (surrounding land use) relative proportion of agricultural land in a 150-m buffer around the marsh, relative proportion of natural land in a 150-m buffer around the marsh, relative proportion of natural land in a 1-km buffer around the marsh; (ratio of open water area to vegetation area) ratio of open water to emergent herbaceous wetlands within the marsh; (marsh surface elevation) elevation referenced to NAVD88 in a representative area of the marsh; (tidal range and groundwater level) percent of time the marsh surface is flooded during deployment of a continuous water-level monitor at a representative marsh location, mean depth of surface flooding as measured by a continuous water-level monitor at a representative location; (salinity) salinity measured in surface water; (vegetation community) vegetation species richness using the point-intercept method in 100-m diameter survey plots, percent cover of various marsh community types within 100-m diameter survey plots; (invasive species abundance) percent cover of invasive plant species measured using the point-intercept method in 100-m diameter survey plots; (nekton community) nekton density, nekton species richness, length of Fundulus heteroclitus ; (breeding bird community) abundance of Willets counted per point during standard call-broadcast surveys, summed abundance of tidal marsh obligate species (Clapper Rail, Willet, Saltmarsh Sparrow, Seaside Sparrow) counted per point during standard call-broadcast surveys. Metrics describing the historical condition, geomorphic setting, and broad landscape features can be assessed using existing GIS databases. Our results support use of rapid methods to assess the majority of field metrics; only those used to describe the nekton community must be measured using intensive methods (throw traps or ditch nets, dependant on habitat configuration). Implementation of these metrics for quantitative assessment of NWRS salt marsh integrity in FWS Region 5 requires developing sampling designs for each refuge. Additionally, it is important to determine how the monitoring information will be used within a management context. SDM should be used to complete the analysis of salt marsh management decisions. The next steps would involve 1) prioritizing and weighting the management objectives; 2) predicting responses to individual management actions in terms of objectives and metrics; 3) using multiattribute utility theory to convert all measurable attributes to a common utility scale; 4) determining the total management benefit of each action by summing utilities across objectives; and 5) maximizing the total management benefits within cost constraints for each refuge. This process would allow the optimum management decisions for NWRS salt marshes to be selected and implemented based directly on monitoring data and current understanding of marsh responses to management actions. Monitoring the outcome of management actions would then allow new monitoring data to be incorporated into subsequent decisions.

Connecticut, Delaware, Maine, Massachusetts, New J↗

Sediment transport and deposition in the lower Missouri River during the 2011 flood

Floodwater in the Missouri River in 2011 originated in upper-basin regions and tributaries, and then travelled through a series of large flood-control reservoirs, setting records for total runoff volume entering all six Missouri River main-stem reservoirs. The flooding lasted as long as 3 months. The U.S Geological Survey (USGS) examined sediment transport and deposition in the lower Missouri River in 2011 to investigate how the geography of floodwater sources, in particular the decanting effects of the Missouri River main-stem reservoir system, coupled with the longitudinal characteristics of civil infrastructure and valley-bottom topography, affected sediment transport and deposition in this large, regulated river system. During the flood conditions in 2011, the USGS, in cooperation with the U.S. Army Corps of Engineers, monitored suspended-sediment transport at six primary streamgages along the length of the lower Missouri River. Measured suspended-sediment concentration (SSC) in the lower Missouri River varied from approximately 150 milligrams per liter (mg/L) to 2,000 mg/L from January 1 to September 30, 2011. Median SSC increased in the downstream direction from 355 mg/L at Sioux City, Iowa, to 490 mg/L at Hermann, Missouri. The highest SSCs were measured downstream from Omaha, Nebraska, in late February when snowmelt runoff from tributaries, which were draining zones of high-sediment production, was entering the lower Missouri River, and releases of water at Gavins Point Dam were small. The combination of dilute releases of water at Gavins Point Dam and low streamflows in lower Missouri River tributaries caused sustained lowering of SSC at all streamgages from early July through late August. Suspended-sediment ranged from 5 percent washload (PW; percent silt and clay) to as much as 98 percent in the lower Missouri River from January 1 to September 30, 2011. Median PW increased in the downstream direction from 24 percent at Sioux City, Iowa, to 78 percent at Hermann, Missouri. Measurements made in early January, when SSC was low, indicate that suspended sediment mostly was composed of bed material, but by mid-February, runoff from the plains caused PW to increase at most streamgages. Total suspended-sediment discharge (SSD) during water year 2011 at the selected streamgages in the lower Missouri River ranged from approximately 29 to 64 million tons. Total estimated SSD had the lowest exceedance frequencies in the reaches between Gavins Point Dam and Nebraska City, Nebraska, but exceedance frequencies increased substantially downstream. In 2011, total SSD with low exceedance frequencies were reported at Sioux City, Iowa, Omaha, Nebraska, and Nebraska City, Nebraska, despite moderate-to-high exceedance frequencies for annual average SSC, indicating that the duration of high-magnitude flooding was the primary driver of total SSD. Comparison of median SSC for samples from water year 2011 with samples in the 20 years prior indicated that median SSC for high-action streamflows (streamflows likely to produce a stage exceeding the National Weather Service&rsquo;s &ldquo;action stage&rdquo;) in 2011 were lower than those typical for high-action streamflows. Multiple-comparison analysis indicated that median SSC values for low-action streamflows (streamflows likely to produce stages lower than the National Weather Service&rsquo;s &ldquo;action stage&rdquo;) and high-action streamflows sampled in 2011 at 4 of 6 streamgages were not significantly distinguishable from median SSC values for low-action streamflows in the previous 20 years. Longitudinal comparison of streamflow and SSD exceedance frequencies for 2011 with corresponding frequencies for 2008 and 1993 indicated the important role of tributary contributions to total SSD in the lower Missouri River. In 1993 and 2008, tributaries were the primary source of floodwater in the lower Missouri River, which resulted in a 20-fold increase in total SSD from Sioux City, Iowa, to Hermann, Missouri. In 2011, releases at Gavins Point Dam were the primary source of floodwater in the lower Missouri River, and total SSD at Hermann, Missouri, was only twice that estimated for Sioux City, Iowa. Sand deposition was estimated using analysis of multispectral satellite imagery collected in October and November 2011. Distributions of sand in the flood plain of the lower Missouri River also were quantified in relation to distance from the banks of the main channel for seven discrete river segments bounded by Gavins Point Dam and selected downstream tributaries. The areal extent of overbank flooding and flood-plain sand deposits increased downstream from Sioux City, Iowa to a broad peak near Rulo, Nebraska, and then decreased to levels near the lower limit of quantification downstream from Kansas City, Missouri. Most of the flood plain inundation and sediment-deposition damage to agricultural fields was observed between river miles 480 and 700, where 2011 peak streamflows had low exceedance frequencies, and the lower Missouri River channel was less incised or had aggraded recently. As channel capacity increased in the downstream direction, the relative magnitude of the flood decreased downstream, and overbank flooding was less extensive. In the constricted reaches, flood-plain sand deposits mainly were observed in association with levee breaks.

Iowa, Kansas, Missouri, Montana, Nebraska, North D↗

Drinking water health standards comparison and chemical analysis of groundwater for 72 domestic wells in Bradford County, Pennsylvania, 2016

Pennsylvania has the second highest number of residential wells of any state in the Nation with approximately 2.4 million residents that depend on groundwater for their domestic water supply. Despite the widespread reliance on groundwater in rural areas of the state, publicly available data to characterize the quality of private well water are limited. In Bradford County, more than half of the residents use groundwater from private domestic-supply wells as their primary drinking source. The quality of private well water is influenced by the regional and local setting, including the surrounding soil, geology, land use, household plumbing, and well construction. The groundwater used for domestic water supply in Bradford County is obtained primarily from shallow bedrock and from unconsolidated (glacial) deposits that overlie the bedrock. Historical land use has been predominately forested, agricultural, and residential, but more recently unconventional oil/gas development has been distributed throughout the landscape. Pennsylvania is one of only two states in the Nation without statewide water-well construction standards. To better assess the quality of groundwater used for drinking water supply in Bradford County, data for 72 domestic wells were collected and analyzed for a wide range of constituents that could be evaluated in relation to drinking water health standards, geology, land use, and other environmental factors. Groundwater samples were collected from May through August 2016 and analyzed for physical and chemical properties, including major ions, nutrients, trace elements, volatile organic compounds, ethylene and propylene glycol, alcohols, gross-alpha/beta-particle activity, uranium, radon-222, and dissolved gases. A subset of samples was analyzed for radium isotopes (radium-226 and -228) and for the isotopic composition of methane. This study was conducted by the U.S. Geological Survey in cooperation with the Northern Tier Regional Planning and Development Commission and is part of a regional effort to characterize groundwater in rural areas of Pennsylvania. Results of the 2016 study show that groundwater quality generally met most drinking-water standards. However, a percentage of samples failed to meet maximum contaminant levels (MCLs) for total coliform bacteria (49.3 percent), Escherichia coli (8.5 percent), barium (2.8 percent), and arsenic (2.8 percent); and secondary maximum contaminant levels (SMCL) for sodium (48.6 percent), manganese (30.6 percent), gross alpha and beta activity (16.7 percent), iron (11.1 percent), pH (8.3 percent), total dissolved solids (5.6 percent), chloride (1.4 percent), and aluminum (1.4 percent). Radon-222 activities exceeded the proposed drinking-water standard of 300 picocuries per liter (pCi/L) in 70.4 percent of the samples. There were no exceedances of drinking water health standards for any volatile organic compounds, and the only detections were for three trihalomethanes in one sample. The pH of the groundwater had a large influence on chemical characteristics and ranged from 6.18 to 9.31. Generally, the higher pH samples had higher potential for elevated concentrations of several constituents, including total dissolved solids, sodium, lithium, chloride, fluoride, boron, arsenic, and methane. For the Bradford County well-water samples, calcium/bicarbonate type waters were most abundant, with others classified as sodium/bicarbonate or mixed water types including calcium-sodium/bicarbonate, calcium-sodium/bicarbonate-chloride, sodium/bicarbonate-chloride, sodium/bicarbonate-sulfate, or sodium/chloride types. Six principal components (pH, redox, hardness, chloride-bromide, strontium-barium, and molybdenum-arsenic) explained nearly 78.3 percent of the variance in the groundwater dataset. Groundwater from 12.5 percent of the wells had concentrations of methane greater than the Pennsylvania action level of 7 milligrams per liter (mg/L); detectable methane concentrations ranged from 0.01 to 77 mg/L. In addition, low levels of ethane (as much as 0.13 mg/L) were present in seven samples with the highest methane concentrations. The isotopic composition of methane in five of these groundwater samples was consistent with the isotopic compositions reported for mud-gas logging samples from these geologic units and a thermogenic source. Isotopic composition from a sixth sample suggested the methane in that sample may be of microbial origin. Well-water samples with the higher methane concentrations also had higher pH values and elevated concentrations of sodium, lithium, boron, fluoride, arsenic, and bromide. Relatively elevated concentrations of some other constituents, such as barium and chloride, commonly were present in, but not limited to, those well-water samples with elevated methane. Four of the six groundwater samples with the highest methane concentrations had chloride/bromide ratios that indicate mixing with a small amount of brine (0.02 percent or less) similar in composition to those reported for gas and oil well brines in Pennsylvania. In several other eastern Pennsylvania counties where gas drilling is absent, groundwater with comparable chloride/bromide ratios and chloride concentrations have been reported, implying a potential natural source of brine. Most of Bradford County well-water samples have chloride concentrations less than 20 mg/L, and those with higher chloride concentrations have chloride/bromide ratios that indicate anthropogenic sources (such as road-deicing salt and septic effluent) or brine. Brines that are naturally present may originate from deeper parts of the aquifer system, whereas anthropogenic sources are more likely to affect shallow groundwater because they occur on or near the land surface. The available data for this study indicate that no one physical factor, such as the topographic setting, well depth, or altitude at the bottom of the well, was particularly useful for predicting those well locations with an elevated dissolved concentration of methane. The 2016 assessment of groundwater quality in Bradford County shows groundwater is generally of good quality, but methane and some constituents that occur in high concentration in naturally occurring brine and also in produced waters may be present at low to moderate concentrations in groundwater in various parts of the aquifer.

Pennsylvania↗

Evaluation of agricultural best-management practices in the Conestoga River headwaters, Pennsylvania: Characterization of surface-runoff and ground-water quantity and quality in a small carbonate basin near Churchtown, Pennsylvania, prior to terracing and implementation of nutrient management: Water-quality study of the Conestoga River headwaters, Pennsylvania

The U.S. Geological Survey, in cooperation with the Pennsylvania Department of Environmental Protection1 , conducted a study as part of the U.S. Department of Agriculture's Rural Clean Water Program to determine the effects of agricultural best-management practices on surface-water and ground-water quality in the Conestoga River headwaters basin. This report describes Field-Site 1 and characterizes the surface-runoff and ground-water quantity and quality at the site from January 1983 through September 1984, before the implementation of terracing and nutrient-management best-management practices. The 22.1-acre site, part of two dairy farms, was cropland used primarily for the production of corn and alfalfa, and is underlain by carbonate rock. During the 21-month study period, 91.2 inches of precipitation fell, of which 66 percent occurred during the 1984 water year. Of the 169 storms of 0.10 inch or larger, 97 produced measurable runoff. The average annual application of nutrients to the 14.4 acres of cornfields was 410 pounds per acre of nitrogen and 110 pounds per acre of phosphorus. About three times more nutrients were applied during 1984 water year than during the 1983 water year. The approximate soil-nitrate concentration as nitrogen for the study period was 8.2 milligrams per kilogram in the top 8 inches of soil. Runoff for the study period totalled 714,000 cubic feet, or 9.8 percent of total precipitation. Eighty-eight percent of the runoff occurred during the 1984 water year. Regression analyses indicate that total runoff was controlled primarily by total precipitation amounts during storms and by antecedent moisture conditions; mean event and maximum instantaneous discharges are controlled by precipitation intensity and by antecedent moisture conditions. Regression analyses also suggest that crop cover on the corn acreage affected runoff amounts and rates. Mean concentrations in runoff ranged from 24 to 73,000 milligrams per liter for suspended sediment, 0.45 to 24 milligrams per liter for total phosphorus, and 1.5 to 55 milligrams per liter for total nitrogen. Of the total nitrogen in instantaneous runoff samples, a median of 72 percent was total organic nitrogen. Results of multiple regression analyses suggest that mean nitrogen concentrations in runoff (1) increased with increased surface-nitrogen applications made prior to runoff, and (2) were diluted with increased precipitation duration. Runoff from storms on frozen ground produced the highest loads of nitrogen in runoff. Over the 21-month study period, an estimated 258 tons of suspended sediment, 314 pounds of nitrogen, and 176 pounds of phosphorus were transported with runoff from the site. Of this, 88 percent of the nitrogen and phosphorus loads and 94 percent of the suspended-sediment load were discharged during the wet 1984 water year. The loads for the study period represent 2.5 percent of the total nitrogen and 5.5 percent of the total phosphorus applied to the site as manure and commercial fertilizer. The ground-water basin at the site is slightly larger than the surface-water basin. Ground-water levels, except at one well, responded quickly to recharge; water levels peaked several hours to a day after precipitation. Surface-applied materials moved rapidly to the water table through macropores and slowly by transport through micropores in the unsaturated zone. Median concentrations of dissolved nitrate ranged from 9.2 to 13 milligrams per liter as nitrogen for ground-water samples collected monthly from five wells and a spring. Dissolved nitrate comprised 93 percent (median percentage) of the total nitrogen in ground-water samples. The effect of recharge on ground-water nitrate concentrations was affected by variations in nutrient availability at the land surface and in the unsaturated zone. A lag time of 1 to 3 months was observed between the time that nitrogen was applied to the land surface and local maximums in nitrate concentrations were detected in ground water unaffected by recharge events. Approximately 3,187,000 cubic feet of ground water and an associated 2,200 pounds of nitrate-nitrogen discharged from the site during the study period. For the study period, 42 percent of the precipitation recharged to ground water, 10 percent became runoff, and 48 percent evapotranspired. Inputs of nitrogen to the study area were estimated to be 93 percent from manure, 5 percent from commercial fertilizer, and 2 percent from precipitation. Nitrogen outputs from the system were estimated to be 38 percent to crop uptake, 39 percent to volatilization, 20 percent to ground-water discharge, and 3 percent to surface runoff.

Pennsylvania↗

Status and understanding of groundwater quality in the Monterey Bay and Salinas Valley Basins, 2005-California GAMA Priority Basin Project

Groundwater quality in the approximately 1,000 square mile (2,590 km2) Monterey Bay and Salinas Valley Basins (MS) study unit was investigated as part of the Priority Basin Project of the Groundwater Ambient Monitoring and Assessment (GAMA) Program. The study unit is located in central California in Monterey, Santa Cruz, and San Luis Obispo Counties. The GAMA Priority Basin Project is being conducted by the California State Water Resources Control Board in collaboration with the U.S. Geological Survey (USGS) and the Lawrence Livermore National Laboratory. The GAMA MS study was designed to provide a spatially unbiased assessment of the quality of untreated (raw) groundwater in the primary aquifer systems (hereinafter referred to as primary aquifers). The assessment is based on water-quality and ancillary data collected in 2005 by the USGS from 97 wells and on water-quality data from the California Department of Public Health (CDPH) database. The primary aquifers were defined by the depth intervals of the wells listed in the CDPH database for the MS study unit. The quality of groundwater in the primary aquifers may be different from that in the shallower or deeper water-bearing zones; shallow groundwater may be more vulnerable to surficial contamination. The first component of this study, the status of the current quality of the groundwater resource, was assessed by using data from samples analyzed for volatile organic compounds (VOC), pesticides, and naturally occurring inorganic constituents, such as major ions and trace elements. This status assessment is intended to characterize the quality of groundwater resources in the primary aquifers of the MS study unit, not the treated drinking water delivered to consumers by water purveyors. Relative-concentrations (sample concentration divided by the health- or aesthetic-based benchmark concentration) were used for evaluating groundwater quality for those constituents that have Federal and (or) California regulatory or non-regulatory benchmarks for drinking-water quality. A relative-concentration greater than (>) 1.0 indicates a concentration greater than a benchmark, and less than or equal to (&le;) 1.0 indicates a concentration less than or equal to a benchmark. Relative-concentrations of organic and special interest constituents [perchlorate, N-nitrosodimethylamine (NDMA), and 1,2,3-trichloropropane (1,2,3-TCP)], were classified as "high" (relative-concentration > 1.0), "moderate" (0.1 < relative-concentration &le; 1.0), or "low" (relative-concentration &le; 0.1). Relative-concentrations of inorganic constituents were classified as "high" (relative-concentration > 1.0), "moderate" (0.5 < relative-concentration &le; 1.0), or "low" (relative-concentration &le; 0.5). Aquifer-scale proportion was used as the primary metric in the status assessment for evaluating regional-scale groundwater quality. High aquifer-scale proportion was defined as the percentage of the area of the primary aquifers with a relative-concentration greater than 1.0 for a particular constituent or class of constituents; percentage is based on an areal rather than a volumetric basis. Moderate and low aquifer-scale proportions were defined as the percentage of the primary aquifers with moderate and low relative-concentrations, respectively. Two statistical approaches-grid-based and spatially weighted-were used to evaluate aquifer-scale proportions for individual constituents and classes of constituents. Grid-based and spatially-weighted estimates were comparable in the MS study unit (within 90-percent confidence intervals). Inorganic constituents with human-health benchmarks were detected at high relative-concentrations in 14.5 percent of the primary aquifers, moderate in 35.5 percent, and low in 50.0 percent. High aquifer-scale proportion of inorganic constituents primarily reflected high aquifer-scale proportions of nitrate (7.9 percent), molybdenum (2.9 percent), arsenic (2.8 percent), boron (1.9 percent), and gross alpha-beta radioactivity (1.5 percent). Relative-concentrations of organic constituents (one or more) were high in 0.2 percent, moderate in 6.6 percent, and low in 93.2 percent (not detected in 48.1 percent) of the primary aquifers. The high aquifer-scale proportion of organic constituents primarily reflected high aquifer-scale proportions of tetrachloroethene (0.1 percent) and methyl tert-butyl ether (0.1 percent). Relative-concentration for inorganic constituents with secondary maximum contaminant levels, manganese, total dissolved solids, iron, sulfate, and chloride were high in 18.6, 8.6, 7.1, 2.9, and 1.4 percent of the primary aquifers, respectively. Of the 205 organic and special-interest constituents analyzed, 32 constituents were detected. One organic constituent, the herbicide simazine, was frequently detected (in 10 percent or more of samples), but was detected at low relative-concentrations. The second component of this study, the understanding assessment, identified the natural and human factors that affect groundwater quality by evaluating land use, physical characteristics of the wells, and geochemical conditions of the aquifer. Results from these evaluations were used to explain the occurrence and distribution of constituents in the study unit. The understanding assessment indicated that most wells that contained nitrate were classified as being in agricultural land-use areas, and depths to the top of perforations in most of the wells were less than 350 ft (76 m). High and moderate relative-concentrations of arsenic may be attributed to reductive dissolution of manganese or iron oxides, or to desorption or inhibition of arsenic sorption under alkaline conditions. Arsenic concentrations increased with increasing groundwater depth and residence time (age). Simazine was detected more often in groundwater from wells with surrounding land use classified as agricultural or urban, and with top of perforation depths less than 200 ft (61 m), than in groundwater from wells with natural land use or with deeper depths. Tritium, helium-isotope, and carbon-14 data were used to classify the predominant age of groundwater samples into three categories: modern (water that has entered the aquifer since 1953), pre-modern (water that entered the aquifer prior to 1953 to tens of thousands of years ago), and mixed (mixtures of modern- and pre-modern-age waters). Arsenic concentrations were significantly greater in groundwater with pre-modern age classification than in groundwater with modern-age classification, suggesting that arsenic accumulates with groundwater residence time.

Scientific Investigations Report↗

Regional and other general factors bearing on evaluation of earthquake and other geologic hazards to coastal communities of southeastern Alaska

The great Alaska earthquake of March 27, 1964, brought into sharp focus the need for engineering geologic studies in seismically active regions. As a result, nine communities in southeastern Alaska were selected for reconnaissance investigations as an integral part of an overall program to evaluate earthquake and other geologic hazards in most of the larger Alaska coastal communities. This report gives background information on the regional and other general factors that bear on these evaluations. Southeastern Alaska, about 525 miles long and averaging about 125 miles in width, consists of a narrow mainland strip and numerous islands. For the most part, it is a region of rugged relief with numerous glaciers capping many of the higher mountainous areas and with long linear fiords forming the inland waterways. A maritime climate prevails with mild winters and cool summers. The southeastern part of the region receives the highest precipitation in the continental United States. Ketchikan, with a population of 6,994 in 1970, is the largest city. Geology and structure of the area are complex. Igneous, metamorphic, and sedimentary rocks crop out and range in age from Paleozoic to Tertiary. Surficial deposits of Pleistocene and Holocene age mantle many areas. All of southeastern Alaska, except probably the highest peaks, was covered by glacier ice advances of late Pleistocene age. Major deglaciation was well advanced by 10,000 years ago--a time which approximately marks the end of the Pleistocene and the beginning of the Holocene. There followed a period of warm climate called the Hypsithermal, which in southeastern Alaska began 7,000-8,000 years ago and ended about 4,800-3,500 years ago. Glaciers in most places receded back of their present positions. The Hypsithermal was followed by an interval (termed Neoglaciation) of cooler climate and resurgence of glacier ice which continues to the present, although most glaciers are now rapidly receding. During the past 10,000 years worldwide sea level has risen about 100 feet, but during the past 4,000 years it has risen only about 10 feet or about 0.03 inch per year. With sea level used as a datum, the amount of sea-level rise must be added to the apparent uplift of land for the time under consideration to determine the actual amount of land uplift. The widespread presence of emergent marine deposits, several hundred feet above sea level, demonstrates that the land in southeastern Alaska has been uplifted since the last major deglaciation. The greatest known has been uplifted since the last major deglaciation. The greatest known uplift is in the vicinity of Juneau where glaciomarine deposits are present 750 feet above present sea level. Part of southeastern Alaska is presently undergoing one of the most rapid rates of uplift of any place in the world. The fastest emergence is occurring in the Glacier Bay area where the land is being uplifted relative to sea level approximately 3.9 cm per year. Most or all of the uplift appears to be due to rebound as a result of deglaciation. Southeastern Alaska lies within the circum-Pacific earthquake belt, one of the world's greatest zones of seismic activity. During historic time, there have been five earthquakes in the region with magnitudes of 8 or greater, three with magnitudes of 7 to 8, eight with magnitudes of 6 to 7, more than 15 with magnitudes of 5 to 6, and about 140 recorded earthquakes with magnitudes smaller than 5 or of unassigned magnitudes. All of the earthquakes with magnitudes 8 or greater, and a large proportion of the others, appear to be related to the active Fairweather- Queen Charlotte Islands fault system or its western extension, the Chugach-St. Elias fault. Earthquake epicenters on the Denali fault system, the other major fault system in southeastern Alaska, are few in comparison. However, because high microearthquake activity has been recorded recently on this system and earthquakes of moderate size have occurred on some of its segments, the Denali fault system probably should not be dismissed as a relict fault system of no current tectonic importance. There are numerous other known faults, as well as lineaments that may be faults of varying degrees of tectonic activity in southeastern Alaska, adjacent Canada, and eastern Alaska. One of these elements is the Totschunda fault system, which connects with the Denali fault system in eastern Alaska; it has been very active during Holocene time but few historical earthquake epicenters appear to be related to it. Both historical seismicity and geologic conditions, such as frequency and recency of faulting, must be considered together to permit an assessment of the future earthquake probability of an area. Data are too few for both factors for an accurate evaluation to be made of earthquake probability in southeastern Alaska. However, information compiled in the form of strain-release and seismic-zone maps permit some generalizations. Thus, it is tentatively concluded that most, if not all, of southeastern Alaska should be placed in seismic zone 3, a zone in which earthquakes of magnitude greater than 6 will occur from time to time and where there may be major damage to manmade structures. Inferred effects from future earthquakes in southeastern Alaska include: (1) surface displacement along faults and other tectonic land-level changes, (2) ground shaking, (3) compaction, (4) liquefaction in cohesionless materials, (5) reaction of sensitive and quick clays, (6) water-sediment ejection and associated subsidence and ground fracturing, (7) earthquake-induced sub aerial slides and slumps, (8) earthquake induced subaqueous slides, (9) effects on glaciers and related features, (10) effects on ground water and stream flow, and (11) tsunamis, seiches, and other abnormal water waves. Because of the reconnaissance nature of our studies in the coastal communities and the sparsity of laboratory data on physical properties of geologic units in each area studied, the inferred effects must be largely empirical and generalized. Therefore, the inferences are based in large part upon the effects of past major earthquakes in Alaska and elsewhere, particularly upon the well-documented effects of the Alaska earthquake of March 27, 1964. Buildings, highways, bridges, tunnels, harbor facilities, pipelines, canals, and other manmade structures may be severely damaged or destroyed by fault displacement or related tectonic land-level changes in southeastern Alaska. Direct damage from fault rupture would be restricted virtually to structures built directly athwart the fault. In California and Nevada, fault rupture almost always accompanies shocks of magnitude 6.5 or greater. The Alaska earthquake of March 27, 1964, and the Chilean earthquake of May 22, 1960, dramatically illustrated the severe adverse effects that can result from uplift or subsidence over a wide area. The variable most responsible for the degree of shaking at any epicentral distance is the type of ground. Generally, shaking is considerably greater in poorly consolidated deposits than in hard bedrock, particularly if the deposits are water saturated. Severe shaking of alluvial deposits and manmade fill, with resultant heavy damage, is well documented from the records of many past earthquakes. Damage commonly has been heavy as a result of ground settlement caused by compaction of loose sediments by shaking during an earthquake. This has been especially true where compaction was accompanied by tectonic downdrop of land, such as occurred during the Chilean earthquake of 1960 and the Alaska earthquake of 1964. Loosely emplaced manmade fill, deltaic deposits, beach deposits, and alluvial deposits may be susceptible to compaction in southeastern Alaska during a severe earthquake. Liquefaction of sand and silt is a fairly common effect of large earthquakes. It was well illustrated at Niigata, Japan, during the earthquake of June 16, 1964, and resulted in extensive damage. When part of a sloping soil mass liquefies, the entire mass can undergo catastrophic failure and can flow as a high-density liquid. In southeastern Alaska, deltaic deposits probably would be most susceptible to liquefaction. Sensitive and quick clays, which lose a considerable part of their strength when shaken, commonly fail during an earthquake and become rapid earthflows. Extensive studies were made of the sensitivity of the Bootlegger Cove Clay at Anchorage because of the marked loss of shear strength and dramatic failures of the deposits during the Alaska earthquake of 1964. If similar sensitive clays are present in some places in southeastern Alaska, they most likely are in some of the emergent fine-grained marine deposits; supporting data to confirm their presence, however, are largely lacking. Records of some 50 major earthquakes show that in at least half of the instances water and sediment have been ejected from surficial deposits Water-sediment ejection and associated subsidence and ground fracturing commonly cause extensive damage to the works of man. Ejecta may fill basements and other low-lying parts of buildings. Agricultural land can be covered with a blanket of infertile soils, and small ponds can be filled or made shallow. In southeastern Alaska these phenomena are most likely to occur on valley floors, deltas, tidal flats, alluvial fans, swamps, and lakeshores. Earthquake-induced sliding on land generally is confined to steep slopes but may take place in fine-grained deposits on moderately to nearly flat surfaces if the deposits are subject to liquefaction. A large rockslide triggered by the Lituya Bay, Alaska, earthquake of July 10, 1958, generated a wave that surged up the opposite wall of the inlet to a record height of 1,740 feet. During the Hebgen Lake, Montana, earthquake of August 17, 1959, a spectacular rockslide plunged into the Madison River canyon, buried 28 people, dammed the river, and created a large lake. Earthquake-records are replete with accounts of sliding of surficial deposits during moderate to large earthquakes. Most or all of the general factors that favor subaerial landsliding are present in southeastern Alaska. Earthquake-induced subaqueous slides can produce adverse effects both nearshore and some distance offshore. Nearshore sliding may progress shoreward and destroy harbor facilities and other structures, commonly with substantial loss of life. Disastrous large submarine slides occurred along the fronts of deltas in Seward and Valdez during the Alaska earthquake of 1964. In similar fashion, the largest submarine slides in southeastern Alaska likely will be triggered along the larger delta fronts. Sliding farther offshore can constitute a threat to navigation because of changes in water depths. Also underwater sliding can break communication cables. Glaciers were not greatly affected by the Alaska earthquake of 1964 despite the fact that about 20 percent of the area that underwent strong shaking is covered by ice. In contrast, the cataclysmic avalanche of ice and rock that fell from a high glacier-covered peak in Peru during the earthquake of May 31, 1970, produced devastating effects downvalley on man and his works in the form of mudflows. Most towns in southeastern Alaska are sufficiently distant from glaciers so as not be to directly affected. Both the Alaska earthquake of 1964 and the Hebgen Lake, Montana, earthquake of 1959 significantly affected ground- and surface-water regimens. Water levels in some wells declined whereas in others flow increased. Some springs discharged at a rate three times as much as normal; flow of others decreased or stopped. Discharge of many streams increased markedly. Most or all of the effects described above could occur in parts of southeastern Alaska during future large earthquakes. Tsunamis, seiches, and other abnormal water waves associated with large earthquakes commonly cause vast property damage and heavy loss of life. Tsunami effects can be devastating to coastal areas as far as many thousands of miles from their generation source. Seiche effects generally are confined to inland bodies of water or to relatively enclosed coastal bodies of water. Abnormal waves generated by submarine sliding or by subaerial sliding into water generally produce only local effects but may be highly devastating. Tsunami waves resulting from the Chilean earthquake of 1960 inflicted extensive damage and loss of life on coastal communities throughout a large part of southern Chile, and significant runups and damage were recorded in many places throughout the Pacific Ocean area. The tsunami waves generated by the Alaska earthquake of 1964 struck with devastating force along a broad stretch of the Alaska coast and produced heavy property damage and loss of life as far away as Crescent City, Calif. Seiche waves generated by that earthquake reached runup heights of 20-30 feet on some lakes in Alaska, and water-level fluctuations were recorded on streams, reservoirs, lakes, and swimming pools in States bordering the Gulf of Mexico. Waves generated by submarine sliding struck violently at a number of places during or immediately after the quake and were the major cause of loss of life and damage to property. Slide-generated waves probably would have a higher destructive potential in southeastern Alaska than either tsunami waves or seiche waves because of their possibly higher local runups and because they can hit the shores almost without warning during or immediately after an earthquake. Nonearthquake-related geologic hazards, although generally far less dramatic than those related to earthquakes, tend to occur so much more frequently or persistently that their aggregate effects can be significant. Three kinds of geologic hazards of this type are discussed: (1) nonearthquake-induced landsliding and subaqueous sliding, (2) flooding, and (3) land uplift. The potential for nonearthquake-triggered landsliding in southeastern Alaska ranges widely from place to place. Past sliding generally furnishes the clue in the prediction of where and in what materials future sliding will occur. Fast-moving rockslides, debris slides, and mudflows can be expected to occur from time to time on steep slopes and be highly destructive to highways, power plants, pipelines, buildings, and other facilities located on a slope or at its base. Present slow downslope movement of talus can be expected to continue at the same general rate unless conditions are changed by man or there are climatic changes. Snow and debris avalanches can be especially hazardous during winter months. Long-inactive landslides may be triggered into renewed activity or new slides may be created by man-induced modifications. Accelerated slope erosion and debris flows may follow large-scale clearing and cutting of timber. Subaqueous sliding can be expected to occur periodically along fronts of deltas and on other oversteepened underwater slopes. Floods have been common in parts of southeastern Alaska because of heavy precipitation and rapid runoff from steep slopes with resulting heavy damage to roads and other facilities. Continued damage can be expected in the future unless more remedial measures are taken. Current uplift of land in southeastern Alaska, although probably not affecting man significantly in a short period of time, may have some adverse long-term effects. These long-term effects should be borne in mind when facilities such as docks and boat harbors are constructed on or near the shore, where there is a critical relation between height of land and water.

Alaska↗

Gear comparison study for sampling nekton in Barataria Basin marshes

This project was funded by the Louisiana Trustee Implementation Group (LA TIG) to support decisions related to investments in long-term monitoring. The LA TIG seeks to ensure long-term monitoring informs coastal restoration activities with the goal of sustaining and improving fisheries impacted by the Deepwater Horizon (DWH) Oil Spill. The project objective was to compare nekton catch across an estuarine gradient using different sampling gear with the goal of identifying trade-offs among nekton sampling approaches. To accomplish this objective, Louisiana Department of Wildlife and Fisheries (LDWF), The Water Institute of the Gulf (the Institute), Dynamic Solutions, LLC, Louisiana State University Agricultural Center (LSU AgCenter), and the U.S. Geological Survey (USGS) completed a field gear comparison study from 2018 to 2019. This work compared electrofisher and seine sampling at 12 fixed stations in Barataria Basin using data collected by LDWF. In addition, and in conjunction with LDWF monthly sampling, the same 12 fixed stations were sampled in May 2019 using a throw trap to compare nekton catch and assemblages collected with the throw trap, seine and electrofisher. LDWF has been conducting seine sampling since 1986, and seine data are used by the State of Louisiana to assess juvenile shrimp, crab and fish abundances, sizes and overall assemblages. In 2018, LDWF began conducting electrofisher sampling at 12 Barataria Basin seine stations in order to determine if the two gear types sample similar species and assemblages for potential future replacement of long-term seine sampling with electrofishing. Throw traps were included as they provide density estimates, which are ultimately the desired statistic used in modeling trophic webs, and are used in assessing habitat restoration outcomes. The project compared the nekton catch and assemblages collected using seine, electrofisher, and throw trap data from marsh edge habitats located across the estuarine gradient in Barataria Basin. Specifically, catch per unit effort (CPUE), species richness, species-specific total length (mm) distribution and nekton assemblages were compared between gear types. The first dataset was collected in May 2019 with throw trap (Appendix A), seine (LDWF data), and electrofisher (LDWF data) gear, and the second dataset (collected by LDWF) spanned 14 months of seine and electrofisher monthly sampling occurring from May 2018 through June 2019 at 12 stations in Barataria Basin. Key findings include that gear bias was not evident across the range of water quality conditions (salinity, temperature, o C, dissolved oxygen, mg L-1 , turbidity, NTU; Appendix B: scatter plots) captured during this pilot study, but differences in nekton catch per unit effort (CPUE) and assemblages were evident between gear types. However, those differences largely depended on the parameter examined. For example, the overall CPUE was highest for electrofishing, followed by seine, and then throw trap. When grass shrimp (the most abundant taxon collected) were removed from CPUE, the electrofisher and seine results were similar in CPUE. When CPUE was corrected for gear efficiency and total area sampled, the throw trap had the highest reported density of nekton sampled, followed by electrofisher and seine results. Electrofishing captured the highest number of species, which included more unique species compared to seine or throw trap catches, though all gear types captured at least one unique species. These highlight a need for caution in interpreting assemblage and density data when comparing datasets derived from different sampling methodologies. These key findings can help inform implementation and interpretation of long-term monitoring data in Louisiana as management decisions are made about coastal restoration projects to sustain and improve fisheries. There are trade-offs in selecting gear types for estuarine nekton monitoring of density, abundance, species richness, and assemblages. The table below (Table 1) summarizes some considerations when selecting gear types for long-term monitoring of estuarine nekton. In addition to biological and ecological considerations, other important considerations include cost, the labor required to conduct sampling, logistics, and potential uncertainties related to how effective each gear type is for sampling the wide variety of conditions found across Louisiana’s coastal habitats. For example, although electrofishing may capture higher CPUE, the equipment is more expensive to obtain and maintain compared to the other gear types. Most importantly, this table highlights differences in the nekton assemblages sampled by each gear type; this consideration is critical when designing the goals of a long-term monitoring program as it will inform how the data can be used and interpreted in the future. This report provides caveats, assumptions, and recommendations that can help support the Louisiana Coastal Protection and Restoration Authority (CPRA), LDWF and the LA TIG in comparing data from different gear types, and in making decisions for future monitoring. Findings from this study are limited to the range of water quality conditions occurring during these data collection events; these data and analyses could benefit from sampling across a wider range of water quality conditions, and collection of habitat structure and bottom type data which are not routinely collected but critically influence nekton. Further investigation examining how relative differences detected in key species abundances between gear types might impact ecosystem indicators and energetics in a modeled food web would provide valuable input to understand outputs of the Comprehensive Aquatic System Model for Barataria Basin, including the potential impacts of nekton monitoring decisions on food web models.

Louisiana↗

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↗