<table><thead><tr><th colspan="2">Famille</th><th rowspan="2">Nom vernaculaire</th><th rowspan="2">Code</th><th rowspan="2">HRH</th><th rowspan="2">Ain</th><th rowspan="2">BRH</th><th rowspan="2">SAO</th><th rowspan="2">ARD</th><th rowspan="2">DUR</th><th rowspan="2">Totaux</th></tr></thead><tbody><tr><td colspan="2">Nom scientifique</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td colspan="2">Tinca tinca</td><td>Tanche</td><td>TAN</td><td>542</td><td>38</td><td>1 300</td><td>212</td><td>9</td><td>284</td><td>2 385</td></tr><tr><td colspan="2">Parachondrostoma toxostoma</td><td>Toxostome</td><td>TOX</td><td>19</td><td>2</td><td>44</td><td></td><td>195</td><td>1 598</td><td>1 858</td></tr><tr><td colspan="2">Carassius gibelio</td><td>Carassin argenté</td><td>CAG</td><td>2</td><td></td><td>1 484</td><td>40</td><td></td><td>16</td><td>1 542</td></tr><tr><td colspan="2">Cyprinus carpio</td><td>Carpe commune</td><td>CCO</td><td>23</td><td>1</td><td>704</td><td>128</td><td>7</td><td>62</td><td>925</td></tr><tr><td colspan="2">Carassius auratus</td><td>Carassin doré</td><td>CAA</td><td>1</td><td></td><td>154</td><td></td><td>15</td><td></td><td>170</td></tr><tr><td colspan="2">Carassius carassius</td><td>Carassin commun</td><td>CAS</td><td>12</td><td></td><td>56</td><td></td><td></td><td>12</td><td>80</td></tr><tr><td colspan="2">Chondrostoma sp.</td><td>Chondrostome</td><td>CHO</td><td></td><td></td><td></td><td></td><td></td><td>68</td><td>68</td></tr><tr><td colspan="2">Cyprinus carpio</td><td>Carpe miroir</td><td>CMI</td><td>9</td><td></td><td>17</td><td></td><td></td><td></td><td>26</td></tr><tr><td colspan="2">Leuciscus idus</td><td>Ide mélanote</td><td>IDE</td><td></td><td></td><td>21</td><td>4</td><td></td><td></td><td>25</td></tr><tr><td colspan="2">Cyprinus carpio</td><td>Carpe cuir</td><td>CCU</td><td>1</td><td></td><td>1</td><td></td><td>1</td><td></td><td>3</td></tr><tr><td colspan="2">Leucaspius delineatus</td><td>Able de Heckel</td><td>ABH</td><td>2</td><td></td><td>1</td><td></td><td></td><td></td><td>3</td></tr><tr><td colspan="2">Esocidae</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td colspan="2">Esox lucius</td><td>Brochet</td><td>BRO</td><td>838</td><td>76</td><td>263</td><td>13</td><td></td><td>56</td><td>1 246</td></tr><tr><td colspan="2">Gasterosteidae</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td colspan="2">Gasterosteus gymnurus</td><td>Epinoche</td><td>EPI</td><td>1 059</td><td>47</td><td>246</td><td></td><td></td><td>2</td><td>1 354</td></tr><tr><td colspan="2">Ictaluridae</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td colspan="2">Ameiurus melas</td><td>Poisson chat</td><td>PCH</td><td>146</td><td>5</td><td>5 086</td><td>252</td><td>1</td><td>11</td><td>5 501</td></tr><tr><td colspan="2">Lotidae</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td colspan="2">Lota lota</td><td>Lote de rivière</td><td>LOT</td><td>105</td><td>1</td><td>1</td><td></td><td></td><td></td><td>107</td></tr><tr><td colspan="2">Moronidae</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td colspan="2">Dicentrarchus labrax</td><td>Loup</td><td>BAR</td><td></td><td></td><td>15</td><td></td><td></td><td></td><td>15</td></tr><tr><td colspan="2">Mugilidae</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td colspan="2">Liza ramada</td><td>Mulet porc</td><td>MUP</td><td></td><td></td><td>294</td><td></td><td></td><td>49</td><td>343</td></tr><tr><td colspan="2">Mugil cephalus</td><td>Mulet cabot</td><td>MUC</td><td></td><td></td><td>14</td><td></td><td></td><td>19</td><td>33</td></tr><tr><td colspan="2">Percidae</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td colspan="2">Perca fluviatilis</td><td>Perche commune</td><td>PER</td><td>2 406</td><td>316</td><td>2 718</td><td>1 092</td><td>19</td><td>312</td><td>6 863</td></tr><tr><td colspan="2">Gymnocephalus cernuus</td><td>Grémille</td><td>GRE</td><td>652</td><td>5</td><td>463</td><td>76</td><td></td><td>46</td><td>1 242</td></tr><tr><td colspan="2">Sander lucioperca</td><td>Sandre</td><td>SAN</td><td>40</td><td></td><td>586</td><td>86</td><td></td><td>39</td><td>751</td></tr><tr><td colspan="2">Zingel asper</td><td>Apron</td><td>APR</td><td></td><td>1</td><td></td><td></td><td></td><td>61</td><td>62</td></tr><tr><td colspan="2">Petromyzontidae</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td colspan="2">Lampetra planeri</td><td>Lamproie de Planer</td><td>LPP</td><td>7</td><td>3</td><td>4</td><td></td><td></td><td></td><td>14</td></tr><tr><td colspan="2">non identifiee</td><td>Lamproie (larve)</td><td>LPA</td><td></td><td></td><td>1</td><td></td><td></td><td></td><td>1</td></tr><tr><td colspan="2">Poeciliidae</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td colspan="2">Gambusia affinis</td><td>Gambusie</td><td>GAM</td><td></td><td></td><td>2</td><td></td><td></td><td>136</td><td>138</td></tr><tr><td colspan="2">Salmonidae</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td colspan="2">Salmo trutta</td><td>Truite commune</td><td>TRF</td><td>971</td><td>415</td><td>85</td><td></td><td>16</td><td>116</td><td>1 603</td></tr><tr><td colspan="2">Thymallus thymallus</td><td>Ombre commun</td><td>OBR</td><td>483</td><td>68</td><td></td><td></td><td></td><td></td><td>551</td></tr><tr><td colspan="2">Oncorhynchus mykiss</td><td>Truite arc en ciel</td><td>TAC</td><td>36</td><td></td><td>27</td><td></td><td>1</td><td>9</td><td>73</td></tr><tr><td colspan="2">Siluridae</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td colspan="2">Silurus glanis</td><td>Silure</td><td>SIL</td><td>454</td><td></td><td>598</td><td>96</td><td></td><td></td><td>1 148</td></tr><tr><td colspan="2">Nombre d'espèces</td><td>51</td><td>Total</td><td>184 671</td><td>24 070</td><td>386 246</td><td>24 696</td><td>26 612</td><td>70 974</td><td>717 269</td></tr></tbody></table>
La famille des Cyprinidae compte 23 espèces, regroupe un peu plus de 90% des effectifs et la majorité des espèces dominantes sur les cours d'eau. Des espèces abondantes de petite taille (ablette, bouvière, pseudorasbora, goujon, spirlin), à croissance rapide et à maturité précoce posent de réel problème d'identification des classes d'âge en l'absence de scalimétrie et/ou d'otolithométrie. Les nombreuses espèces à ponte multiple, parmi lesquelles se trouvent les plus petites espèces, génèrent plusieurs cohortes au cours de la saison de reproduction à l'origine d'une grande dispersion des tailles et des chevauchements avec celles des individus de l'année précédente (RINCHARD & KESTEMONT, 1996; RINCHARD et al., 1998).

Parmi les espèces longévives les plus abondantes (gardon, chevaine, barbeau fluviatile, hotu, brème bordelière et vandoise), toutes n'ont pas une répartition homogène sur l'axe rhodanien. La vandoise n'est trouvée que sur le Haut-Rhône. La brème bordelière est nettement plus abondante sur le Bas-Rhône. La répartition de deux grands rhéophiles (barbeau et hotu) n'est pas uniforme dans le temps sur le Bas-Rhône puisque leur réapparition effective dans les chroniques du Bas-Rhône a été postérieure aux crues 1933-1994.

Finalement, il ne reste que deux espèces euryèces parmi les plus abondantes dans ces pêches faites sur l'axe rhodanien : le gardon et le chevaine (Figure 5). Ces deux espèces étant typiquement associées aux habitats rivulaires à tous les stades de vie (WOLTER & BISCHOFF, 2001), la méthode d'échantillonnage adoptée est globalement représentative de la répartition diurne de toutes les classes d'âge de ces deux populations. Ce ne serait pas le cas pour les brèmes commune et bordelière, deux espèces potamiques préférentiellement inféodées au chenal dès leur plus jeune âge.

Bien que toutes deux euryèces, tolérantes vis-à-vis de la dégradation physique du milieu et ayant une vaste distribution européenne, elles constituent deux modèles biologiques contrastés.

Parmi les aspects importants dans cette étude, leurs traits écologiques associés à la reproduction sont nettement différents, les plaçant parmi deux guildes reproductives généralistes mais distinctes (KRYZHANOVSKII, 1974; BALON, 1975).

Le chevaine est lithophile. Il dépose ses œufs sur des fonds pierreux localisés dans des secteurs courants (PONCIN, 1996). Le substrat abrite les larves photophobes jusqu'au début de leur nutrition externe. Elles ne possèdent pas de glandes cémentaires permettant leur fixation sur le substrat.

Le gardon est phyto-lithophile. La femelle dépose préférentiellement ses œufs sur des macrophytes des eaux calmes, mais peut pondre sur d'autres substrats si nécessaire tels que des débris ligneux, des graviers et des blocs. Les larves ont des glandes cémentaires leur permettant d'adhérer au substrat de ponte et d'éviter les eaux anoxiques proches du fond dans les milieux eutrophes.

Ces deux espèces présentent des stratégies de reproduction différentes. Le gardon est une espèce à ponte unique, réalisée en quelques jours. La fraie dure 2 à 3 jours dans le lac Léman (GILLET, 1989) si les conditions climatiques sont favorables. Si elles se dégradent, la ponte peut être interrompue temporairement, voire définitivement.

Le chevaine peut réaliser plusieurs pontes. Toutefois, parmi les espèces à pontes multiples, ce caractère est moins prononcé que chez l'ablette, la brème bordelière ou le goujon (RINCHARD & KESTEMONT, 1996; RINCHARD et al., 1997; RINCHARD et al., 1998). LIBOSVARSKI (1979) a montré que les ovaires des femelles contenaient plusieurs groupes d'ovules permettant la réalisation potentielle de trois pontes successives. Par ailleurs, PENAZ et al. (1978) ont observé trois fraies successives (du 10 au 12 mai, du 24 au 27 mai et du 2 au 5 juin 1977) dans la rivière Jihlava (République Tchèque). NUNN et al. (2007) signalent un minimum annuel de deux cohortes de chevaines 0+ dans la rivière Trent. Ces pontes multiples étaient plus apparentes en 1999, 2003 et 2004. Par contre, les pontes
La famille des Cyprinidae compte 23 espèces, regroupe un peu plus de 90% des effectifs et la majorité des espèces dominantes sur les cours d'eau. Des espèces abondantes de petite taille (ablette, bouvière, pseudorasbora, goujon, spirlin), à croissance rapide et à maturité précoce posent de réel problème d'identification des classes d'âge en l'absence de scalimétrie et/ou d'otolithométrie. Les nombreuses espèces à ponte multiple, parmi lesquelles se trouvent les plus petites espèces, génèrent plusieurs cohortes au cours de la saison de reproduction à l'origine d'une grande dispersion des tailles et des chevauchements avec celles des individus de l'année précédente (RINCHARD & KESTEMONT, 1996; RINCHARD et al., 1998).

Parmi les espèces longévives les plus abondantes (gardon, chevaine, barbeau fluviatile, hotu, brème bordelière et vandoise), toutes n'ont pas une répartition homogène sur l'axe rhodanien. La vandoise n'est trouvée que sur le Haut-Rhône. La brème bordelière est nettement plus abondante sur le Bas-Rhône. La répartition de deux grands rhéophiles (barbeau et hotu) n'est pas uniforme dans le temps sur le Bas-Rhône puisque leur réapparition effective dans les chroniques du Bas-Rhône a été postérieure aux crues 1933-1994.

Finalement, il ne reste que deux espèces euryèces parmi les plus abondantes dans ces pêches faites sur l'axe rhodanien : le gardon et le chevaine (Figure 5). Ces deux espèces étant typiquement associées aux habitats rivulaires à tous les stades de vie (WOLTER & BISCHOFF, 2001), la méthode d'échantillonnage adoptée est globalement représentative de la répartition diurne de toutes les classes d'âge de ces deux populations. Ce ne serait pas le cas pour les brèmes commune et bordelière, deux espèces potamiques préférentiellement inféodées au chenal dès leur plus jeune âge.

Bien que toutes deux euryèces, tolérantes vis-à-vis de la dégradation physique du milieu et ayant une vaste distribution européenne, elles constituent deux modèles biologiques contrastés.

Parmi les aspects importants dans cette étude, leurs traits écologiques associés à la reproduction sont nettement différents, les plaçant parmi deux guildes reproductives généralistes mais distinctes (KRYZHANOVSKII, 1974; BALON, 1975).

Le chevaine est lithophile. Il dépose ses œufs sur des fonds pierreux localisés dans des secteurs courants (PONCIN, 1996). Le substrat abrite les larves photophobes jusqu'au début de leur nutrition externe. Elles ne possèdent pas de glandes cémentaires permettant leur fixation sur le substrat.

Le gardon est phyto-lithophile. La femelle dépose préférentiellement ses œufs sur des macrophytes des eaux calmes, mais peut pondre sur d'autres substrats si nécessaire tels que des débris ligneux, des graviers et des blocs. Les larves ont des glandes cémentaires leur permettant d'adhérer au substrat de ponte et d'éviter les eaux anoxiques proches du fond dans les milieux eutrophes.

Ces deux espèces présentent des stratégies de reproduction différentes. Le gardon est une espèce à ponte unique, réalisée en quelques jours. La fraie dure 2 à 3 jours dans le lac Léman (GILLET, 1989) si les conditions climatiques sont favorables. Si elles se dégradent, la ponte peut être interrompue temporairement, voire définitivement.

Le chevaine peut réaliser plusieurs pontes. Toutefois, parmi les espèces à pontes multiples, ce caractère est moins prononcé que chez l'ablette, la brème bordelière ou le goujon (RINCHARD & KESTEMONT, 1996; RINCHARD et al., 1997; RINCHARD et al., 1998). LIBOSVARSKI (1979) a montré que les ovaires des femelles contenaient plusieurs groupes d'ovules permettant la réalisation potentielle de trois pontes successives. Par ailleurs, PENAZ et al. (1978) ont observé trois fraies successives (du 10 au 12 mai, du 24 au 27 mai et du 2 au 5 juin 1977) dans la rivière Jihlava (République Tchèque). NUNN et al. (2007) signalent un minimum annuel de deux cohortes de chevaines 0+ dans la rivière Trent. Ces pontes multiples étaient plus apparentes en 1999, 2003 et 2004. Par contre, les pontes
multiples n'étaient pas évidentes chaque année dans deux autres rivières dans lesquelles le recrutement annuel a été plus vraisemblablement le résultat d'une ponte unique.

La ponte multiple est à l'origine d'une plurimodalité marquée des distributions des tailles des jeunes chevaines de l'année (CARREL, 1986; NUNN et al., 2002) et d'une grande disparité des tailles des jeunes de l'année (BOLLAND et al., 2007).

<figure><img src="image_6.png" /><figcaption>Figure 5 - Distribution en classes de taille de 5 mm des gardons et chevaines capturés par pêche électrique (Haut-Rhône, Bas-Rhône) et paramètres synthétiques de la taille.</figcaption></figure>
multiples n'étaient pas évidentes chaque année dans deux autres rivières dans lesquelles le recrutement annuel a été plus vraisemblablement le résultat d'une ponte unique.

La ponte multiple est à l'origine d'une plurimodalité marquée des distributions des tailles des jeunes chevaines de l'année (CARREL, 1986; NUNN et al., 2002) et d'une grande disparité des tailles des jeunes de l'année (BOLLAND et al., 2007).

<figure><img src="image_6.png" /><figcaption>Figure 5 - Distribution en classes de taille de 5 mm des gardons et chevaines capturés par pêche électrique (Haut-Rhône, Bas-Rhône) et paramètres synthétiques de la taille.</figcaption></figure>
## 2.2 Données thermiques

Une convention a été mise en place entre eDF et Irstea Aix-en-Provence afin de disposer des données thermiques moyennes journalières vérifiées, validées et centralisées par le service Ressources en Eau, Branche Environnement, eDF-DTG de Grenoble. Les températures utilisées dans ce travail proviennent des stations permanentes du Rhône [Creys, Bugey amont, Jons, Loire-sur-Rhône, Saint-Alban amont, Seuil de Peyraud (RCC Péage-de-Roussillon), Cruas amont, Tricastin amont, Tricastin aval, Aramon amont], de la Saône (Couzon) et de l'Ain (Pont de Chazey) (LAUTERS & POIREL, 2001a, b).

La température du rejet de Tricastin, fournie par le CNPE de Tricastin, est une valeur calculée, fonction de la température du fleuve à l'amont, du débit du canal de Donzère-Mondragon et de la production énergétique. Les valeurs au droit du rejet sont environ de +10°C par rapport à l'amont. Le panache thermique est plaqué en rive droite du canal, constituant une zone de mélange très hétérogène. Le mélange complet des masses d'eau ne s'effectue qu'au passage de l'usine de Bollène. Dans ce secteur, nous avons deux stations de pêche (rejet et aval immédiat). Afin de considérer une température plus proche des conditions ambiantes locales, nous avons calculé une température moyenne entre la température amont et la température calculée des eaux réchauffées.

L'emplacement des stations de mesure eDF a été adapté pour renseigner au mieux la température de la masse d'eau du chenal principal. Or, de nombreux relevés biologiques sont localisés dans un contexte thermique distinct du chenal principal. Ils sont traités avec les relevés thermiques disponibles les plus proches, une approximation acceptée malgré la forte hétérogénéité spatiale de la température, en particulier des sections fluviales court-circuitées. La température printanière et estivale des eaux de la masse d'eau principale est généralement plus basse que celle de la frange utilisée par les jeunes stades, proche des rives et de la surface.

Les poissons sont aptes à rechercher des conditions thermiques optimales, notamment des eaux plus chaudes chez les jeunes stades. Les larves sont à même de se déplacer dans un gradient thermique pour optimiser leur croissance en fonction de leur stade de développement (GOLOVANOV & SMIRNOV, 2011). Par conséquent, quelles que soient les données thermiques utilisées, même les plus précises et les plus pertinentes, il y aura toujours une grande approximation des conditions réelles.

## 2.3 Données hydrologiques

Les chroniques hydrologiques proviennent des stations de la banque HYDRO, gérées soit par les deux grands exploitants de la ressource hydraulique, eDF et la Compagnie Nationale du Rhône, soit par les services de l'Etat : DREAL Rhône-Alpes, DREAL Provence-Alpes-Côte d'Azur, et le Service de Prévision des Crues Grand Delta. La donnée de base utilisée est le débit moyen journalier (m3/s) relevé au droit des stations limnigraphiques.

## 2.4 Le tri et l'identification des classes d'âge

### 2.4.1 Usage des distributions en classes de taille

En zone tempérée, la croissance des poissons est fortement saisonnalisée, fonction de la température de l'eau, la longueur du jour et la disponibilité alimentaire. Elle se traduit par une augmentation rapide de la taille et du poids au cours de l'été, un arrêt progressif lors de la baisse des températures et la diminution de la longueur du jour, puis une absence virtuelle durant les mois d'hiver avant la reprise printanière. Les données permettant d'établir des courbes de croissance saisonnière sont souvent obtenues par des échantillonnages réguliers au cours de l'année et à partir des histogrammes de taille, une technique couramment adoptée pour le suivi de la croissance des juvéniles dans tous les types de
## 2.2 Données thermiques

Une convention a été mise en place entre eDF et Irstea Aix-en-Provence afin de disposer des données thermiques moyennes journalières vérifiées, validées et centralisées par le service Ressources en Eau, Branche Environnement, eDF-DTG de Grenoble. Les températures utilisées dans ce travail proviennent des stations permanentes du Rhône [Creys, Bugey amont, Jons, Loire-sur-Rhône, Saint-Alban amont, Seuil de Peyraud (RCC Péage-de-Roussillon), Cruas amont, Tricastin amont, Tricastin aval, Aramon amont], de la Saône (Couzon) et de l'Ain (Pont de Chazey) (LAUTERS & POIREL, 2001a, b).

La température du rejet de Tricastin, fournie par le CNPE de Tricastin, est une valeur calculée, fonction de la température du fleuve à l'amont, du débit du canal de Donzère-Mondragon et de la production énergétique. Les valeurs au droit du rejet sont environ de +10°C par rapport à l'amont. Le panache thermique est plaqué en rive droite du canal, constituant une zone de mélange très hétérogène. Le mélange complet des masses d'eau ne s'effectue qu'au passage de l'usine de Bollène. Dans ce secteur, nous avons deux stations de pêche (rejet et aval immédiat). Afin de considérer une température plus proche des conditions ambiantes locales, nous avons calculé une température moyenne entre la température amont et la température calculée des eaux réchauffées.

L'emplacement des stations de mesure eDF a été adapté pour renseigner au mieux la température de la masse d'eau du chenal principal. Or, de nombreux relevés biologiques sont localisés dans un contexte thermique distinct du chenal principal. Ils sont traités avec les relevés thermiques disponibles les plus proches, une approximation acceptée malgré la forte hétérogénéité spatiale de la température, en particulier des sections fluviales court-circuitées. La température printanière et estivale des eaux de la masse d'eau principale est généralement plus basse que celle de la frange utilisée par les jeunes stades, proche des rives et de la surface.

Les poissons sont aptes à rechercher des conditions thermiques optimales, notamment des eaux plus chaudes chez les jeunes stades. Les larves sont à même de se déplacer dans un gradient thermique pour optimiser leur croissance en fonction de leur stade de développement (GOLOVANOV & SMIRNOV, 2011). Par conséquent, quelles que soient les données thermiques utilisées, même les plus précises et les plus pertinentes, il y aura toujours une grande approximation des conditions réelles.

## 2.3 Données hydrologiques

Les chroniques hydrologiques proviennent des stations de la banque HYDRO, gérées soit par les deux grands exploitants de la ressource hydraulique, eDF et la Compagnie Nationale du Rhône, soit par les services de l'Etat : DREAL Rhône-Alpes, DREAL Provence-Alpes-Côte d'Azur, et le Service de Prévision des Crues Grand Delta. La donnée de base utilisée est le débit moyen journalier (m3/s) relevé au droit des stations limnigraphiques.

