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. 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)).

La différenciation entre les juvéniles et les adultes pour chaque espèce n'a pas été réalisée à cause du manque du données de taille. Un travail bibliographique a été effectué pour identifier les adultes reproducteurs (servant de proxi des zones de frayères) des adultes non reproducteurs (en phase de croissance) pour les espèces de poissons.

La différenciation entre adultes reproducteurs et non reproducteurs chez les poissons est indiquée ci-dessous:

Espèce Fonction Sources taille et période de reproduction

Arnoglossus laterna Adulte reproducteur Ilkyaz et al., 2017

Callionymus maculatus Adulte reproducteur Burgeot et al, 1993

Capros aper Adulte reproducteur Yapici & Filiz, 2014

Cepola macrophthalma Adulte reproducteur Stergiou, 1999

Conger conger Adulte reproducteur Vallisneri et al., 2007

Lepidotrigla cavillone Adulte reproducteur Campillo, 1992

Serranus hepatus Adulte reproducteur voir "Source 1"

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. Un variogramme empirique permet d'observer la relation spatiale des points d'échantillonnage : il va calculer l’autocorrélation entre chaque paire de points dans des classes de distances croissantes et identifier si les points ont une organisation particulière dans l’espace. 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.