Abundance indices for long-lived tropical snappers: estimating standardized catch rates from spatially and temporally coarse logbook data

The Northern Demersal Scalefish Fishery has historically comprised a small fleet (≤10 vessels year−1) operating over a relatively large area off the northwest coast of Australia. This multispecies fishery primarily harvests two species of snapper: goldband snapper, Pristipomoides multidens and red emperor, Lutjanus sebae. A key input to age-structured assessments of these stocks has been the annual time-series of the catch rate. We used an approach that combined Generalized Linear Models, spatio-temporal imputation, and computer-intensive methods to standardize the fishery catch rates and report uncertainty in the indices. These analyses, which represent one of the first attempts to standardize fish trap catch rates, were also augmented to gain additional insights into the effects of targeting, historical effort creep, and spatio-temporal resolution of catch and effort data on trap fishery dynamics. Results from monthly reported catches (i.e. 1993 on) were compared with those reported daily from more recently (i.e. 2008 on) enhanced catch and effort logbooks. Model effects of catches of one species on the catch rates of another became more conspicuous when the daily data were analysed and produced estimates with greater precision. The rate of putative effort creep estimated for standardized catch rates was much lower than estimated for nominal catch rates. These results therefore demonstrate how important additional insights into fishery and fish population dynamics can be elucidated from such “pre-assessment” analyses.

A regionally informed abundance index for supporting integrative analyses across butterfly monitoring schemes

1. The rapid expansion of systematic monitoring schemes necessitates robust methods to reliably assess species' status and trends. Insect monitoring poses a challenge where there are strong seasonal patterns, requiring repeated counts to reliably assess abundance. Butterfly monitoring schemes (BMSs) operate in an increasing number of countries with broadly the same methodology, yet they differ in their observation frequency and in the methods used to compute annual abundance indices. 2. Using simulated and observed data, we performed an extensive comparison of two approaches used to derive abundance indices from count data collected via BMS, under a range of sampling frequencies. Linear interpolation is most commonly used to estimate abundance indices from seasonal count series. A second method, hereafter the regional generalized additive model (GAM), fits a GAM to repeated counts within sites across a climatic region. For the two methods, we estimated bias in abundance indices and the statistical power for detecting trends, given different proportions of missing counts. We also compared the accuracy of trend estimates using systematically degraded observed counts of the Gatekeeper Pyronia tithonus (Linnaeus 1767). 3. The regional GAM method generally outperforms the linear interpolation method. When the proportion of missing counts increased beyond 50%, indices derived via the linear interpolation method showed substantially higher estimation error as well as clear biases, in comparison to the regional GAM method. The regional GAM method also showed higher power to detect trends when the proportion of missing counts was substantial. 4. Synthesis and applications. Monitoring offers invaluable data to support conservation policy and management, but requires robust analysis approaches and guidance for new and expanding schemes. Based on our findings, we recommend the regional generalized additive model approach when conducting integrative analyses across schemes, or when analysing scheme data with reduced sampling efforts. This method enables existing schemes to be expanded or new schemes to be developed with reduced within-year sampling frequency, as well as affording options to adapt protocols to more efficiently assess species status and trends across large geographical scales.

Monitoring Butterfly Abundance: Beyond Pollard Walks

Most butterfly monitoring protocols rely on counts along transects (Pollard walks) to generate species abundance indices and track population trends. It is still too often ignored that a population count results from two processes: the biological process (true abundance) and the statistical process (our ability to properly quantify abundance). Because individual detectability tends to vary in space (e.g., among sites) and time (e.g., among years), it remains unclear whether index counts truly reflect population sizes and trends. This study compares capture-mark-recapture (absolute abundance) and count-index (relative abundance) monitoring methods in three species (Maculinea nausithous and Iolana iolas: Lycaenidae; Minois dryas: Satyridae) in contrasted habitat types. We demonstrate that intraspecific variability in individual detectability under standard monitoring conditions is probably the rule rather than the exception, which questions the reliability of count-based indices to estimate and compare specific population abundance. Our results suggest that the accuracy of count-based methods depends heavily on the ecology and behavior of the target species, as well as on the type of habitat in which surveys take place. Monitoring programs designed to assess the abundance and trends in butterfly populations should incorporate a measure of detectability. We discuss the relative advantages and inconveniences of current monitoring methods and analytical approaches with respect to the characteristics of the species under scrutiny and resources availability.