## 2.4 Le tri et l'identification des classes d'âge

### 2.4.1 Usage des distributions en classes de taille

En zone tempérée, la croissance des poissons est fortement saisonnalisée, fonction de la température de l'eau, la longueur du jour et la disponibilité alimentaire. Elle se traduit par une augmentation rapide de la taille et du poids au cours de l'été, un arrêt progressif lors de la baisse des températures et la diminution de la longueur du jour, puis une absence virtuelle durant les mois d'hiver avant la reprise printanière. Les données permettant d'établir des courbes de croissance saisonnière sont souvent obtenues par des échantillonnages réguliers au cours de l'année et à partir des histogrammes de taille, une technique couramment adoptée pour le suivi de la croissance des juvéniles dans tous les types de
milieux aquatiques (FECHHELM et al., 1997; POWER et al., 2000; NUNN et al., 2002; RIBEIRO et al., 2004).

Pour la majorité des Cyprinidés, seules les plus jeunes classes d'âge peuvent être identifiées à partir des distributions de taille, essentiellement les classes d'âge 0+ (jeunes de l'année n), et 1+ (jeunes de l'année n-1). Généralement, un troisième groupe est distinct dans une distribution de taille (individus majoritairement d'âge 2+). Mais il présente un chevauchement considérable avec les classes d'âge proches (1+, 3+ et au-delà).

En l'absence d'une identification précise de l'âge par des techniques spécifiques appropriées (scalimétrie, otolithométrie), ce travail est basé sur une discrimination des âges à partir des histogrammes de taille.

Cette étape, coûteuse en temps, a été rendue possible par l'abondance des juvéniles dans les pêches électriques de rives, l'existence de campagnes pluriannuelles sur de nombreux sites, une relative homogénéité des dates d'échantillonnage dans les longues chroniques des sites industriels, et la complémentarité des informations recueillies par cette approche pluri scalaire, spatiale et temporelle.

Dans l'identification des jeunes de l'année, certaines études considèrent une valeur de taille « butoir » en-deçà de laquelle les individus sont considérés de l'année. Ce choix est discutable dans une étude prenant en compte d'une part l'axe rhodanien français orienté nord-sud, ses nombreux milieux, et des situations extrêmes représentées par les rejets thermiques ; et d'autre part l'usage de chroniques pluri-décennales alors que le contexte hydrologique et thermique est fortement structuré par le changement climatique actuel.

La multiplicité des situations observées n'autorise pas cette simplification extrême d'autant que l'intérêt de ce jeu de données réside dans la variabilité de la taille selon les années et les secteurs. A ce sujet, TEDESCO et al. (2009) préconisent la prise en compte des spécificités locales et notamment le gradient longitudinal, plutôt que de s'appuyer sur des valeurs globales issues de la bibliographie.

Par conséquent, chaque tableau associé à une espèce comporte l'ensemble des pêches par station et par date où l'espèce est présente.

## 2.4.2 Les problèmes de discrimination des classes d'âge

Malgré les atouts du jeu de données, le travail reste complexe du fait (1) de la forte variabilité des effectifs des juvéniles, (2) de la variabilité intra-spécifique de la croissance des juvéniles selon les années et les secteurs et (3) de la stratégie de ponte des espèces.

### 2.4.2.1 Les effectifs de juvéniles

Bien que les effectifs de juvéniles dominent dans les pêches électriques réalisées sur le Rhône, ils sont très variables dans l'espace, au cours de l'année, et d'une année à l'autre. Cette grande amplitude de variabilité a plusieurs origines.

Les conditions de reproduction peuvent ne pas avoir été favorables et les effectifs des pêches automnales sont alors très faibles, voire nuls. Connaissant le synchronisme important à large échelle régionale des réponses aux conditions hydroclimatiques (LE CREN, 1955; HOLMES & HENDERSON, 1990; HUGUENY et al., 1995; FOX et al., 2000; GRENOUILLET et al., 2001), cet effet négatif concerne alors un long linéaire et il s'accroît dans un système fluvial fortement artificialisé dans lequel l'homogénéité et la pauvreté de l'habitat accentue le rôle de l'hydrologie (SCHIEMER et al., 2001). A l'opposé, les conditions peuvent avoir été très favorables pour une ou plusieurs espèces et ce recrutement exceptionnel marquera les années suivantes (MILLS & MANN, 1985).

Dans un milieu fortement dégradé et eutrophe, la compétition intra- et interspécifique peut réduire les capacités de reproduction des adultes au cours des années suivantes (LAMMENS, 1986; TOWNSEND, 1989; TOWNSEND & PERROW, 1989; BRITTON et al.,
milieux aquatiques (FECHHELM et al., 1997; POWER et al., 2000; NUNN et al., 2002; RIBEIRO et al., 2004).

Pour la majorité des Cyprinidés, seules les plus jeunes classes d'âge peuvent être identifiées à partir des distributions de taille, essentiellement les classes d'âge 0+ (jeunes de l'année n), et 1+ (jeunes de l'année n-1). Généralement, un troisième groupe est distinct dans une distribution de taille (individus majoritairement d'âge 2+). Mais il présente un chevauchement considérable avec les classes d'âge proches (1+, 3+ et au-delà).

En l'absence d'une identification précise de l'âge par des techniques spécifiques appropriées (scalimétrie, otolithométrie), ce travail est basé sur une discrimination des âges à partir des histogrammes de taille.

Cette étape, coûteuse en temps, a été rendue possible par l'abondance des juvéniles dans les pêches électriques de rives, l'existence de campagnes pluriannuelles sur de nombreux sites, une relative homogénéité des dates d'échantillonnage dans les longues chroniques des sites industriels, et la complémentarité des informations recueillies par cette approche pluri scalaire, spatiale et temporelle.

Dans l'identification des jeunes de l'année, certaines études considèrent une valeur de taille « butoir » en-deçà de laquelle les individus sont considérés de l'année. Ce choix est discutable dans une étude prenant en compte d'une part l'axe rhodanien français orienté nord-sud, ses nombreux milieux, et des situations extrêmes représentées par les rejets thermiques ; et d'autre part l'usage de chroniques pluri-décennales alors que le contexte hydrologique et thermique est fortement structuré par le changement climatique actuel.

La multiplicité des situations observées n'autorise pas cette simplification extrême d'autant que l'intérêt de ce jeu de données réside dans la variabilité de la taille selon les années et les secteurs. A ce sujet, TEDESCO et al. (2009) préconisent la prise en compte des spécificités locales et notamment le gradient longitudinal, plutôt que de s'appuyer sur des valeurs globales issues de la bibliographie.

Par conséquent, chaque tableau associé à une espèce comporte l'ensemble des pêches par station et par date où l'espèce est présente.

## 2.4.2 Les problèmes de discrimination des classes d'âge

Malgré les atouts du jeu de données, le travail reste complexe du fait (1) de la forte variabilité des effectifs des juvéniles, (2) de la variabilité intra-spécifique de la croissance des juvéniles selon les années et les secteurs et (3) de la stratégie de ponte des espèces.

### 2.4.2.1 Les effectifs de juvéniles

Bien que les effectifs de juvéniles dominent dans les pêches électriques réalisées sur le Rhône, ils sont très variables dans l'espace, au cours de l'année, et d'une année à l'autre. Cette grande amplitude de variabilité a plusieurs origines.

Les conditions de reproduction peuvent ne pas avoir été favorables et les effectifs des pêches automnales sont alors très faibles, voire nuls. Connaissant le synchronisme important à large échelle régionale des réponses aux conditions hydroclimatiques (LE CREN, 1955; HOLMES & HENDERSON, 1990; HUGUENY et al., 1995; FOX et al., 2000; GRENOUILLET et al., 2001), cet effet négatif concerne alors un long linéaire et il s'accroît dans un système fluvial fortement artificialisé dans lequel l'homogénéité et la pauvreté de l'habitat accentue le rôle de l'hydrologie (SCHIEMER et al., 2001). A l'opposé, les conditions peuvent avoir été très favorables pour une ou plusieurs espèces et ce recrutement exceptionnel marquera les années suivantes (MILLS & MANN, 1985).

Dans un milieu fortement dégradé et eutrophe, la compétition intra- et interspécifique peut réduire les capacités de reproduction des adultes au cours des années suivantes (LAMMENS, 1986; TOWNSEND, 1989; TOWNSEND & PERROW, 1989; BRITTON et al.,
2009). Nous ne pouvons exclure cette hypothèse dans les années qui ont suivi le fort
recrutement du gardon en 1989, à laquelle s'ajoute la dégradation progressive de l'état
sanitaire des poissons au cours des années 1990 à 1992.

Des événements hydroclimatiques majeurs, ou des pollutions, marqués par leur
brièveté, sont survenus après la reproduction et ont éliminé les jeunes alors que les
conditions de fraie avaient été satisfaisantes (la pollution industrielle sur le Bas-Rhône en
1985, la crue d'octobre 1993 sur le Bas-Rhône, la crue majeure et les vidanges conjointes
des ouvrages hydro-électriques de la Basse-Isère en juin 2008, ...).

Les stations étudiées peuvent ne pas être propices à la présence de juvéniles du fait d'un habitat physique peu diversifié, ou de leur éloignement des sites de fraie (exemple du canal de Donzère-Mondragon et des stations situées à proximité immédiate du CNPE de Tricastin).

L'échantillonnage peut aussi être à l'origine de ces disparités. Nous rappelons que ces données n'ont pas été acquises pour étudier la taille des juvéniles. Par conséquent, si les effectifs étaient faibles voire nuls, aucun effort n'a été fourni pour obtenir une information minimale sur le niveau de recrutement et la taille potentielle de la cohorte annuelle. Par ailleurs, ces espèces sont grégaires. Les bancs, souvent mono-spécifiques, rassemblent des individus de la même taille. Lorsqu'une année est peu productive en termes de juvéniles, les maigres résultats d'une pêche peuvent provenir d'une capture qui relève du hasard : échantillonnage partiel de l'unique banc existant sur la station dans une prospection en continue de la rive, capture totale des individus présents sur un seul EPA.

Même à l'échelle de ce jeu de données, les problèmes liés aux larges fluctuations d'effectifs et déjà rencontrés par DELAIGUE et al. (2007), restent perceptibles.

2.4.2.2 La variabilité intra-spécifique de la taille

A l'échelle de quelques kilomètres, une hétérogénéité spatiale de la distribution des tailles est largement perceptible (Figure 6, Figure 7). La multiplicité des milieux situés entre deux ouvrages et interconnectés (section court-circuitée et ses annexes, canaux, retenue, affluents) sont autant d'espaces aquatiques présentant des différentiels de ressources à l'origine de taux de croissance distincts. L'existence de cohortes multiples, liées d'une part à la stratégie de reproduction et d'autre part à des reproducteurs d'âge différents pouvant utiliser ces différents milieux, exacerbe la variabilité intra-spécifique de la taille enregistrée à une date donnée. Le fonctionnement hydraulique des sections fluviales court-circuitées peut se résumer à deux états, l'un marqué par la stabilité et l'existence de conditions de reproduction favorables (débit réservé), l'autre par sa variabilité et relative hostilité lors des surverses. Le changement de « statut hydraulique » est propice à une dérive accrue des juvéniles vers l'aval. Il peut alors y avoir des chevauchements au sein des distributions de taille, consécutifs à cette dérive massive depuis ces secteurs « semi-naturels » vers les espaces aquatiques artificialisés situés en aval. Les probabilités d'obtention de distributions hétérogènes issues du « mélange » d'individus de provenances diverses augmentent au cours de la saison froide.

La représentativité de la distribution de tailles pouvant être incertaine du fait des difficultés d'échantillonner un hydrosystème de cette envergure, des méthodes et outils d'échantillonnage utilisés, de la subjectivité du découpage en classe de tailles, nous avons choisi de caractériser statistiquement chaque pêche par la **médiane** des longueurs totales des juvéniles. Celle-ci étant peu sensible aux valeurs extrêmes, elle nous apparaît comme un descripteur robuste de la distribution des tailles.
2009). Nous ne pouvons exclure cette hypothèse dans les années qui ont suivi le fort
recrutement du gardon en 1989, à laquelle s'ajoute la dégradation progressive de l'état
sanitaire des poissons au cours des années 1990 à 1992.

Des événements hydroclimatiques majeurs, ou des pollutions, marqués par leur
brièveté, sont survenus après la reproduction et ont éliminé les jeunes alors que les
conditions de fraie avaient été satisfaisantes (la pollution industrielle sur le Bas-Rhône en
1985, la crue d'octobre 1993 sur le Bas-Rhône, la crue majeure et les vidanges conjointes
des ouvrages hydro-électriques de la Basse-Isère en juin 2008, ...).

Les stations étudiées peuvent ne pas être propices à la présence de juvéniles du fait d'un habitat physique peu diversifié, ou de leur éloignement des sites de fraie (exemple du canal de Donzère-Mondragon et des stations situées à proximité immédiate du CNPE de Tricastin).

L'échantillonnage peut aussi être à l'origine de ces disparités. Nous rappelons que ces données n'ont pas été acquises pour étudier la taille des juvéniles. Par conséquent, si les effectifs étaient faibles voire nuls, aucun effort n'a été fourni pour obtenir une information minimale sur le niveau de recrutement et la taille potentielle de la cohorte annuelle. Par ailleurs, ces espèces sont grégaires. Les bancs, souvent mono-spécifiques, rassemblent des individus de la même taille. Lorsqu'une année est peu productive en termes de juvéniles, les maigres résultats d'une pêche peuvent provenir d'une capture qui relève du hasard : échantillonnage partiel de l'unique banc existant sur la station dans une prospection en continue de la rive, capture totale des individus présents sur un seul EPA.

Même à l'échelle de ce jeu de données, les problèmes liés aux larges fluctuations d'effectifs et déjà rencontrés par DELAIGUE et al. (2007), restent perceptibles.

2.4.2.2 La variabilité intra-spécifique de la taille

A l'échelle de quelques kilomètres, une hétérogénéité spatiale de la distribution des tailles est largement perceptible (Figure 6, Figure 7). La multiplicité des milieux situés entre deux ouvrages et interconnectés (section court-circuitée et ses annexes, canaux, retenue, affluents) sont autant d'espaces aquatiques présentant des différentiels de ressources à l'origine de taux de croissance distincts. L'existence de cohortes multiples, liées d'une part à la stratégie de reproduction et d'autre part à des reproducteurs d'âge différents pouvant utiliser ces différents milieux, exacerbe la variabilité intra-spécifique de la taille enregistrée à une date donnée. Le fonctionnement hydraulique des sections fluviales court-circuitées peut se résumer à deux états, l'un marqué par la stabilité et l'existence de conditions de reproduction favorables (débit réservé), l'autre par sa variabilité et relative hostilité lors des surverses. Le changement de « statut hydraulique » est propice à une dérive accrue des juvéniles vers l'aval. Il peut alors y avoir des chevauchements au sein des distributions de taille, consécutifs à cette dérive massive depuis ces secteurs « semi-naturels » vers les espaces aquatiques artificialisés situés en aval. Les probabilités d'obtention de distributions hétérogènes issues du « mélange » d'individus de provenances diverses augmentent au cours de la saison froide.

La représentativité de la distribution de tailles pouvant être incertaine du fait des difficultés d'échantillonner un hydrosystème de cette envergure, des méthodes et outils d'échantillonnage utilisés, de la subjectivité du découpage en classe de tailles, nous avons choisi de caractériser statistiquement chaque pêche par la **médiane** des longueurs totales des juvéniles. Celle-ci étant peu sensible aux valeurs extrêmes, elle nous apparaît comme un descripteur robuste de la distribution des tailles.
<figure><img src="image_7.png" /><figcaption>Figure 6 - Variabilité de la taille des gardons d'âge 0+ en septembre 2011. Résultats des pêches (PCR) de plusieurs stations proches sur le Bas-Rhône médian. MO (Montélimar, amont et aval du CNPE de Cruas-Meysse), DM (Donzère-Mondragon, amont, rejet et aval immédiat du CNPE du Tricastin), CA (Caderousse, amont et aval de la retenue). DM PSEV (Rhône court-circuité par l'aménagement de Donzère-Mondragon, aval de Pont-Saint-Esprit)</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>Figure 7 - Variabilité de la taille des chevaines d'âge 0+ en septembre 2011. Résultats des pêches (PCR) de plusieurs stations proches sur le Bas-Rhône médian. MO (Montélimar, amont et aval du CNPE de Cruas-Meysse), DM (Donzère-Mondragon, amont, rejet et aval immédiat du CNPE du Tricastin), CA (Caderousse, amont et aval de la retenue). DM PSEV (Rhône court-circuité par l'aménagement de Donzère-Mondragon, aval de Pont-Saint-Esprit)</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>Figure 6 - Variabilité de la taille des gardons d'âge 0+ en septembre 2011. Résultats des pêches (PCR) de plusieurs stations proches sur le Bas-Rhône médian. MO (Montélimar, amont et aval du CNPE de Cruas-Meysse), DM (Donzère-Mondragon, amont, rejet et aval immédiat du CNPE du Tricastin), CA (Caderousse, amont et aval de la retenue). DM PSEV (Rhône court-circuité par l'aménagement de Donzère-Mondragon, aval de Pont-Saint-Esprit)</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>Figure 7 - Variabilité de la taille des chevaines d'âge 0+ en septembre 2011. Résultats des pêches (PCR) de plusieurs stations proches sur le Bas-Rhône médian. MO (Montélimar, amont et aval du CNPE de Cruas-Meysse), DM (Donzère-Mondragon, amont, rejet et aval immédiat du CNPE du Tricastin), CA (Caderousse, amont et aval de la retenue). DM PSEV (Rhône court-circuité par l'aménagement de Donzère-Mondragon, aval de Pont-Saint-Esprit)</figcaption></figure>
### 2.4.2.3 Les stratégies de ponte

Selon la stratégie de ponte, la distribution en classes de taille des juvéniles est très contrastée (Figure 8). La résultante des pontes multiples (cas du chevaine) peut être une succession de cohortes de tailles diverses pour une même pêche, avec parfois un chevauchement des distributions de tailles entre les individus nés tardivement l'année n-1 et de ceux nés précocement l'année n.

<figure><img src="image_9.png" /></figure>

**Figure 8 - Distribution en classes de taille du gardon (en vert) et du chevaine (en bleu) montrant les effets potentiels de leur stratégie de ponte sur la forme de la distribution et son amplitude. Retenue de Montélimar le 17 septembre 2009.**

Les difficultés rencontrées avec le chevaine deviennent problématiques avec de petites espèces comme l'ablette, le goujon, la bouvière, ou le pseudorasbora.

Chez le gardon, la distribution en classe de taille est a priori idéale pour l'identification des âges. Cependant, malgré une stratégie de ponte unique, la distribution uni-modale parfaite est rarement observée à l'échelle des stations de pêche. Elle le sera encore moins lors d'un regroupement sectoriel. La distribution automnale des tailles des individus résulte d'un regroupement progressif de bancs de jeunes poissons dont les tailles sont hétérogènes car la croissance est contrastée à l'échelle locale (CARREL, 1986). Cette hétérogénéité des tailles est favorisée par la dérive active et/ou accidentelle depuis de multiples frayères (JURAJDA, 1998) et aux déplacements des jeunes poissons au cours des premiers mois de leur vie.

### 2.4.3 Vers une automatisation du tri ?

DELAIGUE et al. (2007) ont proposé une procédure de découpage automatisé en classes d'âge à partir des distributions des longueurs totales des individus.

La discrimination des individus d'âges 0+ et 1+ est réalisée à l'aide d'un modèle de mélange gaussien ajusté par maximisation de la vraisemblance. Le modèle est basé sur un mélange de deux lois normales car les pêches contiennent majoritairement des individus ayant moins de deux ans. Pour chaque année et chaque station, les valeurs seuils ont ainsi été définies pour ces deux classes d'âge.

Les limites de cette procédure sont de deux types.

Le modèle ne permet pas de traiter les distributions présentant de faibles effectifs (limite fixée arbitrairement à 30 individus). Or les années marquées par de faibles effectifs
### 2.4.2.3 Les stratégies de ponte

Selon la stratégie de ponte, la distribution en classes de taille des juvéniles est très contrastée (Figure 8). La résultante des pontes multiples (cas du chevaine) peut être une succession de cohortes de tailles diverses pour une même pêche, avec parfois un chevauchement des distributions de tailles entre les individus nés tardivement l'année n-1 et de ceux nés précocement l'année n.

<figure><img src="image_9.png" /><figcaption>Figure 8 - Distribution en classes de taille du gardon (en vert) et du chevaine (en bleu) montrant les effets potentiels de leur stratégie de ponte sur la forme de la distribution et son amplitude. Retenue de Montélimar le 17 septembre 2009.</figcaption></figure>

Les difficultés rencontrées avec le chevaine deviennent problématiques avec de petites espèces comme l'ablette, le goujon, la bouvière, ou le pseudorasbora.

Chez le gardon, la distribution en classe de taille est a priori idéale pour l'identification des âges. Cependant, malgré une stratégie de ponte unique, la distribution uni-modale parfaite est rarement observée à l'échelle des stations de pêche. Elle le sera encore moins lors d'un regroupement sectoriel. La distribution automnale des tailles des individus résulte d'un regroupement progressif de bancs de jeunes poissons dont les tailles sont hétérogènes car la croissance est contrastée à l'échelle locale (CARREL, 1986). Cette hétérogénéité des tailles est favorisée par la dérive active et/ou accidentelle depuis de multiples frayères (JURAJDA, 1998) et aux déplacements des jeunes poissons au cours des premiers mois de leur vie.

### 2.4.3 Vers une automatisation du tri ?

DELAIGUE et al. (2007) ont proposé une procédure de découpage automatisé en classes d'âge à partir des distributions des longueurs totales des individus.

La discrimination des individus d'âges 0+ et 1+ est réalisée à l'aide d'un modèle de mélange gaussien ajusté par maximisation de la vraisemblance. Le modèle est basé sur un mélange de deux lois normales car les pêches contiennent majoritairement des individus ayant moins de deux ans. Pour chaque année et chaque station, les valeurs seuils ont ainsi été définies pour ces deux classes d'âge.

Les limites de cette procédure sont de deux types.

Le modèle ne permet pas de traiter les distributions présentant de faibles effectifs (limite fixée arbitrairement à 30 individus). Or les années marquées par de faibles effectifs
présentent généralement des singularités thermiques et/ou hydrologiques. Dans ce travail, il n'est pas envisageable d'écarter ces années singulières.

