A Retrospective Analysis of Ozone Formation in the Lower Fraser Valley,British Columbia,Canada. Part I: Dynamical Model Evaluation

We conduct a retrospective study of ozone formation in the Lower Fraser Valley (LFV), using numerical models, observations, and emission inventories in order to understand relationships between reductions in local precursor emissions and episodic ozone concentrations. Because there appears to be little or no impact from precursor emissions upwind of the LFV during ozone episodes and because background concentrations of ozone and its precursors are generally from the North Pacific Ocean and quite low, summertime ozone formation in the LFV is mostly caused by local emissions. The observed change in behaviour of ozone formation must, therefore, arise from reductions in local precursor emissions. We exploit the observed changing precursor emission–ozone concentration relationship to perform a dynamical model evaluation. Complicating the analyses are an observed shift in the population patterns within the valley over the last 25 years and a small but documented change in the tropospheric background concentration of ozone. Ozone formation for four episodes, which capture the observed changes in ozone reduction and the different meteorological types that occur during LFV ozone events, are investigated using the Weather Research and Forecasting (WRF)-Sparse Matrix Operator Kernel Emissions (SMOKE)-Community Multiscale Air Quality (CMAQ) modelling system. In order to provide realistic simulations of past events, the SMOKE emission inventory is adjusted to account for temporal changes in the amount of emissions and locations of emission sources. Model output is compared with continuous observations, data collected from field campaigns, and previous modelling efforts. The WRF-SMOKE-CMAQ modelling framework is able to capture the changes in both the magnitude of the ozone concentrations and its spatial behaviour over the period of study. Many of the simulations show that the highest ozone concentrations occur outside the area sampled by the fixed monitoring network and within the LFV's numerous tributary valleys. Not all modelled episodes achieved the same agreement with observations and some of these discrepancies are likely related to shortcomings in the meteorological modelling. The model consistently overpredicts ozone at a number of stations within the City of Vancouver and underpredicts daytime NO_x concentrations there. Both results are consistent with a deficiency in NO_x emissions. The model shows a changing bias over time which also suggests uncertainties in the emission backcasting.

RÉSUMÉ?[Traduit par la rédaction] Nous menons une étude rétrospective de la formation de l'ozone dans le vallée du bas Fraser (VBF) à l'aide de modèles numériques, d'observations et d'inventaires d’émissions dans le but de comprendre les relations entre les réductions dans les émissions locales de polluants précurseurs et les concentrations épisodiques d'ozone. Parce qu'il semble n'y avoir que peu ou pas d'impact des émissions de polluants précurseurs en amont de la VBF durant les épisodes d'ozone et parce que les concentrations de fond de l'ozone et de ses précurseurs proviennent généralement du Pacifique Nord et sont très faibles, la formation d'ozone en été dans la VBF est principalement attribuable aux émissions locales. Le changement de comportement observé dans la formation d'ozone doit par conséquent résulter de réductions dans les émissions locales de précurseurs. Nous exploitons la relation changeante observée entre les émissions de précurseurs et la concentration de l'ozone pour effectuer une évaluation par modèle dynamique. Un déplacement observé dans les configurations de population à l'intérieur de la vallée au cours des 25 dernières années et un changement, petit mais documenté, dans la concentration troposphérique de fond de l'ozone viennent compliquer les analyses. Nous étudions la formation d'ozone lors de quatre épisodes qui capturent les changements observés dans la réduction d'ozone et les différents types météorologiques qui se produisent pendant les événements d'ozone dans la VBF en nous servant du système de modélisation SMOKE (Sparse Matrix Operator Kernel Emissions) – CMAQ (Community Multiscale Air Quality) du WRF (Weather Research and Forecasting). Afin de fournir des simulations réalistes des événements passés, nous avons ajusté l'inventaire d’émissions SMOKE pour tenir compte des changements au cours du temps dans la quantité d’émissions et dans la position des sources d’émissions. Nous comparons la sortie du modèle avec les observations continues, les données recueillies lors d’études sur le terrain et les efforts de modélisation précédents. Le cadre de modélisation WRF-SMOKE-CMAQ est capable de capturer les changements dans la grandeur des concentrations d'ozone ainsi que dans son comportement spatial durant la période de l’étude. Plusieurs des simulations montrent que les plus fortes concentrations d'ozone se produisent en dehors de la région échantillonnée par le réseau fixe de surveillance et à l'intérieur des nombreuses vallées affluentes de la VBF. Tous les épisodes modélisés n'ont pas exhibé le même accord avec les observations et certaines de ces divergences sont vraisemblablement dues à des lacunes dans la modélisation météorologique. Le modèle surprévoit constamment l'ozone à certaines stations dans la ville de Vancouver et sous-prévoit les concentrations de NO_x le jour à cet endroit. Les deux résultats sont cohérents avec un déficit dans les émissions de NO_x. Le modèle montre un biais qui change avec le temps, ce qui porte à croire à des incertitudes dans les émissions rétrospectives.     

Observations dans la couche limite atmosphérique de l'Estuaire du St‐Laurent: Résultats préliminaires

Résumé

Cet article fait état d'une première étude de la structure verticale de la couche limite atmosphérique (CLA) dans l'Estuaire du St‐Laurent, à partir de ballons captifs. Les résultats préliminaires de l'expérience réalisée en août 1977 démontrent que cette CLA est essentiellement formée d'une couche interne très stable (CLI) d'épaisseur variable et d'une couche mélangée correspondant à la partie non perturbée de la CLA côtière. Dans la couche de surface (CS), nous avons observé des gradients de température pouvant atteindre 1°C m^‐1 et des vents généralement très élevés. Le sommet de la CS correspond très souvent à un maximum de vent. Le fait que la direction du vent change considérablement avec l'altitude indique que les effets de vallée et de baroclinité sont très importants dans la CLI.     

Long-Term Variability of Volume and Heat Transport in the Nordic Seas: A Model Study

This article presents results from a model study of interannual and decadal variability in the Nordic Seas. Fifty years of simulations were conducted in an initial condition ensemble mode forced with the National Centers for Environmental Prediction (NCEP) reanalysis. We studied two major events in the interannual and interdecadal variability of the Nordic Seas during the past fifty years: the Great Salinity Anomaly in the 1960s and early 1970s and the warming of the Arctic and subarctic oceans in the late 1990s.

Previous studies demonstrated that the Great Salinity Anomaly observed in the subarctic ocean in 1960 was originally generated by intensified sea-ice and freshwater inflow from the Arctic Ocean. Our model results demonstrate that the increase in the transport of fresh and cold waters through Fram Strait in the 1960s was concurrent with a reduction in the meridional water exchange over the Greenland–Scotland Ridge. The resulting imbalance in salinity and heat fluxes through the strait and over the ridge also contributed to the freshening of the water masses of the Nordic Seas and intensified the Great Salinity Anomaly in the Nordic Seas.

