A general framework for integrating environmental time series into stock assessment models: model description, simulation testing, and example

We present a method to integrate environmental time series into stock assessment models and to test the significance of correlations between population processes and the environmental time series. Parameters that relate the environmental time series to population processes are included in the stock assessment model, and likelihood ratio tests are used to determine if the parameters improve the fit to the data significantly. Two approaches are considered to integrate the environmental relationship. In the environmental model, the population dynamics process (e.g. recruitment) is proportional to the environmental variable, whereas in the environmental model with process error it is proportional to the environmental variable, but the model allows an additional temporal variation (process error) constrained by a log-normal distribution. The methods are tested by using simulation analysis and compared to the traditional method of correlating model estimates with environmental variables outside the estimation procedure. In the traditional method, the estimates of recruitment were provided by a model that allowed the recruitment only to have a temporal variation constrained by a log-normal distribution. We illustrate the methods by applying them to test the statistical significance of the correlation between sea-surface temperature (SST) and recruitment to the snapper (Pagrus auratus) stock in the Hauraki Gulf–Bay of Plenty, New Zealand. Simulation analyses indicated that the integrated approach with additional process error is superior to the traditional method of correlating model estimates with environmental variables outside the estimation procedure. The results suggest that, for the snapper stock, recruitment is positively correlated with SST at the time of spawning

Use of state-space population dynamics models in hypothesis testing: advantages over simple log-linear regressions for modeling survival, illustrated with application to longfin smelt (Spirinchus thaleichthys)

AbstractFactors impacting the survival of individuals between two life stages have traditionally been evaluated using log-linear regression of the ratio of abundance estimates for the two stages. These analyses require simplifying assumptions that may impact the results of hypothesis tests and subsequent conclusions about the factors impacting survival. Modern statistical methods can reduce the dependence of analyses on these simplifying assumptions. State-space models and the related concept of random effects allow the modeling of both process and observation error. Nonlinear models and associated estimation techniques allow for flexibility in the system model, including density dependence, and in error structure. Population dynamics models link information from one stage to the next and over multiple time periods and automatically accommodate missing observations. We investigate the impact of observation error, density dependence, population dynamics, and data for multiple stages on hypothesis testing using data for longfin smelt in the San Francisco Bay-Delta

Using stock assessment information to assess fishing capacity of tuna fisheries

Abstract Arrizabalaga, H., Restrepo, V. R., Maunder, M. N., and Majkowski, J. 2009. Using stock assessment information to assess fishing capacity of tuna fisheries. – ICES Journal of Marine Science, 66: 1959–1966. In tuna and tuna-like fisheries, there is a need for periodic assessments of fishing capacity to aid management. However, the nature and quantity of data needed to apply conventional methodologies for estimating fishing capacity are not usually available for tuna fisheries. We discuss simple alternative approaches to estimate fishing capacity and related quantities (i.e. capacity utilization, excess capacity, and overcapacity) directly from stock assessment inputs and outputs that are usually available for most tuna (and many other) stocks. Sensitivity analyses are performed to assess the effect of different levels of data aggregation and different assumptions made during the stock assessments on estimates of fishing capacity. Main advantages and disadvantages of the proposed methodologies are also illustrated using stock assessment information from different tuna stocks with different historical developments and trends in fishing mortality

Lessons to be learned by comparing integrated fisheries stock assessment models (SAMs) with integrated population models (IPMs)

AEP was partially funded by the Cooperative Institute for Climate, Ocean, & Ecosystem Studies (CICOES) under NOAA Cooperative Agreement NA15OAR4320063, Contribution No. 2023-1331.Integrated fisheries stock assessment models (SAMs) and integrated population models (IPMs) are used in biological and ecological systems to estimate abundance and demographic rates. The approaches are fundamentally very similar, but historically have been considered as separate endeavors, resulting in a loss of shared vision, practice and progress. We review the two approaches to identify similarities and differences, with a view to identifying key lessons that would benefit more generally the overarching topic of population ecology. We present a case study for each of SAM (snapper from the west coast of New Zealand) and IPM (woodchat shrikes from Germany) to highlight differences and similarities. The key differences between SAMs and IPMs appear to be the objectives and parameter estimates required to meet these objectives, the size and spatial scale of the populations, and the differing availability of various types of data. In addition, up to now, typical SAMs have been applied in aquatic habitats, while most IPMs stem from terrestrial habitats. SAMs generally aim to assess the level of sustainable exploitation of fish populations, so absolute abundance or biomass must be estimated, although some estimate only relative trends. Relative abundance is often sufficient to understand population dynamics and inform conservation actions, which is the main objective of IPMs. IPMs are often applied to small populations of conservation concern, where demographic uncertainty can be important, which is more conveniently implemented using Bayesian approaches. IPMs are typically applied at small to moderate spatial scales (1 to 104 km2), with the possibility of collecting detailed longitudinal individual data, whereas SAMs are typically applied to large, economically valuable fish stocks at very large spatial scales (104 to 106 km2) with limited possibility of collecting detailed individual data. There is a sense in which a SAM is more data- (or information-) hungry than an IPM because of its goal to estimate absolute biomass or abundance, and data at the individual level to inform demographic rates are more difficult to obtain in the (often marine) systems where most SAMs are applied. SAMs therefore require more 'tuning' or assumptions than IPMs, where the 'data speak for themselves', and consequently techniques such as data weighting and model evaluation are more nuanced for SAMs than for IPMs. SAMs would benefit from being fit to more disaggregated data to quantify spatial and individual variation and allow richer inference on demographic processes. IPMs would benefit from more attempts to estimate absolute abundance, for example by using unconditional models for capture-recapture data.Publisher PDFPeer reviewe

