Search USGSSearch

SEARCH · Search USGS

Results for “Journal of Applied Probability”

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

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

At least 73 records · Page 4Linked to original sources

Investigations of potential bias in the estimation of lambda using Pradel's (1996) model for capture-recapture data

Pradel's (1996) temporal symmetry model permitting direct estimation and modelling of population growth rate, u i , provides a potentially useful tool for the study of population dynamics using marked animals. Because of its recent publication date, the approach has not seen much use, and there have been virtually no investigations directed at robustness of the resulting estimators. Here we consider several potential sources of bias, all motivated by specific uses of this estimation approach. We consider sampling situations in which the study area expands with time and present an analytic expression for the bias in u i We next consider trap response in capture probabilities and heterogeneous capture probabilities and compute large-sample and simulation-based approximations of resulting bias in u i . These approximations indicate that trap response is an especially important assumption violation that can produce substantial bias. Finally, we consider losses on capture and emphasize the importance of selecting the estimator for u i that is appropriate to the question being addressed. For studies based on only sighting and resighting data, Pradel's (1996) u i ' is the appropriate estimator.

Journal of Applied Statistics

Challenges of DNA-based mark-recapture studies of American black bears

We explored whether genetic sampling would be feasible to provide a region-wide population estimate for American black bears ( Ursus americanus ) in the southern Appalachians, USA. Specifically, we determined whether adequate capture probabilities ( p >0.20) and population estimates with a low coefficient of variation (CV <20%) could be achieved given typical agency budget and personnel constraints. We extracted DNA from hair collected from baited barbed-wire enclosures sampled over a 10-week period on 2 study areas: a high-density black bear population in a portion of Great Smoky Mountains National Park and a lower density population on National Forest lands in North Carolina, South Carolina, and Georgia. We identified individual bears by their unique genotypes obtained from 9 microsatellite loci. We sampled 129 and 60 different bears in the National Park and National Forest study areas, respectively, and applied closed mark&ndash;recapture models to estimate population abundance. Capture probabilities and precision of the population estimates were acceptable only for sampling scenarios for which we pooled weekly sampling periods. We detected capture heterogeneity biases, probably because of inadequate spatial coverage by the hair-trapping grid. The logistical challenges of establishing and checking a sufficiently high density of hair traps make DNA-based estimates of black bears impractical for the southern Appalachian region. Alternatives are to estimate population size for smaller areas, estimate population growth rates or survival using mark&ndash;recapture methods, or use independent marking and recapturing techniques to reduce capture heterogeneity.

Georgia, North Carolina, South Carolina, Tennessee

Hydraulic anisotropy characterization of pneumatic-fractured sediments using azimuthal self potential gradient

The pneumatic fracturing technique is used to enhance the permeability and porosity of tight unconsolidated soils (e.g. clays), thereby improving the effectiveness of remediation treatments. Azimuthal self potential gradient (ASPG) surveys were performed on a compacted, unconsolidated clay block in order to evaluate their potential to delineate contaminant migration pathways in a mechanically-induced fracture network. Azimuthal resistivity (ARS) measurements were also made for comparative purposes. Following similar procedures to those used in the field, compressed kaolinite sediments were pneumatically fractured and the resulting fracture geometry characterized from strike analysis of visible fractures combined with strike data from optical borehole televiewer (BHTV) imaging. We subsequently injected a simulated treatment (electrolyte/dye) into the fractures. Both ASPG and ARS data exhibit anisotropic geoelectric signatures resulting from the fracturing. Self potentials observed during injection of electrolyte are consistent with electrokinetic theory and previous laboratory results on a fracture block model. Visual (polar plot) analysis and linear regression of cross plots show ASPG lobes are correlated with azimuths of high fracture strike density, evidence that the ASPG anisotropy is a proxy measure of hydraulic anisotropy created by the pneumatic fracturing. However, ARS data are uncorrelated with fracture strike maxima and resistivity anisotropy is probably dominated by enhanced surface conduction along azimuths of weak 'starter paths' formed from pulverization of the clay and increases in interfacial surface area. We find the magnitude of electrokinetic SP scales with the applied N2 gas pressure gradient (??PN2) for any particular hydraulically-active fracture set and that the positive lobe of the ASPG anomaly indicates the flow direction within the fracture network. These findings demonstrate the use of ASPG in characterizing the effectiveness of (1) pneumatic fracturing and (2) defining likely flow directions of remedial treatments in unconsolidated sediments and rock. ?? 2008 Elsevier B.V. All rights reserved.

