Le coe ur du mémoire concerne l'enchaînement intégré de la modélisation stochastique des pluies et de la modélisation des débits et autres flux associés. <br />Une première partie traite de l'aspect spatial des transferts, visant à exprimer sous une forme très synthétique la structure arborescente du réseau hydrographique. Elle suffit à transcrire des effets d'échelle et les grands contrastes internes au bassin versant concernant la fonction de production, ce qui asseoit une méthode de régionalisation efficace. Une variante plus réaliste dans la représentation de l'hydraulique permet de représenter l'atténuation hydraulique au sein d'un bassin versant. <br />La modélisation stochastique de la pluie (CECP) procède par désagrégation de la pluie totale d'un épisode en sous épisodes. Cette décomposition \textit{quantité} x \textit{forme} permet d'utiliser des fonctions monotones à forme donnée, et donc d'inverser le modèle ou de calculer facilement des probabilités conditionnelles, opérations généralement impossibles avec d'autres modèles. Le modèle n'est donc pas restreint à une utilisation en simulation, mais permet un calcul explicite (numérique) des fréquences rares. <br />Les probabilités conditionnelles se révèlent par ailleurs une représentation conjointe des aspects dynamiques et probabilistes que l'on propose d'appeler modèle pluie-débit-durée-fréquence, et qui manifeste l'existence de durées critiques dépendantes de la fréquence. On propose par ailleurs une formulation approchée, sous forme de loi dérivée, c'est à dire une distribution de probabilité dont les paramètres sont les paramètres de la transformation pluie-debit. <br />Pour estimer ce modèle dans des conditions variées de disponibilité de l'information, on propose d'utiliser l'information de Kullback-Leibler, avec une adaptation pour les valeurs extrêmes.

Les bassins versants influencés par la présence de barrages-réservoirs constituent les objets d'étude de cette thèse. L'objectif principal est de mettre au point une méthode de prise en compte, au sein d'un modèle pluie-débit global, de l'impact de retenues d'eau artificielles sur le comportement hydrologique du bassin versant. Un tel outil, simple dans sa conception, doit permettre d'améliorer la simulation des débits de bassins versants influencés et d'évaluer l'impact de barrages-réservoirs sur le régime hydrologique des cours d'eau, en particulier lorsque l'on se place loin à l'aval de l'ouvrage. La démarche adoptée est empirique dans le sens où l'on a cherché, sur une base de données de 46 bassins versants situés en France, aux Etats-Unis et au Brésil, à exploiter au mieux, dans la structure du modèle pluie-débit, une information simple sur les ouvrages, à savoir l'évolution temporelle des volumes stockés dans les retenues. La diversité des conditions hydro-climatiques présentes dans la base de données donne aux résultats obtenus un caractère assez général. Au cours de cette thèse, un grand nombre de solutions de prise en compte des barrages-réservoirs a été testé. Les résultats montrent que cette information sur les remplissages des ouvrages peut être utilement exploitée au sein d'une structure de modèle global pour améliorer la simulation des bassins influencés. Les solutions qui se sont révélées être les plus performantes n'engendrent pas de complexification de la structure initiale du modèle puisqu'aucun ajout de paramètre n'a été nécessaire. Des gains de performances significatifs ont été réalisés, en particulier sur la simulation des étiages. Il s'agit d'un résultat important de la thèse car il montre qu'il est possible de prendre explicitement en compte, dans un modèle pluie-débit global, les barrages-réservoirs présents sur le bassin versant. Nous avons ensuite essayé de comprendre le lien entre la méthode de prise en compte des barrages-réservoirs et certaines caractéristiques physiques des bassins versants et des ouvrages de stockage afin de tenter d'adapter la solution au bassin versant considéré. Les résultats de ces tests ont montré que ce lien est difficile à mettre en évidence et qu'à ce stade, une solution générique semble la plus efficace. La solution la plus performante a ensuite été évaluée sur un échantillon complémentaire de 31 bassins versants n'ayant pas servi à son développement. L'évaluation a confirmé l'intérêt de la méthode pour l'amélioration des simulations des débits mais a également montré sa sensibilité à la disponibilité des données de stockages ainsi qu'aux transferts interbassins inconnus.

