Ces travaux de these s'integrent dans les objectifs du projet international d'Analyse de la Mousson Africaine (AMMA). Ils abordent la caracterisation des echanges entre l'ocean et l'atmosphere dans l'Atlantique Equatorial Est (AEE) de l'echelle diurne a interannuelle a partir de donnees d'observation et de la modelisation numerique. Ces echanges sont d'une tr'es grande importance pour le declenchement de la Mousson Africaine qui determine en grande partie la saison des pluies en Afrique de l'Ouest. Nous montrons que les parametres de la couche de melange oceanique dans le GG sont principalement pilotes, a l'echelle diurne, par les flux de chaleur et les processus de subsurface (l'entraınement, le melange vertical) et nous montrons aussi qu'a partir d'une parametrisation de la turbulence oceanique (TKE) fermee a l'ordre 1.5, on parvient a restituer la turbulence oceanique observee dans cette region. Afin de quantifier le role de la tension du vent sur la variabilite des parametres de la couche de melange oceanique, nous utilisons un modele lineaire force par les anomalies de la tension du vent. Les resultats montrent que pres de l'equateur (3°S-3°N), le signal est retrouve ; par contre de part et d'autre de la bande 3°S-3°N, les amplitudes des parametres de la couche de melange sont tres largement sous estimes, ce qui montre que la tension du vent n'est pas le seul parametre a moduler la variabilite de la temperature surface 'a ces latitudes. Les bilans de chaleur integres sur la couche de melange ont permis de montrer le role des flux de chaleur et du melange vertical dans l'AEE. En effet, l'AEE est subdivise en plusieurs boıtes en tenant compte de la dynamique et de la thermodynamique et dans chaque boıte le bilan y est calcule. Les resultats montrent que les flux de chaleur et le m'elange vertical dominent le bilan aux echelles saisonniere et interannuelle. Ce melange vertical, calcule par residu, est compare avec les donnees independantes de turbulence oceanique (Dengler et al., 2010 ; Rhein et al., 2010) collectees lors des campagnes EGEE/AMMA en 2005-2007. Les resultats montrent que le residu se compare tres bien avec ces donnees independantes en terme de variabilite spatiale et temporelle. Le melange vertical est tres fort dans la region sauf dans les boıtes au sud de l'AEE et sa variabilite est largement pilotee par les flux de flottabilite et la tension du vent. En regard de ces resultats avec ceux obtenus sur le cycle diurne, on peut dire que dans l'AEE, la variabilite spatio-temporelle des parametres de la couche de melange est principalement pilotee a toutes les echelles par les flux de chaleur et les processus de subsurface.