La procédure nécessite un contrôle systématique a posteriori, puisque de nombreux cas particuliers ont posé problème. Les distributions de taille unimodales caractérisées par une absence d'individus d'âge 0+ n'ont pas été correctement discriminées, de même que certaines distributions bimodales atypiques (distributions très rapprochées, importance relative des juvéniles 0+ et 1+ inversée, ou poissons d'âge 0+ présentant une distribution multimodale). Les distributions à plus de deux modes, ainsi que celles comportant peu d'effectifs mais largement dispersées, posent d'autres problèmes. Il s'avère que les distributions atypiques sont nombreuses et l'expertise devient obligatoire.

Le tri des classes d'âge à partir des distributions de la taille étant long et délicat, nous avons essayé de compléter ce travail par des approches complémentaires de classification automatique. Cette question a fait l'objet d'une étude spécifique et d'un rapport détaillé rédigé par Jacques VESLOT dans le cadre d'une convention complémentaire de la phase IV.

### 2.4.4 Le choix de l'expertise

Afin de se départir des problèmes rencontrés lors des procédures automatiques, nous avons choisi la discrimination des classes d'âges par expertise.

L'expertise a été principalement faite à partir des tableaux croisés des données brutes comportant les tailles en colonne et les stations-dates en ligne (Figure 9). L'intérêt majeur du tableur Excel dans une procédure de ce type réside dans la possibilité d'utiliser des filtres en colonne permettant d'opérer des regroupements de lignes selon des critères spatiaux et/ou temporels, ou une combinaison de ces critères (dates, mois, années, stations, secteurs, tronçons, ...). Ils facilitent les nombreux allers-retours et la perception visuelle des données dans de grands tableaux.

<table><thead><tr><th>B</th><th>F</th><th>G</th><th>O</th><th>P</th><th>Q</th><th>R</th><th>S</th><th>T</th><th>U</th><th>V</th><th>W</th><th>X</th><th>Y</th><th>Z</th><th>AA</th><th>AB</th><th>AC</th><th>AD</th><th>AE</th><th>AF</th><th>AG</th><th>AH</th><th>AJ</th><th>AK</th><th>AL</th><th>AM</th><th>AN</th><th>AO</th><th>AP</th><th>AQ</th><th>AR</th><th>AS</th><th>AT</th><th>AU</th><th>AV</th><th>AW</th><th>AX</th></tr><tr><th>r_secte</th><th>mois</th><th>annee</th><th>pecl</th><th>spawn_1</th><th>Lmax_1</th><th>Lmin</th><th>spawn_2</th><th>Lmax_2</th><th>spawn_3</th><th>Lmax_3</th><th>10</th><th>14</th><th>15</th><th>16</th><th>17</th><th>18</th><th>19</th><th>20</th><th>21</th><th>22</th><th>23</th><th>24</th><th>25</th><th>26</th><th>27</th><th>28</th><th>29</th><th>30</th><th>31</th><th>32</th><th>33</th><th>34</th><th>35</th><th>36</th><th>37</th><th>38</th><th>39</th><th>40</th><th>41</th></tr></thead><tbody><tr><td>VIVI</td><td>6</td><td>1993</td><td>1993</td><td>31</td><td>0</td><td>1992</td><td>105</td><td>1991</td><td>143</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>6</td><td>1996</td><td>1996</td><td>30</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>6</td><td>2005</td><td>2005</td><td>41</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>7</td><td>2005</td><td>2005</td><td>65</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>7</td><td>2007</td><td>2007</td><td>60</td><td>0</td><td>2006</td><td>127</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>7</td><td>2008</td><td>2008</td><td>58</td><td>0</td><td>2007</td><td>106</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>8</td><td>2005</td><td>2005</td><td>86</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>1994</td><td>1994</td><td>88</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>1995</td><td>1994</td><td>126</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>1996</td><td>1996</td><td>83</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>1997</td><td>1997</td><td>80</td><td>0</td><td>1996</td><td>116</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>2000</td><td>2000</td><td>88</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>2001</td><td>2001</td><td>80</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>2002</td><td>2002</td><td>90</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>2003</td><td>2003</td><td>95</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>2004</td><td>2004</td><td>77</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>2005</td><td>2005</td><td>83</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>2006</td><td>2006</td><td>83</td><td>0</td><td>2005</td><td>133</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>2007</td><td>2007</td><td>77</td><td>0</td><td>2006</td><td>131</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>2008</td><td>2008</td><td>90</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>10</td><td>1998</td><td>1998</td><td>80</td><td>0</td><td>1997</td><td>117</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>11</td><td>1999</td><td>1998</td><td>128</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

Figure 9 - Un extrait du tableau de découpage en classes millimétriques pour la discrimination des classes d'âge des gardons sur le secteur de Viviers. Spawn_1 indique l'année de naissance de la première classe d'âge identifiée et Lmax la longueur maximale associée. L'exercice est poursuivi pour les classes d'âge suivantes, lorsque les données le permettent.

Il est possible de délimiter les tailles maximales et minimales des individus d'âge 0+, et 1+ pour la plupart des stations lorsqu'au moins un individu est présent dans chaque classe d'âge. Les pêches à très faibles effectifs initialement écartées au début de l'étude, ont été réinsérées pour aider les opérations de tri et de discrimination.

Dans une première étape, la délimitation est faite sur les données brutes, accompagnée d'une lecture en parallèle de représentations graphiques (histogrammes des distributions).
présentent généralement des singularités thermiques et/ou hydrologiques. Dans ce travail, il n'est pas envisageable d'écarter ces années singulières.

La procédure nécessite un contrôle systématique a posteriori, puisque de nombreux cas particuliers ont posé problème. Les distributions de taille unimodales caractérisées par une absence d'individus d'âge 0+ n'ont pas été correctement discriminées, de même que certaines distributions bimodales atypiques (distributions très rapprochées, importance relative des juvéniles 0+ et 1+ inversée, ou poissons d'âge 0+ présentant une distribution multimodale). Les distributions à plus de deux modes, ainsi que celles comportant peu d'effectifs mais largement dispersées, posent d'autres problèmes. Il s'avère que les distributions atypiques sont nombreuses et l'expertise devient obligatoire.

Le tri des classes d'âge à partir des distributions de la taille étant long et délicat, nous avons essayé de compléter ce travail par des approches complémentaires de classification automatique. Cette question a fait l'objet d'une étude spécifique et d'un rapport détaillé rédigé par Jacques VESLOT dans le cadre d'une convention complémentaire de la phase IV.

### 2.4.4 Le choix de l'expertise

Afin de se départir des problèmes rencontrés lors des procédures automatiques, nous avons choisi la discrimination des classes d'âges par expertise.

L'expertise a été principalement faite à partir des tableaux croisés des données brutes comportant les tailles en colonne et les stations-dates en ligne (Figure 9). L'intérêt majeur du tableur Excel dans une procédure de ce type réside dans la possibilité d'utiliser des filtres en colonne permettant d'opérer des regroupements de lignes selon des critères spatiaux et/ou temporels, ou une combinaison de ces critères (dates, mois, années, stations, secteurs, tronçons, ...). Ils facilitent les nombreux allers-retours et la perception visuelle des données dans de grands tableaux.

<table><thead><tr><th>B</th><th>F</th><th>G</th><th>O</th><th>P</th><th>Q</th><th>R</th><th>S</th><th>T</th><th>U</th><th>V</th><th>W</th><th>X</th><th>Y</th><th>Z</th><th>AA</th><th>AB</th><th>AC</th><th>AD</th><th>AE</th><th>AF</th><th>AG</th><th>AH</th><th>AJ</th><th>AK</th><th>AL</th><th>AM</th><th>AN</th><th>AO</th><th>AP</th><th>AQ</th><th>AR</th><th>AS</th><th>AT</th><th>AU</th><th>AV</th><th>AW</th><th>AX</th></tr><tr><th>r_secte</th><th>mois</th><th>année</th><th>pecl</th><th>spawn_1</th><th>Ltmx_1</th><th>Ltmx_1</th><th>spawn_2</th><th>Ltmx_2</th><th>spawn_3</th><th>Ltmx_3</th><th>10</th><th>14</th><th>15</th><th>16</th><th>17</th><th>18</th><th>19</th><th>20</th><th>21</th><th>22</th><th>23</th><th>24</th><th>25</th><th>26</th><th>27</th><th>28</th><th>29</th><th>30</th><th>31</th><th>32</th><th>33</th><th>34</th><th>35</th><th>36</th><th>37</th><th>38</th><th>39</th><th>40</th><th>41</th></tr></thead><tbody><tr><td>VIVI</td><td>6</td><td>1993</td><td>1993</td><td>31</td><td>0</td><td>1992</td><td>105</td><td>1991</td><td>143</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>6</td><td>1996</td><td>1996</td><td>30</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>6</td><td>2005</td><td>2005</td><td>41</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>7</td><td>2005</td><td>2005</td><td>65</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>7</td><td>2007</td><td>2007</td><td>60</td><td>0</td><td>2006</td><td>127</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>7</td><td>2007</td><td>2007</td><td>60</td><td>0</td><td>2007</td><td>106</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>8</td><td>2008</td><td>2008</td><td>58</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>2005</td><td>2005</td><td>86</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>1994</td><td>1994</td><td>88</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>1995</td><td>1994</td><td>126</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>1996</td><td>1996</td><td>83</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>1997</td><td>1997</td><td>80</td><td>0</td><td>1996</td><td>116</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>2000</td><td>2000</td><td>88</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>2001</td><td>2001</td><td>80</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>2002</td><td>2002</td><td>90</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>2003</td><td>2003</td><td>95</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>2004</td><td>2004</td><td>77</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>2005</td><td>2005</td><td>83</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>2006</td><td>2006</td><td>83</td><td>0</td><td>2005</td><td>133</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>2007</td><td>2007</td><td>77</td><td>0</td><td>2006</td><td>131</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>9</td><td>2008</td><td>2008</td><td>90</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>10</td><td>1998</td><td>1998</td><td>80</td><td>0</td><td>1997</td><td>117</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>VIVI</td><td>11</td><td>1999</td><td>1998</td><td>128</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

Figure 9 - Un extrait du tableau de découpage en classes millimétriques pour la discrimination des classes d'âge des gardons sur le secteur de Viviers. Spawn_1 indique l'année de naissance de la première classe d'âge identifiée et Ltmx la longueur maximale associée. L'exercice est poursuivi pour les classes d'âge suivantes, lorsque les données le permettent.

Il est possible de délimiter les tailles maximales et minimales des individus d'âge 0+, et 1+ pour la plupart des stations lorsqu'au moins un individu est présent dans chaque classe d'âge. Les pêches à très faibles effectifs initialement écartées au début de l'étude, ont été réinsérées pour aider les opérations de tri et de discrimination.

Dans une première étape, la délimitation est faite sur les données brutes, accompagnée d'une lecture en parallèle de représentations graphiques (histogrammes des distributions).
Dans une seconde étape, les résultats de l'identification des classes d'âge sont repris par des histogrammes dans lesquelles l'ajout de couleurs distinctes attribuées à chaque classe d'âge permet de visualiser et de vérifier les résultats de l'étape initiale (Figure 10).

<figure><img src="image_10.png" /><figcaption>Figure 10 - Exemples de distributions en classes de taille (5 mm) des gardons capturés dans la retenue de Montélimar (regroupement des stations amont et aval) et visualisation colorée des différentes cohortes définies lors de l'expertise des données brutes</figcaption></figure>

Les représentations graphiques ont nécessité de nombreuses adaptations afin de pouvoir utiliser au mieux les potentialités d'affichage des écrans, faciliter l'observation et la rapidité d'affichage dans des zooms avant et arrière sur des figures au format Adobe pdf en haute définition (Figure 11).

L'analyse multivariée - Analyse Factorielle des Correspondances sur les classes de taille (PERSAT & CHESSEL, 1989) – a été utilisée à plusieurs reprises pour la mise en évidence de cohortes atypiques et corriger des erreurs manifestes du découpage.

Des espèces, autres que le gardon et le chevaine, ont fait l'objet de la première étape du tri. La plupart nécessitent de rassembler des connaissances complémentaires pour le découpage en classes d'âge, par des recherches bibliographiques, voire la lecture d'écailles de spécimens récoltés sur le terrain. C'est principalement le cas pour de petites espèces
Dans une seconde étape, les résultats de l'identification des classes d'âge sont repris par des histogrammes dans lesquelles l'ajout de couleurs distinctes attribuées à chaque classe d'âge permet de visualiser et de vérifier les résultats de l'étape initiale (Figure 10).

<figure><img src="image_10.png" /><figcaption>Figure 10 - Exemples de distributions en classes de taille (5 mm) des gardons capturés dans la retenue de Montélimar (regroupement des stations amont et aval) et visualisation colorée des différentes cohortes définies lors de l'expertise des données brutes</figcaption></figure>

Les représentations graphiques ont nécessité de nombreuses adaptations afin de pouvoir utiliser au mieux les potentialités d'affichage des écrans, faciliter l'observation et la rapidité d'affichage dans des zooms avant et arrière sur des figures au format Adobe pdf en haute définition (Figure 11).

L'analyse multivariée - Analyse Factorielle des Correspondances sur les classes de taille (PERSAT & CHESSEL, 1989) – a été utilisée à plusieurs reprises pour la mise en évidence de cohortes atypiques et corriger des erreurs manifestes du découpage.

Des espèces, autres que le gardon et le chevaine, ont fait l'objet de la première étape du tri. La plupart nécessitent de rassembler des connaissances complémentaires pour le découpage en classes d'âge, par des recherches bibliographiques, voire la lecture d'écailles de spécimens récoltés sur le terrain. C'est principalement le cas pour de petites espèces
telles que la bouvière et le pseudorasbora, deux espèces pour lesquelles le chevauchement est tel qu'il nous est actuellement impossible de les délimiter sans des erreurs grossières.

L'expertise présente l'inconvénient de travailler a priori, et d'ajouter une part d'interprétation et une part de subjectivité, rendant la méthode difficilement reproductible, mais en l'absence de méthodes automatiques efficaces, elle s'avère actuellement la plus fiable. Les tableaux expertisés ont servi de base de travail pour une évaluation des méthodes de tri automatique.

<figure><img src="image_11.png" /><figcaption>Figure 11 - Extraits d'un graphique panoramique représentant la chronique complète des captures de gardons dans la station de Cruas Aval (retenue de Montélimar). Les campagnes annuelles sont regroupées par mois. Chaque classe d'âge et ses effectifs sont figurés par des « navettes » de couleur différente dont la taille est proportionnelle au nombre d'individus.</figcaption></figure>

La synthèse par année et par secteur des principaux résultats issus de la discrimination des tailles des poissons d'âge 0+ (Figure 12, Figure 13) permet de visualiser les informations disponibles pour aborder les relations potentielles associant la longueur totale des juvéniles 0+ et les paramètres hydro-climatiques majeurs. Il faut signaler (1) la relative disparité spatiale et temporelle du jeu de données malgré l'abondance de ces deux espèces, (2) des tailles globalement contrastées que ce soit au niveau spatial ou temporel lorsque de longues chroniques existent. Au niveau spatial, de forts contrastes existent entre un affluent froid (Ain) et un affluent chaud (Saône), ainsi qu'entre le Haut-Rhône et le Bas-Rhône. Sur l'axe fluvial, les secteurs sous influence d'un rejet thermique sont marqués par de plus grandes tailles (Bugey - BUGR, Péage - PEAR, Tricastin - DOMR). Sur les longues chroniques, les contrastes signalés dans la Figure 2 se retrouvent en partie dans ces tableaux croisés.
telles que la bouvière et le pseudorasbora, deux espèces pour lesquelles le chevauchement est tel qu'il nous est actuellement impossible de les délimiter sans des erreurs grossières.

L'expertise présente l'inconvénient de travailler a priori, et d'ajouter une part d'interprétation et une part de subjectivité, rendant la méthode difficilement reproductible, mais en l'absence de méthodes automatiques efficaces, elle s'avère actuellement la plus fiable. Les tableaux expertisés ont servi de base de travail pour une évaluation des méthodes de tri automatique.

<figure><img src="image_11.png" /><figcaption>Figure 11 - Extraits d'un graphique panoramique représentant la chronique complète des captures de gardons dans la station de Cruas Aval (retenue de Montélimar). Les campagnes annuelles sont regroupées par mois. Chaque classe d'âge et ses effectifs sont figurés par des « navettes » de couleur différente dont la taille est proportionnelle au nombre d'individus.</figcaption></figure>

La synthèse par année et par secteur des principaux résultats issus de la discrimination des tailles des poissons d'âge 0+ (Figure 12, Figure 13) permet de visualiser les informations disponibles pour aborder les relations potentielles associant la longueur totale des juvéniles 0+ et les paramètres hydro-climatiques majeurs. Il faut signaler (1) la relative disparité spatiale et temporelle du jeu de données malgré l'abondance de ces deux espèces, (2) des tailles globalement contrastées que ce soit au niveau spatial ou temporel lorsque de longues chroniques existent. Au niveau spatial, de forts contrastes existent entre un affluent froid (Ain) et un affluent chaud (Saône), ainsi qu'entre le Haut-Rhône et le Bas-Rhône. Sur l'axe fluvial, les secteurs sous influence d'un rejet thermique sont marqués par de plus grandes tailles (Bugey - BUGR, Péage - PEAR, Tricastin - DOMR). Sur les longues chroniques, les contrastes signalés dans la Figure 2 se retrouvent en partie dans ces tableaux croisés.
<figure><img src="image_12.png" /><figcaption>Figure 12 - Représentation synthétique de la discrimination par expertise des tailles de la base de données. Taille médiane des gardons d'âge 0+ en fin d'année, par année et par secteur.</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>Figure 12 - Représentation synthétique de la discrimination par expertise des tailles de la base de données. Taille médiane des gardons d'âge 0+ en fin d'année, par année et par secteur.</figcaption></figure>
<figure><img src="image_13.png" /><figcaption>Figure 13 - Représentation synthétique de la discrimination par expertise des tailles de la base de données. Taille médiane des chevaines d'âge 0+ en fin d'année, par année et par secteur.</figcaption></figure>
<figure><img src="image_13.png" /><figcaption>Figure 13 - Représentation synthétique de la discrimination par expertise des tailles de la base de données. Taille médiane des chevaines d'âge 0+ en fin d'année, par année et par secteur.</figcaption></figure>
## 2.5 La température

Nous avons utilisé des paramètres simples et/ou synthétiques tels que ceux définis dans des travaux antérieurs (LAUTERS & POIREL, 2001a, b; POIREL et al., 2001; CARREL et al., 2003; POIREL, 2004; BOISNEAU et al., 2008; POIREL et al., 2008). Parmi ces paramètres synthétiques, le plus usuel est la somme des degrés-jours au-dessus d'un seuil, calculée pour un intervalle de temps donné.

La sommation thermique, connue sous le nom de degrés-jours de croissance, est utilisée depuis des décennies par les praticiens agricoles et les entomologistes comme outil de prédiction de la croissance et du développement des plantes et des insectes. C'est une donnée couramment manipulée en aquaculture pour caractériser le « temps physiologique » (combinaison du temps calendaire et de la température) nécessaire pour atteindre un stade de développement (éclosion, résorption de vésicule vitelline, ...). NEUHEIMER et TAGGART (2007) recommandent son usage dans l'étude de la croissance des poissons, montrant par de nombreux exemples la performance de cette variable pour simplifier des modèles de croissance des organismes aquatiques ectothermes. Les relations linéaires entre la taille et la somme de degrés-jours (SDJ) supposent d'être dans la gamme normale de tolérance de l'espèce et de ne pas être proche des limites supérieures de la niche thermique de l'espèce (NEUHEIMER et al., 2011).

Les principales difficultés associées à l'usage de cette variable résident dans le choix de la température seuil et de la date fixant le début de la sommation thermique. Les études dans les eaux tempérées ont montré que les seuils de température pour la croissance effective des cyprinidés sont dans un intervalle de 12-15°C (MANN, 1991). La valeur de 12°C a été largement utilisée, notamment pour le gardon (MANN, 1973; MÜLLER & MENG, 1986; GRENOUILLET et al., 2001; DELAIGUE et al., 2007).

L'usage d'un seuil de 12°C s'avère discutable face à l'extrême variabilité des régimes thermiques rencontrée dans l'aire biogéographique de répartition des deux espèces. Les espèces à large aire de distribution sont capables de s'adapter aux conditions thermiques régionales afin d'assurer la viabilité de leur progéniture (KUCHARCZYK et al., 1997).

Pour le gardon, GILLET (1989) cite un seuil de 20°C pour initier la fraie du gardon dans le lac Léman. Cet auteur signale que l'amplitude thermique journalière peut excéder 5°C à cette époque de l'année, et que la ponte peut alors se dérouler à des températures inférieures (valeurs citées comprises entre 14 et 17°C). ALABASTER et LLOYD (1982) indiquent des valeurs habituellement comprises entre 8 et 19.4°C. TISSOT et SOUCHON (2010) donnent des valeurs comprises entre 10 et 18°C, avec des extrêmes de 7 et 22°C. En 1976, les gardons de la Jihlava (PEŇÁZ et al., 1978) ont pondu deux fois, les 9-10 mai (température moyenne de 15.4°C) et 25-27 mai (mêmes conditions thermiques). Dans une étude estonienne, relative aux effets du changement climatique sur la variabilité des dates de ponte du gardon et de la brème commune, NOGES et JÄRVET (2005) indiquent des températures très basses pendant la fraie, avec des médianes de 8.7°C pour les lacs et de 7.8°C pour les rivières (1951-1990). Selon les auteurs, ces mesures ne sont pas représentatives de la réalité thermique des lieux de ponte. Ces travaux portant sur une longue chronique (40 années) ont montré que malgré la hausse des températures printanières, il n'y avait pas d'évolution significative des dates de ponte du gardon. L'espèce pond en moyenne vers le 30 avril (jour 120) dans les cours d'eau et le 5 mai dans les lacs, si ce n'est que la température moyenne des eaux lors de la fraie s'est accrue de 3 à 3.5°C au cours de la période considérée. La moitié des relevés s'établissent sur 16 jours, entre le 26 avril et la 11 mai.

Les pontes du chevaine sont plus tardives. Leur multiplicité étale la reproduction sur une longue période. NUNN et al. (2007) ont signalé des pontes multiples entre la mi-juin et la mi-juillet dans des cours d'eau anglo-saxons. Dans la rivière Stour (Grande-Bretagne), MANN (1976) a observé les pontes fin mai et courant juin. Dans la rivière Jihlava
## 2.5 La température