Can diagnostic tests help identify model misspecification in integrated stock assessments?

AbstractA variety of data types can be included in contemporary integrated stock assessments to simultaneously provide information on all estimated parameters. Conflicts between data, which are often a symptom of model misspecification and evident as model misfit, can affect the estimates of important parameters and derived quantities. Unfortunately, there are few standard diagnostic tools available for integrated stock assessment models that can provide the analyst with all the information needed to determine if there is substantial model misspecification. In this study, we use simulation methods to evaluate the ability of commonly-used and recently-proposed diagnostic tests to detect model misspecification in the observation model process (i.e., the incorrect form for survey selectivity), systems dynamics (i.e., incorrect assumed values for steepness of the stock-recruitment relationship and natural mortality), and incorrect data weighting. The diagnostic tests evaluated here were: i) residuals analysis (SDNR and runs test); ii) retrospective analysis; iii) the R0 likelihood component profile; iv) the age-structured production model (ASPM); and v) catch-curve analysis (CCA). The efficacy of the diagnostic tests depended on whether the misspecification was in the observation or systems dynamics model. Residual analyses were easily the best detector of misspecification of the observation model while the ASPM test was the only good diagnostic for detecting misspecification of system dynamics model. Retrospective analysis and the R0 likelihood component profile infrequently detected misspecified models, and CCA had a high probability of rejecting correctly-specified models. Finally, applying multiple carefully selected diagnostics can increase the power to detect misspecification without substantially increasing the probability of falsely concluding there is misspecification when the model is correctly specified

Asymmetric Externalities of the Tuna Longline and Tuna Purse-Seine Fisheries in the Eastern Pacific Ocean

Abstract only

A cookbook for using model diagnostics in integrated stock assessments

n/

Review of potential line-transect methodologies for estimating abundance of dolphin stocks in the eastern tropical Pacific

A twelve-year hiatus in fishery-independent marine mammal surveys in the eastern tropical Pacific Ocean (ETP), combined with a mandate to monitor dolphin stock status under international agreements and the need for reliable stock status information to set dolphin bycatch limits in the tuna purse-seine fishery, has renewed debate about how best to assess and monitor ETP dolphin stock status. The high cost of replicating previous ship-based surveys has intensified this debate. In this review, transect methods for estimating animal abundance from dedicated research surveys are considered, with a focus on both contemporary and potential methods suitable for surveying large areas for dolphin species that can form large, multi-species aggregations. Covered in this review are potential improvements to the previous ship-based survey methodology, other ship-based methods, alternative approaches based on high-resolution imagery and passive acoustics, and combinations of ship-based and alternative approaches. It is concluded that for immediate management needs, ship-based surveys, with some suggested modifications to improve precision, are the only reliable option despite their high cost. However, it is recommended that a top research priority should be development of composite methods. Pilot studies on the use of high-resolution imagery and passive acoustics for development of indices of relative abundance to be used in composite methods should be part of any future ship-based survey efforts

The TOP-SCOPE Survey of Planck Galactic Cold Clumps : Survey Overview and Results of an Exemplar Source, PGCC G26.53+0.17

The low dust temperatures (<14 K) of Planck Galactic cold clumps (PGCCs) make them ideal targets to probe the initial conditions and very early phase of star formation. "TOP-SCOPE" is a joint survey program targeting similar to 2000 PGCCs in J = 1-0 transitions of CO isotopologues and similar to 1000 PGCCs in 850 mu m continuum emission. The objective of the "TOP-SCOPE" survey and the joint surveys (SMT 10 m, KVN 21 m, and NRO 45 m) is to statistically study the initial conditions occurring during star formation and the evolution of molecular clouds, across a wide range of environments. The observations, data analysis, and example science cases for these surveys are introduced with an exemplar source, PGCC G26.53+0.17 (G26), which is a filamentary infrared dark cloud (IRDC). The total mass, length, and mean line mass (M/L) of the G26 filament are similar to 6200 M-circle dot, similar to 12 pc, and similar to 500 M-circle dot pc(-1), respectively. Ten massive clumps, including eight starless ones, are found along the filament. The most massive clump as a whole may still be in global collapse, while its denser part seems to be undergoing expansion owing to outflow feedback. The fragmentation in the G26 filament from cloud scale to clump scale is in agreement with gravitational fragmentation of an isothermal, nonmagnetized, and turbulent supported cylinder. A bimodal behavior in dust emissivity spectral index (beta) distribution is found in G26, suggesting grain growth along the filament. The G26 filament may be formed owing to large-scale compression flows evidenced by the temperature and velocity gradients across its natal cloud.Peer reviewe

Magnetic Fields toward Ophiuchus-B Derived from SCUBA-2 Polarization Measurements

We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope in Hawaii. This work is part of the B-fields in Star-forming Region Observations survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over ~0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically toward the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large-scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630 ± 410 μG in the Oph-B2 sub-clump using a Davis–Chandrasekhar–Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio λ = 1.6 ± 1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical