modelisation_hydrologique_corse_pro:start []

Cette page est en lecture seule. Vous pouvez afficher le texte source, mais ne pourrez pas le modifier. Contactez votre administrateur si vous pensez qu'il s'agit d'une erreur.
====== Modélisation hydrologique et le Model builder ======

<WRAP center round box 60%>
Objectif de la séance : reconstruire le réseau hydrographique de la Corse actuel, lors du dernier maximum glaciaire et lors du bas niveau marin messinien.
</WRAP>


<WRAP center round download 60%>
**Données nécessaires à la séance**

{{ :modelisation_hydrologique_corse_pro:modelisationhydrologique.zip |}}

{{ :modelisation_hydrologique_corse_pro:modelisationhydrologique.pdf |La modélisation hydrologique}}


</WRAP>


=====Partie 1. Visualisation du Modèle Numérique de Terrain=====

  * Chargez le MNT de la Corse à la résolution de 100m.
  * Découpez le raster en deux parties avec la fonction "**SetNull**" : la partie sous-marine -> __corse100mer__ et la partie aérienne -> __corse100cont__.

{{:modelisation_hydrologique_corse_pro:capture_d_ecran_2024-10-15_090635.png?300|}}

  * Adaptez la symbologie des deux parties.
  * Créez un ombrage.
 
=====Partie 2. Calcul du réseau hydrologique =====

Vous allez commencer la modélisation hydrologique sur la partie terrestre.
Les outils de modélisation hydrologique sont disponibles dans l'Arctoolbox sous "**Outils Spatial Analyst - Hydrologie**".


====1. Remplissage des cuvette==
Afin de créer des réseaux hydrographiques continus vous devez supprimer les cuvettes du raster.
A partir du raster __corse100cont__ :
  * Supprimer les cuvettes avec la fonction "**Remplissage**" -> __corse100fill__.
  * Vous pouvez essayer de visualiser les cuvettes remplies en utilisant la "**calculatrice raster**". Il suffit de soustraire __corse100cont__ à __corse100fill__.

====2. Création des écoulements====

  * Sur le raster __corse100fill__ déterminez le sens d'écoulement avec la fonction "**Direction de flux**" -> __corse100dir__ 
  * Sur le raster __corse100dir__ calculez les écoulements avec la fonction "**Accumulation de flux**" -> __corse100ecoul__.

====3. Classification des écoulements====
  * Sur le raster __corse100ecoul__ simplifiez et classifiez le réseau hydrologique avec la fonction "**Ordre d'écoulement - Méthode STRAHLER**" -> __corse100class__.
  * Adaptez la symbologie du raster __corse100mclass__ pour que le réseau hydrologique ne soit ni trop simple ni trop complexe à l'affichage.

====4. Transformer le réseau hydrographique en fichier de formes====
  * Transformez le raster __corse100class__ en polylignes avec la fonction "**Ecoulement vers entité**" -> __corse100ecoulements__. Adaptez la symbologie pour que le réseau hydrologique ne soit ni trop simple ni trop complexe à l'affichage.
  
====5. Création des bassins versants====  
  * Déterminez les bassins versants avec la fonction "**Bassin**" -> __corse100bass__.
  * Transformez les Bassins en polygones avec la fonction "**Raster vers polygones**" -> __corse100bassinsversants__.


=====Partie 3. Création d'une chaîne de traitement avec le Model builder=====

Créez votre propre chaîne de traitement avec le model builder. Votre chaîne doit contenir la suite des traitements ci-dessous. En **paramètre d'entrée** du modèle un raster (ici corse100mrgf) et en **paramètres de sortie** le réseau hydrographique en polylignes, les bassins versants en polygones et l'accumulation de flux. 

===Liste des traitements===

  - Remplissage
  - Direction de flux
  - Accumulation de flux
  - Ordre d'écoulement
  - Ecoulement vers entité
  - Bassin
  - Raster vers polygones

{{:modelisation_hydrologique_corse_pro:capture_d_ecran_2024-10-14_190125.png?700|}}

Quand votre chaîne de traitement fonctionne, vous pouvez la faire fonctionner sur le raster initial contenant la partie sous-marine : __corse100mrgf__.

=====Partie 4. Paramètres des bassins versants=====

====1. Connaitre la surface et l'altitude max des bassins versants====

L'outil "**Statistiques zonales (table)**" permet de récupérer les statistiques d'un raster (altitude max, min, ...) située en dessous de polygones et de les stocker dans une table. Utilisez-le pour extraire les caractéristiques des bassins versants. L'option "Couche de jointure en sortie" permet d'intégrer directement sous forme de jointure le résultat dans le fichier de formes des bassins versants.

{{:modelisation_hydrologique_corse_pro:capture_d_ecran_2024-10-15_092542.png?350|}}

