Search USGSSearch

Geology topics

Michael J. Dodrill

Publications and source records attributed to Michael J. Dodrill.

At least 19 recordsLinked to original sources

DRIFteRS: A dataset of drift invertebrate densities in streams and rivers across western North America, 1997–2024

Prey availability is among the most influential and highly variable determinants of fish growth and freshwater habitat carrying capacity, yet it remains understudied compared to physical habitat variables (Ouellet et al., 2025; Rosenfeld et al., 2014; Weber et al., 2017). We often lack a clear understanding of how much food is available to fishes, how it varies spatially and temporally, and how it influences responses to restoration (Ouellet et al., 2025; Rossi et al., 2024; Wipfli et al., 2010). Drift invertebrates—the primary food source for juvenile salmonids and other drift-foraging fishes—play a pivotal role in these dynamics. To better understand the spatiotemporal variability of drift invertebrate abundance and biomass across the freshwater range of drift-feeding salmonids in western North America, we compiled the DRift Invertebrates For salmonids in River Systems (DRIFteRS) dataset. The dataset encompasses 6125 samples of drift invertebrates, and, for a subset of drift samples, associated benthic invertebrate density data, collected from 1360 reaches on 459 unique rivers and streams spanning 55 river basins considered hydrologically independent (i.e., not nested within the same larger watershed) across British Columbia, Canada, and the U.S. states of Alaska, Arizona, California, Colorado, Idaho, Nevada, New Mexico, Oregon, Utah, Washington, and Wyoming. Sample sites represent a diverse array of river and stream habitats (e.g., headwater, mainstem, side channel), in watersheds with diverse land uses (e.g., urban, wilderness, agricultural), and disturbance histories (e.g., fire, restoration). Collected between 1997 and 2024, the data span the full calendar year and capture daily and seasonal patterns in drift abundance and biomass densities. When paired with water quality and quantity data as well as remotely sensed environmental landscape data, such as land use/land cover, climate, and disturbance history, channel morphology, and riparian vegetation composition, the DRIFteRS dataset can aid in identifying key drivers of drift invertebrate densities and mean body size and support predictive modeling in unsampled locations and times. The dataset may also be used to analyze aquatic-terrestrial resource flows, derive prey-encounter rates and profitability (mean prey size), and inform broader investigations of sit-and-wait foraging ecology, especially when paired with data on drift-foraging predators. For salmonid-focused applications, the dataset can be integrated into habitat evaluation models, including bioenergetic (e.g., Naman et al., 2019) and life cycle models (e.g., Beechie et al., 2023), to improve estimates of habitat capacity and population dynamics for river- and stream-rearing salmonids. Understanding prey availability dynamics is increasingly important, because rising water temperatures increase salmonid metabolic demands (Crozier et al., 2010). Flow regime transitions (i.e., snow or glacier dominated to rain dominated; Beechie et al., 2013), wildfire frequency and intensity (Hessburg et al., 2021), as well as plant community and phenology shifts (Cleland et al., 2007; Franklin et al., 2016) are all predicted to change with rising temperatures and are potential drivers of terrestrial and aquatic invertebrate prey quality and availability in lotic systems. These insights can ultimately inform restoration prioritization and design, helping managers consider food resource implications when evaluating restoration priorities and match habitat improvement to food supply. The data are released under the Creative Commons Attribution 4.0 International license and may be reused with attribution.

Alaska, Arizona, British Columbia, California, Col

Restoration is a partial but not complete solution for long-term declines in oyster populations—A case history from the northeastern Gulf of America