The North Sea autumn spawning Herring (Clupea harengus L.) Spawning Component Abundance Index (SCAI)

The North Sea autumn-spawning herring (Clupea harengus) stock consists of a set of different spawning components. The dynamics of the entire stock have been well characterized, but although time-series of larval abundance indices are available for the individual components, study of the dynamics at the component level has historically been hampered by missing observations and high sampling noise. A simple state-space statistical model is developed that is robust to these problems, gives a good fit to the data, and proves capable of both handling and predicting missing observations well. Furthermore, the sum of the fitted abundance indices across all components proves an excellent proxy for the biomass of the total stock, even though the model utilizes information at the individual-component level. The Orkney-Shetland component appears to have recovered faster from historic depletion events than the other components, whereas the Downs component has been the slowest. These differences give rise to changes in stock composition, which are shown to vary widely within a relatively short time. The modelling framework provides a valuable tool for studying and monitoring the dynamics of the individual components of the North Sea herring stock.

The distribution and size composition of finfish, American lobster, and long-finned squid in Long Island Sound based on the Connecticut Fisheries Division Bottom Trawl Survey, 1984–1994

The distribution, abundance, and length composition of marine finfish, lobster, and squid in Long Island Sound were examined relative to season and physical features of the Sound, using Connecticut Department of Environmental Protection trawl survey data collected from 1984 to 1994. The following are presented: seasonal distribution maps for 59 species, abundance indices for 41 species, and length frequencies for 26 species. In addition, a broader view of habitat utilization in the Sound was examined by mapping aggregated catches (total catch per tow, demersal catch per tow, and pelagic catch per tow) and by comparing species richness and mean aggregate catch/tow by analysis of variance (ANOVA) among eight habitat types defined by depth interval and bottom type. For many individual species, seasonal migration patterns and preference for particular areas within Long Island Sound were evident. The aggregate distribution maps show that overall abundance was lower in the eastern Sound than the central and western portions. Demersal and pelagic temporal abundance show opposite trends—demersals were abundant in spring and declined through summer and fall, whereas pelagic abundance was low in spring and increased into fall. The analysis of habitat types revealed significant differences for both species richness and mean catch per tow. Generally, species richness was highest in habitats within the central area of the Sound and lowest in eastern habitats. The aggregate mean catch was highest in the western and central habitats, and declined eastward. (PDF file contains 199 pages.)

Historie und aktueller Stand des Internationalen Heringslarvensurveys in der Nordsee

Since years the International Herring Larvae Survey Program (IHLS) is an important and internationally established survey program in the North Sea. The IHLS serves the calibration of stock abundance estimates based on information from the commercial fishery and the method of Integrated Catch Analysis (ICA) which is a specific derivate of the Virtual Population Analysis (VPA). Meanwhile the IHLS database has been transferred from Aberdeen to Kiel and it has been agreed that the Institut für Meereskunde Kiel should continue to maintain this database and provide the abundance indices to be utilized by the ICES Herring Assessment Working Group as one of the means for assessing the state of the herring stock in the North Sea. For establishing the calculation procedure at Kiel, it was necessary to optimize both, the survey design and the index calculation. This article gives an overview over the survey’s history, it’s geography, the sampling design, the information content of the IHLS data base and the various methods of calculating the different indices necessary for the calibration.

Catches per unit of effort of bigeye tuna: a new analysis with regression trees and simulated annealing