The warming of the Atlantic Waters in the Nordic Seas and Arctic Ocean during the past two decades had an important impact on the variability of these two ocean basins. Some previous observational and model studies demonstrated that the warming of the subpolar Atlantic Ocean in the late 1990s and the meridional transport of the Atlantic Water mass (AW) into the Nordic Seas and Arctic Ocean contributed to this process. At the same time, observations show that the warming of the AW in the Nordic Seas started in the 1980s (i.e., earlier than the warming of the subpolar North Atlantic Ocean). Our model results suggest that this process was triggered by an imbalance in the lateral heat fluxes through Fram Strait and over the Greenland–Scotland Ridge. In the late 1980s the AW transport over the Greenland–Scotland Ridge was stronger than normal while the exchange through Fram Strait was close to normal. The related imbalance in the lateral heat fluxes through the strait and over the ridge warmed the Nordic Seas and caused an increase in the temperature of the AW inflow to the Arctic Ocean in the late 1980s (i.e., about a decade earlier than the warming of the source of the AW in the subpolar North Atlantic Ocean). Thus the model results suggest that the imbalance in lateral heat and salinity fluxes through the strait and over the ridge connecting the Nordic Seas to the North Atlantic and Arctic oceans could amplify the interannual variability in the subarctic ocean.

[Traduit par la rédaction] Cet article présente les résultats d'une étude par modèle de la variabilité interannuelle et décennale dans les mers nordiques. Nous avons effectué des simulations sur une période de cinquante ans en mode d'ensemble de conditions initiales forcé avec les réanalyses des NCEP (National Centers for Environmental Prediction). Nous avons étudié deux événements majeurs survenus dans la variabilité interannuelle et décennale des mers nordiques au cours des cinquante dernières années : la grande anomalie de salinité des années 1960 et du début des années 1970 et le réchauffement des océans Arctique et subarctique vers la fin des années 1990.

Des études précédentes ont démontrées que la grande anomalie de salinité observée dans l'océan subarctique en 1960 a été causée par une intensification de l'apport de glace de mer et d'eau douce depuis l'océan Arctique. Les résultats que nous avons obtenus du modèle montrent que l'accroissement du transport d'eau douce et froide à travers le détroit de Fram dans les années 1960 s'est produit en même temps qu'une réduction dans l’échange méridien d'eau au-dessus de la crête Groenland–Écosse. Le déséquilibre résultant dans les flux de salinité et de chaleur à travers le détroit et au-dessus de la crête a aussi contribué à l'adoucissement des masses d'eau des mers nordiques et a intensifié la grande anomalie de salinité dans les mers nordiques.

