Une incertitude majeure dans la prévision numérique du climat provient des difficultés de représenter convenablement les nuages en phase mixte qui jouent un rôle important dans le processus de formation des précipitations (Wallace et Hobbs, 1977). De nombreuses études ont montré qu'il est nécessaire pour une représentation adéquate des effets indirects des aérosols et des processus microphysiques au sein des nuages, d'utiliser un schéma complet de la microphysique qui inclut d'une manière détaillée la phase glace (Senior et Mitchell, 1993; Lohmann et Feitcher, 1997). Blanchet et Girard (1995) ont montré que la modification des propriétés des nuages par les aérosols pourrait produire un forçage radiatif négatif de 9W/m² dans l'infrarouge en considérant dans leur étude l'impact que pourrait avoir les aérosols de nature anthropique sur la microphysique des nuages mixtes et glacés en Arctique. L'interaction des aérosols avec les nuages de phase mixte et glacée est peu étudiée et donc moins connue. La compréhension de ces nuages dépend du mode d'action des noyaux glaçogènes qui est très complexe (Cooper, 1980), des caractéristiques physiques et chimiques des substances qui peuvent recouvrir ces noyaux. De plus, les techniques de mesures existantes ne permettent pas de bien caractériser ce type de nuage (Bigg, 1990; Vali, 1991; Rogers, 1993; Pruppacher et Klett, 1997). Dans ce travail de recherche, nous avons testé et évalué un nouveau paramétrage de la nucléation hétérogène de la glace élaboré par Khvorostyanov et Curry (2000) avec le modèle local du climat (LCM11) et comparé au paramétrage existant de Meyers et al. (1992) et à ce qui a été observé pendant la campagne de mesure SHEBA (Surface Heat Budget of Arctic). Le modèle local du climat (LCM11) est doté d'un schéma de nuages glacés qui simule les interactions entre les aérosols et les nuages, la formation de cristaux de glace et les précipitations. Le nouveau paramétrage de la nucléation hétérogène de la glace est fonction des caractéristiques des aérosols (paramètre de mouillabilité, sites actifs, tension d'inadaptation élastique). Il est à noter que le paramétrage de Meyers et al. (1992) est empirique, il est fonction soit de la température soit de la sursaturation par rapport à la glace dépendamment du mode de nucléation de la glace. La sensibilité du modèle local du climat (LCM11) a été étudiée en fonction des paramètres de la nouvelle théorie de Khvorostyanov et Curry (2000). Cette étude nous a permis de conclure que l'inadaptation élastique de 1 et 2% et les sites actifs sont des paramètres importants pour simuler de l'eau liquide et la glace dans des proportions raisonnables représentatives de l'Arctique. Les résultats des simulations ont montré aussi que la structure verticale de l'eau liquide, de la glace et de la température est bien représentée par Khvorostyanov et Curry (2000) contrairement à Meyers à condition de faire un choix judicieux des paramètres d'inadaptation élastique et de sites actifs. Cependant, le choix de ces paramètres dépend des propriétés physiques et de la composition chimique des aérosols présents dans la masse d'air qui requièrent encore d'autres mesures et vérifications. Il a été démontré que ce paramétrage peut donner de meilleurs résultats si le type de masse d'air est connu à l'avance ce qui suggère de valider le paramétrage de Khvorostyanov et Curry, (2000) avec ses différents choix de combinaison de paramètre pour d'autres cas observés où la composition des aérosols présents est connue. Il y a un gain associé à l'utilisation du paramétrage de Khvorostyanov et Curry (2000) par rapport à celui de Meyers et al. (1992). Ce dernier peut représenter un cas particulier de la théorie de Khvorostyanov et Curry (2000). ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : Nucléation, Paramétrage de la glace.

Les changements de l'intensité et de la fréquence des extrêmes hydro-climatiques peuvent avoir des impacts significatifs sur les secteurs liés aux ressources en eau. Il est donc nécessaire d'évaluer leur vulnérabilité face aux changements climatiques. Cette étude porte sur l'estimation des changements de la fréquence et l'amplitude des événements de précipitations extrêmes au Canada en utilisant un ensemble de dix simulations de 30 ans effectuées avec le Modèle Régional Canadien du Climat (MRCC), pour une période de référence (1961-1990) et une période future (2040-2071). Les simulations futures utilisent le scénario A2 du SRES. Deux méthodes sont utilisées dans cette étude, avec l'hypothèse de stationnarité tranche de temps (en anglais «time-slice stationarity assumption»): Analyse Fréquentielle Régionale (RFA pour «Regional Frequency Analysis»), qui s'opère à l'échelle des unités statistiquement homogènes des régions climatiques prédéfinies, avec la possibilité de réduction au niveau du point de grille et l'analyse individuelle de point de grille (GBA pour «Grid-box analysis»). Des données d'observations réhabilitées et homogénéisées de 495 stations situées partout au Canada sont utilisées pour vérifier l'homogénéité statistique des régions climatiques canadiennes. Ces données sont également utilisées pour sélectionner la distribution régionale la plus appropriée parmi les cinq distributions candidates aux trois paramètres de la modélisation observée de l, 2, 3, 5, 7 et 10 jours AM (AM pour «Annual Maxima») de la quantité de précipitations (i.e. les extrêmes d'un seul jour et de plusieurs jours), survenues entre les mois d'avril et de septembre, pour la période de 30 ans allant de 1961 à 1990. Les distributions candidates sont les suivantes: la distribution des valeurs extrêmes généralisées (GEV pour «General Extreme Value»), Pareto généralisée (GPA pour «Generalized Pareto»), Logistique généralisée (GLO pour «Generalized Logistic»), Pearson Type 3 (PE3 pour «Pearson Type 3») et Normal généralisée (GNO pour «Generalized Normal»). La validation du modèle de simulation pour les périodes de retour de 20,50 et 100 ans des précipitations extrêmes d'une et de plusieurs journées en comparaison avec les observations durant la période 1961-1990 en utilisant les méthodes de RFA et GBA suggèrent une sous-estimation du MRCC pour une grande partie du Canada. Toutefois, le MRCC a tendance à surestimer légèrement sur la région de YUKON. Les changements de l'amplitude et la fréquence des précipitations extrêmes d'un seul jour et de plusieurs jours au Canada sont estimés en utilisant les deux méthodes, celle de RFA et de GBA. Une estimation d'incertitude sous la forme d'intervalles de confiance des 20, 50 et 100 ans de périodes de retour à l'échelle régionale des cinq paires de simulations de la période de référence et celle du future est effectuée en utilisant la méthode de bootstrap vectoriel (en anglais «nonparametric vector bootstrap resampling method») et ensuite exprimé sous forme d'intervalle de confiance. Les résultats de l'étude ont des implications fortes pour des projets liés à la conception et la gestion des ressources en eau et pour estimer la durabilité des infrastructures existantes dans un climat changeant. ______________________________________________________________________________

Les modèles régionaux de climat (MRC) sont maintenant reconnus comme des outils efficaces pour augmenter la résolution des simulations climatiques à un ordre d'environ 50 km. Cette résolution est jugée nécessaire pour décrire de nombreux impacts liés aux changements climatiques. Un modèle régional de climat (MRC) performe une simulation climatique sur une sous-région du globe. Pour cette raison, un MRC doit utiliser des données pilotes au pourtour de la sous-région pour tenir compte de la circulation atmosphérique qui parcourt le globe. Les données provenant d'un MRC sont de plus en plus utilisées par les décideurs pour prendre des mesures d'adaptation adéquates en fonction des changements climatiques projetés. Pour cette raison, les scientifiques qui livrent les données climatiques aux décideurs doivent fournir les données les plus précises possibles. L'une des incertitudes des projections de changements climatiques produites par les MRC est associée au comportement chaotique du système climatique qui est simulé. En raison de ce comportement, le climat simulé avec un MRC peut diverger d'un autre climat simulé lancé avec des conditions initiales légèrement différentes tout en étant en accord avec les mêmes forçages imposés. La variabilité mesurée entre les simulations lancées avec différentes conditions initiales est appelée variabilité interne. La variabilité interne des MRC est étudiée depuis quelques années. Ces études sont généralement limitées par les ensembles de simulations nécessaires au calcul de la variabilité interne. Ces ensembles comportent généralement peu de simulations qui sont courtes et qui sont intégrées sur de petits domaines. Le premier objectif de cette thèse consistait à vérifier les hypothèses émises dans les articles précédents sur la variabilité interne des MRC en utilisant un ensemble de simulations actuel. Afin de répondre à cet objectif, un ensemble comportant dix simulations de dix ans couvrant l'Amérique du nord a été réalisé avec le modèle régional canadien du climat. L'évolution temporelle de la variabilité interne de cet ensemble pour la pression au niveau moyen de la mer (PNM) et pour la température au niveau de l'abri (TA) révèle un cycle annuel faible avec des valeurs en hiver et au printemps légèrement plus élevées. Ce résultat, qui se distingue de celui des études précédentes; suggère que l'augmentation de la dimension du domaine diminue l'efficacité du contrôle du pilote sur le MRC, ce qui augmente la variabilité interne. Pour les variables étudiées, la variabilité interne augmente généralement d'ouest en est du domaine selon la circulation générale vers l'est. Une comparaison de la variabilité interne avec la variabilité temporelle montre que la variabilité interne pour la PNM est près de son maximum en été dans le nord-est du domaine. Une comparaison du climat de dix ans pour la TA en été de chacune des simulations avec la moyenne d'ensemble montre que la variabilité interne peut avoir un impact important sur le climat simulé. Le deuxième objectif de cette thèse consistait à trouver un outil capable de quantifier le contrôle du pilote sur une simulation effectuée avec un MRC. Cet outil est un traceur qui calcule le temps de résidence des parcelles atmosphériques à l'intérieur du domaine d'un MRC. Ce traceur a été utilisé pour première fois dans l'ensemble de simulations décrit précédemment. Les diagnostiques du climat généré par l'ensemble montre le bon fonctionnement du traceur. Le temps de résidence est plus élevé en été qu'en hiver en liaison avec la circulation atmosphérique qui est plus rapide en hiver. Le temps de résidence augmente d'ouest en est à l'intérieur du domaine aussi en relation avec la circulation atmosphérique générale vers l'est. Un diagramme de disperison généré avec les distributions spatiales de la variabilité interne relative des variables étudiées et du temps de résidence montre que la variabilité interne augmente linéairement avec le temps de résidence similairement en été qu'en hiver. Ce résultat appuie plusieurs études précédentes qui proposaient un lien entre la circulation atmosphérique et la variabilité interne. D'après la relation trouvée, le temps de résidence peut être utilisé comme un indicateur quantitatif du contrôle des données pilotes exercé sur un MRC. On envisage que cet outil sera utile pour les analyses futures effectuées avec les MRC. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : Variabilité interne, Modélisation régionale du climat, Ensemble de simulations, Temps de résidence, Climat nord-américain.

Thèse numérisée par la Direction des bibliothèques de l'Université de Montréal.

Cette thèse vise à rechercher les configurations expérimentales optimales pour étudier la réponse des Modèles Régionaux du Climat à aire limitée (MRC) face à des perturbations de leurs paramètres. Le travail est présenté en deux parties. La première partie aborde le cas d'une comparaison entre les simulations provenant d'un MRC, où un événement météorologique ou saisonnier est mis à l'échelle dynamiquement à partir de données observées (réanalyses). Cette situation implique l'utilisation de périodes d'intégration relativement courtes. Par conséquent, la réponse obtenue dans les moyennes temporelles des simulations par rapport à des modifications aux paramètres a tendance à être noyée dans le bruit quasi-aléatoire provenant de la dynamique chaotique du MRC. La possibilité d'augmenter le rapport signal-bruit par l'application du pilotage spectral ou par une réduction de la taille du domaine est étudiée. L'approche adoptée consiste à analyser la sensibilité des moyennes saisonnières du MRC Canadien (MRCC) face à des perturbations sur deux paramètres variés un à un. Le premier contrôle la convection profonde tandis que le second régit la condensation stratiforme. Les résultats montrent que l'ampleur du bruit diminue avec la réduction de la taille du domaine ainsi que par l'application du pilotage spectral. Toutefois, la réduction de la taille du domaine produit aussi des altérations statistiquement significatives de certains signaux, ce qui favorise l'utilisation de pilotage spectral. La deuxième partie de cette thèse aborde le cas d'une comparaison entre deux simulations d'un MRC en termes du climat simulé. À cet effet, un cadre théorique est développé pour le calcul des statistiques de premier et second ordre sur la différence entre les simulations. Les statistiques de la différence sont décomposées en une composante déterministe et reproductible contrainte par les conditions aux frontières et une composante de bruit provenant de la dynamique interne du MRC. Certaines questions reliées à l'estimation de la différence des moyennes temporelles entre les simulations sont développées en détail. Par exemple, un partage optimal des ressources informatiques entre la taille d'un ensemble et la longueur de la période d'intégration, ou encore l'impact de la t aille du domaine et du pilotage spectral sur l'estimation de la réponse du modèle. Une application de ces considérations théoriques est illustrée à partir de la réponse des simulations du MRCC dont un paramètre lié à la convection profonde a été perturbé. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : modèle régional de climat, perturbation des paramètres, ensemble, composante reproductible, pilotage spectral, taille du domaine, différence de moyennes.

La présente étude a permis d'évaluer le potentiel d'utilisation de deux versions (3.7.1 et 4.1.1) du Modèle Régional Canadien du Climat (MRCC) pilotées en mode réanalyse, afin de reproduire les extrêmes observés en été (sur la période 1961-1990), et de caractériser le régime de précipitation et de températures maximum et minimum au-dessus de trois régions de l'Est de l'Amérique du Nord (Sud et Nord des Grands Lacs et Pennsylvanie). La validation fut réalisée à l'aide de critères diagnostiques multiples liés à la fréquence, l'intensité et la durée des événements extrêmes de précipitation et de température, ceux-ci étant régulièrement utilisés dans l'évaluation des impacts sur les activités humaines et les écosystèmes dans une perspective de changements climatiques. Avec l'utilisation de critères statistiques, les distributions statistiques des indices d'extrêmes de précipitation et de température, les champs moyens saisonniers et la variabilité interannuelle de ces indices et ont été comparés par rapport aux valeurs observées. Les données de référence ont été établies en interpolant sur la grille du MRCC (45km), à l'aide du krigeage ordinaire, les données observées du National Climate Data Center, d'Environnement Canada et du Ministère de l'Environnement du Québec. Les principaux résultats obtenus sur les trois régions montrent que la version 4.1.1 du MRCC simule mieux le régime de température minimum et l'occurrence des jours humides/secs que la version 3.7.1 notamment la moyenne saisonnière, le 10ième centile de température, le nombre maximum de jours secs consécutifs, le nombre de jours de pluie (≥1 mm) ainsi que la variabilité interannuelle. Par contre, peu ou pas de différences existent entre les deux versions quant à la température maximum, l'amplitude thermique quotidienne, le 90ième centile de la température, la précipitation moyenne et le centile de précipitation. Dans tous les cas, la variabilité saisonnière est biaisée (forme et étalement des queues de distribution, i.e. aplatissement et asymétrie), les extrêmes chauds des températures maximum et les précipitations les plus intenses, liées aux phénomènes convectifs, étant largement surestimés ou sous-estimés, respectivement. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : Modèle Régional de Climat, Extrêmes, Krigeage, Variabilité interannuelle.

Cette étude se penche sur l'interaction nuage-rayonnement simulée par le modèle GEM-LAM (modèle Global Environnemental multi-échelle à aire limitée) en évaluant avec des observations provenant de deux sites du programme ARM (Atmospheric Radiation Measurement) les différentes composantes atmosphériques ayant un impact sur le bilan radiatif de surface. Ainsi, le rayonnement vers la surface de courtes et longues longueurs d'ondes est comparé aux observations en fonction de la fraction nuageuse afin d'isoler l'effet de la vapeur d'eau ou de l'eau liquide des nuages sur le rayonnement descendant. À l'aide des cycles diurnes moyens et des distributions de fréquences, le principal biais identifié pour le rayonnement à la surface simulé par GEM-LAM est la surestimation du rayonnement d'ondes courtes incident à la surface vers le milieu de la journée. Ce biais provient, d'une part, d'une sous-estimation de la fraction nuageuse, et d'autre part, d'une trop grande transmissivité du rayonnement solaire des nuages lorsqu'ils sont présents, particulièrement pour les nuages optiquement minces. Le biais radiatif de courtes longueurs d'ondes est responsable d'un biais chaud de température près de la surface pour les saisons d'été aux deux sites. Ceci entraîne un biais positif du rayonnement d'oncles longues pour les conditions de ciel clair qui est toutefois compensé par la sous-estimation de la fraction nuageuse pour donner des biais réduits du rayonnement d'ondes longues pour toutes les conditions. De plus, le biais de courtes longueurs d'ondes pourrait être responsable d'un assèchement excessif de la surface et par conséquent mener à un déficit de vapeur d'eau dans l'atmosphère, particulièrement pour la saison d'été au site SGP. Cette étude illustre l'importance de l'évaluation individuelle des composantes de l'interaction nuage-rayonnement à l'aide de statistiques à hautes fréquences temporelles afin de bien identifier les erreurs compensatoires qui peuvent être présentes. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : Interaction nuage-rayonnement, Modèle régional de climat, Schéma microphysique, Bilan radiatif de surface.

Les nuages de convection restreinte sont omniprésents au-dessus de l'Océan Tropical et durant l'été dans la couche limite des latitudes moyennes . Ils ont un impact sur la structure thermodynamique de la basse atmosphère et, par le biais de leur interaction avec la radiation solaire, ils influencent le budget énergétique à la surface. Dans ce projet, nous introduisons un paramétrage des nuages de convection restreinte dans le Modèle Régional Canadien du Climat (MRCC). Ce paramétrage est couplé au schéma de convection, qui est celui de Bechtold-Kain-Fritsch. La représentation de ces nuages ainsi que du MRCC avec ce paramétrage est évalué au-dessus du Pacific selon le protocole du GEWEX Pacific Cross-section Intercomparaison (GPCI). ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : Convection restreinte, Paramétrage de la couverture nuageuse, Modèles régionaux du climat.

Les Modèles Régionaux du Climat (MRC) ont longtemps été considérés comme des outils performants, de haute résolution à aire limitée permettant une meilleure compréhension du climat passé, présent et futur. Le plus souvent, les MRC sont pilotés à leurs frontières latérales par des Modèles Globaux du Climat (MGC) de basse résolution, et qui couvrent le reste du globe terrestre. Ces modèles ont la particularité de reproduire différentes solutions de l'état de l'atmosphère à cause de leur sensibilité aux conditions initiales (CI). Outre la solution due aux forçages externes (forçage aux frontières latérales et le forçage de surface), les MRC reproduisent une seconde solution associée à la variabilité interne (VI) du modèle du fait de leur sensibilité aux CI. Cette sensibilité est en grande partie causée par la nature non-linéaire de la physique et la dynamique atmosphériques. À l'instar de précédentes études, nous analysons la variabilité interne du Modèle Régional Canadien du Climat (MRCC) en utilisant un ensemble de simulations aux CI différentes. Ce projet de recherche consiste à effectuer un diagnostic quantitatif des termes dynamique et diabatique qui contribuent à la variation temporelle et la distribution spatiale de la VI. L'originalité de ce travail est qu'il propose des équations bilans de la VI pour deux variables atmosphériques: la température potentielle et le tourbillon relatif. Les deux équations établies présentent des termes similaires, notamment les termes relatifs au transport de la VI par l'écoulement de la moyenne d'ensemble et la covariance des fluctuations agissant sur le gradient de la moyenne d'ensemble de la variable considérée. Concrètement, nous avons utilisé un ensemble de 20 simulations aux CI différentes pour analyser les caractéristiques de la VI, afin de déterminer une période et une région d'intérêt caractérisées par une forte croissance de la VI. Ensuite, nous avons validé les équations établies en montrant l'égalité entre les deux parties de chaque équation. Enfin, une étude de bilan a permis d'évaluer la contribution des différents termes au développement et à l'évolution de la VI. Les résultats révèlent que les termes dominants responsables de l'accroissement de la VI sont soit les termes de covariance impliquant les fluctuations de température potentielle et de chauffage diabatique, ou les termes de covariance de fluctuations inter-membres agissent sur le gradient de la moyenne d'ensemble de la variable considérée. Les résultats révèlent également que les épisodes de fortes diminutions de la VI se produisent lorsque les maxima de la VI sont proches de la frontière nord-est, indiquant leur transport en dehors de la zone d'étude par l'écoulement moyen. Enfin, nos résultats ont montré qu'en moyenne, les termes du troisième ordre sont négligeables, mais peuvent devenir importants lorsque la VI est importante. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : Équation de la variabilité interne, modèle régional de climat, ensemble de simulations, climat de l'Amérique du Nord

À cause des non-linéarités dans les processus atmosphériques, les modèles de circulation générale de l'atmosphère (MCG) sont sensibles aux conditions initiales (CI). Autrement dit, les MCG peuvent générer plusieurs solutions de la circulation atmosphérique si un petit changement affecte les CI du modèle. La divergence de ces solutions devient importante en quelques jours de simulations en limitant de cette façon l'habilité des prévisions atmosphériques à une période limitée de temps. Comme une conséquence de l'existence de la variabilité interne, n'importe quel climat peut être décomposé en deux termes: le signal, associé au forçage externe, qui est indépendant des CI, et le «bruit» dû à la variabilité interne qui est générée par cette sensibilité du modèle aux CI. Contrairement aux MCG qui couvrent entièrement le globe et qui possèdent une large variabilité interne, les Modèles Régionaux de Climat (MRC) sont limités et pilotés à leurs frontières latérales par les données qui proviennent d'un MCG. La variabilité interne sera influencée dans ce cas par toutes les contraintes et les conditions qui définissent les MRC comme par exemple, les conditions aux frontières latérales (CFL), la taille ou la position géographique du domaine. Cette étude analyse la variabilité interne du Modèle Régional Canadien de Climat (MRCC) et les possibles conséquences sur les analyses statistiques des simulations saisonnières. Dans ce but, des ensembles de 20 simulations qui diffèrent légèrement dans leurs CI ont été générés avec le MRCC pour différentes tailles du domaine sur une période saisonnière de trois mois. Le degré de variabilité interne du modèle a été apprécié selon le degré d'écart entre les simulations durant la période d'intégration. Les résultats montrent que la variabilité interne du MRCC, estimée statistiquement comme la variance entre les 20 simulations de l'ensemble, dépend fortement des événements synoptiques. Cette dépendance est notamment observée dans les variations de la variance enregistrées durant la période d'intégration. Selon notre étude, ces variations ont une distribution géographique qui dépend de la variable étudiée. L'évolution temporelle des patrons synoptiques montre que les maximums des variances pour les deux champs étudiés, la précipitation (le maximum au sud des États-Unis) et la hauteur du géopotentiel à 850 mb (le maximum vers le nord-est du domaine) sont reliés entre eux. Les évènements de forte précipitation du sud des États-Unis (qui correspondent au maximum de l'écart entre les simulations de la précipitation) engendrent aussi un écart des simulations du champ de la hauteur du géopotentiel. Cet écart entre les simulations du champ de la hauteur du géopotentiel se développe dans les circulations cycloniques et atteint son maximum vers le nord-est (où la variance est maximale), avant de quitter le domaine. Utilisée comme une mesure de la variabilité interne à l'échelle saisonnière, la variance entre les moyennes saisonnières des membres de l'ensemble montre des valeurs importantes pour la statistique saisonnière du champ simulé sur les grands domaines. Cela suggère qu'une seule simulation pourrait produire des erreurs significatives dans l'estimation de la moyenne saisonnière du champ simulé, spécialement pour le champ de la précipitation. Notre expérience montre que les réductions de la taille du domaine diminuent la variabilité interne du modèle dans les aires convectives trouvées comme des aires de forte variabilité interne, mais une importante variation dans la distribution géographique et l'amplitude de la variabilité interne a été détectée à l'intérieur du domaine avec les réductions successives du domaine. Finalement, certains de nos résultats sont reliés au comportement complexe des simulations de l'ensemble durant la période d'intégration, en relevant l'habileté du modèle à générer des solutions bimodales des champs simulés sur les grands domaines d'intégration. La division de l'ensemble en deux groupes différents de simulations a généré pour une période courte de simulation deux estimations différentes du même champ simulé. Les résultats ont montré que cette division de l'ensemble coïncide avec la période durant laquelle la variabilité interne du modèle est maximale. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : modèle régional de climat, ensemble de simulations, variabilité interne.