La modélisation pluie-débit est une discipline des sciences hydrologiques qui s'intéresse principalement au cycle hydrologique afin d'aboutir à un outil de travail exploitable dans différents domaines. Deux dimensions définissent la modélisation pluie-débit. La première est l'échelle d'espace. Pour différentes superficies, allant d'une maille de quelques mètres carrés au bassin versant de quelques dizaines de milliers de kilomètres carrés, il existe une littérature très riche qui va de la proposition de divers types de modèles jusqu'à la critique de ceux ci. La deuxième dimension est l'échelle de temps. Du pas de temps horaire jusqu'au pas de temps pluriannuel, on trouve dans la littérature de nombreux modèles fonctionnant à chaque pas de temps. Cependant, rares sont les auteurs qui se sont intéressés à la modélisation avec une vue d'ensemble de tous les principaux pas de temps. La présente thèse tente de répondre à une question qui nous semble intéressante : peut-on tirer profit de l'échelle de temps dans la modélisation pluie-débit en vue d'obtenir des modèles plus efficaces et plus cohérents entre eux?<br />Nous nous sommes intéressés dans ce travail de recherche, aux pas de temps pluriannuel, annuel et mensuel. Pour chaque échelle de temps, nous avons procédé à une étude comparative des modèles tirés de la littérature en les testant sur un échantillon de 429 bassins versants de caractéristiques hydro-climatiques très variables. Grâce à des allers-retours entre les pas de temps, nous avons largement amélioré les modèles disponibles pour chaque pas de temps. Ainsi nous avons pu proposer trois modèles conceptuels globaux aux pas de temps pluriannuel, annuel et mensuel. De plus, suite à une vue d'ensemble de ces trois pas de temps, en y ajoutant le pas de temps journalier, nous avons décelé une cohérence dans la présence d'un échange souterrain indiquant clairement qu'un bassin versant est un système qui est ouvert sur toutes les frontières et pas seulement sur celle avec l'atmosphère.