Objective We sought to evaluate whether oyster reef restoration is associated with sustained increases in intertidal eastern oyster Crassostrea virginica density in Suwannee Sound, Florida, and to examine how commercial fishing activity and freshwater discharge may modify those patterns. Methods We monitored intertidal oyster reefs for 14 years, including 9 years prior to large-scale restoration (completed in 2018) and 5 years of standardized postrestoration monitoring beginning in 2019 at Lone Cabbage Reef, an offshore intertidal reef complex that was rebuilt to target elevation using dolomite limestone. A Bayesian negative-binomial generalized linear mixed model framework estimated oyster densities while accounting for spatial differences among inshore, nearshore, and offshore reefs; variation in transect length; and temporal variability. We assessed restoration effectiveness via repeated-measures monitoring of fixed restored sites and landscape-scale comparisons of restored versus unrestored reefs. We evaluated covariates of management interest—commercial fishing activity and river discharge—using model comparison and posterior predictions of oyster density at observed covariate values. Results Restoration produced large per-area gains in live-oyster density on offshore reefs. At four fixed restored offshore stations, mean densities rose from near 0 to 91 oysters/m 2 (95% credible interval [CrI] = 43–185) within 6 months and reached 244 oysters/m 2 (95% CrI = 134–433) by year 5 at fixed restoration stations. At the broader landscape scale, restored offshore reefs averaged about 83 oysters/m 2 (95% CrI = 26–157) across sites and postrestoration periods, whereas model-based estimates for restored offshore reefs reached 218 oysters/m 2 (95% CrI = 128–366) in the final sampling period. Conclusions Restoring reef substrate rebuilt the physical capacity for oysters to persist and grow at Lone Cabbage Reef, yielding sustained density gains that degraded sites did not achieve on their own. Harvest status and river discharge showed weak, inconsistent associations once restoration was included, consistent with the interpretation that oyster cultch provides structural capacity, while environmental and fishing conditions modulate around it. Rebuilding and conserving oyster cultch constitute the essential first step to repair and sustain degraded oyster reefs; fishing regulations or favorable environmental conditions alone are unlikely to replace the need for reef structure.

Florida

Assessment of Chinook Salmon smolt survival at offsite release locations in the Sacramento River across wet and dry years

We used acoustic telemetry to estimate survival of tagged release groups, and quantify differences between alternative hatchery release strategies. To assess whether offsite release could increase survival of hatchery fish relative to those released at the hatchery, we compare survival during emigration of hatchery fall Chinook Salmon (FCS, Oncorhynchus tshawytscha ) smolts released onsite, at Coleman National Fish Hatchery (NFH), to those released at alternative offsite locations downstream in the upper Sacramento River, California. Approximately 300 fish in each release group were implanted with acoustic tags in the 3-year study during 2019, 2021, and 2022. Environmental conditions in the emigration corridor varied between study years, with extended high flows in 2019, whereas drought conditions and a flat hydrograph occurred in 2021 and 2022. Survival across years appeared to reflect environmental conditions, with higher overall survival seen in 2019. Survival differences between release sites generally showed higher cumulative survival to Knights Landing or Chipps Island for the downstream release groups, compared to the onsite releases during 2021 and 2022. This information suggests hatchery releases of FCS from Coleman NFH at downstream release sites may benefit emigration survival during years of drought and sub-optimal in-river conditions, compared to standard onsite releases. Hatchery managers can use information from this study to diversify the portfolio of release strategies, which may be useful to buffer against extreme variations of adult abundance by spreading risks across space and time. Ultimately, offsite releases may be a useful tool for adaptive management of these culturally, ecologically, and economically important salmon populations.

California

Who needs closure? Estimating abundance with a Markovian availability model for geographically open removal sampling

Removal sampling is an important method for estimating abundance, but nearly all removal models assume closure during sampling. Yet, closure may be difficult to assume, evaluate, or enforce in many settings. To address situations where populations are geographically open between each removal sample, we incorporated a Markovian availability process into an N-mixture model framework. This model relates local abundance available for sampling to a superpopulation through recruitment of new individuals to the sampling area. To test the model, we (1) conducted parameter identifiability analysis, (2) fit the model to removal data generated from a random walk movement model, and (3) analyzed a case study of empirical removal data. Parameters were increasingly identifiable as capture probability exceeded 0.25 and removal samples increased from 3 to 6. Abundance estimates were unbiased when parameters were identifiable, except for scenarios that simulated a behavioral response to sampling. For our case study, the model estimated negligible recruitment for benthic-oriented fishes, indicating closure, but we found evidence against closure for juvenile Chinook salmon, a highly mobile species. Our removal model allows researchers to formally test closure assumptions, to estimate the degree of closure, and to estimate abundance without bias when closure is violated.