En Amazonie orientale, l'agriculture familiale contribue significativement à la transformation des écosystèmes forestiers en pâturages. La mise en valeur pastorale se fait généralement de manière monospécifique avec une Poacée pérenne d'origine africaine (Brachiaria brizantha cv. Marandu) qui constitue la principale offre fourragère. Mais les pratiques pastorales mises en oeuvre sur une exploitation conduisent souvent à une dégradation des pâturages, qui se manifeste par une prolifération des plantes adventices. Le pâturage a été modifié par l'introduction de deux Fabacées pérennes fourragères originaires d'Amérique Latine, l'une herbacée (Arachis pintoï cv. Amarelo) et l'autre ligneuse (Leucaena leucocephala). Pendant la saison sèche, nous avons étudié les interactions entre B. brizantha et des Fabacées fourragères implantées. Cette étude est la première approche intégrée (écophysiologie) qui étudie les caractéristiques adaptatives et les effets des Fabacées sur la capacité de résistance à la contrainte hydrique des pâturages à partir des évolutions des échanges gazeux foliaires, du potentiel hydrique et de l'état hydrique du sol. Le dispositif expérimental était composé de cinq parcelles de 100 m2 isolées du bétail correspondant à cinq traitements différents. Les trois espèces étudiées ont été cultivées seules et en association. Trois répétitions ont été réalisées pour chaque traitement afin de valider statistiquement les résultats et prendre en compte la variabilité spatiale du sol. En monoculture, les valeurs de conductances stomatiques de B. brizantha sont relativement élevées au regard des conditions climatiques par rapport à des valeurs habituelles de plantes C4 au champ alors que les valeurs d'assimilations nettes sont celles mesurées couramment. A. pintoï et L. leucocephala possèdent des valeurs de conductances stomatiques et d'assimilations mesurées habituellement sur les plantes C3 au champ. En situation de sécheresse, les trois espèces étudiées en monoculture adoptent un mouvement de fermeture stomatique 30 jours après l'arrêt des pluies et ajustent ainsi leur conductance sur l'épuisement de la réserve utile du sol. Cette stratégie permet le maintien de l'activité photosynthétique indispensable à la survie cellulaire et de maintenir une transpiration suffisante pour réguler leur température. Les trois espèces survivent par évitement de la sécheresse, en réduisant leur surface foliaire active et en fermant leurs stomates dès l'abaissement du potentiel hydrique. B. brizantha a le même comportement lorsqu'il est cultivé seul ou associé à A. pintoï que ce soit au niveau de la conductance stomatique, de l'assimilation de CO2 ou de l'efficience instantanée de l'eau. Cependant, son potentiel hydrique est affecté plus précocement au cours d'un stress hydrique lorsqu'il est en concurrence avec A. pintoï. Dans le cas de cette association, notre étude montre que la présence de B. brizantha a un effet négatif sur les activités photosynthétiques et donc sur la production de biomasse d'A. pintoï, espèce considérée pourtant comme sciaphile. Enfin, les fonctions métaboliques de B. brizantha sont réduites lorsqu'il est associé à L. leucocephala du fait de l'ombrage. L. leucocephala évite la sécheresse en réduisant sa surface foliaire. L'augmentation d'insolation en fin de saison sèche affaiblit cependant les plantes associées aux strates inférieures et réduit la quantité de fourrage disponible sur pied. Nous proposons aux agriculteurs de mettre en place une gestion durable de leurs pâturages par la création d'associations végétales fourragères adaptées aux contraintes biotiques et abiotiques. Les résultats de notre étude nous permettent de proposer plusieurs possibilités de gestion des ressources fourragères. Chacune des espèces fourragères montre des aptitudes face aux contraintes environnementales et leur association peut-être profitable si l'implantation et la pression pastorale sont raisonnées selon leurs caractéristiques culturales. Une plus grande diversité de l'offre fourragère permettra de stabiliser ces éleveurs sur leurs terres, de ralentir la déforestation et ainsi conserver la biodiversité existante

Dans ce travail, des données pétro-géochimiques, isotopiques et géochronologiques ont été couplées afin de caractériser et dater les circulations de fluides dans le Cisaillement Sud Armoricain. Ces circulations ont été datées par les méthodes U-Pb sur zircon et monazite et 40Ar/39Ar sur muscovite, dans les mylonites et les granites syntectoniques présentant des traces d'altération hydrothermale. A 315 Ma, le Cisaillement Sud Armoricain localise la circulation de fluides magmato-métamorphiques, responsables de la cristallisation de veines de quartz géantes qui lui sont parallèles. Vers 300 Ma, un afflux de fluides météoriques est enregistré dans les mylonites, certains échantillons granitiques et quelques veines de quartz, obliques au Cisaillement Sud Armoricain, interprétées comme des fentes de tension d'échelle crustale. Les résultats obtenus ont permis d'étudier le comportement des chronomètres isotopiques au cours des interactions fluide-roche. Si le zircon apparaît comme relativement robuste, les circulations de fluides provoquent la recristallisation de la monazite et de la muscovite, et ainsi perturbent, voir remettent à zéro les chronomètres isotopiques U-Pb et KAr. Par exemple, dans le massif syntectonique de Questembert, les dates 40Ar/39Ar obtenues sur six échantillons s'étalent sur 16 Ma, du fait de l'infiltration pervasive de fluides. Ces résultats illustrent donc la nécessité d'associer la géochronologie à une étude pétrogéochimique détaillée, depuis la caractérisation des roches jusqu'à l'étude texturale des monazites ou la documentation des zonations chimiques des muscovites, par exemple.