L'objectif de notre recherche est d'améliorer les performances des modèles pluie-débit par une meilleure prise en compte de l'ETP au sein de ces modèles. Pour cela, l'échantillon de bassins versants étudiés devait représenter diverses conditions climatiques et hydrologiques. L'assemblage des différentes données collectées (pluie, débit et paramètres climatiques) a permis l'élaboration d'un échantillon de 308 bassins versants dont 221 français, 79 nord-américains et 8 australiens. L'avantage de travailler sur un échantillon assez vaste est de tirer des conclusions relativement libres de toute dépendance vis à vis des caractéristiques particulières de tel ou tel bassin. De plus, quatre modèles pluie-débit ont été utilisés: le modèle GR4J et des versions globales des modèles HBV, IHACRES et TOPMODEL. Tout d'abord, une étude de sensibilité de ces modèles pluie-débit à l'entrée d'ETP a été entreprise. Deux aspects ont été traités en particulier. La première motivation était de tester l'impact d'une connaissance fine du climat par rapport à une connaissance des normales (moyennes interannuelles). Les résultats ont confirmé le manque d'utilisation judicieuse de l'information datée contenue dans les formules d'ETP par les modèles pluie-débit. Le deuxième aspect était de tester plusieurs formules pour représenter les variations de l'ETP au sein des modèles pluie-débit. Les résultats ont montré là aussi une faible sensibilité des modèles quant au choix fait pour la formule d'ETP. Nous avons porté en particulier notre attention sur la formule de Penman, couramment utilisée en modélisation, car elle est la plus satisfaisante d'un point de vue physique, à l'échelle de la parcelle. En modélisation pluie-débit, elle est dépassée par d'autres formules plus simples (utilisant moins de données climatiques). Ainsi, nous pouvons parler de la pertinence d'autres formules que Penman pour la modélisation: des formules simples faisant intervenir uniquement la température de l'air sont aussi performantes que la formule de Penman. D'un point de vue opérationnel, ces résultats sont extrêmement rassurants puisqu'ils suggèrent qu'une entrée d'ETP simple à obtenir peut être utilisée à la place de formules d'ETP plus demandeuses en données. Cependant, d'un point de vue scientifique, le constat d'apparente insensibilité des modèles pluie-débit à l'entrée d'ETP est préoccupant. En effet, le modèle semble se contenter d'une représentation extrêmement simplifiée de la demande évaporatoire. Pourtant, il est légitime de penser que le modèle serait plus performant s'il prenait mieux en compte cette information supplémentaire. Les hypothèses avancées pour expliquer la faible sensibilité des modèles à l'entrée d'ETP sont multiples : - S'agit-il d'un phénomène naturel que le modèle retranscrit assez fidèlement ? - Faut-il remettre en cause la structure du modèle, notamment la partie traitant l'évaporation ? - Les formules d'ETP traditionnelles sont-elles inadaptées à la modélisation pluie-débit ? Alors que la première interrogation nous semblait difficile à investiguer, les deux autres étaient plus à la portée d'un utilisateur de modèles pluie-débit conceptuels. Pour cela, nous avons utilisé une approche par données synthétiques d'ETP qui permet notamment de lever toute incertitude sur la pertinence des données d'ETP. Les résultats ont montré que les différences entre ETP datées et interannuelles étaient presque totalement absorbées par les réservoirs de production des modèles pluie-débit. Le modèle se comporte donc comme un filtre passe-bas pour les données atmosphériques, atténuant ainsi la variabilité temporelle des entrées climatiques.

Les sécheresses exceptionnelles observées ces dernières années, ont montré la vulnérabilité de certaines activités humaines face à une pénurie de la ressource en eau. Le retour d'expérience de la sécheresse de 2003 a montré que cette situation aurait pu être mieux anticipée, incitant les services qui ont à charge la gestion des cours d'eau à se doter d'outils de prévision des étiages. Cette recherche porte sur l'analyse des étiages des cours d'eau de l'est de la France situés dans les bassins du Rhin et de la Meuse. L'objectif est de proposer une modélisation conceptuelle applicable pour la prévision opérationnelle. Un des intérêts réside dans le choix d'un nombre élevé de bassins, situés dans un contexte climatique homogène, mais aux caractéristiques géologiques diversifiées. Cet échantillonnage permet d'aborder une étude comparée du comportement des bassins et d'en expliquer la variabilité spatiale des écoulements en étiage. Dans un souci de transposabilité nous avons cherché à orienter notre démarche vers des méthodes uniques, applicables et reproductibles pour l'ensemble des bassins. Au final, malgré la diversité des milieux géologiques des bassins étudiés, des modèles uniques, aux paramètres descripteurs d'une réalité de terrain sont adoptés. Enfin, les incertitudes liées à la démarche conceptuelle ont pu être réduites grâce à l'analyse statistique, en proposant un compromis entre les fondements conceptuels et les applications employant les approches probabilistes.

