<table><thead><tr><th></th><th></th><th>Juin</th><th>Juillet</th><th>Août</th><th>Sept.</th><th>Déficit étage</th></tr></thead><tbody><tr><td rowspan="3">Hérault amort Arre (BV H1)</td><td>Prélèvés</td><td>56 000</td><td>68 000</td><td>53 000</td><td>25 000</td><td></td></tr><tr><td>Vpn</td><td>618 736</td><td>327 434</td><td>231 948</td><td>231 260</td><td></td></tr><tr><td>Bilan</td><td>562 736</td><td>259 434</td><td>178 948</td><td>206 260</td><td></td></tr><tr><td></td><td>Réduction prélèvements</td><td>0,0%</td><td>0,0%</td><td>0,0%</td><td>0,0%</td><td></td></tr><tr><td rowspan="3">Arre (BV A)</td><td>Prélèvés</td><td>116 500</td><td>120 500</td><td>108 500</td><td>60 500</td><td></td></tr><tr><td>Vpn</td><td>846 464</td><td>352 082</td><td>108 500</td><td>150 726</td><td></td></tr><tr><td>Bilan</td><td>729 964</td><td>231 582</td><td>0</td><td>90 226</td><td></td></tr><tr><td></td><td>Réduction prélèvements</td><td>0,0%</td><td>0,0%</td><td>0,0%</td><td>0,0%</td><td></td></tr><tr><td rowspan="3">Amort confluence Hérault-Vis (BV H2)</td><td>Prélèvés</td><td>20 000</td><td>27 000</td><td>19 000</td><td>8 000</td><td></td></tr><tr><td>Vpn</td><td>1 221 060</td><td>488 027</td><td>125 252</td><td>237 923</td><td></td></tr><tr><td>Bilan</td><td>1 201 060</td><td>461 027</td><td>106 252</td><td>229 923</td><td></td></tr><tr><td></td><td>Réduction prélèvements</td><td>0,0%</td><td>0,0%</td><td>0,0%</td><td>0,0%</td><td></td></tr><tr><td rowspan="3">Vis (BV V)</td><td>Prélèvés</td><td>137 000</td><td>137 000</td><td>138 000</td><td>137 000</td><td></td></tr><tr><td>Vpn</td><td>5 333 597</td><td>2 616 324</td><td>1 335 457</td><td>1 155 918</td><td></td></tr><tr><td>Bilan</td><td>5 196 597</td><td>2 479 324</td><td>1 197 457</td><td>1 018 918</td><td></td></tr><tr><td></td><td>Réduction prélèvements</td><td>0,0%</td><td>0,0%</td><td>0,0%</td><td>0,0%</td><td></td></tr><tr><td rowspan="3">Hérault entre H2 et H3 Moulin Bertrand (BV H3)</td><td>Prélèvés</td><td>122 500</td><td>117 500</td><td>117 500</td><td>82 500</td><td></td></tr><tr><td>Vpn</td><td>9 465 629</td><td>4 901 682</td><td>2 416 975</td><td>2 156 100</td><td></td></tr><tr><td>Bilan</td><td>9 343 129</td><td>4 784 182</td><td>2 299 475</td><td>2 073 600</td><td></td></tr><tr><td></td><td>Réduction prélèvements</td><td>0,0%</td><td>0,0%</td><td>0,0%</td><td>0,0%</td><td></td></tr><tr><td rowspan="3">Hérault entre H3 et H4 Pont du Diable (BV H4)</td><td>Prélèvés</td><td>1 556 070</td><td>1 814 780</td><td>1 769 590</td><td>1 526 410</td><td></td></tr><tr><td>Vpn</td><td>10 412 658</td><td>4 819 245</td><td>2 147 860</td><td>2 059 024</td><td></td></tr><tr><td>Bilan</td><td>8 856 588</td><td>3 004 465</td><td>378 270</td><td>532 614</td><td></td></tr><tr><td></td><td>Réduction prélèvements</td><td>0,0%</td><td>0,0%</td><td>0,0%</td><td>0,0%</td><td></td></tr><tr><td rowspan="3">Hérault entre H4 et H5 (amont Lergue BV H5)</td><td>Prélèvés</td><td>134 500</td><td>137 500</td><td>140 500</td><td>133 639</td><td></td></tr><tr><td>Vpn</td><td>9 640 078</td><td>3 304 622</td><td>944 020</td><td>1 429 639</td><td></td></tr><tr><td>Bilan</td><td>9 505 578</td><td>3 167 122</td><td>803 520</td><td>1 296 000</td><td></td></tr><tr><td></td><td>Réduction prélèvements</td><td>0,0%</td><td>0,0%</td><td>0,0%</td><td>0,0%</td><td></td></tr><tr><td rowspan="3">Lergue amort (BV L1)</td><td>Prélèvés</td><td>238 378</td><td>305 040</td><td>241 719</td><td>137 303</td><td></td></tr><tr><td>Vpn</td><td>2 146 090</td><td>1 231 766</td><td>565 805</td><td>336 887</td><td></td></tr><tr><td>Bilan</td><td>1 907 712</td><td>926 726</td><td>324 086</td><td>199 584</td><td></td></tr><tr><td></td><td>Réduction prélèvements</td><td>0,0%</td><td>0,0%</td><td>0,0%</td><td>0,0%</td><td></td></tr><tr><td rowspan="3">Lergue aval (BV L2)</td><td>Prélèvés</td><td>230 315</td><td>336 349</td><td>240 838</td><td>117 479</td><td></td></tr><tr><td>Vpn</td><td>4 544 886</td><td>3 549 191</td><td>2 142 448</td><td>1 390 264</td><td></td></tr><tr><td>Bilan</td><td>4 314 571</td><td>3 212 842</td><td>1 901 610</td><td>1 272 785</td><td></td></tr><tr><td></td><td>Réduction prélèvements</td><td>0,0%</td><td>0,0%</td><td>0,0%</td><td>0,0%</td><td></td></tr><tr><td rowspan="3">Hérault entre confluence Lergue et H6 (BV H6)</td><td>Prélèvés</td><td>287 616</td><td>601 437</td><td>551 978</td><td>227 050</td><td></td></tr><tr><td>Vpn</td><td>16 676 982</td><td>7 581 007</td><td>3 959 705</td><td>3 008 085</td><td></td></tr><tr><td>Bilan</td><td>16 389 366</td><td>6 979 570</td><td>3 407 727</td><td>2 781 035</td><td></td></tr><tr><td></td><td>Réduction prélèvements</td><td>0,0%</td><td>0,0%</td><td>0,0%</td><td>0,0%</td><td></td></tr><tr><td rowspan="3">Hérault amort Thongue (BV H7)</td><td>Prélèvés</td><td>649 544</td><td>902 424</td><td>806 529</td><td>301 454</td><td></td></tr><tr><td>Vpn</td><td>18 298 194</td><td>8 612 107</td><td>4 694 381</td><td>3 514 367</td><td></td></tr><tr><td>Bilan</td><td>17 648 650</td><td>7 709 683</td><td>3 887 852</td><td>3 212 913</td><td></td></tr><tr><td></td><td>Réduction prélèvements</td><td>0,0%</td><td>0,0%</td><td>0,0%</td><td>0,0%</td><td></td></tr><tr><td rowspan="3">Hérault aval (BV H8)</td><td>Prélèvés</td><td>2 864 239</td><td>4 094 683</td><td>4 102 831</td><td>2 423 921</td><td></td></tr><tr><td>Vpn</td><td>16 054 915</td><td>6 074 164</td><td>3 165 391</td><td>4 789 177</td><td></td></tr><tr><td>Bilan</td><td>13 190 676</td><td>1 979 482</td><td>-937 440</td><td>2 365 256</td><td></td></tr><tr><td></td><td>Réduction prélèvements</td><td>0,0%</td><td>0,0%</td><td>23,0%</td><td>0,0%</td><td>-937 000</td></tr></tbody></table>

Département de l'Hérault

## Bassin versant HERAULT

### Tâbleau 1

Volumes prélevables mensuels nets à l'étage (tous usages) pour le fleuve Hérault

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

Synthèse des données issues de l'étage d'élaboration du schéma directeur de la ressource en eau / Détermination des volumes prélevables du bassin versant du fleuve Hérault (CEREG, novembre 2015)

<table><tr><td>Prélèvés :</td><td>Volume net actuellement prélevé (m³) sur le sous-bassin contrôlé par le point nodal (estimation étude sur la période 2007-2011)</td></tr><tr><td>Vpn :</td><td>Volume prélevable net (m3) théorique au point nodal : écart entre le débit renaturisé en sécheresse quinquennale et la somme du débit biologique et des prélèvements actuels à l'amont du sous bassin versant. Le Vpn du sous bassin diminue quand les prélèvements sur l'amont augmentent. Pour un même mois, les Vpn des différents sous bassins ne s'ajoutent pas car ils sont liés.</td></tr><tr><td>Bilan (m3) entre la ressource disponible sur le sous bassin versant et les prélèvements actuels, négatif en cas de déficit sur le mois. Le bilan s'améliore quand celui du sous bassin amont s'améliore. Pour un même mois, les bilans des différents sous bassins ne s'ajoutent pas car ils sont liés</td><td>Bilan (m3) entre la ressource disponible sur le sous bassin versant et les prélèvements actuels, négatif en cas de déficit sur le mois. Le bilan s'améliore quand celui du sous bassin amont s'améliore. Pour un même mois, les bilans des différents sous bassins ne s'ajoutent pas car ils sont liés</td></tr><tr><td>réduction prélèvements :</td><td>Taux de réduction (%) nécessaire des prélèvements nets actuels pour les mois structurellement déficitaires</td></tr><tr><td>Déficit étage :</td><td>Cumul des volumes nets marquants (m³) des mois déficitaires de chaque sous bassin versant. Le déficit à l'étage cumule le déficit des sous bassins versants situés à l'amont</td></tr></table>

DREAL LR - décembre 2015
<table><thead><tr><th colspan="10">Points nodaux sur le fleuve Hérault</th></tr><tr><th colspan="10">SDAGE</th></tr></thead><tbody><tr><td>Hérault<br>amont<br>Aire (H1)</td><td>Arre<br>(A)</td><td>Hérault-<br>amont<br>Vis<br>(H2)</td><td>Vis (V)</td><td>Hérault<br>Moulin<br>Bertrand<br>(H3)</td><td>Gorges<br>Hérault, amont<br>prise canal<br>Gignac</td><td>Hérault Pont<br>du Diable (H4)</td><td>Hérault<br>amont<br>Ligue<br>(H5)</td><td>Ligue<br>amont<br>(BV L1)</td><td>Ligue<br>aval (BV<br>L2)</td><td>Hérault (H6)<br>contrôlé par<br>la station<br>d'Aspiran</td><td>Hérault<br>amont<br>Thongue<br>(H7)</td><td>Hérault aval<br>(H8) contrôlé<br>par la station<br>d'Agde</td></tr><tr><td>QMNA5 naturel</td><td>0,33</td><td>0,55</td><td>0,90</td><td>1,40</td><td>3,30</td><td>***</td><td>3,60</td><td>3,8</td><td>0,80</td><td>1,20</td><td>4,10</td><td>4,200</td><td>3,60</td></tr><tr><td>QMNA5 influencé</td><td>0,30</td><td>0,50</td><td>0,82</td><td>1,40</td><td>3,30</td><td>***</td><td>2,10</td><td>2,80</td><td>0,70</td><td>1,10</td><td>4,00</td><td>4,00</td><td>1,90</td></tr><tr><td>1/10 module</td><td>0,27</td><td>0,22</td><td>0,53</td><td>1,15</td><td>2,34</td><td>***</td><td>2,53</td><td>2,99</td><td>0,53</td><td>0,93</td><td>4,04</td><td>4,40</td><td>4,54</td></tr></tbody></table>

<table><tbody><tr><td rowspan="3">QMNA5 naturel</td><td>Juillet</td><td>0,35</td><td>0,5</td><td>0,95</td><td>1,1</td><td>2,5</td><td>3,0</td><td>2**</td><td>2,5</td><td>0,65</td><td>0,8</td><td>3,5</td><td>3,5</td><td>2,25**</td></tr><tr><td>août</td><td>0,30</td><td>0,5</td><td>0,90</td><td>1,1</td><td>2,5</td><td>3,0</td><td>2**</td><td>2,5</td><td>0,65</td><td>0,8</td><td>3,5</td><td>3,5</td><td>2,25**</td></tr><tr><td>Sept.</td><td>0,25</td><td>0,5</td><td>0,85</td><td>1,1</td><td>2,5</td><td>3,0</td><td>2**</td><td>2,5</td><td>0,65</td><td>0,8</td><td>3,5</td><td>3,5</td><td>2,25**</td></tr></tbody></table>

<table><tbody><tr><td rowspan="3">Débit biologique (DB)</td><td>Juillet</td><td>0,35</td><td>0,5</td><td>0,95</td><td>1,1</td><td>2,5</td><td>3,0</td><td>2**</td><td>2,5</td><td>0,65</td><td>0,8</td><td>3,5</td><td>3,5</td><td>2,25**</td></tr><tr><td>août</td><td>0,30</td><td>0,5</td><td>0,90</td><td>1,1</td><td>2,5</td><td>3,0</td><td>2**</td><td>2,5</td><td>0,65</td><td>0,8</td><td>3,5</td><td>3,5</td><td>2,25**</td></tr><tr><td>Sept.</td><td>0,25</td><td>0,5</td><td>0,85</td><td>1,1</td><td>2,5</td><td>3,0</td><td>2**</td><td>2,5</td><td>0,65</td><td>0,8</td><td>3,5</td><td>3,5</td><td>2,25**</td></tr></tbody></table>

<table><tbody><tr><td rowspan="3">DOE et DG d'étage</td><td>Juillet</td><td>≥ 0,35</td><td>≥ 0,5</td><td>0,95 *</td><td>≥ 1,1</td><td>≥ 2,5</td><td>≥ 3</td><td>≥ 2**</td><td>≥ 2,5</td><td>≥ 0,65</td><td>≥ 0,8</td><td>≥ 3,5</td><td>≥ 3,5</td><td>2,25**</td></tr><tr><td>août</td><td>≥ 0,30</td><td>≥ 0,5</td><td>0,9 *</td><td>≥ 1,1</td><td>≥ 2,5</td><td>≥ 3</td><td>≥ 2**</td><td>≥ 2,5</td><td>≥ 0,65</td><td>≥ 0,8</td><td>≥ 3,5</td><td>≥ 3,5</td><td>2,25**</td></tr><tr><td>Sept.</td><td>≥ 0,25</td><td>≥ 0,5</td><td>0,85 *</td><td>≥ 1,1</td><td>≥ 2,5</td><td>≥ 3</td><td>≥ 2**</td><td>≥ 2,5</td><td>≥ 0,65</td><td>≥ 0,8</td><td>≥ 3,5</td><td>≥ 3,5</td><td>2,25**</td></tr></tbody></table>

* en lien avec l'hypothèse de déconnexion entre l'amont et l'aval de H2 justifiée par d'importants apports karstiques localisés entre H2 et H4 => DG = DB

** débits cibles tenant compte des prélèvements structurants du canal de Gignac et du Syndicat du Bas Languedoc à réexaminer à l'échéance du prochain SDAGE

*** chronique des données hydrométriques en cours d'acquisition

<table><tr><td>Débit biologique (DB) :</td><td>Indicateur moyen mensuel (m's) de soutien du bon état des eaux aux points nodaux (fonctionnement biologique satisfaisant du milieu aquatique sur le tronçon amort) – permet de dimensionner les prélèvements admissibles sur les sous-bassins concernés (volumes prélevables nets)</td></tr></table>

<table><tr><td>Débit Objectif d'étage (DOE) et débit de gestion (DG) aux points nodaux :</td><td>Indicateur moyen mensuel (m's) de la gestion structurelle équilibrée aux points nodaux de fermeture des sous bassins en période d'étage, toujours supérieur ou égal au Débit Biologique (DB). Il dépendra de la répartition concertée des volumes prélevables entre les acteurs</td></tr></table>

<table><tr><td>Etiage :</td><td>Mois pour lesquels les basses eaux coïncident avec le pic annuel de prélèvements (Juin à septembre pour le bassin versant du fleuve Hérault)</td></tr><tr><td>QMNA5 :</td><td>Débit moyen mensuel minimum quinquennal (m3/s)</td></tr></table>

Département de l'Hérault

Bassin versant du fleuve Hérault

Débits d'étage de référence (m3/s) aux
points nodaux

Tableau 2

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

Synthèse des données issues de l'étude d'élaboration du schéma directeur de la ressource en eau /Détermination des volumes prélevables du bassin versant du fleuve Hérault (CEREG, novembre 2015)
# Gestion quantitative des ressources en eaux superficielles du bassin versant du Vidourle

## Ecart aux volumes prélevables nets (étude de détermination des volumes prélevables Ginger, octobre 2012)

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

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

DREAL LR / SN / EMA - avril 2014
Bassin versant du VIDOURLE

Tableau 1

Volumes prélevables mensuels nets en étiage,
tous usages (en milliers de m3)

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

Liberté • Égalité • Fraternité
RÉPUBLIQUE FRANÇAISE

Détection régionale
de l'environnement,
et du logement
LANGUAGOD
ROISSILLOU

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

Synthèse des données issues de l'étude de détermination
de volumes prélevables pour le bassin versant du
Vidourle (GINGER rapports FL34 A 0055 – octobre
2012)

DREAL LR / SN / EMA - avril 2014

<table><thead><tr><th></th><th></th><th>Juillet</th><th>Août</th><th>Septembre</th></tr></thead><tbody><tr><td rowspan="3">Bilan bassin au point nodal V6</td><td>Prélévé</td><td>510</td><td>426</td><td>180</td></tr><tr><td>Vpn</td><td>321</td><td>80</td><td>544</td></tr><tr><td>Bilan</td><td>-189</td><td>-346</td><td>364</td></tr><tr><td></td><td>Réduction prélèvements</td><td>37%</td><td>81%</td><td>0%</td></tr></tbody></table>

<table><thead><tr><th></th><th></th><th>Juillet</th><th>Août</th><th>Septembre</th></tr></thead><tbody><tr><td rowspan="4">Bilan bassin au point nodal V5<br>SDAGE</td><td>Prélévé</td><td>369</td><td>330</td><td>179</td></tr><tr><td>Vpn</td><td>295</td><td>80</td><td>492</td></tr><tr><td>Prélévé sous BV6</td><td>141</td><td>96</td><td>&lt;1</td></tr><tr><td>Vpn avec besoins BV 6</td><td>153</td><td>0</td><td>494</td></tr><tr><td></td><td>Bilan</td><td>-216</td><td>-330</td><td>315</td></tr><tr><td></td><td>Réduction prélèvements</td><td>58%</td><td>100%</td><td>0%</td></tr></tbody></table>

<table><thead><tr><th></th><th></th><th>Juillet</th><th>Août</th><th>Septembre</th></tr></thead><tbody><tr><td rowspan="3">Bilan bassin au point nodal V1</td><td>Prélévé</td><td>109</td><td>133</td><td>108</td></tr><tr><td>Vpn</td><td>118</td><td>37</td><td>166</td></tr><tr><td>Bilan</td><td>9</td><td>-95</td><td>58</td></tr><tr><td></td><td>Réduction prélèvements</td><td>0%</td><td>72%</td><td>0%</td></tr></tbody></table>

<table><thead><tr><th rowspan="3">Sous-bassins<br>aux points<br>nodaux des<br>affluents</th><th>CRE<br>Crespenou</th><th>Prélévé</th><th>Août</th><th>Septembre</th></tr></thead><tbody><tr><td rowspan="2">CRI Crieulon</td><td>Vpn</td><td>68</td><td>63</td><td>38</td></tr><tr><td>Bilan</td><td>-54</td><td>-50</td><td>1</td></tr><tr><td></td><td>Réduction prélèvements</td><td>80%</td><td>79%</td><td>3%</td></tr><tr><td></td><td>Prélévé</td><td>&lt;1</td><td>&lt;1</td><td>&lt;1</td></tr><tr><td></td><td>Vpn</td><td>36</td><td>11</td><td>88</td></tr><tr><td></td><td>Bilan</td><td>36</td><td>11</td><td>88</td></tr><tr><td></td><td>Réduction prélèvements</td><td>0%</td><td>0%</td><td>0%</td></tr></tbody></table>

<table><tr><td>Prélévé (en milliers de m3) :</td><td>Volume net cumulé prélévé actuellement sur l'ensemble du bassin aux points nodaux CRE, CRI, V1, V5 et V6 (estimation étude)</td></tr></table>

<table><tr><td>Vpn (en milliers de m3) :</td><td>Volume prélevable net déterminé pour l'ensemble du bassin versant aux points nodaux CRE, CRI, V1, V5 et V6 (écart entre débit naturel et débit biologique en hydrologie quinquennale sèche)</td></tr></table>

<table><tr><td>Bilan (en milliers de m3) :</td><td>Bilan (en milliers m3) entre la ressource disponible et les prélèvements actuels, négatif en cas de déficit sur le mois</td></tr></table>

<table><tr><td>Réduction prélèvements (%):</td><td>Taux de réduction (%) nécessaire des prélèvements actuels pour les mois structurellement déficitaires</td></tr></table>
<table><thead><tr><th rowspan="2"></th><th colspan="6">Points nodaux sur l'axe du Vidourle</th><th colspan="3">Points nodaux affluents</th></tr><tr><th>V0</th><th>V1</th><th>V2</th><th>V3</th><th>V4</th><th>V5<br>SDAGE</th><th>V6</th><th>Crespenou<br>CRE</th><th>Crieulon<br>CRI</th></tr></thead><tbody><tr><th>QMNA5 naturel (l/s)</th><td>30</td><td>100</td><td>150</td><td>150</td><td>190</td><td>200</td><td>200</td><td>20</td><td>40</td></tr><tr><th>QMNA5 influencé (l/s)</th><td>12</td><td>47</td><td>40</td><td>34</td><td>75</td><td>77</td><td>42</td><td>0</td><td>41</td></tr><tr><th>1/10 module naturel (l/s)</th><td>32</td><td>82</td><td>282</td><td>285</td><td>511</td><td>688</td><td>859</td><td>12</td><td>52</td></tr></tbody></table>

<table><tr><th>Débit biologique étiage (l/s)</th><td>36</td><td>86</td><td>135</td><td>140</td><td>170</td><td>170</td><td>170</td><td>15</td><td>36</td></tr><tr><th>DOE juillet (l/s)</th><td></td><td></td><td></td><td></td><td></td><td>220</td><td></td><td></td><td></td></tr><tr><th>DOE août (l/s)</th><td></td><td></td><td></td><td></td><td></td><td>200</td><td></td><td></td><td></td></tr><tr><th>DOE septembre (l/s)</th><td></td><td></td><td></td><td></td><td></td><td>170</td><td></td><td></td><td></td></tr><tr><th>DCR étiage (l/s)</th><td></td><td></td><td></td><td></td><td></td><td>170</td><td></td><td></td><td></td></tr></table>

