SEARCH · Geology Reports

Results for “Biology Direct”

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

SourceResearch connectionFrom yearThrough year

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

At least 865 records · Page 48Linked to original sources

USGS1999

Time course of salinity adaptation a strongly euryhaline estuarine teleost, Fundulus Heteroclitus: A multivariable approach

Freshwater-adapted killifish ( Fundulus heteroclitus ) were transferred directly from soft fresh water to full-strength sea water for periods of 1 h, 3 h, 8 h and 1, 2, 7, 14 and 30 days. Controls were transferred to fresh water for 24 h. Measured variables included: blood [Na + ], osmolality, glucose and cortisol levels, basal and stimulated rates of ion transport and permeability of in vitro opercular epithelium, gill Na + /K + -ATPase and citrate synthase activity and chloride cell ultrastructure. These data were compared with previously published killifish cystic fibrosis transmembrane conductance regulator (kfCFTR) expression in the gills measured over a similar time course. Plasma cortisol levels peaked at 1 h, coincident with a rise in plasma [Na + ]. At 8 h after transfer to sea water, a time at which previous work has shown kfCFTR expression to be elevated, blood osmolality and [Na + ] were high, and cortisol levels and opercular membrane short-circuit current ( I sc ; a measure of Cl − secretion rate) were low. The 24 h group, which showed the highest level of kfCFTR expression, had the highest plasma [Na + ] and osmolality, elevated plasma cortisol levels, significantly lower opercular membrane resistance, an increased opercular membrane ion secretion rate and collapsed tubule inclusions in mitochondria-rich cells, but no change in gill Na + /K + -ATPase and citrate synthase activity or plasma glucose levels. Apparently, killifish have a rapid (<1 h) cortisol response to salinity coupled to subsequent (8–48 h) expression of kfCFTR anion channel proteins in existing mitochondria-rich cells that convert transport from ion uptake to ion secretion.

Journal of Experimental Biology↗

USGS2025

Integrating the resist–accept–direct framework into natural resource decision-making processes for climate adaptation

The resist–accept–direct (RAD) framework for climate adaptation is a useful tool, particularly when conservation practitioners know they need to address climate change but do not know where to start or when they struggle to implement conservation actions that are outside the status quo. Some conservation practitioners may view RAD as a decision process that will lead them through selecting climate adaptation actions to meet their objectives; however, RAD may be better suited for use with existing decision processes. RAD can improve adaptation planning processes by helping conservation practitioners examine a broader portfolio of climate adaptation actions. Choosing the actions that meet organizational objectives and long-term goals relies on a sound decision process through which to select RAD-generated actions to implement across space and time. The good news for conservation practitioners is that RAD can easily be integrated into the decision processes they are already using—that is, a new decision process for RAD is not necessary. We examined 6 commonly used decision frameworks in the context of RAD for climate adaptation: structured decision-making, adaptive management, conservation standards, climate-smart conservation, strategic habitat conservation, and scenario planning. Each decision framework can easily incorporate RAD methodologies to create climate adaptation menus and portfolios within existing structures, allowing conservation practitioners to continue to utilize existing frameworks for their valuable decision-making tools. Integrating the RAD framework can enhance them by explicitly considering future uncertainties and providing a menu of climate adaptation actions for navigating changing ecosystems.

Conservation Biology↗

USGS2020

Effects of increased temperature on arctic slimy sculpin Cottus cognatus is mediated by food availability: Implications for climate change

Lakes are vulnerable to climate change, and warming rates in the Arctic are faster than anywhere on Earth. Fishes are sensitive to changing temperatures, which directly control physiological processes. Food availability should partly dictate responses to climate change because energetic demands change with temperature, but few studies have simultaneously examined temperature and food availability. We used a fully factorial experiment to test effects of food availability and temperature (7.6, 12.7, and 17.4°C; 50 days) on growth, consumption, respiration, and excretion, and effects of temperature (12 and 19.3°C; 27 days) on habitat use and growth of a common, but understudied, mid-level consumer, slimy sculpin Cottus cognatus , in arctic lakes. We also used bioenergetics modelling to predict consumptive demand under future warming scenarios. Growth rates were 3.4× higher at 12.7°C in high food compared to low food treatments, but the magnitude of differences depended on temperature. Within low food treatments, there was no statistical difference in growth rates among temperatures, suggesting food limitation. Consumption, respiration, and nitrogen excretion increased with temperature independent of food availability. Lower growth rates coincided with lower phosphorus excretion at the highest temperature, suggesting that fish selectively retained phosphorus at high temperatures and low food. In habitat choice experiments, fish were more likely to use the 12°C side of the tank, closely matching their optimal temperature. We predicted a 9% increase in consumption is required to maintain observed growth under a 4°C warming scenario. These results highlight considering changes in food resources and other associated indirect effects (e.g. excretion) that accompany changing temperatures with climate change. Depending on how food webs respond to warming, fish may cope with predicted warming if density-dependent feedback maintains population sizes.

Alaska↗

USGS2004

Coping with mist-net capture-rate bias: Canopy height and several extrinsic factors