California

U.S. Geological Survey Grand Canyon Monitoring and Research Center: Proceedings of the fiscal year 2023 annual reporting meeting to the Glen Canyon Dam Adaptive Management Program

This proceedings is prepared for the USBR and Glen Canyon Dam Adaptive Management Program (GCDAMP) to account for work conducted and products delivered in FY 2023 by SBSC's Grand Canyon Monitoring and Research Center (GCMRC) and to inform the Technical Work Group of science conducted by GCMRC and its cooperators in support of the GCDAMP. It includes a summary of accomplishments, modifications to work plans, and results related to projects included in GCMRC’s FY 2021-23 Triennial Work Plan. This work was done to support the 11 resource goals identified in the Glen Canyon Dam Long-Term Experimental and Management Plan (LTEMP) Environmental Impact Statement and Record of Decision.

Arizona, Nevada, Utah

Insectivorous bat foraging tracks the availability of aquatic flies (Diptera)

Rivers and their adjacent riparian zones are model ecosystems for observing cross-ecosystem energy transfers. Aquatic insects emerging from streams, for example, are resource subsidies that support riparian consumers such as birds, spiders, lizards, and bats. We collaborated with recreational river runners in Grand Canyon, Arizona, USA, to record acoustic bat activity and sample riparian insects using light traps at dusk from April through October 2017–2020. River runners collected these data on 1,428 events over 611 sampling nights at 410 sites throughout a 470-km segment of the Colorado River. We documented 71 insect taxa in light traps and recorded 19 bat species with acoustic detectors. We hypothesized that bat activity along this highly regulated river segment would be influenced primarily by variation in prey availability, as compared to other habitat descriptors. We predicted that bat activity would be positively related to aquatic insect catch rates and unrelated to terrestrial insect abundance. We fit Bayesian regression models to test these hypotheses and to quantify the relationship between bat activity and a suite of environmental variables: time of year, time of day, distance from perennial tributaries, distance from rapids, channel width, geomorphic reach, tall vegetation cover, air temperature, and lunar phase. Bat activity was positively related to the abundance of aquatic flies (Diptera), which outcompeted all other prey categories and structural habitat descriptors in our models. Within our dusk sampling period, activity of small myotis (California myotis [ Myotis californicus ] and Yuma myotis [ M. yumanensis ]) was high late in the evening and canyon bats ( Parastrellus hesperus ), conversely, were more active early in the evening. Activity of canyon bats varied seasonally, with peak activity in August. Our results support the hypothesis that aquatic prey, specifically aquatic flies, are key predictors of bat activity along a large, regulated river corridor and demonstrate the power of community science as a tool for ecosystem monitoring.

Journal of Wildlife Management

Declines in prey production during the collapse of a tailwater Rainbow Trout population are associated with changing reservoir conditions

Objective Understanding how energy moves through food webs and limits productivity at various trophic levels is a central question in aquatic ecology and can provide insight into drivers of fish population dynamics since many fish populations are food limited. In this study, we seek to better understand what factors drove a decline of >85% in the number of Rainbow Trout Oncorhynchus mykiss found in the tailwater portion of the Colorado River below Glen Canyon Dam during 2012–2016. Methods We estimate the production of dominant prey using data from previously published studies of Rainbow Trout abundance and growth alongside drift and diet samples. We test how prey production correlates to both proximate (e.g., nutrients) and distal (e.g., limnological conditions in the upriver reservoir) drivers. Result Results suggest that gross consumption of invertebrate prey by the Rainbow Trout population declined from an annual mean of 423 to 69 kg/d. Daily production rates of dominant prey in aggregate declined from a high of 0.173 to 0.018 g·m −2 ·d −1 . Chironomids accounted for 70% of the decline in prey production. Foraging efficiency by Rainbow Trout (range, 0.99–0.67) was high across the range of prey production rates. After the Rainbow Trout population had declined by ~90%, prey consumption saturated at higher rates of prey production and the gross quantity of daily drift exported from the reach increased from 8.9 to 12.7 kg/d. Conclusion Rainbow Trout population dynamics are largely influenced by changes in prey production, which is itself driven by soluble reactive phosphorus ( SRP ) concentrations in the reservoir. The SRP model predicted that prey production would increase by 32 kg/d (SE, 9) for each 1 μg/L increase in SRP . These concentrations were indirectly influenced by reservoir hydrology and biogeochemistry, linkages that may extend far beyond the confines of this tailwater fishery and into the downstream reaches of the Grand Canyon's Colorado River ecosystem. Impact Statement We combined Rainbow Trout diet, growth, and abundance estimates with concentrations of drifting invertebrates to estimate the biomass of Rainbow Trout prey produced over time. Trends in prey biomass production track trends in phosphorous concentrations in the river.

