Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Engineering Analysis”

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.

1,052 records · Page 59Linked to original sources

Compounds of emerging concern detected in water samples from potable water and wastewater treatment plants and detected in water and bed-sediment samples from sites on the Trinity River, Dallas, Texas, 2009–13

The population in the Dallas-Fort Worth metropolitan area in northern Texas is rapidly growing, resulting in a rapid increase in the demand for potable water and an increase in the discharge of wastewater treatment plant effluent. An assessment of compounds of emerging concern (CECs) in samples collected at potable water and wastewater treatment plants in Dallas and downstream from Dallas in the Trinity River was completed by the U.S. Geological Survey in cooperation with the City of Dallas, Dallas Water Utilities. CECs are synthetic or naturally occurring chemicals that are not commonly monitored in the environment but can enter the environment and cause known or suspected adverse ecological or human health effects. CECs can enter the environment through nonpoint sources (for example, runoff) and point sources (for example, concentrated animal feeding operations and treated-effluent discharge from wastewater treatment plants), which can increase concentrations of CECs especially in highly populated areas. CECs include pharmaceuticals (prescription and nonprescription), steroidal hormones, stanols, sterols, detergents and detergent metabolites (hereinafter referred to as “detergents”), personal-use products, pesticides, polycyclic aromatic hydrocarbons (PAHs), flame retardants, plasticizers, and other organic compounds used in everyday domestic, agricultural, and industrial applications. The release of CECs to the environment went largely unrecognized until relatively recently. Increased loading of certain CECs to the environment, combined with advancements in laboratory analysis methods that resulted in appreciably lower detection levels, brought greater attention to the release of CECs. In addition, synthesis of new chemicals or changes in use and disposal of existing chemicals can create new sources of CECs. Some CECs are endocrine disrupting compounds (EDCs), which can elicit adverse effects on development, behavior, and reproduction of wildlife and can cause dysfunction of human and wildlife endocrine (hormone) systems. Results of studies in the United States and Europe indicate that CECs, their metabolites, and industrial, agricultural, and household wastewater products are present in the aquatic environment, water treatment plants, and septic systems. CECs, especially pharmaceuticals, are of interest because of their persistence, widespread use, and potential to cause adverse effects in humans and nontargeted organisms. There is also concern that some CECs and EDCs resist degradation of water treatment processes at potable water treatment plants (PWTPs) and wastewater treatment plants (WWTPs) and that treated-effluent discharge could contain compounds that negatively affect biota living in receiving waters. Therefore, CECs and EDCs are more likely to be detected in environmental samples collected near areas of high population density where treated effluent from WWTPs can contribute substantially to receiving waters. The U.S. Geological Survey, in cooperation with the City of Dallas, Dallas Water Utilities, evaluated the occurrence and concentrations of selected CECs in samples collected at PWTPs and WWTPs in Dallas and downstream from the Dallas-Fort Worth metropolitan area in the Trinity River, Texas, from August 2009 to December 2013. Water samples were collected at three PWTP sites, two WWTP sites, and five study sites on the Trinity River; all sites where samples were collected were in or downstream from Dallas. These water samples were analyzed for 120 CECs, including human-health pharmaceuticals (prescription and nonprescription), antibiotics, steroidal hormones, stanols, sterols, detergents, personal-use products (flavors and fragrances), pesticides and repellents, industrial wastewater compounds, disinfection compounds, PAHs, flame retardants, and plasticizers. Additionally, bed-sediment samples were collected at each of the five Trinity River sites. The bed-sediment samples were analyzed for 57 CECs. In general, the water treatment processes at PWTPs and WWTPs were effective at reducing detections and concentrations of CECs to undetectable levels or transforming the compounds into degradates that were not analyzed. There were 14 and 73 CECs detected in raw water and in untreated-influent water at PWTPs and WWTPs, respectively. Of these, 11 of the 14 CECs detected in raw-water samples and 44 of the 73 CECs detected in untreated-influent samples were not detected in finished water or in treated-effluent water samples, respectively, indicating that these compounds were removed or degraded to compounds that were not analyzed. Some CECs, however, are resistant to degradation and were detected in untreated and treated water at PWTPs and at WWTPs. The three CECs detected at PWTPs in raw-water and finished-water samples were tris(dichloroisopropyl)phosphate, benzophenone, and methyl salicylate. At WWTPs, 29 CECs were detected, including carbamazepine, sulfamethoxazole, 4-androstene-3,17- dione , 3-beta-coprostanol, acetyl-hexamethyl-tetrahydronaphthalene (AHTN), hexahydro-hexamethyl-cyclopenta-benzopyran (HHCB), 1,4-dichlorobenzene, tribromomethane, benzophenone, and tris(dichloroisopropyl)phosphate, in untreated and treated water, indicating that treatment processes likely did not remove or degrade these compounds. Of the 23 CECs detected in stream-water samples collected at 5 sites on the Trinity River in or near Dallas, 10 CECs (carbamazepine, sulfamethoxazole, caffeine, 3-beta-coprostanol, cholesterol, HHCB, benzophenone, triethyl citrate, tributyl phosphate, and tris(dichloroisopropyl)phosphate) were detected at all 5 sites. The 10 CECs detected in water samples collected at all 5 sites on the Trinity River were also detected in treated-effluent water at WWTPs. Eleven of the 57 targeted CECs were detected in bed-sediment samples collected at study sites on the Trinity River. Of these 11 CECs, only 2 (beta-sitosterol and cholesterol) were detected in bed-sediment samples at all 5 sites on the Trinity River. Nine of these 11 CECs were not detected in any water-column sample, likely because of the strong hydrophobic characteristics of these compounds. Results from water treatment plants indicate that the water treatment process is less effective for removing or degrading compounds that are engineered to be resistant to degradation. These results also indicate the presence of CECs and EDCs at locations upstream from PWTPs in Dallas. Results from Trinity River main-stem sites indicate that some compounds are naturally attenuated during transport, but a few are persistent throughout the study reach. Many CECs and EDCs are hydrophobic and were only detected in bed sediment, indicating multiple pathways through which CECs can persist in the environment. In general, concentrations of CECs in the Dallas-Fort Worth metropolitan area were similar to those found in metropolitan areas nationwide.

Texas↗

Evaluation of some software measuring displacements using GPS in real-time

