Search

An unprecedented heat wave occurred over the Pacific Northwest and southwest Canada on 25–30 June 2021, resulting in all-time temperature records that greatly exceeded previous record maximum temperatures. The impacts were substantial, including several hundred deaths, thousands of hospitalizations, a major wildfire in Lytton, British Columbia, Canada, and severe damage to regional vegetation. Several factors came together to produce this extreme event: a record-breaking midtropospheric ridge over British Columbia in the optimal location, record-breaking midtropospheric temperatures, strong subsidence in the lower atmosphere, low-level easterly flow that produced downslope warming on regional terrain and the removal of cooler marine air, an approaching low-level trough that enhanced downslope flow, the occurrence at a time of maximum insolation, and drier-than-normal soil moisture. It is shown that all-time-record temperatures have not become more frequent and that annual high temperatures only increased at the rate of baseline global warming. Although anthropogenic warming may have contributed as much as 1°C to the event, there is little evidence of further amplification from increasing greenhouse gases. Weather forecasts were excellent for this event, with highly accurate predictions of the extreme temperatures. Significance Statement: This paper describes the atmospheric evolution that produced an extreme heat wave over the Pacific Northwest during June 2021 and puts this event into historical perspective. [ABSTRACT FROM AUTHOR]

Accurately modeling the interactions between inland water bodies and the atmosphere in meteorological and climate models is crucial, given the marked differences with surrounding landmasses. Modeling surface heat fluxes remains a challenge because direct observations available for validation are rare, especially at high latitudes. This study presents a detailed evaluation of the Canadian Small Lake Model (CSLM), a one-dimensional mixed-layer dynamic lake model, in reproducing the surface energy budget and the thermal stratification of a subarctic reservoir in eastern Canada. The analysis is supported by multiyear direct observations of turbulent heat fluxes collected on and around the 85-km2 Romaine-2 hydropower reservoir (50.7°N, 63.2°W) by two flux towers: one operating year-round on the shore and one on a raft during ice-free conditions. The CSLM, which simulates the thermal regime of the water body including ice formation and snow physics, is run in offline mode and forced by local weather observations from 25 June 2018 to 8 June 2021. Comparisons between observations and simulations confirm that CSLM can reasonably reproduce the turbulent heat fluxes and the temperature behavior of the reservoir, despite the one-dimensional nature of the model that cannot account for energy inputs and outputs associated with reservoir operations. The best performance is achieved during the first few months after the ice break-up (mean error = −0.3 and −2.7 W m−2 for latent and sensible heat fluxes, respectively). The model overreacts to strong wind events, leading to subsequent poor estimates of water temperature and eventually to an early freeze-up. The model overestimated the measured annual evaporation corrected for the lack of energy balance closure by 5% and 16% in 2019 and 2020. Significance Statement: Freshwater bodies impact the regional climate through energy and water exchanges with the atmosphere. It is challenging to model surface energy fluxes over a northern lake due to the succession of stratification and mixing periods over a year. This study focuses on the interactions between the atmosphere of an irregular shaped northern hydropower reservoir. Direct measurements of turbulent fluxes using an eddy covariance system allowed the model assessment. Turbulent fluxes were successfully predicted during the open water period. Comparison between observed and modeled time series showed a good agreement; however, the model overreacted to high wind episodes. Biases mostly occur during freeze-up and breakup, stressing the importance of a good representation of the ice cover processes. [ABSTRACT FROM AUTHOR]

The purpose of this study is to present the impacts of a fully interactive coupling between an atmospheric and a sea ice model over the Gulf of St. Lawrence, Canada. The impacts are assessed in terms of the atmospheric and sea ice forecasts produced by the coupled numerical system. The ocean-ice model has been developed at the Maurice Lamontagne Institute, where it runs operationally at a horizontal resolution of 5 km and is driven (one-way coupling) by atmospheric model forecasts provided by the Meteorological Service of Canada (MSC). In this paper the importance of two-way coupling is assessed by comparing the one-way coupled version with a two-way coupled version in which the atmospheric model interacts with the sea ice model during the simulation. The impacts are examined for a case in which the sea ice conditions are changing rapidly. Two atmospheric model configurations have been studied. The first one has a horizontal grid spacing of 24 km, which is the operational configuration used at the Canadian Meteorological Centre. The second one is a high-resolution configuration with a 4-km horizontal grid spacing. A 48-h forecast has been validated using satellite images for the ice and the clouds, and also using the air temperature and precipitation observations. It is shown that the two-way coupled system improves the atmospheric forecast and has a direct impact on the sea ice forecast. It is also found that forecasts are improved with a fine resolution that better resolves the physical events, fluxes, and forcing. The coupling technique is also briefly described and discussed. [ABSTRACT FROM AUTHOR]

