Nous étudions les flux turbulents de chaleur sensible et de chaleur latente, qui sont mal compris et mal mesurés sur les glaciers, à l'aide de campagnes de terrain déployées dans la zone d'ablation du glacier tropical du Zongo (16°S, Bolivie, 4900-6000 m) durant la saison sèche de l'hiver austral et sur le glacier alpin de Saint-Sorlin (Alpes Françaises, 45°N, 2600-3400 m) durant l'été boréal. Un mât de 6 m permettait la mesure des profils verticaux de vitesse de vent et de température de l'air et des mâts de 2 m comportaient des systèmes de covariances turbulentes (CT). Nous étudions l'évolution temporelle des flux turbulents et l'applicabilité de la méthode aérodynamique des profils en terrain complexe de montagne. Les hypothèses sont discutées via la caractérisation des régimes de vent et de la turbulence. Nous calculons ensuite les flux et les erreurs associées. Sur le Zongo, sous forçage synoptique faible, un écoulement catabatique s'installe de la fin d'après-midi jusqu'au matin, avec un maximum de vitesse de vent à environ 2 m de hauteur. Les forçages synoptiques forts s'alignent approximativement avec le glacier, provoquant un intense écoulement descendant, et dans ce cas nous n'observons pas de maximum de vitesse de vent. Souvent, autour de midi, des vents ascendants sont observés. Sur le glacier de Saint-Sorlin les forçages associés à des épisodes de Foehn ou à des dépressions se déplaçant depuis l'ouest, s'alignent approximativement avec le glacier, générant des vents forts descendants. Quand le forçage synoptique est modéré, un maximum de vitesse de vent est observé nuit et jour 50% du temps. Des vents ascendants sont observés 15% du temps, quand le forçage synoptique est faible. La couche de surface est perturbée par des tourbillons de couches externes sous vent fort, ou de lentes oscillations en écoulement catabatique. Ces perturbations influencent les flux turbulents. Les erreurs aléatoires sur la méthode des profils sont dues principalement à des incertitudes sur la température. L'erreur reste faible sur les flux moyens. La couche de surface est rarement plus épaisse que 2 m et la méthode des profils appliquée à l'aide des mesures plus obtenus plus haut sous-estime les flux de surface de 20% à 70% . Quand un maximum de vitesse de vent est observé, les flux sont sous-estimés même à 2 m. L'influence des perturbations de la couche de surface n'est pas capturée par la méthode des profils, et les flux sont environ 40% inférieurs à ceux mesurés par CT. Ces derniers sont affectés par d'importantes erreurs aléatoires, en raison d'un échantillonnage statistique insuffisant des grands tourbillons. La méthode sous-estime probablement les flux à cause d'une sous-estimation de la vitesse verticale (~15%) et de la divergence verticale des flux. Sur le glacier du Zongo, l'air de haute altitude est très sec et la sublimation (quelques mm d'eau par jour) est un important puits d'énergie à la surface. Le flux de chaleur sensible est un important gain d'énergie la nuit sous l'influence de vents forts (de 30 à 50 W m-2), car l'inversion de température est marquée. Quand un maximum de vitesse de vent est observé, les flux sont faibles (de 5 à 20 W m-2) car la vitesse du vent est faible. La somme des flux turbulents est faible dans ces deux cas car ils sont opposés et les biais se compensent. En vent ascendant, le flux de chaleur sensible est faible (<5 W m-2) car la stratification est neutre, mais le flux de chaleur latente reste important (de -25 à -35 W m-2), le flux net est donc important et les biais ne se compensent plus. Sur le glacier de Saint-Sorlin, le flux de chaleur latente est faible car l'air est humide, et le flux de chaleur sensible peut être intense (~25 W m-2) quand la vitesse du vent est élevée. Le flux net est fort par vent fort et les biais sur les flux calculés par la méthode des profils peuvent être élevés. / We study turbulent fluxes of sensible and latent heat, that are a poorly-known and difficult term to measure over glaciers, with the help of two field campaigns deployed over the ablation zone of Zongo glacier (16°S, Bolivia, 4900-6000 m.a.s.l.) during the austral winter dry season and over the Saint-Sorlin glacier (French Alps, 45°N, 2600-3400 m.a.s.l.) during the boreal summer. A 6-m mast allowing for wind speed and air temperature vertical profile measurements was installed, along with 2-m masts holding eddy-covariance systems. The focus is on the temporal evolution of turbulent fluxes and the applicability of the aerodynamic profile method in the complex terrain of high mountains. The assumptions of the method are discussed by characterizing the wind regimes and the turbulence. We then compute fluxes and associated errors. Above Zongo glacier, under weak synoptic forcing, katabatic flows are observed from late afternoon to early morning, with a wind-speed maximum at around 2 m. Strong synoptic forcing roughly aligns with the glacier, leading to strong downslope flows for which no wind-speed maximum is observed. Most of the days around noon, upslope flows are observed. On Saint-Sorlin glacier in summer, flows associated with low-pressure systems coming from the west or Foehn events roughly align with the glacier, leading to strong downslope winds. Wind-speed maxima are observed night and day, ~50% of the time, when synoptic forcing is moderate. Upslope flows are observed 15% of the time, when synoptic forcing is weak. The surface layer is disturbed by outer-layer eddies in strong flows and by slow oscillations if katabatic flow prevails. These disturbances influence turbulent fluxes. Random errors on the fluxes derived from the profile method are mainly due to temperature uncertainties. Errors remain small on the mean fluxes. The surface layer is rarely deeper than 2 m on both glaciers and the profile method with measurements made above that height underestimates the surface fluxes by 20% to 70%. When a wind-speed maximum is observed, fluxes are underestimated even at 2 m. The influence on the fluxes of the surface-layer disturbances is not captured by the profile method, and fluxes are about 40% smaller than the eddy-covariance fluxes. The latter are affected by large random errors due to inadequate statistical sampling of large-scale eddies and are probably underestimated, mainly due to vertical wind speed underestimation (~15%) and to vertical flux divergence. Above Zongo glacier, due to the dry high-elevation air, sublimation (a few millimeters w. e. per day) is a large energy loss for the surface. Sensible heat flux is a large energy gain in strong nocturnal downslope flows (from 30 to 50 W m-2) and strong winds, due to a marked temperature inversion. When a wind-speed maximum is observed, low wind speeds cause small turbulent fluxes (from 5 to 20 W m-2). The sum of turbulent fluxes is small in those two cases because the fluxes are opposed in sign and the biases mostly compensate. In upslope flows, the sensible heat flux is small (<5 W m-2) due to near-neutral stratification, but latent heat losses remain large (around -25 to -35 W m-2), so that the net turbulent flux is large and the biases do not compensate. Above Saint-Sorlin glacier, the latent heat flux remains small because the air is generally humid, whereas the sensible heat flux can be large (~25 W m-2) when wind speed is high. The net flux is large when wind speed is high, and the biases on net turbulent fluxes derived from the profiles can be significant.

La simulation de réacteurs à bulles en régime industriel est un grand défi. L'objectif principal de ce travail est la prédiction de la taille des bulles à l’aide d’un modèle numérique de bilan de population, basé sur la modélisation des phénomènes de brisure et de coalescence, et pouvant être couplé aux conditions hydrodynamiques présentes dans les réacteurs. Différentes données expérimentales sont obtenues pour valider le modèle. La taille des bulles est mesurée à l'aide d'une technique innovante de corrélation croisée. Les essais, réalisés en eau du réseau (partiellement contaminée) et en eau déminéralisée avec ajout éventuel d'éthanol, montrent que les additifs réduisent la coalescence et diminuent la taille moyenne des bulles. Deux distributeurs du gaz différents sont utilisés pour découpler l'étude de la brisure et de la coalescence. Les données expérimentales sont utilisées initialement pour valider des simulations CFD 3D transitoires Eulériennes-Eulériennes. La loi de traînée est corrigée par un facteur de swarm pour intégrer l’effet d’une fraction de gaz élevée. Différents modèles de turbulence sont testés. La contribution de la turbulence induite par les sillages de bulles au mélange de scalaires est évaluée. Enfin, pour prédire la taille des bulles, un bilan de population est couplé au modèle hydrodynamique préalablement validé et est résolu par la méthode de quadrature des moments (QMOM). Un set original de kernels de brisure et coalescence est proposé, capable de prédire la taille des bulles pour différentes conditions opératoires. Le comportement du modèle lors de l’extrapolation des réacteurs est également examiné / The simulation of bubble column reactors under industrial operating conditions is an exciting challenge. The main objective of this work is to predict the bubble size, in turn interconnected to the reactor hydrodynamic conditions, with computational models, by modelling bubble breakage and coalescence. Experimental data is collected for model validation, including bubble size measurements with an innovative cross-correlation technique. Experiments are carried out with tap water and demineralized water, with or without the addition of ethanol, and gathered results show that additives reduce coalescence and lower the mean bubble size. Two different spargers are used, in order to decouple the investigation of breakage and coalescence. The experimental data set is used to validate out unsteady three-dimensional Eulerian-Eulerian CFD simulations. A drag law for oblate bubbles is considered, together with a swarm factor, that accounts for the swarm effect. Several turbulence models are tested. The contribution of bubble induced turbulence (BIT) to scalar mixing is assessed. To predict bubble size, a population balance model is coupled to the hydrodynamic model and is solved with the quadrature method of moments. A set of breakage and coalescence kernels is proposed, capable of predicting the bubble size for different operating conditions. Scale-up effects are also investigated

L'objectif ultime de ce travail de thèse consiste à établir un bilan énergétique à l'échelle du grain afin de caractériser et de vérifier la cohérence thermodynamique de modèles de comportement utilisés pour rendre compte du développement de la plasticité cristalline dans les matériaux métalliques. La première partie de ce travail a consisté à mettre en place un protocole d'élaboration du matériau permettant d'obtenir la microstructure souhaitée, compatible avec des moyens d'observations macroscopiques. Les échantillons d'aluminium à très gros grains (centimétriques) ainsi obtenus sont utilisés pour effectuer des essais cycliques durant lesquels les champs cinématiques et thermiques sont mesurés au moyen de techniques de Corrélation d'Images Numériques et de Thermographie Infra-rouge. Deux techniques de traitement d'image spécifiques ont été proposées. Elles permettent d'introduire des hypothèses sur les champs cinématiques et thermiques qui soient adaptées à la microstructure (ici continuité intra-granulaire du déplacement, de la température et du flux). Ces méthodes permettent d'accéder à des mesures complètement indépendantes d'un grain à l'autre tout en améliorant la robustesse des méthodes de mesure. Ces méthodes ont été validées numériquement en utilisant des images de synthèse sur lesquelles ont été appliqués des champs hétérogènes. Une campagne d'essais cycliques a enfin été menée sur les multi-cristaux d'aluminium élaborés. Les méthodes développées ont permis d'observer le développement de la plasticité intra-granulaire et le développement de la fissuration inter-granulaire. / The main objective of this PhD thesis is to establish an energy balance at the grain scale in order to assess the thermomechanical consistency of material models used to predict the development of crystal plasticity of metals.The first part of this work consists in setting a protocol allowing the material elaboration with the desired microstructure which is to be compatible with the use of classical macroscopic observation devices. The obtained coarse-grained aluminum samples (with centimeter grains) are used in cyclic tensile tests. During these tests, the kinematic and thermal fields are recorded with Digital Image Correlation and Infra-Red Thermography techniques.Two specific imaging techniques were developed. They allow introducing ad hoc hypotheses (i.e. consistent with microstructure) on the kinematic and the thermal fields. In this work, these hypotheses consist in intra-granular continuity conditions on the displacement, temperature and heat flux fields. These methods give independent measures on each grain while improving the robustness of the measurement methods. These methods were numerically validated using computer-generated images heterogeneously loaded.Cyclic tests were finally performed on the processed aluminum multi-crystals. The developed methods allowed the observation of the development of intra-granular plasticity and the development of inter-granular cracking.

L’agrivoltaïque, association sur une même surface d’une production agricole et d’une production d’énergie d’origine photovoltaïque, apparait aujourd’hui comme une solution innovante pour atténuer les effets du changement climatique, notamment sur le secteur agricole. Déjà imaginé en 1982, les premières expérimentations débutées à Montpellier (France) en 2010 ont montré la pertinence de cette association à travers un maintien des rendements sous certaines conditions, une meilleure efficience d’utilisation du sol ainsi qu’une diminution de la consommation en eau des cultures sous ombrage. Suite à ces travaux pionniers réalisés sous des panneaux photovoltaïques à inclinaison fixe, l’utilisation de panneaux dynamiques à inclinaison variable est apparue nécessaire pour réduire l’hétérogénéité du rayonnement disponible pour la culture mais aussi, adapter la stratégie d’ombrage aux besoins radiatifs de la culture au cours de son développement.La présente thèse avait pour objectif de caractériser et de modéliser l’impact de la présence de panneaux solaires sur le bilan hydrique d’une culture conduite sous dispositif agrivoltaïque dynamique en vue d’optimiser la stratégie d’irrigation et le pilotage de l’inclinaison des panneaux. Les expérimentations menées sur une culture de laitues ont mis en évidence l’intérêt de panneaux dynamiques afin de réduire les hétérogénéités radiatives. Le développement d’un modèle de redistribution des pluies par les panneaux photovoltaïques a permis la mise en œuvre d’une stratégie de pilotage de l’inclinaison en temps réel a permis d’homogénéiser les cumuls observés au sol. L’adaptation d’un modèle de bilan hydrique et de développement de culture par l’ajout d’un module décrivant la dynamique stomatique sous ombrage fluctuant a permis de décrire la consommation en eau de la culture et son développement sous différentes stratégies de pilotage de l’ombrage. Enfin, différentes stratégies de pilotage de l’inclinaison des panneaux photovoltaïque ont pu être simulées et évaluées à l’aide d’un indicateur globale intégrant l’efficience d’utilisation du sol, la productivité de l’eau, le décalage de maturité et les sources d’hétérogénéités pouvant affecter la production de laitues. / Agrivoltaism, defined as the association on the same land of agricultural and photovoltaic energy production, appears as an innovating concept to dampen some of the effects of climate change, in the agricultural sector. Although the concept was already imagined in 1982, the first experimentations started in 2010 at Montpellier (France) and showed the relevance of this combination by the maintenance of crop yield under certain conditions, the increase of land use efficiency and a reduction of water consumption for the tested crops. Following this pioneering work done under fixed (but not horizontal) photovoltaic panels, the use of "dynamic" panels, i.e. panels with a variable tilting angle, appears necessary to reduce the spatial heterogeneity of the transmitted radiation but also to adapt the shading strategy to the radiation amount required for crop growth.This thesis aims to characterize and to model the impact of the photovoltaic panels on the water budget of the cultivated plot and to progress towards the optimization of irrigation strategies in such systems controlled by the variations in time of the tilting angle of the panels. Experimentations conduced on lettuces showed the benefits of "dynamic" photovoltaic panels to reduce the radiative heterogeneity. Accounting for rain redistribution by the solar panels permits the implementation of a real time strategy to reduce rainfall heterogeneity on the ground surface. The derivation of a water budget and crop development model which describes the dynamics of stomatal conductance under fluctuating shading allows a better simulation of water consumption and crop development for different shading strategies. Finally, various strategies for the piloting of the solar panels could be tested and evaluated by a new, global index combining land use efficiency, water productivity, maturity delays and heterogeneities (in rain and radiation) which can impact production.

Le modèle de neige est souvent dépendant du modèle hydrologique avec lequel il est couplé, ce qui peut favoriser la représentation du débit au détriment de celle de la neige. L'objectif est de rendre le calage du modèle de neige plus indépendant de celui du modèle hydrologique en restant facilement utilisable en opérationnel. Dans cette optique, un modèle contraint sur des données d'observations de la neige permettrait d'améliorer d'une part la robustesse des paramètres du modèle de neige et d'autre part la simulation de l'état du manteau neigeux. Dans la première partie de cette thèse, nous avons étudié et modifié le modèle degrés-jour semi-distribué CemaNeige afin qu'il puisse simuler de manière plus réaliste la variable de surface d'enneigement du bassin versant. Cette modification, couplée au calage du modèle sur des données de surface enneigée et sur le débit, a permis d'améliorer la simulation de l'enneigement par le modèle sans détériorer significativement les performances en débits. Nous alors ensuite débuté le développement d'un nouveau modèle de neige à l'échelle ponctuelle. Celui-ci se compose d'un modèle de rayonnements, simulant les rayonnements incidents à partir de données d'amplitude de températures journalières, et d'un modèle de manteau neigeux. Le modèle de manteau neigeux résout les équations de la chaleur au sein du manteau neigeux à l'aide d'une représentation spectrale du profil de température. Cette représentation permet de simuler les profils et gradients de températures en utilisant moins de variables d'état qu'une discrétisation verticale par couches. Pour mieux prendre en compte les mesures ponctuelles de neige, ce modèle devra être distribué. / Snow models are often dependent on the hydrological model they are coupled with, which can promote higher performance on runoff simulation at the expense of snow state simulations performances. The objective of this thesis is to make the calibration of the snow model more independent from the calibration of the hydrological model, while remaining easily usable for runoff forecasting. Calibrating snow model on observed snow data would on one hand improve the robustness of the snow model parameters and on the other hand improve the snowpack modelling. In the first part of this manuscript, we modified the semi-distributed CemaNeige degree-day model so that it can explicitly simulate the watershed snow cover area. This modification coupled with the calibration of the model on snow cover area data and on river runoff data significantly improved the simulation of the snow cover area by the model without significantly deteriorating the runoff performances. Then we started the development of a new point scale snow model. It is based on a radiation model, which simulates incoming radiations from daily temperature range data, and a snowpack model. The snowpack model solves the heat equations within the snowpack by using a spectral representation of the temperature profile. This representation simulates the temperature profile and gradients using fewer state variables than a vertical discretization of the snowpack. In order to be able to use point scale snow observations in the model, it should be distributed on the watershed.

L'étude de la dynamique et de la thermodynamique à petite meso-échelle est importante pour l'évaluation de la climatologie urbaine et de la qualité de l'air. Dans les modèles numériques, les hétérogénéités de la canopée urbaine sont intégrées sur chaque maille au sol ; il convient donc de connaître l'influence du degré de description du sol sur les résultats. Dans ce but on a créé un modèle de sol urbain SM2-U détaillant les flux de chaleur à l'interface canopée-atmosphère à l'échelle du quartier, en vue de simuler la Couche Limite Urbaine (CLU) à haute résolution avec le modèle SUBMESO. La première partie décrit les caractéristiques morphologiques et climatologiques d'une zone urbaine. La deuxième partie aborde la notion d'hétérogénéité de surface sur l'ensemble des échelles spatio-temporelles de la troposphère, puis analyse les méthodes les prenant en compte dans les simulations. On présente alors une étude numérique de l'influence du détail aérodynamique d'une surface hétérogène, de type damier de rugosités, sur la dynamique de l'atmosphère. La troisième partie présente SM2-U. Afin d'évaluer les flux de chaleur à l'échelle du quartier, SM2-U modélise les surfaces urbaines, et prend en compte de manière originale les effets de la canopée urbaine. Le modèle a été testé dans deux configurations avec un forçage météorologique sans rétroaction entre le sol et l'atmosphère : - comparaison du comportement hydrique de SM2-U avec celui du modèle hydrologique urbain EHU du LCPC, et avec des mesures effectuées sur un site péri-urbain, - étude du comportement de SM2-U sur cinq quartiers urbains. Dans la quatrième partie, on présente les premières simulations de la structure de la CLU avec SUBMESO couplé à SM2-U (avec rétroaction sol-atmosphère). Une première étude de sensibilité met en relief l'impact de la représentation du sol urbain sur la structure de la CLU, sur les bilans d'énergie et sur les températures de surface.

La consommation importante de nos réserves énergétiques mondiale ainsi que la sensibilisation actuelle aux problèmes d'écologie imposent aux fabricants d'augmenter le rendement de leurs systèmes électromécaniques. L'étude menée dans ce rapport s'intéresse au calcul des pertes magnétiques, méconnues bien que représentant une part non négligeable dans les bilans de puissance. Pour ce faire, nous avons astucieusement instrumenté magnétiquement et thermiquement une machine asynchrone. Ceci nous permet de localiser et d'étudier l'évolution des inductions lors de conditions de fonctionnements variées. A cette fin, un banc d'essai nous permettant d'agir sur divers paramètres (tension, fréquence, couple) a été réalisé. Dans le même temps, une modélisation numérique fine du système est menée avec le logiciel Flux2D. Elle intègre le maximum de phénomène électrique et mécanique : alimentation, cage d'écureuil, effets de bord, température des conducteurs, mouvement du rotor. La cartographie magnétique ainsi obtenue est comparée et validée par rapport aux résultats issus de l'instrumentation. Un modèle de calcul d'hystérésis dynamique basé sur la connaissance de l'induction et de sa dérivée temporelle est ensuite appliqué à chaque élément du maillage. Ce dernier permet une représentation fidèle des phénomènes magnétiques statiques et dynamiques des tôles. Validé sur des geometries simples pour différentes formes d'ondes, il est ici appliqué au cas complexe d'une machine tournante. Les résultats obtenus à vide comme en charge confirment ses performances vis-à-vis des modèles utilisés actuellement dans l'industrie.

Sur la Cordillère Royale en Bolivie, nous avons reconstitué les bilans de masse de 21 glaciers représentatifs, par photogrammétrie numérique. Les MNT's ont été construit à partir des clichés aériens de l'IGM-Bolivie de 1956, 1963, 1975, 1983, 1997 et 2006, avec une précision décimétrique ou métrique.<br />L'analyse des résultats sur le glacier du Zongo montre qu'un tiers de la surface du glacier (zone d'ablation) contrôle 80% du bilan de masse annuel. La série du Zongo constitue la série de bilan de masse glaciaire la plus longue et la plus précise des Andes tropicales.<br />Nous avons analysé les bilans de masse des 21 glaciers depuis 1956, lorsqu'on s'affranchit des tendances de ces glaciers, les variations temporelles du bilan de masse montrent un signal commun d'origine climatique.<br />Nous avons montré ainsi que 88% de la variance des bilans cumulés sont expliqués par l'exposition et l'altitude. Cette corrélation et les restitutions photogrammétriques ont permis de calculer les bilans de masse et la perte de surfaces de 376 glaciers au cours de la période 1963-2006 (43% de perte de leur volume et 48% de perte de leur surface depuis 1975).<br />Enfin, notre étude évalue la contribution des glaciers au système d'alimentation en eau potable de la ville de La Paz. En moyenne sur la période 1975-2006, la fonte glaciaire a contribué à hauteur de 18% et 15% sur la période récente 1997-2006 (12% en saison humide et 27% en saison sèche). Nous avons enfin calculé les débits attendus en cas de disparition totale des glaciers (12% en moins à l'année, 9% en moins en saison humide et 25% en moins en saison sèche).

Les objectifs de ce travail de thèse étaient : (i) identifier les principaux compromis et covariations d'un ensemble de traits foliaires relatifs à l'acquisition du carbone et (ii) caractériser la variabilité de ces traits en réponse à l'éclairement (plasticité phénotypique) et la diversité interspécifique de cette plasticité chez des espèces de forêt tropicale humide. Les mesures ont été réalisées sur des semis de 13 espèces d'arbres élevés en conditions contrôlées dans 3 traitements d'éclairements relatifs distincts. Une analyse coût-bénéfice des traits du bilan de carbone foliaire a été centrée sur les coûts de construction du métamère (CCmetm), le temps de retour sur investissement (PBT) et la mesure des durées de vie des feuilles (LLS). Une forte diversité interspécifique a été observée pour l'ensemble des traits étudiés. Les espèces se répartissent le long d'un axe opposant durée de vie des feuilles et LMA (masse surfacique) d'un côté, aux capacités photosynthétiques (Asat), PNUE (Asat/N), respiration (Rd) et teneurs en azote (Nm) de l'autre. Ces résultats sont en adéquation avec le schéma universel d'acquisition des ressources proposé par Wright et al. (2004). Les relations entre LLS et PBT sont faibles et l'ensemble des espèces amortissent très largement leur CCmetm. De manière générale, le classement des espèces est maintenu d'un traitement à l'autre et, à l'exception des coûts de construction, seuls de faibles effets d'interaction ont été observés entre espèces et traitements. Le calcul d'un indice de plasticité a également permis de montrer qu'en dépit d'une certaine diversité interspécifique de la plasticité phénotypique, aucune espèce n'était globalement plus « plastique » que les autres sur l'ensemble des traits considérés. Enfin, si les espèces pionnières se distinguent nettement par leurs attributs, leur degré de plasticité n'est pas différent de celui des espèces de sous-bois.

Aux moyennes latitudes, les lacs sont le lieu de formation et de préservation d'archives paléoclimatiques nécessaires à la compréhension des changements du climat passé. Par ailleurs, les lacs d'eau douce ont également un rôle important dans l'économie régionale du fait des activités de pêche, des activités touristiques et de leur utilisation pour l'alimentation en eau potable.<br />Durant cette thèse, nous avons cherché à développer un outil de modélisation des lacs conceptuellement simple et nécessitant un minimum de données d'entrées et de paramètres spécifiques au lac. Ce modèle a été mis au point dans le cas de deux lacs européens (Lac d'Annecy, France, et Ammersee, Allemagne) représentatifs de la majorité des lacs de moyennes latitudes. Le modèle inclut la représentation au pas de temps journalier des processus thermiques (température de l'eau en surface et au fond et fonctionnement de mélange de la colonne d'eau), hydrologiques et isotopiques (δ18O) à l'échelle du lac et de son bassin-versant. Le modèle a été mis au point sur l'actuel grâce aux observations hydro-isotopiques réalisées en parallèle. Les processus sont correctement reproduits aux échelles saisonnière (simulations à l'équilibre) et inter-annuelle (simulations transitoires au cours du 20ième siècle).<br />Des expériences de sensibilité ont été effectuées à l'aide de ce modèle pour quantifier les incertitudes des reconstructions paléoclimatiques basées sur les tests carbonatés d'ostracodes benthiques (Candona sp.) dues (i) aux fonctionnements thermiques des lacs, (ii) aux changements de la saisonnalité des précipitations et, (iii) à la déforestation sur les bassin-versants.<br />Le modèle a été ensuite forcé par un scénario de changement climatique au 21ième siècle. Ces simulations transitoires prévoient des changements très importants des écosystèmes des deux lacs étudiés. Ainsi, le lac d'Annecy subirait un réchauffement de l'eau profonde de 2°C au cours du 21ième siècle, ce qui mettrait en danger la faune d'eau froide. Ammersee, celui-ci, pourrait encourir des arrêts brutaux de mélange à partir de 2020 avec un arrêt de la ventilation des eaux profondes et donc l'asphyxie des faunes benthiques. Ce modèle est la première étape en vue d'un modèle intégré physique, chimie et biologie des écosystèmes lacustres nécessaire à la prévision de la vulnérabilité de la biocénose lacustre aux changements climatiques à venir.