Journal of Contaminant Hydrology

Falcon versus grouse: flight adaptations of a predator and its prey

Several falcons were trained to fly along a 500 m course to a lure. The air speeds of the more consistent performers averaged about 1.5 times their calculated minimum power speeds, and occasionally reached 2.1 times the minimum power speed. Wing beat frequencies of all the falcons were above those estimated from earlier field observations, and the same was true of wild Sage Grouse Centrocercus urophasianus, a regular falconer's quarry in the study area. Measurements of grouse killed by falcons showed that their wings were short, with broad slotted tips, whereas the falcons' wings were longer in relation to their body mass, and tapered. The short wings of grouse result in fast flight, high power requirements, and reduced capacity for aerobic flight. Calculations indicated that the grouse should fly faster than the falcons, and had the large amount of flight muscle needed to do so, but that the falcons would be capable of prolonged aerobic flight, whereas the grouse probably would not. We surmise that Sage Grouse cannot fly continuously without incurring an oxygen debt, and are therefore not long-distance migrants, although this limitation is partly due to their large size, and would not apply to smaller galliform birds such as ptarmigan Lagopus spp. The wing action seen in video recordings of the falcons was not consistent with the maintenance of constant circulation. We call it 'chase mode' because it appears to be associated with a high level of muscular exertion, without special regard to fuel economy. It shows features in common with the 'bounding' flight of passerines.

Journal of Avian Biology

Peclet number as affected by molecular diffusion controls transient anomalous transport in alluvial aquifer-aquitard complexes

This study evaluates the role of the Peclet number as affected by molecular diffusion in transient anomalous transport, which is one of the major knowledge gaps in anomalous transport, by combining Monte Carlo simulations and stochastic model analysis. Two alluvial settings containing either short- or long-connected hydrofacies are generated and used as media for flow and transport modeling. Numerical experiments show that 1) the Peclet number affects both the duration of the power-law segment of tracer breakthrough curves (BTCs) and the transition rate from anomalous to Fickian transport by determining the solute residence time for a given low-permeability layer, 2) mechanical dispersion has a limited contribution to the anomalous characteristics of late-time transport as compared to molecular diffusion due to an almost negligible velocity in floodplain deposits, and 3) the initial source dimensions only enhance the power-law tail of the BTCs at short travel distances. A tempered stable stochastic (TSS) model is then applied to analyze the modeled transport. Applications show that the time-nonlocal parameters in the TSS model relate to the Peclet number, P e . In particular, the truncation parameter in the TSS model increases nonlinearly with a decrease in P e due to the decrease of the mean residence time, and the capacity coefficient increases with an increase in molecular diffusion which is probably due to the increase in the number of immobile particles. The above numerical experiments and stochastic analysis therefore reveal that the Peclet number as affected by molecular diffusion controls transient anomalous transport in alluvial aquifer&ndash;aquitard complexes.

Journal of Contaminant Hydrology

Estimating detection and density of the Andean cat in the high Andes

The Andean cat ( Leopardus jacobita ) is one of the most endangered, yet least known, felids. Although the Andean cat is considered at risk of extinction, rigorous quantitative population studies are lacking. Because physical observations of the Andean cat are difficult to make in the wild, we used a camera-trapping array to photo-capture individuals. The survey was conducted in northwestern Argentina at an elevation of approximately 4,200 m during October-December 2006 and April-June 2007. In each year we deployed 22 pairs of camera traps, which were strategically placed. To estimate detection probability and density we applied models for spatial capture-recapture using a Bayesian framework. Estimated densities were 0.07 and 0.12 individual/km 2 for 2006 and 2007, respectively. Mean baseline detection probability was estimated at 0.07. By comparison, densities of the Pampas cat ( Leopardus colocolo ), another poorly known felid that shares its habitat with the Andean cat, were estimated at 0.74–0.79 individual/km 2 in the same study area for 2006 and 2007, and its detection probability was estimated at 0.02. Despite having greater detectability, the Andean cat is rarer in the study region than the Pampas cat. Properly accounting for the detection probability is important in making reliable estimates of density, a key parameter in conservation and management decisions for any species.

High Andes

Support of total maximum daily load programs using spatially referenced regression models