Arizona

Simulating post-dam removal effects of hatchery operations and disease on juvenile Chinook salmon (Oncorhynchus tshawytscha) production in the Lower Klamath River, California

Executive Summary The Federal Energy Regulatory Commission has been considering the approval to breach four dams on lower Klamath River in southern Oregon and northern California. Approval of this application would allow for Strikeouts indicate text deletion hereafter. decommissioning and dam removal, beginning as early as 2023. This action would affect Klamath River salmon ( Oncorhynchus ssp.) populations, a critical food source for federally endangered Southern Resident Killer Whales ( Orcinus orca ). In the long run, reintroduction of salmon populations to the upper Klamath River Basin may increase salmon abundance available to Southern Resident Killer Whales, but in the near term, it is uncertain how changes in hatchery management and disease-caused mortality by the myxosporean parasite Ceratonova shasta will influence abundance of salmon populations entering the ocean. To assess this uncertainty, we used the Stream Salmonid Simulator (S3) to simulate population dynamics of juvenile Chinook salmon ( Oncorhynchus tshawytscha ) for nine different population sources that rear and migrate through the Klamath River. S3 is a spatially explicit population model that runs on a daily time-step and simulates daily growth, survival, and movement of juvenile Chinook salmon from the time of spawning through ocean entry. The key features of this model relevant to this report include (1) a C. shasta disease submodel; (2) a temperature-dependent bioenergetics model that calculates daily growth rates; (3) size-dependent movement; (4) density-dependent dynamics that are influenced by the effect of flow on suitable habitat area; and (5) habitat, river flow, and water temperature specific to each scenario. We constructed and ran four scenarios: two scenarios for dams in place (Dams In) and dams removed (Dams Out), and given these dam-removal conditions, a low- and high-spore scenario for C. shasta . Each scenario was run for nine water years representing a range of conditions from dry to wet. Previously published daily river flows and water temperatures for Dams In and Dams Out provided physical inputs for each scenario. Daily spore concentrations were simulated using a three-part mechanistic model that used river discharge, water temperature, and the prevalence of infection (POI) of hatchery-origin Chinook salmon juveniles with C. shasta in the previous year 3 . We constructed two spore scenarios for each Dams In and Dams Out scenario, a “Low Spore” scenario and a “High Spore” scenario resulting in four scenarios for comparison. Spore scenarios were established by setting the prior-year POI of hatchery fish to 0.15 and 0.75 in the estimation of spore concentrations. Hatchery releases under Dams Out differed from those under the current Dams In scenario. Hatchery releases under the Dams Out scenario were modified to emulate changes in hatchery production that would occur under Dams Out conditions. This included moving hatchery production and releases from Iron Gate Dam to a proposed hatchery at Fall Creek, which would be located about 11 kilometers (km) upstream of Iron Gate Dam. It is anticipated that the Fall Creek hatchery would produce fewer fish at smaller and larger sizes at different release timings. For salmon inputs, we used observed historical abundance of main-stem spawners from brood year 2009 and juvenile salmon entering from tributaries in water year 2010, which represented an average return year for the 2005–18 period. Main-stem spawning was allowed to shift upstream from Iron Gate Dam under the Dams Out scenario. We also included hatchery-origin fish as natural spawners that would have otherwise returned to Iron Gate Hatchery in the first 3 years following dam removal. The S3 model simulated considerably higher total abundance for Dams Out relative to the respective Dams In scenarios, and higher abundance for the Low Spore scenario relative to the High Spore scenario. The difference in abundance between the four combinations of the dam-removal and spore scenarios varied among population groups. For main-stem natural production, juvenile abundance at ocean entry was 2–3 times higher for Dams Out scenarios than for Dams In scenarios, and juvenile abundance for High Spore scenarios was lower than that for the Dams Out Low Spores scenario. For hatchery releases, abundance at ocean entry was similar between Dams In and Dams Out scenarios for most water years, despite lower release sizes from Fall Creek Hatchery under Dams Out. For tributary populations, abundance for the High Spore scenarios was consistently lower than for the Low Spore scenarios, but differences between dam-removal scenarios varied among water years, with Dams Out scenarios having similar or higher abundance than Dams In scenarios, and dry water years having the largest difference between Dams In and Dams Out scenarios. We determined that different factors affected the response of each population group. For main-stem natural production, survival from fry emergence to ocean entry was higher under Dams Out scenarios compared to Dams In scenarios because juveniles emerged later and tended to arrive at the ocean sooner and at larger sizes, causing the population to have less time-dependent in-river mortality. Owing to their late release timing, hatchery populations had high disease-caused mortality in Dams In and Dams Out High Spore scenarios. Furthermore, a high proportion of infected fish (those that would be expected to die at some future point) survived to the ocean. Iron Gate Hatchery fish had lower survival rates than releases from Fall Creek Hatchery because the last mid-June release group from the 2010 Iron Gate Hatchery release incurred nearly total mortality in most water years owing to water temperatures exceeding 24 degrees Celsius. Our analysis shows how the S3 model was able to track different populations and provide insights on how the differential response of each population combined to influence the simulated number of juvenile Chinook salmon arriving at the Pacific Ocean where they become available as a food source for Southern Resident Killer Whales.