The performance of the Weather Research and Forecasting (WRF) Model is evaluated in predicting the meteorological conditions over a complex open-pit mining facility in northern Canada in support of more accurate operational reporting of area-fugitive greenhouse gas emission fluxes from such facilities. WRF is studied in a series of sensitivity tests by varying topography, land use, and horizontal and vertical grid spacings to arrive at optimum configurations for reducing modeling biases in comparison with field meteorological observations. Overall, WRF shows a better performance when accounting for the mine topography and modified land use. As a result, the model biases reduce from 1.10 to 0.08 m s−1, from 1.04 to 0.50 m s−1, from 0.98 to 0.32 K, and from 45.7 to 17.3 W m−2, for near-surface wind speed, boundary layer wind speed, near-surface potential temperature, and turbulent sensible heat flux, respectively. Refining the model horizontal and vertical grid spacings results in bias reductions from 3.31 to 0.08 and from 0.80 to −0.11 m s−1 for near-surface and boundary layer wind speeds, respectively. The simulation results also agree with previous observations of meteorological effects on enclosed Earth depressions, characterized by formation of a cool pool of air, reduced wind speeds, and horizontal wind circulations at the bottom of the depression under thermally stable conditions. The results suggest that such configurations for WRF are necessary to arrive at more accurate meteorological predictions over complex open-pit mining terrains with similar features. [ABSTRACT FROM AUTHOR]

A new ensemble-based land surface data assimilation (DA) system is coupled with the atmospheric four-dimensional ensemble-variational data assimilation (4D-EnVar) system with the goal of improving the analyses within Environment and Climate Change Canada's Global Deterministic Prediction System. Since 2001, the sequential assimilation of surface variables is used to generate the initial conditions to launch the Global Environmental Multiscale (GEM) coupled forecast model. The work presented here is to replace the sequential DA with an independent surface DA system, the Canadian Land Data Assimilation System (CaLDAS) assimilating screen-level observations, and to compare assimilation experiments with CaLDAS run in uncoupled and weakly coupled modes. In the uncoupled mode, CaLDAS is used to initialize the forecast without interacting with the 4D-EnVar system. In the coupled mode, the analyses generated from CaLDAS and 4D-EnVar are used to initialize the forecast model. The analyses and forecasts from uncoupled and coupled runs are evaluated against surface and radiosonde observations over different subdomains to conclude the impact of coupling CaLDAS with 4D-EnVar. Results indicate a statistically significant reduction in bias and standard deviation at the surface for screen-level temperature and dewpoint temperature on the order of 0.1 K, and in the lower troposphere between 1000 and 500 hPa on the order of 0.1 dam for geopotential height and 0.1 K for air temperature and dewpoint depression in the coupled DA runs. The positive impact persists up to 5 days over some subdomains. It is concluded that the coupled DA approach generally performs better than the uncoupled version. [ABSTRACT FROM AUTHOR]

Slant satellite-viewing geometry is investigated for the simulation and assimilation of radiances in Environment Canada's weather forecast system. The standard approach is to extract from a short-term forecast (trial field) a one-dimensional vertical profile, located at the ground footprint of the observation location, to compute the model equivalent of the observation. Since in general, the lines of sight are not vertical, the observation operator can be improved by interpolating the trial field to the slant path. The interpolation relies on the estimation of horizontal gradients of atmospheric variables, extracted at each altitude layer, and representative of a neighborhood of approximately 100 km. Using these slant profiles of the trial fields yields significant improvements in the simulation and assimilation of radiance observations. The impact of this modification in assimilation and forecasts is evaluated over two periods of two months. Both short-term and long-term effects are found, obtained not only from the more accurate interpolation of the trial fields, but also from a revision of the observation error suggested by the lower observation minus trial field difference statistics. The error standard deviation of the geopotential height shows, in the upper troposphere and lower stratosphere, in 12–48-h forecasts, up to 10% reduction in the zonal value at higher latitudes, and 1%–2% reduction in the global value. Given these improved statistics, a revision of the radiance error statistics was explored. Applying these revised statistics suggests that the impact of a slant operator extends into longer range. [ABSTRACT FROM AUTHOR]