Deux volets distincts de l'hydrologie sont abordés, prévision et prédétermination, au travers d'une problématique commune : le transfert à l'exutoire des bassins versants d' une information hydrologique distribuée. Dans le domaine de la prévision des crues, la technologie radar procure une information pluviométrique spatialement continue. Les hydrologues disposent ainsi en temps réel de la connaissance des champs de pluie, atout indéniable pour l'anticipation des crues notamment sur des petits bassins versants par le biais de la modélisation de la pluie en débit. Dans le cadre de la méthode AIGA d'alerte crues, développée au Cemagref, une modélisation mise en oeuvre à l'échelle du pixel de pluie fournit une cartographie des contributions de débit des pixels. Dans le domaine de la prédétermination, le Cemagref a développé la méthode SHYREG qui associe un modèle régionalisé de simulation de pluies horaires à une modélisation de la pluie en débit. Une estimation statistique régionale des pluies et des débits spécifiques de différentes durées, dans une large plage de fréquence (du courant à l'exceptionnel) peut ainsi être proposée et cartographiée. L'objectif du travail présenté est d'étudier et d'élaborer des méthodologies simples de transfert de ces deux informations débitmétriques discrétisées information temps réel pour le volet prévision et information statistique pour le volet prédétermination - à l'information débit à l'exutoire du bassin versant. La méthodologie met en oeuvre des informations spatiales et une modélisation de la pluie en débit. Pour répondre à l'objectif fixé, trois axes de travail sont développés. Le premier est l'étude du comportement d'un modèle pluie-débit simple développé pour être mis en oeuvre à la maille du km². On examine en particulier s'il satisfait les caractéristiques d'invariance et de parcimonie souhaitée pour une utilisation à la fois en reconstitution de crues et en simulation. Le second axe de travail concerne l'agglomération, en prédétermination, de l'information débit statistique connue au km² pour l'estimation des quantiles de débit à l' exutoire de bassins versants de superficie plus importante dans le cadre de la méthode SHYREG. Il s'agit de tenir compte de deux phénomènes hydrologiques distincts : l'abattement spatial de la pluie et le transfert dans le réseau hydraulique. Le troisième axe de travail concerne l'agglomération de l'information hydrologique distribuée pour la reconstitution des crues dans le cadre de l'outil AIGA d'alerte crues. Différentes modélisations sont proposées pour transférer à l'exutoire les contributions des débits modélisées aux pixels. / We address the routing of distributed hydrological information to the outlet of watersheds, in the fields of flood forecasting and flood prediction on ungauged watersheds in the French Mediterranean area.Flood forecasting can benefit of areal rainfall data provided in real-time by radar networks. This data used as an input to rainfall runoff models gives access to flood anticipation on small ungauged watersheds. Within the framework of the AIGA method, developed by CEMAGREF to provide floods alert, a rainfall-runoff model is implemented at the spatial resolution of the radar data, thus providing a map of the 1 km² pixel contributions to the runoff at the catchment outlet.Flood prediction consists of assessing the frequency of occurrence of floods of different given magnitude without reference to the times at which they would occur. The SHYREG flood prediction method, developed by Cemagref associates a regionalized rainfall model with a rainfall-runoff model. It provides grids of statistical estimates of rain and runoff for various duration and return periods. Our purpose is to study and work out simple methodologies to aggregate these two gridded hydrological data - real time information for the AIGA forecasting method and statistical data for the SHYREG prediction method to the catchments outlets. Our methodology implements distributed information and a rainfall-runoff model. We have first studied the behaviour of a simple rainfall-runoff model developed to be implemented in a gridded resolution (1 km² cells) for prediction as well as for forecasting purposes. We have checked that the model parameters show no redundancy and no link with the characteristics of the rainfall events. We have then addressed the question of the aggregation of gridded hydrological data. Within the SHYREG method, it consists of assessing statistical flow estimates at catchments outlets, knowing simulated flow distributions in each cell of the catchments. This aggregation would combine two distinct hydrological phenomena: areal reduction of rainfall and discharge attenuation in the channel network. Within the AIGA method, we have focused on the routing function of the rainfall-runoff model at the 1 km² cell scale, this scale being the first step of the runoff routing from the production area to the outlet of the catchment. We have then produced streamflow hindcasts for selected observed events using different routing function, within our rainfall-runoff model.

