A meta-analytic approach to quantifying scientific uncertainty in stock assessments

Quantifying scientific uncertainty when setting total allowable catch limits for fish stocks is a major challenge, but it is a requirement in the United States since changes to national fisheries legislation. Multiple sources of error are readily identifiable, including estimation error, model specification error, forecast error, and errors associated with the definition and estimation of reference points. Our focus here, however, is to quantify the influence of estimation error and model specification error on assessment outcomes. These are fundamental sources of uncertainty in developing scientific advice concerning appropriate catch levels and although a study of these two factors may not be inclusive, it is feasible with available information. For data-rich stock assessments conducted on the U.S. west coast we report approximate coefficients of variation in terminal biomass estimates from assessments based on inversion of the assessment of the model’s Hessian matrix (i.e., the asymptotic standard error). To summarize variation “among” stock assessments, as a proxy for model specification error, we characterize variation among multiple historical assessments of the same stock. Results indicate that for 17 groundfish and coastal pelagic species, the mean coefficient of variation of terminal biomass is 18%. In contrast, the coefficient of variation ascribable to model specification error (i.e., pooled among-assessment variation) is 37%. We show that if a precautionary probability of overfishing equal to 0.40 is adopted by managers, and only model specification error is considered, a 9% reduction in the overfishing catch level is indicated.

Use of genetic data to assess the uncertainty in stock assessments due to the assumed stock structure: the case of albacore (Thunnus alalunga) from the Atlantic Ocean

Stock assessments can be problematic because of uncertainties associated with the data or because of simplified assumptions made when modeling biological processes (Rosenberg and Restrepo, 1995). For example, the common assumption in stock assessments that stocks are homogeneous and discrete (i.e., there is no migration between the stocks) is not necessarily true (Kell et al., 2004a, 2004b).

Simulierte Auswirkungen des mehrjährigen Bewirtschaftungsplans für die Dorschfischerei in der westlichen Ostsee und die Möglichkeit einer technischen Alternative

The cod stock in the Western Baltic Sea is assessed to be overfished regarding the definitions of the UN World Summit on Sustainable Development at Johannesburg in 2002. Thus, the European Fisheries Council enforced a multi-annual management plan in 2007. Our medium term simulations over the future 10 years assume similar stock productivity as compared with the past four decades and indicate that the goals of the management plan can be achieved through TAC and consistent effort regulations. Taking account of the uncertainty in the recruitment patterns, the target average fishing mortality of age groups 3 – 6 years of F = 0.6 per year as defined in the management plan is indicated to exceed sustainable levels consistent with high long term yields and low risk of depletion. The stipulated constraint of the annual TAC variations of ±15% will dominate future fisheries management and implies a high recovery potential of the stock through continued reductions in fishing mortality. The scientific assessment of sustainable levels of exploitation and consideration in the plan is strongly advised, taking account of uncertainties attributed to environmental and biological effects. We recommend our study to be complemented with economic impact assessments including effects on by-catch species, which have been disregarded in this study. It is further demonstrated, that the goals of the management plan can alternatively be achieved by mesh size adaptations. An alternative technical option of mesh size increases to realize the required reductions in fishing mortality provides avoidance of discards of undersized fish after a few years by means of improved selectivity, another important element of the Common Fisheries Policy. However, it is emphasized that technical regulations since 1990 failed to affect the by-catch and discards of juvenile cod. In any way, the meaningful implementation of the multiannual management plan through stringent control and enforcement appears critical.

Improving reporting of uncertainties in sea level rise assessments

Sea level rise (SLR) assessments are commonly used to identify the extent that coastal populations are at risk to flooding. However, the data and assumptions used to develop these assessments contain numerous sources and types of uncertainty, which limit confidence in the accuracy of modeled results. This study illustrates how the intersection of uncertainty in digital elevation models (DEMs) and SLR lead to a wide range of modeled outcomes. SLR assessments are then reviewed to identify the extent that uncertainty is documented in peer-reviewed articles. The paper concludes by discussing priorities needed to further understand SLR impacts. (PDF contains 4 pages)