Many factors other than a species' actual abundance can affect mist-net capture rates. We used ANCOVA models to quantify some potential biases and control their effects, producing adjusted estimates of capture rates that are more directly comparable among mist-net stations. Data came from 46 two-day mist-net sessions from September 1990 to May 1992 at six subtropical forest stations in the Rio Bravo Conservation and Management Area, northwest Belize. Factors evaluated included canopy height at net sites, long-term net shyness (days elapsed between first and last netting day of the entire study period), season (wet vs. dry), total rainfall during a netting session, and temperature. Number of individuals and species captured/10 net-h declined at each net with increasing canopy height above the net. Capture rates differed significantly among some of the stations. Elapsed days and rainfall caused significant bias in capture rates, which were statistically controlled within the ANCOVA, whereas season and temperature did not. Capture rates varied among sessions, but there was a slight and significant decline over the entire study period for all stations combined. Rainfall significantly depressed capture rates somewhat on a daily basis, but capture rates did not differ between wet and dry seasons. When we replaced the station variable in the ANCOVA with mean canopy height, the model was still highly significant, but did not explain as much of the variation in capture rates. Statistical analysis provides an objective means of interpreting data and estimating reliability, but only if statistical assumptions of the analyses are met. We discuss the need for including randomization in the experimental design, standardizing netting protocol, and quantifying sources of bias in the field, before ANCOVA or other parametric statistical techniques can be used to partition effects of biases.

Studies in Avian Biology↗

USGS2010

Multi-species occurrence models to evaluate the effects of conservation and management actions

Conservation and management actions often have direct and indirect effects on a wide range of species. As such, it is important to evaluate the impacts that such actions may have on both target and non-target species within a region. Understanding how species richness and composition differ as a result of management treatments can help determine potential ecological consequences. Yet it is difficult to estimate richness because traditional sampling approaches detect species at variable rates and some species are never observed. We present a framework for assessing management actions on biodiversity using a multi-species hierarchical model that estimates individual species occurrences, while accounting for imperfect detection of species. Our model incorporates species-specific responses to management treatments and local vegetation characteristics and a hierarchical component that links species at a community-level. This allows for comprehensive inferences on the whole community or on assemblages of interest. Compared to traditional species models, occurrence estimates are improved for all species, even for those that are rarely observed, resulting in more precise estimates of species richness (including species that were unobserved during sampling). We demonstrate the utility of this approach for conservation through an analysis comparing bird communities in two geographically similar study areas: one in which white-tailed deer (Odocoileus virginianus) densities have been regulated through hunting and one in which deer densities have gone unregulated. Although our results indicate that species and assemblage richness were similar in the two study areas, point-level richness was significantly influenced by local vegetation characteristics, a result that would have been underestimated had we not accounted for variability in species detection.

Biological Conservation↗

USGS2008

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↗

USGS2016

Fire effects on wildlife in Central Hardwoods and Appalachian regions

Fire is being prescribed and used increasingly to promote ecosystem restoration (e.g., oak woodlands and savannas) and to manage wildlife habitat in the Central Hardwoods and Appalachian regions, USA. However, questions persist as to how fire affects hardwood forest communities and associated wildlife, and how fire should be used to achieve management goals. We provide an up-to-date review of fire effects on various wildlife species and their habitat in the Central Hardwoods and Appalachians. Documented direct effects (i.e., mortality) on wildlife are rare. Indirect effects (i.e., changes in habitat quality) are influenced greatly by light availability, fire frequency, and fire intensity. Unless fire intensity is great enough to kill a portion of the overstory, burning in closed-canopy forests has provided little benefit for most wildlife species in the region because it doesn’t result in enough sunlight penetration to elicit understory response. Canopy reduction through silvicultural treatment has enabled managers to use fire more effectively. Fire intensity must be kept low in hardwoods to limit damage to many species of overstory trees. However, wounding or killing trees with fire benefits many wildlife species by allowing increased sunlight to stimulate understory response, snag and subsequent cavity creation, and additions of large coarse woody debris. In general, a fire-return interval of 2 yr to 7 yr benefits a wide variety of wildlife species by providing a diverse structure in the understory; increasing browse, forage, and soft mast; and creating snags and cavities. Historically, dormant-season fire was most prevalent in these regions, and it still is when most prescribed fire is implemented in hardwood systems as burn-days are relatively few in the growing season of May through August because of shading from leaf cover and high fuel moisture. Late growing-season burning increases the window for burning, and better control on woody composition is possible. Early growing-season fire may pose increased risk for some species, especially herpetofauna recently emerged from winter hibernacula (April) or forest songbirds that nest in the understory (May to June). However, negative population-level effects are unlikely unless the burned area is relatively large and early growing-season fire is used continually. We did not find evidence that fire is leading to population declines for any species, including Endangered Species Act (ESA)-listed species (e.g., Indiana bat [ Myotis sodalis Mill. Allen] or northern long-eared bat [ M. septentrionalis Trouess.]). Instead, data indicate that fire can enhance habitat for bats by increasing suitability of foraging and day-roost sites. Similarly, concern over burning and displacement of woodland salamanders (Plethodontidae), another taxa of heightened conservation concern, is alleviated when fire is prescribed along ecologically appropriate aspect and slope gradients and not forced into mesic, high site index environments where salamanders are most common. Because topography across the Central Hardwoods and Appalachians is diverse, we contend that applying fire on positions best suited for burning is an effective approach to increase regional landscape heterogeneity and biological diversity. Herein, we offer prescriptive concepts for burning for various wildlife species and guilds in the Central Hardwoods and Appalachians.

Appalachian region, Central Hardwoods region↗

USGS2002

Potential effects of structural controls and street sweeping on stormwater loads to the lower Charles River, Massachusetts

The water quality of the lower Charles River is periodically impaired by combined sewer overflows (CSOs) and non-CSO stormwater runoff. This study examined the potential non-CSO load reductions of suspended solids, fecal coliform bacteria, total phosphorus, and total lead that could reasonably be achieved by implementation of stormwater best management practices, including both structural controls and systematic street sweeping. Structural controls were grouped by major physical or chemical process; these included infiltration-filtration (physical separation), biofiltration-bioretention (biological mechanisms), or detention-retention (physical settling). For each of these categories, upper and lower quartiles, median, and average removal efficiencies were compiled from three national databases of structural control performance. Removal efficiencies obtained indicated a wide range of performance. Removal was generally greatest for infiltration-filtration controls and suspended solids, and least for biofiltration-bioretention controls and fecal coliform bacteria. Street sweeping has received renewed interest as a water-quality control practice because of reported improvements in sweeper technology and the recognition that opportunities for implementing structural controls are limited in highly urbanized areas. The Stormwater Management Model that was developed by the U.S. Geological Survey for the lower Charles River Watershed was modified to simulate the effects of street sweeping in a single-family land-use basin. Constituent buildup and washoff variable values were calibrated to observed annual and storm-event loads. Once calibrated, the street sweeping model was applied to various permutations of four sweeper efficiencies and six sweeping frequencies that ranged from every day to once every 30 days. Reduction of constituent loads to the lower Charles River by the combined hypothetical practices of structural controls and street sweeping was estimated for a range of removal efficiencies because of their inherent variability and uncertainty. This range of efficiencies, with upper and lower estimates, provides reasonable bounds on the load that could be removed by the practices examined. The upper estimated load reduction from combined street sweeping and structural controls, as a percentage of the total non-CSO load entering the lower Charles River downstream of Watertown Dam, was 44 percent for suspended solids, 34 percent for total lead, 14 percent for total phosphorus, and 17 percent for fecal coliform bacteria. The lower estimated load reduction from combined street sweeping and structural controls from non-CSO sources downstream of Watertown Dam, was 14 percent for suspended solids, 11 percent for total lead, 4.9 percent for total phosphorus, and 7.5 percent for fecal coliform bacteria. Load reductions by these combined management practices can be a small as 1.4 percent for total phosphorus to about 4 percent for the other constituents if the total load above Watertown Dam is added to the load from below the dam. Although the reductions in stormwater loads to the lower Charles River from the control practices examined appear to be minor, these practices would likely provide water-quality benefits to portions of the river during those times that they are most impaired-during and immediately after storms. It should also be recognized that only direct measurements of changes in stormwater loads before and after implementation of control practices can provide definitive evidence of the beneficial effects of these practices on water-quality conditions in the lower Charles River.

Massachusetts↗

USGS2012

Tracking climate impacts on the migratory monarch butterfly

Understanding the impacts of climate on migratory species is complicated by the fact that these species travel through several climates that may be changing in diverse ways throughout their complete migratory cycle. Most studies are not designed to tease out the direct and indirect effects of climate at various stages along the migration route. We assess the impacts of spring and summer climate conditions on breeding monarch butterflies, a species that completes its annual migration cycle over several generations. No single, broad-scale climate metric can explain summer breeding phenology or the substantial year-to-year fluctuations observed in population abundances. As such, we built a Poisson regression model to help explain annual arrival times and abundances in the Midwestern United States. We incorporated the climate conditions experienced both during a spring migration/breeding phase in Texas as well as during subsequent arrival and breeding during the main recruitment period in Ohio. Using data from a state-wide butterfly monitoring network in Ohio, our results suggest that climate acts in conflicting ways during the spring and summer seasons. High spring precipitation in Texas is associated with the largest annual population growth in Ohio and the earliest arrival to the summer breeding ground, as are intermediate spring temperatures in Texas. On the other hand, the timing of monarch arrivals to the summer breeding grounds is not affected by climate conditions within Ohio. Once in Ohio for summer breeding, precipitation has minimal impacts on overall abundances, whereas warmer summer temperatures are generally associated with the highest expected abundances, yet this effect is mitigated by the average seasonal temperature of each location in that the warmest sites receive no benefit of above average summer temperatures. Our results highlight the complex relationship between climate and performance for a migrating species and suggest that attempts to understand how monarchs will be affected by future climate conditions will be challenging.

Global Change Biology↗

USGS2020

Broad-scale impacts of an invasive native predator on a sensitive native prey species within the shifting avian community of the North American Great Basin

