A spatial population dynamics simulation model of tropical tunas using a habitat index based on environmental parameters

We are developing a spatial, multigear, multispecies population dynamics si mulation model for tropical tunas in the Pacific Ocean. The model is age-st ructured to account for growth and gear selectivity. It includes a tuna mov ement model based on a diffusion-advection equation in which the advective term is proportional to the gradient of a habitat index. The monthly geogra phical distribution of recruitment is defined by assuming that spawning occ urs in areas where sea surface temperature is above 25 degrees C. During th e first 3 months of their life, simulated tunas are transported by oceanic currents, after which movement is conditioned by gradients in the habitat i ndex. Independent estimates of natural mortality rates and population size from large-scale tagging experiments carried out by the Secretariat of the Pacific Community are used in the simulations. The habitat index consists o f components due to forage density and sea surface temperature, both of whi ch are suspected to play major roles in determining tuna distribution. Beca use direct observations of forage are not available on a basin scale, we de veloped a submodel to simulate the surface tuna forage production (Lehodey et al., 1998). At present, only skipjack (Katsuwonus pelamis; a surface tun a species caught by purse seine and by pole-and-line) is considered, at a 1 degrees-square resolution and on a monthly climatological time series. Des pite the simplicity of the model and the limitations of the data used, the simulation model is able to predict a distribution of skipjack catch rates, of the different fleets involved in the fishery, that is fairly consistent with observations.