<table><tr><th>Débit biologique :</th><td>Indicateur moyen mensuel de soutien du bon état des eaux aux points nodaux (fonctionnement biologique satisfaisant du milieu aquatique sur le tronçon amont) – permet de dimensionner les prélèvements admissibles sur les sous-bassins concernés (volumes prélevables nets)</td></tr></table>

<table><tr><th>DOE (Débit Objectif d'Etiage) :</th><td>Indicateur moyen mensuel de la gestion structurelle équilibrée au Point Stratégique de Référence du SDAGE 2016-2021 (satisfaction du bon état des eaux et, en moyenne 8 années sur 10, des usages), déterminé sur la base des volumes prélevables et de l'identification des réductions de prélèvements envisageables (scénario 1 de l'étude) – contrôlé a posteriori par l'Etat à la station hydrométrique associée</td></tr></table>

<table><tr><th>DCR (Débit de Crise Renforcée) :</th><td>Indicateur seuil de gestion de crise au Point Stratégique de Référence du SDAGE 2016-2021 – contrôlé en valeur journalière par l'Etat à la station hydrométrique associée lorsque les actions nécessaires au respect des DOE auront été mises en place</td></tr></table>

<table><tr><th>Etiage :</th><td>Mois pour lesquels les basses eaux coincident avec le pic annuel de prélèvement (juillet à septembre pour le bassin versant du Vidourle)</td></tr><tr><th>QMNA5 :</th><td>Débit moyen mensuel minimum quinquennal</td></tr></table>

Tableau 2

Débits d'étiage de référence aux points nodaux (en l/s)
Annexe 4 : Liste des stations d'épuration

<table><thead><tr><td>Nom du système</td><td>Code SANDRE</td><td>Commune</td><td>Taille de l'agglomération</td><td>Débit (m3/j)</td><td>Milieu récepteur</td><td>Nom du milieu</td></tr></thead><tbody><tr><td>VIOLS LE FORT</td><td>60934343002</td><td>Viols-le-Fort</td><td>850</td><td>105</td><td>Eau douce de surface</td><td>L'Arnède</td></tr><tr><td>SAUTEYRARGUES (Bourg)</td><td>60934297002</td><td>Sauteyrargues</td><td>180</td><td>23</td><td>Sol</td><td>Sol</td></tr><tr><td>CAUSSE DE LA SELLE</td><td>60934060001</td><td>Causse-de-la-Selle</td><td>500</td><td>25</td><td>Sol</td><td>Sol</td></tr><tr><td>VILLENEUVE LES MAGUELONE</td><td>60934337001</td><td>Villeuneuves-les-Maguelone</td><td>10270</td><td>1611</td><td>Eau douce de surface</td><td>La Mosson</td></tr><tr><td>VALFAUNES - LANCYRE</td><td>60934322003</td><td>Valflaunès</td><td>80</td><td>12</td><td>Sol</td><td>Sol</td></tr><tr><td>SAUTEYRARGUES (Vabre)</td><td>60934297001</td><td>Sauteyrargues</td><td>60</td><td>9</td><td>Sol</td><td>Sol</td></tr><tr><td>CAZEVIEILLE</td><td>60934066001</td><td>Cazevieille</td><td>70</td><td>13</td><td>Sol</td><td>Sol</td></tr><tr><td>MURVIEL LES MONTPELIER</td><td>60934179002</td><td>Murviel-lès-Montpelier</td><td>1680</td><td>291</td><td>Eau douce de surface</td><td>Ruisseau de Lassedéron</td></tr><tr><td>ST JEAN DE CORNIES</td><td>60934265001</td><td>Saint-Jean-de-Cornies</td><td>600</td><td>85</td><td>Eau douce de surface</td><td>Le Budel</td></tr><tr><td>MAS DE LONDRES</td><td>60934152001</td><td>Mas-de-Londres</td><td>480</td><td>56</td><td>Eau douce de surface</td><td>Ruisseau de la Fontaine d'Azémar</td></tr><tr><td>ST GENIES DES MOURGUES</td><td>60934256001</td><td>Saint-Geniès-des-Mourgues</td><td>1110</td><td>570</td><td>Eau douce de surface</td><td>Canal de Lansargues</td></tr><tr><td>VALFLAUNES-VALCYRE</td><td>60934322002</td><td>Valflaunès</td><td>100</td><td>29</td><td>Sol</td><td>Le Rieufrêche</td></tr><tr><td>LAURET</td><td>60934131001</td><td>Lauret</td><td>300</td><td>34</td><td>Sol</td><td>Le Brestalou</td></tr><tr><td>BEAULIEU-RESTIN-CLIERES</td><td>60934227002</td><td>Restinclières</td><td>3460</td><td>479</td><td>Rejet diffus</td><td>Le Dardaillon-Ouest</td></tr><tr><td>ST CLEMENT(Rouargues)</td><td>60934247004</td><td>Saint-Clément-de-Rivière</td><td>4896</td><td>978</td><td>Eau douce de surface</td><td>Le lez</td></tr><tr><td>Bonnel-Plagnol-ZAE-CresseStMartin</td><td>60934088002</td><td>Cournonsec</td><td>240</td><td>47</td><td>Sol</td><td>Ruisseau de Vire</td></tr><tr><td>LAVERUNE</td><td>60934134003</td><td>Laverune</td><td>4039</td><td>513</td><td>Eau douce de surface</td><td>Ruisseau de Lassedéron</td></tr></tbody></table>
<table><tbody><tr><td>SAINT D'ORQUES</td><td>GEORGES</td><td>60934259002</td><td>Saint-Georges-d'Orques</td><td>6345</td><td>766</td><td>Eau douce de surface</td><td>Ruisseau de Lassedéron</td></tr><tr><td>VACQUIERES</td><td></td><td>60934318001</td><td>Vacquières</td><td>220</td><td>271</td><td>Eau douce de surface</td><td>Le Brestalou</td></tr><tr><td>MURLES</td><td></td><td>60934177001</td><td>Murles</td><td>120</td><td>19</td><td>Eau douce de surface</td><td>Ruisseau de Saint-Jean</td></tr><tr><td>FONTANES - bourg</td><td></td><td>60934102001</td><td>Fontanès</td><td>80</td><td>29</td><td>Eau douce de surface</td><td>La Bénovie</td></tr><tr><td>LES MATELLES</td><td></td><td>60934153002</td><td>Les Matelles</td><td>1927</td><td>492</td><td>Eau douce de surface</td><td>Le Lirou</td></tr><tr><td>MAERA (Ex CEREIREDE)</td><td></td><td>60934172001</td><td>Montpellier</td><td>375145</td><td>93365</td><td>Eau côtière</td><td>Le lez</td></tr><tr><td>SAINT-DREZERY</td><td></td><td>60934249002</td><td>Saint-Drézéry</td><td>2180</td><td>303</td><td>Eau douce de surface</td><td>Le Bérange</td></tr><tr><td>BUZIGNARGUES</td><td></td><td>60934043001</td><td>Buzygnargues</td><td>140</td><td>21</td><td>Eau douce de surface</td><td>La Bénovie</td></tr><tr><td>SAINT HILAIRE DE BAUVOR- bourg</td><td></td><td>60934263002</td><td>Saint-Hilaire-de-Beauvoir</td><td>300</td><td>41</td><td>Eau douce de surface</td><td>Le Budel</td></tr><tr><td>Nom du système</td><td></td><td>Code SANDRE</td><td>Commune</td><td>Taille de l'agglomération</td><td>Débit (m3/j)</td><td>Milieu récepteur</td><td>Nom du milieu</td></tr><tr><td>ST VINCENT DE BARBEYRARGUES-bourg</td><td></td><td>60934290002</td><td>Saint-Vincent-de-Barbeyrargues</td><td>360</td><td>1224</td><td>Eau douce de surface</td><td>Ruisseau de Cassagnoles</td></tr><tr><td>COURNONTERRAL</td><td></td><td>60934088001</td><td>Cournonterral</td><td>10800</td><td>45</td><td>Eau douce de surface</td><td>Ruisseau du Coulazou</td></tr><tr><td>GUZARGUES</td><td></td><td>60934118001</td><td>Guzargues</td><td>330</td><td>36</td><td>Eau douce de surface</td><td>La Cadoule</td></tr><tr><td>NOTRE-DAME-DE-LONDRES</td><td></td><td>60934185002</td><td>Notre-Dame-De-Londres</td><td>350</td><td>10</td><td>Eau douce de surface</td><td>La Tourguille</td></tr><tr><td>ST JEAN DE CUCULLES</td><td></td><td>60934266001</td><td>Saint-Jean-de-Cuculles</td><td>150</td><td>21</td><td>Eau douce de surface</td><td>Ruisseau de Yorgues</td></tr><tr><td>ST JEAN DE BUEGES</td><td></td><td>60934264001</td><td>Saint-Jean-de-Buèges</td><td>140</td><td>123</td><td>Eau douce de surface</td><td>La Buèges</td></tr><tr><td>CLARET</td><td></td><td>60934078001</td><td>Claret</td><td>840</td><td>0</td><td>Eau douce de surface</td><td>Le Brestalou</td></tr><tr><td>PEGAIROLLES DE BUEGES MEJANEL</td><td></td><td>60934195001</td><td>Pégairolles-de-Buèges</td><td>60</td><td>664</td><td>Eau douce de surface</td><td>Ruisseau des Vignes Coudoulières</td></tr><tr><td>SAINT MATHIEU DE TREVIERS</td><td></td><td>60934276001</td><td>Saint-Mathieu-de-Tréviers</td><td>4472</td><td>1448</td><td>Eau douce de surface</td><td>Le Terrieu</td></tr><tr><td>BAILLARGUES - ST BRES</td><td></td><td>60934022002</td><td>Baillargues</td><td>16656</td><td>71</td><td>Eau douce de surface</td><td>Ruisseau d'Aigues-Vives</td></tr><tr><td>VAILHAUQUES-Bel air</td><td></td><td>60934320003</td><td>Vailhauquès</td><td>270</td><td>39</td><td>Eau douce de surface</td><td>Le Rieu de Querelle</td></tr><tr><td>VALFLAUNES-bourg</td><td></td><td>60934322004</td><td>Valflaunès</td><td>190</td><td>2763</td><td>Eau douce de surface</td><td>Le Terrieu</td></tr></tbody></table>
<table><tr><td>FABREGUES</td><td>60934095003</td><td>Fabrègues</td><td>18503</td><td>116</td><td>Eau douce de surface</td><td>Ruisseau du Coulazou</td></tr><tr><td>MONTAUD</td><td>60934164001</td><td>Montaud</td><td>450</td><td>171</td><td>Eau douce de surface</td><td>Ruisseau de la Lequette</td></tr><tr><td>COMBAILLAUX-LES SAJOLLES</td><td>60934082003</td><td>Combaillaux</td><td>1161</td><td>102</td><td>Eau douce de surface</td><td>Ruisseau de Miege Sole</td></tr><tr><td>SAINTE CROIX DE QUIN-TILLARGUES (nouvelle)</td><td>60934248002</td><td>Sainte-Croix-de-Quintillargues</td><td>480</td><td>540</td><td>Eau douce de surface</td><td>La Bénovie</td></tr><tr><td>SAINT MARTIN DE LONDRES</td><td>60934274001</td><td>Saint-Martin-de-Londres</td><td>2541</td><td>64</td><td>Eau douce de surface</td><td>Ruisseau de Rieutort</td></tr><tr><td>ST BAUZILLE DE MONT-MEL-BOURG</td><td>60934242001</td><td>Saint-Bauzille-de-Montmel</td><td>280</td><td>2175</td><td>Eau douce de surface</td><td>Le Brau</td></tr><tr><td>SAINT GELY DU FESC</td><td>60934255001</td><td>Saint-Gély-du-Fesc</td><td>13400</td><td>370</td><td>Eau douce de surface</td><td>Ruisseau de Pézouillet</td></tr><tr><td>VAILHAUQUES</td><td>60934320002</td><td>Vailhauquès</td><td>2254</td><td>41</td><td>Eau douce de surface</td><td>La Mosson</td></tr><tr><td>TRIADOU</td><td>60934314001</td><td>Le Triadou</td><td>270</td><td>21</td><td>Eau douce de surface</td><td>Le Terrieu</td></tr><tr><td>SAINT BAUZILLE DE MONTMEL (Favas)</td><td>60934242002</td><td>Saint-Bauzille-de-Montmel</td><td>100</td><td>22</td><td>Eau douce de surface</td><td>Le Brau</td></tr><tr><td>CLARET (les Embrus-calles)</td><td>60934078002</td><td>Claret</td><td>840</td><td>105</td><td>Eau douce de surface</td><td>Ruisseau de la Font</td></tr></table>
# Annexe 5 : Point sur la réglementation autour de la REUT

Sur le plan réglementaire, c'est l'arrêté du 2 août 2010 relatif à l'utilisation d'eaux issues du traitement d'épuration des eaux résiduaires urbaines pour l'irrigation de cultures ou d'espaces verts modifié par l'arrêté du 25 juin 2014 qui représente le texte de référence. L'instruction interministérielle N° DGS/EA4/DEB/DGPE/2016/135 du 26 avril 2016 en précise les modalités d'application.

Par la suite sont précisés les points ayant une importance dans l'élaboration de la méthodologie d'analyse spatiale du projet SoPoLo (tâche 1). Enfin nous évoquerons les évolutions potentielles de la réglementation sur la base de questions posées à l'agence régionale de santé (ARS).

## A. Stations concernées

Selon le premier article dudit arrêté, les eaux usées traitées concernées par la REUT sont celles issues de systèmes de traitements collectifs et non-collectifs (la charge en pollution organique devant être supérieure à 1.2kg de DBO5/j)

Sont interdites à la réutilisation les eaux usées traitées issues de stations de traitement reliées à des « établissement de collecte, d'entreposage, de manipulation après collecte ou de transformation des sous-produits animaux de catégorie 1 ou 2 au sens du règlement (CE) n° 1069/2009 [...] à l'exception des cas où les eaux sont, préalablement à leur rejet dans le réseau de collecte, traitées thermiquement à 133° C pendant 20 minutes sous une pression de 3 bars ».

Pour rappel, la classification des sous-produits animaux obéit à une logique de risque à l'encontre de la santé humaine ou animale. Les catégories 1 et 2 présentent des risques tandis que la catégorie 3 rassemble les « sous-produits animaux issus d'animaux sains, ou des sous-produits animaux au travers desquels ne peuvent être véhiculées de maladies transmissibles à l'homme ou à l'animal » (Guide de classification des sous-produits animaux et de leurs devenirs du MAAF)

Les stations produisant des boues ne respectant pas les contraintes énoncées au titre de l'annexe I de l'arrêté du 8 janvier 1998 relatif aux prescriptions techniques applicables aux épandages de boues sur les sols agricoles ne peuvent être retenues dans le cadre d'un projet de REUT. Pour rappel cet arrêté fixe des seuils de concentration en éléments traces métalliques (ETM) à ne pas dépasser dans les boues d'épandages ou les sols. Ainsi les sols pollués ne peuvent pas non plus être sujet à de la REUT. Toutefois ne disposant pas d'informations suffisantes sur ces pollutions (délimitation des sols pollués), ce critère n'est pas retenu.

## B. Limitations selon les zones irriguées

Comme mentionné auparavant, les sols pollués par des ETM au titre de l'arrêté du 8 janvier 1998 ne peuvent être retenus pour un projet de REUT.

Les zones situées à l'intérieur de périmètres de protection rapprochés de captages d'eau potables ne peuvent pas être retenues pour un projet de REUT. Toutefois elles peuvent être sujettes à une dérogation dans le cadre d'un captage de surface ou karstique après avis d'un hydrogéologue agréé en matière d'hygiène publique.

Les maires ou préfets peuvent également délimiter des périmètres dans lesquels la REUT est interdite car elle aurait un « impact sanitaire sur un usage sensible de l'eau, tel qu'un captage public utilisé pour la consommation humaine, un site de conchyliculture, de pisciculture, de cressiculture, de pêche à pied, de baignade ou d'activités nautiques et, en cas d'absence de réseau public d'eau potable, un puits
ou un forage réalisé à des fins domestiques de l'eau et ayant fait l'objet d'une déclaration auprès du maire de la commune concerné conformément aux dispositions de l'article L 2224-9 du CGCT ».

A ces délimitations réglementaires se rajoutent des contraintes de terrains dont les modalités sont synthétisées dans le tableau suivant :

<table><thead><tr><th rowspan="2"></th><th rowspan="2"></th><th rowspan="2"></th><th colspan="2">Saturé</th><th colspan="3">Non saturé</th></tr><tr><th></th><th>&lt;3% de pente</th><th>Entre 3% et 7% de pente</th><th>Plus de 7% de pente</th></tr></thead><tbody><tr><td rowspan="4">Sans couvert végétal</td><td rowspan="2">Sol fin</td><td>Karstique</td><td>Non</td><td>Non</td><td>Non</td><td>Non</td><td>Non</td></tr><tr><td>Non karstique</td><td>Non</td><td>Oui</td><td>Oui</td><td>Localisé</td><td>Localisé</td></tr><tr><td rowspan="2">Sol épais (&gt;1m)</td><td>Karstique</td><td>Non</td><td>Non</td><td>Non</td><td>Non</td><td>Non</td></tr><tr><td>Non karstique</td><td>Non</td><td>Oui</td><td>Oui</td><td>Localisé</td><td>Localisé</td></tr><tr><td rowspan="4">Avec couvert végétal</td><td rowspan="2">Sol fin</td><td>Karstique</td><td>Non</td><td>Non</td><td>Non</td><td>Non</td><td>Non</td></tr><tr><td>Non karstique</td><td>Non</td><td>Oui</td><td>Oui</td><td>Oui</td><td>Oui</td></tr><tr><td rowspan="2">Sol épais (&gt;1m)</td><td>Karstique</td><td>Non</td><td>Eau A ou B</td><td>Eau A ou B (localisé)</td><td>Eau A ou B (localisé)</td><td>Eau A ou B (localisé)</td></tr><tr><td>Non karstique</td><td>Non</td><td>Oui</td><td>Oui</td><td>Oui</td><td>Oui</td></tr></tbody></table>

N.B. : Eau A ou B désigne le niveau de qualité de traitement des eaux usées traitées selon les critères énoncés ci-dessous.

### C. Définition de la qualité des eaux usées traitées

Ci-suivent les tableaux directement issus de l'arrêté du 2 août 2010 :

Les eaux usées traitées sont classées dans le niveau de qualité qui correspond au classement du paramètre le plus défavorable.

Les abattements sont mesurés entre les eaux brutes, en entrée de la station de traitement des eaux usées, et les eaux usées traitées, en sortie de la station de traitement des eaux usées ou de la filière de traitement complémentaire, le cas échéant.

Quatre niveaux de qualité sanitaire des eaux usées traitées (A, B, C et D) sont définis comme suit :

<table><thead><tr><th rowspan="2">PARAMÈTRES</th><th colspan="4">NIVEAU DE QUALITÉ SANITAIRE DES EAUX USÉES TRAITÉES</th></tr><tr><th>A</th><th>B</th><th>C</th><th>D</th></tr></thead><tbody><tr><td>Matières en suspension (mg/L)</td><td>&lt; 15</td><td colspan="3">Conforme à la réglementation des rejets d'eaux usées traitées pour l'exutoire de la station hors période d'irrigation</td></tr><tr><td>Demande chimique en oxygène (mg/L)</td><td>&lt; 60</td><td colspan="3"></td></tr><tr><td>Escherichia coli (UFC/100mL)</td><td>≤ 250</td><td>≤ 10 000</td><td>≤ 100 000</td><td>-</td></tr><tr><td>Entérocoques fécaux (abattement en log)</td><td>≥ 4</td><td>≥ 3</td><td>≥ 2</td><td>≥ 2</td></tr><tr><td>Phages ARN F-spécifiques (abattement en log)</td><td>≥ 4</td><td>≥ 3</td><td>≥ 2</td><td>≥ 2</td></tr><tr><td>Spores de bactéries anaérobies sulfito-réductrices (abattement en log)</td><td>≥ 4</td><td>≥ 3</td><td>≥ 2</td><td>≥ 2</td></tr></tbody></table>

Pour les stations traitant des eaux brutes dont les charges en entérocoques fécaux, en phages ARN F-
spécifiques ou en spores de bactéries anaérobies sulfito-réductrices sont inférieures ou égales à 4 log, les critères suivants sont appliqués :

<table><thead><tr><th rowspan="2">Concentration en micro-organismes (*) dans les eaux brutes</th><th colspan="4">Concentration en micro-organismes (*) dans les eaux usées traitées</th></tr><tr><th>10<sup>4</sup></th><th>10<sup>3</sup></th><th>10<sup>2</sup></th><th>&lt;10</th></tr></thead><tbody><tr><th>10<sup>4</sup></th><td>-</td><td>-</td><td>Qualité C ou D (**)</td><td>Qualité A, B, C ou D (**)</td></tr><tr><th>10<sup>3</sup></th><td>-</td><td>-</td><td>-</td><td>Qualité A, B, C ou D (**)</td></tr><tr><th>10<sup>2</sup></th><td>-</td><td>-</td><td>-</td><td>Qualité A, B, C ou D (**)</td></tr><tr><th>10</th><td>-</td><td>-</td><td>-</td><td>Qualité A, B, C ou D (**)</td></tr></tbody></table>

(**) Le niveau de qualité sanitaire des eaux usées traitées doit être déterminé en prenant également en compte les micro-organismes de concentration supérieure à 10⁴ dans les eaux brutes ainsi que les autres paramètres: Escherichia coli, matières en suspension et demande chimique en oxygène, en respectant les valeurs définies à l'annexe II de l'arrêté.

A savoir que les mesures en phages ARN F-spécifiques et en spores de bactéries anaérobies ne sont réalisées que dans le cadre d'un projet de REUT (demande & suivi). Il est donc impossible de prévoir le niveau de qualité sanitaire des EUT dans le cadre d'une étude sauf si des mesures sont mises en place (proximité d'une zone de baignade par exemple).

## D. Distances réglementaires et activités à protéger

Les distances minimales à respecter (en mètres) entre l'irrigation par des eaux usées traitées et les activités à protéger figurent dans les tableaux suivants :

<table><thead><tr><th>CARACTÉRISTIQUES DE L'ASPERSEUR</th><th>DISTANCE ASPERSEUR À ZONE SENSIBLE (1)</th><th></th></tr></thead><tbody><tr><td>Portée</td><td>Avec écran 2 et basse pression (2)</td><td>Dans les autres cas</td></tr><tr><td>Faible portée : &lt; 10 m</td><td>5 m (3)</td><td></td></tr><tr><td>Moyenne portée : 10 à 20 m</td><td>10 m (3)</td><td rowspan="2">Deux fois la portée</td></tr><tr><td>Grande portée : &gt; 20 m</td><td>10 m (3)</td></tr></tbody></table>

(1) Habitations, cours et jardins attenants aux habitations, voies de circulation, lieux publics de passage et de loisir, bâtiments publics et bâtiments d'entreprise, quels que soient le sens et la vitesse du vent dominant.

(2) Dispositif végétalisé arbustif ou écrans fixes ou mobiles tels que murs, brise-vents, canisses, panneaux d'occultation, etc., dont la hauteur doit être au moins égale à celle de l'apogée de l'asperseur.

(3) Cette valeur est augmentée de la portée pour le secteur couvert par l'arrosage.

## NIVEAU DE QUALITÉ SANITAIRE DES EAUX USÉES TRAITÉES

<footnote>(*) Micro-organismes nécessitant un abattement : entérocoques fécaux, phages ARN F-spécifiques et spores de bactéries anaérobies sulfito-réductrices.</footnote>
<table><thead><tr><th>NATURE DES ACTIVITÉS À PROTÉGER</th><th>A</th><th>B</th><th>C ET D</th></tr></thead><tbody><tr><td>Plan d'eau (1)</td><td>20 m</td><td>50 m</td><td>100 m</td></tr><tr><td>Bassin aquacole (à l'exception des coquillages filtreurs)</td><td>20 m</td><td>50 m</td><td>100 m</td></tr><tr><td>Pisciculture y compris pêche de loisir</td><td>50 m</td><td>200 m</td><td>300 m</td></tr><tr><td>Conchyliculture</td><td>50 m</td><td>100 m</td><td>200 m</td></tr><tr><td>Pêche à pied des coquillages filtreurs</td><td>50 m</td><td>100 m</td><td>200 m</td></tr><tr><td>Baignades et activités nautiques</td><td>50 m</td><td>100 m</td><td>200 m</td></tr><tr><td>Abreuvement du bétail</td><td>50 m</td><td>100 m</td><td>200 m</td></tr><tr><td>Cressiculture</td><td>50 m</td><td>200 m</td><td>300 m</td></tr></tbody></table>

<footnote>(1) A l'exception du plan d'eau servant d'exutoire au rejet de la station d'épuration et des plans d'eau privés où l'accès est réglementé et où aucune activité telle que baignade, sport nautique et aquatique, pêche ou abreuvement du bétail n'est pratiquée.</footnote>
## E. Usages autorisés selon le niveau de qualité des EUT

<table><thead><tr><th rowspan="2">TYPE D'USAGE</th><th colspan="4">NIVEAU DE QUALITÉ SANITAIRE DES EAUX USÉES TRAITÉES</th></tr><tr><th>A</th><th>B</th><th>C</th><th>D</th></tr></thead><tbody><tr><td>Cultures maraîchères, fruitières et légumières non transformées par un traitement thermique industriel adapté (excepté cressiculture (1))</td><td>+</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Cultures maraîchères, fruitières, légumières transformées par un traitement thermique industriel adapté</td><td>+</td><td>+</td><td>-</td><td>-</td></tr><tr><td>Pâturage (2)</td><td>+</td><td>+ (3)</td><td>-</td><td>-</td></tr><tr><td>Espaces verts ouverts au public (4)</td><td>+ (5)</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Fleurs vendues coupées</td><td>+</td><td>+ (6)</td><td>-</td><td>-</td></tr><tr><td>Pépinières et arbustes et autres cultures florales</td><td>+</td><td>+</td><td>+ (6)</td><td>-</td></tr><tr><td>Fourrage frais</td><td>+</td><td>+ (3)</td><td>-</td><td>-</td></tr><tr><td>Autres cultures céréalières et fourragères</td><td>+</td><td>+</td><td>+ (6)</td><td>-</td></tr><tr><td>Arboriculture fruitière</td><td>+</td><td>+ (7)</td><td>+ (8)</td><td>-</td></tr><tr><td>Taillis à courte rotation ou à très courte rotation, avec accès contrôlé du public</td><td>+</td><td>+</td><td>+ (6)</td><td>+ (6)</td></tr><tr><td>Forêt, hors taillis à courte rotation avec accès contrôlé du public</td><td>-</td><td>-</td><td>-</td><td>-</td></tr></tbody></table>

+ autorisée, - : interdite.

(1) La réutilisation d'eaux usées traitées est interdite pour la cressiculture.

(2) En cas d'aspersion, les animaux ne doivent pas être au champ au moment de l'opération et les abreuvoirs, au cas où ils seraient arrosés, doivent être rincés avant utilisation.

(3) Sous réserve du respect d'un délai après irrigation de 10 jours en l'absence d'abattoir relié à la station de traitement des eaux usées et de 21 jours dans le cas contraire.

(4) On entend par espace vert, notamment : les aires d'autoroutes, cimetières, golfs, hippodromes, parcs, jardins publics, parties communes de lotissements, ronds-points (pas envisageable d'après l'ARS) et autres terre-pleins, squares, stades, etc.

(5) Irrigation en dehors des heures d'ouverture au public, ou fermeture aux usagers pendant l'irrigation et deux heures suivant l'irrigation dans le cas d'espaces verts fermés ; irrigation pendant les heures de plus faible fréquentation et interdiction d'accès aux passants pendant l'irrigation et deux heures suivant l'irrigation dans le cas d'espaces verts ouverts de façon permanente.

(6) Uniquement par irrigation localisée, telle que définie à l'article 2.

(7) Interdite pendant la période allant de la floraison à la cueillette pour les fruits non transformés, sauf en cas d'irrigation au goutte à goutte.

(8) Uniquement par goutte à goutte.

Dans le cas d'une culture sous serre, seule l'irrigation localisée, telle que définie à l'article 2, est autorisée, en cas de micro-aspersion l'entrée dans les serres est interdite au cours et une heure après l'arrosage.

A titre indicatif ont été calculées les surfaces agricoles compatibles avec les eaux usées de qualité A, B et C sur le territoire d'étude :
<figure><img src="image_47.png" /></figure>

## Précisions apportées par l'ARS

Dans le cadre du projet SoPoLo nous avons été amenés à poser quelques questions (en bleu) sur la réglementation autour de la REUT, les réponses sont en vert.

« L'annexe III de l'arrêté du 2 août 2010 relatif à l'utilisation d'eaux issues du traitement d'épuration des eaux résiduaires urbaines pour l'irrigation de cultures ou d'espaces verts mentionne la possibilité d'arrosage d'espaces verts ouverts au public (selon des conditions précisées dans la même annexe). De quel type d'irrigation est-il question ici?

o S'il s'agit d'une irrigation par aspersion, n'y-a-t-il pas un conflit avec l'annexe I du même arrêté qui impose au moins 5 à 10m de distance par rapport aux lieux publics? (Le lieu de l'irrigation est donc également le lieu à éviter). »

« Non pas obligatoirement, ça dépend du type d'espace public, de sa taille et de son type de fréquentation. C'est à voir au cas par cas

Le conflit peut se résoudre de 2 façons :

- Par une bande de non aspersion dans espaces verts vis-à-vis des lieux fréquentés par le public

- Aspersion pendant fermeture accès au public

Dans tous les cas l'objectif est d'éviter le contact entre l'eau et le public ainsi que la plante arrosée pendant quelques heures. »

« Est-il possible d'utiliser les eaux usées traitées dans le cadre de nettoyage, par exemple de matériel agricole (cuves) ou de la voirie?

o Si oui, selon quelles conditions?

o Si non, peut-on s'attendre à une évolution de la législation en ce sens? »

« Actuellement la réglementation ne traite pas de cet usage, l'utilisation d'eaux usées traitées pour des usages de nettoyage dans l'espace public (voiries, véhicules, etc.) ne fait pas l'objet actuellement d'un encadrement réglementaire.
L'utilisation d'eaux usées traitées pour des usages de nettoyage (voiries, ouvrages, véhicules etc.) interne à une STEU relève d'une maîtrise des risques par l'exploitant et des obligations de protection des travailleurs fixés par le code du travail (articles L. 4121-1 à 4), mais n'est pas, en l'état actuel de la réglementation, soumise à autorisation préfectorale préalable ;

Le ministère de la santé a saisi l'ANSES à ce sujet qui a répondu en 2012 :

« Lavage des voiries :

En l'absence de données d'exposition, notamment sur les caractéristiques des particules d'eaux émises par les engins de nettoyage des espaces publics, de données sur les scenarii d'exposition des travailleurs et des passants (seules catégories de population ayant été identifiées comme potentiellement exposées), il n'a pas été possible de réaliser une ERS pour ces deux catégories de population.

Cependant, considérant que les travailleurs, notamment les opérateurs de lance, pourraient être particulièrement exposés à des particules d'EUT, et dans l'optique de pouvoir mener une ERS, le CES recommande la réalisation d'une étude permettant de caractériser leurs expositions en fonction du matériel utilisé.

De même que pour l'irrigation, d'autres ressources en eau de qualité généralement non contrôlées ou contrôlées partiellement sont utilisées pour le lavage des voiries, le CES recommande la création d'une base de données afin de recenser et compiler les qualités d'eau utilisées. »

« Avez-vous connaissance d'autres débouchés possibles pour l'utilisation des eaux usées traitées hormis le rejet en milieu naturel et l'irrigation agricole/espaces verts ?

o Si oui, lesquels et selon quelles conditions ?

o Si non, peut-on s'attendre à une évolution de la législation pour introduire de nouveaux débouchés? »

« Des réflexions sont en cours au niveau national, portant sur des possibilités d'évolution de la réglementation (niveaux de qualité des eaux usées traitées A+ sans limitation d'usage, simplifications administratives...).