Human enterprise has modified ecosystem processes through direct and indirect alteration of native predators’ distribution and abundance. For example, human activities subsidize food, water, and shelter availability to generalist predators whose subsequent increased abundance impacts lower trophic-level prey species. The common raven (Corvus corax; hereafter, raven) is an avian predator, native to the northern hemisphere, that can become invasive when subsidized. Raven populations are increasing at unprecedented rates in many regions globally. Information regarding scale of impact and potential ecological thresholds is needed to guide conservation actions aimed at reducing adverse effects on sensitive prey. We conducted a multi-part analysis to investigate broad-scale variation in raven densities and impacts on nesting greater sage-grouse (Centrocercus urophasianus), an indicator species for sagebrush ecosystems in western North America. We estimated raven densities using 16,000 point surveys over 10 years within the Great Basin, USA, and examined associations with anthropogenic and environmental covariates. Average density was 0.54 ravens km-2 (95% CI: 0.42–0.70), with higher densities at lower relative elevations comprising increased agriculture and development. We then used a reduced dataset to estimate the effect of raven density on sage-grouse nest survival (nests = 737). We identified negative impacts to nesting sage-grouse, especially where raven density exceeded ~ 0.40 km-2, a potential ecological threshold. We mapped regions where elevated raven densities were predicted to depress sage-grouse population growth in the absence of compensatory demographic responses from other sage-grouse life-history stages, and found ~ 64% of sage-grouse breeding areas were adversely impacted by high raven density.

California, Idaho, Nevada, Oregon, Utah↗

USGS2014

2013 update on sea otter studies to assess recovery from the 1989 Exxon Valdez oil spill, Prince William Sound, Alaska

On March 24, 1989, the tanker vessel Exxon Valdez ran aground in Prince William Sound, Alaska, spilling an estimated 42 million liters of Prudhoe Bay crude oil. Oil spread in a southwesterly direction and was deposited on shores and waters in western Prince William Sound (WPWS). The sea otter (Enhydra lutris) was one of more than 20 nearshore species considered to have been injured by the spill. Since 1989, the U.S. Geological Survey has led a research program to evaluate effects of the spill on sea otters and assess progress toward recovery, as defined by demographic and biochemical indicators. Here, we provide an update on the status of sea otter populations in WPWS, presenting findings through 2013. To assess recovery based on demographic indicators, we used aerial surveys to estimate abundance and annual collections of sea otter carcasses to evaluate patterns in ages-at-death. To assess recovery based on biochemical indicators, we quantified transcription rates for a suite of genes selected as potential indicators of oil exposure in sea otters based on laboratory studies of a related species, the mink (Mustela vison). In our most recent assessment of sea otter recovery, which incorporated results from a subset of studies through 2009, we concluded that recovery of sea otters in WPWS was underway. This conclusion was based on increasing abundance throughout WPWS, including increasing numbers at northern Knight Island, an area that was heavily oiled in 1989 and where the local sea otter population had previously shown protracted injury and lack of recovery. However, we did not conclude that the WPWS sea otter population had fully recovered, due to indications of continuing reduced survival and exposure to lingering oil in sea otters at Knight Island, at least through 2009. Based on data available through 2013, we now conclude that the status of sea otters—at all spatial scales within WPWS—is consistent with the designation of recovery from the spill as defined by the Exxon Valdez Oil Spill Trustee Council. The support for this conclusion is based primarily on demographic data, including (1) a return to estimated pre-spill abundance of sea otters at northern Knight Island, and (2) a return to pre-spill mortality patterns. Gene transcription rates in 2012 were similar in sea otters from oiled, moderately oiled and unoiled areas, suggesting abatement of exposure effects in 2012. However, because 2012 gene transcription rates generally were low for sea otters from all areas relative to 2008, we cannot fully interpret these observations without data from a wider panel of genes. This slight uncertainty with respect to the data from the biochemical indicator is outweighed by the strength of the data for the demographic indicators. The return to pre-spill numbers and mortality patterns suggests a gradual dissipation of lingering oil over the past two decades, to the point where continuing exposure is no longer of biological significance to the WPWS sea otter population.

Alaska↗

USGS2019

Effects of 21st century climate, land use, and disturbances on ecosystem carbon balance in California

Terrestrial ecosystems are an important sink for atmospheric carbon dioxide (CO 2 ), sequestering ~30% of annual anthropogenic emissions and slowing the rise of atmospheric CO 2 . However, the future direction and magnitude of the land sink is highly uncertain. We examined how historical and projected changes in climate, land use, and ecosystem disturbances affect the carbon balance of terrestrial ecosystems in California over the period 2001–2100. We modeled 32 unique scenarios, spanning 4 land use and 2 radiative forcing scenarios as simulated by four global climate models. Between 2001 and 2015, carbon storage in California's terrestrial ecosystems declined by −188.4 Tg C, with a mean annual flux ranging from a source of −89.8 Tg C/year to a sink of 60.1 Tg C/year. The large variability in the magnitude of the state's carbon source/sink was primarily attributable to interannual variability in weather and climate, which affected the rate of carbon uptake in vegetation and the rate of ecosystem respiration. Under nearly all future scenarios, carbon storage in terrestrial ecosystems was projected to decline, with an average loss of −9.4% (−432.3 Tg C) by the year 2100 from current stocks. However, uncertainty in the magnitude of carbon loss was high, with individual scenario projections ranging from −916.2 to 121.2 Tg C and was largely driven by differences in future climate conditions projected by climate models. Moving from a high to a low radiative forcing scenario reduced net ecosystem carbon loss by 21% and when combined with reductions in land‐use change (i.e., moving from a high to a low land‐use scenario), net carbon losses were reduced by 55% on average. However, reconciling large uncertainties associated with the effect of increasing atmospheric CO 2 is needed to better constrain models used to establish baseline conditions from which ecosystem‐based climate mitigation strategies can be evaluated.