A-SCALA: an age-structured statistical catch-at-length analysis for assessing tuna stocks in the eastern tropical Pacific Ocean

English: We describe an age-structured statistical catch-at-length analysis (A-SCALA) based on the MULTIFAN-CL model of Fournier et al. (1998). The analysis is applied independently to both the yellowfin and the bigeye tuna populations of the eastern Pacific Ocean (EPO). We model the populations from 1975 to 1999, based on quarterly time steps. Only a single stock for each species is assumed for each analysis, but multiple fisheries that are spatially separate are modeled to allow for spatial differences in catchability and selectivity. The analysis allows for error in the effort-fishing mortality relationship, temporal trends in catchability, temporal variation in recruitment, relationships between the environment and recruitment and between the environment and catchability, and differences in selectivity and catchability among fisheries. The model is fit to total catch data and proportional catch-at-length data conditioned on effort. The A-SCALA method is a statistical approach, and therefore recognizes that the data collected from the fishery do not perfectly represent the population. Also, there is uncertainty in our knowledge about the dynamics of the system and uncertainty about how the observed data relate to the real population. The use of likelihood functions allow us to model the uncertainty in the data collected from the population, and the inclusion of estimable process error allows us to model the uncertainties in the dynamics of the system. The statistical approach allows for the calculation of confidence intervals and the testing of hypotheses. We use a Bayesian version of the maximum likelihood framework that includes distributional constraints on temporal variation in recruitment, the effort-fishing mortality relationship, and catchability. Curvature penalties for selectivity parameters and penalties on extreme fishing mortality rates are also included in the objective function. The mode of the joint posterior distribution is used as an estimate of the model parameters. Confidence intervals are calculated using the normal approximation method. It should be noted that the estimation method includes constraints and priors and therefore the confidence intervals are different from traditionally calculated confidence intervals. Management reference points are calculated, and forward projections are carried out to provide advice for making management decisions for the yellowfin and bigeye populations. Spanish: Describimos un análisis estadístico de captura a talla estructurado por edad, A-SCALA (del inglés age-structured statistical catch-at-length analysis), basado en el modelo MULTIFAN- CL de Fournier et al. (1998). Se aplica el análisis independientemente a las poblaciones de atunes aleta amarilla y patudo del Océano Pacífico oriental (OPO). Modelamos las poblaciones de 1975 a 1999, en pasos trimestrales. Se supone solamente una sola población para cada especie para cada análisis, pero se modelan pesquerías múltiples espacialmente separadas para tomar en cuenta diferencias espaciales en la capturabilidad y selectividad. El análisis toma en cuenta error en la relación esfuerzo-mortalidad por pesca, tendencias temporales en la capturabilidad, variación temporal en el reclutamiento, relaciones entre el medio ambiente y el reclutamiento y entre el medio ambiente y la capturabilidad, y diferencias en selectividad y capturabilidad entre pesquerías. Se ajusta el modelo a datos de captura total y a datos de captura a talla proporcional condicionados sobre esfuerzo. El método A-SCALA es un enfoque estadístico, y reconoce por lo tanto que los datos obtenidos de la pesca no representan la población perfectamente. Además, hay incertidumbre en nuestros conocimientos de la dinámica del sistema e incertidumbre sobre la relación entre los datos observados y la población real. El uso de funciones de verosimilitud nos permite modelar la incertidumbre en los datos obtenidos de la población, y la inclusión de un error de proceso estimable nos permite modelar las incertidumbres en la dinámica del sistema. El enfoque estadístico permite calcular intervalos de confianza y comprobar hipótesis. Usamos una versión bayesiana del marco de verosimilitud máxima que incluye constreñimientos distribucionales sobre la variación temporal en el reclutamiento, la relación esfuerzo-mortalidad por pesca, y la capturabilidad. Se incluyen también en la función objetivo penalidades por curvatura para los parámetros de selectividad y penalidades por tasas extremas de mortalidad por pesca. Se usa la moda de la distribución posterior conjunta como estimación de los parámetros del modelo. Se calculan los intervalos de confianza usando el método de aproximación normal. Cabe destacar que el método de estimación incluye constreñimientos y distribuciones previas y por lo tanto los intervalos de confianza son diferentes de los intervalos de confianza calculados de forma tradicional. Se calculan puntos de referencia para el ordenamiento, y se realizan proyecciones a futuro para asesorar la toma de decisiones para el ordenamiento de las poblaciones de aleta amarilla y patudo.