For the past decade, the USGS has been monitoring deformation at various locations in the western United States using continuous GPS. The main focus of these measurements are estimates of displacement averaged over one day. Essentially, these consist of recording at 30 seconds intervals the carrier-frequency phase-data (equivalent to travel-time) between a GPS receiver and the GPS satellite network. In turn, these observations, which are converted to pseudo—ranges, are processed using one of the “research grade” programs (GIPSY, Zumberge et al., or GAMIT, wwwgpsg.mit.edu/~simon/gtgk) to estimate the position of the GPS receiver averaged over 24 hours. However, it is possible and desirable to estimate the position of the receiver (actually the antenna) more frequently and to do this within a few seconds of the time actual measurement (known as real-time). A recent example, the 2004 Magnitude 6, Parkfield, California earthquake, demonstrated that having GPS estimates of position more frequently than simply a daily average is required if one requires discrimination between co-seismic and post-seismic deformation (Langbein et al., 2006). The high-rate estimates of position obtained at Parkfield show that post-seismic deformation started less than one-hour after the mainshock and that this deformation was roughly the same magnitude as the co-seismic deformation. The high-rate solutions for Parkfield were done by others including Yehuda Bock at UCSD and Kristine Larson at U. of Colorado, but not the USGS. The Parkfield experience points out the need for an in-house capability by the USGS to be able to accurately measure co-seismic displacements and other rapid, deformation signals using GPS. This applies to both the Earthquake and Volcano Hazard programs. Although at many locations where we monitor deformation, we have strainmeters and tiltmeters in addition to GPS which, in principle, are far more sensitive to rapid deformation over periods of less than a day (Langbein and Bock, 2004). But, not all locales include strain and tiltmeters. Thus, having the capability to extract signals with periods of less than a day is desirable since the distribution of GPS is more extensive than strain and tilt. At both Parkfield and Long Valley, the USGS has been using other software packages to process the GPS data at sub-daily intervals and in real-time. The underlying goal of these types of measurements is to detect any deformation event as it evolves; the 24 hour processing might not provide timely results if such a deformation event is precursory to a geologic hazard (an earthquake for Parkfield and either a volcanic event or an earthquake for Long Valley). In Long Valley, We use the software package called 3DTracker (http://www.3dtracker.com, http://www.condorearth.com) to estimate the changes of in position of a remote site relative to a “fixed” site. The 3DTracker software uses double difference GPS code measurements and receiversatellite-time triple differences from one epoch to the next of the GPS phase data (a proxy for travel-time measurements) and employs a Kalman filter to obtain stability in the estimate of position. That is, the estimate of the current position depends upon the estimate of the prior position. Hence, a time series of position looks fairly smooth depending upon the coefficient selected for the Kalman filter. With triple differences, the sometimes troublesome initial integer cycle ambiguity terms cancel (number of wavelengths between the receiver and each satellite), but only the incremental change in position is calculated. This triple difference Kalman filter solution is slow to converge and less accurate than a double difference (e.g., RTD, Track) solution, but it is robust and computationally efficient (Remondi and Brown, 2000). 3D-Tracker allows use of various single-frequency and dual-frequency GPS phase and code observables including the ionospheric-free combinations (known as LC or L3 and P(L3)) formed from an linear combination of the L1 and L2 carrier phase and code data. The lowest noise observable is the L1 carrier, but it is biased by ionospheric refraction that has amplitudes of about 1 to 10 ppm. This results in a systematic scale error in the relative positions. The L3 phase noise is about 3 times greater than the L1 phase noise, but it is generally used to solve for all but the shortest baselines (< 5 km). In addition, the software does output the position changes is a standard format that can be used for other analysis. At Parkfield, we use the software package called RTD (http://www.geodetics.com). The RTD software has been described in the literature (Bock et al., 2000) but basically, it estimates the position without the constraint of a Kalman filter. It uses double differences (in our studies the LC or ionospheric free observable is used) and the integer ambiguities are resolved independently for each 1-second measurement; Most GPS software that use double-differences require several epochs of measurements to resolve the integer ambiguities. The data files use a proprietary format and can not be read by me or others; rather, Yehuda Bock at UCSD (and author of RTD) translates these files into a standard format that can be read by me. Recently, Tom Herring of MIT has modified the GAMIT software to process kinematically GPS data (www-gpsg.mit.edu/~simon/gtgk/tutorial/Lecture_13.pdf). At this time, the software, known as TRACK, does not process the observations in real-time. Consequently, the latency between the time of the observation and the time when a position estimate is available depends upon the frequency that the data are downloaded and the speed of actually processing the observations; there could be a delay of an hour or two before the a position estimates are available. Unlike RTD and 3DTracker, TRACK comes with GAMIT (which is distributed freely) and is currently operating in a test mode at the USGS office in Pasadena. The LC or ionosphere free observable is used in our TRACK solutions. JPL has a version of their GIPSY software called “Real-time GIPSY (RTG)” (gipsy.jpl.nasa.gov/orms/rtg), which, like TRACK, can process the pseudo-range data “off—line”. However, this software is not freely distributed. Instead, at least one company, NAVCOM, has teamed with JPL to integrate RTG with GPS receivers and telemetry that yields positions in realtime. Kristine Larson of University of Colorado has modified the original GIPSY to estimate positions kinematically. Again, like TRACK, the positions are estimated off—line. Much of her research is described in Larson et al. (2003), and Choi et al. (2004). For Long Valley, out of the 17 GPS sites, we monitor 5 baselines within the caldera at 5 second intervals relative to the Bald Mountain site at the edge of the caldera using 3DTracker. The baseline measurement using 3DTracker consists of determination of the 3 dimensional positions of the 5 remote points (GPS receivers) relative to a GPS site at Bald. A second, independent system collects and downloads once a day the 30-second data used for the 24-hour solutions for the 12 sites not monitored with 3DTracker. For the sites monitored with 3DTracker, the pseudo—range data are decimated to 30 seconds and converted to a form used for the 24-hour solutions. Both sets of telemetry employ 900 MHz spread spectrum radios which require line of site between all of the links. The telemetry for the 3DTracker sites require a dedicated radios at each end and intermediate repeaters as needed, while the telemetry required for the other sites use a single master radio, repeaters as needed, and a radio at each remote site. (The 5 sites being monitored with 3DTracker require 13 radios.) At Parkfield, RTD is used to measure the position changes all 12 baselines at 1 second intervals relative to a site, Pomm, adjacent to the San Andreas Fault. The complete RTD package (hardware and software) collects all of the data and determines the position of each site relative to Pomm. In addition, the system stores both the 1-second and 30-second pseudo-range data for later downloading which are ultimately used in the 24-hour solutions. To do this, each site has a 2.4 GHz radio and a telemetry buffer. The telemetry buffer holds 24-hours of data (in the event that the telemetry link is broken) and converts the RS232 data stream from the GPS receiver into a form compatible with an IP (Internet protocol) network connection. In contrast with the Long Valley system, the telemetry link for GPS at Parkfield consists of a single radio at each remote sites and a single radio at the central site. Although position estimates are produced within 1-second of the observations, these results are not immediately available because there is no high speed Internet connection to Parkfield. Instead, the data are stored on a removable disk and sent to UCSD once per month. Below, I describe the results of a simple experiment to examine the response of some of these systems to simulated deformation that could be an analogue of a tectonic or volcanic event. In many engineering applications, the system response is tested by inputting a step to the system and measuring the output of the system. Essentially, this is what I've done. The experiment described below moves the GPS antenna from its original position to a new position within 1 second; the software tracks the translation. These measurements were conducted in August 2004 with the RTD software at Parkfield, and twice in Long Valley. The first Long Valley test was conducted in September 2004 using 3DTracker on a single baseline. The test was repeated in September 2005 using 3DTracker on two baselines and, importantly, saving the RINEX files of the data so that the data could be replayed through 3DTracker using other options in the program and, using other software packages including TRACK. In addition, we observed a short-term event at the Three Sisters volcano in Oregon. This event was snow melt at a remote GPS site which gave an apparent 15 cm displacement in vertical in less than one-day. 3DTracker is used to monitor this site, and the event was captured with this software. In addition, with the assistance of others, I got additional estimates of position using other software packages; those results are presented. Finally, the precision of both 3DTracker and RTD are compared using a power spectrum. Those results would suggest that 3DTracker using appropriate Kalman filter coefficients would have better precision than RTD; instead, the lower noise level from 3DTracker is a result of smoothing from the Kalman filter. Given the results described in this report, high-rate GPS is certainly capable of accurately measuring displacements of 1 centimeter with a high degree of statistical confidence. Plotting these results show that the time of the displacement can be visually determined to that of the sampling interval of the data. However, especially with small amplitude signals, any of the software packages can yield erroneous deformation “signals” that are either due excess travel-time of the GPS carrier frequency from multipath or a limitation in the software. Thus, the time series of displacements must be viewed with caution and knowledge of external circumstances that might cause a change in position. The casual reader should continue with the next section describing the methods then jump to the last two sections for the discussion and conclusions. I have made some recommendations there.

Open-File Report↗

Davis Pond freshwater prediversion biomonitoring study: freshwater fisheries and eagles

In January 2001, the construction of the Davis Pond freshwater diversion structure was completed by the U.S. Army Corps of Engineers. The diversion of freshwater from the Mississippi River is intended to mitigate saltwater intrusion from the Gulf of Mexico and to lessen the concomitant loss of wetland areas. In addition to the freshwater inflow, Barataria Bay basin would receive nutrients, increased flows of sediments, and water-borne and sediment-bound compounds. The purpose of this biomonitoring study was, therefore, to serve as a baseline for prediversion concentrations of selected contaminants in bald eagle ( Haliaeetus leucocephalus ) nestlings (hereafter referred to as eaglets), representative freshwater fish, and bivalves. Samples were collected from January through June 2001. Two similarly designed postdiversion studies, as described in the biological monitoring program, are planned. Active bald eagle nests targeted for sampling eaglet blood (n = 6) were generally located southwest and south of the diversion structure. The designated sites for aquatic animal sampling were at Lake Salvador, at Lake Cataouatche, at Bayou Couba, and along the Mississippi River. Aquatic animals representative of eagle prey were collected. Fish were from three different trophic levels and have varying feeding strategies and life histories. These included herbivorous striped mullet ( Mugil cephalus ), omnivorous blue catfish ( Ictalurus furcatus ), and carnivorous largemouth bass ( Micropterus salmoides ). Three individuals per species were collected at each of the four sampling sites. Freshwater Atlantic rangia clams ( Rangia cuneata ) were collected at the downstream marsh sites, and zebra mussels ( Dreissena spp.) were collected on the Mississippi River. The U.S. Geological Survey (USGS) Biomonitoring of Environmental Status and Trends (BEST) protocols served as guides for fish sampling and health assessments. Fish are useful for monitoring aquatic ecosystems because they accumulate pesticides and other contaminants. Biomarker data on individual fish, generated at the USGS National Wetlands Research Center (Lafayette, La.), included percent white blood cells in whole blood, spleen weight to body weight ratio, liver weight to body weight ratio, condition factor, splenic macrophage aggregates, and liver microsomal 7-ethoxyresorufin-o-deethylase (EROD) activity. Fish age was estimated by comparing total lengths with values from the same species in the Southeast United States as determined from the literature. Contaminant analyses were coordinated by the U.S. Fish and Wildlife Service (USFWS) Analytical Control Facility (Laurel, Md.), where residues of organochlorine (OC) pesticides, total polychlorinated biphenyls (PCBs), polycyclic aromatic hydrocarbons (PAHs), aliphatic hydrocarbons (AHs), and trace elements were determined. The organic contaminant data were generated at the Mississippi State University Chemical Lab (Mississippi State, Miss.), and the inorganic contaminant data were generated by the Texas A&M University Geochemical and Environmental Research Group (College Station, Tex.). Statistical tests were performed to assess relationships among contaminants, fish age, fish species, and collection sites. Trends in interspecific differences among fish in concentrations of contaminants were noted. Striped mullet (hereafter mullet) frequently displayed the highest chemical concentrations. Levels of contaminants were generally higher in samples obtained from the Mississippi River than in those collected from the diversion area and were higher in mussels and clams (hereafter bivalves) than in fish. Because the Mississippi River sampling site for mullet and largemouth bass was downriver of the structure and south of New Orleans and the catfish site was upriver, the downriver data may not be directly reflective of the results from the receiving waters at the Davis Pond structure. Compared to the Caernarvon freshwater prediversion study in 1990 that assessed possible influx of contaminants with the freshwater diversion, contaminant levels in fishes and bivalves in this study were generally lower, yet three nontoxic inorganic elements in Davis Pond fish samples exhibited ranges of concentrations that were more than two times higher than did those from Caernarvon. Levels in bivalves were different between diversions but about equal in the numbers of trace elements showing high levels per location. Contaminant values were compared to those listed in various literature and agency sources, both regional and national, including the National Contaminant Biomonitoring Program (NCBP), in which the 85th percentile and above represents what is considered to be an elevated contaminant concentration and cause for concern. Generally, bivalves were at the high end of their ranges for both organic and inorganic contaminants. In this study, OCs were detectable in 67 percent of fish from the Mississippi River site, ranging from 0.15 to 1.09 μg/g wet weight (ww) or fresh weight (fw), and in 11 percent of the fish from the marsh sites, ranging from 0.06 to 0.612 μg/g ww. Bivalves from the Mississippi River had OC levels of 0.096 μg/g ww, whereas none were detectable in bivalves at the marsh sites. In this study, p,p ’-dichlorodiphenyldichloroethylene ( p,p’ -DDE) (a biodegradation product of DDT [dichlorodiphenyl trichloroethane]) and total PCBs were the most frequently detected OCs and were primarily from the Mississippi River. For total OC content, using adjusted least squares means, some significant interactions were noted between fish species and sites. PAHs were detected in aquatic animals at all sites (range of 0.017–17.534 μg/g ww), as were AHs (range of 0.423–4.549 μg/g ww); the highest levels of PAHs and AHs were found in bivalves from the Mississippi River. When analysis of variance (α = 0.05) was performed with data from aquatic animals, there were only two significant relationships between PAHs, AHs, and OCs between species, site, and age or the interaction among these variables. There was an interaction between fish species and n-decane (an AH) in that mullet and largemouth bass had significantly higher levels than did catfish ( P = 0.0175). When general linear means were used to investigate associations of inorganic contaminants among fish species, site, and age or any interactions among these variables, no significant results were noted for arsenic, cadmium, lead, beryllium, boron, molybdenum, or nickel. The range of mercury in fish in this study was 0.04–0.14 μg/g ww (0.14– 0.48 μg/g dry weight [dw]), with the most elevated levels detected in predatory largemouth bass at the sampling point farthest downstream from the structure and within the marsh area. Mercury was positively correlated with fish age ( P = 0.0152), where levels were estimated to increase 0.0253 parts per million (ppm) dw per year. In the Mississippi River, catfish showed significantly higher levels of mercury than did mullet or largemouth bass ( P = 0.00167). Among fish species, mullet displayed the highest levels in fish of aluminum, barium, manganese, and iron, all considered to have low toxicity in hydrologic systems. An interaction between fish and site was seen with aluminum ( P = 0.0031), where concentrations in mullet were significantly higher in the Mississippi River than at the other sites, as was also seen with barium ( P = 0.0009), chromium ( P = <0.0001), manganese ( P = 0.0004), strontium ( P = 0.0074), vanadium ( P = 0.0156), and zinc ( P = 0.0059). For iron ( P = 0.0.0001), mullet and largemouth bass at both the Mississippi River and Lake Salvador showed higher levels than did catfish, and these two species showed higher levels at two of the four sites. An interaction between fish and site was also seen with chromium ( P = <0.0001) in that concentrations in mullet were significantly higher in the Mississippi River than at the other sites, as was also seen with strontium ( P = 0.0074), vanadium ( P = 0.0156), and zinc ( P = 0.0059), metals for which deleterious effects have been demonstrated in other ecosystems. The NCBP program lists the 85th percentile for zinc at 34.2 μg/g fw (117.9 μg/g dw). In the Davis Pond prediversion biomonitoring study (hereafter the current study), one fish (MUL31RIVER, fish ID 8) showed values higher than that (125.4 μg/g dw or 37.54 μg/g ww), and the Mississippi River bivalve sample (MUSSRIVER) had a value of 140 μg/g dw (41.2 μg/g ww). In the current study, approximately 86 percent of the fish had measurable selenium levels, yet none reached the 85th percentile. The 85th percentile for selenium from the NCBP was 0.73 μg/g ww. Significantly higher levels of selenium were seen in mullet than in largemouth bass and catfish ( P = 0.0023). The NCBP 85th percentile for lead is 0.22 μg/g ww (0.76 μg/g dw). In the current study, the range of concentrations of lead was as much as 18.3 ppm dw (MUL31RIVER, fish ID 8), with the three most elevated values (range of 3.46–5.31 μg/g ww) coming from mullet from the Mississippi River. Biomarker data are measurable and directly reflect the condition of the animal, and measuring more than one biomarker in an individual increases confidence in health assessments. In the current study, biomarkers included macrophage aggregates (MAs), liver (hepatosomatic index [HSI]) and spleen (splenosomatic index [SSI]) weight to body weight ratios, percent white blood cells (WBCs) in whole blood, and condition factor. Few significant differences were noted with any of the biomarkers between sites, and there were no relationships between species and sites. For improved use of biomarker assessments, an increase in fish sample size would be useful for postdiversion sampling, as would comparisons of fish of the same sex and reproductive condition. During the current study, success for eagle nests in the diversion area and reference sites was similar as determined by numbers of nestlings fledged. When temperatures were below average during winter 2000, nests in both regions similarly failed. At each nest, the primary evidence of food items was small mammals. Eaglets (n = 6) generally appeared healthy, and whole blood concentrations of organic contaminants exceeded detection limits with three incidences of p,p’ -DDE (0.002–0.006 μg/L ww) and one incidence of oxychlordane (0.002 μg/L ww). The levels of p,p’ -DDE were well below those that have been inversely correlated with productivity and success rates of nesting bald eagles on a regional scale. The low values found in the whole blood samples for OC pesticides and PCBs were even lower when corrected for plasma volume. Aluminum values were 3.66 and 5.75 μg/L in two samples, zinc ranged from 5.21 to 6.77 μg/L ww in six samples, and silicon ranged from 1.7 to 4.6 μg/L in four samples. Selenium was detectable in each bird with the range at 0.332–0.566 μg/L ww, and strontium ranged from 0.0581 to 0.0975 μg/L ww. Mercury was detectable in blood samples from each bird and ranged from 0.0254 to 0.0845 μg/L ww, whereas lead was detectable in four samples and ranged from 0.0042 to 0.0136 μg/L ww. Although no detectable levels of total PCBs were found (also correlated with decreased reproductive productivity), 70 percent of the aquatic animals from the Mississippi River contained total PCBs (range 0.13–0.79 μg/L), whereas only about 7 percent of the aquatic animals sampled from the marsh area contained PCBs. Suggestions for postdiversion sampling include lowering the analytical detection limit for some metals, sampling aquatic animals over the course of a single season, obtaining a higher sample number of mature fish of one species (for example, blue catfish) within a range of total lengths for biomarker analyses, obtaining otoliths for estimating fish ages, assessing dioxins in eaglet blood, examining triazines in water, and obtaining all Mississippi River fish samples as close to the Davis Pond structure intake as possible. Because contaminants found in blood of eaglets reflect their prey species and because of the contaminant levels found in fish in the current study, eaglets may not be consuming primarily these species; therefore, obtaining juvenile nutria ( Myocastor coypus ) or turtle species for contaminant analyses might be considered, as well as collecting greater blood volume and using plasma to measure OCs and PCBs. Data obtained postdiversion will be compared with prediversion data to monitor changes.

Louisiana↗

Evaluation of groundwater levels in the South Platte River alluvial aquifer, Colorado, 1953-2012, and design of initial well networks for monitoring groundwater levels

The South Platte River and underlying alluvial aquifer form an important hydrologic resource in northeastern Colorado that provides water to population centers along the Front Range and to agricultural communities across the rural plains. Water is regulated based on seniority of water rights and delivered using a network of administration structures that includes ditches, reservoirs, wells, impacted river sections, and engineered recharge areas. A recent addendum to Colorado water law enacted during 2002-2003 curtailed pumping from thousands of wells that lacked authorized augmentation plans. The restrictions in pumping were hypothesized to increase water storage in the aquifer, causing groundwater to rise near the land surface at some locations. The U.S. Geological Survey (USGS), in cooperation with the Colorado Water Conservation Board and the Colorado Water Institute, completed an assessment of 60 years (yr) of historical groundwater-level records collected from 1953 to 2012 from 1,669 wells. Relations of "high" groundwater levels, defined as depth to water from 0 to 10 feet (ft) below land surface, were compared to precipitation, river discharge, and 36 geographic and administrative attributes to identify natural and human controls in areas with shallow groundwater. Averaged per decade and over the entire aquifer, depths to groundwater varied between 24 and 32 ft over the 60-yr record. The shallowest average depth to water was identified during 1983-1992, which also recorded the highest levels of decadal precipitation. Average depth to water was greatest (32 ft) during 1953-1962 and intermediate (30 ft) in the recent decade (2003-2012) following curtailment of pumping. Between the decades 1993-2002 and 2003-2012, groundwater levels declined about 2 ft across the aquifer. In comparison, in areas where groundwater levels were within 20 ft of the land surface, observed groundwater levels rose about 0.6 ft, on average, during the same period, which demonstrated preferential rise in areas with shallow groundwater. Approximately 29 percent of water-level observations were identified as high groundwater in the South Platte River alluvial aquifer over the 60-yr record. High groundwater levels were found in 17 to 33 percent of wells examined by decade, with the largest percentages occurring over three decades from 1963 to 1992. The recent decade (2003-2012) exhibited an intermediate percentage (25 percent) of wells with high groundwater levels but also had the highest percentage (30 percent) of high groundwater observations, although results by observations were similar (26-29 percent) over three decades prior, from 1963 to 1992. Major sections of the aquifer from north of Sterling to Julesburg and areas near Greeley, La Salle, and Gilcrest were identified with the highest frequencies of high groundwater levels. Changes in groundwater levels were evaluated using Kendal line and least trimmed squares regression methods using a significance level of 0.01 and statistical power of 0.8. During 2003-2012, following curtailment of pumping, 88 percent of wells and 81 percent of subwatershed areas with significant trends in groundwater levels exhibited rising water levels. Over the complete 60-yr record, however, 66 percent of wells and 57 percent of subwatersheds with significant groundwater-level trends still showed declining water levels; rates of groundwater-level change were typically less than 0.125 ft/yr in areas near the South Platte River, with greater declines along the southern tributaries. In agreement, 58 percent of subwatersheds evaluated between 1963-1972 and 2003-2012 showed net declines in average decadal groundwater levels. More areas had groundwater decline in upgradient sections to the west and rise in downgradient sections to the east, implying a redistribution of water has occurred in some areas of the aquifer. Precipitation was identified as having the strongest statistically significant correlations to river discharge over annual and decadal periods (Pearson correlation coefficients of 0.5 and 0.8, respectively, and statistical significance defined by p-values less than 0.05). Correlation coefficients between river discharge and frequency of high groundwater levels were statistically significant at 0.4 annually and 0.6 over decadal periods, indicating that periods of high river flow were often coincident with high groundwater conditions. Over seasonal periods in five of the six decades examined, peak high groundwater levels occurred after spring runoff from July to September when administrative structures were most active. Between 1993-2002 and 2003-2012, groundwater levels rose while river discharge decreased, in part from greater reliance on surface water and curtailed pumping from wells without augmentation plans. Geographic attributes of elevation and proximity to streams and rivers showed moderate correlations to high groundwater levels in wells used for observing groundwater levels (correlation coefficients of 0.3 to 0.4). Local depressions and regional lows within the aquifer were identified as areas of potential shallow groundwater. Wells close to the river regularly indicated high groundwater levels, while those within depleted tributaries tended to have low frequencies of high groundwater levels. Some attributes of administrative structures were spatially correlated to high groundwater levels at moderate to high magnitudes (correlation coefficients of 0.3 to 0.7). The number of affected river reaches or recharge areas that surround a well where groundwater levels were observed and its distance from the nearest well field showed the strongest controls on high groundwater levels. Influences of administrative structures on groundwater levels were in some cases local over a mile or less but could extend to several miles, often manifesting as diffuse effects from multiple surrounding structures. A network of candidate monitoring wells was proposed to initiate a regional monitoring program. Consistent monitoring and analysis of groundwater levels will be needed for informed decisions to optimize beneficial use of water and to limit high groundwater levels in susceptible areas. Finalization of the network will require future field reconnaissance to assess local site conditions and discussions with State authorities.

Colorado↗

Effect of storms on barrier island dynamics, Core Banks, Cape Lookout National Seashore, North Carolina, 1960-2001

The effect of storms on long-term dynamics of barrier islands was evaluated on Core Banks, a series of barrier islands that extend from Cape Lookout to Okracoke Inlet in the Cape Lookout National Seashore, North Carolina. Shoreline and elevation changes were determined by comparing 77 profiles and associated reference markers established by the U.S. Army Corps of Engineers (USACE) on Core Banks from June 1960 to July 1962 to a follow-up survey by Godfrey and Godfrey (G&G) in 1971 and a survey by the Department of Geology at East Carolina University (ECU) in 2001, in which 57 of the original 77 profiles were located. Evaluation of the baseline data associated with the USACE study supplies an important record of barrier island response to two specific storm events—Hurricane Donna in September 1960 and the Ash Wednesday extra-tropical cyclone in March 1962. The 1962 USACE survey was followed by 9 years characterized by no major storms; this low-energy period was captured by the G&G survey in 1971. The G&G survey was followed by 22 years characterized by occasional small to moderate storms. Starting in 1993, however, and continuing through 1999, the North Carolina coast experienced a major increase in storm activity, with seven major hurricanes impacting Core Banks. Both the USACE 1960–1962 and G&G 1962–1971 surveys produced short-term data sets that reflected very different sets of weather conditions. The ECU 2001 survey data were then compared with the USACE 1960 survey data to develop a long-term (41 years) data set for shoreline erosion on Core Banks. Those resulting long-term data were compared with the long-term (52 years) data sets by the North Carolina Division of Coastal Management (NCDCM) from 1940–1992 and 1946–1998; a strong positive correlation and very similar rates of average annual erosion resulted. However, the ECU and NCDCM long-term data sets did not correlate with either of the USACE and G&G short-term survey data and had very different average annual erosion rates. The average annual long-term rate of shoreline erosion for all of Core Banks and for both the ECU 1960–2001 and the NCDCM 1946–1998 surveys was -5 feet per year (ft/yr). These long-term rates of shoreline recession are in strong contrast with the short-term, storm-dominated rates of shoreline erosion for all of Core Banks developed by the USACE 1960–1961 and USACE 1961–1962 surveys, which have average annual erosion rates of -40 ft/yr and -26 ft/yr, respectively, and range from -226 feet (ft) to +153 ft. The combined short-term, storm-dominated shoreline erosion rate for the USACE surveys (1960–1962) was -36 ft/yr. In contrast, the average annual short-term, non-stormy period G&G 1962–1971 survey demonstrated shoreline accretion for all of Core Banks with an average annual rate of +12 ft/yr. In general, North Core Banks has higher erosion and accretion rates than South Core Banks. In the 1961 survey, the USACE installed 231 reference markers (RM-0 is closest to the ocean and RM-2 is farthest from the ocean) along the 77 profiles, as well as 33 reference markers labeled RM-4, RM-6, and RM-8 in the wider portions of the islands. The G&G survey recovered a total of 141 reference markers (61 percent), and the ECU survey recovered a total of 83 reference markers (36 percent) of the RM-0, RM-1, and RM-2 markers. The average ground elevation measured by the USACE in 1961 was RM-0 = +5.8 ft, RM-1 = +5.2 ft, and RM-2 = +4.8 ft. The G&G 1970 survey measured average ground elevations of RM-0 = +6.7 ft, RM-1 = +6.4 ft, and RM-2 = +6.1 ft, and the average ground elevation measured by ECU in 2001 was RM-0 = +10.1 ft, RM-1 = +9.1 ft, and RM-2 = +8.5 ft. The latter numbers represent approximately an overall 72-percent increase in island elevation from 1961 to 2001. Based on aerial photographic time-slice analyses, it is hypothesized that this increase in island elevation occurred during the post-1962 period with storm overwash systematically raising the island elevation through time, which in turn led to decreased numbers of overwash events. The latter processes and responses in turn led to a substantial increase in vegetative growth on the barrier island, as well as submerged aquatic vegetation on the back-barrier sand shoals. Integration of the USACE, G&G, ECU, and NCDCM shoreline erosion data for Core Banks shows several important points about shoreline recession. (1) The ECU and NCDCM data sets demonstrate that there is an ongoing net, long-term, but small-scale shoreline recession associated with Core Banks; (2) the USACE short-term data sets demonstrate that processes associated with individual storm events or sets of events produce extremely large-scale changes that include both erosion and accretion; (3) the short-term, non-stormy period data set of G&G demonstrates that if given enough time between storm events, barriers can rebuild to their pre-storm period conditions; and (4) the post-storm response generally tends to approach the pre-storm location, but rarely reaches it before the next storm or stormy period sets in. The result is the net long-term change documented by both the ECU 1960–2001 and NCDCM 1946–1998 Core Banks data sets that resulted in erosion rates ranging from 0 to -30 ft/yr with net annual average recession rates of -5 ft/yr. Analysis and comparison of these data sets supply important information for understanding the dynamics and responses of barrier island systems through time. In addition, the results of the present study on Core Banks supply essential process-response information that can be used to design and implement management plans for the Cape Lookout and Cape Hatteras National Seashores and for other seashores in the U.S. National Park Service system.

North Carolina↗

Evaluation of the Bushy Park Reservoir three-dimensional hydrodynamic and water-quality model, South Carolina, 2012–15

The Bushy Park Reservoir is a relatively shallow impoundment in southeastern South Carolina. The reservoir, located under a semi-tropical climate, is the principal water supply for the city of Charleston, South Carolina, and the surrounding areas including the Bushy Park Industrial Complex. Although there was an adequate supply of freshwater in the reservoir in 2022, water-quality concerns are present over taste-and-odor and saltwater-intrusion issues. From 2013 to 2015, the U.S. Geological Survey (USGS), in cooperation with the Charleston Water System, engaged in a multi-year study of the hydrology and hydrodynamics of Bushy Park Reservoir to better understand factors affecting water-quality conditions in the reservoir. As part of this study, Charleston Water System worked with Tetra Tech, Inc., a consulting and engineering firm, to develop a Bushy Park Reservoir hydrodynamic and water-quality modeling framework, built upon earlier efforts by both Tetra Tech and the U.S. Army Corps of Engineers. At the completion of the new modeling framework, the USGS was requested to evaluate the calibrated hydrodynamic and water-quality model. The Bushy Park Reservoir Environmental Fluid Dynamics Code (EFDC) model was calibrated for the time period from January 1, 2012, to December 31, 2015. The general modeling approach for the newly revised modeling framework, as briefly detailed in this report, was developed with EFDC. The EFDC is a grid-based modeling package that can simulate three-dimensional flow, transport, and water quality in surface-water systems. This report evaluated the capacity of Tetra Tech’s Bushy Park EFDC model to simulate water discharge, water circulation, surface elevations, temperature, salinity, and other water-quality parameters. The USGS model review focused specifically on the following criteria: (1) determine if the model, with additional effort, could be developed into an adequate planning tool for Bushy Park Reservoir; (2) assess the capacity of the model to specifically address water-quality issues in the reservoir related to taste-and-odor and saltwater intrusions; and, (3) evaluate three preliminary water-management scenarios related to reduced water withdrawals in the reservoir and the effect on saltwater intrusion. Overall, the model was able to simulate discharge, flow velocity, and water-surface elevations with generally good agreement between the simulated and measured values. Specifically, the model was able to demonstrate good agreement for discharge at two USGS continuous discharge locations (USGS station 02172002; USGS station 02172040), with Wilmott index of agreements of 0.86 and 0.75, respectively. A total of seven USGS streamgages, located on the West Branch of the Cooper River, Durham Canal, and the Cooper River, were available for water-surface elevations, with index of agreements ranging from 0.74 to 0.99. However, model-simulated water-surface elevation ranges were appreciably high (compared to measured ranges) for two locations near Pinopolis Dam, farthest upstream on the West Branch of the Cooper River. This result may indicate that too much simulated tidal energy propagated through the model domain. For water temperature, 16 calibration stations were available for at least part of the 4-year simulation. The index of agreement range for temperature comparisons was from 0.95 to 1.00, indicating excellent agreement between the measured and simulated results. One of the primary future applications for the Bushy Park Reservoir EFDC model is to determine the extent of saltwater intrusions. A wide range in the salinity prediction quality was simulated with the model. The prediction quality ranged from an index of agreement of 0.15 at Cooper River approximately 2.75 miles southeast of the Tee, South Carolina, to 0.92 at West Branch Cooper River near Moncks Corner, South Carolina. Although the model did not accurately simulate some of the larger salinity deviations resulting during individual hydrologic events, the seasonal salinity trends were adequately simulated with the model during the study period (2012–15). Therefore, it may be difficult to simulate extreme hydrologic events, such as during large storms, where high salinity water is exchanged with Bushy Park Reservoir. There was agreement in model simulation with the measured data either on the quantitative index of agreement values or qualitative agreement in the seasonal salinity data trends. For water quality, the index of agreement values were generally low for total nitrogen, ammonia, nitrate, total Kjeldahl nitrogen, total phosphorus, and orthophosphate. Although general trends were adequately simulated at specific stations, particularly for Bushy Park Reservoir, the model-simulated fit was low across all the constituents described above with index of agreements usually below 0.50. A limitation for simulating nutrient concentrations across the model domain was the lack of characterization for the constituents directly entering Bushy Park Reservoir, or the lack of data directly attributed to the boundary condition (for example, the Cooper River). The other two calibrated water-quality constituents (besides the nutrients mentioned above) were dissolved oxygen and chlorophyll a . Dissolved oxygen varied from index of agreement values from 0.58 to 0.94 for 11 stations, generally indicating agreement with the available measured data. Chlorophyll a , calibrated for seven stations, had a wider range from 0.11 to 0.74 for the index of agreement. With the current modeling framework, taste-and-odor events, related to cyanobacterial blooms, cannot be directly simulated. However, indirect estimates of cyanobacteria concentrations may be obtained by using the chlorophyll a model outputs, which represent total phytoplankton biomass, and the phytoplankton biovolume data by group (diatoms, green algae, cyanobacteria and others) collected from 2012 to 2015. For the Bushy Park Reservoir modeling framework to be used directly for taste-and-odor issues, cyanobacteria must be simulated and calibrated based on observations of cyanobacteria biomass concentrations. In addition to the cyanobacteria sampling conducted within the reservoir between 2012 and 2015, the new model calibration would also require new algae biomass data-collection efforts to characterize the external sources of cyanobacteria entering the Bushy Park Reservoir from tributaries, as well as the internal cycling, production, and decay of cyanobacteria in the hydrologic system. Further improvements to the EFDC model would include expanding the collection of boundary condition datasets, such as water-quality monitoring to determine improved nutrient loads into the model domain. Along with improved water-quality monitoring for the major boundary conditions, continuous discharge, for both Foster Creek and the Back River, would further constrain the flow balance and the loads into Bushy Park Reservoir. In addition to better boundary-condition characterization, it is important to better characterize possible shortcomings specifically to the model domain, such as the grid resolution, bathymetry, and numerical hydrodynamic errors. Further consideration of the model may involve a sensitivity analysis to determine if errors in the simulation outputs, such as discharge, water-surface elevations, and salinity, were more likely caused by poor boundary condition characterization or, specifically, the model setup. Three model scenarios were run with the revised Bushy Park Reservoir model: (1) reduced withdrawals from one of the large intake-discharge locations for Bushy Park Reservoir, the Williams Station; (2) elevated (above background levels) ocean water level causing saltwater intrusion from the ocean through Durham Canal into Bushy Park Reservoir; and (3) overtopping of the Back River Dam at the southernmost end of Bushy Park Reservoir. For the reduced withdrawals scenarios, the largest shift in flow resulted near the Williams Station intake, with the next largest flow change at the southern end of Bushy Park Reservoir, and a net increase in flow out of the Bushy Park Reservoir to the Cooper River by way of the Durham Canal. The effect resulting from scenario 3 on water quality and salinity was small, with larger increases for dissolved oxygen than other constituents at several monitoring stations. For the two scenarios related to saltwater intrusion (including dam overtopping), the changes in salinity generally were found to dissipate in the following 2 weeks and generally back to baseline salinity conditions within 3 months. This result did vary depending on the severity of the storm or length of the dam overtopping event.

South Carolina↗

Least Bell's Vireos and Southwestern Willow Flycatchers at the San Luis Rey Flood Risk Management Project Area in San Diego County, California—Breeding activities and habitat use—2022 annual report

Executive Summary We completed four protocol surveys for Least Bell’s Vireos ( Vireo bellii pusillus ; vireo) during the breeding season, supplemented by weekly territory monitoring visits. We identified a total of 133 territorial male vireos; 114 were confirmed as paired, and 3 were confirmed as single males. For the remaining 16 territories, we were unable to confirm breeding status. Two transient vireos were detected in 2022. The vireo population in the Project Area increased by 9 percent from 2021 to 2022. The vireo population at Marine Corps Base Camp Pendleton also increased (4 percent), whereas the population at Marine Corps Air Station remained relatively stable (decreased from 10 pairs to 9) and the Otay River population decreased by 10 percent (2 territories). We used an index of treatment (Treatment Index) to evaluate the effect of on-going vegetation clearing on the Project Area vireo population. The Treatment Index measures the cumulative effect of vegetation treatment within a territory (since 2005) by using the percentage area treated weighted by the number of years since treatment. We determined that the Treatment Index for unoccupied habitat was more than four times that of occupied habitat, indicating that vireos selected habitat that was less treated in which to settle. We monitored vireo nests at three general site types: (1) within the flood channel where exotic and native vegetation removal has occurred regularly (Channel), (2) three sites near the flood channel where limited exotic and native vegetation removal has occurred (Off-channel), and (3) three sites that have been actively restored by planting native vegetation (Restoration). Nesting activity was monitored in 80 territories, 3 of which were occupied by single males and 1 by a male whose breeding status could not be confirmed. Overall, 38 percent of completed nests were successful and nest success did not differ among the three sites. In 2022, there were no differences with regard to clutch size, hatching, or fledging success among Channel, Off-channel and Restoration sites. Overall breeding success and productivity were slightly higher in 2022 than in 2021, with 72 percent of pairs fledgling at least one young and pairs fledging an average of 2.2±1.7 young. To investigate if the cumulative years of treatment had an effect on vireo reproductive effort, we looked at the effects of the Treatment Index on reproductive parameters. Results from generalized linear models indicated that treatment did not have an effect on vireo nesting effort or the number of vireo fledglings per pair produced in 2022. Similarly, we did not detect an effect of Treatment Index on daily survival rate (DSR) of nests. Analysis of vegetation data collected at vireo nests from 2006 to 2022 did not reveal an effect of vegetation cover at the nest on DSR. We did find, however, that Channel nests were placed higher in the host plant than Off-channel nests. In the Channel and Off-channel sites, successful nests were placed closer to the edge of the host plants than unsuccessful nests. Additionally, successful Off-channel nests were placed lower in the vegetation, in shorter host plants, and closer to the edge of the vegetation clump than unsuccessful nests. Red/arroyo willow ( Salix laevigata or Salix lasiolepis ) were the species most commonly selected for nesting by vireos in all three site types. Black willow ( Salix gooddingii ) and mule fat ( Baccharis salicifolia ) also were commonly used. Vireos used a wider variety of species for nesting in Channel and Off-channel sites (eight and six species, respectively) compared to Restoration sites (two species), although there was limited nesting in Restoration sites in 2022. There were 43 vireos banded before the 2022 breeding season that were resighted and identified at the Project Area in 2022, all of which were originally banded in the Project Area. Adult birds of known age ranged from 1 to 7 years old. A total of 146 vireos were newly banded in 2022. There were 8 adult vireos banded with a unique color combination, and 138 nestlings were banded with a single dark blue numbered federal band on the left leg. Between 2006 and 2022, survivorship of males (66±11 percent) was consistently higher than that of females (59±12 percent). First-year birds from 2006 to 2022 had an average annual survivorship of 15±6 percent. First-year dispersal in 2022 averaged 6.7±7.4 kilometers (km), with the longest dispersal (15.3 km) by a male that was recaptured at Fallbrook Creek, Fallbrook Naval Weapons Station (FNWS). From 2007 to 2011, most returning first-year vireos returned to the Project Area, whereas from 2014 to 2016, the majority of returning birds dispersed to areas outside of the Project Area. From 2018 to 2021, the trend shifted, and more first-year vireos returned to the Project area. In 2022, only one first-year vireo returned to the project area and two dispersed to sites outside the Project Area (upstream to the middle San Luis Rey River and to Fallbrook Creek, FNWS). However, the total number of identified first-year vireos was low and the trend in 2022 will likely shift as additional returning first-year vireos are identified in subsequent years. Most of the returning adult male vireos showed strong between-year site fidelity to their previous territories. Seventy-three percent of males (27/37) occupied a territory in 2022 that they had defended in 2021 (within 100 meters [m]). There were no females (0/4) detected in 2022 that returned to a territory they occupied in 2021; however, 50 percent of females (2/4) detected in 2022 returned to areas adjacent to their previous territories (within 300 m). The average between-year movement for returning adult vireos was 0.3±0.7 km. The amount of treatment at adults’ 2021 territories did not affect the distance adults moved to their 2022 territories. We completed four protocol surveys for the endangered Southwestern Willow Flycatcher ( Empidonax traillii extimus ; flycatcher) at the Project Area between May 16 and July 25, 2022. Four transient Willow Flycatchers were detected in the Project Area in 2022. Two transients were detected in Reach 1, one in Reach 3a, and one in Pilgrim Pond. There were not any resident flycatchers documented in the Project Area in 2022. A total of 46 vegetation transects (528 points) were sampled at the Project Area in 2022. Seventy-one percent (378/528) of points were located in the Channel, and 22 percent (115/528) were in Upper Pond. The remaining 7 percent (35/528) of points were at the Whelan Restoration site. Foliage cover below 2 m was higher at the Channel points compared to Upper Pond and Whelan Restoration, which can be attributed to the dense herbaceous vegetation that grows after mowing. Above 2 m, foliage cover was similar at the Channel and Whelan Restoration sites and was higher than at Upper Pond. Average canopy height was higher in the Channel (5.6±3.4 m) compared to Upper Pond (4.7±2.9 m) and Whelan Restoration (4.6±1.9 m). From 2006 to 2022, total foliage cover declined above 2 m in the Channel, in contrast to Upper Pond and Whelan Restoration, where little directional change in vegetation cover has occurred and where vegetation cover has largely recovered to 2006 levels. Within the Channel, the steepest declines occurred between 2009 and 2013 and between 2014 and 2016. Since 2016, we observed an increase in foliage cover, largely herbaceous, between 0 and 2 m within the Channel. The percent cover remained below levels detected before 2009 for other height classes. We sampled vegetation at 44 vireo nests and 44 random plots (“territory” plots) within territories in the Channel and Upper Pond after the 2022 breeding season. Vireos in the Channel established territories in areas with significantly more cover from 3 to 6 m but less cover below 2 m relative to the available habitat. Within territories, Channel vireos selected nest sites with significantly more foliage cover from 2 to 3 m. Vireos at Upper Pond established territories in areas with significantly more foliage cover from 5 to 6 m and below 1 m relative to available habitat. However, within territories, Upper Pond vireos selected nest sites with significantly less foliage cover from 5 to 6 m and below 1 m. Data either are not available or have limited availability owing to restrictions of the funding entity (U.S. Army Corps of Engineers). Please contact Christopher Chabot, Planning Division, Los Angeles District, U.S. Army Corps of Engineers, for more information.

Califoria↗

Reservoir evolution, downstream sediment transport, downstream channel change, and synthesis of geomorphic responses of Fall Creek and Middle Fork Willamette River to water years 2012–18 streambed drawdowns at Fall Creek Lake, Oregon

Executive Summary Chapter A. Introduction Fall Creek Dam impounds Fall Creek Lake, a 10-kilometer-long reservoir in western Oregon and is operated by the U.S. Army Corps of Engineers (USACE) primarily for flood-risk management (or flood control) in late autumn through early spring months, as well as for water quality, irrigation, recreation, and habitat in late spring through early autumn. Since 2011 (water year [WY] 2012), Fall Creek Lake has been temporarily drawn down each year to facilitate downstream passage of juvenile spring Chinook salmon ( Oncorhynchus tshawytscha ) through the 55-meter (m) high dam. This annual dam operation is temporary, typically lasting about 1–2 weeks from WY 2012 through 2020 (drawdown operations in WY 2022–24 have increased to more than 6 weeks). Drawdown of the reservoir results in lake levels being lowered to the elevation near the historical, pre-dam streambed. The annual streambed drawdowns of WY 2012–18 have improved fish passage and led the USACE to formally adopt streambed drawdowns as part of annual operations at Fall Creek Dam. However, temporarily lowering the lake to streambed creates free-flowing conditions in the reservoir that result in the erosion and episodic export of predominantly sand and finer-grained sediments (less than 2 millimeters [mm]) to the lower gravel-bed reaches of Fall Creek and the Middle Fork Willamette River. The introduction of large volumes of sand and finer-grain sediment into the dam-regulated reaches downstream from Fall Creek Dam prompted questions about the geomorphic responses to annual streambed drawdowns within Fall Creek Lake and downstream reaches along Fall Creek and the Middle Fork Willamette River. The U.S. Geological Survey (USGS) in partnership with USACE initiated a comprehensive geomorphic and sediment transport investigation to assess the coupled processes of reservoir erosion, sediment evacuation from Fall Creek Lake, and patterns of sediment transport and deposition in reaches downstream from the Fall Creek Dam that have resulted from annual streambed drawdowns. The purpose of this report is to systematically describe the processes of sediment erosion, transport, and deposition at Fall Creek Lake and geomorphic interactions between reaches upstream and downstream from Fall Creek Dam that relate to dam operations. Specifically, this report focuses on evaluating geomorphic responses to streambed drawdowns from WY 2012 through 2018 and placing drawdown-induced geomorphic responses within the broader context of physiographic and historical conditions and dam operations of Fall Creek and Middle Fork Willamette Rivers. Key objectives for this study were to characterize changes in reservoir morphology and substrate at Fall Creek Lake, describe the character and temporal pattern of sediment transport downstream from Fall Creek Dam, characterize geomorphic changes in channel reaches downstream from the Fall Creek Dam, and relate these data to the annual streambed drawdowns of WY 2012–18. This study uses multiple independent monitoring and measurement approaches to assess site, reach, and river-scale geomorphic responses to drawdowns to inform dam and reservoir management. Patterns and processes of reservoir evolution were assessed with geomorphic mapping and volumetric analyses of topography through comparison of multiple digital surface models (DSMs). Just downstream from Fall Creek Dam, analyses of sediment export from the reservoir focused on suspended sediment but also incorporated bedload analyses to assess sediment sizes. Geomorphic assessments downstream from the dam used reach-scale and site-scale approaches to document changes in channel morphology and substrate, including site measurements of sand and finer-grained sediment deposition and in-channel bed-material, volumetric change analyses from comparison of digital elevation models (DEMs), and repeat geomorphic mapping. Findings from this study inform river management and dam operations by providing an understanding of (1) coupled upstream-downstream geomorphic responses to the Fall Creek Lake streambed drawdowns, (2) geomorphic responses of Fall Creek Lake streambed drawdowns in comparison to drawdowns at other large dams, (3) controls on reservoir erosion and downstream geomorphic responses, and (4) implications for future hydrogeomorphic changes that may result from continued drawdowns and monitoring activities to assess those changes. Chapter B. Reservoir Morphology and Evolution Related to Dam Operations at Fall Creek Lake To understand the volume and distribution of sediment accumulation in Fall Creek Lake since dam closure in 1965, decadal-scale sedimentation patterns (spanning approximately 1965–2016) are evaluated using a combination of storage curve analyses and geomorphic mapping. Short-term (drawdown event-scale) patterns of erosion, sedimentation, and sediment export downstream are evaluated using a combination of geomorphic mapping and change detection analyses that quantify the distribution and total volume of sediment erosion and deposition within Fall Creek Lake. Geomorphic mapping of reservoir topography and analyses of historical datasets reveals four categories of landforms and sediment processes within Fall Creek Lake related to lake level operations: lacustrine sedimentation expressed in the reservoir floor, fluvial erosion and deposition within historical stream channels during streambed drawdowns, channel-like features created by erosion within the reservoir floor during streambed drawdowns, and erosion on reservoir hillslopes. Where the reservoir floor is mapped for this study as pelagic (deep water), deposition up to 3 meters (m) thick by lacustrine processes and burial of pre-dam topography with deposits thinning toward the edges of the valley floor and upstream areas of reservoir are observed. Despite over 50 years of sediment accumulation since dam construction, the main stream channels of Fall and Winberry Creeks (or reservoir thalwegs) through the reservoir are well defined, though their distinct morphology is likely influenced by a long history of recurring historical drawdowns to or near streambed since dam construction. Unregulated streamflow and sediment transport through the reservoir primarily are confined to these channels during the streambed drawdown periods. Erosional channel-like features created by drawdowns are carved through underlying, unconsolidated reservoir floor sediments and are most prominent in the lower reservoir below minimum conservation pool (the low pool elevation during winter flood season); sediment generated from the formation of these drawdown channels is more likely to be transported through and out of the reservoir than sediment deposits along the reservoir hillslopes at the valley margins that are separated from main channels by areas of low-gradient reservoir floor. Morphologic changes in the lower reservoir topography between January 2012 and November 2016 indicate overall net erosion of about 129,500 cubic meters (m 3 ). The most prominent geomorphic changes occurred along the main channels of Fall and Winberry Creeks near the Fall Creek Dam where incision, lateral migration, and slumping banks resulted in vertical and lateral adjustments to channel position, whereas most changes fell below the detectable limit on higher-elevation reservoir floor surfaces except where erosion occurred along features mapped as drawdown channels. Chapter C. Sediment Delivery from Fall Creek Lake and Transport through Downstream Reaches USGS implemented a sediment monitoring program in WY 2013–18 to evaluate the quantity and character of reservoir sediment exported from Fall Creek Lake during streambed drawdowns. Turbidity and suspended sediments were monitored annually autumn through spring to span the WY 2013–18 streambed drawdowns; however, unequal monitoring timeframes each year reduced the ability to compare results and factors affecting sediment export from the reservoir difficult between years. These data were originally measured to develop regressions and compute suspended-sediment loads (SSL). Bedload sediment monitoring from a cableway at the Fall Creek streamgage was completed in the autumn-winter of WY 2013 and 2017. The limited number of samples and presumed variability in sediment supply from the reservoir precluded construction of streamflow and bedload discharge relations to compute more than instantaneous bedload. Sand and finer-grained silts and clays were transported from the reservoir in suspension, though some coarser grains (up to 32 mm) were also mobilized and transported downstream from the dam as bedload. Observations of increased sediment transport downstream from Fall Creek Dam coincided with lake levels approaching about 3 m (10 feet [ft] or elevation 690 ft) above the streambed regulating outlets. Suspended-sediment loads computed for the full monitoring periods WY 2013–18 at the Fall Creek streamgage, located 1.4 kilometers (km) downstream from Fall Creek Dam, range from 54,700 metric tons (t) in WY 2013 to 13,900 t in WY 2018. Although the total annual SSL varied from year to year, the overall seasonal patterns of suspended sediment transport throughout each year were similar during monitoring in WY 2013-18. Suspended-sediment loads were low prior to the drawdown, then increased rapidly as lake levels lowered and approached the streambed. In the weeks following the drawdown period, as pool levels were increased, SSL remained slightly elevated above pre-drawdown levels but generally declined through the following winter and spring except during streamflow-driven pulses of suspended-sediment transport. WY 2013 had the greatest total computed SSL for each streambed drawdown and partial-year monitoring period. SSL computed for the partial-year period have generally decreased since WY 2013 and have varied by about 6,800 t with the exception of WY 2014. WY 2014 SSL reflects anomalously low sediment export due to low streamflows and freezing conditions that stabilized reservoir floor deposits. Bedload measurements in the short 1.4-km reach between Fall Creek Dam and the Fall Creek streamgage showed an inverse correlation between bedload transport rates and discharge, which probably reflects diminishing supply of coarse-sized sediment. Sand was more abundant (60–100 percent) than gravel in bedload samples confirming sand and finer-grained sediment dominated sediment evacuated from the reservoir during streambed drawdowns at Fall Creek Lake. Chapter D. Geomorphic Responses to Fall Creek Lake Streambed Drawdowns Downstream from Fall Creek Dam In the days, weeks, and months following streambed drawdown operations at Fall Creek Dam through WY 2018, sites downstream from the dam displayed a variety of geomorphic responses to reservoir sediment delivery within the main channel and overbank areas. Evaluation of streambed elevations at two streamgages located 1.4 km downstream from the dam on Fall Creek and 16.3 km downstream from the dam on the Middle Fork Willamette River indicated the effects of drawdown sediment on bed elevations were modest and transient. Repeat particle size measurements (October 2015 and September 2016) at five sites along Fall Creek and the Middle Fork Willamette River showed similar grain-sized distributions that do not reveal substantial deposition of fine-grained sediment related to the WY 2016 streambed drawdown. Altogether, these findings indicate that transport capacity in the main, low-flow channels of Fall Creek and Middle Fork Willamette River during WY 2012–18 was sufficient to mobilize and evacuate reservoir sediments from streambed drawdowns or other bank material and tributary sources. However, other monitoring for this study indicate low-velocity zones in off-channel areas are prime locations for sand and finer-grain sediment deposition. Patterns of overbank sediment accumulation indicate that the magnitude and timing of overbank deposition on bars and low-elevation floodplain varies with proximity to the dam, geomorphic setting, streamflows, and other factors. Sand and finer-grained reservoir sediments carried as suspended-sediment load in the reaches downstream from Fall Creek Dam were deposited in overbank areas as observed with clay-horizon markers during WY 2016–17. Overbank deposition quantified with Geomorphic Change Detection (GCD) software evaluated landform-scale patterns of erosion and deposition using repeat light detection and ranging (lidar) surveys at two sites in the Upper Fall Creek reach and one site in the Jasper reach for 3 years (2012–15) and one site in the Clearwater reach for 6 years (2009–15). Deposition thickness and spatial patterns from the GCD analysis were variable; some sites had dispersed but measurable deposition while at others, deposition was highly localized and exceeded 1 m in depth. Patterns of overbank deposition illustrate interactions among bar morphology, local hydraulics, and suspended-sediment transport dynamics that can create patches of highly localized deposition. The measured deposition at the two Fall Creek GCD sites likely resulted from reservoir sediments released from Fall Creek Lake during streambed drawdowns in WY 2016 and 2017 because the limited sediment inputs from bank material (geomorphically laterally stable reach) or tributaries (no significant tributaries) provided few other sediment sources. On the Middle Fork Willamette River, observed patterns of overbank deposition could reflect sediment sourced from upstream tributaries, bank erosion, or Fall Creek Lake streambed drawdown operations. Despite the introduction of several thousand tons of reservoir sediment delivered from the Fall Creek Lake streambed drawdowns to below-dam river corridors, reach-scale mapping of channel features downstream from Fall Creek Dam shows minimal evidence of changes in channel planform or landforms that can be attributed to a drawdowns in WY 2012–16. On Upper Fall Creek reach, widespread increases in gravel bars or other in-channel sediment did not result from the five streambed drawdowns. The main changes attributable to sediment releases from Fall Creek Lake were localized increases in vegetated bar area, particularly on channel margin areas where sand and finer-grain sediment was deposited and rapidly colonized by vegetation. The area of mapped secondary water features decreased between 2005 and 2016, but that may be due to lower discharges depicted in the 2016 aerial photographs and less mapped area of inundation. Primary changes along the Lower Fall Creek reach include a 6.4 percent decrease in area of secondary water features between 2011 and 2016, and a nearly twofold increase in the area of unvegetated bars. Immediately downstream from the Fall Creek confluence, there were negligible changes in the location and areas of vegetated bars and the main wetted channel between 2005 and 2016, and local increases in bar area cannot be attributed solely to deposition of reservoir sediments from Fall Creek Lake because (1) areas along the Middle Fork Willamette River just upstream from the Fall Creek confluence display similar type and magnitude of changes and (2) some of the increases at the confluence area pre-date the drawdowns. The cumulative effect of sediment releases from Fall Creek Lake streambed drawdowns from WY 2012 to 2016 on downstream channel planform and landforms are modest compared to the river-scale transformations and planform changes that occurred in the decades following dam construction. Chapter E. Discussion of Geomorphic Responses of Fall Creek and Middle Fork Willamette River to Streambed Drawdowns at Fall Creek Lake Multiple aspects of Fall Creek Dam infrastructure and operations exert first-order controls on the magnitudes of reservoir erosion that occur during the streambed drawdowns and ultimately determine the sediment delivery to downstream reaches. Key aspects of the dam and its operations that are most relevant to assessing geomorphic responses to streambed drawdowns include the (1) dam infrastructure, including configuration and size of regulating outlets and their proximity to the streambed which dictates the capacity and competence of the river to deliver sediment to downstream reaches and mode of sediment transport as suspended-sediment load or bedload; (2) frequency of historical drawdowns and long-term, year-round dam operations and lake level management, which partly dictate reservoir morphology and locations and magnitudes of readily erodible materials; (3) dam operations and hydroclimatic conditions during the streambed drawdown (including length of the drawdown and streamflows entering the reservoir), which directly control the timing, duration and magnitude of reservoir erosion and sediment evacuation; and (4) dam operations following the streambed drawdown operation that regulate streamflows (and thereby sediment transport conditions) downstream of Fall Creek Dam which primarily reflect interactions between hydroclimatic conditions and flood control operations. Patterns of sediment erosion and evacuation observed in this study at Fall Creek Lake from WY 2012–18 suggest that reservoir erosion during annual streambed drawdowns may remain similar or decrease in future years assuming (1) annual streambed drawdown operations are implemented in similar manner as the WY 2012–18 drawdowns (in terms of duration, late autumn or early winter implementation, rate of pool-level lowering to reach streambed, and other factors), (2) streambed drawdowns coincide with similar conditions as were observed WY 2012–18 (similar sediment yield into reservoir, low reservoir inflows, limited precipitation, moderate air temperature), and (3) no major geomorphic changes in the main reservoir channels of Fall and Winberry Creeks occur (for example, channel avulsion). Under such conditions, it is hypothesized that the stream channel within the reservoir would achieve a quasi-equilibrium state with respect to annual influx and export of sediment and aided by the substantial amount of in-channel bedrock, will remain laterally stable without erosion across reservoir deposits. Patterns of sediment transport measured at the Fall Creek streamgage downstream from Fall Creek Dam provide insight into the potential effects of future streambed drawdowns at Fall Creek Lake. Analyses of suspended sediment measured in WY 2013–18 show a major reduction in suspended-sediment loads between WY 2013 and later years, indicating streamflows transporting sediment through the reservoir to downstream reaches during streambed drawdowns have become supply limited. The 6-year suspended-sediment monitoring and sampling program is insufficient to make predictions about future sediment transport conditions because of uneven monitoring periods and varying controls on reservoir sediment erosion. It is likely that future suspended-sediment loads will be variable but similar to those observed in WY 2015–18 if operational, climatic, and geomorphological factors remain similar to those monitored WY 2015–18. Suspended-sediment loads downstream from Fall Creek Lake will likely remain highest when regulating outlets are fully open and Fall Creek is free flowing with the reservoir fully drained with little to no residual pool. Over time, it is possible that the suspended-sediment loads would reflect mobilization of reservoir sediment deposited in the previous year rather than erosion of sediment deposited years or decades earlier. Bedload is likely to remain a small fraction of the total sediment load evacuated from the reservoir and is relatively modest compared with pre-dam bedload transport rates because most coarse sediment remains trapped by the dam. If sediment releases from Fall Creek Lake and ensuing streamflow conditions follow a similar pattern in the future as was assessed in this study spanning WY 2012–18, near-term geomorphic adjustments downstream of the dam are expected to be modest. Barring major operational, climatic, and geomorphological changes, local site-scale deposition on bars, overbank areas, or off-channel features that persists several months after the streambed drawdown will likely continue to be highly variable, ranging from negligible to several centimeters of deposition. At the landform-scale, low velocity areas nearest to Fall Creek Dam will likely continue to undergo rapid deposition immediately during and after a streambed drawdown event, similar to patterns observed for WY 2012–18. Some of the sediment entering these off-channel features and margin areas may be temporarily stored, then later remobilized and dispersed farther downstream. But if newly deposited sediment persists through the following spring, there is a greater likelihood that local vegetation will establish, reinforce deposited material, and trap sediment during later drawdowns. The reach-scale geomorphic changes may become more apparent if (1) streambed drawdowns continued for several decades, and geomorphic changes were measured at decadal scales or (2) the amount of sediment introduced to downstream reaches substantially increased and (or) sediment transport capacity decreased. The continued streamflow regulation of Fall Creek Dam after sediment releases provides an opportunity to strategically manage streamflows during and after the streambed drawdowns to minimize downstream sediment impacts and ensure other operational thresholds are satisfied. This study provides a comprehensive foundation of datasets and geomorphic analyses to inform dam operations at Fall Creek Lake, monitor sediment transport downstream, and consider operational schemes for future drawdowns. The datasets from this study also provide baselines of sediment transport and geomorphic conditions to assess future changes in reservoir and downstream environments. Future monitoring could be tailored to address specific questions regarding the long-term geomorphic effects of streambed drawdowns on fluvial habitats, flood hazards, cultural resources, or downstream water quality. Future monitoring activities could focus on the relevant geomorphic processes and spatial domains within the three categories used for this study: (1) reservoir erosion and net sediment evacuation, (2) sediment delivery to downstream reaches, including magnitude and temporal pattern of sediment transport, and (3) geomorphic responses of downstream reaches to sediment delivery. Specifically, high priority future monitoring activities could include: Repeat topographic or photographic surveys in the reservoir to characterize changes occurring within individual drawdowns, to quantify sediment export, to determine temporal changes in reservoir storage, and to identify locations of erosion and deposition. Continuous, year-round turbidity monitoring supplemented with suspended-sediment measurements at a streamflow-gaging station immediately downstream from the dam to quantify sediment export. Repeat geomorphic monitoring, mapping, or modeling in downstream reaches to track changes in channel and over bank features using a combination of site- and reach-scale monitoring approaches. This could support assessments of sediment deposition and ensuing vegetation encroachment on flood hazards and habitats and examine how sediment transport and depositional processes may be affected by different sediment supply, streamflow, or dam management scenarios.

Oregon↗