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 possédant une occurrence moyenne supérieure ou égale à 20% sur la période 1994-2019 ont bénéficié d'une cartographie pluriannuelle (distribution moyenne sur l’ensemble de la 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). 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:

Chelidonichthys cuculus : 18.0 cm, Juvénile/Adulte reproducteur (Campillo, 1992; voir "Source 1")

Boops boops : 13.0 cm, Juvénile/Adulte reproducteur (Campillo, 1992)

Eledone moschata : 9.0 cm (Krstulović-Šifner & Vrgoc, 2009 et Campillo, 1992)

Illex coindetii : 14.9 cm (Données MEDITS; Campillo, 1992)

Micromesistius poutassou : 18.9 cm, Juvénile/Adulte non reproducteur (Serrat et al., 2019; Mir-Arguimbau et al., 2020)

Nephrops norvegicus : 3.4 cm (Données MEDITS; voir "Source 2")

Pagellus acarne : 19.9 cm, Juvénile/Adulte non reproducteur (Velasco et al., 2011 et Campillo, 1992)

Pagellus bogaraveo : 30.0 cm, Juvénile/Adulte non reproducteur (Mytilineou et al., 2013; Campillo, 1992)

Pagellus erythrinus : 16.5 cm, Juvénile/Adulte reproducteur (Busalacchi et al., 2014; Campillo, 1992)

Parapenaeus longirostris : 4.0 cm (Données MEDITS; Arculeo et al., 2014)

Phycis blennoides : 19.5 cm, Juvénile/Adulte non reproducteur (Rotllant et al., 2002)

Scomber scombrus : 20.7 cm, Juvénile/Adulte non reproducteur (Campillo, 1992)

Trachurus mediterraneus : 23.0 cm, Juvénile/Adulte reproducteur (Campillo, 1992)

Zeus faber : 40.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. 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, en utilisant les valeurs multi-annuelles, les valeurs de densité 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 la carte de densité pluriannuelle interpolée par espèce.