California↗

USGS2020

Wildfire and the ecological niche: Diminishing habitat suitability for an indicator species within semi‐arid ecosystems

Globally accelerating frequency and extent of wildfire threatens the persistence of specialist wildlife species through direct loss of habitat and indirect facilitation of exotic invasive species. Habitat specialists may be especially prone to rapidly changing environmental conditions because their ability to adapt lags behind the rate of habitat alteration. As a result, these populations may become increasingly susceptible to ecological traps by returning to suboptimal breeding habitats that were dramatically altered by disturbance. We demonstrate a multistage modeling approach that integrates habitat selection and survival during the key nesting life‐stage of a bird species of high conservation concern, the greater sage‐grouse ( Centrocercus urophasianus ; hereafter, sage‐grouse). We applied these spatially explicit models to a spatiotemporally robust dataset of sage‐grouse nest locations and fates across wildfire‐altered sagebrush ecosystems of the Great Basin ecoregion, western United States. Female sage‐grouse exhibited intricate habitat selection patterns that varied across regional gradients of ecological productivity among sagebrush communities, but often selected nest sites that disproportionately resulted in nest failure. For example, 23% of nests occurred in wildfire‐affected habitats characterized by reduced sagebrush cover and greater composition of invasive annual grasses. We found survival of nests was negatively associated with wildfire‐affected areas, but positively associated with higher elevations with increased ruggedness and overall shrub cover. Strong site fidelity likely drove sage‐grouse to continue nesting in habitats degraded by wildfire. Hence, increasing frequency and extent of wildfire may contribute disproportionately to reduced reproductive success by creating ecological traps that act as population sinks. Identifying such habitat mismatches between selection and survival facilitates deeper understanding of the mechanisms driving reduced geographic niche space and population decline at broad spatiotemporal scales, while guiding management actions to areas that would be most beneficial to the species.

Global Change Biology↗

USGS2022

Microbiome assembly in thawing permafrost and its feedbacks to climate

The physical and chemical changes that accompany permafrost thaw directly influence the microbial communities that mediate the decomposition of formerly frozen organic matter, leading to uncertainty in permafrost–climate feedbacks. Although changes to microbial metabolism and community structure are documented following thaw, the generality of post-thaw assembly patterns across permafrost soils of the world remains uncertain, limiting our ability to predict biogeochemistry and microbial community responses to climate change. Based on our review of the Arctic microbiome, permafrost microbiology, and community ecology, we propose that Assembly Theory provides a framework to better understand thaw-mediated microbiome changes and the implications for community function and climate feedbacks. This framework posits that the prevalence of deterministic or stochastic processes indicates whether the community is well-suited to thrive in changing environmental conditions. We predict that on a short timescale and following high-disturbance thaw (e.g., thermokarst), stochasticity dominates post-thaw microbiome assembly, suggesting that functional predictions will be aided by detailed information about the microbiome. At a longer timescale and lower-intensity disturbance (e.g., active layer deepening), deterministic processes likely dominate, making environmental parameters sufficient for predicting function. We propose that the contribution of stochastic and deterministic processes to post-thaw microbiome assembly depends on the characteristics of the thaw disturbance, as well as characteristics of the microbial community, such as the ecological and phylogenetic breadth of functional guilds, their functional redundancy, and biotic interactions. These propagate across space and time, potentially providing a means for predicting the microbial forcing of greenhouse gas feedbacks to global climate change.

Global Change Biology↗

USGS2012

Foreword: Contributions of Arctic PRISM to monitoring western hemispheric shorebirds

Long-term monitoring of populations is of paramount importance to understanding responses of organisms to global environmental change and to evaluating whether conservation practices are yielding intended results through time (Wiens 2009). The population status of many shorebird species, the focus of this volume, remain poorly known. Long-distance migrant shorebirds have proven particularly difficult to monitor, in part because of their highly inaccessible regions. As migrant shorebirds travel the length of the hemisphere, the congregate and disperse in ways that vary among species, locations, and years, presenting serious challenges to designing and implementing monitoring programs. Rigorous field and quantitative methods that estimate population size and monitor trends are vitally needed to direct and evaluate effective conservation measures. Many management efforts depend on unbiased population size estimates; for examples, the shorebird conservation plans for both Canada and the United States seek to restore populations to levels calculated for the 1970s based on the best information available from existing surveys. Further, federal wildlife agencies within the United States and Canada have mandates to understand the state of their nations' resources under various conventions for the protection of migratory birds. Accurate estimates of population size are vital statistics for a variety of conservation activities, such as prioritizing species for conservation action and setting management targets. Areas of essential habitat, such as those designated under the Western Hemisphere Shorebird Reserve Network, the Important Bird Areas program of BirdLife Internationals and the National Audubon Society, or Canada's National Wildlife Areas program, are all evaluated on the basis of proportions of species' populations which they contain. The size, and trends in size, of a species' population are considered key information for assessing its vulnerability and subsequent listing under the U.S. Endangered Species Act and the Canadian Species at Risk Act. To meet the need for information on population size and trends, shorebird biologists from Canada and the United States proposed a shared blueprint for shorebird monitoring across the Western Hemisphere in the late 1990s; this effort was undertaken in concert with the development of the Canadian and the U.S. Shorebird Conservation Plans. Soon thereafter, partners in the monitoring effort adopted the name "Program for Regional and International Shorebird Monitoring" (PRISM). Among the primary objectives of PRISM were to estimate the population sizes and trends of breeding North American shorebirds and describe their distributions. PRISM members evaluated ongoing and potential monitoring approached to address 74 taxa (including subspecies) and proposed a combination of arctic and boreal breeding surveys, temperate breeding and non-breeding surveys, and neotropical surveys.