Le réchauffement des eaux atlantiques dans les mers nordiques et dans l'océan Arctique au cours des deux dernières décennies a eu un impact important sur la variabilité de ces deux bassins océaniques. Des études observationnelles et par modèle précédentes ont établi que le réchauffement de l'océan Atlantique subpolaire dans les années 1990 et le transport méridien de la masse d'eau atlantique dans les mers nordiques et dans l'océan Arctique ont contribué à ce processus. En même temps, les observations montrent que le réchauffement des eaux atlantiques dans les mers nordiques a commencé dans les années 1980 (c.–à–d. plus tôt que le réchauffement de l'océan Nord-Atlantique subpolaire). Les résultats du modèle suggèrent que ce processus a été déclenché par un déséquilibre dans les flux de chaleur latéraux à travers le détroit de Fram et au-dessus de la crête Groenland–Écosse. À la fin des années 1980, le transport des eaux atlantiques au-dessus de la crête Groenland–Écosse était plus fort que la normale alors que l’échange à travers le détroit de Fram était près de la normale. Le déséquilibre résultant dans les flux de chaleur latéraux à travers le détroit et au-dessus de la crête a réchauffé les mers nordiques et causé une augmentation de la température des eaux atlantiques parvenant à l'océan Arctique à la fin des années 1980 (c.-à-d. environ une décennie avant le réchauffement de la source d'eaux atlantiques dans l'océan Nord-Atlantique subpolaire). Donc, les résultats du modèle suggèrent que le déséquilibre dans les flux de chaleur et de salinité latéraux à travers le détroit et au-dessus de la crête reliant les mers nordiques à l'Atlantique Nord et à l'Arctique pourrait amplifier la variabilité interannuelle dans l'océan subarctique.     

Influence des rejets atmosphériques d'anhydride sulfureux sur la qualité de l'eau des lacs de la région de Rouyn‐Noranda*

Résumé

En mars 1982, le ministère de l'Environnement du Québec a entrepris un échantillonnage physico‐chimique de 64 lacs de tête de la région de Rouyn‐Noranda. Les objectifs de l'étude étaient de déterminer la sensibilité des lacs de cette région et de délimiter l'aire d'influence de la plus importante source anthropique d'émission d'anhydride sulfureux au Québec.

Le traitement graphique et statistique des résultats, ainsi que l'application de modèles empiriques d'acidification ont permis de dégager les zones les plus vulnérables à l'acidification. L'application de ces techniques d'analyse a montré, entre autres, que les concentrations de sulfates en lac étaient beaucoup plus élevées à proximité de Rouyn‐Noranda où est située l'usine Mines Noranda Itée. De fait, la concentration de sulfates dans les lacs de cette région passe de 16 mg L^?1, aux alentours de Rouyn‐Noranda, à moins de 2 mg L^?1 en direction du nord‐est, alors qu'elle baisse à 8 mg L^?1 en bordure de la rivière des Outaouais. De plus, en faisant abstraction de la distribution en forme de cloche des concentrations de sulfates autour de Rouyn‐Noranda, l'étude a permis de faire ressortir la décroissance générale des sulfates du sud‐ouest vers le nord‐est.

Les points majeurs qui ressortent de cette étude sont la présence de lacs très acides à proximité de Rouyn‐Noranda et V extrême vulnérabilité des plans d'eau situés au sud et sud‐est de la ville. En effet, l'eau des lacs de la région du Témiscamingue est en général moins minéralisée et plus acide que celle des lacs au nord de cette région. Cette distinction est due à la nature géologique des substrats dans les bassins versants. Ainsi, les lacs sensibles à l'acidification sont principalement situés sur des bassins versants associés à des roches granitiques et gneissiques, comparativement aux roches volcaniques et aux dépôts d'argiles lacustres qui caractérisent les lacs aux eaux plus minéralisées.     

La rationalisation du réseau météorologique du Québec: Stratégie d'intervention et méthodes d'analyse des données

Resume

Le Service de météorologie du ministère québécois de l'Environnement a voulu savoir comment ses réseaux météorologiques devraient évoluer pour répondre le plus adéquatement possible aux besoins réels, actuels et futurs du Québec.

C'est dans cette perspective que l'on définit une stratégie d'intervention caractérisée par (a) une conception dynamique des réseaux, (b) une régionalisation des besoins et des caractéristiques des variables, (c) des critères réalistes, (d) une étude du réseau région par région, (e) une planification visant à intégrer les réseaux entre eux, (f) une considération des méthodes d'analyse ou de simulation selon une approche régionale, (g) une durée et une densité d'observation adaptées aux besoins et, enfin, (h) une intégration rationnelle des diverses sources de données météorologiques. L'étude proprement dite des réseaux, selon une approche régionale, a été réalisée en faisant appel à l'analyse des données en composantes principales et à l'interpolation optimale. Les notions théoriques particulières à chacune de ces deux méthodes d'analyse sont présentées. On explique, par la suite, avec exemples à l'appui, comment ces méthodes sont appliquées, dans le cadre de la rationalisation du réseau météorologique du Québec, pour la détermination des régions statistiquement homogènes, l'estimation des fonctions de structure et des écarts types d'interpolation régionaux, en fonction de la distance entre les stations, et la cartographie des écarts types réels du réseau en opération, à une date donnée. On indique, aussi, comment les modifications projetées à ce réseau peuvent être évaluées.     

Simulation de l'énergie solaire au sol

Résumé

Le régime journalier du rayonnement solaire global au sol a été reconstitué aux six stations du réseau radiométrique de base canadien situées au Québec à l'aide d'un modèle numérique simulant le transfert de l'énergie solaire à travers l'atmosphère. Les données d'humidité atmosphérique, de trouble atmosphérique et d'albédo au sol ont été interpolées dans le temps et l'espace à partir de la climatologie mensuelle connue de ces champs. La nébulosité a été tirée des séquences horaires des durées d'ensoleillement observées à ces mêmes stations et échelonnées sur une même période de 974 jours. On a estimé les limites de précision du modèle en déterminant l'erreur quadratique moyenne entre la valeur calculée et mesurée pour 1, 2, 3, … et 20 jours consécutifs. On a ainsi trouvé que l'erreur quadratique moyenne diminue rapidement pour les regroupements de 1 à 10 jours et se stabilise pour les périodes plus longues. L'amplitude de l'erreur quadratique moyenne varie saisonnièrement et latitudinalement étant en général plus faible en été et au sud. La variabilité des valeurs journalières simulées n'est que légèrement supérieure à celle des données journalières mesurées. Il est démontré finalement que l'utilisation journalière des résultats du modèle est préférable à l'utilisation de la climatologie du régime solaire.     

La mesure de la durée d'insolation: Notes de recherche

A catastrophic rain event occurred in early June 2002 during a major drought over the Canadian Prairies. The storm brought record-breaking rainfall and major flooding to many locations in the region. Given the importance of this event, this study's overall objectives are to characterize and to understand the physical nature of the rainstorm better. The event was associated with a major extratropical cyclone which acted in concert with the Great Plains low-level jet to transport a tremendous amount of moisture into the eastern Prairies producing intense diurnal mesoscale convective systems over the region. At the same time, moisture was transported to the western Prairies by a strong easterly low-level jet which produced heavy and long-lived orographic precipitation near the foothills. Several working hypotheses were developed to explain the severity and longevity of the rainstorm; it was found that the Rockies played a central role in the organization and development of the system.

Atmospheric features that are critical to the development of an important class of extreme rain events in the Canadian Prairies were also identified. The severity of the June 2002 system is partly a result of the rare co-occurrence of these features during the period. Results from a preliminary analysis show that the atmospheric conditions associated with the extreme background drought enhanced the likelihood of the co-occurrence of these features during spring 2002, hence facilitating the development of the extreme rain event. In return, the tremendous precipitation from the storm alleviated the drought conditions in the southern Prairies.

R ésumé ?[Traduit par la rédaction] Un événement de pluie catastrophique s'est produit au début de juin 2002 pendant une sécheresse importante dans les Prairies canadiennes. La tempête a produit des chutes de pluie records et des inondations importantes à plusieurs endroits dans la région. Étant donné l'importance de cet événement, les objectifs généraux de la présente étude sont de mieux caractériser et de mieux comprendre la nature physique de la tempête de pluie. L'événement a été causé par une importante dépression extratropicale qui s'est associée au courant-jet à basse altitude des Grandes Plaines pour transporter une énorme quantité d'humidité dans l'est des Prairies et produire d'intenses systèmes convectifs de mésoéchelle diurne dans la région. En même temps, de l'humidité transportée dans l'ouest des Prairies par un fort courant-jet à basse altitude soufflant de l'est a produit d'intenses précipitations orographiques persistantes près des contreforts. Nous avons formulé plusieurs hypothèses de travail pour expliquer l'intensité et la persistance de la tempête de pluie; il ressort que les Rocheuses ont joué un rôle central dans l'organisation et le développement du système.

Nous avons aussi identifié les caractéristiques atmosphériques essentielles à la formation d'une catégorie importante d'événements de pluie extrêmes dans les Prairies canadiennes. L'intensité du système de juin 2002 est en partie le résultat de la rare présence simultanée de ces caractéristiques durant la période. Les résultats d'une analyse préliminaire montrent que les conditions atmosphériques liées à la sécheresse extrême qui régnait en toile de fond ont augmenté la probabilité d'une présence simultanée de ces caractéristiques au cours du printemps 2002, ce qui a favorisé la formation de l'événement de pluie extrême. En revanche, les précipitations diluviennes produites par la tempête ont atténué les conditions de sécheresse dans le sud des Prairies.     

Streamflow Modelling: A Primer on Applications,Approaches and Challenges

This article examines the current practice of streamflow modelling, a field under development for over a century. A sample of the wide range of assessment and planning applications of streamflow models is presented. The diversity in the use of these models is mirrored in the diversity of model complexity, and modelling approaches ranging from empirical to physically based and from lumped to fully distributed are described with examples. Predictions derived from hydrological models are subject to many sources of error; these are discussed along with methods for error minimization or anticipation. Model error is generally quantified using an ensemble of forecasts meant to sample the range of predictive uncertainty. This ensemble can be used to generate reliable probabilistic forecasts of hydrological quantities if all sources of error are accounted for. To date, applications of ensemble methods in streamflow forecasting have typically focused on only one or two error sources. A challenge will be to develop ensemble streamflow forecasts that sample a wider range of predictive uncertainty.

[Traduit par la rédaction] Le présent article examine la pratique actuelle en modélisation d’écoulement fluvial, un domaine qui évolue depuis plus d'un siècle. Nous présentons un échantillon de la vaste gamme d'applications d’évaluation et de planification des modèles d’écoulement fluvial. La diversité dans l'utilisation de ces modèles est le reflet de la diversité dans la complexité des modèles, et nous décrivons à l'aide d'exemples les approches de modélisation qui peuvent être empiriques ou basées sur la physique ou encore localisées ou entièrement réparties. Plusieurs sources d'erreur peuvent affecter les prévisions issues des modèles hydrologiques; nous discutons de ces sources d'erreur de même que des méthodes de réduction ou d'anticipation des erreurs. L'erreur du modèle est généralement quantifiée à l'aide d'un ensemble de prévisions servant à échantillonner la grandeur de l'incertitude prévisionnelle. Cet ensemble peut servir à produire des prévisions probabilistes fiables des grandeurs hydrologiques si toutes les sources d'erreur sont prises en compte. Jusqu’à maintenant, les applications des méthodes d'ensemble à la prévision des écoulements fluviaux n'ont généralement tenu compte que d'une ou deux sources d'erreur. Ce sera un défi de mettre au point des prévisions d'ensemble d’écoulement fluvial qui échantillonnent un plus large éventail d'incertitude prévisionnelle.     

Evidence of Change in the Winter Mixed Layer in the Northeast Pacific Ocean: A Problem Revisited

This paper presents a reanalysis of trends in water properties around the Gulf of Alaska. The datasets used are the lighthouse time series and the Line-P time series augmented with data from the Argo global float array. It is shown that though trends have changed, the essential story that was told 15 years ago has not changed in any meaningful way. Sea-surface temperatures are rising over a large part of the Gulf of Alaska, and sea-surface salinities are declining. Both of these lead to a decrease in surface density. We show that the temperature change extends deep into the Gulf of Alaska, and the freshening trend extends to a depth of about 100 db. Potential energy in the water column is decreasing, and the result of the decreasing potential energy and the increasing stratification implies a decrease in mid-winter mixed layer depths.

RÉSUMÉ?[Traduit par la rédaction] Cet article présente une réanalyse des tendances des propriétés de l'eau dans le golfe d'Alaska. Les ensembles de données utilisés sont les séries chronologiques des phares et les séries chronologiques de la ligne P augmentés des données du réseau mondial de flotteurs Argo. Il ressort que même si les tendances ont changé, à peu près tout ce que l'on a pu dire il y a 15 ans demeure valable. Les températures de la surface de la mer augmentent dans une grande partie du golfe d'Alaska et les salinités de la surface diminuent. Ces deux facteurs entraînent une diminution de la densité de la surface. Nous montrons que le changement de température se produit jusqu’à une bonne profondeur dans le golfe d'Alaska et que la tendance à l'adoucissement s’étend jusqu’à une profondeur d'environ 100 db. L’énergie potentielle dans la colonne d'eau diminue et cette diminution, de pair avec l'accroissement de la stratification, implique une diminution des profondeurs de la couche de mélange au milieu de l'hiver.     

Circulation Around Luzon Strait in September as Inferred from CTD,Argos and Argo Measurements and a Generalized Topography-Following Ocean Model

The hydrography, wind, Argos and Argo measurements in the areas surrounding Luzon Strait were collected. Based on the hydrographic data obtained during September 1994, the improved Princeton Ocean Model using a generalized topography-following coordinate system together with a modified inverse method was applied to study the circulation in September. Observations and the diagnostic simulation produce a consistent circulation pattern, and the main dynamical features can be summarized as follows. (1) The Kuroshio has two branches with the main Kuroshio existing above 800 m depth and the western part existing above 400 m depth. The western branch of the Kuroshio leaves the main current near 20.5°N, then flows northwestward through Luzon Strait and finally flows out of the northern boundary southwest of Taiwan, consistent with the trajectory of Argos drifters. (2) The non-linear term is important and cannot be neglected in the momentum equations in the northern part of Luzon Strait under the baroclinicity field. Using non-linear dynamics, the westward intrusion of the Kuroshio into the northern part of Luzon Strait is more curved than when using linear dynamics. However, the non-linear term is smaller and so negligible around Luzon Strait under the homogeneous density field. (3) In the area from 117°E to 119°E and from 20.2°N to 21.7°N, an anticyclonic eddy appears east of Dongsha Islands. (4) At depths above 400 m, the circulation is mainly dominated by the basin-scale cyclonic gyre. (5) In the computational domain west of 121°E, the circulation below 800 m is mainly dominated by the basin-scale anticyclonic gyre. (6) The South China Sea water flows eastward across Luzon Strait in the middle layers, then turns cyclonically, finally flowing northward into the region east of Taiwan Island, which is qualitatively in agreement with the trajectories of Argo floats at about 1000 m depth in the area east of 121°E.

RÉSUMÉ?[Traduit par la rédaction] Nous avons rassemblé l'hydrographie, le vent, les Argo et les mesures Argo dans les régions avoisinant le détroit de Luçon. Sur la base des données hydrographiques relevées durant le mois de septembre 1994, nous avons appliqué le modèle océanique amélioré de Princeton utilisant un système de coordonnées « suivant le terrain » généralisé de pair avec une méthode inverse modifiée pour étudier la circulation en septembre. Les observations et la simulation diagnostique produisent une configuration de circulation cohérente dont on peut résumer comme suit les principales caractéristiques dynamiques. (1) Le Kuroshio possède deux branches : le Kuroshio principal existant au-dessus d'une profondeur de 800 m et la partie ouest existant au-dessus d'une profondeur de 400 m. La branche ouest du Kuroshio quitte le courant principal aux environs de 20,5°N et se dirige vers le nord-ouest à travers le détroit de Luçon pour finalement franchir la limite nord du détroit au sud-ouest de Taïwan, ce qui correspond à la trajectoire des bouées dérivantes Argos. (2) Le terme non linéaire dans les équations est important et ne peut pas être considéré négligeable de quantité de mouvement dans la partie nord du détroit de Luçon dans le champ de baroclinité. Lorsqu'on utilise la dynamique non linéaire, l'intrusion vers l'ouest du Kuroshio dans la partie nord du détroit de Luçon est plus incurvée que lorsqu'on utilise la dynamique linéaire. Cependant, le terme non linéaire est plus petit et donc négligeable autour du détroit de Luçon dans le champ de densité homogène. (3) Dans la région allant de 117°E à 119°E et de 20,2°N à 21,7°N, un remous anticyclonique apparaît à l'est des îles Dongsha. (4) Aux profondeurs inférieures à 400 m, la circulation est principalement dominée par le gyre cyclonique à l’échelle du bassin. (5) Dans le domaine de calcul à l'ouest de 121°E, la circulation sous 800 m est principalement dominée par le gyre anticyclonique à l’échelle du bassin. (6) Les eaux de la mer de Chine méridionale s’écoulent vers l'est à travers le détroit de Luçon dans les couches intermédiaires, puis tournent de façon cyclonique pour finalement s’écouler vers le nord dans la région à l'est de l’île Taïwan, ce qui, qualitativement, s'accorde avec la trajectoire des flotteurs Argo à environ 1000 m de profondeur dans la région à l'est de 121°E.     

Extreme Waves and Wave Run-up in Halifax Harbour under Climate Change Scenarios

The tendency for rising sea levels, combined with changes in the frequency and intensity of extreme storm events raises the potential for flooding and inundation of coastal locations such as Halifax Harbour, in the context of climate change. In this study, we consider three scenarios for extreme high water levels based on climate change scenarios and estimates for land subsidence, rising mean sea levels, and return period analysis for extreme events (from previous studies). We also investigate the effect of ocean waves on these estimates for extreme high water levels. Because the most damaging storm to make landfall in Nova Scotia over the last century was Hurricane Juan (2003), it was chosen to simulate the extreme case of storm-generated waves and wave run-up in Halifax Harbour. To simulate waves generated by Hurricane Juan, a nested-grid system consisting of two modern state-of-the-art operational wave models was used. High quality winds were used to drive the wave models, and the simulations used recently updated high resolution coastal bathymetry. Observed water elevation changes in Halifax Harbour were used in wave model simulations of Hurricane Juan. The Federal Emergency Management Agency's (FEMA) run-up model is used to estimate wave run-up elevation, which is validated with recorded high water observations along the coastline. Simulated waves and wave run-up elevations for Hurricane Juan suggest that the maximum significant wave heights at the mouth of the Harbour were 9.0 m, and the wave run-up was as high as 2.0 m along the shoreline of Halifax Harbour. In this way, we estimated the impact of waves and wave run-up on extreme high water elevations for three climate change scenarios in Halifax Harbour, under worst-case conditions. The sensitivity of these estimates is analyzed for different water level variations, wave propagation directions and shore slope profiles.

[Traduit par la rédaction] La tendance à la hausse du niveau de la mer de pair avec les changements dans la fréquence et l'intensité des événements de tempêtes extrêmes augmentent le risque d'inondation et de submersion à des emplacements côtiers comme le port d'Halifax dans le contexte du changement climatique. Dans cette étude, nous examinons trois scénarios de niveaux de hautes eaux extrêmes basés sur des scénarios de changement climatique et estimations de subsidence du terrain, l’élévation du niveau moyen de la mer ainsi que l'analyse de la période de retour d’événements extrêmes (faite lors d’études antérieures). Nous examinons aussi les effets des vagues de l'océan sur ces estimations de niveaux de hautes eaux extrêmes. Étant donné que la tempête la plus dévastatrice à avoir touché terre en Nouvelle–Écosse au cours du dernier siècle a été l'ouragan Juan (2003), c'est celle-ci que nous avons choisie pour simuler le cas extrême de vagues produites par une tempête et de remontée de vagues dans le port d'Halifax. Pour simuler les vagues produites par Juan, nous nous sommes servis d'un système à grilles imbriquées consistant en deux modèles de vagues opérationnels à la fine pointe de la technologie. Les modèles de vagues étaient pilotés par des données de vent de haute qualité et les simulations disposaient d'une bathymétrie côtière à haute résolution récemment mise à jour. Les changements observés d’élévation de l'eau dans le port d'Halifax ont été utilisés dans les simulations de l'ouragan Juan par les modèles de vagues. Nous utilisons le modèle de remontée des vagues de la Federal Emergency Management Agency (FEMA) pour estimer la hauteur des remontées, qui est validée par rapport aux observations de hautes eaux enregistrées le long de la côte. Les vagues et les hauteurs de remontée des vagues simulées pour l'ouragan Juan suggèrent que la hauteur significative maximale des vagues à l'entrée du port était de 9,0 m et la remontée des vagues atteignait 2,0 m le long de la côte du port d'Halifax. Nous avons de cette manière estimé l'impact des vagues et de la remontée des vagues sur l’élévation des hautes eaux extrêmes pour trois scénarios de changement climatique dans le port d'Halifax dans les conditions les plus défavorables. Nous analysons la sensibilité de ces estimations à divers changements de niveau d'eau, directions de propagation et profils de pente côtière.     

A Time-Averaged Covariance Method in the EnKF for Argo Data Assimilation

In the ensemble Kalman filter (EnKF), ensemble size is one of the key factors that significantly affects the performance of a data assimilation system. A relatively small ensemble size often must be chosen because of the limitations of computational resources, which often biases the estimation of the background error covariance matrix. This is an issue of particular concern in Argo data assimilation, where the most complex state-of-the-art models are often used. In this study, we propose a time-averaged covariance method to estimate the background error covariance matrix. This method assumes that the statistical properties of the background errors do not change significantly at neighbouring analysis steps during a short time window, allowing the ensembles generated at previous steps to be used in present steps. As such, a joint ensemble matrix combining ensembles of previous and present steps can be constructed to form a larger ensemble for estimating the background error covariance. This method can enlarge the ensemble size without increasing the number of model integrations, and this method is equivalent to estimating the background error covariance matrix using the mean ensemble covariance averaged over several assimilation steps. We apply this method to the assimilation of Argo and altimetry datasets with an oceanic general circulation model.

Experiments show that the use of this time-averaged covariance can improve the performance of the EnKF by reducing the root mean square error (RMSE) and improving the estimation of error covariance structure as well as the relationship between ensemble spread and RMSE.

RÉSUMÉ [Traduit par la rédaction] Dans le filtre de Kalman d'ensemble (EnKF), la taille de l'ensemble est l'un des facteurs clés qui ont une influence importante sur la performance d'un système d'assimilation de données. Il faut souvent choisir une taille d'ensemble assez petite à cause des limites des ressources informatiques, ce qui biaise souvent l'estimation de la matrice de covariance de l'erreur de fond. Cette question revêt une importance particulière pour l'assimilation des données Argo, qui fait souvent appel à des modèles de pointe très complexes. Dans cette étude, nous proposons une méthode de covariance moyennée dans le temps pour estimer la matrice de covariance de l'erreur de fond. Cette méthode suppose que les propriétés statistiques des erreurs de fond ne changent pas de façon importante d'une étape d'analyse à la suivante durant un court laps de temps, ce qui permet d'utiliser dans les étapes courantes les ensembles générés aux étapes précédentes. Ainsi, on peut construire une matrice d'ensembles conjoints combinant les ensembles des étapes précédentes et courantes pour former un plus grand ensemble dans le but d'estimer la covariance de l'erreur de fond. Cette méthode peut accroître la taille de l'ensemble sans augmenter le nombre d'intégrations du modèle; elle équivaut à estimer la matrice de covariance de l'erreur de fond en utilisant la covariance moyenne de l'ensemble calculée sur plusieurs étapes d'assimilation. Nous appliquons cette méthode à l'assimilation des ensembles de données Argo et d'altimétrie avec un modèle de circulation océanique générale.

Des essais montrent que l'emploi de cette covariance moyennée dans le temps peut améliorer la performance de l'EnKF en réduisant l’écart-type et en améliorant l'estimation de la structure de la covariance de l'erreur de même que la relation entre l'étalement et l'écart-type l'ensemble.     

Assessing the Application of Argo Profiling Float Data to the Study of the Seasonal Variation of the Hydrological Parameters and the Current Field East of Luzon Strait

The seasonal characteristics of water masses, mixed-layer depth and the current field east of Luzon Strait are studied using Argo profiling float data. Based on data from March 2006 to November 2010, the temperature-salinity relation indicates that the seasonal variation of water masses is not obvious, except that the surface temperature is lower and the salinity is higher in spring and winter; the water masses in other layers change little. The seasonal variation of the mixed-layer depth, which is deepest at over 170 m in winter and shallowest at about 25 m in summer, is also discussed. Based on the P-vector method, multi-year seasonal mean Argo data and Levitus data are used to calculate the respective current fields. The surface geostrophic current derived from altimetric data is compared with those computed from Argo data and Levitus data. It is shown that, the current field computed from Argo data, which shows the Kuroshio and eddies clearly, is similar to that obtained from altimetric data, except that the current speed is less in the former case than in the latter one; this may be related to the scarcity and uneven distribution of Argo floats, the subsequent interpolation and extrapolation errors, the shortcomings of the P-vector method, the lack of a barotropic component in the calculation or the altimetry-derived current error associated with the marine geoid, but the resultant current field is much better than that obtained from Levitus data. Moreover, Argo data have an advantage in that they can be used to derive a three-dimensional current field, whereas the altimetric data can only be used to derive the surface current field. The seasonal variations in the vertical current structure and the Kuroshio transport east of Taiwan Island are also discussed; it is found that the thickness of the northward Kuroshio is about 700 m, and the seasonal zonal variation in the main axis of the Kuroshio is obvious.

RÉSUMÉ?[Traduit par la rédaction] Nous étudions les caractéristiques saisonnières des masses d'eau, la profondeur de la couche de mélange et le champ de courant à l'est du détroit de Luçon à l'aide des données des flotteurs profileurs Argo. D'après les données recueillies entre mars 2006 et novembre 2010, la relation température-salinité indique que la variation saisonnière des masses d'eau n'est pas évidente, sauf que la température de surface est plus basse et la salinité plus élevée au printemps et en hiver; les masses d'eau dans les autres couches changent peu. Nous discutons aussi de la variation saisonnière de la profondeur de la couche de mélange, qui est maximale à plus de 170 m en hiver et minimale à environ 25 m en été. En nous basant sur la méthode des vecteurs P, nous utilisons les données Argo et les données Levitus pluriannuelles de moyennes saisonnières pour calculer les champs de courant respectifs. Nous comparons le courant géostrophique de surface déduit des données altimétriques avec ceux calculés à l'aide des données Argo et des données Levitus. Il apparait que le champ de courant calculé d'après les données Argo, qui montre clairement le Kuroshio et les remous, est semblable à celui obtenu des données altimétriques, sauf que la vitesse du courant est moindre dans le premier cas que dans le dernier; cela peut être lié à la rareté et à la distribution inégale des flotteurs Argo, aux erreurs subséquentes d'interpolation et d'extrapolation, aux faiblesses de la méthode des vecteurs P, au manque d'une composante barotropique dans les calculs ou à l'erreur dans le courant déduit des données altimétriques associée au géoïde marin, mais le champ de courant résultant est bien meilleur que celui obtenu à l'aide des données Levitus. De plus, les données Argo présentent l'avantage de pouvoir être utilisées pour déduire un champ de courant tridimensionnel, alors que les données altimétriques ne permettent que de déduire le champ de courant de surface. Nous discutons aussi des variations saisonnières dans la structure verticale du courant et dans le transport par le Kuroshio à l'est de l’île de Taïwan; nous trouvons que la profondeur du Kuroshio circulant vers le nord est d'environ 700 m et la variation saisonnière zonale dans l'axe principal du Kuroshio est évidente.     

Drought and Associated Cloud Fields over the Canadian Prairie Provinces

Little is known about clouds during drought. From 1999 to 2005 the Canadian Prairies experienced one of the most severe and prolonged droughts in the historical record. This study characterizes clouds during drought in the Canadian Prairie provinces with a particular focus on this recent drought.

Drought severity was determined using the Standardized Precipitation Index (SPI) based on monthly precipitation on a 1° × 1° grid. Cloud fields from the National Aeronautics and Space Administration/Global Energy and Water Experiment's (NASA/GEWEX) Surface Radiation Budget database were used to examine overall cloud amount, optical thickness, and top-of-the-atmosphere albedo. Anomalies in monthly precipitation in the satellite record from 1984 to 2004, with an emphasis on the recent drought from 1999 to 2004, were related to anomalies in cloud fields.

During drought, a decrease in cloud amount was observed. During the spring and summer months of the 1999–2004 drought, for example, the observed cloud cover fraction decreased by approximately 7% between severely wet and severely dry conditions. There was, however, large month-to-month and spatial variability, and the correlation of cloud cover fraction anomaly with precipitation was weak. A higher correlation was found between the top-of-the-atmosphere albedo and precipitation. The occurrence of thick clouds and clouds of medium thickness did decrease with drought severity. These trends also applied to subregions within the overall domain. These observations further the understanding of the role of clouds in feedback mechanisms during drought.

R ésumé ?[Traduit par la rédaction] On sait peu de choses à propos des nuages durant les sécheresses. De 1999 à 2005, les Prairies canadiennes ont connu l'une des périodes de sécheresse les plus longues et les plus intenses depuis que l'on tient des relevés. La présente étude porte sur les nuages durant les périodes de sécheresse dans les Prairies canadiennes, en mettant l'accent sur les sécheresses récentes.

Nous avons déterminé l'intensité des sécheresses au moyen de l'indice de précipitations normalisé (SPI) fondé sur les précipitations mensuelles sur une grille de 1° × 1°. Nous avons étudié les champs de nuages obtenus de la base de données du bilan radiatif de la surface du GEWEX (Expérience mondiale sur les cycles de l'énergie et de l'eau) de la NASA pour examiner l'étendue générale des nuages, l'épaisseur optique et l'albédo du sommet de l'atmosphère. Nous avons mis en relation les anomalies dans les précipitations mensuelles d'après les données satellitaires de 1984 à 2004, en mettant l'accent sur la sécheresse récente de 1999 à 2004, avec les anomalies dans les champs de nuages.

Pendant les sécheresses, nous avons noté une diminution de l'étendue des nuages. Durant les mois de printemps et d'été de la période de sécheresse de 1999 à 2004, par exemple, la fraction de couverture du ciel observée a diminué d'environ 7 % entre les conditions très humides et très sèches. Cependant, la variabilité intermensuelle et spatiale est grande, et la corrélation de l'anomalie de la fraction de couverture nuageuse avec les précipitations est faible. Nous avons trouvé une meilleure corrélation entre l'albédo du sommet de l'atmosphère et les précipitations. La quantité de nuages épais et de nuages d'épaisseur moyenne diminuait quand l'intensité de la sécheresse augmentait. Ces tendances s'observent aussi dans les sous-régions à l'intérieur domaine général. Ces observations permettent de mieux comprendre le rôle des nuages dans les mécanismes de rétroaction au cours des sécheresses.     

Étude de la circulation verticale dans l'estuaire du Saint‐Laurent au moyen de la modélisation mathématique

Résumé

Un modèle bidimensionnel dans le sens vertical est utilisé afin d'étudier et de simuler la circulation marine engendrée par la marée dans l'estuaire du Saint‐Laurent. A partir d'une schématisation de l'Estuaire et de conditions en niveaux d'eau et salinités aux extrémités aval et amont du domaine d'étude, ce modèle génère la répartition des hauteurs d'eau et les distributions longitudinales et verticales de vitesses et salinités. Les résultats permettent alors de déduire les circulations instantanées obtenues aux différentes heures de la marée. Celles‐ci mettent en évidence des phénomènes physiques caractéristiques de l'estuaire du Saint‐Laurent et qui ont été plusieurs fois observés. Parmi ceux‐ci, les processus les plus importants concernent les marées internes, la remontée d'eau profonde et les effets des gradients horizontaux de salinités sur les profils verticaux de vitesse.     

Current Observations of Internal Tides and Parametric Subharmonic Instability in Luzon Strait

Acoustic Doppler Current Profiler (ADCP) measurements of the velocity structure in the thermocline in Luzon Strait are presented. The statistics for current, vertical variation of the inertia-gravity waveband, parametric subharmonic instability (PSI), and current shear are analyzed. It was found that 1) barotropic flow primarily consists of a nearly circular mixed tide. Diurnal tides are strongest and show smooth variance with a fortnightly spring-neap cycle, indicative of the astronomical tide-generating force. However, the semi-diurnal band power exhibits a high-frequency oscillation as a result of non-linear interactions. The high-frequency band power with high values during the spring tide oscillates with the tidal cycle. Near-inertial wave motions showing random variance may be caused by changes in the wind forcing at the sea surface or by random forcing. 2) Baroclinic velocities exhibit strong shear structure. The observed large changes in the amplitude of the baroclinic velocity and the limited vertical extent of the high-velocity cores may be interpreted as internal wave beams that pass through the observed water column. Semi-diurnal tides are dominant in the baroclinic velocity. Kinetic energy spectra also revealed that additional peaks were centred at sum-tidal-inertial interaction frequencies (such as M₂ ?+?f) and difference-interaction frequencies (such as M₂ ???f). The spectral exponent of the baroclinic velocity is ω^?α (1?<?α?<?3). 3) Strong non-linear interactions among internal waves exist, and the semi-diurnal (M₂) component plays a key role in these interactions. Bicoherence analysis showed that M₂/2 waves were non-linearly coupled with the dominant M ₂ internal tide. 4) The polarization relations were used to diagnose observational internal tidal motions. Diurnal waves propagate to the east-northeast, and the semi-diurnal–diurnal waves propagate westward. In the case of diurnal tides, the minor to major axis ratio is different from the expected value of f/ω_{K 1} because of the deviation of inclinations, whereas, for semi-diurnal tides, it is close to the expected value of f/ω _M2 at depths from 30 to 150 m.