ENGLISH: We analyzed catches per unit of effort (CPUE) from the Japanese longline fishery for bigeye tuna (Thunnus obesus) in the central and eastern Pacific Ocean (EPO) with regression tree methods. Regression trees have not previously been used to estimate time series of abundance indices fronl CPUE data. The "optimally sized" tree had 139 parameters; year, month, latitude, and longitude interacted to affect bigeye CPUE. The trend in tree-based abundance indices for the EPO was similar to trends estimated from a generalized linear model and fronl an empirical model that combines oceanographic data with information on the distribution of fish relative to environmental conditions. The regression tree was more parsimonious and would be easier to implement than the other two nl0dels, but the tree provided no information about the nlechanisms that caused bigeye CPUEs to vary in time and space. Bigeye CPUEs increased sharply during the mid-1980's and were more variable at the northern and southern edges of the fishing grounds. Both of these results can be explained by changes in actual abundance and changes in catchability. Results from a regression tree that was fitted to a subset of the data indicated that, in the EPO, bigeye are about equally catchable with regular and deep longlines. This is not consistent with observations that bigeye are more abundant at depth and indicates that classification by gear type (regular or deep longline) may not provide a good measure of capture depth. Asimulated annealing algorithm was used to summarize the tree-based results by partitioning the fishing grounds into regions where trends in bigeye CPUE were similar. Simulated annealing can be useful for designing spatial strata in future sampling programs. SPANISH: Analizamos la captura por unidad de esfuerzo (CPUE) de la pesquería palangrera japonesa de atún patudo (Thunnus obesus) en el Océano Pacifico oriental (OPO) y central con métodos de árbol de regresión. Hasta ahora no se han usado árboles de regresión para estimar series de tiempo de índices de abundancia a partir de datos de CPUE. EI árbol de "tamaño optimo" tuvo 139 parámetros; ano, mes, latitud, y longitud interactuaron para afectar la CPUE de patudo. La tendencia en los índices de abundancia basados en árboles para el OPO fue similar a las tendencias estimadas con un modelo lineal generalizado y con un modelo empírico que combina datos oceanográficos con información sobre la distribución de los peces en relación con las condiciones ambientales. EI árbol de regresión fue mas parsimonioso y seria mas fácil de utilizar que los dos otros modelos, pero no proporciono información sobre los mecanismos que causaron que las CPUE de patudo valiaran en el tiempo y en el espacio. Las CPUE de patudo aumentaron notablemente a mediados de los anos 80 y fueron mas variables en los extremos norte y sur de la zona de pesca. Estos dos resultados pueden ser explicados por cambios en la abundancia real y cambios en la capturabilidad. Los resultados de un arbal de regresión ajustado a un subconjunto de los datos indican que, en el OPO, el patudo es igualmente capturable con palangres regulares y profundos. Esto no es consistente con observaciones de que el patudo abunda mas a profundidad e indica que clasificación por tipo de arte (palangre regular 0 profundo) podría no ser una buena medida de la profundidad de captura. Se uso un algoritmo de templado simulado para resumir los resultados basados en el árbol clasificando las zonas de pesca en zonas con tendencias similares en la CPUE de patudo. El templado simulado podría ser útil para diseñar estratos espaciales en programas futuros de muestreo. (PDF contains 45 pages.)

Using experiments and expert judgment to model catchability of Pacific rockfishes in trawl surveys, with application to bocaccio (Sebastes paucispinis) off British Columbia