Analysis of sex-specific spawning biomass per recruit of the sailfish (Istiophorus platypterus) in the waters off eastern Taiwan

Analyses of sex-specific yield per recruit and spawning stock biomass per recruit were conducted to evaluate the current status of the sailfish (Istiophorus platypterus) fishery in the waters off eastern Taiwan. Natural mortality rates estimated from Pauly’s empirical equation were 0.26/yr for females and 0.27/yr for males. The current fishing mortality rates were estimated as 0.24/yr and 0.43/yr for females and males, respectively, which are much lower than the estimated F0 .1 (0.62/yr and 0.79/yr for females and males, respectively) and FSSB40 (0.46/yr for females) which are commonly used as target reference points in fisheries management. The effects of the fishing mortality, natural mortality, and age at first capture on the estimates of biological reference points were evaluated by using the Monte Carlo simulation. The results indicate that failure to consider the uncertainty in parameters such as natural mortality or age at first capture may lead to the improper estimation of biological reference points. This study indicates the possibility of current fishing mortality exceeding the target biological reference points may be negligible for sailfish in the waters off eastern Taiwan. However, in view of the recent rapid increase in fishing effort, it is evident that the stock status and development of the fishery need to be closely monitore

Morphometric differentiation in small juveniles of the pink spotted shrimp (Farfantepenaeus brasiliensis) and the southern pink shrimp (F. notialis) in the Yucatan Peninsula, Mexico

The morphometric and morphological characters of the rostrum have been widely used to identify penaeid shrimp species (Heales et al., 1985; Dall et al., 1990; Pendrey et al., 1999). In this setting, one of the constraints in studies of penaeid shrimp populations has been the uncertainty in the identification of early life history stages, especially in coastal nursery habitats, where recruits and juveniles dominate the population (Dall et al., 1990; Pérez-Castañeda and Defeo, 2001). In the western Atlantic Ocean, Pérez-Farfante (1969, 1970, 1971a) described diagnostic characters of the genus Farfantepenaeus that allowed identification of individuals in the range of 8−20 mm CL (carapace length) on the basis of the following morphological features: 1) changes in the structure of the petasma and thelycum; 2) absence or presence of distomarginal spines in the ventral costa of the petasma; 3) the ratio between the keel height and the sulcus width of the sixth abdominal somite; 4) the shape and position of the rostrum with respect to the segments and flagellum of the antennule; and 5) the ratio between rostrum length (RL) and carapace length (RL/CL). In addition, she classified Farfantepenaeus into two groups according to the shape and position of the rostrum with respect to the segments and flagellum of the antennule and the ratio RL/CL: 1) F. duorarum and F. notialis: short rostrum, straight distally, and the proximodorsal margin convex, usually extending anteriorly to the end of distal antennular segment, sometimes reaching to proximal one-fourth of broadened portion of lateral antennular flagellum, with RL/CL <0.75; and 2) F. aztecus, F. brasiliensis, F. paulensis, and F. subtilis: long rostrum, usually almost straight along the entire length, extending anteriorly beyond the distal antennular segment, sometimes reaching to the distal one-third of broadened portion of lateral antennular flagellum, with RL/CL >0.80. Pérez-Farfante stressed that, for the recognition to species level of juveniles <10 mm CL, all the characters listed above should be considered because occasionally one alone may not prove to be diagnostic. However, the only characters that could be distinguished for small juveniles in the range 4−8 mm CL are those defined on the rostrum. Therefore, it has been almost impossible to identify and separate small specimens of Farfantepenaeus (Pérez-Farfante, 1970, 1971a; Pérez-Farfante and Kensley, 1997).