RÉSUMÉ?[Traduit par la rédaction] Nous présentons les mesures de la structure des vitesses dans la thermocline faites par profileur de courant à effet Doppler (ADCP) dans le détroit de Luçon. Nous analysons les statistiques sur le courant, la variation verticale de la gamme d'ondes d'inertie-gravité, l'instabilité subharmonique paramétrique (PSI) et le cisaillement du courant. Il ressort que 1) l’écoulement barotrope consiste principalement en une marée mixte presque circulaire. Les marées diurnes sont les plus fortes et présentent une variance régulière dans un cycle vives-eaux mortes-eaux de deux semaines qui révèle la nature astronomique de la force qui produit les marées. Cependant, la puissance de bande semi-diurne affiche une oscillation de haute fréquence causée par des interactions non linéaires. La puissance de bande de haute fréquence avec des valeurs élevées durant les vives-eaux oscille avec le cycle des marées. Les mouvements des ondes quasi inertielles affichant une variance aléatoire peuvent être dus à des variations dans le forçage par le vent à la surface de la mer ou à un forçage aléatoire. 2) Les vitesses baroclines exhibent une structure fortement cisaillée. Les changements marqués observés dans l'amplitude des vitesses baroclines et l'extension verticale limitée des c?urs de vitesse élevée peuvent être interprétés comme des faisceaux d'ondes internes qui traversent la colonne d'eau observée. Les marées semi-diurnes sont dominantes dans les vitesses baroclines. Les spectres d’énergie cinétique ont aussi révélé que des crêtes additionnelles étaient centrées aux fréquences d'interaction cumulative marée-inertie (comme à M₂ ?+?f) et aux fréquences d'interaction différentielle (comme à M₂ ???f). L'exposant spectral de la vitesse barocline est ω^?α (1?<?α?<?3). 3) Il existe de fortes interactions non linéaires entre les ondes internes et la semi-diurne (M₂) joue un rôle clé dans ces interactions. L'analyse de bicohérence a montré que les ondes M₂/2 étaient non-linéairement couplées avec la marée interne dominante M₂ . 4) Nous avons utilisé les relations de polarisation pour diagnostiquer les mouvements de marée internes déduits des observations. Les ondes diurnes se propagent vers l'est-nord-est et les ondes semi-diurnes se propagent vers l'ouest. Dans le cas des marées diurnes, le rapport de l'axe secondaire à l'axe principal diffère de la valeur attendue (f/ω_{K 1} ) à cause de l’écart des inclinaisons alors que pour les marées semi-diurnes, il est proche de la valeur attendue (f/ω _M2 ) à des profondeurs de 30 à 150 m.     