Studies in Avian Biology↗

USGS2015

Southern Salish Sea Habitat Map Series: Admiralty Inlet

In 2010 the Environmental Protection Agency, Region 10 initiated the Puget Sound Scientific Studies and Technical Investigations Assistance Program, designed to support research in support of implementing the Puget Sound Action Agenda. The Action Agenda was created in response to Puget Sound having been designated as one of 28 estuaries of national significance under section 320 of the U.S. Clean Water Act, and its overall goal is to restore the Puget Sound Estuary's environment by 2020. The Southern Salish Sea Mapping Project was funded by the Assistance Program request for proposals process, which also supports a large number of coastal-zone- and ocean-management issues. The issues include the recommendations of the Marine Protected Areas Work Group to the Washington State Legislature (Van Cleve and others, 2009), which endorses a Puget Sound and coast-wide marine conservation needs assessment, gap analysis of existing Marine Protected Areas (MPA) and recommendations for action. This publication is the first of four U.S. Geological Survey Scientific Investigation Maps that make up the Southern Salish Sea Mapping Project. The remaining three map blocks to be published in the future, located south of Admiralty Inlet, are shown in figure 1. Puget Sound is a deep, fjord-type estuary covering an area of 2,330 km 2 in the Pacific Northwest region of the United States (fig. 1). It is connected to the ocean by the Strait of Juan de Fuca, a turbulent passage approximately 160 km in length and 22 km wide at its west end, expanding to over 40 km wide at its east end (Thomson, 1994). During the Pleistocene, the area was occupied several times by lobes of continental ice, resulting in a complex basin-fill of glacial and interglacial deposits that are locally as thick as 1100 m (Johnson and others, 2001). The last glaciation, called the Fraser glaciation, began after 28,800±740 14 C yr B.P. when ice started a slow expansion (Clague, 1981). At peak advance the westward Juan de Fuca lobe reached the edge of the continental shelf through the Juan de Fuca Strait shortly before 14,460±200 14 C yr B.P. (Herzer and Bornhold, 1982). The southward Puget lobe advanced to its terminal position in Puget Sound by around 14,150 14 C yr B.P. (Porter and Swanson, 1998). Ice retreated from its maximum to northern Whidbey Island by 13,650±350 14 C yr B.P. (Dethier and others, 1995). Retreating glaciers resulted in a thick sequence of ice-contact, glacial-marine sediment, and early post-glacial sediments (Linden and Schurrer, 1988). These deposits have experienced the effects of a marine transgression followed by regression, resulting in a sea-level several tens of meters lower than the present day (Linden and Schurrer, 1988). A second transgression brought sea level to about the present level by around 5,470±120 14 C yr B.P. (Clague and others, 1982) establishing the present oceanographic and geologic environment Puget Sound is separated into four interconnected basins; Whidbey, Central (Main), Hood Canal, and South (Thomson, 1994). The Whidbey, Central, and Hood Canal basins are the three main branches of the Puget Sound estuary and are separated from the Strait of Juan de Fuca by a double sill at Admiralty Inlet. The Admiralty Inlet map area includes the Inlet and a portion of the Whidbey Basin (fig. 1). The shallower South Basin is separated by a sill at Tacoma Narrows and is highly branched with numerous finger inlets. Flow within Puget Sound is dominated by tidal currents of as much as 1 m/s at Admiralty Inlet, reducing to approximately 0.5 m/s in the Central Basin (Lavelle and others, 1988). The lack of silt and clay-sized sediments in the Admiralty Inlet map area is likely a result of the strong currents (see Ground-Truth Studies for the Admiralty Inlet Map Area, sheet 3). The subtidal component of flow reaches approximately 0.1 m/s and is driven by density gradients arising from the contrast in salty ocean water at the entrance and freshwater inputs from stream flow (Lavelle and others, 1988). The total freshwater input to Puget Sound is approximately 3.4 x 10 6 m 3 /day, primarily from the Skagit River (Cannon, 1983). The subtidal circulation mostly consists of a two-layered flow in the basins with fresher water exiting at the surface and saltier water entering at depth (Ebbesmeyer and Cannon, 2001). In general, surface waters flow north and deeper waters flow south; variations arise from wind effects that can drive a surface current in the same direction as the wind, and a baroclinic response in the lower layer to about 100-m depth (Matsuura and Cannon, 1997). Oceanographic properties are influenced by temporal forcing parameters such as reduced stream flow during the 2000-01 drought that increased surface salinity and decreased differences between surface and bottom waters (Newton and others, 2003). On offshore seismic-reflection profiles, Pleistocene strata (excluding latest Pleistocene glacial and post-glacial deposits) form a distinct seismic unit, bounded below by pre-Tertiary or Tertiary basement and above by typically flat-lying latest Pleistocene to Holocene deposits that fill in erosional or depositional relief (Johnson and others, 2001). Cores from central Puget Sound have accumulation rates that range from 85 to 1200 mg/cm 2 /yr, or 0.12 to 2.4 cm/yr; the highest accumulation rates are near the southern end of central Puget Sound (Carpenter and others, 1985). Carpenter and others (1985) un-weighted arithmetic mean of accumulation rates for central Puget Sound deeper stations is 480±340 (± one standard deviation) mg/cm 2 /yr. Lavelle and others (1985) also found rates as high as 1200 mg/cm 2 /yr over the past approximately 70 years in cores in the Central Basin off of and north and south of Elliott Bay. Puget Sound basin rates are comparable to rates in midshelf silt deposits on the Washington coast north of the Columbia River (Nittrouer and others, 1979). The deep subtidal (in other words, below SCUBA depths) habitats of Puget Sound are relatively poorly known. A few subtidal surveys exist for several habitat types from the 1960s and 1970s (reviewed in Dethier, 1990), using grab and box core data. The Dethier (1990) review divides habitat up into Coast and Marine Ecological Classification Standard (CMECS) substrate, water column energy, and depth zones but does not attempt to map these habitats, rather it is an inventory of habitats found in the area and the flora and fauna associated with each habitat. The approach of the Southern Salish Sea Mapping project is to create highly detailed seafloor maps through collection, integration, interpretation, and visualization of swath sonar data (the undersea equivalent of satellite remote-sensing data in terrestrial mapping), acoustic backscatter, seafloor video, seafloor photography, and bottom-sediment sampling data. This approach is based in part on methods presented and data collection and product needs identified at the Washington State Seafloor Mapping Workshop (Washington State Seafloor Mapping Workshop Steering Committee, 2008), attended by coastal and marine managers and scientists. The map products display seafloor geomorphology and substrate, and identify potential marine benthic habitats. It is emphasized that the more interpretive habitat and geology maps rely on the integration of multiple, new high-resolution datasets and that mapping at small scales would not be possible without such data. Oceanographic current and wave data is not included in this analysis, however, the accompanying geographic information system (GIS) data set is designed and intended to be combined with oceanographic and biologic data sets assembled by others in the future and some of the GIS data has already been incorporated in the unpublished Nature Conservancy Benthic Habitats of Puget Sound database. This publication includes four map sheets, explanatory text, and a descriptive pamphlet. Each map sheet is published as a portable document format (PDF) file. ESRI ArcGIS compatible geotiffs (for example, bathymetry) and shapefiles (for example video observation points) will be available for download in the data catalog associated with this publication (Cochrane, 2015). An ArcGIS Project File with the symbology used to generate the map sheets is also provided. For those who do not own the full suite of ESRI GIS and mapping software, the data can be read using ESRI ArcReader, a free viewer that is available at http://www.esri.com/software/arcgis/arcreader/index.html .