California

Extending the Stream Salmonid Simulator to accommodate the life history of coho salmon (Oncorhynchus kisutch) in the Klamath River Basin, Northern California

In this report, we apply the stream salmonid simulator (S3) to coho salmon ( Oncorhynchus kisutch ) in the Klamath River Basin by extending the original model to account for life history and disease dynamics specific to coho salmon. This version of S3 includes tracking of three separate life-history strategies representing the different time periods and ages at which fish leave natal tributaries such as the Scott and Shasta Rivers (age-0 spring, age-0 fall, or age-1 smolt). Once fish leave their natal tributaries and enter the Klamath River, the deterministic life-stage-structured population model simulates daily growth, movement, and survival. We extend the model to include non-natal tributary dynamics, where spring age-0 fish entry to non-natal tributaries is simulated based on environmental conditions in the main-stem Klamath River. Fish that use non-natal tributaries then reenter the Klamath River during the winter or spring as smolts and actively migrate downstream. We also consider the life history strategy where fish rear in natal tributaries and enter the Klamath River as age-1 smolts. In addition to simulating different life history pathways that coho salmon may take, we model disease dynamics, incorporating new information on Ceratonova shasta related infection and mortality. We incorporate competitive interactions between juvenile coho and Chinook salmon ( Oncorhynchus tshawytscha ) by simulating density-dependent movement dynamics in response to Chinook salmon abundance. Model simulations suggest that total abundance and survival to the ocean differed between life-history strategies. In general, spring age-0 fish that leave their natal tributaries in their first spring had lower survival compared with fish that remained in natal tributaries and out-migrated later. Spring age-0 fish also had higher disease related mortality, owing to their residence in the main-stem Klamath River overlapping with periods of elevated C. shasta spore concentrations. Age-0 fish leaving their natal tributaries in the fall had near-zero disease related mortality. Most non-natal tributary use occurred at upstream tributary locations and was variable between the brood years depending on passage timing and environmental conditions. The inclusion of Chinook salmon in simulations resulted in decreased abundance and survival of Coho salmon reaching the ocean. In addition, we developed an R package to facilitate use of and continued development of S3 as a tool to guide management of juvenile salmonid populations.

California

Experimental reductions in sub-daily flow fluctuations increased gross primary productivity for 425 river kilometers downstream