Lapse Rate Adjustments of Gridded Surface Temperature Normals in an Area of Complex Terrain: Atmospheric Reanalysis versus Statistical Up-Sampling

The applicability of elevation-regression based interpolation methods for long-term temperature normals, for example the Parameter-elevation Regressions on Independent Slopes Model (PRISM), becomes increasingly limited in data sparse, complex terrain such as that found in mountainous British Columbia (BC), Canada. Recent methods to improve both the resolution and accuracy of interpolation models have focused on the development of “up-sampling” algorithms based on local lapse rate adjustments to the original interpolated surfaces. Lapse rates can be derived from statistical models (e.g., elevation-based polynomial regression equations) or dynamical models (e.g., vertical temperature profiles from numerical weather prediction (NWP) models). This study compares a widely used statistical up-sampling algorithm, ClimateBC, with two NWP reanalysis products, the National Centers for Environmental Prediction/National Corporation for Atmospheric Research, Reanalysis 1 (NCEP1) and the more modern European Centre for Medium-range Weather Forecasts (ECMWF) Reanalysis Interim (ERA-Interim). Thirty-year climate normals for maximum and minimum temperatures were calculated using statistical up-sampling and NWP lapse rate adjustments to existing PRISM-based climate normals at a subset of stations in BC. Specifically, up-sampling model evaluation was performed using 1951–80 climate normals from an independent set of 54 surface stations (1 m to 2347 m) which were not included in the PRISM interpolation or assimilated into the NWP reanalysis products. All models performed similarly for minimum temperature, which showed only a slight improvement over PRISM. For maximum temperature, ClimateBC, NCEP1 and ERA-Interim all performed significantly better than PRISM, in particular during spring and summer. The ERA-Interim reanalysis outperformed NCEP1 in almost all months. The results suggest that lapse rate adjustment algorithms based on reanalysis products will have greater potential as progress continues on developing NWP components.