Washington↗

USGS2014

Synthesis of studies in the fall low-salinity zone of the San Francisco Estuary, September-December 2011

In fall 2011, a large-scale investigation (fall low-salinity habitat investigation) was implemented by the Bureau of Reclamation in cooperation with the Interagency Ecological Program to explore hypotheses about the ecological role of low-salinity habitat in the San Francisco Estuary—specifically, hypotheses about the importance of fall low-salinity habitat to the biology of delta smelt Hypomesus transpacificus , a species endemic to the San Francisco Estuary and listed as threatened or endangered under federal and state endangered species legislation. The Interagency Ecological Program is a consortium of 10 agencies that work together to develop a better understanding of the ecology of the Estuary and the effects of the State Water Project and Federal Central Valley Project operations on the physical, chemical, and biological conditions of the San Francisco Estuary. The fall low-salinity habitat investigation constitutes one of the actions stipulated in the Reasonable and Prudent Alternative issued with the 2008 Biological Opinion of the U.S. Fish and Wildlife Service, which called for adaptive management of fall Sacramento-San Joaquin Delta outflow following “wet” and “above normal” water years to alleviate jeopardy to delta smelt and adverse modification of delta smelt critical habitat. The basic hypothesis of the adaptive management of fall low-salinity habitat is that greater outflows move the low-salinity zone (salinity 1–6), an important component of delta smelt habitat, westward and that moving the low-salinity zone westward of its position in the fall of recent years will benefit delta smelt, although the specific mechanisms providing such benefit are uncertain. An adaptive management plan was prepared to guide implementation of the adaptive management of fall low-salinity habitat and to reduce uncertainty. This report has three major objectives: • To provide a summary of the results from the first year of coordinated fall low-salinity habitat studies and monitoring. • To provide a synthesis of the results of the fall low-salinity habitat studies and other ongoing research and monitoring, to determine if the available information supports the hypotheses behind the adaptive management of fall low-salinity habitat as set forth in the adaptive management plan. • To begin to put the results from the fall low-salinity habitat studies into context within the larger body of knowledge regarding the San Francisco Estuary and, in particular, the upper San Francisco Estuary, including the Sacramento-San Joaquin Delta, Suisun Bay, and associated embayments. The basic approach of this report is to evaluate predictions derived from the hypotheses included in the conceptual model developed within the adaptive management plan. All available data from studies and monitoring conducted in fall 2011 and similar data from fall 2006, which was the most recent wet year preceding 2011, were considered. Data from 2005 and 2010 were also considered, to include the conditions antecedent to those years. Many of the predictions either could not be evaluated with the data available, or the needed data were not collected. Most of the predictions that could be addressed involved either the abiotic habitat components (that is, the physical environment) or delta smelt responses. In general, the fall low-salinity habitat investigation has been largely inconclusive as of the writing of this report. This is not to be unexpected in the first year of what is intended to be a multi-year adaptive-management effort. This report can be viewed as the first chapter of a “living document” that is to be continually updated as part of the adaptive management cycle. The results of this report, especially predictions with insufficient data for evaluation, indicate a number of science-based approaches to improve the fall low-salinity habitat investigations: • Develop a method of measuring “hydrodynamic complexity.” This concept is central to a number of the predictions that could not be evaluated. • Determine if wind speed warrants a stand-alone prediction. The wind-speed prediction is directly related to the turbidity predictions, and wind is only one of several factors important for determining turbidity. • Determine the correct spatial and temporal scale or scales necessary for monitoring and for studies to address the predicted abiotic and biotic responses. Many of the assessments in this report were based on monthly sampling of dynamic habitat components, such as phytoplankton and zooplankton populations, that can change on daily scales. • Address the nutrient predictions as part of developing a phytoplankton production model that includes nutrient cycling and other important processes, if feasible. At a minimum develop a mechanistic conceptual model to support more processed-based interpretations of data or design of new studies, rather than making simple predictions of increase or decrease. • Determine if studies of predation rates are feasible in areas where there are delta smelt.