Le cadre réglementaire national pourrait être amené à évoluer afin de prendre en compte les éléments qui auront pu être validés aux niveaux européen et national.

Un texte devrait sortir sur la recharge de nappe. »

« Auriez-vous des éléments concernant l'état d'avancement des réflexions de la Commission Européenne autour de la REUT, notamment vis-à-vis des usages précédemment mentionnés? »

« Il n'existe pas à l'heure actuelle de cadre commun concernant la REUT au niveau européen. Cependant, la Commission européenne mène des réflexions sur cette thématique avec comme objectif d'aboutir à une proposition d'outil commun à l'ensemble des États membres concernant la REUT. Je n'ai aucune information sur l'état d'avancement de la réflexion

Une Consultation publique sur les possibilités d'introduire des exigences minimales de qualité concernant l'eau réutilisée dans l'Union européenne a été faite début 2017 »
Annexe 6 : Construction d'une carte d'occupation du sol

Définir l'occupation précise du sol est un travail nécessaire puisque cette dernière contient des informations sur les possibles usages de la REUT mais également sur les zones de contraintes existantes et obstacles à franchir. Bien que certains acteurs publics aient déjà effectués un certain nombre de travaux sur cette occupation des sols, ce n'est pas nécessairement le cas de l'intégralité des communes ou intercommunalités, et même si cette condition était remplie, rien ne garantit que les nomenclatures et précisions des données soient identiques et comparables. L'objectif ici est donc de proposer une méthode de définition de l'occupation du sol qui soit suffisamment précise pour être appliquée à l'échelle d'une commune mais potentiellement applicable à des ensembles plus vastes. Pour information l'IGN travaille actuellement sur la définition d'un référentiel d'occupation du sol à grande échelle (OCS GE®) qui devrait satisfaire les critères de la méthodologie. Mais tant que cette solution n'est pas utilisable il a été choisi de proposer une alternative.

Pour les objectifs du projet SoPoLo il a été choisi d'utiliser 23 classes d'occupation du sol au regard de l'utilisation future de cette information (calcul de chemin favorables, localisation des parcelles, etc.), la nomenclature est définie dans le tableau suivant :

<table><thead><tr><td>Code</td><td>Libellé</td><td>Source</td></tr></thead><tbody><tr><td>1</td><td>Inconnu</td><td>-</td></tr><tr><td>2</td><td>Cultures</td><td>RPG, CCGPSL, 3M</td></tr><tr><td>3</td><td>Vigne</td><td>RPG, BD Topo®, CCGPSL, 3M</td></tr><tr><td>4</td><td>Chemins et sentier</td><td>BD Topo®</td></tr><tr><td>5</td><td>Route empierrée</td><td>BD Topo®</td></tr><tr><td>6</td><td>Route goudronnée petite</td><td>BD Topo®</td></tr><tr><td>7</td><td>Route goudronnée grande</td><td>BD Topo®, 3M (A9)</td></tr><tr><td>8</td><td>Chemin de fer</td><td>BD Topo®, CCGPSL, 3M</td></tr><tr><td>9</td><td>Landes</td><td>BD Forêt®, RPG, CCGPSL, 3M</td></tr><tr><td>10</td><td>Foret ouverte</td><td>BD Forêt®</td></tr><tr><td>11</td><td>Peupleraie</td><td>BD Forêt®, CCGPSL, 3M</td></tr><tr><td>12</td><td>Taillis et futaies</td><td>BD Forêt®</td></tr><tr><td>13</td><td>Autres boisements</td><td>BD Topo®, CCGPSL, 3M</td></tr><tr><td>14</td><td>Zone urbaine</td><td>BD Topo®, CCGPSL, 3M</td></tr><tr><td>15</td><td>Décharges et carrières</td><td>CCGPSL, 3M</td></tr><tr><td>16</td><td>Affleurements</td><td>CCGPSL, 3M</td></tr><tr><td>17</td><td>Masse d'eau</td><td>BD Topo®, CCGPSL, 3M</td></tr><tr><td>18</td><td>Prairie</td><td>RPG, CCGPSL, 3M</td></tr><tr><td>19</td><td>Espaces verts urbains</td><td>CCGPSL, 3M</td></tr><tr><td>20</td><td>Sables</td><td>CCGPSL, 3M</td></tr><tr><td>21</td><td>Très petit cours (0 à 15m)</td><td>BD Carthage®</td></tr><tr><td>22</td><td>Petit cours (15 à 50m)</td><td>BD Carthage®</td></tr><tr><td>23</td><td>Cours moyen (Plus de 50m)</td><td>BD Carthage®</td></tr></tbody></table>

A noter qu'il existe des classes d'occupation du sol qui ne sont pas représentées (ex : tourbes...) et qui ont peut-être leur importance sur d'autres territoires. Le principe est d'obéir à une logique d'impédance pour le choix des classes : Quelles sont les type de couverture les plus représentatifs du territoire traduisant la facilité de déplacement d'un point A à un point B? Néanmoins certaines catégories d'occupation bien particulières ne peuvent être omises comme les sables ou les carrières.

Comme nous pouvons le voir, 5 sources de données principales sont utilisées pour construire cette
couche d'occupation du sol : BD Topo, BD Forêt, RPG, ainsi que des données d'occupation du sol et de plans d'urbanisme transmises par Montpellier Méditerranée Métropole (3M) et la Communauté de Communes du Grand-Pic-Saint-Loup. Concernant ces dernières sources, elles sont utilisées pour compléter l'occupation du sol dans le cas où celle-ci est inconnue.

En effet la BD Topo, la BD Forêt et le RPG, bien que couvrant la majeure partie de la zone d'étude (80%), ne conviennent pas aux besoins de la méthodologie, la couverture devant être intégrale pour pouvoir notamment étudier les chemins de moindre contrainte. Les zones non couvertes (inconnu) sont donc complétées avec les informations transmises par les intercommunalités.

Les modalités de croisement seront présentées dans la suite ainsi que des alternatives.

Parmi les données présentées plus haut, l'une d'entre-elles a fait l'objet d'un travail en particulier : la zone urbaine.

Il existe en effet plusieurs sources d'information sur le contour des espaces urbaine : Corine Land Cover, la BD CARTO, les documents d'urbanismes ou encore la BD Topo. Malheureusement les trois premières sources ont tendance à être soit trop grossières, soit à omettre les bâtis en zone rurales. A contrario la BD Topo fournit le contour des bâtiments, mais ne délimite pas d'aire urbaine à proprement parler. Nous avons donc besoin d'une information intermédiaire pour pouvoir répondre aux besoins de la méthodologie.

> L'INSEE définit une unité urbaine comme une « commune ou un ensemble de communes qui comporte sur son territoire une zone bâtie d'au moins 2 000 habitants où aucune habitation n'est séparée de la plus proche de plus de 200 mètres. En outre, chaque commune concernée possède plus de la moitié de sa population dans cette zone bâtie. »

Ce seuil de 200 mètres est intéressant car il permet de définir une distance d'agrégation du bâti, permettant de passer d'un ensemble de bâtiment à une aire urbaine proprement dite.

Deux différentes méthodes ont été testées en utilisant ce seuil : Une méthode par dilatation-érosion et une méthode par agrégation de polygone.