Aquatic primary production is the foundation of many river food webs. Dams change the physical template of rivers, often driving food webs toward greater reliance on aquatic primary production. Nonetheless, the effects of regulated flow regimes on primary production are poorly understood. Load following is a common dam flow management strategy that involves sub-daily changes in water releases proportional to fluctuations in electrical power demand. This flow regime causes an artificial tide, wetting and drying channel margins and altering river depth and water clarity, all processes that are likely to affect primary production. In collaboration with dam operators, we designed an experimental flow regime whose goal was to mitigate negative effects of load following on ecosystem processes. The experimental flow contrasted steady-low flows on weekends with load following flows on weekdays. Here, we quantify the effect of this experimental flow on springtime gross primary production (GPP) 90-to-425 km downstream of Glen Canyon Dam on the Colorado River, AZ, USA. GPP during steady-low flows was 41% higher than during load following flows, mostly owing to non-linear reductions in sediment-driven turbidity. The experimental flow increased weekly GPP even after controlling for variation in weekly mean discharge, demonstrating a negative effect of load following on GPP. We estimate that this environmental flow increased springtime carbon fixation by 0.27 g C m –2 d –1 , which is ecologically meaningful considering median C fixation in 356 U.S. rivers of 0.44 g C m –2 d –1 and the fact that native fish populations in this river are food-limited.

PNAS Nexus

Quantifying the effects of tides, river flow, and barriers on movements of Chinook Salmon smolts at junctions in the Sacramento–San Joaquin River Delta using multistate models

Successful migration of Chinook Salmon ( Oncorhynchus tshawytscha ) smolts seaward in the Sacramento – San Joaquin River Delta (hereafter, Delta) requires navigating a network of numerous branching channels. Within the Delta, several key junctions route smolts either towards more direct paths to the ocean or towards the interior Delta, an area associated with decreased survival. Movements within these junctions that determine route choice can be influenced by numerous behavioral and environmental factors, including the complex interplay between tidal and riverine hydraulics. Here, we apply continuous time multistate Markov models to examine the influence of tidal and riverine hydraulics, behavioral factors, and management actions on smolt movements. These models incorporate more information from acoustic telemetry data compared with previous approaches to modeling smolt movements in the Delta. By decomposing modeled flows into tidal and net flow signals we elucidate how each component influences movements into and out of distributary channels. Increasing net flows generally increased movement rates, while flood tides decreased seaward movement rates. Similarly, ebb tides increased downstream movements as fish go with the flow. We found less support for diel movement behaviors compared to flow metrics. Additionally, we quantify the effects of a large management action, the placement of a physical barrier, which was effective at decreasing entrainment into the interior Delta. Together, these results help inform management of Chinook Salmon and increase our understanding of the major factors driving smolt movements within these key junctions.

California

As the prey thickens: Rainbow trout select prey based upon width not length

Drift-feeding fish are typically considered size-selective predators. Yet, few studies have explicitly tested which aspect of prey “size” best explains size selection by drift-foraging fish. Here, we develop a Bayesian discrete choice model to evaluate how attributes of both prey and predator simultaneously influence size-selective foraging. We apply the model to a large dataset of paired invertebrate drift ( n = 784) and rainbow trout ( Oncorhynchus mykiss ) diets ( n = 1028). We characterized prey “size” using six metrics (length, width, area, hemispherical area, volume, mass) and used pseudo- R 2 to determine which metric best explained observed prey selection across seven taxa. We found that rainbow trout are positively size-selective, they are selecting prey based upon differences in prey width, and size-selectivity increases with fish length. Rainbow trout demonstrated strong selection for the adult and pupae stages of aquatic insects relative to their larval stages. Our study provides strong empirical evidence for size-selective foraging in rainbow trout and demonstrates prey selection is based primarily upon width, not length or area as has been widely reported.

Arizona

Changes in prey, turbidity, and competition reduce somatic growth and cause the collapse of a fish population