The spatially referenced regressions on watershed attributes modeling approach, as applied to predictions of total nitrogen flux in three North Carolina river basins, addresses several information needs identified by a National Research Council evaluation of the total maximum daily load program. The model provides reach-level predictions of the probability of exceeding water-quality criteria, and estimates of total nitrogen budgets. Model estimates of point- and diffuse-source contributions and nitrogen loss rates in streams and reservoirs compared moderately well with literature estimates. Maps of reach-level predictions of nutrient inputs and delivery provide an intuitive and spatially detailed summary of the origins and fate of nutrients within a basin.

Journal of Water Resources Planning and Management

Bayesian inversion of data from effusive volcanic eruptions using physics-based models: Application to Mount St. Helens 2004--2008

Physics-based models of volcanic eruptions can directly link magmatic processes with diverse, time-varying geophysical observations, and when used in an inverse procedure make it possible to bring all available information to bear on estimating properties of the volcanic system. We develop a technique for inverting geodetic, extrusive flux, and other types of data using a physics-based model of an effusive silicic volcanic eruption to estimate the geometry, pressure, depth, and volatile content of a magma chamber, and properties of the conduit linking the chamber to the surface. A Bayesian inverse formulation makes it possible to easily incorporate independent information into the inversion, such as petrologic estimates of melt water content, and yields probabilistic estimates for model parameters and other properties of the volcano. Probability distributions are sampled using a Markov-Chain Monte Carlo algorithm. We apply the technique using GPS and extrusion data from the 2004–2008 eruption of Mount St. Helens. In contrast to more traditional inversions such as those involving geodetic data alone in combination with kinematic forward models, this technique is able to provide constraint on properties of the magma, including its volatile content, and on the absolute volume and pressure of the magma chamber. Results suggest a large chamber of >40 km 3 with a centroid depth of 11–18 km and a dissolved water content at the top of the chamber of 2.6–4.9 wt%.

Washington

Detection of deer at remote camera sites in relation to snow conditions

In the rain-snow transition zone of the Pacific Northwest, climate change is expected to alter the incidence of rain-on-snow and freeze-thaw events, which will change snow density and hardness dynamics. In winter, the ability of economically and ecologically important wildlife species, such as deer ( Odocoileus spp.), to efficiently move through the landscape and access forage is mediated by snow conditions. Therefore, snow properties such as density and hardness can directly affect how energetically costly it is for these animals to survive. However, little is known about whether and how ungulates use habitats based on snow density and hardness. We deployed a stratified network of remote camera stations in complex forested terrain in Latah County, Idaho, USA, to remotely measure snow depth and detect deer. We also collected snow density and hardness measurements throughout the winter. We used these data to determine the degree to which the probability of deer presence at cameras could be explained by snow conditions and air temperature. Snow depth and density had negative relationships with the probability of deer presence, while ram resistance (a proxy for snow hardness) had a marginal positive effect. We were able to estimate snow conditions important to deer in winter 2020–2021 primarily using data obtained from cameras. This provides an important proof-of-concept that can be applied at different sites and climate conditions to gain a deeper understanding of how deer are affected by snowpack properties. These methods can be used by managers to determine how ungulates are affected by snow depth, density, and hardness collectively and subsequently inform ungulate management in a changing climate.

Idaho

Applying mark-resight, count, and telemetry data to estimate effective sampling area and fish density with stationary underwater cameras

Accurate estimates of abundance and density for geographically open populations must account for the effective sampling area (ESA) of sampling gears. We describe a Marked N-Mixture model to estimate ESA and density (number of individuals/unit area) from repeated counts of unmarked and marked individuals, integrating mark-resight, camera counts, and telemetry data of red snapper ( Lutjanus campechanus ) at a 1.6 km 2 reef off North Carolina, USA. Cameras recorded observations of unmarked and marked individuals, whereas telemetry data indicated the number of tagged fish present on the reef. We estimated density (95 individuals/km 2 , 95%CI: 58–149), ESA (which was lower when current direction was towards the camera), detection probability (0.06, 95%CI: 0.03–0.09), and covariate relationships. Simulation studies under different scenarios of data quality and space use identified positive bias in density estimates from N-mixture models due to fish movement. In contrast, the Marked N-Mixture model returned unbiased estimates of density, ESA, and detection parameters, and appears to be a more robust method for modeling density given the data available for this analysis. This approach can be applied to other populations where count and telemetry data overlap in space and time.

North Carolina

Biocrust and seed placement influence emergence rates of perennial grass Elymus elymoides across five North American deserts

