Accounting for individual variability in the von Bertalanffy growth model

Estimation of von Bertalanffy growth parameters has received considerable attention in fisheries research. Since Sainsbury (1980, Can. J. Fish. Aquat. Sci. 37: 241-247) much of this research effort has centered on accounting for individual variability in the growth parameters. In this paper we demonstrate that, in analysis of tagging data, Sainsbury's method and its derivatives do not, in general, satisfactorily account for individual variability in growth, leading to inconsistent parameter estimates (the bias does not tend to zero as sample size increases to infinity). The bias arises because these methods do not use appropriate conditional expectations as a basis for estimation. This bias is found to be similar to that of the Fabens method. Such methods would be appropriate only under the assumption that the individual growth parameters that generate the growth increment were independent of the growth parameters that generated the initial length. However, such an assumption would be unrealistic. The results are derived analytically, and illustrated with a simulation study. Until techniques that take full account of the appropriate conditioning have been developed, the effect of individual variability on growth has yet to be fully understood.

A note on the von Bertalanffy growth function concerning the allocation of surplus energy to reproduction

We propose an extended form of the von Bertalanffy growth function (VBGF), where the allocation of surplus energy to reproduction is considered. Any function can be used in our model to describe the ratio of energy allocation for reproduction to that for somatic growth. As an example, two models for energy allocation were derived: a step-function and a logistic function. The extended model can jointly describe growth in adult and juvenile stages. The change in growth rate between the two stages can be either gradual or steep; the latter gives a biphasic VBGF. The results of curve fitting indicated that a consideration of reproductive energy is meaningful for model extension. By controlling parameter values, our comprehensive model gives various growth curve shapes ranging from indeterminate to determinate growth. An increase in the number of parameters is unavoidable in practical applications of this new model. Additional information on reproduction will improve the reliability of model estimates.

The von Bertalanffy growth curve: when a good fit is not good enough

The von Bertalanffy growth function is used for length based analysis of growth and mortality patterns for management of fisheries. However, certain fish have growth patterns that the VBGF may not be able to describe adequately.e.g. the Acanthurus lineatus in Samoa. In such cases a two phase VBGF may be a useful approach.

Food conversion efficiency and the von Bertalanffy growth function 1: a modification of Pauly's model

A simple modification of Pauly's model for relating food conversion efficiency (K sub(1)) and body weight is proposed. The key parameter is an index to how efficiently food can be absorbed; the other parameter is related to the surface-limiting growth, an important component of von Bertalanff's and Pauly's theories of fish growth.

Food conversion efficiency and the von Bertalanffy growth function. Part II and conclusion: Extension of the new model to the generalized von Bertalanffy growth function

The simple model relating food conversion efficiency (K sub(1)) to body weight derived from the theoretical concepts behind von Bertalanffy's growth model, is extended here in the context of Pauly's generalization of that model. The exponent, which was fixed to 1/3 in the simple model, is in the extended model equivalent to 1-d, with d being the weight exponent of the anabolism term in Pauly's growth model. This makes the model applicable to fish for which the assumptions of the original (special) version of von Bertalanffy's growth model are violated.

A von Bertalanffy Based Model for the Estimation of Oyster (Crassostrea virginica) Growth on Restored Oyster Reefs in Chesapeake Bay

A model to estimate the mean monthly growth of Crassostrea virginica oysters in Chesapeake Bay was developed. This model is based on the classic von Bertalanffy growth function, however the growth constant is changed every monthly timestep in response to short term changes in temperature and salinity. Using a dynamically varying growth constant allows the model to capture seasonal oscillations in growth, and growth responses to changing environmental conditions that previous applications of the von Bertalanffy model do not capture. This model is further expanded to include an estimation of Perkinsus marinus impacts on growth rates as well as estimations of ecosystem services provided by a restored oyster bar over time. The model was validated by comparing growth estimates from the model to oyster shell height observations from a variety of restoration sites in the upper Chesapeake Bay. Without using the P. marinus impact on growth, the model consistently overestimates mean oyster growth. However, when P. marinus effects are included in the model, the model estimates match the observed mean shell height closely for at least the first 3 years of growth. The estimates of ecosystem services suggested by this model imply that even with high levels of mortality on an oyster reef, the ecosystem services provided by that reef can still be maintained by growth for several years. Because larger oyster filter more water than smaller ones, larger oysters contribute more to the filtration and nutrient removal ecosystem services of the reef. Therefore a reef with an abundance of larger oysters will provide better filtration and nutrient removal. This implies that if an oyster restoration project is trying to improve water quality through oyster filtration, it is important to maintain the larger older oysters on the reef.

Big bang bifurcations in von Bertalanffy’s dynamics with strong and weak Allee effects

The main purpose of this work was to study population dynamic discrete models in which the growth of the population is described by generalized von Bertalanffy's functions, with an adjustment or correction factor of polynomial type. The consideration of this correction factor is made with the aim to introduce the Allee effect. To the class of generalized von Bertalanffy's functions is identified and characterized subclasses of strong and weak Allee's functions and functions with no Allee effect. This classification is founded on the concepts of strong and weak Allee's effects to population growth rates associated. A complete description of the dynamic behavior is given, where we provide necessary conditions for the occurrence of unconditional and essential extinction types. The bifurcation structures of the parameter plane are analyzed regarding the evolution of the Allee limit with the aim to understand how the transition from strong Allee effect to no Allee effect, passing through the weak Allee effect, is realized. To generalized von Bertalanffy's functions with strong and weak Allee effects is identified an Allee's effect region, to which is associated the concepts of chaotic semistability curve and Allee's bifurcation point. We verified that under some sufficient conditions, generalized von Bertalanffy's functions have a particular bifurcation structure: the big bang bifurcations of the so-called box-within-a-box type. To this family of maps, the Allee bifurcation points and the big bang bifurcation points are characterized by the symmetric of Allee's limit and by a null intrinsic growth rate. The present paper is also a significant contribution in the framework of the big bang bifurcation analysis for continuous 1D maps and unveil their relationship with the explosion birth and the extinction phenomena.