Dans un contexte où des millions de personnes dépendent de la ressource en eau exposée au caprice de la mousson en Afrique de l'Ouest, ce travail de thèse vise à mieux appréhender les processus hydrologiques en zone de socle, et notamment à évaluer le rôle de la redistribution latérale de l'eau dans le sol par les interactions entre les réservoirs souterrains, la végétation et l'atmosphère, par la caractérisation et la modélisation à l'échelle de la parcelle et du versant. Ce travail s'appuie sur le dispositif expérimental déployé dans le petit bassin versant de l'Ara dans le cadre de l'observatoire AMMA-CATCH. La mise en oeuvre du modèle ParFlow-CLM permet de simuler les transferts dans la zone saturée et la zone non saturée par la résolution de l'équation de Richards en 3D, en étant conditionné par un forçage atmosphérique en surface. Après avoir identifié les paramètres influents pour les transferts verticaux, une caractérisation spatiale de ces paramètres a été menée. La configuration du modèle a ensuite été évaluée en 1D. Il est montré que le modèle reproduit de manière pertinente les séries temporelles du bilan d'énergie et la distribution de l'eau dans le sol. L'effet de la variabilité spatiale des paramètres hydrodynamiques est ensuite étudié à l'échelle de la parcelle. Enfin, en incluant les géométries de socle connues et une distribution de végétation, les transferts horizontaux souterrains sont mis en évidence avec la formation de zones sèches ou humides en relation avec des distributions spatiales d'évapotranspiration.

Le réchauffement climatique est lié aux augmentations des concentrations de gaz à effet de serre dans l'atmosphère terrestre et en particulier aux émissions anthropiques de CO2. La séquestration géologique a la capacité et la longévité potentielles pour diminuer de façon significative les émissions anthropiques de CO2. Cette séquestration à grande profondeur induit des risques de fuite des réservoirs géologiques. Parmi les scénarios de fuite envisagés, celui d'une fuite diffuse est le plus inquiétant puisque sans surveillance, cette fuite pourrait perdurer et entrainer des séquelles sur l'environnement ainsi que des risques pour les populations. Des outils et protocoles de surveillance doivent donc être mis au point pour la surveillance en proche surface. Ce travail de thèse s'inscrit dans le cadre de cette problématique. Il a pour objectif la caractérisation, la quantification et la modélisation du transport et des interactions du CO2 dans une zone non saturée carbonatée. Ce travail a suivi une approche expérimentale sur un site pilote naturel à Saint-Emilion (Gironde, France), avec la réalisation de fuites diffuses en ZNS carbonatée. Cette étude aborde plusieurs thématiques: la description et l'instrumentation du site pilote naturel ; la caractérisation physico-chimique de l'hétérogénéité du réservoir carbonaté ; l'étude du fonctionnement naturel de la ZNS carbonatée et en particulier la mise en place d'une ligne de base des concentrations en CO2 ; la caractérisation de l'extension des panaches de gaz suite à des expériences de fuite diffuse dans la ZNS carbonatée et l'étude par simulation numérique des interactions gaz-eau-roche lors d'une fuite diffuse de CO2 dans une ZNS carbonatée. Les résultats de ces travaux montrent l'importance de la caractérisation de l'hétérogénéité du réservoir carbonaté ainsi que des techniques d'échantillonnage et d'analyse des différentes phases en présence. L'établissement de la ligne de base a une importance particulière pour permettre de distinguer les variations naturelles de celles induites par une fuite diffuse de CO2 dans la ZNS carbonatée. Les modes de transport du CO2 vont évoluer en fonction des paramètres physico-chimiques. Ce transport se fait par advection et/ou par diffusion. L'utilisation de gaz inertes au niveau du site de séquestration géologique est très importante puisque la détection de ces traceurs permettrait de prédire les arrivées de panaches de CO2 en proche surface. Par ailleurs, les interactions chimiques doivent être prises en compte dans les modèles de transport afin de pouvoir définir les facteurs de retard et l'impact d'une fuite diffuse de CO2 sur une ZNS carbonatée.

