La campagne MEDITS est une campagne de chalutage de fond standardisée et coordonnée à l'échelle internationale ciblant les ressources démersales et benthiques en mer Méditerranée et qui se déroule tous les ans au printemps depuis 1994. Chaque année, 65 stations fixes sont échantillonnées en moyenne dans le Golfe du Lion. L’identification de chaque individu capturé dans les traits de chalut est réalisée au niveau spécifique puis les poids et le nombre d'individus par espèce sont déterminés. La taille, le sexe et la maturité sont récoltés quand cela est prévu en fonction du protocole défini pour chaque espèce.

Les espèces correctement échantillonnées et suffisamment abondantes, c'est à dire possédant une occurrence minimum supérieure ou égale à 20% entre 1994 et 2019, ont bénéficié d'une cartographie annuelle sur cette période.

Une préparation des données a été préalablement effectuée avant la création des cartes de distribution. Le calcul de la surface balayée par trait de chalut a été effectué afin de calculer la densité totale (nombre individus/km²) par station, par année et par espèce. La densité a ensuite été moyennée par station et par espèce sur la période d'étude. Les moyennes des latitudes et longitudes ont été calculées par trait de chalut. De façon à travailler dans un système orthonormé, la longitude a été transformée (Lon*cos(Lat*pi/180)).

Pour les espèces bénéficiant d'au moins 80% d’observations de la distribution en taille sur l’ensemble des observations disponibles, la différenciation entre les juvéniles et les adultes a été effectuée. Les stades de maturité sont échantillonnés pour certaines espèces lors de la campagne MEDITS ce qui nous a permis de distinguer dans les données les individus matures et non matures. Dans le but de déterminer la taille à laquelle les espèces atteignent leur maturité sexuelle, une régression logistique a été utilisée pour étudier la relation entre le stade de maturité et la taille des individus. Ce type de modèle relie le statut de maturité (un facteur bimodal) à la taille et permet de déterminer la longueur à laquelle 50% des individus ont atteint la maturité sexuelle (LM50). Les données biologiques sur les stades de maturité collectées dans le cadre de la DFMAP «Data Collection Multi-Annual Programme » ont été utilisées pour déterminer la période de reproduction et identifier la taille à maturité sexuelle du merlu Merluccius merluccius et du rouget de vase Mullus barbatus. Pour les espèces ne bénéficiant pas de données de maturité, un travail bibliographique a été réalisé afin de déterminer la taille seuil. La taille à maturité sexuelle mise en évidence a servi à identifier la proportion des juvéniles présents dans chaque trait sur la base des spectres de taille disponibles. Cette proportion a été rapportée à la densité totale pour obtenir la densité des juvéniles et des adultes dans chaque trait de chalut par an. Un travail bibliographique a été effectué pour distinguer les adultes reproducteurs (servant de proxi des zones de frayères), les adultes non reproducteurs (en phase de croissance) pour les espèces de poissons avec les juvéniles représentant les zones de nourriceries.

La taille à maturité sexuelle de chaque espèce, la différenciation entre adultes reproducteurs et non reproducteurs chez les poissons, et les sources, sont indiqués ci-dessous:

Citharus linguatula : 10.5 cm, Juvénile/Adulte reproducteur (Garcia-Rodriguez & Esteban, 2000)

Eledone cirrhosa : 9.0 cm (Donnaloia et al., 2010; Orsi Relini et al., 2006)

Eutrigla gurnardus : 21.0 cm, Juvénile/Adulte reproducteur (Campillo, 1992)

Lepidorhombus boscii : 12.2 cm, Juvénile/Adulte non reproducteur (Vassilopoulou & Ondrias, 1999; Campillo, 1992)

Lophius budegassa : 40.8 cm, Juvénile/Adulte reproducteur (Colmenero et al., 2013)

Lophius piscatorius : 51.3 cm, Juvénile/Adulte reproducteur (Colmenero et al., 2017)

Merluccius merluccius : 29.0 cm, Juvénile/Adulte reproducteur (Evaluation des stocks : Jadaux et al., 2018b; Campillo, 1992)

Mullus barbatus : 8.6 cm, Juvénile/Adulte reproducteur (Evaluation des stocks : Jadaux et al., 2018a)

Octopus vulgaris 12.0 cm (González et al., 2011; Quetglas et al., 1998)

Scyliorhinus canicula : 40.5 cm, Juvénile/Adulte reproducteur (Données MEDITS; Capapé et al., 2008)

Trachurus trachurus : 16.0 cm, Juvénile/Adulte reproducteur (Abaunza et al., 2003)

Trisopterus capelanus : 13.0 cm, Juvénile/Adulte reproducteur (Campillo, 1992)

Concernant la distribution des espèces, les données ont d'abord été transformées par logarithme décimal log10(x+1) puis une analyse variographique et un krigeage ordinaire ont été réalisés par stade de vie pour créer les cartes de densité interpolée. Le variogramme, basé sur les données multi-annuelles (moyenne des densités par station), a été ajusté par une des fonctions continues autorisées dans le cadre de la géostatistique : effet pépitique, exponentiel, circulaire, sphérique, puissance ou gaussien, choisie sur la base du meilleur ajustement statistique. Les géostatistiques ne peuvent être en principe appliquées que si les données ont un comportement isotropique (varient de la même façon dans toutes les directions). De façon à vérifier cette condition, un variogramme directionnel dans quatre directions (0°, 45°, 90°, 135°) a également été produit pour détecter une possible anisotropie des données, ce qui n’a pas été le cas. Enfin, de façon à vérifier qu’aucune tendance longue portée ne vienne perturber significativement l’estimation krigée, une cross-validation a été réalisée a posteriori pour vérifier l’utilité de sa prise en compte, ce qui a été rejeté dans tous les cas. Un krigeage utilisant les paramètres du variogramme ajusté précédemment et avec une résolution de 0.1° va estimer les valeurs de densité annuelles aux emplacements non échantillonnés ainsi que l’estimation de la variance (ou erreur) de krigeage permettant d’avoir une mesure de l’incertitude liée à l’estimation. L’erreur de krigeage correspond ici à l’écart-type (racine carré) de la variance de krigeage. Une rastérisation avec une résolution de 0.05° et une distance maximum d’extrapolation de 0.2° est ensuite effectuée pour obtenir les cartes de densité interpolée annuelles par espèce. Une carte de variance de krigeage a été calculée pour chaque espèce et chaque pas de temps.