Generalised growth models for aquatic species with an application to blacklip abalone (Haliotis rubra)

This paper presents a maximum likelihood method for estimating growth parameters for an aquatic species that incorporates growth covariates, and takes into consideration multiple tag-recapture data. Individual variability in asymptotic length, age-at-tagging, and measurement error are also considered in the model structure. Using distribution theory, the log-likelihood function is derived under a generalised framework for the von Bertalanffy and Gompertz growth models. Due to the generality of the derivation, covariate effects can be included for both models with seasonality and tagging effects investigated. Method robustness is established via comparison with the Fabens, improved Fabens, James and a non-linear mixed-effects growth models, with the maximum likelihood method performing the best. The method is illustrated further with an application to blacklip abalone (Haliotis rubra) for which a strong growth-retarding tagging effect that persisted for several months was detected

Improved estimation of size-transition matrices using tag-recapture data

We derive a new method for determining size-transition matrices (STMs) that eliminates probabilities of negative growth and accounts for individual variability. STMs are an important part of size-structured models, which are used in the stock assessment of aquatic species. The elements of STMs represent the probability of growth from one size class to another, given a time step. The growth increment over this time step can be modelled with a variety of methods, but when a population construct is assumed for the underlying growth model, the resulting STM may contain entries that predict negative growth. To solve this problem, we use a maximum likelihood method that incorporates individual variability in the asymptotic length, relative age at tagging, and measurement error to obtain von Bertalanffy growth model parameter estimates. The statistical moments for the future length given an individual's previous length measurement and time at liberty are then derived. We moment match the true conditional distributions with skewed-normal distributions and use these to accurately estimate the elements of the STMs. The method is investigated with simulated tag-recapture data and tag-recapture data gathered from the Australian eastern king prawn (Melicertus plebejus).

Latitudinal and seasonal effects on growth of the Australian eastern king prawn (Melicertus plebejus)

The growth of the Australian eastern king prawn (Melicertus plebejus) is understood in greater detail by quantifying the latitudinal effect. The latitudinal effect is the change in the species' growth rate during migration. Mark-recapture data (N = 1635, latitude 22.21 degrees S-34.00 degrees S) presents northerly movement of the eastern king prawn, with New South Wales prawns showing substantial average movement of 140 km (standard deviation: 176 km) north. A generalized von Bertalanffy growth model framework is used to incorporate the latitudinal effect together with the canonical seasonal effect. Applying this method to eastern king prawn mark-recapture data guarantees consistent estimates for the latitudinal and seasonal effects. For M. plebejus, it was found that growth rate peaks on 25 and 29 January for males and females, respectively; is at a minimum on 27 and 31 July, respectively; and that the shape parameter, k (per year), changes by -0.0236 and -0.0556 every 1 degree of latitude south increase for males and females, respectively.

Study of Pinna nobilis growth from inner record: How biased are posterior adductor muscle scars estimates?

Previous studies have shown that the external growth records of the posterior adductor muscle scar (PAMS) of the bivalve Pinna nobilis are incomplete and do not produce accurate age estimations. We have developed a new methodology to study age and growth using the inner record of the PAMS, which avoids the necessity of costly in situ shell measurements or isotopic studies. Using the inner record we identified the positions of PAMS previously obscured by nacre and estimated the number of missing records in adult specimens with strong abrasion of the calcite layer in the anterior portion of the shell. The study of the PAMS and inner record of two shells that were 6 years old when collected showed that only 2 and 3 PAMS were observed, while 6 inner records could be counted, thus confirming our working methodology. Growth parameters of a P. nobilis population located in Moraira, Spain (western Mediterranean) were estimated with the new methodology and compared to those obtained using PAMS data and in situ measurements. For the comparisons, we applied different models considering the data alternatively as length-at-age (LA) and tag-recapture (TR). Among every method we tested to fit the Von Bertalanffy growth model, we observed that LA data from inner record fitted to the model using non-linear mixed effects and the estimation of missing records using the calcite width was the most appropriate. The equation obtained with this method, L = 573*(1 - e(-0.16(t-0.02))), is very similar to that calculated previously from in situ measurements for the same population.

A maximum likelihood approach for estimating growth from tag–recapture data

The Fabens method is commonly used to estimate growth parameters k and l infinity in the von Bertalanffy model from tag-recapture data. However, the Fabens method of estimation has an inherent bias when individual growth is variable. This paper presents an asymptotically unbiassed method using a maximum likelihood approach that takes account of individual variability in both maximum length and age-at-tagging. It is assumed that each individual's growth follows a von Bertalanffy curve with its own maximum length and age-at-tagging. The parameter k is assumed to be a constant to ensure that the mean growth follows a von Bertalanffy curve and to avoid overparameterization. Our method also makes more efficient use nf thp measurements at tno and recapture and includes diagnostic techniques for checking distributional assumptions. The method is reasonably robust and performs better than the Fabens method when individual growth differs from the von Bertalanffy relationship. When measurement error is negligible, the estimation involves maximizing the profile likelihood of one parameter only. The method is applied to tag-recapture data for the grooved tiger prawn (Penaeus semisulcatus) from the Gulf of Carpentaria, Australia.