Cette recherche porte sur les interactions entre la tectonique, le climat et les processus géomorphologiques de la surface de la Terre à travers l'étude des terrasses fluviales au Venezuela et l'Albanie. Ces deux domaines ont été soumis à une surrection modérée, aux variations climatiques du Quaternaire et fournissent un large registre de terrasses fluviales. Ces contextes donnent l'opportunité d'étudier la dynamique de la formation des terrasses fluviales à l'échelle de temps 102-105 ans. Ainsi, une approche morpho-chronologique a été appliquée afin de parvenir à une meilleure compréhension de cette dynamique dans les rivières du Venezuela et de l'Albanie.Dans le système des rivières Pueblo Llano et Santo Domingo situé sur le flanc sud-est des Andes de Mérida au Venezuela, douze terrasses fluviales ont été identifiées pour les derniers 200 ka. Les mesures des concentrations du 10Be réalisées sur les dépôts fluviatiles et glaciaires ont permis d'estimer l'âge d'exposition de sept terrasses et d'un complexe de moraines frontales. Un modèle dynamique de terrasses soutenues par ces datations et des données géomorphologiques, stratigraphiques et sédimentologiques, indique que la formation de terrasses est principalement contrôlée par des variations climatiques de haute fréquence (103 à 104 ans) qui induisent un déséquilibre entre capacité de transport de l'eau et apport des sédiments. Néanmoins, le type de réponse a été très lié à l'altitude du site et l'influence des glaciers en amont. En effet, dans la partie supérieure du système, la succession des phases d'aggradation et d'incision a été synchronisée avec la succession des périodes froides-sèches et chaudes-humides, tandis que la partie inférieure du système montre une tendance opposée. Basé sur la restauration temporelle du taux d'incision de la partie inférieure du système, un taux de surrection de 1.1 mm/an au cours des dernières 70 ka a été estimé pour le flanc sud-est des Andes de Mérida. En outre, l'identification et la datation d'un complexe de moraines frontales situé à l'altitude de 2300 m au dessus du niveau de la mer met en évidence le fait que l'avancée glaciaire au cours du Dernier Maximum Glaciaire, dans d'autres domaines des Andes de Mérida, pourrait également avoir atteint des élévations inférieures à celles déjà signalés entre 2900 et 3500 m au dessus du niveau de la mer.En Albanie, les six principales rivières ont été analysées. De nouvelles données géomorphologiques et géochronologiques ont été intégrées aux données publiées afin de proposer un cadre stratigraphique et chronologique homogène à l'échelle régional pour les derniers 200 ka. Sur la base de ce cadre, l'âge de la formation de onze terrasses fluviales régionales a été établi. En Albanie, les processus de formation des terrasses ont également été principalement contrôlés par des variations climatiques de haute fréquence (103 à 104 ans). Néanmoins, les résultats montrent également que la réponse géomorphologique du système fluvial était probablement modulée par la taille des bassins versants et par les variations eustatiques. En effet, pour la période précédant le MIS 2, les terrasses de remplissage et d'abrasion ont été formées respectivement dans les grands et petits bassins versants. Pour la période après le début du MIS 2, une relation complexe entre les variations climatiques et eustatiques a favorisé le développement de terrasses d'abrasion partout en Albanie. Malgré des différences entre les réponses des rivières, la succession de phases d'aggradation ou d'abrasion latérale et de phases d'incision a été synchronisée avec la succession des périodes froides-sèches et chaudes-humides dans toutes les rivières de l'Albanie.