Nous avons utilisé des paramètres simples et/ou synthétiques tels que ceux définis dans des travaux antérieurs (LAUTERS & POIREL, 2001a, b; POIREL et al., 2001; CARREL et al., 2003; POIREL, 2004; BOISNEAU et al., 2008; POIREL et al., 2008). Parmi ces paramètres synthétiques, le plus usuel est la somme des degrés-jours au-dessus d'un seuil, calculée pour un intervalle de temps donné.

La sommation thermique, connue sous le nom de degrés-jours de croissance, est utilisée depuis des décennies par les praticiens agricoles et les entomologistes comme outil de prédiction de la croissance et du développement des plantes et des insectes. C'est une donnée couramment manipulée en aquaculture pour caractériser le « temps physiologique » (combinaison du temps calendaire et de la température) nécessaire pour atteindre un stade de développement (éclosion, résorption de vésicule vitelline, ...). NEUHEIMER et TAGGART (2007) recommandent son usage dans l'étude de la croissance des poissons, montrant par de nombreux exemples la performance de cette variable pour simplifier des modèles de croissance des organismes aquatiques ectothermes. Les relations linéaires entre la taille et la somme de degrés-jours (SDJ) supposent d'être dans la gamme normale de tolérance de l'espèce et de ne pas être proche des limites supérieures de la niche thermique de l'espèce (NEUHEIMER et al., 2011).

Les principales difficultés associées à l'usage de cette variable résident dans le choix de la température seuil et de la date fixant le début de la sommation thermique. Les études dans les eaux tempérées ont montré que les seuils de température pour la croissance effective des cyprinidés sont dans un intervalle de 12-15°C (MANN, 1991). La valeur de 12°C a été largement utilisée, notamment pour le gardon (MANN, 1973; MÜLLER & MENG, 1986; GRENOUILLET et al., 2001; DELAIGUE et al., 2007).

L'usage d'un seuil de 12°C s'avère discutable face à l'extrême variabilité des régimes thermiques rencontrée dans l'aire biogéographique de répartition des deux espèces. Les espèces à large aire de distribution sont capables de s'adapter aux conditions thermiques régionales afin d'assurer la viabilité de leur progéniture (KUCHARCZYK et al., 1997).

Pour le gardon, GILLET (1989) cite un seuil de 20°C pour initier la fraie du gardon dans le lac Léman. Cet auteur signale que l'amplitude thermique journalière peut excéder 5°C à cette époque de l'année, et que la ponte peut alors se dérouler à des températures inférieures (valeurs citées comprises entre 14 et 17°C). ALABASTER et LLOYD (1982) indiquent des valeurs habituellement comprises entre 8 et 19.4°C. TISSOT et SOUCHON (2010) donnent des valeurs comprises entre 10 et 18°C, avec des extrêmes de 7 et 22°C. En 1976, les gardons de la Jihlava (PEŇÁZ et al., 1978) ont pondu deux fois, les 9-10 mai (température moyenne de 15.4°C) et 25-27 mai (mêmes conditions thermiques). Dans une étude estonienne, relative aux effets du changement climatique sur la variabilité des dates de ponte du gardon et de la brème commune, NOGES et JÄRVET (2005) indiquent des températures très basses pendant la fraie, avec des médianes de 8.7°C pour les lacs et de 7.8°C pour les rivières (1951-1990). Selon les auteurs, ces mesures ne sont pas représentatives de la réalité thermique des lieux de ponte. Ces travaux portant sur une longue chronique (40 années) ont montré que malgré la hausse des températures printanières, il n'y avait pas d'évolution significative des dates de ponte du gardon. L'espèce pond en moyenne vers le 30 avril (jour 120) dans les cours d'eau et le 5 mai dans les lacs, si ce n'est que la température moyenne des eaux lors de la fraie s'est accrue de 3 à 3.5°C au cours de la période considérée. La moitié des relevés s'établissent sur 16 jours, entre le 26 avril et la 11 mai.

Les pontes du chevaine sont plus tardives. Leur multiplicité étale la reproduction sur une longue période. NUNN et al. (2007) ont signalé des pontes multiples entre la mi-juin et la mi-juillet dans des cours d'eau anglo-saxons. Dans la rivière Stour (Grande-Bretagne), MANN (1976) a observé les pontes fin mai et courant juin. Dans la rivière Jihlava
(République Tchèque), en 1976, PENAZ et al. (1978) ont observé trois pontes successives réparties entre le 10 mai et le 5 juin 1976, avec des températures moyennes de 15.8, 15.4 et 15.1°C. Les auteurs indiquent que la première ponte a été initiée lorsque le premier pic thermique a dépassé 15°C. La ponte s'est interrompue à chaque chute de la température, reprenant dès que ce seuil était atteint de nouveau. Dans la rivière Spree (Allemagne), les chevaines ont migré deux fois de suite sur les mêmes lieux de ponte en 1995 et 1996, respectivement les 21-22 mai et 17-18 juin. En 1995, la première ponte a eu lieu lors d'une hausse rapide de la température, de 13 à 21°C. La seconde s'est déroulée alors que la température montait de 16 à 19°C (FREDRICH et al., 2003).

Sur l'axe rhodanien, nous n'avons pas d'informations quant aux dates de ponte des deux espèces, hormis une estimation faite en 1983 à partir des courbes de croissance des larves et juvéniles de deux annexes fluviales du Haut-Rhône (CARREL, 1986). Plusieurs cohortes ont été identifiées au cours du suivi : cinq pour le gardon, six pour le chevaine. Les pontes ont été estimées entre le 23 avril et le 7 juillet pour le gardon, entre le 7 et le 25 juin pour le chevaine. La première ponte du gardon se serait déroulée avec des températures plutôt basses (9°C) ; celle du chevaine lors d'une hausse thermique de 16 à 19°C, en accord avec les travaux cités précédemment.

Ces quelques informations nous renseignent indirectement sur la température approximative et la période initiale du développement embryonnaire. Ces deux informations pourraient être une base du calcul des degrés-jours dans la mesure où nous aurions à la fois le seuil et la date de début de la sommation. Mais ces données sont trop limitées au regard de la taille du bassin versant.

En l'absence de ces données, il nous fallait établir des valeurs arbitraires et réalistes vis-à-vis des espèces étudiées et permettant de vérifier notre hypothèse, à savoir qu'il existerait une relation linéaire entre le cumul thermique (SDJ) et taille (LT) en fin d'année des jeunes poissons, alors que la période de croissance saisonnière est terminée. Un ajustement logistique sur l'ensemble des pêches disponibles montre que la croissance initialement rapide des juvéniles 0+ diminue progressivement en fin de saison chaude, puis se stabilise vers une valeur maximale entre l'automne et le printemps suivant (Figure 14). Nous utilisons des pêches réalisées entre la fin août et la fin de l'année pour tester cette hypothèse en considérant que la taille maximale des individus 0+ est atteinte alors que la courbe thermique annuelle est décroissante (Figure 15). La majorité des données utilisées proviennent de pêches effectuées en septembre (Figure 16).

La taille des juvéniles en fin d'année est intimement liée au cumul d'énergie reçue au cours de la période de croissance (MANN, 1991). La durée réelle de cette période nous est inconnue, et elle est variable (1) d'une année à l'autre selon la précité et/ou la tardiveté thermique, et (2) selon la localisation spatiale des captures sur le gradient latitudinal et le type de milieu.

Pour les deux espèces, nous avons défini un seuil thermique inférieur pour le calcul de degrés-jours (DJ), et deux dates fixes entre lesquelles nous avons calculé le cumul des températures (SDJ).
(République Tchèque), en 1976, PENAZ et al. (1978) ont observé trois pontes successives réparties entre le 10 mai et le 5 juin 1976, avec des températures moyennes de 15.8, 15.4 et 15.1°C. Les auteurs indiquent que la première ponte a été initiée lorsque le premier pic thermique a dépassé 15°C. La ponte s'est interrompue à chaque chute de la température, reprenant dès que ce seuil était atteint de nouveau. Dans la rivière Spree (Allemagne), les chevaines ont migré deux fois de suite sur les mêmes lieux de ponte en 1995 et 1996, respectivement les 21-22 mai et 17-18 juin. En 1995, la première ponte a eu lieu lors d'une hausse rapide de la température, de 13 à 21°C. La seconde s'est déroulée alors que la température montait de 16 à 19°C (FREDRICH et al., 2003).

Sur l'axe rhodanien, nous n'avons pas d'informations quant aux dates de ponte des deux espèces, hormis une estimation faite en 1983 à partir des courbes de croissance des larves et juvéniles de deux annexes fluviales du Haut-Rhône (CARREL, 1986). Plusieurs cohortes ont été identifiées au cours du suivi : cinq pour le gardon, six pour le chevaine. Les pontes ont été estimées entre le 23 avril et le 7 juillet pour le gardon, entre le 7 et le 25 juin pour le chevaine. La première ponte du gardon se serait déroulée avec des températures plutôt basses (9°C) ; celle du chevaine lors d'une hausse thermique de 16 à 19°C, en accord avec les travaux cités précédemment.

Ces quelques informations nous renseignent indirectement sur la température approximative et la période initiale du développement embryonnaire. Ces deux informations pourraient être une base du calcul des degrés-jours dans la mesure où nous aurions à la fois le seuil et la date de début de la sommation. Mais ces données sont trop limitées au regard de la taille du bassin versant.

En l'absence de ces données, il nous fallait établir des valeurs arbitraires et réalistes vis-à-vis des espèces étudiées et permettant de vérifier notre hypothèse, à savoir qu'il existerait une relation linéaire entre le cumul thermique (SDJ) et taille (LT) en fin d'année des jeunes poissons, alors que la période de croissance saisonnière est terminée. Un ajustement logistique sur l'ensemble des pêches disponibles montre que la croissance initialement rapide des juvéniles 0+ diminue progressivement en fin de saison chaude, puis se stabilise vers une valeur maximale entre l'automne et le printemps suivant (Figure 14). Nous utilisons des pêches réalisées entre la fin août et la fin de l'année pour tester cette hypothèse en considérant que la taille maximale des individus 0+ est atteinte alors que la courbe thermique annuelle est décroissante (Figure 15). La majorité des données utilisées proviennent de pêches effectuées en septembre (Figure 16).

La taille des juvéniles en fin d'année est intimement liée au cumul d'énergie reçue au cours de la période de croissance (MANN, 1991). La durée réelle de cette période nous est inconnue, et elle est variable (1) d'une année à l'autre selon la précocité et/ou la tardiveté thermique, et (2) selon la localisation spatiale des captures sur le gradient latitudinal et le type de milieu.

Pour les deux espèces, nous avons défini un seuil thermique inférieur pour le calcul de degrés-jours (DJ), et deux dates fixes entre lesquelles nous avons calculé le cumul des températures (SDJ).
<figure><img src="image_14.png" /></figure>

**Figure 14 -** Taille médiane des chevaines et des gardons juvéniles d'âge 0+ en fonction du jour de pêche. Tous les individus pêchés depuis le début de la chronique sont représentés (points gris). En rouge et vert, les tailles des juvéniles de chaque pêche. Le modèle logistique (courbe orange) a été ajusté sur les tailles médianes des pêches de l'année civile en cours (points verts).

<figure><img src="image_15.png" /></figure>

**Figure 15 -** Profils thermiques du Rhône à l'amont du CNPE du Tricastin. Ensemble de la chronique et quelques années remarquables. Les lissages des valeurs moyennes journalières ont été obtenus pas une régression locale (loess) sur 73 jours.

<figure><img src="image_16.png" /></figure>
<figure><img src="image_14.png" /></figure>

**Figure 14 -** Taille médiane des chevaines et des gardons juvéniles d'âge 0+ en fonction du jour de pêche. Tous les individus pêchés depuis le début de la chronique sont représentés (points gris). En rouge et vert, les tailles des juvéniles de chaque pêche. Le modèle logistique (courbe orange) a été ajusté sur les tailles médianes des pêches de l'année civile en cours (points verts).

<figure><img src="image_15.png" /></figure>

**Figure 15 -** Profils thermiques du Rhône à l'amont du CNPE du Tricastin. Ensemble de la chronique et quelques années remarquables. Les lissages des valeurs moyennes journalières ont été obtenus pas une régression locale (loess) sur 73 jours.

<figure><img src="image_16.png" /></figure>
Figure 16 - Nombre total de pêches par mois dans le jeu de données (en gris). Nombre de pêches contenant des gardons (en vert) et des chevaines (en bleu).

## 2.5.1 Les seuils thermiques

Face à la variabilité des températures de ponte et de croissance des cyprinidés, leur capacité d'adaptation au contexte thermique régional et à des périodes favorables de courte durée (PIVNICKA, 1983; MILLS, 1988; NOGES & JARVET, 2005), le seuil de 12°C a été considéré comme trop élevé. Nous ne pouvons pas exclure la possibilité de pontes très précoces suite à des conditions favorables conjoncturelles, et ignorer l'apport énergétique par le choix d'une valeur seuil élevée.

Nous avons donc défini une valeur seuil nettement plus basse (8.8°C), correspondant à la moyenne des valeurs journalières à l'équinoxe du printemps (20 mars), et calculée sur l'ensemble des chroniques thermiques à notre disposition (données eDF-DTG). Ce seuil a été appliqué dans le calcul des SDJ pour les deux espèces bien qu'il puisse être considéré comme bas pour le chevaine.

Nous n'avons pas fixé de seuil limite supérieur, considérant que les juvéniles des espèces étudiés ne se trouvaient pas en situation de stress thermique au plus fort de l'été (ALABASTER & LLOYD, 1982).

## 2.5.2 Les dates pour le calcul des SDJ

Le calcul des SDJ a été réalisé depuis l'équinoxe de printemps (20 mars), jusqu'au jour moyen du maximum thermique<fnref n="3" /> (Cmax) calculé pour le bassin versant du Rhône (date du 8 août).

La période de croissance est donc définie entre l'équinoxe de printemps et le 8 août, soit entre les jours juliens 80 et 220.

## 2.5.3 Le calcul des SDJ.

Les SDJ ont été calculées par année-secteur. Chaque secteur regroupe un ensemble de stations de pêches supposées avoir le même régime thermique ; régime caractérisé par la station d'enregistrement de la température la plus proche.

La SDJ est calculée par :

$$ \sum_{j=80}^{220} (T_{syj} - 8,8^{\circ}\text{C}) $$

Avec *j* jour julien compris entre 80 et 220 de l'année *y* dans le secteur *s* où la température *T* est supérieure au seuil de 8.8°C (Figure 17).

Une valeur unique de SDJ est attribuée à toutes les pêches d'un même secteur pour une même année (Figure 17).

<footnote n="3">³ Les courbes de températures des chroniques complètes ont été lissées par régression de Lowess sur 73 jours. Le 8 août correspond à la moyenne du jour du maximum thermique (Cmax) de toutes ces courbes lissées.</footnote>
Figure 16 - Nombre total de pêches par mois dans le jeu de données (en gris). Nombre de pêches contenant des gardons (en vert) et des chevaines (en bleu).

## 2.5.1 Les seuils thermiques

Face à la variabilité des températures de ponte et de croissance des cyprinidés, leur capacité d'adaptation au contexte thermique régional et à des périodes favorables de courte durée (PIVNICKA, 1983; MILLS, 1988; NOGES & JARVET, 2005), le seuil de 12°C a été considéré comme trop élevé. Nous ne pouvons pas exclure la possibilité de pontes très précoces suite à des conditions favorables conjoncturelles, et ignorer l'apport énergétique par le choix d'une valeur seuil élevée.

Nous avons donc défini une valeur seuil nettement plus basse (8.8°C), correspondant à la moyenne des valeurs journalières à l'équinoxe du printemps (20 mars), et calculée sur l'ensemble des chroniques thermiques à notre disposition (données eDF-DTG). Ce seuil a été appliqué dans le calcul des SDJ pour les deux espèces bien qu'il puisse être considéré comme bas pour le chevaine.

Nous n'avons pas fixé de seuil limite supérieur, considérant que les juvéniles des espèces étudiés ne se trouvaient pas en situation de stress thermique au plus fort de l'été (ALABASTER & LLOYD, 1982).

## 2.5.2 Les dates pour le calcul des SDJ

Le calcul des SDJ a été réalisé depuis l'équinoxe de printemps (20 mars), jusqu'au jour moyen du maximum thermique<fnref n="3" /> (Cmax) calculé pour le bassin versant du Rhône (date du 8 août).

La période de croissance est donc définie entre l'équinoxe de printemps et le 8 août, soit entre les jours juliens 80 et 220.

## 2.5.3 Le calcul des SDJ.

Les SDJ ont été calculées par année-secteur. Chaque secteur regroupe un ensemble de stations de pêches supposées avoir le même régime thermique ; régime caractérisé par la station d'enregistrement de la température la plus proche.

La SDJ est calculée par :

$$ \sum_{j=80}^{220} (T_{syj} - 8,8^{\circ}\text{C}) $$

Avec *j* jour julien compris entre 80 et 220 de l'année *y* dans le secteur *s* où la température *T* est supérieure au seuil de 8.8°C (Figure 17).

Une valeur unique de SDJ est attribuée à toutes les pêches d'un même secteur pour une même année (Figure 17).

<footnote n="3">³ Les courbes de températures des chroniques complètes ont été lissées par régression de Lowess sur 73 jours. Le 8 août correspond à la moyenne du jour du maximum thermique (Cmax) de toutes ces courbes lissées.</footnote>
<figure><img src="image_17.png" /><figcaption>Température (°C) de l'année 2011 à Tricastin amont</figcaption></figure>

Figure 17 - Illustration du calcul de la sommation thermique (SDJ) en 2011 à partir des relevés de la station eDF-DTG située à l'amont du CNPE du Tricastin. La courbe représente les moyennes journalières des températures. La SDJ est égale au cumul des écarts thermiques par rapport au seuil inférieur 8.8°C (température moyenne du Rhône à l'équinoxe de mars) entre les jours juliens 80 et 220.
<figure><img src="image_17.png" /><figcaption>Température (°C) de l'année 2011 à Tricastin amont</figcaption></figure>

Figure 17 - Illustration du calcul de la sommation thermique (SDJ) en 2011 à partir des relevés de la station eDF-DTG située à l'amont du CNPE du Tricastin. La courbe représente les moyennes journalières des températures. La SDJ est égale au cumul des écarts thermiques par rapport au seuil inférieur 8.8°C (température moyenne du Rhône à l'équinoxe de mars) entre les jours juliens 80 et 220.
## 2.6 L'hydrologie

### 2.6.1 Généralités

L'influence des débits sur la croissance des juvéniles est difficile à évaluer dans la mesure où le débit et la température de l'eau sont des paramètres intimement liés et négativement corrélés à certaines périodes de l'année. Au cours de la période de reproduction et de croissance des poissons, des débits forts associés à la fonte des neiges et/ou des évènements pluvieux se traduisent par des températures basses, et inversement des débits faibles résultant d'un déficit hydrique et d'un ensoleillement important favorisent la hausse des températures.

L'hydrologie a des effets essentiellement indirects sur la taille des jeunes poissons en augmentant ou en réduisant le volume des nurseries, des espaces lentiques ou stagnants propices à la croissance des jeunes poissons (GARNER, 1996), et en modifiant les teneurs de matières en suspension donc de la turbidité des eaux. La transparence permet le bon développement des macrophytes et des herbiers qui servent d'abris et de zones préférentielles de productions phyto- et zooplanctoniques (GRENOUILLET et al., 2000; GRENOUILLET et al., 2001). Elle facilite la sélectivité et la capturabilité des proies par les poissons.

La variabilité hydrologique est à l'origine d'une dispersion plus grande des individus en augmentant la dérive accidentelle (PAVLOV et al., 2008). Des charges minérales élevées augmentent la mortalité des jeunes stades (REICHARD et al., 2002; REICHARD & JURAJDA, 2004). Les chasses sédimentaires peuvent être à l'origine d'une élimination totales des cohortes de juvéniles de l'année (ROUX, 1984).

Le débit est le facteur abiotique majeur qui contrôle le succès du recrutement d'une classe d'âge (Year-Class Strength, YCS) en agissant directement sur la mortalité par une action mécanique (CATTANEO et al., 2001), ou indirectement via une réduction de la croissance par des températures plus basses, une augmentation des dépenses énergétiques ou une disponibilité alimentaire réduite (NUNN et al., 2003).

### 2.6.2 Choix d'une métrique pour la période de croissance

L'hydrologie est vraisemblablement le facteur abiotique le plus complexe à appréhender, que les hydrologues analysent par un grand nombre de descripteurs plus ou moins redondants (OLDEN & POFF, 2003) et généralement peu satisfaisants pour nos préoccupations.

NUNN et al. (2003) ont opté pour une variable synthétique sensiblement équivalente au SDJ en calculant un cumul de débits journaliers (débit-jour) au-dessus d'un débit de base pour la période d'avril à septembre inclus. Le débit seuil correspond au débit moyen journalier calculé pour chaque mois sur 15 années.

Nous avons opté pour une démarche similaire. Afin d'éliminer des effets spatiaux liés à l'augmentation progressive des débits de l'amont vers l'aval, le descripteur est basé sur les centiles des débits moyens journaliers classés de chaque station limnigraphique. Il a été calculé pour la période allant du 20 mars au 8 août, en utilisant les chroniques complètes de débits disponibles dans la banque HYDRO<fnref n="4" />.

<footnote n="4">⁴ Nous avons utilisé les relevés des stations limnigraphiques suivantes : Rhône (Pougny, Lagnieu, Perrache, Ternay, Viviers, Beaucaire), Saône (Couzon, Macon), Ain (Chazey). Nous avons estimé les débits du Rhône à Caderousse en sommant les débits du fleuve à Viviers et ceux de l'Ardèche à Sauze.</footnote>
## 2.6 L'hydrologie

### 2.6.1 Généralités

L'influence des débits sur la croissance des juvéniles est difficile à évaluer dans la mesure où le débit et la température de l'eau sont des paramètres intimement liés et négativement corrélés à certaines périodes de l'année. Au cours de la période de reproduction et de croissance des poissons, des débits forts associés à la fonte des neiges et/ou des évènements pluvieux se traduisent par des températures basses, et inversement des débits faibles résultant d'un déficit hydrique et d'un ensoleillement important favorisent la hausse des températures.

L'hydrologie a des effets essentiellement indirects sur la taille des jeunes poissons en augmentant ou en réduisant le volume des nurseries, des espaces lentiques ou stagnants propices à la croissance des jeunes poissons (GARNER, 1996), et en modifiant les teneurs de matières en suspension donc de la turbidité des eaux. La transparence permet le bon développement des macrophytes et des herbiers qui servent d'abris et de zones préférentielles de productions phyto- et zooplanctoniques (GRENOUILLET et al., 2000; GRENOUILLET et al., 2001). Elle facilite la sélectivité et la capturabilité des proies par les poissons.