Cette thèse a pour but de contribuer à l’étude des aquifères karstiques en (1) identifiant les facteurs influençants le stockage et l’écoulement de l’eau souterraine et en (2) localisant et estimant de manière quantitative les ressources en eau au sein d’une série carbonatée. Cette étude a été réalisée sur l’aquifère de Dardennes, près de Toulon, en France. L’eau de l’aquifère est utilisée pour l’alimentation en eau potable et alimente le fleuve Las, susceptible de provoquer de fortes crues. Le bassin versant de l’aquifère est estimé à 70 km². Des mesures de porosité et de géochimie isotopique sur les roches carbonatées sont mises en relation dans le but de localiser et d’expliquer les vides créés au sein de la série sédimentaire. Des coupes équilibrées et des analyses structurales ont été réalisées pour évaluer la géométrie de l’aquifère. L’acquisition de données en continu (CTD) et de données hydrochimiques permet d’étudier le dynamisme et le fonctionnement de l'aquifère. Une approche de modélisation pluie-débit avec la plateforme de modélisation KarstMod, basée sur les chroniques de précipitations et de débits sur trois années a été utilisée afin d’estimer la part de l’écoulement disponible pour l’eau potable et celle générant des crues éclair. Une déconvolution du débit des sources a été effectuée à l’échelle annuelle et à l’échelle des crues. Les résultats du modèle ont été comparés avec la gestion actuelle de l’aquifère, montrant une exploitation actuelle équilibrée, et ouvrant des perspectives pour la mise en place d’une gestion active, combinant une exploitation plus importante de l’eau, un écrêtage des crues et un respect de l’écologie du cours d’eau en aval. / The aim of this PhD is to contribute to the karst aquifers understanding by (1) identifying the factors influencing the groundwater storage and flow and by (2) locating and quantitatively estimating water resources within the carbonate sedimentary series. This study was conducted on the Dardennes aquifer near Toulon, France. The aquifer water is used for the supply of drinking water and feeds the Las River, which can cause important floods. The recharge area of the aquifer is estimated at 70 km². Measurements of porosity and isotopic geochemistry on carbonate rocks are related in order to locate and explain the voids created within the sedimentary series. Balanced cross-sections and structural analyses were performed to evaluate the geometry of the aquifer at depth. Continuous data acquisition (specific conductivity, temperature and water depth-pressure) and hydrochemical data (major ions, stables isotopes of water and carbon-13) allow to study the aquifer dynamism and functioning. A daily rainfall-discharge modeling approach was used with the KarstMod platform to estimate the flow part available for drinking water and the part that generates karst flash-floods. A hydrograph separation of springs discharge was carried out at annual and floods scales. Model results have been compared with the current groundwater management, showing a balanced current exploitation. In new perspective, it would be interesting to establish active management, combining a greater groundwater exploitation, floods mitigation and stream ecology respect.

