An environmentally based growth model that uses finite difference calculus with maximum likelihood method: its application to the brackish water bivalve Corbicula japonica in Lake Abashiri, Japan

We present a growth analysis model that combines large amounts of environmental data with limited amounts of biological data and apply it to Corbicula japonica. The model uses the maximum-likelihood method with the Akaike information criterion, which provides an objective criterion for model selection. An adequate distribution for describing a single cohort is selected from available probability density functions, which are expressed by location and scale parameters. Daily relative increase rates of the location parameter are expressed by a multivariate logistic function with environmental factors for each day and categorical variables indicating animal ages as independent variables. Daily relative increase rates of the scale parameter are expressed by an equation describing the relationship with the daily relative increase rate of the location parameter. Corbicula japonica grows to a modal shell length of 0.7 mm during the first year in Lake Abashiri. Compared with the attain-able maximum size of about 30 mm, the growth of juveniles is extremely slow because their growth is less susceptible to environmental factors until the second winter. The extremely slow growth in Lake Abashiri could be a geographical genetic variation within C. japonica.

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.

Integrated conceptual ecosystem model development for the Florida Keys/Dry Tortugas coastal marine ecosystem: MARine Estuarine goal Setting (MARES) for South Florida

The overall goal of the MARine and Estuarine goal Setting (MARES) project for South Florida is “to reach a science-based consensus about the defining characteristics and fundamental regulating processes of a South Florida coastal marine ecosystem that is both sustainable and capable of providing the diverse ecosystem services upon which our society depends.” Through participation in a systematic process of reaching such a consensus, science can contribute more directly and effectively to the critical decisions being made by both policy makers and by natural resource and environmental management agencies. The document that follows briefly describes the MARES project and this systematic process. It then describes in considerable detail the resulting output from the first two steps in the process, the development of conceptual diagrams and an Integrated Conceptual Ecosystem Model (ICEM) for the first subregion to be addressed by MARES, the Florida Keys/Dry Tortugas (FK/DT). What follows with regard to the FK/DT is the input received from more than 60 scientists, agency resource managers, and representatives of environmental organizations beginning with a workshop held December 9-10, 2009 at Florida International University in Miami, Florida.

Integrated conceptual ecosystem model development for the southeast Florida coastal marine ecosystem: MARine Estuarine goal Setting (MARES) for South Florida

The overall goal of the MARES (MARine and Estuarine goal Setting) project for South Florida is “to reach a science-based consensus about the defining characteristics and fundamental regulating processes of a South Florida coastal marine ecosystem that is both sustainable and capable of providing the diverse ecosystem services upon which our society depends.” Through participation in a systematic process of reaching such a consensus, science can contribute more directly and effectively to the critical decisions being made both by policy makers and by natural resource and environmental management agencies. The document that follows briefly describes MARES overall and this systematic process. It then describes in considerable detail the resulting output from the first step in the process, the development of an Integrated Conceptual Ecosystem Model (ICEM) for the third subregion to be addressed by MARES, the Southeast Florida Coast (SEFC). What follows with regard to the SEFC relies upon the input received from more than 60 scientists, agency resource managers, and representatives of environmental organizations during workshops held throughout 2009–2012 in South Florida.

Integrated conceptual ecosystem model development for the southwest Florida shelf coastal marine ecosystem: MARine Estuarine goal Setting (MARES) for South Florida

The overall goal of the MARine and Estuarine goal Setting (MARES) project for South Florida is “to reach a science-based consensus about the defining characteristics and fundamental regulating processes of a South Florida coastal marine ecosystem that is both sustainable and capable of providing the diverse ecosystem services upon which our society depends.” Through participation in a systematic process of reaching such a consensus, science can contribute more directly and effectively to the critical decisions being made by both policy makers and by natural resource and environmental management agencies. The document that follows briefly describes the MARES project and this systematic process. It then describes in considerable detail the resulting output from the first two steps in the process, the development of conceptual diagrams and an Integrated Conceptual Ecosystem Model (ICEM) for the second subregion to be addressed by MARES, the Southwest Florida Shelf (SWFS). What follows with regard to the SWFS is the input received from more than 60 scientists, agency resource managers, and representatives of environmental organizations beginning with a workshop held August 19-20, 2010 at Florida Gulf Coast University in Fort Myers, Florida.

Study of the development and evolution of seasonal thermocline and turbulent processes in the northern parts of Persian Gulf using numerical modeling