La variabilité hydrologique est à l'origine d'une dispersion plus grande des individus en augmentant la dérive accidentelle (PAVLOV et al., 2008). Des charges minérales élevées augmentent la mortalité des jeunes stades (REICHARD et al., 2002; REICHARD & JURAJDA, 2004). Les chasses sédimentaires peuvent être à l'origine d'une élimination totales des cohortes de juvéniles de l'année (ROUX, 1984).

Le débit est le facteur abiotique majeur qui contrôle le succès du recrutement d'une classe d'âge (Year-Class Strength, YCS) en agissant directement sur la mortalité par une action mécanique (CATTANEO et al., 2001), ou indirectement via une réduction de la croissance par des températures plus basses, une augmentation des dépenses énergétiques ou une disponibilité alimentaire réduite (NUNN et al., 2003).

### 2.6.2 Choix d'une métrique pour la période de croissance

L'hydrologie est vraisemblablement le facteur abiotique le plus complexe à appréhender, que les hydrologues analysent par un grand nombre de descripteurs plus ou moins redondants (OLDEN & POFF, 2003) et généralement peu satisfaisants pour nos préoccupations.

NUNN et al. (2003) ont opté pour une variable synthétique sensiblement équivalente au SDJ en calculant un cumul de débits journaliers (débit-jour) au-dessus d'un débit de base pour la période d'avril à septembre inclus. Le débit seuil correspond au débit moyen journalier calculé pour chaque mois sur 15 années.

Nous avons opté pour une démarche similaire. Afin d'éliminer des effets spatiaux liés à l'augmentation progressive des débits de l'amont vers l'aval, le descripteur est basé sur les centiles des débits moyens journaliers classés de chaque station limnigraphique. Il a été calculé pour la période allant du 20 mars au 8 août, en utilisant les chroniques complètes de débits disponibles dans la banque HYDRO<fnref n="4" />.

<footnote n="4">⁴ Nous avons utilisé les relevés des stations limnigraphiques suivantes : Rhône (Pougny, Lagnieu, Perrache, Ternay, Viviers, Beaucaire), Saône (Couzon, Macon), Ain (Chazey). Nous avons estimé les débits du Rhône à Caderousse en sommant les débits du fleuve à Viviers et ceux de l'Ardèche à Sauze.</footnote>
L'indice de débit est calculé par année et par station limnigraphique. Chaque jour, le débit moyen est remplacé par le score centré-réduit du centile auquel il appartient. La même opération est réalisée sur l'ensemble des jours de la période, puis les scores journaliers sont sommés pour obtenir l'indice de débit caractéristique de l'année et de la station (Figure 18).

Un indice positif indique des débits plutôt élevés entre le 20 mars et le 8 août par rapport à l'ensemble de la chronique et inversement pour un indice négatif. L'hypothèse sous-jacente est qu'une période caractérisée par de forts débits serait défavorable à la croissance des juvéniles et qu'une période d'étiage serait plutôt favorable.

<figure><img src="image_18.png" /></figure>