Somatic growth exerts strong control on patterns in the abundance of animal populations via effects on maturation, fecundity, and survival rates of juveniles and adults. In this paper, we quantify abiotic and biotic drivers of rainbow trout growth in the Colorado River, AZ, and the resulting impact on spatial and temporal variation in abundance. Inferences are based on approximately 10,000 observations of individual growth grates obtained through an intensive mark‐recapture effort conducted over five years (2012‐2016) in a 130 km‐long study segment downstream of Glen Canyon Dam. Prey availability, turbidity‐driven feeding efficiency, and intra‐specific competition were the dominant drivers of rainbow trout growth. Discharge, water temperature, and solar insulation were also evaluated but had a smaller influence. Mixed‐effect models explained 79‐82% of the variability in observed growth rates, with fixed covariate effects explaining 79‐87% of the total variation in growth parameters across five reaches and 18 quarterly sampling intervals. Reductions in growth owing in part to a phosphorous‐driven decline in prey availability, led to substantive weight loss and poor fish condition. This in turn lowered survival rates and delayed maturation, which led to a rapid decline in abundance and later recruitments. Reductions in feeding efficiency, due to episodic inputs of fine sediment from tributaries, and warmer water temperatures, contributed to reduced growth in downstream reaches, which led to more severe declines in abundance. Somatic growth rates increased following the population collapse due to reduced competition, and in the absence of substantive increases in prey availability. Our study elucidates important linkages between abiotic and biotic factors, somatic growth, and vital rates, and demonstrates how variation in somatic growth influences temporal and spatial patterns in abundance.

Arizona

Trends in oyster populations in the northeastern Gulf of Mexico: An assessment of river discharge and fishing effects over time and space

Within the Big Bend region of the northeastern Gulf of Mexico, one of the least developed coastlines in the continental USA, intertidal and subtidal populations of eastern oyster Crassostrea virginica (hereafter referred to as “oyster”) are a critical ecosystem and important economic constituent. We assessed trends in intertidal oyster populations, river discharge, and commercial fishing activity in the Suwannee River estuary within the Big Bend region using fisheries‐independent data from irregular monitoring efforts and publicly available environmental data. We used generalized linear models to evaluate counts of oysters from line‐transect surveys over time and space. We assessed model performance using simulation to understand potential bias and then evaluated whether these counts were related to freshwater inputs from the Suwannee River and commercial oyster fishing effort and landings at different time lags. We found that intertidal oyster counts have declined over time and that most of these declines are found in inshore intertidal oyster bars, which are becoming degraded. We also found a significant relationship between oyster counts and a 1‐year lag on mean daily Suwannee River discharge, but including commercial fishery trips or landings did not improve model fit. It is unclear whether declines in intertidal oyster bars are offset by formation of new oyster reefs elsewhere. These results quantify rapid declines in intertidal oyster reefs in a region of coastline with high conservation value that can be used to inform ongoing and proposed restoration projects in the region.

Florida

A need for speed in Bayesian population models: A practical guide to marginalizing and recovering discrete latent states

Bayesian population models can be exceedingly slow due, in part, to the choice to simulate discrete latent states. Here, we discuss an alternative approach to discrete latent states, marginalization, that forms the basis of maximum likelihood population models and is much faster. Our manuscript has two goals: 1) to introduce readers unfamiliar with marginalization to the concept and provide worked examples, and 2) to address topics associated with marginalization that have not been previously synthesized and are relevant to both Bayesian and maximum likelihood models. We begin by explaining marginalization using a Cormack-Jolly-Seber model. Next, we apply marginalization to multistate capture-recapture, community occupancy, and integrated population models and briefly discuss random effects, priors, and pseudo-R2. Then, we focus on recovery of discrete latent states, defining different types of conditional probabilities and showing how quantities such as population abundance or species richness can be estimated in marginalized code. Lastly, we show that occupancy and site abundance models with auto-covariates can be fit with marginalized code with minimal impact on parameter estimates. Marginalized code was anywhere from five to >1000 times faster than discrete code. Differences in inferences were minimal using marginalized code. Discrete latent states and fully conditional approaches provide the best estimates of conditional probabilities for a given site or individual. However, estimates for parameters and derived quantities such as species richness and abundance were minimally affected by marginalization and use of imperfect estimates of conditional probabilities. The results applied even when auto-covariates based on imperfect estimates of conditional probabilities were used. Understanding how marginalization works shrinks the divide between Bayesian and maximum likelihood approaches to population models. Some models that have only been presented in a Bayesian framework can easily be fit in maximum likelihood. On the other hand, factors such as informative priors, random effects, or pseudo-R2 values may motivate a Bayesian approach in some applications. An understanding of marginalization allows users to minimize the speed that is sacrificed when switching from a maximum likelihood approach. Widespread application of marginalization in Bayesian population models will facilitate more thorough simulation studies, comparisons of alternative model structures, and faster learning.