La densification actuelle de l'urbanisation conduit à un changement d'occupation du sol et du réseau de drainage en zone péri-urbaine. L'écoulement est concentré dans des fossés ou des réseaux d'assainissement et ainsi accéléré. Ces phénomènes peuvent avoir des conséquences importantes pour les nappes et les cours d'eau, comme par exemple l'aggravation des crues et des sécheresses et l'altération de la qualité chimique et/ou biologique du milieu. Sous la pression réglementaire (Directive Cadre Européenne sur l'Eau) et la demande sociale, les gestionnaires sont donc confrontés à des choix complexes en terme d'aménagement. Ainsi, il est nécessaire de mettre au point des méthodes et des modèles capables de quantifier l'impact de l'augmentation de l'urbanisation sur la vulnérabilité des hydro-systèmes péri-urbains. Dans ce travail, nous formulons l'hypothèse qu'une modélisation hydrologique spatialisée continue, prenant explicitement en compte les objets des paysages périurbains (parcelles urbaines, agricoles, forestières, haies,..) et les éléments déterminant les chemins de l'eau (topographie, mais aussi réseaux d'assainissement ou de routes) peut aider à comprendre et hiérarchiser le rôle des différents objets du paysage sur la réponse hydrologique. Pour ce faire, le modèle distribué PUMMA (Peri-Urban Model for landscape Management) adapté aux bassins versants péri-urbains, a été développé dans le cadre de cette thèse. Il consiste en une intégration du modèle existant URBS, décrivant des processus hydrologiques urbains à l'échelle d'une parcelle cadastrale, dans la plate-forme de modélisation LIQUID, qui contient déjà des modules représentants des processus hydrologiques en zone rurale. PUMMA a également été complété par de nouveaux modules simulant les déversoirs d'orage, les bassins de rétention et le transfert du ruissellement de surface. Le modèle suit une approche orientée objet dans laquelle le paysage est divisé en mailles irrégualières, correspondant aux parcelles cadastrales en zone urbaine et aux unités de réponse hydrologiques (HRUs) en zone rurale. Afin de pouvoir appliquer le modèle à l'échelle d'un bassin versant, des méthodes automatiques pour la préparation des données géographiques ont été mises au point. De plus, une méthode a été développée pour déterminer les contours de bassins versants péri-urbains, en distinguant les surfaces contributives de temps sec et humide. Le modèle a été appliqué au bassin versant de la Chaudanne, un sous-bassin de l'Yzeron, situé en zone péri-urbaine lyonnaise pour deux années en conditions contrastées (sèche et humide). Les paramètres du modèle ont été spécifiés à partir des observations disponibles et des données de la littérature. Les résultats montrent un comportement du modèle réaliste et une aptitude à représenter les comportements différents en période sèche et humide, en lien avec des degrés de saturation des sols différents, même si les pics de débits d'été sont en général surestimés. Différents tests de sensibilité sur certains processus/paramètres montrent l'importance des processus urbains sur la réponse hydrologique du bassin, comme en particulier la génération de ruissellement de surface par les surfaces imperméables et naturelles urbaines, le drainage de l'eau du sol par les réseaux d'assainissement et les connexions entre les îlots urbains et le réseau hydrographique naturel et artificiel. L'épaisseur des sols et la conductivité hydraulique à saturation latérale jouent aussi un rôle important sur la dynamique du débit de base. Nous montrons aussi le potentiel du modèle pour tester différents scenarii d'aménagement ou de gestion des eaux pluviales.

Il est aujourd'hui avéré que la composition chimique de l'atmosphère des régions enneigées - et notamment des régions polaires - est sensiblement affectée par les échanges d'espèces chimiques réactives entre l'air et la neige. En effet, le manteau neigeux constitue un véritable réacteur photochimique multiphasique, mais les mécanismes physico-chimiques à l'œuvre en son sein sont encore mal connus. Une compréhension détaillée des processus s'y déroulant est indispensable pour modéliser correctement la composition et la réactivité de l'atmosphère au-dessus des régions enneigées. De plus, la reconstitution de l'évolution post-dépôt des composés chimiques stables de la neige est également un préalable indispensable pour permettre l'interprétation paléoclimatique de leurs profils de concentration enregistrés dans les carottes de glace.Le nitrate (NO3-) présent dans la neige joue un rôle fondamental, car sa photolyse induit notamment l'émission d'oxydes d'azote (NOx = NO + NO2) par le manteau neigeux, qui modifient la capacité oxydante de l'atmosphère via la production d'ozone. L'objet de cette thèse a donc été d'étudier par modélisations les processus physico-chimiques intervenants dans l'évolution de la concentration du nitrate dans la neige.Une première approche, prolongeant des études préexistantes, a visé à identifier un mécanisme réactionnel pour la photochimie du nitrate dans la neige, en postulant notamment l'existence d'une couche quasi-liquide à la surface des grains de neige. Néanmoins, les propriétés exactes de l'interface air - glace sont, à l'heure actuelle, encore mal caractérisées, et il est apparu que cette démarche présentait de trop larges incertitudes pour être poursuivie.Une discussion approfondie a alors été menée afin d'évaluer les tentatives actuelles de modélisation de la chimie de la neige, et dans le but de proposer une nouvelle approche plus réaliste au regard du niveau de connaissance actuel.Ainsi, dans une seconde partie, l'ensemble des processus d'échange physico-chimiques du nitrate entre l'air et la neige ont été étudiés puis modélisés : adsorption à la surface, diffusion en phase solide et co-condensation. Parmi les résultats obtenus, il est apparu que les paramétrisations actuelles de la couverture surfacique en nitrate étaient incapables de reproduire les concentrations mesurées, dans le cas de la neige de surface à Dome C, et révèlent d'importantes surestimations. A contratio, la prise en compte conjointe de la diffusion en phase solide ainsi que d'un processus de co-condensation permet de bien reproduire qualitativement les séries temporelles de plus d'un an, couvrant donc à la fois l'été et l'hiver austral, qui présentent chacun des caractéristiques distinctes en terme de concentration mesurées.Cette étude révèle ainsi l'importance de ces processus physico-chimiques d'échange dans la modélisation de la chimie de la neige, et pose les bases des mécanismes à prendre en compte dans le cadre de développements futurs.