1. Dryland vascular plant emergence is often limited by water availability especially with projected climate and precipitation changes. Biological soil crusts (biocrusts) can promote water retention and nutrient availability that benefit germination, and emergence yet can also act as a surface barrier preventing critical seed soil contact and hindering emergence. Alongside these factors, dryland fire frequency has increased with the inclusion of invasive species and vegetation structural changes. With enhanced continuous fine fuel distribution there is a high potential to disrupt biocrust-plant interactions and influence associated management actions. 2. This study explores the dynamic relationship between biocrusts and fire-related heating effects on seedling emergence across five North American deserts: the Chihuahuan Desert, Colorado Plateau, Great Basin, Mojave Desert and Sonoran Desert. We conducted a greenhouse-based seedling emergence experiment using Elymus elymoides (bottlebrush squirreltail), a common perennial grass, with biocrust and bare soil mesocosms in which half were heated to mimic the effects of wildfire temperature. 3. The variables that had the greatest influence on germination rate and germination timing were the presence of biocrust and seed placement (on top of vs within the biocrust/soil matrix). Emergence rate was greatest atop bare soil followed by seeds inserted into biocrust. Emergence timing was faster with biocrust present, but only when seeds were inserted into it. Desert origin of biocrust and soil collection also influenced germination where the probability of any one seed emerging was highest in the Chihuahuan and Mojave desert sites relative to the Sonoran desert site which showed the lowest probability. Heating had mixed effects whereby it increased the likelihood of emergence but did not affect the overall rate or timing. 4. Synthesis and applications . This study underscores the importance of healthy and impaired biocrusts in dryland systems and suggests ways in which the combination of biocrust and seed placement can influence plant establishment, in addition to providing insight into seeding strategies for managers and restoration practitioners working in dryland sites.

Arizona, California, Colorado, Nevada, New Mexico

Population dynamics of hispid cotton rats (Sigmodon hispidus) across a nitrogen-amended landscape

We conducted a mark-recapture experiment to examine the population dynamics of hispid cotton rats (Sigmodon hispidus) in response to low-level nitrogen amendments (16.4 kg nitrogen/ha per year) and exclosure fencing in an old-field grassland. The experimental design consisted of sixteen 0.16-ha plots with 4 replicates of each treatment combination. We predicted that densities, reproductive success, movement probabilities, and survival rates of cotton rats would be greater on nitrogen-amended plots because of greater aboveground biomass and canopy cover. Population densities of cotton rats tended to be highest on fenced nitrogen plots, but densities on unfenced nitrogen plots were similar to those on control and fenced plots. We observed no distinct patterns in survival rates, reproductive success, or movement probabilities with regard to nitrogen treatments. However, survival rates and reproductive success tended to be higher for cotton rats on fenced plots than for those on unfenced plots and this was likely attributable to decreased predation on fenced plots. As low-level nitrogen amendments continue to be applied, we predict that survival, reproduction, and population-growth rates of cotton rats on control plots, especially fenced plots with no nitrogen amendment, will eventually exceed those on nitrogen-amended plots as a result of higher plant-species diversity, greater food availability, and better quality cover.

Canadian Journal of Zoology

Relation of surflcial erosion on hillslopes to profile geometry

Differences in erosion rates in relation to the geometry of hillslope profiles were hypothesized by R. E. Horton in 1945. Experimental testing of this hypothesis is described for natural hillslopes in western Colorado using a rainfall simulator to apply five "storms" each having an intensity of 1.85 inches per hour (47 millimeters per hour). Erosion from straight hillslope segments was about twice the amount from convex and concave segments along nine hillslope transects. Standard deviations of erosion measurements for each segment show that the greatest variation occurs in concave slope segments. This is probably because the concave segment is the only place on the profile where deposition occurs. Analysis of variance indicates that the amount of erosion is significant at the 1-percent level of confidence between convex and straight slope segments and at the 5-percent level of confidence between straight and convex slope segments.

Colorado

Broad‐scale occurrence of a subsidized avian predator: reducing impacts of ravens on sage‐grouse and other sensitive prey

