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 biomasse totale (poids/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 de données de taille.

Concernant la distribution des invertébrés benthiques, les données de 2012 à 2019 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 pour créer les cartes de biomasse 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 biomasses 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 biomasse 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 biomasse 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.