R ésumé ?[Traduit par la rédaction] L'application des techniques d'interpolation par régression en fonction de l'altitude pour les normales de température à long terme, comme le Parameter-elevation Regressions on Independent Slopes Model (PRISM), devient très difficile dans les régions accidentées pour lesquelles on dispose de données insuffisantes, par exemple les secteurs montagneux de la Colombie-Britannique (C.-B.) au Canada. Les toutes dernières méthodes destinées à augmenter le degré de résolution des modèles d'interpolation et leur précision reposent sur la conception d'algorithmes d’échantillonnage vertical fondés sur l'ajustement des surfaces interpolées originales au moyen du gradient vertical local. Nous pouvons établir les gradients verticaux à partir de modèles statistiques (p. ex., des équations de régression polynomiales en fonction de l'altitude) ou de modèles dynamiques (p. ex., des profils verticaux de température à partir de modèles de prévision numérique du temps (PNT)). Dans la présente étude, nous comparons un algorithme d’échantillonnage vertical statistique communément utilisé, le programme ClimateBC, à deux produits de réanalyse de PNT, celle des National Centres for Environmental Prediction/National Corporation for Atmospheric Research Reanalysis 1 (NCEP1), et la réanalyse provisoire (ERA-Interim) du Centre européen pour les prévisions météorologiques à moyen terme (ECMWF). Les normales climatiques de trente ans pour les températures maximums et minimums ont été calculées en appliquant la méthode d’échantillonnage vertical statistique et l'ajustement du gradient obtenu par PNT aux normales climatiques établies à partir du PRISM pour un sous-ensemble de stations en Colombie-Britannique. Plus particulièrement, nous avons procédé à l’évaluation du modèle d’échantillonnage vertical en nous servant des normales climatiques (1951–1980), pour un ensemble de 54 stations d'observation en surface indépendantes (1?m à 2347?m), exclues du modèle d'interpolation PRISM et des produits de réanalyse de PNT. Pour tous les modèles, nous avons obtenu des résultats comparables pour la température minimum, soit une légère amélioration seulement par rapport au PRISM. Pour la température maximum, nous avons obtenu avec ClimateBC, NCEP1 et ERA-Interim, des résultats nettement plus probants qu'avec PRISM, notamment au printemps et en été. Les réanalyses ERA-Interim ont donné de meilleurs résultats que NCEP1 pour pratiquement tous les mois. D'après ces résultats, le potentiel des algorithmes d'ajustements des gradients verticaux de température, établis à partir de produits de réanalyse se renforcera à mesure que les composantes de PNT se développeront.     

Downscaling of Precipitation over Vancouver Island using a Synoptic Typing Approach

A statistical downscaling technique is employed to link atmospheric circulation produced by an ensemble of global climate model (GCM) simulations over the twenty-first century to precipitation recorded at weather stations on Vancouver Island. Relationships between the different spatial scales are established with synoptic typing, coupled with non-homogeneous Markov models to simulate precipitation intensity and occurrence. Types are generated from daily precipitation observations spanning 1971 to 2000. Atmospheric predictors used to influence the Markov models are derived from two versions of GCM output: averages of GCM grid cells selected by correlation maps of circulation and precipitation data and an approach involving common Empirical Orthogonal Functions (EOFs) calculated from GCM output over the northeast Pacific Ocean. Projections for 2081 to 2100 made using averaged grid cells find that winter (November–February) precipitation anomalies produce modestly positive values, with gains of 7.5% in average precipitation, typical increases of 9.0% rising to 20% in the case of high-intensity precipitation, and little spatial dependence. In contrast, average and high-intensity summer precipitation (June–September) decline negligibly at most island weather stations with the exception of those in the southwestern sections, which experience reductions of 15% relative to 1971 to 2000. Projections made using common EOFs display a strong spatial dependence. Future winter precipitation is expected to increase only on the west coast of the island by 11%, on average, while the southeastern coast will experience decreases of 5% to 10%. The same pattern repeats in summer, though with negligible increases on the west coast and declines of 12% to 16% on the southeastern coast. The reliability of this novel EOF method remains to be confirmed definitively, however. In both seasons precipitation occurrence decreases slightly at all stations with declines in the total days with measurable precipitation ranging from 2% to 8%.

RÉSUMÉ [Traduit par la rédaction] Nous employons une technique statistique de réduction d’échelle pour lier la circulation atmosphérique produite par un ensemble de simulations du GCM (Global Climate Model) durant le XXI^e siècle aux précipitations enregistrées à des stations météorologiques sur l’île de Vancouver. Les relations entre les différentes échelles spatiales sont établies au moyen d'un typage synoptique couplé avec des modèles markoviens non homogènes pour simuler l'intensité et la fréquence des précipitations. Les types sont générés à partir des observations quotidiennes de précipitations au cours de la période 1971–2000. Les prédicteurs atmosphériques utilisés pour influencer les modèles markoviens sont dérivés de deux versions de sorties du GCM : les moyennes de mailles du GCM sélectionnées par tables de corrélation des données de circulation et de précipitations et une approche fondée sur les fonctions orthogonales empiriques (EOF) communes calculées d'après la sortie du GCM pour le nord-est du Pacifique. Les projections pour la période 2081–2100 basées sur des moyennes de mailles montrent que les anomalies de précipitations hivernales (novembre–février) produisent de faibles valeurs positives, avec des gains de 7.5% dans les précipitations moyennes, des accroissements caractéristiques de 9.0% augmentant à 20% dans le cas des précipitations de forte intensité, et peu de dépendance spatiale. En revanche, les précipitations estivales (juin–septembre) moyennes et de forte intensité diminuent de façon négligeable à la plupart des stations météorologiques de l’île, à l'exception de celles situées dans secteur sud-ouest qui subissent une réduction de 15% par rapport à 1971–2000. Les projections faites à l'aide des fonctions orthogonales empiriques communes exhibent une forte dépendance spatiale. Les précipitations hivernales futures devraient augmenter seulement sur la côte ouest de l’île de 11% en moyenne alors que la côte sud-est connaîtra des diminutions de 5 à 10%. La même configuration se répète en été, bien qu'avec des accroissements négligeables sur la côte ouest et des diminutions de 12 à 16% sur la côte sud-est. La fiabilité de cette nouvelle méthode EOF reste toutefois à établir. Dans les deux saisons, la fréquence des précipitations diminue légèrement à toutes les stations, les diminutions du nombre total de jours avec précipitations mesurables variant entre 2 et 8%.     

Considérations sur l'évaluation des bénéfices à attendre d'une opération de lutte contre la grêle

Résumé

Dans la mesure où l'on entrevoit la possibilité de réduire, sinon de sup‐primer les chutes de grêle dévastatrices, il est nécessaire d'établir un bilan coût /bénéjîce ajin dejuger l'intérêt de telles opérations. La connaissance des relations entre caractéristiques physiques des chutes de grêle et les dégâts qu'elles occasion‐nent à la végétation est une élément indispensable à ce calcul. La méthode analytique, étant donné la diversité des cas de figures, risque d'être longue à mettre en ?uvre. Il peut être intéressant de rapprocher les résultats d'une réseau de grêlimètres aux constats des Compagnies d'assurances pour déterminer des seuils de façon statistique.     

Chaleur Urbaine à Montréal

Par des sondages de température effectués en hélicoptère à 12 sites de l'île de Montréal, on a pu mesurer la chaleur générée par la ville en 37 jours assez bien distribués à travers l'année. Cette chaleur fut mise en relation avec deux fonctions exprimant l'intensité du chauffage, la première, selon le nombre de degrés sous 65°F de l'air ambiant, et la seconde fonction, selon la même variable multipliée par la vitesse du vent. Une relation linéaire fut établie dans l'un et l'autre cas entre la chaleur urbaine et l'intensité du chauffage, mais avec une dispersion minimale dam le second cas. Un calcul à partir de la meilleure droite de régression permet de déterminer la chaleur totale génétrée par la ville dans une année, soit 40.2×10¹⁶ calories‐gramme Cette quantité de chaleur est de beaucoup supérieure à la chaleur artificielle générée par les combustions, soit 8.40 × 10¹⁶ calories. On en conclut que la chaleur urbaine ne peut être expliquée sans faire appel au rayonnement so‐Iaire.