Expanding human enterprise across remote environments impacts numerous wildlife species. Anthropogenic resources provide subsidies for generalist predators that can lead to cascading effects on prey species at lower trophic levels. A fundamental challenge for applied ecologists is to disentangle natural and anthropogenic influences on species occurrence, and subsequently develop spatially explicit models to help inform management and conservation decisions. Using Bayesian hierarchical occupancy models, we mapped the broad‐scale occurrence of common ravens Corvus corax as a function of natural and anthropogenic landscape covariates using >15,000 point count surveys performed during 2007–2016 within the Great Basin region, USA. Raven abundance and distribution is substantially increasing across the American west due to unintended anthropogenic resource subsidies. Importantly, ravens prey on eggs and chicks of numerous species including greater sage‐grouse Centrocercus urophasianus , an indicator species whose decline is at the centre of national conservation strategies and land use policies.Anthropogenic factors that contributed to greater raven occurrence were: increased road density, presence of transmission lines, agricultural activity, and presence of roadside rest areas. Natural landscape characteristics included lower elevations with greener vegetation (NDVI), greater stream and habitat edge densities, and lower percentages of big sagebrush A. tridentata spp . Interactions between anthropogenic sources of nesting substrate and food subsidies suggested that raven occurrence increased multiplicatively when these resource subsidies co‐occurred. Overall, the average probability of raven occurrence estimated within sagebrush ecosystems of the study area was ~0.83. Synthesis and applications . We demonstrate how anthropogenic factors can be disentangled from natural effects when making spatially‐explicit predictions of subsidized predators occurring across expansive landscapes. This approach can guide management decisions where subsidized predators overlap sensitive prey habitats. For example, we identify areas where elevated raven occurrence coincides with breeding sage‐grouse concentration areas and appears to be largely driven by anthropogenic factors. Management applications could focus on reducing raven access to anthropogenic subsidies in these areas, while prioritizing habitat improvements for sage‐grouse elsewhere. Our approach is applicable to other species where widespread survey data are available.

Journal of Applied Ecology

Adaptive monitoring in action: Reconsidering design-based estimators reveals underestimation of whitebark pine disease prevalence in the Greater Yellowstone Ecosystem

Identifying and understanding status and trends in ecological indicators motivates continual monitoring over decades. Many programs rely on probability surveys and their companion design-based estimators for status assessments (e.g. Horvitz–Thompson). Design-based estimators do not easily extend to trend estimation nor situations with observation errors. Field-based monitoring efforts inevitably have turnover of field crew members which may affect consistency and accuracy of data collection over time. Additionally, design-based estimators ignore the complexities of spatial and temporal heterogeneity in an ecological indicator and how this variability may be linked to environmental or biological dynamics. We propose monitoring programs should re-evaluate their prescribed statistical methods, consider model-based approaches and adapt their sampling designs as needed to improve inferences. The Greater Yellowstone Ecosystem, home to two of the most iconic U.S. National Parks, has experienced significant declines in whitebark pine Pinus albicaulis communities due to forest pathogens, insect outbreaks, wildland fires and drought. Whitebark pine is a keystone species found in mountainous environments throughout the Western U.S. and Canada. We assessed the design-based ratio estimator originally recommended for estimating prevalence of white pine blister rust Cronartium ribicola . We compared the design-based estimator to a model-based approach that accounts for the sampling design, imperfect detection and allows for infection probabilities to vary over space and time. Ignoring observation errors led to lower estimated prevalence of white pine blister rust in the general population. Using model-based approaches, we found that the probability of infection has increased since 2004. However, overall prevalence likely has not changed because of the mountain pine beetle Dendroctonus ponderosae -induced shift towards smaller diameter trees that have a lower probability of infection compared to their larger cohorts. Synthesis and Applications . Using a design-based approach to detect change in ecological indicators falls short because of the inability to account for observation errors or to explore environmental or biological factors explaining temporal dynamics. Inherently understanding the mechanisms leading to changes in an ecological indicator over time informs potential management actions. Our assessment underscores the need for continued evaluation and updating of a monitoring program's sampling design and analytical procedures to maintain relevancy.

Wyoming

Sampling characteristics and calibration of snorkel counts to estimate stream fish populations