The time series of abundance indices for many groundfish populations, as determined from trawl surveys, are often imprecise and short, causing stock assessment estimates of abundance to be imprecise. To improve precision, prior probability distributions (priors) have been developed for parameters in stock assessment models by using meta-analysis, expert judgment on catchability, and empirically based modeling. This article presents a synthetic approach for formulating priors for rockfish trawl survey catchability (qgross). A multivariate prior for qgross for different surveys is formulated by using 1) a correction factor for bias in estimating fish density between trawlable and untrawlable areas, 2) expert judgment on trawl net catchability, 3) observations from trawl survey experiments, and 4) data on the fraction of population biomass in each of the areas surveyed. The method is illustrated by using bocaccio (Sebastes paucipinis) in British Columbia. Results indicate that expert judgment can be updated markedly by observing the catch-rate ratio from different trawl gears in the same areas. The marginal priors for qgross are consistent with empirical estimates obtained by fitting a stock assessment model to the survey data under a noninformative prior for qgross. Despite high prior uncertainty (prior coefficients of variation ≥0.8) and high prior correlation between qgross, the prior for qgross still enhances the precision of key stock assessment quantities.

Removing observational noise from fisheries-independent time series data using ARIMA models

Abundance indices derived from fishery-independent surveys typically exhibit much higher interannual variability than is consistent with the within-survey variance or the life history of a species. This extra variability is essentially observation noise (i.e. measurement error); it probably reflects environmentally driven factors that affect catchability over time. Unfortunately, high observation noise reduces the ability to detect important changes in the underlying population abundance. In our study, a noise-reduction technique for uncorrelated observation noise that is based on autoregressive integrated moving average (ARIMA) time series modeling is investigated. The approach is applied to 18 time series of finfish abundance, which were derived from trawl survey data from the U.S. northeast continental shelf. Although the a priori assumption of a random-walk-plus-uncorrelated-noise model generally yielded a smoothed result that is pleasing to the eye, we recommend that the most appropriate ARIMA model be identified for the observed time series if the smoothed time series will be used for further analysis of the population dynamics of a species.

Age composition and distribution of red drum (Sciaenops ocellatus) in offshore waters of the north central Gulf of Mexico: an evaluation of a stock under a federal harvest moratorium

Because of a lack of fishery-dependent data, assessment of the recovery of fish stocks that undergo the most aggressive form of management, namely harvest moratoriums, remains a challenge. Large schools of red drum (Sciaenops ocellatus) were common along the northern Gulf of Mexico until the late 1980s when increased fishing effort quickly depleted the stock. After 24 years of harvest moratorium on red drum in federal waters, the stock is in need of reassessment; however, fisherydependent data are not available in federal waters and fishery-independent data are limited. We document the distribution, age composition, growth, and condition of red drum in coastal waters of the north central Gulf of Mexico, using data collected from a nearshore, randomized, bottom longline survey. Age composition of the fishery-independent catch indicates low mortality of fish age 6 and above and confirms the effectiveness of the federal fishing moratorium. Bottom longline surveys may be a cost-effective method for developing fishery-independent indices for red drum provided additional effort can be added to nearshore waters (<20 m depth). As with most stocks under harvest bans, effective monitoring of the recovery of red drum will require the development of fishery-independent indices. With limited economic incentive to evaluate non-exploited stocks, the most cost-effective approach to developing such monitoring is expansion of existing fishery independent surveys. We examine this possibility for red drum in the Gulf of Mexico and recommend the bottom longline survey conducted by the National Marine Fisheries Service expand effort in nearshore areas to allow for the development of long-term abundance indices for red drum.

Biology and population dynamics of cowcod (Sebastes levis) in the southern California Bight