La dilatation-érosion est une méthode classiquement utilisée pour agréger des entités, le principe est de réaliser une zone tampon autour de chaque polygone (ici de 200m), de fusionner ces zones tampon puis de les éroder de la même distance (200m). Si la distance entre les entités de base est suffisamment faible, alors les entités seront reliées (voir captures d'écran suivantes).

<figure><img src="image_48.png" /><figcaption>Dilatation/érosion, en noir le bâti, en rouge la zone tampon de 200m et en vert la zone après érosion de 200m</figcaption></figure>
Le défaut de cette méthode est dans la morphologie des aires résultantes qui sont arrondies et ne correspondent pas à la représentation d'une aire urbaine (qui devrait être plutôt polygonale).

Pour pallier à ce défaut la méthode d'agrégation de polygones (proposée dans ArcGIS, outil *Agrégation de polygones*) permet de relier par une droite les entités séparées par une distance inférieure au seuil défini (200m ici). De plus il est possible de définir des entités d'interruption qui empêchent de relier deux bâtiments si ceux-ci sont séparés par une entité d'interruption, typiquement il peut s'agir d'un fleuve, d'un champ, d'une falaise, etc. Donnant plus de sens au découpage de l'aire urbaine :

<figure><img src="image_49.png" /><figcaption>Aires urbaines par agrégation de 200m (gris), en vert clair le résultat par dilatation-érosion, en noir le bâti, et en vert foncé le boisement (entité d'interruption)</figcaption></figure>

Néanmoins ce seuil de 200m peut poser problème en zone peu densément bâtie comme l'illustre la situation suivante :

<figure><img src="image_50.png" /><figcaption>Agrégation à 200m et à 100m</figcaption></figure>

Sur la photo ci-dessus le résultat de l'agrégation à 200m (en gris) traverse clairement des zones boisées et de cultures, un problème nettement moins visible si le seuil est abaissé à 100m (en vert).

Par soucis de cohérence, il a donc été choisi de conserver cette valeur de 100m. Il s'agit ici typiquement d'une valeur qui peut être choisie manuellement suivant le contexte (grande ville, espace rural, etc.).
Toutefois il subsiste deux défauts au croisement effectué : En ayant choisi les boisements comme entité d'interruption, il existe un certain nombre de trous dans les tâches urbaines qui correspondent en réalisé à des bois, espaces verts inclus dans les villes mais qui doivent être traités comme faisant partie de l'ensemble urbain. Le second point est celui de la portée temporelle de ces zones urbaines. En effet un projet de REUT se planifie sur le long terme, afin de pouvoir prendre en compte ce point, notamment dans l'étude des chemins, il peut être important d'écarter les zones à urbaniser des plans d'urbanisme.

Bien qu'il soit impossible de ne pas prendre en compte les boisements en entité d'interruptions (apparition d'aires urbaines dans des blocs forestiers en milieu rural) les deux défauts précédemment cités se règlent si l'on superpose le résultat de l'agrégation à 100m aux aires urbaines et à urbaniser des PLU, dans lesquelles sont incluses les espaces verts.

<figure><img src="image_51.png" /><figcaption>Définition des aires urbaines (gris)</figcaption></figure>

Le résultat final permet de localiser des grands ensembles urbains actuels et avenirs sans négliger le bâti peu dense et isolé. Le détail des opérations est illustré en annexe 2, diagramme 1.

Une fois les aires urbaines délimitées, il convient de délimiter l'emprise des voies de communication terrestres. Différent type de voies se distinguent : chemins et sentiers, routes empierrées, et les routes « classiques » (goudronnées) suivant la facilité avec laquelle elles sont franchissables pour des travaux. Il a été décidé qu'une route ayant une emprise supérieure à 10 mètres était infranchissable, on distingue donc les routes « classiques » moyennes et grandes.

L'emprise elle a été calculée en utilisant les travaux du Cerema (Centre d'études et d'expertise sur les risques, l'environnement, la mobilité et l'aménagement) sur la mesure de la consommation d'espace<fnref n="12" />.

Les aires urbaines dessinées et l'emprise routière calculée, il convient de construire la couche d'occupation du sol. Cette dernière est construire par opération de *mise à jour* successive (opération permettant d'insérer un polygone dans une couche effaçant les entités ainsi recouvertes).

<footnote n="12">¹²http://piece-jointe-carto.developpement-durable.gouv.fr/REG074B/FONCIER SOL/N OCCUPATION SOL/L EMPRISE ROUTE R74/Fiche1-7-1.pdf</footnote>
<figure><img src="image_52.png" /></figure>

Résultat de l'opération de mise à jour

La mise à jour permet de partir d'un ensemble de couches pour arriver à une couche unique par « superposition » (notion d'ordre sous entendue).

L'ordre de mise à jour est le suivant : Cultures, vignes, forêt (BD Topo), forêt (BD Forêt), zones urbaines, voie ferrées, routes, chemins, masses d'eau. Un champ « catégorie » est préalablement créé dans chacune des couches et prends la valeur du code indiquée dans la nomenclature. Ce champ est conservé et mis à jour lui aussi pendant l'opération. Le détail de cette opération figure en annexe 2, diagramme 3.

On réalise ensuite une opération de fusion (association de polygones adjacents qui possèdent une ou plusieurs propriété identique, ici l'occupation du sol) et on élimine les polygones résiduels dont la surface est inférieure à 25m² (taille de la maille du MNT et donc de la carte de coût, l'élimination fusionne le polygone avec le polygone adjacent ayant la surface de contact la plus importante).

Une fois cette opération effectuée, on extrait les polygones dont l'occupation est inconnue. Ceux-ci sont croisés (Intersection) avec les couches d'occupation du sol de 3M et CCGPSL). Les deux couches résultantes sont ensuite réinjecté dans la couche d'occupation du sol par mise à jour. On ré-effectue le nettoyage de couche (fusion/élimination). La couche résultante ne possède ainsi plus aucune zone inconnue.

Les forêts issues de la BD Topo permettent de localiser des surfaces boisées plus petites que celles de la BD Forêt, en revanche cette dernière est plus intéressante sur le plan sémantique (différentiation de landes, forêt ouvertes, peupleraies, taillis) pour l'utilisation que l'on souhaite en faire vis-à-vis de la carte de coût. Elle est donc superposée à la BD Topo.

Le résultat final est ensuite converti en raster d'un pas de 5m. Le diagramme 4 en annexe 4 synthétise le processus. La carte 2 en annexe 9 présente le résultat.
# Annexe 7 : Précisions techniques sur le bilan hydrique

La présente annexe a pour but de donner des informations complémentaires au calcul de déficit hydrique présenté en section **Erreur ! Source du renvoi introuvable**. Il s'agit notamment de présenter la manière dont les résultats ont été extraits de la base de données BDSol L-R, des tables Safran et calculées sous Excel.

## A. Extraction de données pédologiques :

La base de données BDSol L-R est construite selon le modèle DoneSol (version 3) commun à l'ensemble des RRP. Notamment les données surfaciques, organisées en strates, UTS et UCS (cf. section 2.5.1.2) sont reliées de la façon suivante :

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

Relations entre les tables de données surfaciques de la BDSol L-R, d'après le dictionnaire de données DoneSol Version 3

Au cours du calcul de déficit hydrique, certains paramètres du sol sont recherchés :

- Pour les strates : taux de sables, d'argile et de limon, pourcentage en éléments grossier, leur nature et épaisseur de la strate

- Pour les UTS, les strates les composant

- Pour les UCS, les UTS les composant

Pour obtenir ces différentes informations est utilisé le système de gestion de base de données Access dans lequel on réalise des requêtes en langage SQL pour extraire les informations recherchées. L'exemple suivant donne la requête permettant d'extraire le pourcentage d'UTS dans chaque UCS ainsi que des informations de nom et d'identifiant :

```sql
SELECT info_ucs.pourcent, uts.id_uts, uts.rp_95_nom, info_ucs.nom_ucs, info_ucs.no_int_ucs, info_ucs.no_ucs
```

```sql
FROM (l_ucs_uts INNER JOIN uts ON l_ucs_uts.id_uts = uts.id_uts) INNER JOIN info_ucs ON
(l_ucs_uts.pourcent = info_ucs.pourcent) AND (l_ucs_uts.id_ucs = info_ucs.id_ucs);
```
La capture d'écran suivante présente le résultat de cette requête :

<table><thead><tr><th>pourcent</th><th>no_int_ucs</th><th>id_uts</th><th>no_ucs</th><th>nom_ucs</th><th>rp_95_nom</th></tr></thead><tbody><tr><td>2</td><td>505F</td><td>9117</td><td>50506</td><td>Pentes faibles et dépressions des t</td><td>CALCOSOL</td></tr><tr><td>5</td><td>551A</td><td>9117</td><td>55101</td><td>Petits domes et cuvettes en pentes</td><td>CALCOSOL</td></tr><tr><td>5</td><td>551A</td><td>9117</td><td>55101</td><td>Petits domes et cuvettes en pentes</td><td>CALCOSOL</td></tr><tr><td>3</td><td>551B</td><td>9117</td><td>55102</td><td>Dépressions cultivées associées à c</td><td>CALCOSOL</td></tr><tr><td>3</td><td>551B</td><td>9117</td><td>55102</td><td>Dépressions cultivées associées à c</td><td>CALCOSOL</td></tr><tr><td>10</td><td>519A</td><td>9119</td><td>51901</td><td>Plateaux de calcaire lacustre massif</td><td>RENDOSOL</td></tr><tr><td>30</td><td>309W</td><td>9136</td><td>30923</td><td>Alluvions villafranchiennes d'origi</td><td>FERSIALSOL luvique rédoxique</td></tr><tr><td>3</td><td>505F</td><td>9138</td><td>50506</td><td>Pentes faibles et dépressions des t</td><td>CALCOSOL rédoxique</td></tr><tr><td>3</td><td>551B</td><td>9138</td><td>55102</td><td>Dépressions cultivées associées à c</td><td>CALCOSOL rédoxique</td></tr><tr><td>3</td><td>551B</td><td>9138</td><td>55102</td><td>Dépressions cultivées associées à c</td><td>CALCOSOL rédoxique</td></tr><tr><td>2</td><td>557A</td><td>9145</td><td>55701</td><td>Cuvettes et bordures peu déclives.</td><td>LUVIOSOL hydromorphe</td></tr><tr><td>1</td><td>551B</td><td>9157</td><td>55102</td><td>Dépressions cultivées associées à c</td><td>CALCOSOL</td></tr><tr><td>2</td><td>309W</td><td>9208</td><td>30923</td><td>Alluvions villafranchiennes d'origi</td><td>BRUNISOL</td></tr><tr><td>100</td><td>104E</td><td>9221</td><td>10405</td><td>Plaines alluviales des petites rivière</td><td>FLUVIOSOL</td></tr><tr><td>60</td><td>106B</td><td>9223</td><td>10602</td><td>Vistrenque, Uzégeois (Gard), Lunel</td><td>FLUVIOSOL calcaire à horizons réc</td></tr><tr><td>100</td><td>171B</td><td>9224</td><td>17102</td><td>Corbières. Bordures cévenoles. Val</td><td>FLUVIOSOL calcaire</td></tr><tr><td>40</td><td>106B</td><td>9225</td><td>10602</td><td>Vistrenque, Uzégeois (Gard), Lunel</td><td>FLUVIOSOL calcaire à horizon pro</td></tr><tr><td>100</td><td>106A</td><td>9233</td><td>10601</td><td>Couloir rhodanien (Gard) et vallée</td><td>COLLUVIOSOL à horizon réductique</td></tr><tr><td>20</td><td>173B</td><td>9236</td><td>17302</td><td>Contacts plaines fluviatiles et littor</td><td>FLUVIOSOL calcaire rédoxique</td></tr><tr><td>40</td><td>173B</td><td>9237</td><td>17302</td><td>Contacts plaines fluviatiles et littor</td><td>FLUVIOSOL à horizons réductique</td></tr><tr><td>40</td><td>173B</td><td>9237</td><td>17302</td><td>Contacts plaines fluviatiles et littor</td><td>FLUVIOSOL à horizons réductique</td></tr></tbody></table>

Résultats d'une requête SQL sous Access : Obtention des proportions relatives de chaque UTS par UCS

Ces données sont ensuite exportées sous Excel pour traitement. L'utilisation de tableaux croisés dynamiques est vivement conseillée pour agréger les informations extraites à l'échelle d'UTS ou d'UCS.

Le calcul de réserve utile est ensuite réalisé selon les modalités présentées en section 2.5.1.2.

## B. Extraction des données SAFRAN

Les données SAFRAN sont fournies sous deux formats, un fichier forme (Shapefile) localisant les centres des grilles SAFRAN. Attention, le système de coordonnées utilisé pour les données SAFRAN est le Lambert II étendu. L'information est ici re-projetée en Lambert 93 (outil *Projeter*). Enfin les points sont utilisés pour calculer les grilles (polygones) avec l'outil *Créer un quadrillage* :

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

Ces grilles peuvent ensuite être intersectées avec le parcellaire agricole pour pouvoir associer les données climatiques aux parcelles.

Ces données climatiques constituent le second format de données SAFRAN. Elles sont stockées sous
format .csv dans le cadre de cette étude. Un fichier .csv est propre à chaque grille et indique l'ensemble des paramètres du climat (ETP, températures, pluviométries, etc.) de 1958 à 2016 avec un pas de temps journalier. Un script VBA a été écrit pour importer l'ensemble des grilles en lien avec le territoire d'étude :

Sub ImportPLEtextFiles()
Dim i As Integer "initial file name"
Dim k As Integer "final file name"
Dim DQ As String " Double quotation marks"
Dim Folder As String
Dim File As String
k = 2
DQ = ""
Folder = Worksheets("Feuil1").Cells(3, 3).Value
Do While k <= 32
i = Worksheets("Feuil1").Cells(k, 1).Value
Worksheets.Add.Name = i
File = "TEXT;" & Folder & i & "_SAFRAN.txt"
Worksheets("Feuil1").Cells(6, 6).Value = File
With ActiveSheet.queryTables.Add(
    Connection:=File,
    Destination:=Range("$A$1"))
    .FieldNames = True
    .RowNumbers = False
    .FillAdjacentFormulas = False
    .PreserveFormatting = True
    .RefreshOnFileOpen = False
    .RefreshStyle = xlInsertDeleteCells
    .SavePassword = False
    .SaveData = True
    .AdjustColumnWidth = False
    .RefreshPeriod = 0
    .TextFilePromptOnRefresh = False
    .TextFilePlatform = 850
    .TextFileStartRow = 31
    .TextFileParseType = xlDelimited
    .TextFileTextQualifier = xlTextQualifierDoubleQuote
    .TextFileConsecutiveDelimiter = False
    .TextFileTabDelimiter = False
    .TextFileSemicolonDelimiter = True
    .TextFileCommaDelimiter = False
    .TextFileSpaceDelimiter = False
    .TextFileColumnDataTypes = Array(1)
    .TextFileTrailingMinusNumbers = True
    .Refresh BackgroundQuery:=False
End With
k = k + 1
Loop
End Sub

Script VBA permettant l'importation de l'ensemble des informations des grilles SAFRAN du territoire sous Excel
Le script est basé sur une boucle (while) lisant le nom du fichier à importer (File = "TEXT;" & Folder & i & "SAFRAN.txt") et qui extrait les valeurs dans une feuille correspondant à la grille SAFRAN importée (ActiveSheet.QueryTables.Add). Les 30 premières lignes étant des données descriptives ne sont pas importées (TextFileStartRow = 31).

Les données antérieures à 1986 sont supprimées, les informations d'ETP et de pluviométrie sont moyennées pour chaque jour sur la période 1986-2016 à l'aide d'un tableau croisé dynamique. Le résultat est synthétisés sous la forme de deux tableaux à double entrée chacun présentant la valeur d'un paramètre (ETP et pluviométrie respectivement) en fonction du jour et de la grille. Les données sont stockées sous le format suivant (exemple pour la pluie, les colonnes correspondent aux grilles):

<table><thead><tr><th>jour_mois</th><th>n°</th><th>8497</th><th>8498</th><th>8586</th><th>8587</th></tr></thead><tbody><tr><td>1_1</td><td>1</td><td>1,9</td><td>2,5</td><td>4,0</td><td>2,1</td></tr><tr><td>2_1</td><td>2</td><td>1,8</td><td>2,1</td><td>2,2</td><td>2,8</td></tr><tr><td>3_1</td><td>3</td><td>2,8</td><td>2,3</td><td>4,6</td><td>6,1</td></tr><tr><td>4_1</td><td>4</td><td>1,3</td><td>1,3</td><td>2,6</td><td>2,3</td></tr><tr><td>5_1</td><td>5</td><td>2,2</td><td>2,3</td><td>2,2</td><td>1,3</td></tr><tr><td>6_1</td><td>6</td><td>3,3</td><td>3,1</td><td>4,8</td><td>2,7</td></tr><tr><td>7_1</td><td>7</td><td>0,5</td><td>0,6</td><td>1,8</td><td>0,6</td></tr><tr><td>8_1</td><td>8</td><td>3,2</td><td>2,6</td><td>3,9</td><td>2,8</td></tr><tr><td>9_1</td><td>9</td><td>1,7</td><td>1,4</td><td>3,4</td><td>3,3</td></tr><tr><td>10_1</td><td>10</td><td>2,9</td><td>3,9</td><td>4,5</td><td>2,1</td></tr><tr><td>11_1</td><td>11</td><td>1,3</td><td>0,9</td><td>2,6</td><td>1,8</td></tr><tr><td>12_1</td><td>12</td><td>3,9</td><td>3,7</td><td>4,2</td><td>2,9</td></tr><tr><td>13_1</td><td>13</td><td>5,1</td><td>5,6</td><td>7,5</td><td>4,7</td></tr><tr><td>14_1</td><td>14</td><td>6,0</td><td>4,7</td><td>7,4</td><td>5,6</td></tr></tbody></table>

Tableau synthétique à double entrée présentant les valeurs moyennes de précipitation en fonction du jour et des grilles SAFRAN.
C. Données de coefficients culturaux

Les Kc utilisés dans cette étude sont issus pour la majorité des cas du Mémento BRL Irrigation, le Kc utilisé pour l'arboriculture est issu des données FAO pour le citronnier tandis que les valeurs de Kc pour la vigne sont issues de la bibliographie (Rousseau, 2009) :

<table><thead><tr><th></th><th>no_mois</th><th>1</th><th>1</th><th>1</th><th>2</th><th>2</th><th>2</th><th>3</th><th>3</th><th>3</th><th>4</th><th>4</th><th>4</th></tr><tr><th></th><th>no_decade</th><th>1</th><th>2</th><th>3</th><th>1</th><th>2</th><th>3</th><th>1</th><th>2</th><th>3</th><th>1</th><th>2</th><th>3</th></tr></thead><tbody><tr><td>BLE TENDRE</td><td>1</td><td>0,5</td><td>0,5</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1,2</td><td>1,2</td><td>1,2</td><td>1,2</td></tr><tr><td>MAIS GRAIN ET ENSILAGE</td><td>2</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>ORGE</td><td>3</td><td>0,675</td><td>0,675</td><td>0,675</td><td>0,825</td><td>0,825</td><td>0,825</td><td>1</td><td>1</td><td>1</td><td>1,125</td><td>1,125</td><td>1,125</td></tr><tr><td>AUTRES CEREALES</td><td>4</td><td>0,5</td><td>0,5</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1,2</td><td>1,2</td><td>1,2</td><td>1,2</td></tr><tr><td>COLZA</td><td>5</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td></tr><tr><td>TOURNESOL</td><td>6</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0,5</td></tr><tr><td>PROTEAGINEUX</td><td>8</td><td>0,5</td><td>0,5</td><td>0,5</td><td>0,7</td><td>0,7</td><td>0,7</td><td>1</td><td>1</td><td>1</td><td>1,1</td><td>1,1</td><td>1,1</td></tr><tr><td>SEMENCES</td><td>10</td><td>0,5</td><td>0,5</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1,2</td><td>1,2</td><td>1,2</td><td>1,2</td></tr><tr><td>AUTRES GELS</td><td>13</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>LEGUMINEUSES A GRAINS</td><td>15</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0,4</td><td>0,4</td><td>0,025</td><td>0,025</td><td>0,025</td><td>0,025</td></tr><tr><td>FOURRAGE</td><td>16</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td></tr><tr><td>ESTIVES LANDES</td><td>17</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>PRAIRIES PERMANENTES</td><td>18</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td><td>1,1</td><td>1,1</td><td>1,1</td><td>1,1</td><td>1,1</td><td>1,1</td></tr><tr><td>PRAIRIES TEMPORAIRES</td><td>19</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td><td>1,1</td><td>1,1</td><td>1,1</td><td>1,1</td><td>1,1</td><td>1,1</td></tr><tr><td>VERGERS</td><td>20</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0,5</td><td>0,5</td><td>0,5</td></tr><tr><td>VIGNES</td><td>21</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0,2</td><td>0,2</td><td>0,2</td></tr><tr><td>FRUITS A COQUE</td><td>22</td><td>0,5</td><td>0,5</td><td>0,5</td><td>0,5</td><td>0,5</td><td>0,5</td><td>0,055</td><td>0,055</td><td>0,055</td><td>0,055</td><td>0,055</td><td>0,055</td></tr><tr><td>OLIVIERS</td><td>23</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0,55</td><td>0,55</td><td>0,55</td></tr><tr><td>LEGUMES-FLEURS</td><td>25</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0,6</td><td>0,6</td></tr><tr><td>ARBORICULTURE</td><td>27</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td></tr><tr><td>DIVERS</td><td>28</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr></tbody></table>
<table><thead><tr><th></th><th>no_mois</th><th>5</th><th>5</th><th>5</th><th>6</th><th>6</th><th>6</th><th>7</th><th>7</th><th>7</th><th>8</th><th>8</th><th>8</th></tr><tr><th></th><th>no_decade</th><th>1</th><th>2</th><th>3</th><th>1</th><th>2</th><th>3</th><th>1</th><th>2</th><th>3</th><th>1</th><th>2</th><th>3</th></tr></thead><tbody><tr><td>BLE TENDRE</td><td>1</td><td>1,2</td><td>1,2</td><td>1,1</td><td>1,1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>MAIS GRAIN ET ENSILAGE</td><td>2</td><td>0,4</td><td>0,4</td><td>0,4</td><td>0,8</td><td>0,8</td><td>1,1</td><td>1,1</td><td>1,1</td><td>1,1</td><td>1,1</td><td>1</td><td>1</td></tr><tr><td>ORGE</td><td>3</td><td>1,15</td><td>1,15</td><td>1,15</td><td>0,475</td><td>0,475</td><td>0,475</td><td>0,25</td><td>0,25</td><td>0,25</td><td>0</td><td>0</td><td>0</td></tr><tr><td>AUTRES CEREALES</td><td>4</td><td>1,2</td><td>1,2</td><td>1,1</td><td>1,1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>COLZA</td><td>5</td><td>1</td><td>1</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>TOURNESOL</td><td>6</td><td>0,5</td><td>0,8</td><td>0,8</td><td>1,1</td><td>1,1</td><td>1,1</td><td>1,1</td><td>1,1</td><td>0,6</td><td>0,6</td><td>0</td><td>0</td></tr><tr><td>PROTEAGINEUX</td><td>8</td><td>1,1</td><td>1,1</td><td>1,1</td><td>1,1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>SEMENCES</td><td>10</td><td>1,2</td><td>1,2</td><td>1,1</td><td>1,1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>AUTRES GELS</td><td>13</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>LEGUMINEUSES A GRAINS</td><td>15</td><td>1,15</td><td>1,15</td><td>1,15</td><td>0,35</td><td>0,35</td><td>0,35</td><td>0,35</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>FOURRAGE</td><td>16</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td></tr><tr><td>ESTIVES LANDES</td><td>17</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>PRAIRIES PERMANENTES</td><td>18</td><td>1,1</td><td>1,1</td><td>1,1</td><td>1,1</td><td>1,1</td><td>1,1</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td></tr><tr><td>PRAIRIES TEMPORAIRES</td><td>19</td><td>1,1</td><td>1,1</td><td>1,1</td><td>1,1</td><td>1,1</td><td>1,1</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td></tr><tr><td>VERGERS</td><td>20</td><td>0,85</td><td>0,85</td><td>0,85</td><td>0,85</td><td>0,85</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>0,7</td><td>0,7</td></tr><tr><td>VIGNES</td><td>21</td><td>0,3</td><td>0,3</td><td>0,3</td><td>0,4</td><td>0,4</td><td>0,4</td><td>0,4</td><td>0,4</td><td>0,4</td><td>0,5</td><td>0,5</td><td>0,5</td></tr><tr><td>FRUITS A COQUE</td><td>22</td><td>1,05</td><td>1,05</td><td>1,05</td><td>1,05</td><td>1,05</td><td>1,05</td><td>1,05</td><td>1,05</td><td>1,05</td><td>1,05</td><td>1,05</td><td>1,05</td></tr><tr><td>OLIVIERS</td><td>23</td><td>0,55</td><td>0,55</td><td>0,55</td><td>0,55</td><td>0,55</td><td>0,55</td><td>0,55</td><td>0,55</td><td>0,55</td><td>0,55</td><td>0,55</td><td>0,55</td></tr><tr><td>LEGUMES-FLEURS</td><td>25</td><td>0,6</td><td>0,6</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,5</td><td>0,5</td><td>0,5</td><td>0</td><td>0</td><td>0</td></tr><tr><td>ARBORICULTURE</td><td>27</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td></tr><tr><td>DIVERS</td><td>28</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr></tbody></table>
<table><thead><tr><th></th><th>no_mois</th><th>9</th><th>9</th><th>9</th><th>10</th><th>10</th><th>10</th><th>11</th><th>11</th><th>11</th><th>12</th><th>12</th><th>12</th></tr><tr><th></th><th>no_decade</th><th>1</th><th>2</th><th>3</th><th>1</th><th>2</th><th>3</th><th>1</th><th>2</th><th>3</th><th>1</th><th>2</th><th>3</th></tr></thead><tbody><tr><td>BLE TENDRE</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0,5</td><td>0,5</td><td>0</td></tr><tr><td>MAIS GRAIN ET ENSILAGE</td><td>2</td><td>1</td><td>0,7</td><td>0,7</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>ORGE</td><td>3</td><td>0</td><td>0</td><td>0</td><td>0,3</td><td>0,3</td><td>0,3</td><td>0,35</td><td>0,35</td><td>0,35</td><td>0,5</td><td>0,5</td><td>0,5</td></tr><tr><td>AUTRES CEREALES</td><td>4</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0,5</td><td>0,5</td><td>0</td></tr><tr><td>COLZA</td><td>5</td><td>0</td><td>0</td><td>0,3</td><td>0,3</td><td>0,3</td><td>0,7</td><td>0,7</td><td>0,7</td><td>1</td><td>1</td><td>1</td><td>1</td></tr><tr><td>TOURNESOL</td><td>6</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>PROTEAGINEUX</td><td>8</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0,5</td><td>0,5</td></tr><tr><td>SEMENCES</td><td>10</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0,5</td><td>0,5</td><td>0</td></tr><tr><td>AUTRES GELS</td><td>13</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>LEGUMINEUSES A GRAINS</td><td>15</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>FOURRAGE</td><td>16</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td></tr><tr><td>ESTIVES LANDES</td><td>17</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>PRAIRIES PERMANENTES</td><td>18</td><td>1,1</td><td>1,1</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td></tr><tr><td>PRAIRIES TEMPORAIRES</td><td>19</td><td>1,1</td><td>1,1</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td><td>0,6</td></tr><tr><td>VERGERS</td><td>20</td><td>0,7</td><td>0,7</td><td>0,7</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>VIGNES</td><td>21</td><td>0,5</td><td>0,5</td><td>0,5</td><td>0,25</td><td>0,25</td><td>0,25</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>FRUITS A COQUE</td><td>22</td><td>1,05</td><td>1,05</td><td>1,05</td><td>0,65</td><td>0,65</td><td>0,65</td><td>0,5</td><td>0,5</td><td>0,5</td><td>0,5</td><td>0,5</td><td>0,5</td></tr><tr><td>OLIVIERS</td><td>23</td><td>0,55</td><td>0,55</td><td>0,55</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>LEGUMES-FLEURS</td><td>25</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>ARBORICULTURE</td><td>27</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td><td>0,8</td></tr><tr><td>DIVERS</td><td>28</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr></tbody></table>
## D. Calcul du déficit hydrique :

Le calcul du déficit hydrique est ensuite réalisé en croisant les facteurs précédemment rassemblés (RU, ETP, pluviométrie, Kc).

Sous Excel, il est réalisé en trois temps. Tout d'abord la réserve réelle est calculée, elle correspond au calcul suivant :

$$ Réserve réelle_{t+1} = Réserve relative_{t} + Précipitation_{t+1} - Consommation_{t+1} $$

Elle peut prendre des valeurs négatives mais ne peut dépasser la valeur de réserve utile maximale du sol. A t=0 la réserve est supposée maximale.

La réserve relative vérifie les conditions suivantes :

$$ Si Réserve réelle_{t} \geq 0 alors Réserve relative_{t} = Réserve réelle_{t} $$

$$ Sinon Réserve relative_{t} = 0 $$

Enfin le déficit hydrique vérifie les conditions suivantes :

$$ Si Réserve réelle_{t} \geq 0 alors Déficit_{t} = 0 $$

$$ Sinon Déficit_{t} = 0 $$

Pour faciliter le calcul, les fonctions INDEX et EQUIV ont été combinées pour rechercher les paramètres souhaités (ETP, Kc, pluviométrie) à chaque date, pour chaque grille et chaque UCS. Le processus de calcul du déficit hydrique est schématisé en fin d'annexe.

## E. Données complémentaires :

Teneur en eau des sols selon les classes texturales (Jamagne et al, 1977)

<table><thead><tr><th>Classe de texture</th><th>Humidité % à la capacité au champ (pF=2.5)</th><th>Humidité % au point de flétrissement permanent (pF=4.2)</th><th>Eau utile %</th><th>Densité apparente</th><th>Réserve Utile (mm/cm) RU</th></tr></thead><tbody><tr><td>S</td><td>8</td><td>3</td><td>5</td><td>1.35</td><td>0.7</td></tr><tr><td>SL</td><td>12</td><td>5</td><td>7</td><td>1.40</td><td>1.0</td></tr><tr><td>SA</td><td>19</td><td>10</td><td>9</td><td>1.50</td><td>1.35</td></tr><tr><td>LIS</td><td>15</td><td>7</td><td>8</td><td>1.50</td><td>1.20</td></tr><tr><td>LS</td><td>19</td><td>9</td><td>10</td><td>1.45</td><td>1.45</td></tr><tr><td>LmS</td><td>20</td><td>9</td><td>11</td><td>1.45</td><td>1.60</td></tr><tr><td>LSA</td><td>22</td><td>11</td><td>11</td><td>1.50</td><td>1.65</td></tr><tr><td>LAS</td><td>24</td><td>12</td><td>12</td><td>1.45</td><td>1.75</td></tr><tr><td>LI</td><td>17</td><td>8</td><td>9</td><td>1.45</td><td>1.30</td></tr><tr><td>Lm</td><td>23</td><td>10</td><td>13</td><td>1.35</td><td>1.75</td></tr><tr><td>LA</td><td>27</td><td>13</td><td>14</td><td>1.40</td><td>1.95</td></tr><tr><td>AS</td><td>33</td><td>22</td><td>11</td><td>1.55</td><td>1.70</td></tr><tr><td>A</td><td>37</td><td>25</td><td>12</td><td>1.45</td><td>1.75</td></tr><tr><td>AL</td><td>32</td><td>19</td><td>13</td><td>1.40</td><td>1.80</td></tr><tr><td>Alo (sédimentaire)</td><td>29</td><td>18</td><td>11</td><td>1.50</td><td>1.65</td></tr><tr><td>Alo (altération)</td><td>38</td><td>25</td><td>13</td><td>1.30</td><td>1.70</td></tr></tbody></table>
<figure><img src="image_55.png" /><figcaption>Triangle textural de l'Aisne (Jamagne et al, 1977)</figcaption></figure>

## Surfaces (en ha) par classe de déficit hydrique.

<table><thead><tr><td>Déficit en m3/ha</td><td>Surfaces (ha)</td></tr></thead><tbody><tr><td>Plus de 3 000</td><td>195</td></tr><tr><td>Entre 2 000 et 3 000</td><td>950</td></tr><tr><td>Entre 1 000 et 2 000</td><td>1235</td></tr><tr><td>Entre 500 et 1000</td><td>2086</td></tr><tr><td>Moins de 500</td><td>3172</td></tr><tr><td>Nul</td><td>17967</td></tr></tbody></table>
<figure><img src="image_56.png" /></figure>
# Annexe 8 : Reclassement des critères pour l'intérêt zonal

Pour les calculs d'aptitude zonale, les classes suivantes ont été déterminées pour chacun des facteurs. Sont également mentionné les valeurs de reclassement utilisées pour l'AHP (pour rappel plus la valeur est haute plus la classe est « favorable »).

## Déficit hydrique :

<table><thead><tr><td>Classe de déficit (m³/ha)</td><td>Reclasse-ment</td></tr></thead><tbody><tr><td>Nul</td><td>0</td></tr><tr><td>Entre 0 et 500</td><td>2</td></tr><tr><td>Entre 500 et 1000</td><td>5</td></tr><tr><td>Entre 1000 et 2000</td><td>7</td></tr><tr><td>Entre 2000 et 3000</td><td>8</td></tr><tr><td>Plus de 3000</td><td>9</td></tr></tbody></table>

## Diversité d'usage :

<table><thead><tr><td>Nombre d'usages</td><td>Reclasse-ment</td></tr></thead><tbody><tr><td>1</td><td>1</td></tr><tr><td>2</td><td>2</td></tr><tr><td>3</td><td>5</td></tr><tr><td>4</td><td>8</td></tr><tr><td>5</td><td>9</td></tr><tr><td>6</td><td>9</td></tr></tbody></table>

## Nombre de contraintes :

<table><thead><tr><td>Niveau de contrainte</td><td>Reclasse-ment</td></tr></thead><tbody><tr><td>0</td><td>9</td></tr><tr><td>1</td><td>7</td></tr><tr><td>2</td><td>5</td></tr><tr><td>3</td><td>3</td></tr><tr><td>4</td><td>0</td></tr></tbody></table>

## Diversité d'agriculteurs :

<table><thead><tr><td>Nombre d'agriculteurs</td><td>Reclasse-ment</td></tr></thead><tbody><tr><td>1</td><td>1</td></tr><tr><td>2</td><td>2</td></tr><tr><td>3</td><td>5</td></tr><tr><td>4</td><td>7</td></tr><tr><td>5</td><td>9</td></tr><tr><td>6</td><td>9</td></tr><tr><td>7</td><td>9</td></tr><tr><td>8</td><td>9</td></tr></tbody></table>
# Annexe 9 : Cartes

Carte 1 : Comparaison des masses d'eau selon la source de donnée considérée

Carte 2 : Occupation du sol sur le territoire d'étude

Carte 3 : Chemin de coût entre station et parcelle d'intérêt

Carte 4 : Chemin de coût entre station et meilleures parcelles
Comparaison des masses d'eau selon la source de donnée considérée

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

Réalisation : F. Hubert, IRSTEA Montpellier, UMR G-EAU, 2017.

Sources : IGN, WISE, OpenStreetMap
Occupation du sol sur le territoire d'étude

<figure><img src="image_58.png" /></figure>
Chemin de moindre coût
entre station et meilleures parcelles

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

Réalisation : F. Hubert, IRSTEA Montpellier, UMR G-EAU, 2017

Sources : BDSol L-R, IGN, Ministère de la Transition

Ecologique et Solidaire, SAFRAN, ASP, OpenStreetMap.

Fond : BD Ortho®

<figure><img src="image_60.png" /></figure>
Chemin de moindre coût entre station et parcelles d'intérêt

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

Réalisation : F. Hubert; IRSTEA Montpellier, UMR G-EAU, 2017

Sources : BDSol L-R, IGN, Ministère de la Transition

Ecologique et Solidaire, SAFRAN, ASP OpenStreetMap.

Fond : BD Ortho®

<figure><img src="image_62.png" /></figure>
# Annexe 10 : Guide Méthodologique d'Analyse spatiales des potentialités de projets de REUT
# PROJET SOPOLO : LA REUTILISATION DES EAUX USEES TRAITEES : UNE SOLUTION POSSIBLE LOCALEMENT ?

## PHASE 1 : ANALYSE SPATIALES DES POTENTIALITES DE PROJETS DE REUT

### Guide Méthodologique (version QGIS 2.18.XX)

- Mars 2019 -

Florent Hubert - Myriam Campardon
Le présent guide méthodologique a pour but d'offrir une présentation détaillée de la méthode développée dans le cadre du projet SoPoLo, relative à l'identification de site d'intérêt pour des projets de réutilisation d'eaux usées traitées à des fins d'irrigation. Ce document fait écho au premier rapport du projet SoPoLo (La REUT, une solution possible localement?) concernant une " Méthodologie d'analyse spatiale des potentialités locales de réutilisation des eaux usées traitées ". Contrairement au rapport, ce guide ne présentera pas de point de contexte ni d'analyse de résultats, en revanche seront présenté et illustrés les sources de données, principales définitions, étapes et outils utilisés. Cette version du guide présente le déroulé de la méthode au travers des logiciel QGIS (en version 2.18). Sont également utilisés (dans une moindre mesure) les logiciels R (3.4.X) et Microsoft Excel.

## GLOSSAIRE

**STEU** : Système de Traitement des Eaux Usées

**REUT** : Réutilisation des Eaux Usées Traitées (à des fins d'arrosage ou d'irrigation)

**Site** : Dans le cadre présent un site représente une parcelle agricole ou un espace vert/sportif potentiel pour un projet de REUT.

**Couche** : On entend ici par couche un ensemble logique d'entités définis par leur forme, position et attributs. Une couche peut être vectorielle ou raster.

**Entité** : Unité indivisible d'une couche, pour les couches vectorielles ce peut être un polygone, une ligne ou un point, pour les couches raster c'est un pixel.

**Raster** : Couche « image » géoréférencée. Elle est composée d'un nombre entier de pixels de taille définie, traditionnellement carrés.

**Vecteur** : Objet géographique (surfacique, linéaire ou ponctuel) délimité par des points (nodes, sommets) géo-localisés et reliés entre eux (polygones et lignes). Ces objets présentent des attributs rassemblés au sein d'une table attributaire.

**Système de coordonnées** : Système de référence permettant la localisation des objets géographiques selon des axes (3D ou 2D pour les systèmes projetés).

**Table attributaire** : Support de l'information sémantique des objets géographiques. Contient l'ensemble des attributs des entités des couches géographiques (ex : Aire, Périmètre, etc.).

**Critères, facteurs, contraintes** : Dans son ouvrage « Analyse spatiale de l'information géographique », Régis Caloz définit le critère comme « l'élément de base d'une décision relative à l'objectif de l'analyse décisionnelle. [...] ces critères se distribuent en deux groupes bien distincts : les facteurs et les contraintes. » (Caloz, Collet, 2011). Un facteur se comprend comme un critère influençant l'aptitude d'un objet (relative à l'atteinte d'un objectif) tandis qu'une contrainte discrédite ou non ce dernier en fonction de la modalité du critère. Si pour une analyse plusieurs critères sont pris en compte, on parle alors d'analyse multicritère.
## SOMMAIRE

<table><tr><td>Présentation du guide</td><td>2</td></tr><tr><td>Glossaire</td><td>2</td></tr><tr><td>Sommaire</td><td>3</td></tr><tr><td>Objet d'étude</td><td>4</td></tr><tr><td>Présentation QGIS</td><td>5</td></tr><tr><td>Présentation générale</td><td>5</td></tr><tr><td>Modeleur graphique</td><td>7</td></tr><tr><td>Présentation de la méthodologie</td><td>9</td></tr><tr><td>Concept général et schémas</td><td>9</td></tr><tr><td>Sources de données</td><td>10</td></tr><tr><td>Préparation des données</td><td>11</td></tr><tr><td>Modèles</td><td>15</td></tr><tr><td>Limites et critiques de l'exercice</td><td>32</td></tr><tr><td>Et ensuite ?</td><td>33</td></tr><tr><td>Glossaire des algorithmes</td><td>36</td></tr><tr><td>Métadonnées</td><td>42</td></tr></table>
En assainissement, les eaux usées traitées (EUT), généralement restituées au milieu naturel en sortie de station, peuvent-être sous certaines conditions réutilisées à des fins alternatives telles que l'irrigation de parcelles agricoles, d'espaces verts ou encore de golfs. Les stations de traitement des eaux usées (STEU) deviennent alors une source d'eau alternative, concept connu sous le nom de réutilisation des eaux usées traitées (REUT ou wastewater reuse en anglais).

En tant que source d'eau alternative, la REUT suscite un intérêt croissant de la part d'institutions publiques, nationales et internationales, des gestionnaires de l'eau ou encore des organismes de recherche. En effet l'eau propre et l'assainissement font partie des 17 objectifs du développement durable de l'agenda 2030 et l'UNESCO a consacré son rapport mondial sur la mise en valeur des ressources en eau 2017 sur la question de la réutilisation des eaux usées traitées, qu'elle identifie notamment comme « source abordable et durable d'eau, d'énergie, de nutriments, de matières organiques et d'autres sous-produits utiles. » dans un contexte où 80 % des eaux usées produites dans le monde ne sont pas traitées (Programme mondial pour l'évaluation des ressources en eau, 2017).

Motivée par les questions d'économies en eau et de préservation des milieux sensibles, l'Agence de l'Eau Rhône-Méditerranée-Corse (AE RMC) a lancé en 2015 un appel à projet afin de consolider les connaissances sur le sujet.

C'est dans ce contexte que l'UMR G-EAU (IRSTEA Montpellier), en association avec les Communautés de Communes du Grand Pic Saint-Loup et de Montpellier Méditerranée Métropole, la Chambre d'Agriculture de l'Hérault et le Conseil Départemental de l'Hérault a répondu à cet appel en se portant pilote du projet SoPoLo (« La réutilisation des eaux usées traitées : une Solution Possible Localement ? ») visant à « mieux caractériser la demande en EUT, les attitudes et perceptions des décideurs, usagers et consommateurs finaux » sur un territoire test composé des deux communautés de communes. En particulier le manque de méthodologie d'étude globale de ces différents critères est souvent identifié comme un frein majeur pour le développement des projets de réutilisation des eaux usées traitées puisqu'un tel outil faciliterait la compréhension et l'évaluation des différents projets par les décideurs.

Parmi les différentes questions identifiées celle de l'analyse des potentialités spatiales de projet de REUT constitue le premier axe du projet SoPoLo. L'objectif étant d'aboutir à une méthodologie permettant de répondre aux différentes interrogations soulevées en lien avec la dimension spatiale de projets de REUT. La zone d'étude utilisée pour le développement de la méthode est constituée des territoires respectifs de la Communauté de Communes du Grand Pic Saint-Loup et de Montpellier Méditerranée Métropole.

Le présent document a pour objet de présenter les principaux résultats méthodologiques au travers de modèles de traitement (algorithmes) développés sous QGIS 2.18.
L'objet de cette section est de fournir une présentation succincte de QGIS (interface, représentation cartographique, modeleur). Elle se destine donc aux utilisateurs novices, les « habitués » peuvent donc ignorer cette section pour passer à la suivante.

## N.B : La méthode et les algorithmes tels que présentés sont prévus pour fonctionner avec la version 2.18 de QGIS avec GRASS 7.

### PRESENTATION GENERALE

QGIS est un logiciel libre de SIG (Système d'Information Géographique), permettant par définition la visualisation, manipulation et création de données géographiques (géométrique, sémantique et topologique). L'objectif de la première tâche du projet SoPoLo étant de développer une méthode pour l'identification de sites potentiellement intéressants pour la REUT, QGIS est un outil particulièrement adapté pour cet exercice. Ce dernier peut être téléchargé sur son site officiel à l'adresse suivante : https://www.qgis.org/fr/site/.

### DOCUMENTATION

QGIS dans sa version 2.18 dispose d'une documentation relativement conséquente néanmoins pas intégralement traduite en français : https://docs.qgis.org/2.18/fr/docs/index.html

Notamment sont très utile pour tous novices :

<table><tr><td>Premiers pas</td><td>https://docs.qgis.org/2.18/fr/docs/user_manual/introduction/getting_started.html</td></tr><tr><td>Présentation de l'interface</td><td>https://docs.qgis.org/2.18/fr/docs/user_manual/introduction/qgis Gui.html</td></tr><tr><td>Outils généraux</td><td>https://docs.qgis.org/2.18/fr/docs/user_manual/introduction/general_tools.html</td></tr><tr><td>Gestion des données</td><td>https://docs.qgis.org/2.18/fr/docs/user_manual/managing_data_source/index.html</td></tr><tr><td>Boîte à outils</td><td>https://docs.qgis.org/2.18/fr/docs/user_manual/processing/toolbox.html</td></tr></table>

Il est recommandé de prendre en main le logiciel avant de se lancer dans les manipulations directes de la méthode. Pour ce faire il existe de nombre tutoriels sur le net, en anglais ou en français avec des jeux de données test (ex : http://cours-fad-public.ensg.eu/, site de l'ENSG - http://sigea.educagri.fr/tutoriels-de-logiciels-sig/qgis/, site de SIGEA).

## CARTES : STYLE ET MISE EN PAGE

QGIS permet de produire des cartes à l'aide d'un composeur d'impression, permettant la mise en forme des données, l'ajout d'éléments (tables, échelle, cadre, texte, rose des vents, etc.) et leur export (format pdf, png,svg et autres). Celui-ci est accessible à partir de l'onglet Projet ou en raccourci avec les touches Ctrl+P.

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

Mais avant d'explorer le composeur il est important de savoir que l'on peut sous QGIS modifier l'apparence des jeux de données, qu'ils soient vecteur ou raster, selon leur propriétés (table attributaire ou valeur des pixels). Cette apparence peut être modifiée dans l'onglet Style des propriétés des couches (accessible dans les propriétés par clic droit sur la couche) ou dans le
panneau Style de Couche (ouvert par défaut dans QGIS, sinon dans l'onglet Vue -> Panneaux -> Style de Couche).

## Onglet Style pour une couche vectorielle :

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

Pour les vecteurs le principe du Style est de générer des symboles qui seront l'apparence des entités à l'écran. Ce symbole est classiquement composé d'un fond et d'un contour, mais il est possible d'y ajouter de nombreux éléments (lignes, points, étoile, etc) ou de modifier les couleurs (dégradé, transparence, épaisseur de contour, etc). On trouve également des symboles préconfigurés. N'hésitez pas à explorer les possibilités, c'est le meilleur moyen de comprendre le fonctionnement de l'outil.

Il est possible de définir des styles catégorisés (pour des attributs nominaux ou ordinaux, par exemple un type de sol, un nombre de culture) et des styles gradués (pour des attributs continus, comme une hauteur de précipitation annuelle).

## Onglet Style pour une couche raster :

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

Le style des couches raster se base sur les valeurs des pixels d'une ou de plusieurs bandes (type photo aérienne avec bandes RVB). Il est possible de choisir un rendu en bande grise unique
(dégradé noir/blanc), en bandes multiples (RVB), en palette (si les couches raster possède une information de palette intégrée), ombrage ou pseudo-couleur à bande unique (carte de chaleur, occupation du sol) où il est possible de régler les couleurs par valeur ou intervalle.

Il est possible de modifier les options de rendu de couche (teinte, saturation, contraste, etc) et de
ré-échantillonnage (interpolation entre pixels, moyenne) permettant de « lisser » plus ou moins
le raster.

MODELEUR GRAPHIQUE

Si vous avez explorés différents tutoriels de QGIS vous avez certainement dû être confrontés à des traitements, ces derniers étant des algorithmes manipulant des données en entrée pour produire un résultat (sortie). Par exemple l'algorithme Pente permet de calculer la pente (en % ou °) à chaque pixel à partir d'un MNT (Modèle Numérique de Terrain, raster « d'altitude ») ou encore Distance tampon fixe permettant de créer une zone tampon d'une taille déterminée autour d'un vecteur.

Si dans le cas général l'utilisateur de QGIS utilise ces traitements de manière anecdotique il existe des situations où il faut répéter plusieurs traitements un grand nombre de fois. Auquel cas la manipulation se révèle longue et répétitive.

Pour pallier à cette éventualité QGIS met à disposition des utilisateurs un modeleur graphique,
outil permettant de créer des modèles où l'on définit entrées, algorithmes et sorties. Pour ceux
connaissant le terme, ce modeleur se rapproche d'un ETL (Extract, Transform, Load). Ce dernier
ce présente dans la forme suivante (y accéder en cliquant sur Create new model dans la boîte à
outils de traitements) :

<figure><img src="image_66.png" /></figure>
On y trouve une fenêtre de visualisation de l'algorithme, à gauche un onglet où l'on trouve les différents traitements disponibles sous QGIS ainsi que les entrées définissables (passer d'un onglet à l'autre en cliquant sur les noms en bas). Sur la partie haute les différents boutons permettent l'enregistrement des modèles, leur export en image (pour documentation) ou en script python (pour retravailler avec plus de flexibilité moyennant des connaissances python). Le bouton permet d'éditer l'aide du modèle (qui s'affiche sur la droite en lançant le modèle et explique le but du modèle ainsi que les entrées/sorties prévues). Enfin lance le modèle.

Dans le modèle les entrées et traitements présentent des icônes , la croix permettant de supprimer l'entrée/traitement et le crayon pour la/le modifier. L'aide de QGIS pour le modeleur est disponible à l'adresse suivante : https://docs.qgis.org/2.18/fr/docs/user_manual/processing/modeler.html

Dans le cadre de la première phase du projet SoPoLo, un certain nombre de modèles ont été
créés pour faciliter la mise en place et l'application de la méthode.
Cette section présente la méthodologie développée dans le cadre du projet SoPoLo, en abordant tout d'abord le concept sous-jacent avec l'aide de schémas, quelques remarques, les sources de données utilisées, leur préparation et les modèles.

## CONCEPT GENERAL ET SCHEMAS

L'objectif de la première phase du projet SoPoLo est de permettre l'identification de sites potentiellement intéressant pour la REUT. Ces sites, surtout représentés par le parcellaire agricole, doivent présenter des caractéristiques intéressantes pour favoriser les chances de réussite de projets de REUT. Ils doivent de plus rester accessibles par rapport à la station d'épuration pour limiter les coûts liés au transport de l'eau entre le STEU et la parcelle.

La méthode repose donc sur deux axes (ou indicateurs) principaux :

- L'intérêt (construit à partir de critères liés à la culture) reflétant les caractéristiques liées au site pouvant être recherchées dans un projet de REUT (besoin en eau, usages, gains, inégalités spatiales d'accès à l'eau)

- L'accessibilité au STEU, reflétant la facilité avec laquelle on peut rejoindre le site à partir de la station, dépend de la topographie et de l'occupation du sol.

Ces deux indices sont utilisés pour évaluer une potentialité (reliant les propriétés de la parcelle dans l'absolu à son positionnement par rapport à la STEU) à partir de laquelle on opère une sélection de site et identifie un chemin de moindre coût entre STEU et site :

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

On distingue ainsi plusieurs phases au sein de cette méthode :

- La construction des critères, qui passe par une normalisation des données en vue de leur croisement ainsi que de l'exploitation des données de bases pour leur calcul.

- Le calcul des deux indicateurs d'intérêt et d'accessibilité

- Leur combinaison pour obtenir l'indicateur de potentialité
- La sélection de site et la construction d'un chemin reliant site(s) à STEU.

Afin de faciliter la mise en place de la méthode la majeure partie de la méthode est sous forme d'algorithmes construits avec le modeleur graphique de QGIS. Ces algorithmes sont multiples pour pouvoir visualiser les résultats intermédiaires et faciliter le débogage en cas de problèmes.

Les différentes phases de la méthode seront précisées au fur et à mesure de la présentation des modèles. Les couches seront indiqués en **souligné**, les traitements en **italique**, les champs « entre guillemets » et les modèles **en gras**.

## SOURCES DE DONNEES

Différentes sources de données ont été exploitées dans le cadre de la méthode, la majeure partie étant gratuite ou produite par des organismes publics (IGN, INRA, etc.) :

<table><thead><tr><td>Couche</td><td>Type</td><td>Année</td><td>Echelle</td><td>Intérêt</td><td>Source/fournisseur</td></tr></thead><tbody><tr><td>Registre Parcellaire Graphique (niveau 4)</td><td>Vecteurs</td><td>2014</td><td>1 : 5 000</td><td>Données sur les cultures à l'échelle parcellaire</td><td>ASP</td></tr><tr><td>BD Topo®</td><td>Vecteurs</td><td>2016</td><td>1 : 5 000</td><td>Entités géographiques et administratives du territoire</td><td>IGN</td></tr><tr><td>Station d'épuration des eaux usées</td><td>Points</td><td>2016</td><td>-</td><td>Localisation et informations sur les stations d'épuration</td><td>Ministère en charge de l'Environnement</td></tr><tr><td>Occupation du sol 3M</td><td>Vecteurs</td><td>2012</td><td>1 : 5 000</td><td>Espaces verts et structures sportives</td><td>3M</td></tr><tr><td>Occupation du sol CCGPSL</td><td>Vecteurs</td><td>2012</td><td>1 : 5 000</td><td>Espaces verts et structures sportives</td><td>CCGPSL</td></tr><tr><td>BD Carthage®</td><td>Vecteurs</td><td>2014</td><td>1 : 50 000</td><td>Référentiel national pour le repérage des entités hydrographiques de surfaces</td><td>IGN</td></tr><tr><td>Référentiel masse d'eau DCE</td><td>Vecteurs</td><td>2016</td><td>1 : 5 000 à 1 : 50 000</td><td>Référentiel d'évaluation de la DCE en France</td><td>Sandre</td></tr><tr><td>Grilles SAFRAN</td><td>Vecteurs</td><td>2016</td><td>-</td><td>Données météorologies interpolées sur le territoire d'étude</td><td>Météo France</td></tr><tr><td>BDSOL L-R</td><td>Vecteurs BDD Access</td><td>2012</td><td>1 : 250 000</td><td>Données pédologiques</td><td>INRA Orléans (diffusion SIG L-R)</td></tr><tr><td>Documents d'urbanismes</td><td>Vecteurs</td><td></td><td>1 : 5 000 à 1 : 50 000</td><td>Zones urbaines et à urbaniser</td><td>3M, CCGPSL</td></tr><tr><td>Indice de Qualité des Sols</td><td>Vecteurs</td><td>2015</td><td>1 : 250 000</td><td>Indicateur de potentiel agronomique des sols pour les cultures</td><td>DRAAF Occitanie/ BDSol L-R</td></tr><tr><td>Périmètres de Protection Rapprochés</td><td>Vecteurs</td><td>-</td><td></td><td>Délimitation des périmètres de protections rapprochés des captages AEP</td><td>3M</td></tr><tr><td>MNT</td><td>Raster</td><td>2013</td><td>5m</td><td>Données sur l'altitude</td><td>RGE Alti® IGN</td></tr><tr><td>BD Forêt V1</td><td>Vecteur</td><td>2014</td><td>1 : 50 000</td><td>Délimitation des espaces boisés</td><td>IGN</td></tr><tr><td>OSO CESBIO</td><td>Raster</td><td>2016</td><td>20m</td><td>Occupation du sol (télédétection)</td><td>CESBIO</td></tr></tbody></table>

Une présentation plus détaillée des couches est disponible dans la section Métadonnées.
# FICHIER TEXTES ET EXCELS

Au cours de la méthode qui va être exposée un certain nombre de fichiers textes et Excel seront utilisés, la plupart d'entre eux sont formatés et certains contiennent des scripts permettant d'automatiser certains processus longs.

Les principales tables et fichiers (avec format) sont décrits ci-dessous :

<table><thead><tr><td>Nom du fichier</td><td>Description</td></tr></thead><tbody><tr><td>Interet_relatif_R.xlsm</td><td>Tableau à double entrée pour l'intérêt relatif</td></tr><tr><td>Access_reclass.txt</td><td>Table de reclassement de l'accessibilité</td></tr><tr><td>Cesbio_reclass.txt</td><td>Table de reclassement nomenclature CESBIO vers occupation du sol</td></tr><tr><td>Contr_IQS.txt</td><td>Table de jointure pour les contraintes de l'indice de qualité des sols</td></tr><tr><td>IQS_reclass.txt</td><td>Table de reclassement de l'indice de qualité des sols</td></tr><tr><td>Occup_reclass.txt</td><td>Table de reclassement de l'occupation du sol</td></tr><tr><td>Pente_reclass.txt</td><td>Table de reclassement des pentes</td></tr><tr><td>Table_deficit.txt</td><td>Table de reclassement du déficit</td></tr><tr><td>Table_exp.txt</td><td>Table de reclassement du nombre d'exploitant</td></tr><tr><td>Table_usage.txt</td><td>Table de reclassement du type d'usage</td></tr><tr><td>Rpg_usage.csv</td><td>Table de jointure entre RPG et type d'usage agricole</td></tr><tr><td>Largeur_route.csv</td><td>Table de jointure contenant les informations d'emprise routière</td></tr></tbody></table>

# PREPARATION DES DONNEES

# REMARQUES GENERALES ET PRECAUTIONS

Veillez à lancer le logiciel en tant qu'administrateur sans quoi les modèles ne fonctionneront pas. Avant l'application de la méthode, veuillez à décocher l'option « Ignorer la déclaration interne d'encodage des shapefiles » dans les options QGIS :

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

Désactivez également la reprojection à la volée dans les options SCR. Veillez à ce que les fournisseurs de traitement GRASS, SAGA, GDAL/ORG et R soient activés et paramétrés.
Les traitements utilisés sous QGIS fonctionnant sous Python et PyQt, certaines erreurs peuvent être occasionnées si les fichiers ou leur chemin d'accès contiennent accents et/ou espaces. Il est recommandé d'éviter ces derniers lors de la création/manipulation de données.

## FOURNISSEURS DE TRAITEMENT

Dans sa version 2.18 QGIS intègre différents fournisseurs de traitement dans son environnement (GRASS, SAGA, GDAL, OGR pour citer les principaux). Ceux-ci doivent être activés dans le menu Traitement > Options. Rechercher les fournisseurs GRASS GIS 7, SAGA, R script s et GDAL/OGR et cliquer sur Activate :

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

## EXTENSIONS

QGIS permet l'installation d'extensions permettant l'ajout de nouvelles fonctionnalités et/ou traitements pouvant être utiles. Dans le cadre de cette méthode a été fait appel à l'extension *Concave Hull* permettant le calcul d'une enveloppe concave autour d'entités avec possibilité d'utiliser la table attributaire pour définir différentes enveloppes.

L'installation d'extensions se fait par l'onglet Extension > Installer/Gérer les extensions. Rechercher Concave Hull et cliquer sur installer.

Note : L'extension *Group Stats* est également utile pour l'analyse des résultats produits en permettant le calcul de statistiques à partir des données de table attributaire (sous forme de table pivot comme sous Excel).
## CHARGER LES MODELES ET SCRIPTS

Les différents modèles et scripts utilisés pour cette méthode doivent être ajoutés aux outils de traitement de QGIS pour pouvoir les exécuter.

Dans la fenêtre boîte à outils de traitements il est possible d'ajouter modèles et scripts python à partir des onglets « Modèles » et « Scripts » en cliquant sur « Add model (script) from file ».

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

Les scripts R doivent eux être copiés dans le répertoire de scripts R (format .rsx) :
« C:\Users\Utilisateur\.qgis2\processing\rscripts »

Une fois chargés les modèles apparaissent de la façon suivante dans l'onglet « Modèles » :

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

## TABLES DE RECLASSEMENT

Différentes tables de reclassement (permettant la modification des valeurs d'une couche suivant des règles prédéfinies) seront utilisées tout au long des modèles (voir section ci-dessus). Celles-ci présentent un formatage particulier, exemple pour reclassement de l'occupation du sol :

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

La syntaxe des règles est la suivante : `VALEUR_MIN:VALEUR_MAX:VALEUR_NOUVELLE`
On définit des intervalles (valeurs minimales et maximales) pour la couche d'entrée et indique la valeur à renseigner en sortie. A noter que les intervalles peuvent contenir une seule valeur unique. Ou au contraire l'utilisation d'une « * » indique une borne infinie.

La syntaxe suivante existe aussi : VALEUR_MIN : VALEUR_MAX : NOUVEAU_MIN : NOUVEAU_MAX

Cette syntaxe permet par exemple de changer d'unité ou d'échelle dans les données raster (ex : Passer de °C en °F pour une carte de température moyenne)

Des tables de reclassement « par défaut » sont fournies dans les livrables, mais il est tout à fait possible d'en créer en suivant la syntaxe décrite auparavant.

## CALCUL DU DEFICIT HYDRIQUE

Parmi les critères utilisés lors du développement de la méthode, le déficit hydrique des cultures est sans doute le plus complexe à obtenir.

Son calcul nécessite trois types de données : La réserve utile en eau des sols, les données de précipitation et d'évapotranspiration foliaire et le coefficient cultural des plants cultivés. Ces données sont fournies par défaut pour la zone test.

Le formatage est le suivant :

### Table de réserve utile :

<table><thead><tr><td>Nom de l'UCS</td><td>Valeur de RU (mm)</td></tr></thead><tbody><tr><td>104A</td><td>250</td></tr></tbody></table>

### Données météo :

<table><thead><tr><td>ETP (mm)</td><td>Pluie (mm)</td></tr></thead><tbody><tr><td>4.2</td><td>3.1</td></tr></tbody></table>

### Kc (coefficient cultural) :

<table><thead><tr><td>Code culture/n°jour</td><td>1 (1er janvier)</td><td>2</td><td>3</td></tr></thead><tbody><tr><td>21 (Vigne)</td><td>0</td><td>0.1</td><td>0.1</td></tr><tr><td>3 (Orge)</td><td>0.675</td><td>0.675</td><td>0.675</td></tr></tbody></table>

Les tables doivent être sous format .csv (séparateur point-virgule).

Exécuter alors sous R le script déficit en renseignant répertoire où se trouvent les tables, et le nom des tables.

Chemin (modifier à convenance) : setwd("C:/SoPoLo/tables/Deficit") (set working directory).

Nom des tables (modifier à convenance)

- Kc <- read.csv2(file = "kc_jour.csv"

- RU <- read.csv2(file = "RU.csv"

- Meteo <- read.csv2(file="grid.csv"
Sélectionner tout le code et cliquer sur « Run ». Un fichier contenant les déficits par combinaison de type de sol et de culture est créé dans le fichier renseigné. Un identifiant unique type 1 UCS CULTURE est indiqué.

## MODELES

Cette sous-section fait la présentation des différents modèles utilisés dans le cadre de la méthode. Le schéma suivant illustre de manière générale la façon dont les modèles s'enchaînent et se complémentent :

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

### POINT PREALABLE : EXECUTION PAR LOT

Pour chacun des modèles présentés dans cette section il est demandé de définir des entrées (couches, tables, champs, valeurs, etc...) et des sorties. Néanmoins il peut être souhaitable d'appliquer la méthode sur différents sites ou pour différents critères. Pour ce faire il est possible de lancer le modèle en « processus de lot », une opération qui itère l'algorithme autant de fois que l'on souhaite et permet de définit entrées et sortie pour chacune des itérations. Il suffit de cliquer sur le bouton « Exécuter comme processus de lot... » en haut à droite de la fenêtre de paramétrage du traitement :
<figure><img src="image_74.png" /></figure>

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

**Note:** Lors du paramétrage d'une exécution en traitement par lot pour l'algorithme **Calcul de potentialité**, veillez à remplacer les « \ » par des « / » pour le chemin d'accès à la table de croisement et seulement pour celui-ci. Ceci est dû à une spécificité du langage R qui est utilisé. Ex : « C:\Documents\table_potentialite.csv » vers « C:/Documents/table_potentialite.csv »

## OCCUPATION DU SOL

Comme mentionné dans la présentation du concept général de la méthode, l'occupation du sol est une composante requise dans le calcul d'indicateur d'accessibilité. Bien qu'il existe de nombreuses sources de données sur l'occupation du sol, la présente méthode requiert une information précise (échelle communale ou intercommunale) et qui puisse retranscrire une accessibilité, autrement dit la nomenclature doit présenter des classes permettant d'associer facilement une valeur de « friction » ou de contrainte de passage pour l'établissement d'un réseau d'eau brute par exemple. De plus l'information doit être facilement accessible et couvrir le territoire national.

Face au manque de données répondant à ces différents critères il a été choisi d'élaborer des modèles exploitant les données BD Topo®, BD Forêt® et BD Carthage®, des PLU (Plans Locaux d'Urbanisme) et d'un raster d'occupation du sol (OSO CESBIO) comme base pour construire une occupation du sol théorique pouvant être exploitée par la suite.

La nomenclature par défaut utilisée lors des tests est la suivante :

<table><thead><tr><td>Code</td><td>Libellé</td><td>Source</td></tr></thead><tbody><tr><td>1</td><td>Inconnu</td><td>-</td></tr><tr><td>2</td><td>Cultures</td><td>OSO CESBIO</td></tr><tr><td>3</td><td>Vigne</td><td>OSO CESBIO, BD Topo®</td></tr><tr><td>4</td><td>Vergers</td><td>OSO CESBIO</td></tr><tr><td>5</td><td>Chemins et sentier</td><td>BD Topo®</td></tr><tr><td>6</td><td>Route empierrée</td><td>BD Topo®</td></tr><tr><td>7</td><td>Route goudronnée petite</td><td>BD Topo®</td></tr><tr><td>8</td><td>Route goudronnée grande</td><td>OSO CESBIO, BD Topo®</td></tr><tr><td>9</td><td>Chemin de fer</td><td>BD Topo®</td></tr><tr><td>10</td><td>Landes</td><td>BD Forêt®</td></tr></tbody></table>
<table><tr><td>11</td><td>Forêt ouverte</td><td>BD Forêt®</td></tr><tr><td>12</td><td>Peupleraie</td><td>BD Forêt®</td></tr><tr><td>13</td><td>Taillis et futaies</td><td>BD Forêt®</td></tr><tr><td>14</td><td>Autres boisements</td><td>OSO CESBIO</td></tr><tr><td>15</td><td>Zone urbaine</td><td>OSO CESBIO, PLU</td></tr><tr><td>16</td><td>Affleurements</td><td>OSO CESBIO</td></tr><tr><td>17</td><td>Masse d'eau</td><td>BD Topo®</td></tr><tr><td>18</td><td>Prairies et pelouses</td><td>OSO CESBIO</td></tr><tr><td>19</td><td>Sables</td><td>OSO CESBIO</td></tr><tr><td>20</td><td>Très petit cours (0 à 15m)</td><td>BD Carthage®</td></tr><tr><td>21</td><td>Petit cours (15 à 50m)</td><td>BD Carthage®</td></tr><tr><td>22</td><td>Cours moyen (Plus de 50m)</td><td>BD Carthage®</td></tr><tr><td>23</td><td>Neige et glaciers</td><td>OSO CESBIO</td></tr></table>

## A) RESEAU SURFACIQUE

Le réseau routier constitue un élément important dans l'analyse d'accessibilité dans la mesure où il peut constituer une contrainte forte (autoroute) ou au contraire une aubaine (chemin, sentier). Néanmoins l'information fournie dans la BD Topo® est linéaire, or pour pouvoir comparer l'ensemble des sources de données et calculer l'accessibilité il nous faut une information surfacique. De plus sous le nom « route » se cache une grande variabilité de voies de communication de revêtement et/ou de largeur variable. Il faut donc pouvoir exploiter ces deux informations.

Le modèle **Route surfacique** utilise en entrée le réseau routier, qu'il découpe à la zone d'étude voulue, puis joint une table contenant des informations de largeur. Une zone tampon variable est ensuite réalisée en utilisant la largeur renseignée par la table. On sépare enfin les routes des chemins. Est ajouté également un champ « cat » pour la nomenclature d'occupation permettant de différencier certaines sous catégories (chemin, route emperrée, route large, etc.).

<table><thead><tr><td>Entrées</td><td>Description</td></tr></thead><tbody><tr><td>Routes</td><td>Description du réseau routier par la BD Topo® (ROUTE.shp)</td></tr><tr><td>Chemin de fer</td><td>Description du réseau de voies de fer par la BD Topo® (TRONCONVOIE_FERREE.shp)</td></tr><tr><td>Zone d'étude</td><td>Zone d'étude</td></tr><tr><td>Table buffer</td><td>Table de jointure contenant les largeurs de route corrigées par le CEREMA (largeur route.csv)</td></tr></tbody></table>
<figure><img src="image_76.png" /></figure>

<table><thead><tr><td>Sorties</td><td>Description</td></tr></thead><tbody><tr><td>Routes surfaciques</td><td>Surfaces routières</td></tr><tr><td>Chemins surfaciques</td><td>Surfaces de chemins, sentiers et routes empierrées</td></tr><tr><td>Chemin de fer surfacique</td><td>Surfaces de voies ferrées</td></tr></tbody></table>

On notera qu'à ce stade le chemin de fer n'est pas découpé à la zone d'étude, en effet si l'on réalisait cette opération dans le modèle on risquerait de générer une sortie vide (toutes les zones d'études n'ont pas de voies ferrées !) et donc de faire planter le modèle.

## B) OCCUPATION DU SOL

Une fois les précédentes couches générées, il est possible de créer la couche d'occupation du sol. Le modèle **Occupation du sol** s'occupe de ce traitement. Celui-ci utilise une couche raster d'occupation du sol développée par le CESBIO via un algorithme de télédétection. Disposant d'une résolution de 20m, d'une couverture nationale et de statistiques de validation (confiance et validité), cette couche représente une très bonne base de travail pour le calcul d'accessibilité. Il faut en effet compléter cette couche avec des données vectorielles (parcellaires, forêt) pour affiner
la couche (qui sera utilisée conjointement à un MNT d'une résolution de 5m) et faire apparaître les structures importantes pour l'accessibilité (routes, chemins, zones urbaines). Bien que le schéma du modèle ci-dessous soit difficile de lecture le principe est le suivant :

- Reclassement du raster d'occupation du sol CESBIO pour le faire correspondre à la nomenclature finale

- Chargement des couches vecteurs et ajout, si nécessaire, d'un champ « cat » contenant l'identifiant d'occupation selon la nomenclature finale

- Rastérisation des couches vectorielles

- Fusion des différents rasters les uns sur les autres selon un ordre prédéfini

L'ordre de fusion des couches est le suivant (superposition) : Chemins, routes, voies ferrées, réseau hydrographique, surfaces en eau, zones urbaines, vignes (BD Topo®), zones boisées (BD Forêt®), occupation du sol CESBIO.

<table><thead><tr><td>Entrées</td><td>Description</td></tr></thead><tbody><tr><td>Zone d'étude</td><td>Zone d'étude</td></tr><tr><td>OCS Raster</td><td>Occupation du sol raster (20m, CESBIO)</td></tr><tr><td>Table de reclassement OCS</td><td>Table de reclassement de l'occupation du sol (CESBIO) vers la nomenclature finale</td></tr><tr><td>Zone urbaine</td><td>Délimitation des zones urbaines (PLU)</td></tr><tr><td>Végétation Topo</td><td>Couvert végétal selon la BD Topo® (ZONE VEGETATION.shp)</td></tr><tr><td>Voies ferrées</td><td>Surfaces en voies de fer (créé par modèle)</td></tr><tr><td>Routes</td><td>Surfaces en route (créé par modèle)</td></tr><tr><td>Chemins</td><td>Surfaces en chemins et sentiers (créé par modèle)</td></tr><tr><td>Surfaces en eau</td><td>Délimitation des surfaces en eau selon la BD Topo® (SURFACE_EAU.shp)</td></tr><tr><td>Réseau hydrographique</td><td>Tracé du réseau hydrographique selon la BD Carthage® (TRONCON_HYDROGRAPHIQUE.shp)</td></tr><tr><td>Forêt</td><td>Couvert végétale selon la BD Forêt® v1</td></tr></tbody></table>

<figure><img src="image_77.png" /></figure>
<table><thead><tr><td>Sorties</td><td>Description</td></tr></thead><tbody><tr><td>OCS Final</td><td>Occupation du sol raster sur la zone d'étude</td></tr></tbody></table>

## ACCESSIBILITE

Une fois l'occupation du sol rastérisée elle peut être utilisée pour le calcul d'accessibilité. Celle-ci se définit ici par la facilité avec laquelle on rejoint le STEU à partir de n'importe quel point de la zone d'étude. L'objectif est ici de produire une carte de coût (dont chaque pixel représente le coût minimal pour rejoindre le STEU) et une carte d'antécédence (dont chaque pixel représente la direction à prendre pour rejoindre le STEU suivant une logique de moindre coût).

Le processus est le suivant :

- Découpage d'un MNT et de l'occupation du sol à la zone d'étude

- Calcul des pentes à partir du MNT

- Reclassement des pentes pour isoler les zones de pentes fortes (supérieures à 45° par défaut)

- Ajout des zones de pentes fortes à l'occupation du sol

- Reclassement de l'occupation du sol obtenue en carte de friction (chaque pixel possède une valeur de coût)

- Calcul des cartes de coût et d'antécédente (traitement r.walk.rast de GRASS)

<table><thead><tr><td>Entrées</td><td>Description</td></tr></thead><tbody><tr><td>Zone d'étude</td><td>Zone d'étude</td></tr><tr><td>OCS</td><td>Couche d'occupation du sol rastérisée</td></tr><tr><td>STEU</td><td>Système de traitement des eaux usées (point)</td></tr><tr><td>MNT</td><td>Modèle Numérique de Terrain (5m)</td></tr><tr><td>Table de reclassement des pentes</td><td>Table de reclassement des pentes fortes</td></tr><tr><td>Table de reclassement OCS</td><td>Table de reclassement de l'occupation en friction</td></tr><tr><td>Distance buffer STEU</td><td>Taille de la zone tampon (« emprise » de la STEU)</td></tr><tr><td>Table de reclassement du coût</td><td>Table de reclassement du coût en indicateur d'accessibilité</td></tr></tbody></table>

<figure><img src="image_78.png" /></figure>
<table><thead><tr><td>Sorties</td><td>Description</td></tr></thead><tbody><tr><td>Friction</td><td>Carte de friction</td></tr><tr><td>Coût</td><td>Carte de coût</td></tr><tr><td>Antécédence</td><td>Carte d'antécédence</td></tr><tr><td>Accessibilité</td><td>Carte d'accessibilité</td></tr></tbody></table>

RESTRICTION ZONE D'ETUDE

Cette étape, purement facultative, permet de définir une zone d'étude plus restreinte à partir de l'accessibilité. L'objectif étant de pouvoir éventuellement compléter les données disponibles (parcellaire notamment) en se focalisant sur les zones à priori plus accessibles.

Le processus est le suivant :

- **Reclassement** de la **carte de coût** en **indicateur** (1, très forte, à 6, très faible)

- **Extraction** des zones d'indice inférieur à un seuil d'accessibilité.

- **Polygonisation** des **zones d'accessibilité voulue**

- **Dissolution/simplification** du résultat, **comblement** des trous, **enveloppe concave**

<table><thead><tr><td>Entrées</td><td>Description</td></tr></thead><tbody><tr><td>Coût</td><td>Carte de coût</td></tr><tr><td>Seuil</td><td>Seuil d'accessibilité</td></tr><tr><td>Table de reclassement coût</td><td>Table de reclassement de la carte de coût en indicateur</td></tr></tbody></table>
<figure><img src="image_79.png" /></figure>

<table><thead><tr><td>Sorties</td><td>Description</td></tr></thead><tbody><tr><td>Zone restreinte</td><td>Zone d'étude restreinte</td></tr></tbody></table>

CRITERE

Une fois l'accessibilité calculée et mise en forme d'indicateur, il convient de préparer les données permettant le calcul de l'intérêt. Celui-ci se fait par le biais du croisement de différents critères dont nous allons exposer les modalités ici.

Pour rappel l'intérêt reflétant les caractéristiques liées au site pouvant être recherchées dans un projet de REUT (besoin en eau, usages, gains, inégalités spatiales d'accès à l'eau). Il est donc une somme pondérée des critères le définissant, eux-mêmes ordonnés selon les objectifs recherchés (ex : favoriser les besoins en eau moyens, éviter les zones avec peu d'usage). Il y a donc une double balance : au sein des critères et entre critères.

Si certains critères sont propre au parcellaire, d'autre font appel à des caractéristiques d'environnement (proximité) comme par exemple le nombre d'exploitants autour de chaque parcelle.
Le script critère permet pour chaque critère d'opérer un reclassement de 1 à 9 selon les situations plus ou moins préférentielles et l'export du résultat en raster pour le croisement entre critères. Il est également possible de choisir entre critère de proximité ou non. L'option de critère de proximité permet de réaliser un compte des valeurs uniques dans un voisinage de taille défini (option Distance en pixel de raster, par défaut 61 pour 150m : 1 pixel de centre et 30 pixels de 5m de chaque côté).

<table><thead><tr><td>Entrées</td><td>Description</td></tr></thead><tbody><tr><td>Zone d'étude</td><td>Zone d'étude</td></tr><tr><td>Critère</td><td>Couche vecteur contenant l'information du critère</td></tr><tr><td>Champ du critère</td><td>Champ de la couche vecteur contenant l'information du critère</td></tr><tr><td>Table de reclassement du critère</td><td>Table de reclassement des valeurs du critère (ou de leur diversité) vers une note (1 à 9).</td></tr><tr><td>Proximité</td><td>Boolean (critère de proximité ou non)</td></tr><tr><td>Distance (cellules)</td><td>Si le critère est de proximité, quel diamètre (en nombre de cellules) est utilisé pour le calcul de diversité ? (nb_cellule = 2*rayon(m)/resolution_raster+1)</td></tr></tbody></table>

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

<table><thead><tr><td>Sortie</td><td>Description</td></tr></thead><tbody><tr><td>Critère final</td><td>Critère reclassé de 1 à 9</td></tr></tbody></table>

Ex : Pour le critère « Nombre d'exploitants », on renseigne le parcellaire, le champ voulu est un identifiant d'exploitation, le critère est de type « proximité ». On définit également une table de reclassement qui pour chaque compte (nombre d'exploitant) donne une note de 1 à 9.

L'algorithme rastérise alors le parcellaire, chaque pixel ayant une valeur par identifiant d'exploitant, puis compte dans un voisinage de taille définit le nombre de valeurs uniques. Le résultat est ensuite reclassé de 1 à 9 suivant les règles de la table de reclassement. C'est cette dernière couche qui est ensuite sortie sous QGIS.

## ADAPTABILITE : NOUVEAUX CRITERES

Dans le cadre de ce guide 4 critères sont mentionnés : l'IQS (Indice de Qualité des Sols », le déficit en eau, les usages (agricoles et espaces verts/sportifs) et le nombre d'exploitants. Toutefois le
modèle Critère est susceptible de prendre en compte d'autres critères, mais cela sous certaines conditions.

Tout d'abord il faut distinguer deux types de critères : les contraintes et les facteurs. Une contrainte représente un critère binaire (oui/non), il permet l'exclusion ou le maintien d'une solution (au sens de possibilité de site). Un facteur quant à lui influence l'aptitude d'une entité à répondre à un objectif voulu (ex : Une parcelle à être un site intéressant pour la REUT). Ceux-ci ne sont donc pas traités de la même manière. De plus le modèle Critère part du principe que les critères en question sont des facteurs et stockés dans une table attributaire d'une couche vecteur. Dans le détail il peut exister plusieurs cas selon que l'on traite une contrainte, un facteur, qu'il soit vectoriel ou raster. Dans le cas d'une contrainte stockée dans une couche vectorielle on peut également distinguer les cas où la couche vecteur est le parcellaire (les sites potentiels) ou d'autres couches vecteurs (des périmètres de protection par exemple). Le schéma ci-dessous présente de manière générale les différentes situations rencontrées et comment les prendre en compte au plus simple :

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

Les facteurs sous forme vectorielle sont à traiter dans le modèle critère, ceux sous forme raster sont à reclasser (r.recode de GRASS, avec une table de reclassement).

Pour les contraintes raster le plus simple est de vectoriser la couche pour la traiter comme une contrainte vectorielle.
Si une contrainte vectorielle concerne le parcellaire, effectuer une extraction par attribut ou par expression (si plusieurs champ de contrainte) pour ne retenir que les parcelles pour lesquelles la **contrainte (ou combinaison de contraintes) ne s'exprime pas**.

Enfin pour les autres contraintes de type vectoriel, le plus simple est d'extraite par attribut/expression les entités pour lesquelles la **contrainte s'exprime**. S'il existe plusieurs couches de contrainte, les fusionner (merge) et fusionner les entités (dissolve). Le but est d'obtenir une couche unique dont les entités représentent les zones de contraintes. Utiliser ensuite le modèle **Restriction du parcellaire** (voir section Autres modèles). Ce modèle réalise l'intersection des entités du parcellaire avec la couche de contrainte, puis calcule le ratio entre la surface des intersections dans une parcelle (il peut y en avoir plusieurs) et la surface de la parcelle entière. Si ce ratio dépasse un certain seuil (10% par défaut, modifiable) alors la parcelle est écartée.

**Attention** : Veillez à ce que la couche utilisée pour le parcellaire ait une topologie propre (pas de superposition, de doublons, etc) auquel cas le ratio de surface peut être faussé.

<table><thead><tr><td>Entrées</td><td>Description</td></tr></thead><tbody><tr><td>Parcellaire</td><td>Couche du parcellaire</td></tr><tr><td>Contrainte</td><td>Zones de contraintes</td></tr><tr><td>Seuil</td><td>Seuil pour le ratio de surface (entre 0 et 1)</td></tr></tbody></table>

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

<table><thead><tr><td>Sorties</td><td>Description</td></tr></thead><tbody><tr><td>Parcellaire retenu</td><td>Parcelles retenues</td></tr><tr><td>Parcellaire écarté</td><td>Parcelles écartées</td></tr></tbody></table>

Le parcellaire retenu peut désormais être utilisé dans le modèle **Critère** si les informations contenues dans sa table peuvent constituer un facteur (diversité d'exploitants par exemple).
Une fois les différents critères traités il convient de les croiser. Une somme pondérée des différents critères peut être réalisée à cette fin.

Si le nombre de critères est relativement faible (3 à 4), ces pondérations peuvent être déterminées relativement aisément, néanmoins lorsque ce nombre augmente il peut être difficile de donner priorité à tel ou tel critère par rapport à l'ensemble.

On peut faire alors appel à une méthode nommée AHP calculant pour chacun des critères un poids en se basant sur une comparaison par paires de chacun des critères. Exemple : Si l'on dispose de trois critères A, B et C, la détermination des poids de chacun se base sur la comparaison de A à B, A à C et B à C. Il existait sous QGIS des extensions permettant de réaliser ce calcul hélas celles-ci ne sont plus maintenues à jour. On trouve sous R des packages (FuzzyAHP par exemple) permettant de réaliser cette AHP. Sur internet également certains sites offrent cette possibilité : https://bpmsg.com/academic/ahp_CALC.php

Pour le site du bpmsg :

- Renseigner le nombre de critères et cliquer sur OK :

AHP Criteria

Select number of criteria:

Input number and names (2 - 20) 3

Go

OK

- Renseigner le nom des critères :

AHP Criteria Names

Please fill out project title and name of criteria

AHP priorities

Name of criteria

1

A

2

B

3

C

max. 45 character ea.

OK

- Réaliser les comparaisons par paire des critères (1 pour égale importance, 9 pour extrême prédominance) :
Pairwise Comparison AHP priorities

3 pairwise comparisons. Please do the pairwise comparison of all criteria. When completed, click Check Consistency to get the priorities.

Which criterion with respect to AHP priorities is more important, and how much more on a scale 1 to 9?

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

AHP Scale: 1- Equal Importance, 3- Moderate importance, 5- Strong Importance, 7- Very strong importance, 9- Extreme importance (2,4,6,8 values in-between).

- Cliquer sur Check Consistency, le CR (Consistency Ratio) est calculé et vérifie la cohérence des comparaisons, le commentaire à la droite invite éventuellement à recommencer les comparaisons (si supérieur à 10%).

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

Les tableaux en bas renseignent les poids des critères à utiliser :

Priorities

These are the resulting weights for the criteria
based on your pairwise comparisons

<table><thead><tr><th>Category</th><th></th><th>Priority</th><th>Rank</th></tr></thead><tbody><tr><td>1</td><td>A</td><td>69.4%</td><td>1</td></tr><tr><td>2</td><td>B</td><td>13.2%</td><td>3</td></tr><tr><td>3</td><td>C</td><td>17.4%</td><td>2</td></tr></tbody></table>

Number of comparisons = 3

Consistency Ratio CR = 8.4%

Le résultat est exportable en fichier texte (.csv).

Lancer alors la calculatrice raster SAGA et renseigner les critères en rasters d'entrée (un raster principal et des rasters additionnels) :
<figure><img src="image_85.png" /></figure>

Renseigner alors la formule suivante avec les poids calculés :

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

« a » correspond au raster principal, « b » et « c » aux rasters additionnels dans l'ordre de la fenêtre de sélection multiple (b est iqs et c usage dans l'exemple). Adapter le nombre de critères et les poids selon le contexte. Les poids sont des pourcentages, attention la calculatrice prends les points « . » comme séparateurs décimaux et non les virgules ! La fonction « int » corresponds à la partie entière, le +1 sert à arrondir à l'entier le plus proche.

Cliquer sur Run, le résultat obtenu est l'indicateur d'intérêt.

## CALCUL DE POTENTIALITE

Les indicateurs d'accessibilité et d'intérêt étant désormais calculés, il convient désormais de les combiner pour obtenir l'indicateur de potentialité.

Le modèle **Calcul de potentialité** réalise cette opération au travers des étapes suivante :

- Formatage des deux couches d'indicateurs

- Identification des combinaisons uniques de valeurs entre les deux couches

- Utilisation d'une table de reclassement indiquant pour chaque combinaison la valeur de potentialité correspondante

- Création de la couche de potentialité

- Formatage pour sortie

La table de reclassement utilisée est générée à l'aide du tableur Excel *Interet_relatif_r.xlsx*. Il s'agit d'un tableur avec macro (qu'il faut activer) où l'on trouve une feuille « table » et une feuille « Export ». La feuille table comprends un tableau à double entrée renseignant la valeur de potentialité pour chaque combinaison d'accessibilité/intérêt. Remplir ce tableau comme voulu (ne pas modifier les valeurs NoData) puis aller dans la table « Export » pour cliquer sur le bouton
« Exporter vers csv », il suffit alors de renseigner un nom de fichier et la table de reclassement est créée dans le fichier où se trouve le tableur Excel :

<table><thead><tr><th rowspan="2"></th><th colspan="10">ACCESSIBILITE</th></tr><tr><th>1 (Très forte)</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th><th>8</th><th>9 (Nul)</th><th>NoData</th></tr></thead><tbody><tr><th>0 (Nul)</th><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>9999</td></tr><tr><th>1</th><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>9999</td></tr><tr><th>2</th><td>2</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>9999</td></tr><tr><th>3</th><td>3</td><td>2</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>9999</td></tr><tr><th>4</th><td>4</td><td>3</td><td>2</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>9999</td></tr><tr><th>5</th><td>5</td><td>4</td><td>3</td><td>2</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>9999</td></tr><tr><th>6</th><td>6</td><td>5</td><td>4</td><td>3</td><td>2</td><td>1</td><td>0</td><td>0</td><td>0</td><td>9999</td></tr><tr><th>7</th><td>7</td><td>6</td><td>5</td><td>4</td><td>3</td><td>2</td><td>1</td><td>0</td><td>0</td><td>9999</td></tr><tr><th>8</th><td>8</td><td>7</td><td>6</td><td>5</td><td>4</td><td>3</td><td>2</td><td>1</td><td>0</td><td>9999</td></tr><tr><th>9 (Très fort)</th><td>9</td><td>8</td><td>7</td><td>6</td><td>5</td><td>4</td><td>3</td><td>2</td><td>1</td><td>9999</td></tr><tr><th>NoData</th><td>9999</td><td>9999</td><td>9999</td><td>9999</td><td>9999</td><td>9999</td><td>9999</td><td>9999</td><td>9999</td><td>9999</td></tr></tbody></table>

Exporter vers csv

<table><thead><tr><th>Entrées</th><th>Description</th></tr></thead><tbody><tr><td>Accessibilité</td><td>Raster d'accessibilité</td></tr><tr><td>Intérêt</td><td>Raster d'intérêt</td></tr><tr><td>Table de croisement indicateurs</td><td>Table de croisement des indicateurs</td></tr><tr><td>Extent</td><td>Etendue de traitement (laisser par défaut)</td></tr></tbody></table>

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

<table><thead><tr><th>Sorties</th><th>Description</th></tr></thead><tbody><tr><td>Potentialité</td><td>Raster de potentialité</td></tr></tbody></table>

## SELECTION PARCELLE
Une fois l'indicateur de potentialité déterminé, il est possible d'opérer une sélection des sites selon une valeur seuil. Pour ce faire utiliser le modèle **Identification parcelle**.

Son fonctionnement est le suivant :

- Calcul de la valeur médiane de potentialité par **parcelle**.

- Extraction des **parcelles** suivant une valeur seuil de potentialité.

- Calcul de **centroïdes** à partir desquels sont tracés des chemins de moindre coût jusque le STEU de départ.

<table><thead><tr><td>Entrées</td><td>Description</td></tr></thead><tbody><tr><td>Potentialité</td><td>Raster de potentialité</td></tr><tr><td>Parcellaire</td><td>Description du parcellaire agricole</td></tr><tr><td>Seuil</td><td>Seuil minimal de potentialité à retenir</td></tr><tr><td>Coût</td><td>Raster de coût (créée par modèle)</td></tr><tr><td>Antécédence</td><td>Raster d'antécédence (créé par modèle)</td></tr></tbody></table>

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

<table><thead><tr><td>Sorties</td><td>Description</td></tr></thead><tbody><tr><td>Chemin</td><td>Chemin de moindre coût entre parcelle et STEU (vecteur)</td></tr><tr><td>Parcelles</td><td>Parcelles de potentialité supérieure ou égale au seuil fixé</td></tr></tbody></table>

A noter que si la valeur médiane de potentialité calculée par **parcelle** est systématiquement inférieure au seuil, c'est la valeur maximale qui est retenue.

## AUTRES MODELES

### COMBINAISON PARCELLE/UCS/METEO

Ce modèle a pour objet la préparation du critère déficit à partir de données parcellaires, pédologiques et climatiques. Il procède de la manière suivante :

- Intersection entre les couches de parcellaire, de type de sol et de zonage météo

- Calcul d'un identifiant à partir des identifiants des couches précédentes (contexte agro/pédo-climatique)
- Jointure entre la couche résultante et une table contenant les informations de déficit pour chaque identifiant calculé.

A noter que si l'utilisateur ne dispose pas de zonage météo (type grille) mais d'une station unique, il est possible de remplacer la couche par la zone d'étude dans le paramétrage du modèle (en choisissant un champ numérique adéquat : id ayant comme valeur 1). L'important est que la valeur du champ se retrouve dans l'identifiant calculé, à la fois dans la couche et dans la table (champ de jointure).

<table><thead><tr><td>Entrées</td><td>Description</td></tr></thead><tbody><tr><td>Zone d'étude</td><td>Zone d'étude</td></tr><tr><td>RPG</td><td>Description du parcellaire agricole (RPG ou équivalent)</td></tr><tr><td>UCS</td><td>Délimitation des unités cartographiques de sol (BDSOL/DONESOL)</td></tr><tr><td>Zone météo</td><td>Délimitation du zonage météo</td></tr><tr><td>Table de jointure déficit</td><td>Table de jointure entre parcelle et déficit par l'identifiant grille/UCS/culture</td></tr></tbody></table>

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

<table><thead><tr><td>Sorties</td><td>Description</td></tr></thead><tbody><tr><td>Déficit</td><td>Parcellaire agricole possédant une information de déficit (Critère)</td></tr></tbody></table>

## CRITERES : DISTANCE A UNE ENTITE

Parmi les différents facteurs imaginables certains peuvent être liés à la distance (euclidienne) séparant une parcelle d'une entité (réseau d'irrigation, masse d'eau, zone urbaine, etc). Si l'on se rappelle la section Adaptabilité : Nouveaux critères, sont mentionnés des facteurs sous forme raster, or c'est exactement le cas présent ici. Le modèle **Distance à une entité** permet de créer un raster de distance à l'entité la plus proche. Bien que la fonction existe par défaut dans QGIS elle demande à ce que l'entité soit sous forme raster, le modèle présent fait office de facilitation en :
- Découpant la couche d'entités à la zone d'étude

- Rasterise l'information (5m)

- Effectue le calcul de proximité (dans la même étendue que la zone d'étude)

La possibilité d'exécuter le processus en lot permet un gain de temps si les couches d'entités sont multiples.

<table><thead><tr><td>Entrées</td><td>Description</td></tr></thead><tbody><tr><td>Zone d'étude</td><td>Zone d'étude</td></tr><tr><td>Entité</td><td>Couche d'entités à partir desquelles on souhaite calculer la distance la plus faible en tout point</td></tr></tbody></table>

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

<table><thead><tr><td>Entrées</td><td>Description</td></tr></thead><tbody><tr><td>Distance</td><td>Raster de distance vers l'entité la plus proche</td></tr></tbody></table>

Il suffit ensuite de reclasser le raster (r.recode) à l'aide d'une table de reclassement pour obtenir un critère utilisable en calculatrice raster.

## LIMITES ET CRITIQUES DE L'EXERCICE

Cette méthodologie, bien que pouvant prendre en compte un nombre variable de critères, demeure avant tout relativement théorique. Sa validité reposant sur le choix des critères, leurs classifications et pondérations. Il y a donc un travail en amont et une « calibration » terrain qui peut avoir un impact important sur les résultats offerts.
Le caractère spatial de la méthode rend difficile l'étude de l'impact qualitatif et quantitatif d'un projet de REUT, les données hydrologiques étant partielles (autant sur le plan spatial que temporel).

Egalement les sources de données utilisées sont définies à des échelles assez variables et leurs croisements entrainent une incertitude sur les résultats obtenus. Outre la variabilité spatiale il existe des précautions à prendre sur la validité temporelle des résultats, le parcellaire agricole pouvant fortement varier d'une année à une autre (moins vrai pour la vigne et les vergers).

Enfin il faut bien comprendre qu'il s'agit là d'une étude de potentialité et non de faisabilité, l'identification même d'un site d'intérêt n'assure en rien la réussite d'un projet de REUT, interviennent également (mais pas que) des dimensions sociales, économiques et/ou environnementales.

## ET ENSUITE ?

Avec l'arrivée de QGIS 3 début 2018 et le remplacement progressif de QGIS 2.18 la question de la pérennité de la méthode est soulevée. En effet QGIS 3 marquera une rupture nette par rapport à la version précédence puisque les langages de programmation intervenant changent (passage de Python 2 à Python 3, PyQt 4 à PyQt 5) ce qui rend les modèles et fichiers projet incompatibles d'une version à l'autre. Egalement les extensions disponibles pour la version 2.18 ne seront utilisables dans la version 3 que si un travail de mise à jour a été effectué.

Un travail de réadaptation des modèles est donc à prévoir. En toute logique cela revient à reconstruire le modèle comme une copie de la version 2.18, mais il est possible que certains traitements aient changés (nom/fournisseur de traitement/paramétrage/etc...). Point positif toutefois le modeleur graphique passe d'un code en Python à un code en C++, ce qui implique une accélération très nette de la vitesse d'exécution. Egalement bien que QGIS 3 sorte début 2018 la version 2.18 est prévue pour rester opérationnelle encore pour un certain temps.
## Sites potentiels pour la STEU de Combaillaux
### Données test

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

### Légende

● STEU

■ Territoire CCGPSL/3M

Potentialité

■ Nulle

■ Intérêt faible

■ Intérêt moyen

■ Intérêt élevé

<table><tr><td>Volume total rejeté :</td><td>58400</td></tr><tr><td>Volume stockable (80% ) :</td><td>50656</td></tr><tr><td>Volume estival (Mai-Sept)</td><td>19680</td></tr></table>

<table><thead><tr><th>Potentiel</th><th>Total volume (m3)</th><th>% volume annuel</th><th>%volume stockable</th><th>%volume estival</th></tr></thead><tbody><tr><td>7</td><td>7002</td><td>12%</td><td>14%</td><td>36%</td></tr><tr><td>6 (Moyen)</td><td>2280</td><td>16%</td><td>18%</td><td>47%</td></tr><tr><td>5</td><td>63064</td><td>136%</td><td>157%</td><td>403%</td></tr><tr><td>4</td><td>156843</td><td>432%</td><td>498%</td><td>1283%</td></tr><tr><td>3 (Faible)</td><td>74436</td><td>724%</td><td>834%</td><td>2147%</td></tr><tr><td>2</td><td>256033</td><td>1758%</td><td>2027%</td><td>5217%</td></tr><tr><td>1</td><td>153949</td><td>3341%</td><td>3852%</td><td>9915%</td></tr><tr><td>0 (Nul)</td><td>281504</td><td>6901%</td><td>7956%</td><td>20478%</td></tr></tbody></table>

Attention : Les résultats présentés sur cette carte sont issus de données test à des fins d'illustration et ne peuvent en aucun cas faire l'objet d'une interprétation quelconque.

Auteur : Florent Hubert, IRSTEA Monpellier. Réalisé le : 16/2/2018
Occupation du sol (station de Combaillaux)

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

Légende :

Cultures

Autres boisements

Vignes

Zone urbaine

Vergers

Affleurements

Chemins et sentiers

Masse d'eau

Route empirée

Prairies et pelouses

Route petite

Sables

Route grande

Très petit cours

Chemin de fer

Petit cours

Landes

Cours moyen

Foret ouverte

Neige et glaciers

Peupleraie

Accessibilité

Taillis et futaies

Zone d'accessibilité forte

Auteur : Florent Hubert
UMR G-EAU, IRSTEA Montpellier

Réalisé le : 16/2/2018
Cette section regroupe une description synthétique des différents traitements employés dans le cadre de la méthode.

## DEFINIR UNE PROJECTION

L'outil *Définir la projection courante* permet de spécifier ou remplacer le système de coordonnées d'un jeu de données. Cet outil s'utilise si le système de projection est inconnu ou erroné.

## PROJETER

L'outil *Reprojeter une couche* permet de passer d'un système de coordonnées d'un jeu de données à un autre.

## DECOUPER (VECTEUR)

L'outil *Couper* (Clip en anglais) permet l'extraction de polygone d'une couche suivant le contour d'une couche de découpage :

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

Seule la table attributaire de la couche en entrée est conservée.

## ZONE TAMPON

L'outil *Distance tampon fixe* permet de créer une zone tampon d'une distance spécifiée autour d'entités d'une couche vecteur. Il est possible de fusionner par défaut ces zones :

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

L'outil *Distance tampon variable* permet d'utiliser un champ de la table attributaire comme distance de buffer, permettant de faire varier celle-ci selon le souhait de l'utilisateur.
## STATISTIQUES DE CHAMP

Le script *R Summary Statistics* pouvant être installé sous QGIS permet de calculer l'occurrence de valeurs uniques (un champ) ou combinaisons de valeurs uniques (plusieurs champs) d'une table attributaire. Il est également possible de calculer pour chacun des cas identifiés une somme des attributs d'un ou plusieurs autres champs (dit champ de récapitulation), par exemple la surface par type de parcelle et par commune avec l'occurrence de chacun des cas ainsi que des statistiques pour chacun des champs (min, max, moyenne, somme, écart-type, premier, dernier, intervalle).

Le moyen le plus simple de parvenir au même résultat sans R est l'utilisation d'un tableau croisé dynamique sous Excel.

## UNION

L'outil *Union* calcule l'union géométrique de couches de polygones en conservant la table attributaire de chacune des couches. L'outil *Merge* permet l'union de plus de deux couches en simultané (i.e. $A \cup B$) :

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

## INTERSECTER

L'outil *Intersection* calcule l'intersection de multiples couches vectorielles (i.e. $A \cap B$) :

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

## IDENTITE

L'outil *Polygon Identity* (SAGA) permet de réaliser l'union entre une couche en entrée et le résultat de l'intersection entre cette couche et une couche d'identité (i.e. $A \cup (A \cap B)$ :
<figure><img src="image_97.png" /></figure>

Les entités issues de l'intersection reçoivent les informations attributaires de la couche d'identité.

## JOINTURE

L'outil *Jointure* (Join en anglais) permet de joindre à une couche la table attributaire d'une autre couche ou d'une table partageant un champ en commun (champ de jointure). Le champ ne doit pas avoir nécessairement le même nom ni le même nombre d'entrées.

## METTRE A JOUR

L'outil *Polygon update* (SAGA) permet la superposition géométrique d'une couche en entrée par une couche de mise à jour. L'ordre des couches importe donc. Si les deux couches possèdent le même champ dans leurs tables attributaires respectives, ces champs seront fusionnés. (i.e. $(A\B)\cup B)$ :

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

## ELIMINER

L'outil *Supprimer les débordements* permet de supprimer des entités sélectionnées d'une couche en les fusionnant aux entités adjacentes. L'entité à supprimer est fusionnée à l'entité adjacente de plus grande taille ou de plus grande surface adjacente.

<figure><img src="image_99.png" /></figure>
Cet outil permet de supprimer les micro-polygones ou artéfacts issus de croisements successifs.

## FUSIONNER

L'outil *Dissolve* permet de fusionner les entités adjacentes d'une même couche. Il est possible de définir des conditions de fusion selon les attributs des entités :

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

## POLYGONE VERS RASTER

L'outil *Rastérisation* permet de convertir une couche de polygones vers une couche raster de taille de cellule spécifiée. Les valeurs du raster sont celles du champ de l'entité vecteur spécifiée.

## POLYGONE VERS LIGNE

L'outil *De polygones à lignes ou v.to.lines* (GRASS) permet de créer une couche de lignes correspondant aux contours des polygones d'une couche en entrée.

## ENTITES VERS POINT

L'outil *Centroïdes de polygones* permet de créer une couche de point correspondant aux barycentres des polygones d'une couche en entrée.

## RASTER VERS POLYGONES

L'outil *Polygoniser* (GDAL/OGR) permet de créer une couche de polygones à partir d'un raster, chaque polygone étant délimité par un ensemble de cellules de même valeur.

## RASTER VERS POLYLIGNES

L'outil *v.to.vect* (GRASS) permet de créer une couche de polylignes à partir d'un raster, chaque polyligne étant délimitée par un ensemble de cellules de même valeur.

## CREER UN QUADRILLAGE

L'outil *Grille vecteur* permet la création d'un quadrillage régulier sur une étendue définie (manuellement ou à partir d'une autre couche). Le quadrillage peut être composé de lignes ou de polygones. Cet outil est utile pour découper un ensemble en zones régulières ou pour créer des points disposés régulièrement.

## DECOUPER (RASTER)
L'outil *Découper un raster selon une couche de masque (GDAL/ODR)* permet l'extraction de données d'un jeu raster selon une étendue spécifiée (manuellement ou par une entité vecteur). Si l'étendue est spécifiée par une entité vecteur il est possible de l'utiliser comme contour pour le raster sortant (tout pixel n'étant pas en intersection avec le vecteur ne sera pas inclus dans la couche en sortie) :

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

## RECLASSIFICATION

L'outil *Reclassify value (SAGA)* ou *r.recode (GRASS)* permet de modifier ou classer les valeurs d'un raster selon une table de classification à rentrer ou à charger.

## MOSAÏQUE VERS NOUVEAU RASTER

L'outil *Mosaic raster layer (SAGA)* permet de combiner plusieurs jeux de données raster vers un nouveau selon l'ordre de superposition.

## CALCULATRICE RASTER

La *Calculatrice raster* est un outil permettant l'écriture et l'exécution d'expression pour la création et modification de raster. Il existe différentes calculatrice sous QGIS (GRASS, SAGA, GDAL/OGR)

## STATISTIQUE FOCALE

L'outil *r.neighbors (GRASS)* permet de calculer pour chaque cellule une statistique de voisinage (min, max, moyenne, diversité, écart-type, etc). La taille et forme (rectangle, cercle) du voisinage sont modifiables.

## DISTANCE DE CHEMIN

L'outil *r.walk.rast* créé un raster pour lequel chaque cellule présente la valeur de coût cumulé minimal pour rejoindre la cellule depuis une source. Il dépend donc d'un raster de coût et d'un raster de surface (MNT). L'outil créé aussi un raster d'antécédence indiquant la direction vers la source.

## CHEMIN DE COUT

L'outil *r.drain* permet de créer un raster correspondant au chemin de moindre coût entre une source et une destination. Ce chemin est créé à partir des rasters créés par l'outil *r.walk.rast*.

## SUPPRIMER UN CHAMP
L'outil *Refactoriser les champs* permet la suppression rapide d'un ou plusieurs champs de la table attributaire d'une couche en entrée.

## PROCHE

L'outil *v.distance* (GRASS) identifie pour chaque entité d'une couche vecteur en entrée l'entité la plus proche d'une autre couche vecteur et calcule la distance séparant les deux (élément d'une couche A le plus proche d'un autre d'une couche B).

## MULTI-PARTIES VERS UNE PARTIE

L'outil *De morceaux multiples vers morceaux uniques* permet de créer une couche d'entités à partir d'une autre en séparant les entités multi-parties.

## PENTE

L'outil *Pente* permet de calculer une pente en degré ou pourcent à partir d'un modèle numérique de terrain (altitude).
Registre Parcellaire Graphique (RPG)

En 2000, la réglementation européenne a institué dans tous les Etats Membres, de localiser et d'identifier les parcelles agricoles sur une base graphique réputée plus simple et mieux à même de prendre en compte les enjeux territoriaux et environnementaux de l'agriculture.

Les îlots et leur occupation culturale sont déclarés annuellement par les exploitants agricoles pour bénéficier des aides PAC. Un îlot est un ensemble de parcelles culturales :

-contiguës, entières ou partielles, portant une ou plusieurs cultures, exploitées par le déclarant ;

- limité par des éléments facilement repérables et permanents, comme des voies de communications, les cours d'eau ou d'autres exploitations ;

- stable d'une année sur l'autre

Description

Format d'objet :

Shapefile (polygones)

Objets :

Ilots de parcelles

Attributs :

Groupement de cultures exploité, surface, irrigation selon le niveau de détail (4 niveaux d'information)

Intérêts pour le diagnostic régional

Origine de l'information géographique

Source :

Déclaration PAC sur orthophotographies

Mode de production :

Saisie à 1 :5 000 initiale effectuée par l'exploitant sur fonds A3.

L'ASP traite ensuite les données afin d'en supprimer les données attributaires et nominatives des îlots. Enfin un identifiant numérique non significatif est généré par îlot afin de faire le lien avec les tables attributaires

Qualité des données graphiques

Echelle :

Départementale

Précision / Qualité :

1 :5 000

Territoire disponible :

France entière

Système de projection :

Lambert 93

Validité des données

Date de validité :

2014

Délai de mise à jour :

Annuel

Accès et diffusion

Producteur :

ASP (Agence de Services et de Paiement s) et Ministère de l'Alimentation, de l'Agriculture et de la Pêche.

Modalité de diffusion :

Convention de mise à disposition entre la CA34 et IRSTEA

Modalité d'utilisation :

Information réutilisable conformément aux clauses de la convention signée

Format des fichiers :

.shp
# Référentiel Régional Pédologique RRP (BDSOL L-R)

Le programme IGCS (Inventaire, Gestion et Conservation des Sols) intègre des opérations d'inventaire des sols réalisées à différentes échelles sur un même territoire. Le produit élaboré, appelé "Référentiel Régional Pédologique", est une base de données géographique régionale, documentée et mise à jour.

Il est structuré pour permettre une manipulation aisée des informations par des logiciels d'analyses cartographiques (système d'informations géographique), et son intégration dans des domaines d'application variés.

Sa précision correspond au minimum à celle d'une représentation cartographique au 1/250 000. Certains petits secteurs de références, choisis pour la représentativité de leurs sols sont eux cartographiés au 1 : 10 000. Ils servent de sites expérimentaux pour les domaines agricoles et environnementaux

<table><thead><tr><th></th><th>Description</th></tr></thead><tbody><tr><td>Format d'objet :<br>Objets :<br>Attributs :</td><td>Vecteur et base de données<br>Unité Cartographique de Sols (UCS), unités typologiques de sols (UTS), profils, horizons.<br>En principe doivent respecter les spécifications de la base de données DONESOL établie par l'INRA : Code et libellé d'UCS, nombre code libellé et pourcentage des UTS (Unité Typologique de Sol)/UCS, contraintes/UTS*, régime hydrique/UTS*, régime submersion/UTS*, excès d'eau/UTS*, discontinuité (type, profondeur d'apparition, conséquence), texture (pierrosité, argile, limons fins et totaux, sable, taux MO : note de 1 à 10) /UTS*, pH/UTS*<br>*pas forcément renseigné pour toutes les UTS d'une même UCS</td></tr><tr><td colspan="2">Intérêts pour le diagnostic régional</td></tr><tr><td colspan="2">Informations pédologiques géolocalisées les plus complètes à ce jour pour cette échelle</td></tr><tr><td colspan="2">Origine de l'information géographique</td></tr><tr><td>Source :<br>Mode de production :</td><td>INRA Orléans + Partenaires techniques et financiers régionaux<br>Numérisation de cartes levées au 1 : 250 000</td></tr><tr><td colspan="2">Qualité des données graphiques</td></tr><tr><td>Echelle :</td><td>1 : 250 000</td></tr><tr><td>Précision / Qualité :</td><td>30 m</td></tr><tr><td>Territoire disponible :</td><td>38 référentiels régionaux disponibles. RRP du Languedoc-Roussillon utilisé dans le cadre de cette étude</td></tr><tr><td>Système de projection :</td><td>Lambert II étendu</td></tr><tr><td colspan="2">Accès et diffusion</td></tr><tr><td>Producteur :</td><td>Unité INFOSOL de l'INRA Orléans, diffusion par SIG LR</td></tr><tr><td>Modalité de diffusion :</td><td>Adhésion IRSTEA à l'association SIG L-R</td></tr><tr><td>Format des fichiers :</td><td>Shapefile, BDD Access</td></tr></tbody></table>
SAFRAN

Le Système d'Analyse Fournissant des Renseignements Adaptés à la Nivologie est un système d'analyse conçu pour l'analyse prévisionnelle des risques d'avalanches et les impacts du changement climatique. Sa couverture nationale et la fréquence d'actualisation horaire en font une source de données climatiques de grande valeur pour les études à petite et moyenne échelle.

Description

Format d'objet :

Grilles de 8 x 8 km au pas de temps journalier

Objets :

Points

Attributs :

Données évapotranspiration, température et pluie

Intérêts pour le diagnostic régional

Origine de l'information géographique

Source :

Météo France

Mode de production :

Le modèle Safran s'appuie à la fois sur l'ensemble des observations disponibles du réseau d'observations (précipitation, température...) ainsi que sur les analyses des modèles atmosphériques (Arpège ou CEPMMT) qui servent d'ébauche pour l'analyse.

Qualité des données graphiques

Echelle :

--

Précision / Qualité :

8x8 km

Territoire disponible :

France entière

Système de projection :

Lambert II étendu

Validité des données

Date de validité :

--

Délai de mise à jour :

--

Accès et diffusion

Producteur :

Météo- France (Centre National de Recherches
Météorologiques, CNRM)

Modalité de diffusion :

IRSTEA

Modalité d'utilisation :

Extraction via R ou Excel

Format des fichiers :

Binaire ou txt
# Données Agence de l'Eau Rhône-Méditerranée et Corse

Afin de répondre aux exigences réglementaires de diffusion des données environnementales, l'agence de l'Eau Rhône-Méditerranée et Corse met à disposition sur son site internet, des données brutes ou valorisées obtenues par les réseaux de mesures de la qualité de l'eau qu'elle pilote et les études qu'elle mène

Elle propose notamment des données sur l'état des lieux des masses d'eau et les objectifs du SDAGE 2021.

## Description

Format d'objet :

Tables

Objets :

Etat des masses d'eau, RNABE, RNAOE

Attributs :

Nom, état chimique, état biologique

## Intérêts pour le diagnostic régional

Ces données permettent d'attribuer à chaque exploitation, le nombre de forage quelle possède ainsi que les volumes quelle prélève.

## Origine de l'information géographique

Source :

Agence de l'Eau Rhône-Méditerranée et Corse

Mode de production :

## Qualité des données graphiques

Echelle :

--

Précision / Qualité :

--

Territoire disponible :

Bassins versants Rhône-Méditerranée et Corse

Système de projection :

Lambert 93

## Validité des données

Date de validité :

--

Délai de mise à jour :

--

## Accès et diffusion

Producteur :

Agence de l'eau Rhône-Méditerranée et Corse

Modalité de diffusion :

https://www.eaurmc.fr/

Modalité d'utilisation :

La réutilisation des informations publiques suppose que celles-ci ne soient pas altérées, que leur sens ne soit pas dénaturé et que leurs sources et la date de leur dernière mise à jour soient mentionnées.

Format des fichiers :

Excel
BD LISA

La BD LISA constitue le référentiel géographique national qui découpe le territoire national en entités hydrogéologiques (formations géologiques aquifères ou non). Une entité hydrogéologique est une partie de l'espace géologique :

* délimitée à une certaine échelle géographique : le "niveau",

* rattachée à un type de formation géologique : le "thème",

* définie par ses potentialités aquifères : la "nature",

* caractérisée par un type de porosité (qui permet de distinguer les principaux modes de circulation de l'eau) : le "milieu",

* caractérisée par la présence ou non d'une nappe, qui peut être libre et/ou captive : l'"état".

Description

Format d'objet :

Shapefile (polygones)

Objets :

Entités hydrogéologiques

Attributs :

Dénomination, niveau, état, milieu, ordre, nature, thème,

Intérêts pour le diagnostic régional

Localisation des formations karstiques

Origine de l'information géographique

Source :

Sandre (diffusion)

Mode de production :

Forages, modélisations, études

Qualité des données graphiques

Echelle :

Du 1 :1 000 000ème au 1 :50 000

Précision / Qualité :

--

Territoire disponible :

France

Système de projection :

Lambert 93

Validité des données

Date de validité :

--

Délai de mise à jour :

--

Accès et diffusion

Producteur :

BRGM

Modalité de diffusion :

Open Data

Modalité d'utilisation :

Open Data

Format des fichiers :

shp
BD TOPO®

La BD TOPO® contient une description des éléments du paysage sous forme de vecteurs de précision métrique, classés selon une thématique adaptée.

Elle couvre de manière cohérente l'ensemble des entités géographiques et administratives du territoire national. Elle permet la visualisation, le positionnement, la simulation au service de l'analyse et de la gestion opérationnelle du territoire.

Description

Format d'objet :

Vectoriel (points, lignes, polygones)

Objets :

Bâtiments, limites administratives, végétation, toponymes, orographie, points d'intérêts

Attributs :

Nature, hauteur, altitude, nom, etc.

Intérêts pour le diagnostic régional

Ces données sont utiles pour la construction d'une occupation du sol et l'identification de structures clés

Origine de l'information géographique

Source :

Institut Géographique National (IGN-F)

Mode de production :

Photogrammétrie, GPS, orthophotographie, autres sources

Qualité des données graphiques

Echelle :

1: 5000 à 1 : 50 000

Précision / Qualité :

Métrique

Territoire disponible :

France entière

Système de projection :

Lambert 93

Validité des données

Date de validité :

2016

Délai de mise à jour :

2 ans

Accès et diffusion

Producteur :

IGN

Modalité de diffusion :

Convention IRSTEA

Modalité d'utilisation :

Base de données soumise aux conditions générales d'utilisation des fichiers IGN

Format des fichiers :

shp
BD FORET V1®

La BD Forêt® est une base de données vecteur de référence pour l'espace forestier et les milieux semi-naturels. Elaborée par photo-interprétation d'images en infrarouge couleurs de la BD ORTHO®, la BD Forêt® est réalisée par emprises départementales sur le territoire métropolitain.

Elle attribue à chaque plage cartographiée de plus de 5000m² un type de formation végétale. Ses principales caractéristiques sont les suivantes :

- une nomenclature nationale de 32 postes qui repose sur une décomposition hiérarchique des critères, distinguant par exemple les peuplements purs des principales essences forestières de la forêt françaises

- un type de formation végétale attribué à chaque plage cartographiée supérieure ou égale à 50 ares (5 000 m²)

- une couche géométriquement compatible avec le RGE® et donc en parfaite cohérence avec la couche végétation de la BDTOPO®

Description

Format d'objet :

Vectoriel (polygones)

Objets :

Contours forestiers et des espaces semi-naturels

Attributs :

Type de végétation (nomenclature en 32 postes)

Intérêts pour le diagnostic régional

Ces données sont utiles pour la construction d'une occupation du sol.

Origine de l'information géographique

Source :

Institut Géographique National (IGN-F)

Mode de production :

Photo-interprétation

Qualité des données graphiques

Echelle :

---

Précision / Qualité :

10m

Territoire disponible :

France entière (2014)

Système de projection :

Lambert 93

Validité des données

Date de validité :

---

Délai de mise à jour :

---

Accès et diffusion

Producteur :

IGN

Modalité de diffusion :

Convention IRSTEA

Modalité d'utilisation :

Base de données soumise aux conditions générales d'utilisation des fichiers IGN

Format des fichiers :

shp
RGE ALTI® 5M

Le RGE ALTI® décrit la forme et l'altitude normale de la surface du sol à grande échelle.

Description

Format d'objet :

Raster

Objets :

Modèle Numérique de Terrain

Attributs :

Altitude, précision

Intérêts pour le diagnostic régional

Fournit la donnée d'altitude permettant les calculs de pentes.

Origine de l'information géographique

Source :

Institut Géographique National (IGN-F)

Mode de production :

Données LIDAR, images aériennes

Qualité des données graphiques

Echelle :

1 : 3 000 et le 1 : 20 000

Précision / Qualité :

5m

Territoire disponible :

France entière (2016)

Système de projection :

Lambert 93

Validité des données

Date de validité :

---

Délai de mise à jour :

---

Accès et diffusion

Producteur :

IGN

Modalité de diffusion :

Convention IRSTEA

Modalité d'utilisation :

Base de données soumise aux conditions générales d'utilisation des fichiers IGN

Format des fichiers :

shp
# BD CARTHAGE®

La BD CARTHAGE® est la base de données complète du réseau hydrographique français.

## Description

Format d'objet :

Vecteur (ligne, polygone)

Objets :

Entités hydrographiques, découpage administratif en quatre niveaux du territoire

Attributs :

Nom, code, attributs différents selon l'objet

## Intérêts pour le diagnostic régional

Fournit la donnée d'altitude permettant les calculs de pentes.

## Origine de l'information géographique

Source :

SANDRE

Mode de production :

Elle est produite par les Agences de l'eau à partir de la base de données cartographiques BD CARTO® de l'IGN qui décrit sous forme numérique l'ensemble des informations présentes sur le territoire national.

## Qualité des données graphiques

Echelle :

1 : 5 000 et le 1 : 50 000

Précision / Qualité :

Métrique à décimétrique

Territoire disponible :

France entière

Système de projection :

Lambert 93

## Validité des données

Date de validité :

2016

Délai de mise à jour :

Annuel

## Accès et diffusion

Producteur :

IGN

Modalité de diffusion :

Licence ouverte Etalab

Modalité d'utilisation :

Licence ouverte Etalab

Format des fichiers :

shp
# OSO CESBIO

La BD CARTHAGE® est la base de données complète du réseau hydrographique français.

## Description

Format d'objet :

Raster

Objets :

Occupation du sol, statistiques de confiance et validité

Attributs :

-

## Intérêts pour le diagnostic régional

Fournit les données de bases pour l'obtention d'un raster d'occupation du sol satisfaisant

## Origine de l'information géographique

Source :

CESBIO

Mode de production :

Produit issu d'algorithmes de télédétection, développé par le CESBIO depuis 2013 à partir de données Sentinel-2.

## Qualité des données graphiques

Echelle :

1 : 5000

Précision / Qualité :

20m

Territoire disponible :

France entière

Système de projection :

Lambert 93

## Validité des données

Date de validité :

2016

Délai de mise à jour :

Annuel

## Accès et diffusion

Producteur :

CESBIO

Modalité de diffusion :

Open Data Commons Attribution Licence, visualisation et téléchargement en ligne.

Modalité d'utilisation :

Open Data Commons Attribution Licence

Format des fichiers :

GeoTiff<
<figure><img src="image_102.png" /></figure>

Étude réalisée avec le soutien financier
De l'Agence de l'Eau Rhône-Méditerranée-Corse