California↗

USGS2020

Modeling Escherichia coli in the Missouri River near Omaha, Nebraska, 2012–16

The city of Omaha, Nebraska, has a combined sewer system in some areas of the city. In Omaha, Nebr., a moderate amount of rainfall will lead to the combination of stormwater and untreated sewage or wastewater being discharged directly into the Missouri River and Papillion Creek and is called a combined sewer overflow (CSO) event. In 2009, the city of Omaha began the implementation of their Long Term Control Plan (LTCP) to mitigate the effects of CSOs on the Missouri River and Papillion Creek. As part of the LTCP, the city partnered with the U.S. Geological Survey (USGS) in 2012 to begin monitoring in the Missouri River. Since 2012, monthly discrete water-quality samples for many constituents have been collected from the Missouri River at four sites. At 3 of the 4 sites, water quality has been monitored continuously for selected constituents and physical properties. These discrete water-quality samples and continuous water-quality monitoring data (from July 2012 to 2020) have been collected to better understand the water quality of the Missouri River, how it is changing with time, how it changes upstream from the city of Omaha to downstream, and how it varies during base-flow conditions and during periods of runoff. The purpose of this report is to document the development of Escherichia coli ( E. coli ) concentration models for these four Missouri River sites. Analysis was completed using the first 5 years of data (through 2016) to determine if the current approach is sufficient to meet future analysis goals and to understand if proposed models such as Load Estimator (LOADEST) models will be able to represent water-quality changes in the Missouri River. Multiple linear regression models were developed to estimate E. coli concentration using LOADEST as implemented in the rloadest package in the R statistical software program. A set of explanatory variables, including streamflow and streamflow anomalies, precipitation, information about CSOs, and continuous water quality, were evaluated for potential inclusion in regression models. The best model at Missouri River at NP Dodge Park at Omaha, Nebr. (USGS station 412126095565201; hereafter “NP Dodge”) included basin explanatory variables of upstream antecedent precipitation index measured at Tekamah, Nebr.; decimal time; season; and turbidity. The best model at Missouri River at Freedom Park Omaha, Nebr. (USGS station 411636095535401; hereafter “Freedom Park”) included the same explanatory variables as the NP Dodge model with the addition of turbidity anomalies and flow anomalies. The best models at the two downstream sites (Missouri River near Council Bluffs, Iowa, USGS station 06610505 and Missouri River near La Platte, Nebr., USGS station 410333095530101) included the same explanatory variables as the Freedom Park model with the addition of local antecedent precipitation index as measured at Eppley Airport in Omaha, Nebr., and additional turbidity and flow anomalies. The final selected models were the best models given our modeling design constraint in which explanatory variables included in the model for the upstream site were included in the downstream models. Explanatory variables currently (2020) being collected and included in the selected models through 2016 explained 64–75 percent of the variability of E. coli concentration in the Missouri River. Explaining 64–75 percent of the variability might be considered low when working with physical constituents (total nitrogen or sediment), but with the natural variability of biological constituents such as E. coli , the uncertainty of E. coli laboratory measurements, and the added complexity of modeling in a large drainage basin with multiple sources, these results are adequate and indicate that the explanatory variables being collected and models such as LOADEST can represent water-quality changes in the Missouri River for E. coli concentration from 2012 to 2016.

Nebraska↗