The Status of the United States Population of Night Shark, Carcharhinus signatus

Night sharks, Carcharhinus signatus, are an oceanic species generally occurring in outer continental shelf waters in the western North Atlantic Ocean including the Caribbean Sea and Gulf of Mexico. Although not targeted, night sharks make up a segment of the shark bycatch in the pelagic longline fishery. Historically, night sharks comprised a significant proportion of the artisanal Cuban shark fishery but today they are rarely caught. Although information from some fisheries has shown a decline in catches of night sharks, it is unclear whether this decline is due to changes in fishing tactics, market, or species identification. Despite the uncertainty in the decline, the night shark is currently listed as a species of concern due to alleged declines in abundance resulting from fishing effort, i.e. overutilization. To assess their relevance to the species of concern list, we collated available information on the night shark to provide an analysis of its status. Night shark landings were likely both over- and under-reported and thus probably did not reflect all commercial and recreational catches, and overall they have limited relevance to the current status of the species. Average size information has not changed considerably since the 1980’s based on information from the pelagic longline fishery when corrected for gear bias. Analysis of biological information indicates night sharks have intrinsic rates of increase (r) about 10% yr–1 and have moderate rebound potential and an intermediate generation time compared to other sharks. An analysis of trends in relative abundance from four data sources gave conflicting results, with one series in decline, two series increasing, and one series relatively flat. Based on the analysis of all currently available information, we believe the night shark does not qualify as a species of concern but should be retained on the prohibited species list as a precautionary approach to management until a more comprehensive stock assessment can be conducted.

Quantifying annual variation in catchability for commercial and research fishing

Biomass indices, from commercial catch per unit of effort (CPUE) or random trawl surveys, are commonly used in fisheries stock assessments. Uncertainty in such indices, often ex-pressed as a coefficient of variation (CV), has two components: observation error, and annual variation in catchability. Only the former can be estimated directly. As a result, the CVs used for these indices either ignore the annual-variation component or assume a value for it (often implicitly). Two types of data for New Zealand stocks were examined: 48 sets of residuals and catchability estimates from stock assessments using either CPUE or trawl survey indices; and biomass estimates from 17 time series of trawl surveys with between 4 and 25 species per time series. These data show clear evidence of significant annual variation in catchability. With the trawl survey data, catchability was detectably extreme for many species in about one year in six. The assessment data suggest that this annual variability typically has a CV of about 0.2. For commercial CPUE the variability is slightly less, and a typical total CV (including both components) of 0.15 to 0.2. This is much less than the values of 0.3 to 0.35 that have commonly been assumed in New Zealand. Some estimates of catchability are shown to be implausible.

Small scale turbulence and mixing in semi-enclosed seas (Persian Gulf)

Turbulence and internal waves are probably important in generating layered structures in frontal region of marine environments (e.g. near river plumes outflow into the sea). Here we investigate the role of normal modes of internal waves in generation of layered structure in a part of Persian Gulf where river plume inters and in some laboratory experiments. The model prediction and observations show that layers so formed have a thickness of about 2m based on salinity variations with depth, but layers (about 5m) based on horizontal velocity profiles. Laboratory experiments with a plume outflow in a Filling Box profile also generate normal mode layered structure with 21H=0.5 (where A is layer thickness and H is the plume depth). In these experiments as Re of the flow is smaller than the Re of field flow. The normal modes are substantially dissipated with depth. Typical values of dissipation factor is about 0(100). This factor for field observation is 0(10) which is still substantial. Qualitative comparison between layered structure in field and laboratory is good. It should be emphasized that field observation is for semi-enclosed seas but the laboratory experiments are for enclosed region. Hence some of the discrepancies in the results of two cases are inevitable. Layered structures in marine environments are also produced by double diffusive convection. In this region this should be studied separately.

A probability-based approach to setting annual catch levels.

