A coupled movement and bioenergetics model to explore the spawning migration of anchovy in the Bay of Biscay

Adult anchovies in the Bay of Biscay perform north to south migration from late winter to early summer for spawning. However, what triggers and drives the geographic shift of the population remains unclear and poorly understood. An individual-based fish model has been implemented to explore the potential mechanisms that control anchovy's movement routes toward its spawning habitats. To achieve this goal, two fish movement behaviors – gradient detection through restricted area search and kinesis – simulated fish response to its dynamic environment. A bioenergetics model was used to represent individual growth and reproduction along the fish trajectory. The environmental forcing (food, temperature) of the model was provided by a coupled physical–biogeochemical model. We followed a hypothesis-testing strategy to actualize a series of simulations using different cues and computational assumptions. The gradient detection behavior was found as the most suitable mechanism to recreate the observed shift of anchovy distribution under the combined effect of sea-surface temperature and zooplankton. In addition, our results suggested that southward movement occurred more actively from early April to middle May following favorably the spatio-temporal evolution of zooplankton and temperature. In terms of fish bioenergetics, individuals who ended up in the southern part of the bay presented better condition based on energy content, proposing the resulting energy gain as an ecological explanation for this migration. The kinesis approach resulted in a moderate performance, producing distribution pattern with the highest spread. Finally, model performance was not significantly affected by changes on the starting date, initial fish distribution and number of particles used in the simulations, whereas it was drastically influenced by the adopted cues.

Oceanic control of multidecadal variability in an idealized coupled GCM

Idealized ocean models are known to develop intrinsic multidecadal oscillations of the meridional overturning circulation (MOC). Here we explore the role of ocean–atmosphere interactions on this low-frequency variability. We use a coupled ocean–atmosphere model set up in a flat-bottom aquaplanet geometry with two meridional boundaries. The model is run at three different horizontal resolutions (4°, 2° and 1°) in both the ocean and atmosphere. At all resolutions, the MOC exhibits spontaneous variability on multidecadal timescales in the range 30–40 years, associated with the propagation of large-scale baroclinic Rossby waves across the Atlantic-like basin. The unstable region of growth of these waves through the long wave limit of baroclinic instability shifts from the eastern boundary at coarse resolution to the western boundary at higher resolution. Increasing the horizontal resolution enhances both intrinsic atmospheric variability and ocean–atmosphere interactions. In particular, the simulated atmospheric annular mode becomes significantly correlated to the MOC variability at 1° resolution. An ocean-only simulation conducted for this specific case underscores the disruptive but not essential influence of air–sea interactions on the low-frequency variability. This study demonstrates that an atmospheric annular mode leading MOC changes by about 2 years (as found at 1° resolution) does not imply that the low-frequency variability originates from air–sea interactions.