Ecological Applications

Colorado River fish monitoring in Grand Canyon, Arizona; 2002–14 humpback chub aggregations

The humpback chub ( Gila cypha ) is an endangered cyprinid species endemic to the Colorado River. The largest remaining population of the species spawns and rears in the Little Colorado River in Grand Canyon. Construction and operation of Glen Canyon Dam has altered the main-stem Colorado River in Glen and Grand Canyons. Cold, clear water releases from the dam result in a river that is generally unsuitable for successful humpback chub reproduction. During the early 1990s, nine locations within the main-stem Colorado River were identified as humpback chub aggregations—areas with a consistent and disjunct group of fish with no significant exchange of individuals with other aggregations. We monitored main-stem Colorado River aggregations of humpback chub in Grand Canyon during 2010 to 2014 and compared our results to previous investigations. Relative abundance, as described by catch per unit effort (fish per hour) of adult humpback chub at most main-stem aggregations, generally increased from the 1990s to 2014. In addition, distribution of humpback chub in the main-stem Colorado River has increased since the 1990s. Movement of humpback chub between the Little Colorado River and other aggregations likely adds fish to those aggregations. There is clear evidence of reproduction near the 30-Mile aggregation, and reproduction at Middle Granite Gorge and downstream seems likely based on catches of gravid fish and captures of very young fish, especially during relatively warm water releases from Glen Canyon Dam, 2004 to 2011. Humpback chub relative abundance at Shinumo and Havasu Creek inflows increased following translocations of young humpback chub starting in 2009. In light of this information, we modify the original nine aggregations, combining two previously separate aggregations and dropping two locations to form six distinct aggregations of humpback chub. Trends in humpback chub abundance at main-stem aggregations, relative to management actions (for example, translocations) or changing environmental conditions (for example, river warming), informs management of the species across a riverscape scale within the Colorado River.

Colorado

Population expansion of Humpback chub in western Grand Canyon and hypothesized mechanisms

Humpback chub, Gila cypha , is an endangered warm water fish endemic to the Colorado River basin of southwestern North America. In Grand Canyon National Park, cold hypolimnetic water-release temperatures from Glen Canyon Dam have largely precluded successful spawning and recruitment of humpback chub in the mainstem Colorado River. Therefore, the species has utilized the warmer, more saline, and free-flowing Little Colorado River for its primary spawning habitat and continued existence. Based on long-term fish sampling efforts, we document local recruitment and population expansion of humpback chub in the western Grand Canyon and hypothesize that this is a result of recent warmer mainstem water temperatures. Continued recruitment and population expansion of humpback chub in the western Grand Canyon could potentially reduce extinction risk by providing population redundancy and less reliance upon the Little Colorado River for the species survival in the Grand Canyon.

Arizona

Nonlinear relationships can lead to bias in biomass calculations and drift-foraging models when using summaries of invertebrate drift data

Drift-foraging models offer a mechanistic description of how fish feed in flowing water and the application of drift-foraging bioenergetics models to answer both applied and theoretical questions in aquatic ecology is growing. These models typically include nonlinear descriptions of ecological processes and as a result may be sensitive to how model inputs are summarized because of a mathematical property of nonlinear equations known as Jensen’s inequality. In particular, we show that the way in which continuous size distributions of invertebrate prey are represented within foraging models can lead to biases within the modeling process. We begin by illustrating how different equations common to drift-foraging models are sensitive to invertebrate inputs. We then use two case studies to show how different representations of invertebrate prey can influence predictions of energy intake and lifetime growth. Greater emphasis should be placed on accurate characterizations of invertebrate drift, acknowledging that inferences from drift-foraging models may be influenced by how invertebrate prey are represented.

Environmental Biology of Fishes