Ce document constitue une synthèse de ma dizaine d'années d'activité scientifique passées sous les tropiques, en tant que physicien adjoint affecté à l'IPSL et détaché à la Réunion, au LPA qui devint LACY en s'associant à Météo-France. Ce document contient donc une présentation de l'essor instrumental de la station de l'OPAR et sa place dans les réseaux de recherche internationaux, qui reflète mon activité dans le cadre de ma mission de participation aux services d'observation. Une synthèse des principaux résultats obtenus par rapport au bilan de l'ozone troposphérique tropical est ensuite présentée. En effet, cette dizaine d'années d'observation de l'ozone troposphérique par lidar et radiosondages nous a permis d'approfondir des cas d'études pour détailler les mécanismes source de l'ozone troposphérique à la Réunion, d'établir la climatologie de l'ozone et de la tropopause et de dégager les premières estimations de tendances de l'ozone troposphérique.

L'importance climatique des hautes latitudes est exacerbée par le contexte actuel du réchauffement climatique, de part de leur forte sensibilité à ces changements et en raison des rétroactions globales majeures qu'elles sont susceptibles d'engendrer. La modélisation offre un moyen d'estimer ces impacts dans les temps passés, présents et futurs, tout en quantifiant les incertitudes procédant des imperfections de notre connaissance de ces environnements et de leur représentation. L'amélioration et l'évaluation de la représentation des hautes latitudes dans les modèles de climat globaux répondent donc à de forts enjeux scientifiques et sociétaux : c'est dans ce cadre précis que s'inscrit mon travail de thèse. Le gel des sols est une spécificité majeure des régions circum-arctiques, porteuse d'implications climatiques aux plans thermiques, hydrologiques et biogéochimiques. Une paramétrisation des impacts hydrologiques du gel des sols a été introduite dans le schéma hydrologique multi-couches du modèle de surfaces continentales ORCHIDEE : ses effets sur le contenu en eau des sols et le régime hydrologique des principaux bassins de l'Arctique ont été évalués par comparaison à des données de terrain, révélant la plus-value d'une telle représentation mais aussi les lacunes résiduelles de la modélisation, qui touchent à l'absence de représentation des réservoirs temporaires d'eau de surface et, sans doute, d'une paramétrisation sous-maille du gel des sols. Parallèlement, une représentation des effets thermiques du gel des sols développée pour un modèle antérieur à ORCHIDEE a été révisée et évaluée à différentes échelles spatiales par comparaison à des données observationnelles : si la représentation de l'énergie de chaleur latente augmente la température des sols soumis au gel saisonnier, un biais froid subsiste dans la modélisation, imputable à une représentation imparfaite de la neige. Une étude de sensibilité conduite sur cette variable en confirme les implications thermiques mais aussi biogéochimiques à l'échelle des régions circum-arctiques, sous-tendues par les importantes quantités de matière organique que ces régions renferment. Alors que les caractéristiques de la neige sont souvent représentées comme spatialement uniformes dans les modèles de climat globaux, la simple prise en compte du caractère particulièrement isolant de la neige de taïga engendre des changements importants dans le cycle du carbone aux hautes latitudes, et souligne les incertitudes entachant notre représentation actuelle de ces écosystèmes. Les propriétés thermiques de la neige n'en sont pas l'unique vecteur, mais une évaluation détaillée de notre modélisation sur un site de permafrost arctique (station de Bayelva, Svalbard) désigne la neige comme une source majeure des incertitudes associées à notre modélisation des hautes latitudes, au travers de représentations inadaptées de son albédo, sa rugosité de surface, son contenu variable en eau liquide pouvant accommoder de l'eau de pluie. En termes hydrologiques, l'absence de représentation spécifique des zones de montagne, des caractéristiques hydrauliques des sols à granulométrie grossière du Haut-Arctique, et des nombreuses étendues d'eau libre des régions circum-arctiques, limite notre capacité à représenter raisonnablement des principales caractéristiques de l'hydrologie de surface de ces régions. Le diagnostique de ces limites définit autant de potentiels d'amélioration de la modélisation des hautes latitudes, sources possibles de développements futurs.