Figure 18 - Visualisation du calcul de l'indice hydrologique, descripteur de la période de croissance. Représentation de deux années contrastées au cours de la dernière décennie (2001 et 2011). En orange, la somme des degrés-jours (SDJ) pour la période de croissance de l'équinoxe de printemps à la date du maximum thermique moyen sur le Rhône. Le débit moyen journalier du Rhône à Viviers est représenté dans la partie inférieure de chaque cadre (tracé bleu) ainsi que le tracé du lissage du Quantile 90% (pointillé). En rouge (+) et en vert (-), les valeurs journalières de l'indice de débit dont le cumul donne l'indice global de la période considérée. Un tracé et une trame bleu clair indique le dépassement du quantile (cas de l'année 2001).

### 2.6.3 Les caractéristiques stationnelles

Afin d'évaluer l'influence des caractéristiques stationnelles sur la taille des jeunes poissons, mise en évidence par GRENOUILLET et al. (2001), nous avons proposé des métriques descriptives des stations. Ces métriques ont été construites avec deux contraintes : (1) représenter au mieux la gamme de milieux existants dans le jeu de données, (2) être les plus simples et objectives possible. Le Tableau X présente les métriques utilisées.
L'indice de débit est calculé par année et par station limnigraphique. Chaque jour, le débit moyen est remplacé par le score centré-réduit du centile auquel il appartient. La même opération est réalisée sur l'ensemble des jours de la période, puis les scores journaliers sont sommés pour obtenir l'indice de débit caractéristique de l'année et de la station (Figure 18).

Un indice positif indique des débits plutôt élevés entre le 20 mars et le 8 août par rapport à l'ensemble de la chronique et inversement pour un indice négatif. L'hypothèse sous-jacente est qu'une période caractérisée par de forts débits serait défavorable à la croissance des juvéniles et qu'une période d'étiage serait plutôt favorable.

<figure><img src="image_18.png" /></figure>

Figure 18 - Visualisation du calcul de l'indice hydrologique, descripteur de la période de croissance. Représentation de deux années contrastées au cours de la dernière décennie (2001 et 2011). En orange, la somme des degrés-jours (SDJ) pour la période de croissance de l'équinoxe de printemps à la date du maximum thermique moyen sur le Rhône. Le débit moyen journalier du Rhône à Viviers est représenté dans la partie inférieure de chaque cadre (tracé bleu) ainsi que le tracé du lissage du Quantile 90% (pointillé). En rouge (+) et en vert (-), les valeurs journalières de l'indice de débit dont le cumul donne l'indice global de la période considérée. Un tracé et une trame bleu clair indique le dépassement du quantile (cas de l'année 2001).

### 2.6.3 Les caractéristiques stationnelles

Afin d'évaluer l'influence des caractéristiques stationnelles sur la taille des jeunes poissons, mise en évidence par GRENOUILLET et al. (2001), nous avons proposé des métriques descriptives des stations. Ces métriques ont été construites avec deux contraintes : (1) représenter au mieux la gamme de milieux existants dans le jeu de données, (2) être les plus simples et objectives possible. Le Tableau X présente les métriques utilisées.
<table><caption>Tableau X - Métriques qualitatives utilisées pour caractériser les stations.</caption><tr><td>Situation</td><td>Localisation géographique par tronçons</td></tr><tr><td>Affluents</td><td>Situation par rapport aux principaux affluents</td></tr><tr><td>Vitesse</td><td>Indice global de vitesse du courant, de nulle à faible (0) à forte (3)</td></tr><tr><td>Macrophytes</td><td>Emprise moyenne des macrophytes, de nulle (0) à forte(3)</td></tr><tr><td>Type de berge</td><td>bétonnée, enrochement, enrochement + naturel(faible, fort, mixte), naturel (faible, fort, mixte), naturel à pente faible, naturel à pente forte</td></tr><tr><td>Rejets thermiques</td><td>0 - absence, 1 - dans le panache thermique, 2 - dans le rejet</td></tr><tr><td>Largeur du chenal</td><td>0 à 50 m, 50-150, 150-300, > à 300 m</td></tr><tr><td>Annexes fluviales</td><td>Niveau d'isolement par rapport au chenal principal (0, 1, 2)</td></tr><tr><td>Complexité hydraulique</td><td>Trois niveaux croissants</td></tr></table>

Les métriques des stations de pêche ont été renseignées par les collaborateurs du projet, en demandant un avis d'expert. Certains descripteurs n'ont pas été utilisés. D'autre part, le descripteur "type de berge" demanderait à être amélioré, car il s'est avéré difficile à renseigner et à interpréter.

Les tests utilisés sont non paramétriques : le test de Kruskal-Wallis pour l'égalité des moyennes entre les groupes, et le test post-hoc Pairwise Wilcoxon rank sum test pour l'égalité des moyennes des groupes pris deux à deux, avec correction de Bonferroni pour les tests multiples dans le second cas.

# 3 Résultats et discussion

## 3.1 Le rôle de la température

Nous rappelons deux questions initiales de ce travail. Existe-t-il une relation générale associant la taille automnale des juvéniles et la température au cours de leur croissance printanière et/ou estivale ? Peux-t-on considérer que les tailles des jeunes de l'année ont sensiblement été modifiées par les évolutions thermiques actuelles ?

### 3.1.1 Une relation linéaire simple

Malgré une dispersion importante des tailles médianes, il est possible d'établir une relation linéaire simple et positive entre la taille des juvéniles d'âge 0+ en fin d'été et la somme des degrés-jours au cours de la période de croissance des jeunes poissons (Figure 19), totalement applicable à l'échelle d'un grand bassin versant.
<table><caption>Tableau X - Métriques qualitatives utilisées pour caractériser les stations.</caption><tr><td>Situation</td><td>Localisation géographique par tronçons</td></tr><tr><td>Affluents</td><td>Situation par rapport aux principaux affluents</td></tr><tr><td>Vitesse</td><td>Indice global de vitesse du courant, de nulle à faible (0) à forte (3)</td></tr><tr><td>Macrophytes</td><td>Emprise moyenne des macrophytes, de nulle (0) à forte(3)</td></tr><tr><td>Type de berge</td><td>bétonnée, enrochement, enrochement + naturel(faible, fort, mixte), naturel (faible, fort, mixte), naturel à pente faible, naturel à pente forte</td></tr><tr><td>Rejets thermiques</td><td>0 - absence, 1 - dans le panache thermique, 2 - dans le rejet</td></tr><tr><td>Largeur du chenal</td><td>0 à 50 m, 50-150, 150-300, > à 300 m</td></tr><tr><td>Annexes fluviales</td><td>Niveau d'isolement par rapport au chenal principal (0, 1, 2)</td></tr><tr><td>Complexité hydraulique</td><td>Trois niveaux croissants</td></tr></table>

Les métriques des stations de pêche ont été renseignées par les collaborateurs du projet, en demandant un avis d'expert. Certains descripteurs n'ont pas été utilisés. D'autre part, le descripteur "type de berge" demanderait à être amélioré, car il s'est avéré difficile à renseigner et à interpréter.

Les tests utilisés sont non paramétriques : le test de Kruskal-Wallis pour l'égalité des moyennes entre les groupes, et le test post-hoc Pairwise Wilcoxon rank sum test pour l'égalité des moyennes des groupes pris deux à deux, avec correction de Bonferroni pour les tests multiples dans le second cas.

# 3 Résultats et discussion

## 3.1 Le rôle de la température

Nous rappelons deux questions initiales de ce travail. Existe-t-il une relation générale associant la taille automnale des juvéniles et la température au cours de leur croissance printanière et/ou estivale ? Peux-t-on considérer que les tailles des jeunes de l'année ont sensiblement été modifiées par les évolutions thermiques actuelles ?

### 3.1.1 Une relation linéaire simple

Malgré une dispersion importante des tailles médianes, il est possible d'établir une relation linéaire simple et positive entre la taille des juvéniles d'âge 0+ en fin d'été et la somme des degrés-jours au cours de la période de croissance des jeunes poissons (Figure 19), totalement applicable à l'échelle d'un grand bassin versant.
<figure><img src="image_19.png" /><figcaption>Figure 19 - Régression de la LT médiane de chaque pêche sur la SDJ pour le gardon (vert) et le chevaine (bleu). La taille des points correspond à leur poids dans le modèle (x0.5 lorsque le nombre d'individus dans la pêche est inférieur à 30, x1 lorsqu'il est supérieur à 30). Les pêches des stations sous influence d'un rejet thermique ont été retirées des analyses.</figcaption></figure>

<figure><img src="image_20.png" /><figcaption>Figure 20 - Tailles médianes (haut) et résidus de la régression des tailles médianes sur la SDJ (bas) pour le gardon (vert) et le chevaine (bleu) par années. Les points indiquent la médiane des tailles médianes par année; les traits fins indiquent les espaces inter-quantiles à 95%; les traits épais donnent les espaces inter-quantiles à 50%.</figcaption></figure>
<figure><img src="image_19.png" /><figcaption>Figure 19 - Régression de la LT médiane de chaque pêche sur la SDJ pour le gardon (vert) et le chevaine (bleu). La taille des points correspond à leur poids dans le modèle (x0.5 lorsque le nombre d'individus dans la pêche est inférieur à 30, x1 lorsqu'il est supérieur à 30). Les pêches des stations sous influence d'un rejet thermique ont été retirées des analyses.</figcaption></figure>

<figure><img src="image_20.png" /><figcaption>Figure 20 - Tailles médianes (haut) et résidus de la régression des tailles médianes sur la SDJ (bas) pour le gardon (vert) et le chevaine (bleu) par années. Les points indiquent la médiane des tailles médianes par année; les traits fins indiquent les espaces inter-quantiles à 95%; les traits épais donnent les espaces inter-quantiles à 50%.</figcaption></figure>
La SDJ explique à elle seule 37% de la variabilité des LT médianes en fin d'année pour le gardon et 35% pour le chevaine, avec une pente significativement différente de 0. Pour le gardon, la taille médiane en fin d'année est plus grande de 0,03 mm pour chaque degré-jour cumulé supplémentaire, et de 0,02 mm pour le chevaine.

Dans les rubriques précédentes, nous avons énuméré un ensemble de problèmes en partie explicatifs des résidus de la régression.

## 3.1.2 Les rejets thermiques

Les stations situées au niveau des rejets thermiques des CNPE fonctionnant en circuit ouvert présentent des spécificités, au-delà d'une simple augmentation de la température. A ce niveau, la température des eaux est élevée (+10°C) et les observations biologiques sont d'un autre ordre en regard de situations normales (Figure 21). En effet, de nombreuses espèces évitent ces secteurs car les températures sont au-dessus des limites tolérables (HADDERINGH, 1975). Chez les cyprinidés, relativement tolérants vis-à-vis des températures élevées, la croissance des juvéniles de gardon et de brème pourrait s'interrompre au-delà de 28°C (ALABASTER & LLOYD, 1982) mais le cycle saisonnier thermique étant très décalé vers les hautes valeurs par rapport à une situation amont, la durée de la période de croissance est notablement étendue. La photopériode redevient alors le facteur de contrôle (KARAS & NEUMAN, 1981). Cet effet « positif » doit être modéré quand les autres effets physiologiques sont considérés (LUKSIENE & SANDSTRÖM, 1994; LUKSIENE & SVEDÄNG, 1997).

La relation linéaire précédente a été établie en ne prenant pas en compte les stations singulières situées dans les panaches thermiques. MATTHEEUWS et al. (1981) ont montré que les gardons de la Meuse soumis aux rejets thermiques de la centrale de Tihange auraient une reproduction avancée de 3 semaines. Cette précocité de la reproduction peut expliquer aisément une taille plus grande dans des eaux chaudes si les ressources alimentaires ne sont pas limitées. Cette croissance rapide est largement mise à profit dans la pisciculture (PEŇÁZ, 1980; PONCIN et al., 1989).

<figure><img src="image_21.png" /><figcaption>Figure 21 - Cycle thermique annuel de la température des eaux à l'amont, au rejet du CNPE du Tricastin et à l'aval de l'usine hydroélectrique de Bollène. Moyennes journalières (points) et quantiles journaliers à 50% entre 1981 et 2011.</figcaption></figure>
La SDJ explique à elle seule 37% de la variabilité des LT médianes en fin d'année pour le gardon et 35% pour le chevaine, avec une pente significativement différente de 0. Pour le gardon, la taille médiane en fin d'année est plus grande de 0,03 mm pour chaque degré-jour cumulé supplémentaire, et de 0,02 mm pour le chevaine.

Dans les rubriques précédentes, nous avons énuméré un ensemble de problèmes en partie explicatifs des résidus de la régression.

### 3.1.2 Les rejets thermiques

Les stations situées au niveau des rejets thermiques des CNPE fonctionnant en circuit ouvert présentent des spécificités, au-delà d'une simple augmentation de la température. A ce niveau, la température des eaux est élevée (+10°C) et les observations biologiques sont d'un autre ordre en regard de situations normales (Figure 21). En effet, de nombreuses espèces évitent ces secteurs car les températures sont au-dessus des limites tolérables (HADDERINGH, 1975). Chez les cyprinidés, relativement tolérants vis-à-vis des températures élevées, la croissance des juvéniles de gardon et de brème pourrait s'interrompre au-delà de 28°C (ALABASTER & LLOYD, 1982) mais le cycle saisonnier thermique étant très décalé vers les hautes valeurs par rapport à une situation amont, la durée de la période de croissance est notablement étendue. La photopériode redevient alors le facteur de contrôle (KARAS & NEUMAN, 1981). Cet effet « positif » doit être modéré quand les autres effets physiologiques sont considérés (LUKSIENE & SANDSTRÖM, 1994; LUKSIENE & SVEDÄNG, 1997).

La relation linéaire précédente a été établie en ne prenant pas en compte les stations singulières situées dans les panaches thermiques. MATTHEEUWS et al. (1981) ont montré que les gardons de la Meuse soumis aux rejets thermiques de la centrale de Tihange auraient une reproduction avancée de 3 semaines. Cette précocité de la reproduction peut expliquer aisément une taille plus grande dans des eaux chaudes si les ressources alimentaires ne sont pas limitées. Cette croissance rapide est largement mise à profit dans la pisciculture (PEŇÁZ, 1980; PONCIN et al., 1989).

<figure><img src="image_21.png" /><figcaption>Figure 21 - Cycle thermique annuel de la température des eaux à l'amont, au rejet du CNPE du Tricastin et à l'aval de l'usine hydroélectrique de Bollène. Moyennes journalières (points) et quantiles journaliers à 50% entre 1981 et 2011.</figcaption></figure>
<figure><img src="image_22.png" /></figure>

Figure 22 - Taille moyenne des gardons et des chevaines d'âge 0+ en fonction de la localisation de la station par rapport aux rejets thermiques. 0 - station sans influence d'un rejet, 1 - station située dans le panache thermique en aval d'un rejet, 3 - station située au niveau même du rejet.

Bien que nous ayons à notre disposition la température du rejet de Tricastin, cette valeur n'est guère représentative du panache thermique dans lequel circulent les poissons. Par ailleurs, bien que ces deux poissons eurythermes puissent momentanément tolérer ces températures extrêmes (ALABASTER & LLOYD, 1982), ils ne séjournent pas en permanence dans les eaux réchauffées comme le montrent les premiers résultats des expérimentations de nos collègues de Irstea Lyon (BERGE et al., 2011).

Les tailles des gardons et chevaines d'âge 0+ sont significativement différentes selon les trois modalités (Figure 22). La différence de taille est cependant plus marquée chez le gardon.

Cette différence de taille peut être associée à la durée de croissance « potentielle » plus longue dans les zones réchauffées, et vraisemblablement une productivité localement plus importante de la flore et de la faune benthique (ROGER et al., 1991).

Il est important de prendre en compte les déplacements possibles des juvéniles. En effet, si la température devient trop élevée pour leur survie, ils ont la possibilité de s'éloigner du rejet thermique, et de retrouver des conditions plus favorables. L'existence des déplacements des juvéniles entre les zones réchauffées et non réchauffées ont pu être mises en évidence en milieu lacustre (MCINERNY & DEGAN, 1991).

Les pêches faites à Tricastin montrent que la taille globalement plus grande des juvéniles n'est toutefois pas en accord avec le cumul thermique enregistré dans les eaux du rejet. Nous avons calculé des SDJ à partir des données thermiques du rejet, puis à partir d'une température moyenne entre l'amont et le rejet afin d'évaluer la pertinence de la relation initiale vis-à-vis de situations extrêmes. L'exercice a été fait pour les deux espèces et il montre la robustesse de cette relation linéaire (Figure 23, Figure 24), mais aussi la nécessité d'avoir des relevés thermiques « au plus près » des habitats thermiques fréquentés par les poissons.
<figure><img src="image_22.png" /></figure>

Figure 22 - Taille moyenne des gardons et des chevaines d'âge 0+ en fonction de la localisation de la station par rapport aux rejets thermiques. 0 - station sans influence d'un rejet, 1 - station située dans le panache thermique en aval d'un rejet, 3 - station située au niveau même du rejet.

Bien que nous ayons à notre disposition la température du rejet de Tricastin, cette valeur n'est guère représentative du panache thermique dans lequel circulent les poissons. Par ailleurs, bien que ces deux poissons eurythermes puissent momentanément tolérer ces températures extrêmes (ALABASTER & LLOYD, 1982), ils ne séjournent pas en permanence dans les eaux réchauffées comme le montrent les premiers résultats des expérimentations de nos collègues de Irstea Lyon (BERGE et al., 2011).

Les tailles des gardons et chevaines d'âge 0+ sont significativement différentes selon les trois modalités (Figure 22). La différence de taille est cependant plus marquée chez le gardon.

Cette différence de taille peut être associée à la durée de croissance « potentielle » plus longue dans les zones réchauffées, et vraisemblablement une productivité localement plus importante de la flore et de la faune benthique (ROGER et al., 1991).

Il est important de prendre en compte les déplacements possibles des juvéniles. En effet, si la température devient trop élevée pour leur survie, ils ont la possibilité de s'éloigner du rejet thermique, et de retrouver des conditions plus favorables. L'existence des déplacements des juvéniles entre les zones réchauffées et non réchauffées ont pu être mises en évidence en milieu lacustre (MCINERNY & DEGAN, 1991).

Les pêches faites à Tricastin montrent que la taille globalement plus grande des juvéniles n'est toutefois pas en accord avec le cumul thermique enregistré dans les eaux du rejet. Nous avons calculé des SDJ à partir des données thermiques du rejet, puis à partir d'une température moyenne entre l'amont et le rejet afin d'évaluer la pertinence de la relation initiale vis-à-vis de situations extrêmes. L'exercice a été fait pour les deux espèces et il montre la robustesse de cette relation linéaire (Figure 23, Figure 24), mais aussi la nécessité d'avoir des relevés thermiques « au plus près » des habitats thermiques fréquentés par les poissons.
<figure><img src="image_23.png" /></figure>

Figure 23 - Régression de la LT médiane des gardons de chaque pêche sur la SDJ, et projection en données supplémentaires des pêches faites au rejet et à l'aval immédiat du CNPE du Tricastin. En vert (1), les pêches réalisées sur l'ensemble du Rhône (rejets exclus) ; en vert cerclé de noir (2) celles réalisées dans le panache thermique du CNPE avec la SDJ calculée sur la thermie amont du site ; en noir (3) les mêmes pêches avec la SDJ calculée avec la thermie du rejet ; en gris cerclé de noir (4) la SDJ est calculée sur la thermie moyenne entre l'amont et le rejet. Les régressions 2, 3 et 4 contiennent donc toutes les pêches du Rhône y compris les pêches du panache thermique du CNPE du Tricastin. Seules les températures utilisées pour le calcul de la SDJ changent entre ces régressions.

<figure><img src="image_24.png" /></figure>

Figure 24 - Régression de la LT médiane des chevaines de chaque pêche sur la SDJ, et projection en données supplémentaires des pêches faites au rejet et à l'aval immédiat du CNPE du Tricastin. En vert (1), les pêches réalisées sur l'ensemble du Rhône (rejets exclus) ; en vert cerclé de noir (2) celles réalisées dans le panache thermique du CNPE avec la SDJ calculée sur la thermie amont du site ; en noir (3) les mêmes pêches avec la SDJ calculée avec la thermie du rejet ; en gris cerclé de noir (4) la SDJ est calculée sur la thermie moyenne entre l'amont et le rejet. Les régressions 2, 3 et 4 contiennent donc toutes les pêches du Rhône y compris les pêches du panache thermique du CNPE du Tricastin. Seules les températures utilisées pour le calcul de la SDJ changent entre ces régressions.

### 3.1.3 Des évolutions temporelles marquées

Il importe de rappeler que les évolutions de la thermie du Rhône consécutives au changement climatique se traduisent par des températures estivales de plus en plus élevées et une accélération de la hausse thermique printanière. La prolongation de la saison chaude vers la fin de l'automne est un phénomène plus récent. Ces changements peuvent être traduits simplement par le suivi inter-annuel de la date de passage d'un seuil thermique, par exemple 12°C (Figure 25).
<figure><img src="image_23.png" /></figure>

**Figure 23** - Régression de la LT médiane des gardons de chaque pêche sur la SDJ, et projection en données supplémentaires des pêches faites au rejet et à l'aval immédiat du CNPE du Tricastin. En vert (1), les pêches réalisées sur l'ensemble du Rhône (rejets exclus) ; en vert cerclé de noir (2) celles réalisées dans le panache thermique du CNPE avec la SDJ calculée sur la thermie amont du site ; en noir (3) les mêmes pêches avec la SDJ calculée avec la thermie du rejet ; en gris cerclé de noir (4) la SDJ est calculée sur la thermie moyenne entre l'amont et le rejet. Les régressions 2, 3 et 4 contiennent donc toutes les pêches du Rhône y compris les pêches du panache thermique du CNPE du Tricastin. Seules les températures utilisées pour le calcul de la SDJ changent entre ces régressions.

<figure><img src="image_24.png" /></figure>

**Figure 24** - Régression de la LT médiane des chevaines de chaque pêche sur la SDJ, et projection en données supplémentaires des pêches faites au rejet et à l'aval immédiat du CNPE du Tricastin. En vert (1), les pêches réalisées sur l'ensemble du Rhône (rejets exclus) ; en vert cerclé de noir (2) celles réalisées dans le panache thermique du CNPE avec la SDJ calculée sur la thermie amont du site ; en noir (3) les mêmes pêches avec la SDJ calculée avec la thermie du rejet ; en gris cerclé de noir (4) la SDJ est calculée sur la thermie moyenne entre l'amont et le rejet. Les régressions 2, 3 et 4 contiennent donc toutes les pêches du Rhône y compris les pêches du panache thermique du CNPE du Tricastin. Seules les températures utilisées pour le calcul de la SDJ changent entre ces régressions.

### 3.1.3 Des évolutions temporelles marquées

Il importe de rappeler que les évolutions de la thermie du Rhône consécutives au changement climatique se traduisent par des températures estivales de plus en plus élevées et une accélération de la hausse thermique printanière. La prolongation de la saison chaude vers la fin de l'automne est un phénomène plus récent. Ces changements peuvent être traduits simplement par le suivi inter-annuel de la date de passage d'un seuil thermique, par exemple 12°C (Figure 25).
<figure><img src="image_25.png" /></figure>

Figure 25 – Date (jour Julien) de passage du seuil de 12°C à Cruas amont au cours de la chronique. Le C12s (vert) est le premier jour de l'année où la courbe de température, lissée par régression locale (loess sur 73 jours), est au-dessus de 12°C. Le C12a (rouge) est le dernier. Evolution de la somme des degrés-jours (SDJ) sur le même secteur.

Lorsque les années de relevés sont projetées sur la relation linéaire, il s'avère que la taille des jeunes poissons a globalement suivi la hausse thermique (Figure 26) de la dernière décennie.

<figure><img src="image_26.png" /></figure>

Figure 26 - Identification des années de pêche dans le nuage de points des LT médianes sur la SDJ pour le gardon (haut) et le chevaine (bas). Utilisation d'un dégradé de couleur allant du bleu pour les plus anciennes vers le jaune pour les plus récentes, en passant par le rouge.

Cependant, l'évolution croissante de la taille ne prend réellement effet qu'à partir des années 1990 (Figure 27, Figure 28).
<figure><img src="image_25.png" /></figure>

Figure 25 – Date (jour Julien) de passage du seuil de 12°C à Cruas amont au cours de la chronique. Le C12s (vert) est le premier jour de l'année où la courbe de température, lissée par régression locale (loess sur 73 jours), est au-dessus de 12°C. Le C12a (rouge) est le dernier. Evolution de la somme des degrés-jours (SDJ) sur le même secteur.

Lorsque les années de relevés sont projetées sur la relation linéaire, il s'avère que la taille des jeunes poissons a globalement suivi la hausse thermique (Figure 26) de la dernière décennie.

<figure><img src="image_26.png" /></figure>

Figure 26 - Identification des années de pêche dans le nuage de points des LT médianes sur la SDJ pour le gardon (haut) et le chevaine (bas). Utilisation d'un dégradé de couleur allant du bleu pour les plus anciennes vers le jaune pour les plus récentes, en passant par le rouge.

Cependant, l'évolution croissante de la taille ne prend réellement effet qu'à partir des années 1990 (Figure 27, Figure 28).
<figure><img src="image_27.png" /><figcaption>Figure 27 - Evolution sectorielle de la taille médiane des gardons au cours du temps.</figcaption></figure>

<figure><img src="image_28.png" /><figcaption>Figure 28 - Evolution sectorielle de la taille médiane des chevaines au cours du temps.</figcaption></figure>
<figure><img src="image_27.png" /><figcaption>Figure 27 - Evolution sectorielle de la taille médiane des gardons au cours du temps.</figcaption></figure>

<figure><img src="image_28.png" /><figcaption>Figure 28 - Evolution sectorielle de la taille médiane des chevaines au cours du temps.</figcaption></figure>
## 3.2 Influence du débit

La régression de la taille médiane sur l'indice de débit (Figure 30) explique 12% de la variabilité pour le gardon et 8% pour le chevaine Toutefois cet indice est très corrélé à la SDJ (Figure 29). Lorsqu'ils sont tous deux intégrés en covariables pour expliquer la LT médiane, la pente associée à l'indice de débit n'est plus significative.

La corrélation négative associant l'indice de débit et la somme des degrés-jours montre parfaitement l'influence du débit sur la température de l'eau pendant la période de printanière.

<figure><img src="image_29.png" /><figcaption>Figure 29 - Corrélation entre la SDJ et l'indice de débit pour la période de croissance. Le coefficient de corrélation de Pearson (0.55) est hautement significatif (avec une p-value inférieure à 0,001).</figcaption></figure>

<figure><img src="image_30.png" /><figcaption>Figure 30 - Régression de la LT médiane de chaque pêche sur l'indice de débit pour le gardon (vert) et le chevaine (bleu). La taille des points correspond à leur poids dans le modèle (×0.5 lorsque le nombre d'individus dans la pêche est inférieur à 30 et ×1 lorsqu'il est supérieur à 30). Les pêches faites à proximité des rejets de centrales ont été retirés des analyses.</figcaption></figure>
## 3.2 Influence du débit

La régression de la taille médiane sur l'indice de débit (Figure 30) explique 12% de la variabilité pour le gardon et 8% pour le chevaine Toutefois cet indice est très corrélé à la SDJ (Figure 29). Lorsqu'ils sont tous deux intégrés en covariables pour expliquer la LT médiane, la pente associée à l'indice de débit n'est plus significative.

La corrélation négative associant l'indice de débit et la somme des degrés-jours montre parfaitement l'influence du débit sur la température de l'eau pendant la période de printanière.

<figure><img src="image_29.png" /><figcaption>Figure 29 - Corrélation entre la SDJ et l'indice de débit pour la période de croissance. Le coefficient de corrélation de Pearson (0.55) est hautement significatif (avec une p-value inférieure à 0,001).</figcaption></figure>

<figure><img src="image_30.png" /><figcaption>Figure 30 - Régression de la LT médiane de chaque pêche sur l'indice de débit pour le gardon (vert) et le chevaine (bleu). La taille des points correspond à leur poids dans le modèle (×0.5 lorsque le nombre d'individus dans la pêche est inférieur à 30 et ×1 lorsqu'il est supérieur à 30). Les pêches faites à proximité des rejets de centrales ont été retirés des analyses.</figcaption></figure>
## 3.3 Le gradient amont-aval

A l'échelle du Rhône, dont le bassin versant est nettement orienté nord-sud, le gradient spatial est marqué. Les jeunes gardons et chevaines sont petits dans les stations du Haut-Rhône et de l'Ain. Leur taille est intermédiaire dans les stations du tiers amont du Bas-Rhône, tandis que les plus grands sont retrouvés en aval. Mais les poissons 0+ de la Saône se démarquent de tous les autres ensembles géographiques par leur grande taille, et une relative similarité avec les stations du Bas-Rhône aval, voire les rejets thermiques. Dans une étude portant sur la taille du chevaine dans le gradient longitudinal et latitudinal de l'axe Saône-Rhône, TEDESCO et al., (2009) ont montré l'accroissement de la taille sur ce gradient sans avoir cependant la donnée thermique nécessaire pour faire la part des effets respectifs liés à la température et à l'habitat.

Nos données couvrent une échelle spatiale et temporelle beaucoup plus importante sur laquelle la structuration fluviale par grands tronçons est bien définie (Figure 31).

<figure><img src="image_31.png" /></figure>

Figure 31 - Identification des tronçons fluviaux dans le nuage de points des LT médianes sur la SDJ pour le gardon (haut) et le chevaine (bas). Haut-Rhône (HR_AMT & HR_AVL), Bas-Rhône (BR_AMT, BR_MED & BR_AVL), Saône (AF_SAO) et Ain (AF_AIN).
## 3.3 Le gradient amont-aval

A l'échelle du Rhône, dont le bassin versant est nettement orienté nord-sud, le gradient spatial est marqué. Les jeunes gardons et chevaines sont petits dans les stations du Haut-Rhône et de l'Ain. Leur taille est intermédiaire dans les stations du tiers amont du Bas-Rhône, tandis que les plus grands sont retrouvés en aval. Mais les poissons 0+ de la Saône se démarquent de tous les autres ensembles géographiques par leur grande taille, et une relative similarité avec les stations du Bas-Rhône aval, voire les rejets thermiques. Dans une étude portant sur la taille du chevaine dans le gradient longitudinal et latitudinal de l'axe Saône-Rhône, TEDESCO et al., (2009) ont montré l'accroissement de la taille sur ce gradient sans avoir cependant la donnée thermique nécessaire pour faire la part des effets respectifs liés à la température et à l'habitat.

Nos données couvrent une échelle spatiale et temporelle beaucoup plus importante sur laquelle la structuration fluviale par grands tronçons est bien définie (Figure 31).

<figure><img src="image_31.png" /><figcaption>Figure 31 - Identification des tronçons fluviaux dans le nuage de points des LT médianes sur la SDJ pour le gardon (haut) et le chevaine (bas). Haut-Rhône (HR_AMT & HR_AVL), Bas-Rhône (BR_AMT, BR_MED & BR_AVL), Saône (AF_SAO) et Ain (AF_AIN).</figcaption></figure>
## 3.4 Le fleuve aménagé

L'usage par défaut de valeurs thermiques identiques pour caractériser plusieurs stations réparties dans des entités fluviales contrastées (lône, RCC, chenal « naturel », canal, retenue) est vraisemblablement à l'origine d'une grande partie de la variance non expliquée. Nous n'avions pas la possibilité de renseigner thermiquement les sections court-circuitées du fleuve et avons utilisé les relevés thermiques « chenal principal » les plus proches.

En distinguant les pêches selon leur appartenance à des entités « semi-naturelles » (anciens tracés du Rhône) ou des sections chenalisées de grandes dimensions (chenal, canal de navigation, retenue, bief sur la Saône), la relation ne change pas entre les deux grands types de milieux avec une pente non significativement différente (Figure 32, Figure 33). Il n'apparaît pas d'effet du milieu, ni d'interaction : p-values respectivement de 0.36 et 0.14 pour le gardon, et 0.06 et 0.53 pour le chevaine. Le R2 montre même un meilleur ajustement là où la thermie est censée être la plus approximative. Faut-il considérer que la thermie utilisée pour le cours principal représente aussi une grande approximation des conditions locales ?

<figure><img src="image_32.png" /></figure>

Figure 32 - Régression de la LT médiane de chaque pêche sur la SDJ pour le gardon en distinguant les milieux « semi-naturels » des sections fluviales de grandes dimensions. Les stations sous influence d'un rejet thermique ont été retirées des analyses.
## 3.4 Le fleuve aménagé

L'usage par défaut de valeurs thermiques identiques pour caractériser plusieurs stations réparties dans des entités fluviales contrastées (lône, RCC, chenal « naturel », canal, retenue) est vraisemblablement à l'origine d'une grande partie de la variance non expliquée. Nous n'avions pas la possibilité de renseigner thermiquement les sections court-circuitées du fleuve et avons utilisé les relevés thermiques « chenal principal » les plus proches.

En distinguant les pêches selon leur appartenance à des entités « semi-naturelles » (anciens tracés du Rhône) ou des sections chenalisées de grandes dimensions (chenal, canal de navigation, retenue, bief sur la Saône), la relation ne change pas entre les deux grands types de milieux avec une pente non significativement différente (Figure 32, Figure 33). Il n'apparaît pas d'effet du milieu, ni d'interaction : p-values respectivement de 0.36 et 0.14 pour le gardon, et 0.06 et 0.53 pour le chevaine. Le R2 montre même un meilleur ajustement là où la thermie est censée être la plus approximative. Faut-il considérer que la thermie utilisée pour le cours principal représente aussi une grande approximation des conditions locales ?

<figure><img src="image_32.png" /></figure>

Figure 32 - Régression de la LT médiane de chaque pêche sur la SDJ pour le gardon en distinguant les milieux « semi-naturels » des sections fluviales de grandes dimensions. Les stations sous influence d'un rejet thermique ont été retirées des analyses.
<figure><img src="image_33.png" /><figcaption>Figure 33 - Régression de la LT médiane de chaque pêche sur la SDJ pour le chevaine en distinguant les milieux « semi-naturels » des sections fluviales de grandes dimensions. Les stations sous influence d'un rejet thermique ont été retirées des analyses.</figcaption></figure>
<figure><img src="image_33.png" /><figcaption>Figure 33 - Régression de la LT médiane de chaque pêche sur la SDJ pour le chevaine en distinguant les milieux « semi-naturels » des sections fluviales de grandes dimensions. Les stations sous influence d'un rejet thermique ont été retirées des analyses.</figcaption></figure>
## 3.5 Les paramètres locaux

Afin de prendre en compte des facteurs abiotiques autres qu'hydroclimatiques, chaque station a fait l'objet d'une description qualitative simplifiée.

En regroupant les stations de pêches selon les différentes catégories définies pour ces facteurs, des différences significatives des LT médianes peuvent apparaître entre les catégories (Figure 34, Figure 36).

Ces différences s'amoindrissent après contrôle de la thermie (par la SDJ) (Figure 35, Figure 37).

Les écarts de taille entre grands secteurs se réduisent, moins cependant chez le gardon. Pour ce dernier, la Saône ainsi que le Bas-Rhône médian reste au-dessus des autres. Les caractéristiques de ces tronçons se retrouvent en partie dans d'autres descripteurs (bief pour la Saône, canal amenée pour le Bas-Rhône médian, nature des berges, la présence de macrophytes). Pour le chevaine, l'écart positif ne persiste que pour la Saône.

Rappel des catégories de la description qualitative des stations utilisées dans les figures suivantes (Figure 34, Figure 35, Figure 36, Figure 37).

- Situation (AF_AIN : Ain, AF_SAO : Saône, BR_AMT : Bas-Rhône amont, BR_AVL : Bas-Rhône aval, BR_MED : Bas-Rhône médian, HR_AMT : Haut-Rhône amont, HR_AVL : Haut-Rhône aval.

- Milieu (BIE : Bief, CAA : Canal d'amenée, CAF : Canal de fuite, CHE : Chenal, RCC : Rhône court-circuité, RET : Retenue)

- Berges (be : béton, en : enrochement, en/na-mi : enrochement/naturel mixte, en/nafa : enrochement/naturel à pente faible, na/fa : naturel à pente faible, na-mi : naturel mixte, na-pefa : naturel à pente faible, na-pefo : naturel à pente forte)

- Macrophytes (0 : absence, 1 : densité faible, 2 : densité moyenne, 3 : densité forte)

- Annexe [Niveau d'isolement par rapport au chenal principal (0, 1, 2)]

- Largeur du chenal (en mètres)

- Chenal [1 : chenal unique, 2 : chenal tressé (plusieurs chenaux sur la station)]

- Complexité hydraulique (1 : homogène, largeur fixe, 2 : semi-hétérogène, largeur variable, 3 : type "plaine d'inondation »)
## 3.5 Les paramètres locaux

Afin de prendre en compte des facteurs abiotiques autres qu'hydroclimatiques, chaque station a fait l'objet d'une description qualitative simplifiée.

En regroupant les stations de pêches selon les différentes catégories définies pour ces facteurs, des différences significatives des LT médianes peuvent apparaître entre les catégories (Figure 34, Figure 36).

Ces différences s'amoindrissent après contrôle de la thermie (par la SDJ) (Figure 35, Figure 37).

Les écarts de taille entre grands secteurs se réduisent, moins cependant chez le gardon. Pour ce dernier, la Saône ainsi que le Bas-Rhône médian reste au-dessus des autres. Les caractéristiques de ces tronçons se retrouvent en partie dans d'autres descripteurs (bief pour la Saône, canal amenée pour le Bas-Rhône médian, nature des berges, la présence de macrophytes). Pour le chevaine, l'écart positif ne persiste que pour la Saône.

Rappel des catégories de la description qualitative des stations utilisées dans les figures suivantes (Figure 34, Figure 35, Figure 36, Figure 37).

- Situation (AF_AIN : Ain, AF_SAO : Saône, BR_AMT : Bas-Rhône amont, BR_AVL : Bas-Rhône aval, BR_MED : Bas-Rhône médian, HR_AMT : Haut-Rhône amont, HR_AVL : Haut-Rhône aval.

- Milieu (BIE : Bief, CAA : Canal d'amenée, CAF : Canal de fuite, CHE : Chenal, RCC : Rhône court-circuité, RET : Retenue)

- Berges (be : béton, en : enrochement, en/na-mi : enrochement/naturel mixte, en/nafa : enrochement/naturel à pente faible, na/fa : naturel à pente faible, na-mi : naturel mixte, na-pefa : naturel à pente faible, na-pefo : naturel à pente forte)

- Macrophytes (0 : absence, 1 : densité faible, 2 : densité moyenne, 3 : densité forte)

- Annexe [Niveau d'isolement par rapport au chenal principal (0, 1, 2)]

- Largeur du chenal (en mètres)

- Chenal [1 : chenal unique, 2 : chenal tressé (plusieurs chenaux sur la station)]

- Complexité hydraulique (1 : homogène, largeur fixe, 2 : semi-hétérogène, largeur variable, 3 : type "plaine d'inondation »)
<figure><img src="image_34.png" /><figcaption>Figure 34 - Représentation des LT médianes des gardons d'âge 0+ classées selon certaines catégories propres aux stations.</figcaption></figure>
<figure><img src="image_34.png" /><figcaption>Figure 34 - Représentation des LT médianes des gardons d'âge 0+ classées selon certaines catégories propres aux stations.</figcaption></figure>
<figure><img src="image_35.png" /><figcaption>Figure 35 - Représentation des résidus de la régression sur la SDJ de la taille des gardons d'âge 0+ classés selon certaines catégories propres aux stations.</figcaption></figure>
<figure><img src="image_35.png" /><figcaption>Figure 35 - Représentation des résidus de la régression sur la SDJ de la taille des gardons d'âge 0+ classés selon certaines catégories propres aux stations.</figcaption></figure>
LTmm

situation

<figure><img src="image_36.png" /></figure>

milieu

<figure><img src="image_37.png" /></figure>

vitesse

<figure><img src="image_38.png" /></figure>

LTmm

Berge

<figure><img src="image_39.png" /></figure>

Macrophytes

<figure><img src="image_40.png" /></figure>

Annexe

<figure><img src="image_41.png" /></figure>

LTmm

Largeur_chenal

<figure><img src="image_42.png" /></figure>

Chenal

<figure><img src="image_43.png" /></figure>

Complexité hydraulique

<figure><img src="image_44.png" /></figure>

**Figure 36 - Représentation des LT médianes des chevaines d'âge 0+ classées selon certaines catégories propres aux stations.**
LTmm

situation

<figure><img src="image_36.png" /></figure>

milieu

<figure><img src="image_37.png" /></figure>

vitesse

<figure><img src="image_38.png" /></figure>

<figure><img src="image_39.png" /></figure>

Macrophytes

<figure><img src="image_40.png" /></figure>

Annexe

<figure><img src="image_41.png" /></figure>

<figure><img src="image_42.png" /></figure>

Largeur_chenal

<figure><img src="image_43.png" /></figure>

Chenal

<figure><img src="image_44.png" /></figure>

Complexité hydraulique

Figure 36 - Représentation des LT médianes des chevaines d'âge 0+ classées selon certaines catégories propres aux stations.
<figure><img src="image_45.png" /><figcaption>Figure 37 - Représentation des résidus de la régression sur la SDJ de la taille des chevaines d'âge 0+ classés selon certaines catégories propres aux stations.</figcaption></figure>
<figure><img src="image_45.png" /><figcaption>Figure 37 - Représentation des résidus de la régression sur la SDJ de la taille des chevaines d'âge 0+ classés selon certaines catégories propres aux stations.</figcaption></figure>
## 3.6 Autres facteurs de variabilité

Une part importante de la littérature relative à l'écologie des jeunes poissons en eaux courantes est consacrée aux effets directs des facteurs hydroclimatiques sur le recrutement et la croissance. Bien que les paramètres hydrologiques et thermiques soient intimement associés, l'hydrologie joue un rôle majeur dans le succès du recrutement et la température, agissant sur le métabolisme des organismes ectothermes, s'avère déterminante pour la croissance des poissons.

Cependant, nombre d'autres facteurs abiotiques et biotiques plus difficiles à appréhender agissent sur la reproduction et la croissance des poissons. Des publications récentes mettent l'accent sur des facteurs peu étudiés, tel que le parasitisme (LONGSHAW et al., 2010) ou une amélioration de la qualité des eaux (BEARDSLEY & BRITTON, 2012).

### 3.6.1 Evolution de la qualité des eaux

Bien que cela ne soit pas encore la panacée, les efforts contemporains réalisés pour améliorer la qualité des eaux, en particulier minimiser la pollution urbaine et industrielle, ont inversé des tendances observées avant épuration des eaux.

A partir de chroniques long terme sur la rivière Wensum en Grande-Bretagne, BEARDSLEY et BRITTON (2012) ont montré qu'un changement des procédés de traitement des eaux avait entraîné une réduction des apports en nutriments et un glissement progressif vers des conditions moins eutrophes. Sur une période de 25 ans, le gardon a montré une baisse significative de sa croissance. Les taux de croissance observés en 2005 et 2008 sont beaucoup plus bas que ceux enregistrés entre 1983 et 1994. Avant la mise en place de ces nouveaux procédés d'épuration, la thermique expliquait l'essentiel de la variabilité de la croissance. Cette relation n'est plus valide depuis la modification des conditions trophiques, et la croissance est désormais expliquée par les apports annuels d'ortho-phosphates.

Nos observations sur le long terme montrent des changements structurels inexpliqués, tels que l'augmentation de la proportion de juvéniles dans les pêches postérieures aux grands épisodes de crues, et à l'arrêt de rejets industriels de phosphates dans le tiers amont du Bas-Rhône.

### 3.6.2 Les facteurs biotiques

#### 3.6.2.1 La densité-dépendance

DELAIGUE et al. (2007) avaient proposé la densité-dépendance pour expliquer la taille des juvéniles. Nous rappelons que les années de faible recrutement avaient été « éliminées » du jeu de données du fait de la modicité des effectifs. Nous pouvons supposer que cette relation négative entre les effectifs et la taille pourrait être effective lorsque les effectifs deviennent importants.

Dans cette étude, aucune relation valide n'a été obtenue entre les effectifs et la taille des juvéniles d'âge 0+ (Figure 38). Des observations ponctuelles sur la chronique suggèrent cependant cet effet lors de recrutement massif et de faible mortalité hivernale (exemple du gardon en 1989).
## 3.6 Autres facteurs de variabilité

Une part importante de la littérature relative à l'écologie des jeunes poissons en eaux courantes est consacrée aux effets directs des facteurs hydroclimatiques sur le recrutement et la croissance. Bien que les paramètres hydrologiques et thermiques soient intimement associés, l'hydrologie joue un rôle majeur dans le succès du recrutement et la température, agissant sur le métabolisme des organismes ectothermes, s'avère déterminante pour la croissance des poissons.

Cependant, nombre d'autres facteurs abiotiques et biotiques plus difficiles à appréhender agissent sur la reproduction et la croissance des poissons. Des publications récentes mettent l'accent sur des facteurs peu étudiés, tel que le parasitisme (LONGSHAW et al., 2010) ou une amélioration de la qualité des eaux (BEARDSLEY & BRITTON, 2012).

### 3.6.1 Evolution de la qualité des eaux

Bien que cela ne soit pas encore la panacée, les efforts contemporains réalisés pour améliorer la qualité des eaux, en particulier minimiser la pollution urbaine et industrielle, ont inversé des tendances observées avant épuration des eaux.

A partir de chroniques long terme sur la rivière Wensum en Grande-Bretagne, BEARDSLEY et BRITTON (2012) ont montré qu'un changement des procédés de traitement des eaux avait entraîné une réduction des apports en nutriments et un glissement progressif vers des conditions moins eutrophes. Sur une période de 25 ans, le gardon a montré une baisse significative de sa croissance. Les taux de croissance observés en 2005 et 2008 sont beaucoup plus bas que ceux enregistrés entre 1983 et 1994. Avant la mise en place de ces nouveaux procédés d'épuration, la thermique expliquait l'essentiel de la variabilité de la croissance. Cette relation n'est plus valide depuis la modification des conditions trophiques, et la croissance est désormais expliquée par les apports annuels d'ortho-phosphates.

Nos observations sur le long terme montrent des changements structurels inexpliqués, tels que l'augmentation de la proportion de juvéniles dans les pêches postérieures aux grands épisodes de crues, et à l'arrêt de rejets industriels de phosphates dans le tiers amont du Bas-Rhône.

### 3.6.2 Les facteurs biotiques

#### 3.6.2.1 La densité-dépendance

DELAIGUE et al. (2007) avaient proposé la densité-dépendance pour expliquer la taille des juvéniles. Nous rappelons que les années de faible recrutement avaient été « éliminées » du jeu de données du fait de la modicité des effectifs. Nous pouvons supposer que cette relation négative entre les effectifs et la taille pourrait être effective lorsque les effectifs deviennent importants.

Dans cette étude, aucune relation valide n'a été obtenue entre les effectifs et la taille des juvéniles d'âge 0+ (Figure 38). Des observations ponctuelles sur la chronique suggèrent cependant cet effet lors de recrutement massif et de faible mortalité hivernale (exemple du gardon en 1989).
<figure><img src="image_46.png" /></figure>

**Figure 38** - Relation entre tailles médianes et nombre d'individus selon une échelle logarithmique, chez le gardon (Cruas aval) et le chevaine (Cruas amont) pour la chronique complète des pêches de ces deux stations.

### 3.6.2.2 L'état sanitaire des populations piscicoles

LONGSHAW et al. (2010) se sont intéressés aux effets du parasitisme sur la croissance au cours de la première année et sur le succès du recrutement du gardon et du chevaine. L'étude a été conduite sur 3 rivières du Yorkshire (Angleterre) entre 1993 et 2006. Elle a montré que des infections parasitaires des muscles avaient des effets négatifs majeurs sur la croissance des juvéniles 0+ et diminuaient le recrutement annuel des deux espèces. Dans ces cours d'eau de plaine, le parasitisme pourrait avoir une influence supérieure à celle des paramètres abiotiques.

Les études spécifiques sur l'état sanitaire des poissons du Rhône n'ont réellement démarré qu'après l'épisode caniculaire 2003. Mais il est probable que les faibles effectifs de juvéniles dans les pêches du début de la chronique aient une part d'explication dans l'épizootie de bucéphalose signalée sur le fleuve au cours des premières études d'impact.

### 3.6.2.3 Les ressources alimentaires

Parmi les métriques thermiques utilisées par WOLTER (2007) dans le cadre des études piscicoles de la Basse-Oder (Allemagne), cet auteur a fait une association entre la durée de la période de développement des macrophytes (une durée fixée par l'application d'un même seuil thermique au printemps et à l'automne) et un cumul de DJ sur cette période<fnref n="5" />. La métrique VegSum est corrélée positivement avec la taille des jeunes gardons. Elle associe les effets saisonniers de plusieurs composantes climatiques (la température et l'ensoleillement) sur la présence, le développement des herbiers, et la croissance des jeunes poissons utilisant les macrophytes comme abris et zones préférentielles d'alimentation. D'une certaine manière, la conjonction de la durée de vie de l'abri végétal et de la température des eaux favorisent la croissance de ses occupants.

Les données montrent que la taille est maximale dans la Saône. L'abondance des herbiers et l'eutrophie de cette grande rivière de plaine pourrait favoriser la croissance. Cet effet est vraisemblablement important dans toutes les grandes retenues, en particulier sur le

<footnote n="5">⁵ VegPeriod : Number of days from first reaching and last having 10°C water temperature, VegSumme : Cumulative degree-days of VegPeriod, VegSum12 : Cumulative degree-days above 12°C of VegPeriod.</footnote>
<figure><img src="image_46.png" /></figure>

**Figure 38** - Relation entre tailles médianes et nombre d'individus selon une échelle logarithmique, chez le gardon (Cruas aval) et le chevaine (Cruas amont) pour la chronique complète des pêches de ces deux stations.

### 3.6.2.2 L'état sanitaire des populations piscicoles

LONGSHAW et al. (2010) se sont intéressés aux effets du parasitisme sur la croissance au cours de la première année et sur le succès du recrutement du gardon et du chevaine. L'étude a été conduite sur 3 rivières du Yorkshire (Angleterre) entre 1993 et 2006. Elle a montré que des infections parasitaires des muscles avaient des effets négatifs majeurs sur la croissance des juvéniles 0+ et diminuaient le recrutement annuel des deux espèces. Dans ces cours d'eau de plaine, le parasitisme pourrait avoir une influence supérieure à celle des paramètres abiotiques.

Les études spécifiques sur l'état sanitaire des poissons du Rhône n'ont réellement démarré qu'après l'épisode caniculaire 2003. Mais il est probable que les faibles effectifs de juvéniles dans les pêches du début de la chronique aient une part d'explication dans l'épizootie de bucéphalose signalée sur le fleuve au cours des premières études d'impact.

### 3.6.2.3 Les ressources alimentaires

Parmi les métriques thermiques utilisées par WOLTER (2007) dans le cadre des études piscicoles de la Basse-Oder (Allemagne), cet auteur a fait une association entre la durée de la période de développement des macrophytes (une durée fixée par l'application d'un même seuil thermique au printemps et à l'automne) et un cumul de DJ sur cette période<fnref n="5" />. La métrique VegSum est corrélée positivement avec la taille des jeunes gardons. Elle associe les effets saisonniers de plusieurs composantes climatiques (la température et l'ensoleillement) sur la présence, le développement des herbiers, et la croissance des jeunes poissons utilisant les macrophytes comme abris et zones préférentielles d'alimentation. D'une certaine manière, la conjonction de la durée de vie de l'abri végétal et de la température des eaux favorisent la croissance de ses occupants.

Les données montrent que la taille est maximale dans la Saône. L'abondance des herbiers et l'eutrophie de cette grande rivière de plaine pourrait favoriser la croissance. Cet effet est vraisemblablement important dans toutes les grandes retenues, en particulier sur le

<footnote n="5">⁵ VegPeriod : Number of days from first reaching and last having 10°C water temperature, VegSumme : Cumulative degree-days of VegPeriod, VegSum12 : Cumulative degree-days above 12°C of VegPeriod.</footnote>
Bas-Rhône. Il est probable que les périodes de déficit hydrologique et les étiages prolongés favorisent le développement massif des herbiers et indirectement la croissance des jeunes poissons.

# 4 Conclusions

Les chroniques hydrobiologiques du Rhône présentent les avantages de leur durée et de leur homogénéité en termes de méthodes et sites d'échantillonnage pour pouvoir évaluer l'évolution spatio-temporelle des assemblages faunistiques. Réalisées dans le cadre des suivis hydrobiologiques permanents des sites industriels, ou celui de la restauration du fleuve, elles peuvent être utiles pour aborder d'autres questions (ROGER & FAESSEL, 1989; ROGER et al., 1991; DESSAIX et al., 1995; DOLEDEC et al., 1996; CATTANEO et al., 2001; FRUGET et al., 2001; GRENOUILLET et al., 2001; DAUFRESNE et al., 2003; DAUFRESNE et al., 2007; DELAIGUE et al., 2007; DESSAIX & FRUGET, 2008; OLIVIER et al., 2009; PIFFADY et al., 2010).

Cependant, les objectifs des suivis hydrobiologiques permanents ont orienté le choix des stations (par ex. amont, rejet, aval d'un site), des méthodes d'échantillonnage (pêches électriques) et leur saisonnalité, la nature (vision globale de l'assemblage) et le niveau de précision dans l'acquisition de ces données (par ex. la constitution de lots lorsque la fraction juvénile est trop importante). Leur usage pour satisfaire de nouveaux objectifs, ou répondre à des questions nouvelles soulevées par les gestionnaires, n'est pas sans problèmes.

Les chroniques piscicoles ont montré la part croissante des juvéniles dans les captures, et la nécessité d'approfondir les connaissances sur le fonctionnement écologique du fleuve à partir de cette fraction juvénile très sensible aux modifications anthropiques. Des propositions dans ce sens ont déjà été faites par des collègues européens (JURAJDA et al., 2010) lorsque l'échantillonnage des adultes pose de réels problèmes.

Une première étude a été effectuée par DELAIGUE et al. (2007) afin d'étudier la variabilité de la taille et du coefficient de condition des juvéniles de l'année en utilisant des pêches collectées sur le Rhône par Irstea Aix-en-Provence. Ces données concernaient essentiellement le Bas-Rhône médian (de la retenue de Montélimar à celle de Caderousse). Les pêches comportant peu d'individus juvéniles (n<30) avaient été écartées des analyses. Or ces pêches correspondent le plus souvent à des années de faible recrutement, caractérisées par des hautes eaux printanières et un printemps froid. Par conséquent, les analyses ont privilégié les années à fort recrutement, plutôt chaudes. L'extension de cette étude préliminaire à l'échelle du bassin versant (183 stations dans le cadre de cette compilation) a permis d'avoir une amplitude thermique plus grande, et des températures froides sans pour autant avoir de faibles effectifs de juvéniles. Grâce à nos partenaires de l'étude Rhône Phase 4, nous avons pu bénéficier de chroniques physiques (relevés des températures eDF-DTG) et piscicoles à l'échelle du Rhône français susceptibles de fournir des situations hydroclimatiques contrastées et biologiquement renseignées.

Après avoir validé les données et les avoir intégrées dans une base de données, il a fallu procéder à l'identification des âges à partir des distributions de taille en l'absence d'informations sur l'âge des poissons capturés. Une part très importante de ce travail a été consacrée à cette identification des classes d'âge, sur la base d'une analyse par expertise des histogrammes de tailles pour avoir entre autres les tailles des juvéniles de l'année (âge 0+).

La complexité de cette tâche est liée (1) à l'extrême variabilité saisonnière des effectifs, (2) à une variabilité du recrutement d'une année sur l'autre, (3) au gradient thermique créé par la répartition longitudinale des stations sur environ 450 km, depuis la Chautagne jusqu'à Vallabrègues, et (4) un décalage progressif des modes des histogrammes de taille des poissons 0+ sur trois décennies, surtout perceptible à partir des années 1990, du fait des évolutions thermiques actuelles.
Bas-Rhône. Il est probable que les périodes de déficit hydrologique et les étiages prolongés favorisent le développement massif des herbiers et indirectement la croissance des jeunes poissons.

# 4 Conclusions

Les chroniques hydrobiologiques du Rhône présentent les avantages de leur durée et de leur homogénéité en termes de méthodes et sites d'échantillonnage pour pouvoir évaluer l'évolution spatio-temporelle des assemblages faunistiques. Réalisées dans le cadre des suivis hydrobiologiques permanents des sites industriels, ou celui de la restauration du fleuve, elles peuvent être utiles pour aborder d'autres questions (ROGER & FAESSEL, 1989; ROGER et al., 1991; DESSAIX et al., 1995; DOLEDEC et al., 1996; CATTANEO et al., 2001; FRUGET et al., 2001; GRENOUILLET et al., 2001; DAUFRESNE et al., 2003; DAUFRESNE et al., 2007; DELAIGUE et al., 2007; DESSAIX & FRUGET, 2008; OLIVIER et al., 2009; PIFFADY et al., 2010).

Cependant, les objectifs des suivis hydrobiologiques permanents ont orienté le choix des stations (par ex. amont, rejet, aval d'un site), des méthodes d'échantillonnage (pêches électriques) et leur saisonnalité, la nature (vision globale de l'assemblage) et le niveau de précision dans l'acquisition de ces données (par ex. la constitution de lots lorsque la fraction juvénile est trop importante). Leur usage pour satisfaire de nouveaux objectifs, ou répondre à des questions nouvelles soulevées par les gestionnaires, n'est pas sans problèmes.

Les chroniques piscicoles ont montré la part croissante des juvéniles dans les captures, et la nécessité d'approfondir les connaissances sur le fonctionnement écologique du fleuve à partir de cette fraction juvénile très sensible aux modifications anthropiques. Des propositions dans ce sens ont déjà été faites par des collègues européens (JURAJDA et al., 2010) lorsque l'échantillonnage des adultes pose de réels problèmes.

Une première étude a été effectuée par DELAIGUE et al. (2007) afin d'étudier la variabilité de la taille et du coefficient de condition des juvéniles de l'année en utilisant des pêches collectées sur le Rhône par Irstea Aix-en-Provence. Ces données concernaient essentiellement le Bas-Rhône médian (de la retenue de Montélimar à celle de Caderousse). Les pêches comportant peu d'individus juvéniles (n<30) avaient été écartées des analyses. Or ces pêches correspondent le plus souvent à des années de faible recrutement, caractérisées par des hautes eaux printanières et un printemps froid. Par conséquent, les analyses ont privilégié les années à fort recrutement, plutôt chaudes. L'extension de cette étude préliminaire à l'échelle du bassin versant (183 stations dans le cadre de cette compilation) a permis d'avoir une amplitude thermique plus grande, et des températures froides sans pour autant avoir de faibles effectifs de juvéniles. Grâce à nos partenaires de l'étude Rhône Phase 4, nous avons pu bénéficier de chroniques physiques (relevés des températures eDF-DTG) et piscicoles à l'échelle du Rhône français susceptibles de fournir des situations hydroclimatiques contrastées et biologiquement renseignées.

Après avoir validé les données et les avoir intégrées dans une base de données, il a fallu procéder à l'identification des âges à partir des distributions de taille en l'absence d'informations sur l'âge des poissons capturés. Une part très importante de ce travail a été consacrée à cette identification des classes d'âge, sur la base d'une analyse par expertise des histogrammes de tailles pour avoir entre autres les tailles des juvéniles de l'année (âge 0+).

La complexité de cette tâche est liée (1) à l'extrême variabilité saisonnière des effectifs, (2) à une variabilité du recrutement d'une année sur l'autre, (3) au gradient thermique créé par la répartition longitudinale des stations sur environ 450 km, depuis la Chautagne jusqu'à Vallabrègues, et (4) un décalage progressif des modes des histogrammes de taille des poissons 0+ sur trois décennies, surtout perceptible à partir des années 1990, du fait des évolutions thermiques actuelles.
Initialement, nous avons travaillé sur les cyprinidés dominants. Parmi les grandes espèces (chevaine, gardon, barbeau, hotu, vandoise, brème bordelière), seuls le chevaine et le gardon présentent une répartition spatiale et temporelle sur l'ensemble des tronçons étudiés. Elles sont très documentées par la littérature scientifique européenne. Ces deux espèces ont donc fait l'objet des nombreux allers-retours dans la boucle « vérification-correction-validation », aidé dans cette opération par des d'outils graphiques conçus pour simplifier l'expertise, et l'usage de l'analyse multivariée (PERSAT & CHESSEL, 1989). Nous avons ainsi pu obtenir les tableaux finaux a priori valides pour la suite des analyses.

Les autres grandes espèces ont été triées mais la validation n'a pas été réalisée. Ces espèces présentent un intérêt pour des approches géographiques plus réduites car leur répartition est spatialement limitée. Certaines d'entre elles sont quasiment absentes des pêches pendant de longues périodes. Les espèces de petite taille sont des espèces à pontes multiples (ablette, spirlin, goujon, bouvière, pseudorasbora). L'identification des classes d'âge est plus complexe encore en raison de la présence de plusieurs modes ; et la petite taille des adultes amplifie le chevauchement des cohortes annuelles. Elles nécessitent des travaux spécifiques sur la croissance, avec identification des âges par des techniques appropriées. C'est surtout le cas de la bouvière et du pseudorasbora, deux espèces dont les effectifs n'ont fait qu'augmenter au cours des dernières années.

Un travail spécifique portant sur des méthodes de tri automatique a été réalisé afin de construire de nouveaux outils permettant de réduire le temps consacré à la discrimination des classes d'âge. Les essais ont été réalisés à partir de la base de tri expertisée (VESLOT, 2012).

L'étude de la croissance saisonnière a été réalisée sur la base de plusieurs pêches échelonnées entre mai et septembre sur quelques sites, avec toutes les réserves liées aux problèmes d'échantillonnage des stades les plus jeunes. Elle a été complétée par des suivis spécifiques réalisés en 1983 sur des annexes fluviales du Haut-Rhône (CARREL, 1986). Pour répondre à notre premier objectif, à savoir l'influence de la température de la période de croissance sur la taille des poissons en fin de leur première année, il fallait s'assurer de la relative stabilité du palier atteint en fin d'été pour ne pas biaiser cette relation par la prise en compte d'individus n'ayant pas achevé leur première saison de croissance. L'étude de la taille est faite sur des poissons capturés après la date moyenne du maximum thermique (8 août).

Une majorité de travaux portant sur l'étude de la croissance des organismes ectothermes utilise la somme des degrés-jours pour caractériser le « temps physiologique » (combinaison du temps calendaire et de la température) nécessaire pour atteindre un stade de développement (éclosion, résorption de vésicule vitelline, ...). Bien que la sommation des degrés-jours puisse être faite avec un seuil fixé à 0°C, des seuils thermiques plus élevés sont utilisés en fonction de l'organisme étudié, considérant qu'il n'y a pas de croissance effective en dessous de ce seuil. Le seuil de 12°C proposé par Richard Mann est le plus communément utilisé pour les cyprinidés. Toutefois, la variabilité des températures de ponte indiquées dans la littérature nous a incités à prendre un seuil thermique beaucoup plus bas (8.8°C) correspondant à la température moyenne du 20 mars (équinoxe de printemps) pour l'ensemble du Rhône. Sachant que les cyprinidés terminent leur gamétogénèse au cours du printemps, nous pouvons considérer que le choix de ce seuil nous permet de prendre en compte une énergie utile pour les adultes reproducteurs, favorable pour leur migration vers les frayères (WOLTER, 2007), et pouvant générer une précocité de la ponte. Nous avons calculé une somme de degrés-jours au dessus de 8.8°C (SDJ) pour la période dite de croissance comprise entre le 20 mars et le 8 août de chaque année et chaque station thermique.

L'hydrologie est contraignante pour la reproduction et la croissance des poissons (CATTANEO et al., 2001; PETERSON & JENNINGS, 2007). Nous avons opté pour un cumul de débits-jours pour la période de croissance. L'indice de débit a été calculé en sommant les
Initialement, nous avons travaillé sur les cyprinidés dominants. Parmi les grandes espèces (chevaine, gardon, barbeau, hotu, vandoise, brème bordelière), seuls le chevaine et le gardon présentent une répartition spatiale et temporelle sur l'ensemble des tronçons étudiés. Elles sont très documentées par la littérature scientifique européenne. Ces deux espèces ont donc fait l'objet des nombreux allers-retours dans la boucle « vérification-correction-validation », aidé dans cette opération par des d'outils graphiques conçus pour simplifier l'expertise, et l'usage de l'analyse multivariée (PERSAT & CHESSEL, 1989). Nous avons ainsi pu obtenir les tableaux finaux a priori valides pour la suite des analyses.

Les autres grandes espèces ont été triées mais la validation n'a pas été réalisée. Ces espèces présentent un intérêt pour des approches géographiques plus réduites car leur répartition est spatialement limitée. Certaines d'entre elles sont quasiment absentes des pêches pendant de longues périodes. Les espèces de petite taille sont des espèces à pontes multiples (ablette, spirlin, goujon, bouvière, pseudorasbora). L'identification des classes d'âge est plus complexe encore en raison de la présence de plusieurs modes ; et la petite taille des adultes amplifie le chevauchement des cohortes annuelles. Elles nécessitent des travaux spécifiques sur la croissance, avec identification des âges par des techniques appropriées. C'est surtout le cas de la bouvière et du pseudorasbora, deux espèces dont les effectifs n'ont fait qu'augmenter au cours des dernières années.

Un travail spécifique portant sur des méthodes de tri automatique a été réalisé afin de construire de nouveaux outils permettant de réduire le temps consacré à la discrimination des classes d'âge. Les essais ont été réalisés à partir de la base de tri expertisée (VESLOT, 2012).

L'étude de la croissance saisonnière a été réalisée sur la base de plusieurs pêches échelonnées entre mai et septembre sur quelques sites, avec toutes les réserves liées aux problèmes d'échantillonnage des stades les plus jeunes. Elle a été complétée par des suivis spécifiques réalisés en 1983 sur des annexes fluviales du Haut-Rhône (CARREL, 1986). Pour répondre à notre premier objectif, à savoir l'influence de la température de la période de croissance sur la taille des poissons en fin de leur première année, il fallait s'assurer de la relative stabilité du palier atteint en fin d'été pour ne pas biaiser cette relation par la prise en compte d'individus n'ayant pas achevé leur première saison de croissance. L'étude de la taille est faite sur des poissons capturés après la date moyenne du maximum thermique (8 août).

Une majorité de travaux portant sur l'étude de la croissance des organismes ectothermes utilise la somme des degrés-jours pour caractériser le « temps physiologique » (combinaison du temps calendaire et de la température) nécessaire pour atteindre un stade de développement (éclosion, résorption de vésicule vitelline, ...). Bien que la sommation des degrés-jours puisse être faite avec un seuil fixé à 0°C, des seuils thermiques plus élevés sont utilisés en fonction de l'organisme étudié, considérant qu'il n'y a pas de croissance effective en dessous de ce seuil. Le seuil de 12°C proposé par Richard Mann est le plus communément utilisé pour les cyprinidés. Toutefois, la variabilité des températures de ponte indiquées dans la littérature nous a incités à prendre un seuil thermique beaucoup plus bas (8.8°C) correspondant à la température moyenne du 20 mars (équinoxe de printemps) pour l'ensemble du Rhône. Sachant que les cyprinidés terminent leur gamétogénèse au cours du printemps, nous pouvons considérer que le choix de ce seuil nous permet de prendre en compte une énergie utile pour les adultes reproducteurs, favorable pour leur migration vers les frayères (WOLTER, 2007), et pouvant générer une précocité de la ponte. Nous avons calculé une somme de degrés-jours au dessus de 8.8°C (SDJ) pour la période dite de croissance comprise entre le 20 mars et le 8 août de chaque année et chaque station thermique.

L'hydrologie est contraignante pour la reproduction et la croissance des poissons (CATTANEO et al., 2001; PETERSON & JENNINGS, 2007). Nous avons opté pour un cumul de débits-jours pour la période de croissance. L'indice de débit a été calculé en sommant les
débits journaliers exprimés par les centiles normés issus des chroniques complètes de débits (données HYDRO). Du fait d'une relation inverse entre le débit et la température, la taille diminue lorsque les débits printaniers sont soutenus car la température des eaux est alors plus froide.

Les tailles des gardons et des chevaines sont globalement croissantes de l'amont vers l'aval. Plus petites sur le Haut-Rhône amont et sur l'Ain, elles augmentent sur le Haut-Rhône aval, puis nettement sur le Bas-Rhône. Les tailles maximales ont été enregistrées sur la Saône, sur le Bas-Rhône aval, et au niveau des rejets thermiques des centrales. Pour ces deux espèces, nous avons obtenu une relation linéaire positive associant la taille des juvéniles mesurée après le maximum thermique et la température printanière (traduite par le cumul de degrés-jours).

Avec la hausse des températures attribuée au changement climatique, la période favorable pour la croissance s'allonge du fait d'une précocité printanière de plus en plus marquée. Une prolongation de la période chaude en automne a été observée plusieurs fois dans les années récentes. Les résultats des pêches faites en 2011, l'année thermique la plus chaude des chroniques météorologiques, ont montré une précocité exceptionnelle, et la poursuite de la croissance en cours d'automne, mise en évidence sur des sites ayant deux pêches successives en fin d'année. A l'heure actuelle, du fait des évolutions climatiques, la taille des juvéniles de l'année de ces deux espèces tend à augmenter. Cette tendance a débuté au cours des années 1990.

Contrairement à de nombreux travaux réalisés sur de grands cours d'eau de plaine, qui ont mis en évidence des effets de compétition sur la taille, nous n'avons pas obtenu de relations significatives négatives (DELAIGUE et al., 2007; WOLTER, 2007) entre la taille et les effectifs de poissons capturés, une diminution de la taille attribuée à la compétition intra-et/ou interspécifique pour l'exploitation de ressources trophiques (DELAIGUE et al., 2007; WOLTER, 2007). Il est probable que ce constat soit en partie attribuable au niveau trophique des hydrosystèmes étudiés (rivières et grands cours d'eau de plaine).

D'autres facteurs ayant des effets majeurs sur la croissance ont été étudiés dans des publications récentes, tels que le parasitisme des jeunes poissons (diminution de la taille) (LONGSHAW et al., 2010), ou l'amélioration de la qualité des eaux (diminution de la taille) (BEARDSLEY & BRITTON, 2012). Ils montrent que des facteurs sans liens directs avec l'hydro-climatologie ne peuvent être sous-estimés.

Nous avons une base de données dans laquelle l'identification des poissons d'âge 0+ de deux espèces peut être considérée comme valide. Le tri préliminaire a été réalisé sur 7 autres espèces (ablette, barbeau, blageon, brème bordelière, goujon, hotu, rotengle, spirlin et vandoise). Le travail accompli pour l'analyse de la taille des gardons et des chevaines mériterait d'être assuré pour ces 7 autres espèces. Un travail préparatoire a été réalisé pour deux petites espèces (bouvière et pseudorasbora), mais la complexité du tri nous a conduis à leur abandon. Ces espèces nécessitent des informations complémentaires (lecture d'otolithes) pour garantir la véracité d'un tri par expertise sur les tailles.
débits journaliers exprimés par les centiles normés issus des chroniques complètes de débits (données HYDRO). Du fait d'une relation inverse entre le débit et la température, la taille diminue lorsque les débits printaniers sont soutenus car la température des eaux est alors plus froide.

Les tailles des gardons et des chevaines sont globalement croissantes de l'amont vers l'aval. Plus petites sur le Haut-Rhône amont et sur l'Ain, elles augmentent sur le Haut-Rhône aval, puis nettement sur le Bas-Rhône. Les tailles maximales ont été enregistrées sur la Saône, sur le Bas-Rhône aval, et au niveau des rejets thermiques des centrales. Pour ces deux espèces, nous avons obtenu une relation linéaire positive associant la taille des juvéniles mesurée après le maximum thermique et la température printanière (traduite par le cumul de degrés-jours).

Avec la hausse des températures attribuée au changement climatique, la période favorable pour la croissance s'allonge du fait d'une précocité printanière de plus en plus marquée. Une prolongation de la période chaude en automne a été observée plusieurs fois dans les années récentes. Les résultats des pêches faites en 2011, l'année thermique la plus chaude des chroniques météorologiques, ont montré une précocité exceptionnelle, et la poursuite de la croissance en cours d'automne, mise en évidence sur des sites ayant deux pêches successives en fin d'année. A l'heure actuelle, du fait des évolutions climatiques, la taille des juvéniles de l'année de ces deux espèces tend à augmenter. Cette tendance a débuté au cours des années 1990.

Contrairement à de nombreux travaux réalisés sur de grands cours d'eau de plaine, qui ont mis en évidence des effets de compétition sur la taille, nous n'avons pas obtenu de relations significatives négatives (DELAIGUE et al., 2007; WOLTER, 2007) entre la taille et les effectifs de poissons capturés, une diminution de la taille attribuée à la compétition intra-et/ou interspécifique pour l'exploitation de ressources trophiques (DELAIGUE et al., 2007; WOLTER, 2007). Il est probable que ce constat soit en partie attribuable au niveau trophique des hydrosystèmes étudiés (rivières et grands cours d'eau de plaine).

D'autres facteurs ayant des effets majeurs sur la croissance ont été étudiés dans des publications récentes, tels que le parasitisme des jeunes poissons (diminution de la taille) (LONGSHAW et al., 2010), ou l'amélioration de la qualité des eaux (diminution de la taille) (BEARDSLEY & BRITTON, 2012). Ils montrent que des facteurs sans liens directs avec l'hydro-climatologie ne peuvent être sous-estimés.

Nous avons une base de données dans laquelle l'identification des poissons d'âge 0+ de deux espèces peut être considérée comme valide. Le tri préliminaire a été réalisé sur 7 autres espèces (ablette, barbeau, blageon, brème bordelière, goujon, hotu, rotengle, spirlin et vandoise). Le travail accompli pour l'analyse de la taille des gardons et des chevaines mériterait d'être assuré pour ces 7 autres espèces. Un travail préparatoire a été réalisé pour deux petites espèces (bouvière et pseudorasbora), mais la complexité du tri nous a conduis à leur abandon. Ces espèces nécessitent des informations complémentaires (lecture d'otolithes) pour garantir la véracité d'un tri par expertise sur les tailles.
# 5 Bibliographie

ALABASTER, J.S. & LLOYD, R. (1982): Water quality criteria for freshwater fish. - Butterworths, London, 361 p.

ALLARDI, J., DUGUET, C. & LEYNAUD, G. (1975): Etude du peuplement ichtyologique d'un milieu fluvial réchauffé artificiellement. - In: WELCOME, R.L. (Ed.), Symposium on the methodology for the survey, monitoring and appraisal of fishery resources in lakes and large rivers. EIFAC Technical Paper, FAO, Rome, pp. 209-227.

BALON, E.K. (1975): Reproductive guilds of fishes: a proposal and definition. - J. Fish. Res. Board Can. **32**: 821-864.

BEARDSLEY, H. & BRITTON, J.R. (2012): Contribution of temperature and nutrient loading to growth rate variation of three cyprinid fishes in a lowland river. - Aquat. Ecol. **46**: 143-152.

BOISNEAU, C., MOATAR, F., BODIN, M. & POISNEAU, P. (2008): Does global warming impact on migration patterns and recruitment of Allis shad (*Alosa alosa* L.) young of the year in the Loire River, France? - Hydrobiologia **602**: 179-186.

BOLLAND, J.D., BRITTON, J.R. & COWX, I.G. (2007): Lifetime consequences of variable 0 year group length in riverine populations of chub *Leuciscus cephalus* (L.). - J. Fish Biol. **71**: 1810-1819.

BRITTON, J.R., DAVIES, G.D. & BRAZIER, M. (2009): Eradication of the invasive *Pseudorasbora parva* results in increased growth and production of native fishes. - Ecol. Freshwat. Fish **18**: 8-14.

CARREL, G. (1981): Etude physico-chimique et piscicole d'une lône du Haut-Rhône français. Analyse de plusieurs méthodes d'observation et de comptage des poissons en eau courantes. - DEA Ecologie fondamentale et appliquée des eaux continentales, Université Lyon 1, Villeurbanne, 73 p.

CARREL, G. (1986): Caractérisation physico-chimique du Haut-Rhône français et de ses annexes ; incidences sur la croissance des populations d'alevins. - Thèse de Doctorat, Université Claude Bernard - Lyon I, Villeurbanne, France, 186 p.

CARREL, G., OLIVIER, J.-M., POIREL, A. & BREIL, P. (2003): Analyse des réponses biologiques à la variabilité climatique par l'examen de chroniques temporelles. Variabilite inter-annuelle de la reproduction des poissons cyprinidae dans le Rhône a l'aval de Lyon : rôle de la température. - In: PONT, D. (Ed.), Conséquences potentielles du changement climatique sur les biocénoses aquatiques et riveraines françaises. Rapport final MEDD - GICC Aquabio, Paris.

CATTANEO, F., BREIL, P. & CARREL, G. (1999): Variabilité inter-annuelle des assemblages piscicoles en relation avec le régime hydrologique : le Bas-Rhône à Montélimar, 4ème Congrès International de Limnologie-Océanographie, Bordeaux, 7-10 Septembre 1999.

CATTANEO, F., CARREL, G., LAMOUROUX, N. & BREIL, P. (2001): Relationship between hydrology and cyprinid reproductive success in the Lower Rhône at Montélimar, France. - Arch. Hydrobiol. **151**: 427-450.

COPP, G.H. (1987): Le rôle et le fonctionnement des milieux aquatiques du Haut-Rhône français comme sites de reproduction et de nurserie pour les poissons du fleuve. - Doctorat de l'Université, Lyon I, 97 p.

COPP, G.H. (1989a): Electrofishing for fish larvae and juveniles: equipment modifications for increased efficiency with short fishes. - Aquacult. Fish. Manage. **20**: 453-462.

COPP, G.H. (1989b): The habitat diversity and fish reproductive function of floodplain ecosystems. - Environ. Biol. Fish. **26**: 1-27.

CRYER, M., PEIRSON, G. & TOWNSEND, C.R. (1986): Reciprocal interactions between roach, *Rutilus rutilus*, and zooplankton in a small lake: prey dynamics and fish growth and recruitment. - Limnol. Oceanogr. **31**: 1022-1038.
# 5 Bibliographie

ALABASTER, J.S. & LLOYD, R. (1982): Water quality criteria for freshwater fish. - Butterworths, London, 361 p.

ALLARDI, J., DUGUET, C. & LEYNAUD, G. (1975): Etude du peuplement ichtyologique d'un milieu fluvial réchauffé artificiellement. - In: WELCOME, R.L. (Ed.), Symposium on the methodology for the survey, monitoring and appraisal of fishery resources in lakes and large rivers. EIFAC Technical Paper, FAO, Rome, pp. 209-227.

BALON, E.K. (1975): Reproductive guilds of fishes: a proposal and definition. - J. Fish. Res. Board Can. **32**: 821-864.

BEARDSLEY, H. & BRITTON, J.R. (2012): Contribution of temperature and nutrient loading to growth rate variation of three cyprinid fishes in a lowland river. - Aquat. Ecol. **46**: 143-152.

BOISNEAU, C., MOATAR, F., BODIN, M. & POISNEAU, P. (2008): Does global warming impact on migration patterns and recruitment of Allis shad (*Alosa alosa* L.) young of the year in the Loire River, France? - Hydrobiologia **602**: 179-186.

BOLLAND, J.D., BRITTON, J.R. & COWX, I.G. (2007): Lifetime consequences of variable 0 year group length in riverine populations of chub *Leuciscus cephalus* (L.). - J. Fish Biol. **71**: 1810-1819.

BRITTON, J.R., DAVIES, G.D. & BRAZIER, M. (2009): Eradication of the invasive *Pseudorasbora parva* results in increased growth and production of native fishes. - Ecol. Freshwat. Fish **18**: 8-14.

CARREL, G. (1981): Etude physico-chimique et piscicole d'une lône du Haut-Rhône français. Analyse de plusieurs méthodes d'observation et de comptage des poissons en eau courantes. - DEA Ecologie fondamentale et appliquée des eaux continentales, Université Lyon 1, Villeurbanne, 73 p.

CARREL, G. (1986): Caractérisation physico-chimique du Haut-Rhône français et de ses annexes ; incidences sur la croissance des populations d'alevins. - Thèse de Doctorat, Université Claude Bernard - Lyon I, Villeurbanne, France, 186 p.

CARREL, G., OLIVIER, J.-M., POIREL, A. & BREIL, P. (2003): Analyse des réponses biologiques à la variabilité climatique par l'examen de chroniques temporelles. Variabilite inter-annuelle de la reproduction des poissons cyprinidae dans le Rhône a l'aval de Lyon : rôle de la température. - In: PONT, D. (Ed.), Conséquences potentielles du changement climatique sur les biocénoses aquatiques et riveraines françaises. Rapport final MEDD - GICC Aquabio, Paris.

CATTANEO, F., BREIL, P. & CARREL, G. (1999): Variabilité inter-annuelle des assemblages piscicoles en relation avec le régime hydrologique : le Bas-Rhône à Montélimar, 4ème Congrès International de Limnologie-Océanographie, Bordeaux, 7-10 Septembre 1999.

CATTANEO, F., CARREL, G., LAMOUROUX, N. & BREIL, P. (2001): Relationship between hydrology and cyprinid reproductive success in the Lower Rhône at Montélimar, France. - Arch. Hydrobiol. **151**: 427-450.

COPP, G.H. (1987): Le rôle et le fonctionnement des milieux aquatiques du Haut-Rhône français comme sites de reproduction et de nurserie pour les poissons du fleuve. - Doctorat de l'Université, Lyon I, 97 p.

COPP, G.H. (1989a): Electrofishing for fish larvae and juveniles: equipment modifications for increased efficiency with short fishes. - Aquacult. Fish. Manage. **20**: 453-462.

COPP, G.H. (1989b): The habitat diversity and fish reproductive function of floodplain ecosystems. - Environ. Biol. Fish. **26**: 1-27.

CRYER, M., PEIRSON, G. & TOWNSEND, C.R. (1986): Reciprocal interactions between roach, *Rutilus rutilus*, and zooplankton in a small lake: prey dynamics and fish growth and recruitment. - Limnol. Oceanogr. **31**: 1022-1038.
DAUFRESNE, M., BADY, P. & FRUGET, J.F. (2007): Impacts of global changes and extreme hydroclimatic events on macroinvertebrate community structures in the French Rhone River. *Oecologia* **151**: 544-559.

DAUFRESNE, M. & BOET, P. (2005): Evolution temporelle de la structure des peuplements piscicoles des grands fleuves français. EDF - Cemagref, Antony, 19 p.

DAUFRESNE, M., LENGFELLNER, K. & SOMMER, U. (2009): Global warming benefits the small in aquatic ecosystems. - Proc. Natl. Acad. Sci. U. S. A. **106**: 12788-12793.

DAUFRESNE, M., ROGER, M.C., CAPRA, H. & LAMOUROUX, N. (2003): Long-term changes within the invertebrate and fish communities of the Upper Rhône River: effects of climatic factors. - Glob. Change Biol. **10**: 124-140.

DE KINKELIN, P., TUFFERY, G., LEYNAUD, G., ARRIGNON, J., POIRIER, M., JAMAIN, M.M., JOANNIS, G., JOLLET, G. & LECLERCQ, J. (1969): La Bucéphalose larvaire à *Bucephalus polymorphus* (Baer 1827) : pathogénie, épizootiologie, possibilités d'intervention. - Bull. fr. Piscic. **234**: 5-20.

DELAIGUE, O., CARREL, G. & PONT, D. (2007): Influence de la température sur la croissance des juvéniles de poissons Cyprinidae. Cemagref, U.R. Hydrobiologie, Aix-en-Provence, France, 44 p.

DESAINT, B. & POIREL, A. (2004): Etude thermique du Rhône - Phase 2 - Rapport final. - D4161/RAP/2004-00314-A. EDF - DTG - Dpt Surveillance - Branche Environnement aquatique, Grenoble, 88 p.

DESSAIX, J., FRUGET, J.-F., OLIVIER, J.-M. & BEFFY, J.-L. (1995): Changes of the macroinvertebrate communities in the dammed and by-passed sections of the French Upper Rhône after regulation. - Regul. Rivers **10**: 265-279.

DESSAIX, J. & FRUGET, J.F. (2008): Évolution des peuplements de crustacés du Rhône Moyen au cours des 20 dernières années, relation avec la variabilité hydroclimatique. - Hydroécol. Appl. **16**: 1-27.

DOLEDEC, S., DESSAIX, J. & TACHET, H. (1996): Changes within the Upper Rhône River macrobenthic communities after the completion of three hydroelectric schemes: anthropogenic effects or natural change? - Arch. Hydrobiol. **136**: 19-40.

FECHHELM, R.G., BRYAN, J.D., GRIFFITHS, W.B. & MARTIN, L.R. (1997): Summer growth patterns of northern Dolly Varden (*Salvelinus malma*) smolts from the Prudhoe Bay region of Alaska. - Can. J. Fish. Aquat. Sci. **54**: 1103-1110.

FOX, C.J., PLANQUE, B.P. & DARBY, C.D. (2000): Synchrony in the recruitment time-series of plaice (*Pleuronectes platessa* L) around the United Kingdom and the influence of sea temperature. - J. Sea Res. **44**: 159-168.

FREDRICH, F., OHMANN, S., CURIO, B. & KIRSCHBAUM, F. (2003): Spawning migrations of the chub in the River Spree, Germany. - J. Fish Biol. **63**: 710-723.

FRUGET, J.-F., CENTOFANTI, M., DESSAIX, J., OLIVIER, J.-M., DRUART, J.-C. & MARTINEZ, P.-J. (2001): Temporal and spatial dynamics in large rivers: examples of a long-term monitoring of the middle Rhône River. - Annales de Limnologie **37**: 237-251.

FRY, F.E.J. (1971): The effect of environmental factors on the physiology of fish. - In: HOAR, W.S. & RANDALL, D.J. (Eds), Fish Physiology. Academic Press, pp. 1-98.

GARNER, P. (1996): Microhabitat use and diet of 0+ cyprinid fishes in a lentic, regulated reach of the River Great Ouse, England. - J. Fish Biol. **48**: 367-382.

GILLET, A. (1989): Le déroulement de la fraie des principaux poissons lacustres. - Hydroécol. Appl. **1/2**: 117-143.

GINOT, V., SOUCHON, Y. & ROGER, P. (1996): Impact de l'élévation artificielle de température induite par le fonctionnement du Centre Nucléaire du Bugey (fleuve Rhône) sur les communautés de poissons. - Hydroécol. Appl. **8**: 1-33.

GOLOVANOV, V. & SMIRNOV, A. (2011): Specific features of thermoregulation behavior in early juveniles of roach "Rutilus rutilus" under thermogridient conditions. - J. Ichthyol. **51**: 466-473.

GRAHAM, C.T. & HARROD, C. (2009): Implications of climate change for the fishes of the British Isles. - J. Fish Biol. **74**: 1143-1205.