The Persian Gulf (PG) is a semi-enclosed shallow sea which is connected to open ocean through the Strait of Hormuz. Thermocline as a suddenly decrease of temperature in subsurface layer in water column leading to stratification happens in the PG seasonally. The forcing comprise tide, river inflow, solar radiation, evaporation, northwesterly wind and water exchange with the Oman Sea that influence on this process. In this research, analysis of the field data and a numerical (Princeton Ocean Model, POM) study on the summer thermocline development in the PG are presented. The Mt. Mitchell cruise 1992 salinity and temperature observations show that the thermocline is effectively removed due to strong wind mixing and lower solar radiation in winter but is gradually formed and developed during spring and summer; in fact as a result of an increase in vertical convection through the water in winter, vertical gradient of temperature is decreased and thermocline is effectively removed. Thermocline development that evolves from east to west is studied using numerical simulation and some existing observations. Results show that as the northwesterly wind in winter, at summer transition period, weakens the fresher inflow from Oman Sea, solar radiation increases in this time interval; such these factors have been caused the thermocline to be formed and developed from winter to summer even over the northwestern part of the PG. The model results show that for the more realistic monthly averaged wind experiments the thermocline develops as is indicated by summer observations. The formation of thermocline also seems to decrease the dissolved oxygen in water column due to lack of mixing as a result of induced stratification. Over most of PG the temperature difference between surface and subsurface increases exponentially from March until May. Similar variations for salinity differences are also predicted, although with smaller values than observed. Indeed thermocline development happens more rapidly in the Persian Gulf from spring to summer. Vertical difference of temperature increases to 9 centigrade degrees in some parts of the case study zone from surface to bottom in summer. Correlation coefficients of temperature and salinity between the model results and measurements have been obtained 0.85 and 0.8 respectively. The rate of thermcline development was found to be between 0.1 to 0.2 meter per day in the Persian Gulf during the 6 months from winter to early summer. Also it is resulted from the used model that turbulence kinetic energy increases in the northwestern part of the PG from winter to early summer that could be due to increase in internal waves activities and stability intensified through water column during this time.

A catch-free stock assessment model with application to goliath grouper (Epinephelus itajara) off southern Florida

Many modern stock assessment methods provide the machinery for determining the status of a stock in relation to certain reference points and for estimating how quickly a stock can be rebuilt. However, these methods typically require catch data, which are not always available. We introduce a model-based framework for estimating reference points, stock status, and recovery times in situations where catch data and other measures of absolute abundance are unavailable. The specif ic estimator developed is essentially an age-structured production model recast in terms relative to pre-exploitation levels. A Bayesian estimation scheme is adopted to allow the incorporation of pertinent auxiliary information such as might be obtained from meta-analyses of similar stocks or anecdotal observations. The approach is applied to the population of goliath grouper (Epinephelus itajara) off southern Florida, for which there are three indices of relative abundance but no reliable catch data. The results confirm anecdotal accounts of a marked decline in abundance during the 1980s followed by a substantial increase after the harvest of goliath grouper was banned in 1990. The ban appears to have reduced fishing pressure to between 10% and 50% of the levels observed during the 1980s. Nevertheless, the predicted fishing mortality rate under the ban appears to remain substantial, perhaps owing to illegal harvest and depth-related release mortality. As a result, the base model predicts that there is less than a 40% chance that the spawning biomass will recover to a level that would produce a 50% spawning potential ratio.

Year-to-year changes of vertical temperature distribution in the California Current region: 1954 to 1986

Studies by Enfield and Allen (1980), McLain et al (1985), and others have shown that anomalously warm years in the northern coastal California Current correspond to El Niño conditions in the equatorial Pacific Ocean. Ocean model studies suggest a mechanical link between the northern coastal California Current and the equatorial ocean through long waves that propagate cyclonically along the ocean boundary (McCreary 1976; Clarke 1983; Shriver et al 1991). However, distinct observational evidence of such an oceanic connection is not extensive. Much of the supposed El Niño variation in temperature and sea level data from the coastal California Current region can be associated with the effects of anomalously intense north Pacific atmospheric cyclogenesis, which is frequently augmented during El Niño years (Wallace and Gutzler 1981; Simpson 1983; Emery and Hamilton 1984). This study uses time series of ocean temperature data to distinguish between locally forced effects, initiated by north Pacific atmospheric changes, and remotely forced effects, initiated by equatorial Pacific atmospheric changes related to El Niño events.

An evaluation of back-calculation methodology using simulated otolith data

I simulated somatic growth and accompanying otolith growth using an individual-based bioenergetics model in order to examine the performance of several back-calculation methods. Four shapes of otolith radius-total length relations (OR-TL) were simulated. Ten different back-calculation equations, two different regression models of radius length, and two schemes of annulus selection were examined for a total of 20 different methods to estimate size at age from simulated data sets of length and annulus measurements. The accuracy of each of the twenty methods was evaluated by comparing the back-calculated length-at-age and the true length-at-age. The best back-calculation technique was directly related to how well the OR-TL model fitted. When the OR-TL was sigmoid shaped and all annuli were used, employing a least squares linear regression coupled with a log-transformed Lee back-calculation equation (y-intercept corrected) resulted in the least error; when only the last annulus was used, employing a direct proportionality back-calculation equation resulted in the least error. When the OR-TL was linear, employing a functional regression coupled with the Lee back-calculation equation resulted in the least error when all annuli were used, and also when only the last annulus was used. If the OR-TL was exponentially shaped, direct substitution into the fitted quadratic equation resulted in the least error when all annuli were used, and when only the last annulus was used. Finally, an asymptotically shaped OR-TL was best modeled by the individually corrected Weibull cumulative distribution function when all annuli were used, and when only the last annulus was used.