Cowcod (Sebastes levis) is a large (100-cm-FL), long-lived (maximum observed age 55 yr) demersal rockfish taken in multispecies commercial and recreational fisheries off southern and central California. It lives at 20–500 m depth: adults (>44 cm TL) inhabit rocky areas at 90–300 m and juveniles inhabit fine sand and clay at 40–100 m. Both sexes have similar growth and maturity. Both sexes recruit to the fishery before reaching full maturity. Based on age and growth data, the natural mortality rate is about M =0.055/yr, but the estimate is uncertain. Biomass, recruitment, and mortality during 1951–98 were estimated in a delay-difference model with catch data and abundance indices. The same model gave less precise estimates for 1916–50 based on catch data and assumptions about virgin biomass and recruitment such as used in stock reduction analysis. Abundance indices, based on rare event data, included a habitat-area–weighted index of recreational catch per unit of fishing effort (CPUE index values were 0.003–0.07 fish per angler hour), a standardized index of proportion of positive tows in CalCOFI ichthyoplankton survey data (binomial errors, 0–13% positive tows/yr), and proportion of positive tows for juveniles in bottom trawl surveys (binomial errors, 0–30% positive tows/yr). Cowcod are overfished in the southern California Bight; biomass during the 1998 season was about 7% of the virgin level and recent catches have been near 20 metric tons (t)/yr. Projections based on recent recruitment levels indicate that biomass will decline at catch levels > 5 t/yr. Trend data indicate that recruitment will be poor in the near future. Recreational fishing effort in deep water has increased and has become more effective for catching cowcod. Areas with relatively high catch rates for cowcod are fewer and are farther offshore. Cowcod die after capture and cannot be released alive. Two areas recently closed to bottom fishing will help rebuild the cowcod stock.

Camarāo de profundidade: relatório dos cruzeiros realizados com o n/i "Ernst Haeckel" em Agosto-Setembro 1980, Novembro-Dezembro 1980, Janeiro-Fevereiro 1981, Janeiro-Fevereiro 1982

Reports of cruises carried out off Mozambique between parallels 21 degree 00 and 26 degree 20'S and at depths of 400 and 800 meters including abundance indices, catch composition, geographic and bathymetric distributions of shrimp catches are presented. The main species found were Hymenopenaeus triarthrus and Aristeomorpha foliacea which represent 64 to 72% of total catches.

Effects of predator control on behaviour of an apex predator and indirect consequences for mesopredator suppression

1. Apex predators can benefit ecosystems through top–down control of mesopredators and herbivores. However, apex predators are often subject to lethal control aimed at minimizing attacks on livestock. Lethal control can affect both the abundance and behaviour of apex predators. These changes could in turn influence the abundance and behaviour of mesopredators.

2. We used remote camera surveys at nine pairs of large Australian rangeland properties, comparing properties that controlled dingoes Canis lupus dingo with properties that did not, to test the effects of predator control on dingo activity and to evaluate the responses of a mesopredator, the feral cat Felis catus.

3. Indices of dingo abundance were generally reduced on properties that practiced dingo control, in comparison with paired properties that did not, although the effect size of control was variable. Dingoes in uncontrolled populations were crepuscular, similar to major prey. In populations subject to control, dingoes became less active around dusk, and activity was concentrated in the period shortly before dawn.

4. Shifts in feral cat abundance indices between properties with and without dingo control were inversely related to corresponding shifts in indices of dingo abundance. There was also a negative relationship between predator visitation rates at individual camera stations, suggesting cats avoided areas where dingoes were locally common. Reduced activity by dingoes at dusk was associated with higher activity of cats at dusk.

5. Our results suggest that effective dingo control not only leads to higher abundance of feral cats, but allows them to optimize hunting behaviour when dingoes are less active. This double effect could amplify the impacts of dingo control on prey species selected by cats. In areas managed for conservation, stable dingo populations may thus contribute to management objectives by restricting feral cat access to prey populations.

6. Synthesis and applications. Predator control not only reduces indices of apex predator abundance but can also modify their behaviour. Hence, indicators other than abundance, such as behavioural patterns, should be considered when estimating a predator's capacity to effectively interact with lower trophic guilds. Changes to apex predator behaviour may relax limitations on the behaviour of mesopredators, providing enhanced access to resources and prey.

Dingoes can help conserve wildlife and our methods can tell