The requirement of setting annual catch limits to prevent overfishing has been added to the Magnuson-Stevens Fishery Conservation and Management Reauthorization Act of 2006 (MSRA). Because this requirement is new, a body of applied scientific practice for deriving annual catch limits and accompanying targets does not yet exist. This article demonstrates an approach to setting levels of catch that is intended to keep the probability of future overfishing at a preset low level. The proposed framework is based on stochastic projection with uncertainty in population dynamics. The framework extends common projection methodology by including uncertainty in the limit reference point and in management implementation, and by making explicit the risk of overfishing that managers consider acceptable. The approach is illustrated with application to gag (Mycteroperca microlepis), a grouper that inhabits the waters off the southeastern United States. Although devised to satisfy new legislation of the MSRA, the framework has potential application to any fishery where the management goal is to limit the risk of overfishing by controlling catch.

The River Exe SAP Final Plan

This is the River Exe Salmon Action Plan Final document produced by the Environment Agency in 2003. This final Salmon Action Plan (SAP) for the River Exe catchment has been produced after consideration of feedback from external consultation. The actions presented within this Salmon Action Plan clarify the important issues and factors currently limiting the salmon stock on the river. An attempt has been made to cost these actions, identify possible sources of funding and to provide a timescale for action. This SAP aims also to promote long term collaboration between the Agency and other interested parties in managing the River Exe salmon stock and fisheries. The River Exe salmon population is currently judged to be passing its Conservation Limit. However, its apparent declining trend in egg deposition in the recent years and the high uncertainty in its stock assessment suggest the following actions as priorities: actions promoting good land management, maximising salmon natural spawning activity and protecting smolts throughout the Exe catchment. Also, the lack of information on salmon stocks and its habitat quality and availability is recognised as the main factor limiting the better management of salmon fisheries. The actions presented in this document are perceived as those required to address the important issues and factors limiting the salmon stock.

Quantification and reduction of unobserved mortality rates for snow, southern Tanner, and red king crabs (Chionoecetes opilio, C. bairdi, and Paralithodes camtschaticus) after encounters with trawls on the seafloor

Unobserved mortalities of nontarget species are among the most troubling and difficult issues associated with fishing, especially when those species are targeted by other fisheries. Of such concern are mortalities of crab species of the Bering Sea, which are exposed to bottom trawling from groundfish fisheries. Uncertainty in the management of these fisheries has been exacerbated by unknown mortality rates for crabs struck by trawls. In this study, the mortality rates for 3 species of commercially important crabs—red king crab, (Paralithodes camtschaticus), snow crab (Chionoecetes opilio) and southern Tanner crab (C. bairdi)—that encounter different components of bottom trawls were estimated through capture of crabs behind the bottom trawl and by evaluation of immediate and delayed mortalities. We used a reflex action mortality predictor to predict delayed mortalities. Estimated mortality rates varied by species and by the part of the trawl gear encountered. Red king crab were more vulnerable than snow or southern Tanner crabs. Crabs were more likely to die after encountering the footrope than the sweeps of the trawl, and higher death rates were noted for the side sections of the footrope than for the center footrope section. Mortality rates were ≤16%, except for red king crab that passed under the trawl wings (32%). Herding devices (sweeps) can expand greatly the area of seafloor from which flatfishes are captured, and they subject crabs in that additional area to lower (4–9%) mortality rates. Raising sweep cables off of the seafloor reduced red king crab mortality rates from 10% to 4%.

Monte Carlo simulation of streamflow of the Salt River, Arizona

EXTRACT (SEE PDF FOR FULL ABSTRACT): Evaluations of the impact of climate change (such as a greenhouse effect) upon water resources should represent both the expected change and the uncertainty in that expectation. Since water resources such as streamflow and reservoir levels depend on a variety of factors, each of which is subject to significant uncertainty, it is desirable to formulate methods of representing that uncertainty in the forcing factors and from this determine the uncertainty in the response variables of interest. We report here progress in the representation of the uncertainty in climate upon the uncertainty in the estimated hydrologic response.