Les modèles pluie-débit sont des outils essentiels pour de nombreuses applications hydrologiques, notamment la prévision des crues. L’objet de cette thèse est d’examiner les performances d’un modèle événementiel distribué, dont l’intérêt est de résumer la représentation des processus à la phase de crue, et la condition initiale à un indice de saturation du bassin facilement observable ou accessible. Ce dernier dispense de modéliser la phase inter-crue, et simplifie la paramétrisation et le calage du modèle. Le modèle étudié combine une fonction de production type SCS et une fonction de transfert type lag and route, appliquées à une discrétisation du bassin en mailles carrées régulières.Le modèle est d’abord testé sur le bassin versant du Real Collobrier. Ce bassin méditerranéen est suivi depuis plus de 50 ans par l’IRSTEA, et dispose d’une exceptionnelle densité de mesures de pluies et de débits. Cet environnement favorable permet de limiter les incertitudes sur l’estimation des pluies et d’évaluer les performances intrinsèques du modèle. Dans ces conditions, les crues sont bien reconstituées à l’aide d’un jeu de paramètres unique pour l’ensemble des épisodes testés, à l’exception de la condition initiale du modèle. Celle-ci apparaît fortement corrélée avec l’humidité du sol en début d’épisode, et peut être prédéterminée de façon satisfaisante par le débit de base ou l’indice w2 fourni par le modèle SIM de Météo-France. Les performances du modèle sont ensuite étudiées en dégradant la densité des pluviomètres, et rendent compte du niveau de performances du modèle dans les cas que l’on rencontre le plus souvent. .La variabilité spatiale des paramètres du modèle est étudiée à l’échelle de différents sous-bassins du Real Collobrier. La comparaison a permis de mettre en évidence et de corriger un effet d’échelle concernant l’un des paramètres de la fonction de transfert. Les relations entre la condition initiale du modèle et les indicateurs d’humidités des sols en début d’épisode restent bonnes à l’échelle des sous-bassins, mais peuvent être significativement différentes selon les sous-bassins. Une seule relation ne permet pas de normaliser l’initialisation du modèle sur l’ensemble des sous-bassins, à une échelle spatiale de quelques km2 ou dizaines de km2. Dans le cas de l’indice d’humidité du sol w2, une explication possible est que cet indice ne prend pas en compte suffisamment finement les propriétés des sols. Enfin, la variabilité spatiale des paramètres du modèle est étudiée à l’échelle d’un échantillon d’une quinzaine de bassins méditerranéens de quelques centaines de km2, associés à des paysages et des fonctionnements hydrologiques divers. La comparaison montre qu’à cette échelle, le lien entre les indicateurs de saturation du bassin et la condition initiale peut rester stable par type de bassin, mais varie significativement d’un type de bassin à l’autre. Des pistes sont proposées pour expliquer cette variation.En conclusion, ce modèle événementiel distribué représente un excellent compromis entre performances et facilité de mise en œuvre. Les performances sont satisfaisantes pour un bassin donné ou pour un type de bassin donné. L’analyse et l’interprétation de la variabilité spatiale des paramètres du modèle apparaît cependant complexe, et doit faire l’objet du test d’autres indicateurs de saturation des bassins, par exemple mesures in situ ou mesures satellitaires de l’humidité des sols. / Rainfall-runoff models are essential tools for many hydrological applications, including flood forecasting. The purpose of this thesis was to examine the performances of a distributed event model for reproducing the Mediterranean floods. This model reduces the parametrization of the processes to the flood period, and estimates the saturation of the catchment at the beginning of the event with an external predictor, which is easily observable or available. Such predictor avoids modelling the inter-flood phase and simplifies the parametrization and the calibration of the model. The selected model combines a distributed SCS production function and a Lag and Route transfer function, applied to a discretization of the basin in a grid of regular square meshes.The model was first tested on the Real Collobrier watershed. This Mediterranean basin has been monitored by IRSTEA for more than 50 years and has an exceptional density of rainfall and flow measurements. This favourable environment made it possible to reduce the uncertainties on the rainfall input and to evaluate the actual performances of the model. In such conditions, the floods were correctly simulated by using constant parameters for all the events, but the initial condition of the event-based model. This latter was highly correlated to predictors such as the base flow or the soil water content w2 simulated by the SIM model of Meteo-France. The model was then applied by reducing the density of the rain gauges, showing loss of accuracy of the model and biases in the model parameters for lower densities, which are representative of most of the catchments.The spatial variability of the model parameters was then studied in different Real Collobrier sub-basins. The comparison made it possible to highlight and correct the scale effect concerning one of the parameters of the transfer function. The catchment saturation predictors and the initial condition of the model were still highly correlated, but the relationships differed from some sub-catchments. Finally, the spatial variability of the model parameters was studied for other larger Mediterranean catchments, of which area ranged from some tenth to hundreds of square kilometres. Once more, the model could be efficiently initialized by the base flow and the water content w2, but significant differences were found from a catchment to another. Such differences could be explained by uncertainties affecting as well the rainfall estimation as the selected predictors. However, the relationships between the initial condition of the model and the water content w2 were close together for a given type of catchment.In conclusion, this distributed event model represents an excellent compromise between performance and ease of implementation. The performances are satisfactory for a given catchment or a given type of catchment. The transposition of the model to ungauged catchment is less satisfactory, and other catchment saturation indicators need to be tested, e.g. in situ measurements or satellite measurements of soil moisture.