Snorkeling is a versatile technique for estimating lotic fish population characteristics; however, few investigators have evaluated its accuracy at population or assemblage levels. We evaluated the accuracy of snorkeling using prepositioned areal electrofishing (PAE) for estimating fish populations in a medium-sized Appalachian Mountain river during fall 2008 and summer 2009. Strip-transect snorkel counts were calibrated with PAE counts in identical locations among macrohabitats, fish species or taxa, and seasons. Mean snorkeling efficiency (i.e., the proportion of individuals counted from the true population) among all taxa and seasons was 14.7% (SE, 2.5%), and the highest efficiencies were for River Chub Nocomis micropogon at 21.1% (SE, 5.9%), Central Stoneroller Campostoma anomalum at 20.3% (SE, 9.6%), and darters (Percidae) at 17.1% (SE, 3.7%), whereas efficiencies were lower for shiners ( Notropis spp., Cyprinella spp., Luxilus spp.) at 8.2% (SE, 2.2%) and suckers (Catostomidae) at 6.6% (SE, 3.2%). Macrohabitat type, fish taxon, or sampling season did not significantly explain variance in snorkeling efficiency. Mean snorkeling detection probability (i.e., probability of detecting at least one individual of a taxon) among fish taxa and seasons was 58.4% (SE, 6.1%). We applied the efficiencies from our calibration study to adjust snorkel counts from an intensive snorkeling survey conducted in a nearby reach. Total fish density estimates from strip-transect counts adjusted for snorkeling efficiency were 7,288 fish/ha (SE, 1,564) during summer and 15,805 fish/ha (SE, 4,947) during fall. Precision of fish density estimates is influenced by variation in snorkeling efficiency and sample size and may be increased with additional sampling effort. These results demonstrate the sampling properties and utility of snorkeling to characterize lotic fish assemblages with acceptable efficiency and detection probability, less effort, and no mortality, compared with traditional sampling methods.

North American Journal of Fisheries Management

Integrating encounter theory with decision analysis to evaluate collision risk and determine optimal protection zones for wildlife

1.Better understanding human‐wildlife interactions and their links with management can help improve the design of wildlife protection zones. One example is the problem of wildlife collisions with vehicles or human‐built structures (e.g. power lines, wind farms). In fact, collisions between marine wildlife and watercraft are among the major threats faced by several endangered species of marine mammals. Natural resource managers are therefore interested in finding cost‐effective solutions to mitigate these threats. 2.We combined abundance estimators with encounter rate theory to estimate relative lethal collision risk of the Florida manatee ( Trichechus manatus latirostris ) from watercraft. We first modeled seasonal abundance of watercraft and manatees using a Bayesian analysis of aerial survey count data. We then modeled relative lethal collision risk in space and across seasons. Finally, we applied decision analysis and Linear Integer Programming to determine the optimal design of speed zones in terms of relative risk to manatees and costs to waterway users. We used a Pareto efficient frontier approach to evaluate the performance of alternative zones, which included additional practical considerations (e.g. spatial aggregation of speed zones) in relation to the optimal zone configurations. 3.Under the various relationships for probability of death given strike speed that we considered, the current speed zones reduced the relative lethal collision risk by an average of 51.5% to 70% compared to the scenario in which all speed regulations were removed (i.e. the no‐protection scenario). We identified optimal zones and near‐optimal zones with additional management considerations that improved upon the current zones in terms of cost or relative risk. 4.Policy Implications : Our analytical framework combines encounter rate theory and decision analysis to quantify the effectiveness of speed zones protecting manatees while accounting for uncertainty. Our approach can be used to optimize the design of protection zones intended to reduce conflicts between human waterborne activity and marine mammals. This framework could be extended to address many other problems of human‐wildlife interactions, such as the optimal placement of wind farms to minimize collisions with wildlife or the optimal allocation of ranger effort to mitigate poaching threats.

Journal of Applied Ecology

Recent changes in burbot growth in Lake Erie

Recruitment of burbot Lota lota in eastern Lake Erie, estimated by catches of age-4 burbot, was high during 1997–2001 and then abruptly declined to low levels during 2002–2007. The invasive round goby Neogobius melanostomus , a benthic species, was first collected in trawl assessments in eastern Lake Erie in 1999, and was first found in stomachs of burbot in 2001. By 2003, round goby became an important prey in the diet of burbot. We hypothesized that the combined effects of low recruitment and consumption of round goby would result in increased size-at-age in burbot. We reasoned that: (i) decreased competition for resources among juveniles should result in larger adults, and (ii) consumption of a benthic prey by a bottom-dwelling predator such as burbot should require less foraging in the water column, and thus less energetic expenditure. We divided our data into two temporal periods: one in which burbot belonged to strong year classes and ate few, if any round goby (i.e., year classes 1989–1997 collected during 1997–2001) and one in which burbot belonged to weak year classes and probably ate round gobies by age 4 (year classes 1998–2003 collected during 2002–2007). Mass and total lengths at ages 4–7 were generally higher during the second period. However, the rates of growth between ages 4 and 7 were not different for the two periods. The results indicate that greater growth at ages 0–4 resulted in larger size at ages 4–7 in the latter period. More information on juvenile diet and growth in burbot is needed for effective conservation of burbot stocks.

Journal of Applied Ichthyology