Management of apex predators is among the most controversial wildlife management issues globally. In Australia, some ecologists have advocated using the dingo, Canis dingo, as a tool for conservation management, due to evidence that they suppress invasive mesopredators. Hayward & Marlow (Journal of Applied Ecology, 51, 2014 and 835) questioned the capacity of dingoes to provide benefits to native biodiversity due to their inability to eradicate foxes and cats. They also argued that indices of abundance commonly used in studies of mesopredator release by dingoes (namely, track-based indices) invalidate the conclusions of the studies. Hayward & Marlow caution conservation practitioners against incorporating dingoes into conservation programmes. Counter to their claims, we summarise research showing that the suppression of invasive mesopredators (cf. eradication) can enhance populations of native species and is therefore a meaningful conservation objective. We highlight literature supporting the hypothesis that dingoes suppress mesopredator abundance and activity, which in turn benefits native biodiversity. We show that Hayward & Marlow overlook many studies of carnivores that show track indices capture a large amount of the variation in the density of medium- and large-sized carnivores. Synthesis and applications. Practitioners cannot afford to wait to act given the perilous state of Australia's mammal species, and we argue that the evidence is sufficiently strong to justify managing dingoes for biodiversity conservation.

End-to-end modelling for an ecosystem approach to fisheries in the Northern Humboldt Current Ecosystem

This work represents an original contribution to the methodology for ecosystem models' development as well as the rst attempt of an end-to-end (E2E) model of the Northern Humboldt Current Ecosystem (NHCE). The main purpose of the developed model is to build a tool for ecosystem-based management and decision making, reason why the credibility of the model is essential, and this can be assessed through confrontation to data. Additionally, the NHCE exhibits a high climatic and oceanographic variability at several scales, the major source of interannual variability being the interruption of the upwelling seasonality by the El Niño Southern Oscillation, which has direct e ects on larval survival and sh recruitment success. Fishing activity can also be highly variable, depending on the abundance and accessibility of the main shery resources. This context brings the two main methodological questions addressed in this thesis, through the development of an end-to-end model coupling the high trophic level model OSMOSE to the hydrodynamics and biogeochemical model ROMS-PISCES: i) how to calibrate ecosystem models using time series data and ii) how to incorporate the impact of the interannual variability of the environment and shing. First, this thesis highlights some issues related to the confrontation of complex ecosystem models to data and proposes a methodology for a sequential multi-phases calibration of ecosystem models. We propose two criteria to classify the parameters of a model: the model dependency and the time variability of the parameters. Then, these criteria along with the availability of approximate initial estimates are used as decision rules to determine which parameters need to be estimated, and their precedence order in the sequential calibration process. Additionally, a new Evolutionary Algorithm designed for the calibration of stochastic models (e.g Individual Based Model) and optimized for maximum likelihood estimation has been developed and applied to the calibration of the OSMOSE model to time series data. The environmental variability is explicit in the model: the ROMS-PISCES model forces the OSMOSE model and drives potential bottom-up e ects up the foodweb through plankton and sh trophic interactions, as well as through changes in the spatial distribution of sh. The latter e ect was taken into account using presence/ absence species distribution models which are traditionally assessed through a confusion matrix and the statistical metrics associated to it. However, when considering the prediction of the habitat against time, the variability in the spatial distribution of the habitat can be summarized and validated using the emerging patterns from the shape of the spatial distributions. We modeled the potential habitat of the main species of the Humboldt Current Ecosystem using several sources of information ( sheries, scienti c surveys and satellite monitoring of vessels) jointly with environmental data from remote sensing and in situ observations, from 1992 to 2008. The potential habitat was predicted over the study period with monthly resolution, and the model was validated using quantitative and qualitative information of the system using a pattern oriented approach. The nal ROMS-PISCES-OSMOSE E2E ecosystem model for the NHCE was calibrated using our evolutionary algorithm and a likelihood approach to t monthly time series data of landings, abundance indices and catch at length distributions from 1992 to 2008. To conclude, some potential applications of the model for shery management are presented and their limitations and perspectives discussed.