Le modèle conceptuel pluie-débit journalier MEDOR a été développé pour le milieu méditerranéen. L'analyse de sensibilité des processus hydrologiques aux caractéristiques du milieu et du climat lui confère une conception spécifique. Ce modèle à quatre paramètres a été élaboré à partir des données de six bassins méditerranéens, français et libanais. Son calage au moyen du critère de Nash pose des problèmes d'equifinalité. Le balayage exhaustif de l'espace critère montre l'existence de relations entre paramètres de production d'une part et ceux de transfert d'autre part. Un critère à échéance variable détermine la relation d'equifinalité de production, à partir des seuls cumuls annuels de pluie et de débit et de la structure stochastique journalière de la pluie.<br />Couplé à un modèle stochastique adimensionnel de pluie, MEDOR génère une surface représentative des coefficients d'écoulement des bassins dans l'espace des paramètres. Cinq zones ont été définies sur le pourtour méditerranéen avec leurs surfaces de référence. MEDOR possède des propriétés d'agglomération spatiale et temporelle. Les paramètres du modèle pour un bassin formé de l'union de plusieurs sous bassins sont déterminés par des lois d'agglomération utilisant les paramètres de chacun d'eux. Un des paramètres de la fonction de production est lié à la profondeur utile du sol, ce qui permet de transformer le modèle global en un modèle semi-distribué. L'extension à un nombre plus important de bassins permettrai de préciser la signification des paramètres, et donc d'envisager une application à des bassins non jaugés.

Les dommages dus aux inondations sont en constante augmentation malgré la réglementation de l'urbanisation. La politique traditionnelle de lutte contre les crues ne donne pas satisfaction, tant du point de vue de la sécurité que des impacts sur les milieux naturels. La gestion du risque d'inondation doit s'inscrire à l'échelle du bassin versant et respecter le fonctionnement naturel des rivières. Il est nécessaire de développer des outils techniques nouveaux pour l'aide à la décision. Nous proposons une méthodologie de modélisation opérationnelle pour la gestion intégrée du risque d'inondation dans le contexte des rivières de moyenne montagne. La recherche de solutions techniques adaptées à chaque rivière peut s'appuyer sur un diagnostic complet du risque et l'exploration de scénarios de prospective. Pour permettre une approche à grande échelle, nous proposons une méthode de production et de sélection d'événements de référence conservant la même fréquence en tout point du bassin versant. Ces événements sont utilisés en simulation pour produire une cartographie précise de l'aléa en intensité et en fréquence. Le diagnostic du risque d'inondation permet d'identifier les enjeux à protéger et les hydrosystèmes à réhabiliter. Pour répondre à ces attentes, nous avons testé en prospective une renaturation de la rivière. L'application à la vallée de la Thur (68) et la simulation à l'aide d'une rivière virtuelle montrent que ce type de projet ne permet pas de réduire sensiblement l'aléa pour les rivières à pente forte. La pertinence des événements de référence choisis fait l'objet d'une discussion. Nous avons par ailleurs montré que les endiguements ne constituent pas une protection durable contre les crues. L'étude de scénarios de changement climatique annonce une aggravation des crues à moyen terme. Ainsi, nous intégrons dans la méthodologie l'étude de stratégies de résilience, encore peu explorées en France. La société doit réapprendre à vivre avec les crues pour concilier développement économique et équilibre des milieux naturels.