===Nettoyage du fichier de forme des bassins versants===

  * Ce fichier possède toujours une jointure, cassez cette jointure en l'exportant vers un autre fichier de forme -> __corse100bassinsversantsstat__.
  * Avec l'outil "Selectionner Lasso" supprimer tous les bassins versants hors de la Corse. Vous pouvez selectionner les BV corses et inverser la sélection avant de supprimer.
  * Dans la table attributaire la surface est calculée en m², ce qui n'est pas pratique. Recalculez la géométrie du champ AREA en km².
  * Supprimez tous les petits bassins versants dont la surface est inférieure à 1 km².
  * Enregistrer les modifications


Pour analyser les bassins versants, vous pouvez créer un graph : surface vs altitude max des BV.

{{:modelisation_hydrologique_corse_pro:capture_d_ecran_2024-10-14_180156.png?700|}}


====2. Connaitre les parties du littoral corse les plus influencées par les apports fluviatiles====

L'objectif est maintenant de connaitre les sections du littoral soumises aux apports des cours d'eaux.
Il faut commencer par localiser l'ensemble des exutoires sur le littoral. Pour cela nous allons utiliser l'outil "**Intersection**" entre le réseau hydro et le trait de côte.

{{:modelisation_hydrologique_corse_pro:capture_d_ecran_2024-10-15_095815.png?300|}}

Le résultat de ce traitement est un fichier de forme de type "multi-points" ce qui est génant pour la suite du traitement, transformez-le en points simples avec le traitement "**Entités vers points**" -> __ExutoiresCorses__.

Nous devons maintenant intégrer la surface du bassin versant au niveau des exutoires (l'idéal aurait été d'avoir directement les débits, mais nous ne les avons pas pour l'ensemble des cours d'eaux corses).
Utilisez l'outil "**Extraction de valeurs multiples vers points**" afin d'intégrer les valeurs de l'accumulation de flux au niveau des exutoires.

{{:modelisation_hydrologique_corse_pro:capture_d_ecran_2024-10-15_100328.png?300|}}

Dans la table attributaire, rajouter un champ "SurfaceBVkm2" de type double.
Avec la "**Calculatrice de champs**" calculer la surface des bassins versants à partir des valeurs d'écoulements. Pour rappel cette valeur correspond au nombre de cellules contenues dans le BV. Si votre raster à une résolution de 100m vous devez donc multiplier cette valeur par 100² pour l'avoir en m² et donc le multiplier par 0.01 pour l'avoir en km².

Le traitement "Densité de noyaux" va vous permettre de cartographier la densité des exutoires sur le littoral mais en prenant en compte la surface des bassins versants.

{{:modelisation_hydrologique_corse_pro:capture_d_ecran_2024-10-15_101958.png?300|}}

{{:modelisation_hydrologique_corse_pro:capture_d_ecran_2024-10-14_183815.png?600|}}


=====Partie 5. Comparaison avec le réseau Sandre=====

Service d'administration nationale des données et référentiels sur l'eau ([[https://www.sandre.eaufrance.fr/v2/|Sandre]]) fourni une cartographie du réseau hydrographique national.

Les cours d'eaux et bassins versants corses sont disponibles : {{ :modelisation_hydrologique_corse_pro:hydrosandrecorse.zip|ici}}. Comparez les données diffusées par le réseau Sandre avec le résultat de votre modélisation hydrologique. 
=====Partie 6. PaléoHydrologie de la Corse au LGM et lors de la crise messinienne=====

Pour la suite de la séance nous allons reconstruire les réseaux hydrographiques de la Corse lors du dernier niveau marin (-120 m lors du Dernier Maximum Glaciaire) et lors de la crise de salinité messinienne (-1500 m il y a 5,5 millions d'années).

  * Découper le raster 100mCorse à -120 m avec la fonction "**SetNull**"
  * Faites fonctionner votre chaîne sur la partie continentale supérieure à -120 m

Faites la même opération pour la basse de niveau marin à -1500 m du Messinien.


===== Corrections =====

Actuel

{{:modelisation_hydrologique_corse_pro:modern.png?600|}}

Dernier maximum glaciaire

{{:modelisation_hydrologique_corse_pro:lgm.png?600|}}

Messinien

{{:modelisation_hydrologique_corse_pro:messinien.png?600|}}



===== Divers =====

<WRAP center round download 100%>
{{:modelisation_hydrologique:geomatique2012.pdf|Modélisation hydrologique des bassins versants du Sahel}}

[[https://www.nature.com/articles/s41598-022-12055-1|Freshwater influx to the Eastern Mediterranean Sea from the melting of the Fennoscandian ice sheet during the last deglaciation]]


</WRAP>