All of the operational scripts and programs, working either on the frontend or the backend machines were converted to ensure year 2000 compliance. A considerable effort has been devoted in 1999 towards that important issue.

The CMC ensemble prediction system for medium range was upgraded to 16 members in August 1999. The additional 8 members are produced using the unified model (GEM).

The mesoscale short range model (HIMAP) was upgraded from 15 km to 10 km resolution, along with improvements in its physics package.

The 30-day forecast production methodology was improved by the use of model’s climatology to get the temperature anomalies.

The blending method of the two numerical models producing the seasonal forecasts has been improved. This method uses a hybrid technique, based on twenty-six years of historical runs of seasonal forecasts for each model.

The following types of observations are presently used at the Centre. For these types, we use all observations that are available from the GTS, on the global scale. The numbers indicate typical amounts received during a 24-hour period :

*** A third of these are used for ice analyses.

b) Products

GRIB ECMF
GRIB KWBC
GRIB EGRR
FDCN KWBC
FDUS KWBC
U.S. Difax products
Significant weather forecasts
Winds/Temperature forecasts for various flight levels

Fully automated.

Various real-time quality control checks are performed for each observation received from the GTS. In particular :

• all reports are checked for gross errors;

• values for main items, such as height, pressure, temperature, dew-point and wind are checked to be inside physical and climatological limits;

• temperature profile check;

• hydrostatic check;

• horizontal check (spatial consistency with neighbours and first guess fields)

These checks are done at, or after, the decoding phase of the bulletins. Canadian observations are put on the GTS before such quality control is performed. However, Canadian observations are subject to quality control at the observing site, before transmission to the national centre.

The information generated by the quality control system inside the objective analysis is fed back into the observations database in order for non real-time monitoring and quality control activities to be performed. This monitoring is done on the global scale. Nationally, we also monitor the bursting altitude of upper-air soundings and results are distributed to data producers on a daily basis and monthly reports are distributed.

Each Canadian synoptic report (manned stations only) is also monitored in real time for completeness and timeliness. Requests to individual stations are made if certain criteria are met. Observing stations send corrections if time permits. These corrections are sent to the GTS for transmission. A monthly summary of errors is produced and distributed to data producers.

Monitoring the availability of observations on the global scale is an inherent portion of operations at the CMC. Information on the current content of the observational databases is available in real time, by observation types and by geographical areas. A chart showing the geographical distribution of observations, by types, used to initialize the numerical models is distributed to forecast centres across the country in real-time. A monthly report describing the availability of upper-air observations is produced and distributed to data producers.

The information on the availability and quality of observations available for use in the final global analyses is assembled each month into the "CMC Global Data Monitoring Report". The statistics presented in the reports are prepared in accordance to the WMO/CBS approved procedures. The reports are sent to the WMO Secretariat as well as other major GDPS centres.

In 1993, CMC was designated by CBS as the lead centre for the monitoring of the quality of land surface observations in WMO RA-IV (North and Central America). In 1994, the CMC began to fulfil its role and since then has regularly produced its 6-monthly reports entitled " Report on the Quality of Land Surface Observations in Region IV". Two such reports were distributed in 1999. Monitoring results are distributed directly to national focal points for most countries within RA-IV.

7.1 System run schedule

The following table summarizes the operational runs at CMC. The core of the operational runs executes in batch on the NEC SX-4/32. Most of the postprocessing jobs, including CMC products, execute on the front end computer (SGI Origin 2000).

Note: There are also runs (not described here) that perform surface objective analyses and update geophysical fields; these are runs G3, G4, G5, G6 and R6.

7.2.1   Data assimilation and objective analysis

7.2.3.1 Analysis

A series of classic analysis products are available in electronic or chart form ( i.e. surface analysis of snow and cover, sea surface temperature, surface MSLP and fronts, upper-air geopotential, winds and temperature at 1000, 850, 700, 500, 250 hPa, etc.).

7.2.3.2 Forecasts

A series of classic forecast products are available in electronic or chart form ( i.e. MSLP and 1000-500 hPa thickness, 500 hPa geopotential height and absolute vorticity, accumulated precipitation and vertical velocity, 700 hPa geopotential height and relative humidity). A wide range of bulletins containing spot forecasts are produced for many locations across the world. As well, other specialized products such as precipitation and probability of precipitation forecasts, temperature and temperature anomaly forecasts, etc., are produced.

The 16 member ensemble prediction system (EPS) runs once a day up to 10 days (Houtekamer et al., 1996; Lefaivre et al., 1997; Plante et al. 1999). Eight perturbed analyses are obtained by running independent assimilation cycles that use perturbed sets of observations and are driven by eight different versions of the spectral global model (SEF model T95, Ritchie, 1991). The number of perturbed analyses is doubled as follows: the mean of the analyses is subtracted to the operational analysis and a fraction of this difference is added to the original perturbed analyses. Every day, at 00 UTC, two separate models are used to produce the 10-day forecasts: the SEF model and the GEM model (resolution of 1.875°, Côté et al., 1998a and 1998b). Each model uses different versions of their physical parameterizations.

Ensemble outputs of the following products are available on the web (http://www.cmc.ec.gc.ca/cmc/htmls/forecasts.html, then select ensemble forecast system): spaghetti plots of the 500 hPa heights; composite MSLP highs and lows; accumulated quantity of precipitation; forecast charts of precipitation amount probability at various thresholds. Ensemble outputs are also used to feed the Perfect Prog statistical package to forecast daily maximum/minimum temperatures.

7.3.3.1 Analysis

A series of classic analysis products are available in electronic or chart form ( i.e. surface analysis of snow and cover, sea surface temperature, surface MSLP and fronts, upper-air geopotential, winds and temperature at 1000, 850, 700, 500, 250 hPa, etc.).

7.3.3.2 Forecasts

A wide variety of forecast products are available in electronic or chart form. These include the classic charts such as MSLP and 1000-500 hPa thickness, 500 hPa geopotential height and absolute vorticity, accumulated precipitation and vertical velocity, 700 hPa geopotential height and relative humidity. A myriad of special charts are produced in the context of the summer or winter severe weather (tropopause, stability indices, wind shear, helicity, wind chill, liquid water content, streamlines, low-level maximum wind, vertical motion, etc.) or in the specific support for aviation forecasting (icing, freezing level, height of cloud ceiling, momentum flux, turbulence, etc.). A wide range of bulletins containing spot forecasts are produced for many locations across North America.

A high resolution model is run once a day for 24 hours over two sub-areas of Canada: western Canada and upstream waters; Great Lakes and eastern Canada. This strategy was given the name of High Resolution Model Applications Project (HIMAP, Pellerin et al., 1998). The model used is the same unified model described in section 7.3.2, except for the following differences:

• the window of uniform grid has a resolution of .09 degrees (~10 km);

• the number of vertical levels has been increased to 35.

The model is started from the 6-hour forecast of the regional model following the 00 UTC run. Outputs of surface fields covering the uniform grid area are transmitted in GRIB formats to Canadian Regions. Series of coloured images (including animation) are also made available through the internal Web.

The CMC is able to provide in real-time air concentrations and surface deposition estimates of airborne pollutants. These fields are obtained from a 3-D long range atmospheric transport/dispersion/deposition model, named the "Canadian Emergency Response Model" or "CANERM". The main applications for this model have been for estimating concentrations of radionuclides and volcanic ash. Based on this operational capability, the CMC has been designed by the WMO as a Regional Specialised Meteorological Centre (RSMC) with specialization in Atmospheric Transport Modelling Products for Environmental Emergency Response. In addition, CMC has been designed by the ICAO as a Volcanic Ash Advisory Centre (VAAC).

Fields of wind, moisture, temperature and geopotential heights must be provided to CANERM. These are obtained either from the Global or the Regional forecast and objective analysis systems. Please refer to the above section 7.2 for more information on these NWP products.

Latitude, longitude and time of the release are necessary input parameters for CANERM. Estimates of intensity and duration of the release are also required. In the absence of actual source data, the standard default values adopted at the WMO's First International Workshop on Users' Requirements for the Provision of Atmospheric Transport Model Products for Environmental Emergency Response (September 1993) would be used. These are :

• uniform vertical distribution up to 500 m above the ground;

• uniform emission rate during the first 6 hours;

• total pollutant release of 1 arbitrary unit;

• type of radionuclide is Caesium 137.

CANERM was developed by Janusz Pudykiewicz of Environment Canada and is described in Pudykiewicz, 1989. The horizontal and vertical advection in the model are performed using the semi-Lagrangian algorithm of Ritchie, 1987. Diffusion is modelled according to K-theory. The diffusivities are constant in the free atmosphere but have a vertical profile in the boundary layer which is dependent on the state of the surface layer; the vertical diffusivity within the surface layer is approximated using the relations provided by the analytical theory of the surface layer. CANERM simulates wet and dry scavenging, wet and dry deposition and radioactive decay for selected tracers. Wet scavenging is modelled by a simplified statistical parameterization based on the relative humidity. The source term is modelled by a narrow gaussian distribution to simulate both the release and subgrid scale mixing.

CANERM can be executed in forecast mode up to day 10, using the operational Global forecast model, and up to 2 days using the operational Regional forecast model. CANERM can also be executed in hindcast mode using Global or Regional objective analyses. Presently, three horizontal resolutions are available : a resolution of 150 km on a quasi-hemispheric domain, a movable continental domain with a resolution of 50 km and a mesoscale domain with a resolution of 25 km. CANERM can be executed in both the Northern and the Southern Hemispheres.

(atmospheric transport/dispersion/deposition) products

Upon request from the appropriate WMO Member Countries Delegated Authorities, the CMC will provide the following standard set of basic products :

• three dimensional trajectories starting at 500, 1500 and 3000 m above the ground, with particle locations indicated at synoptic hours;

• time integrated pollutant concentration within the 500 m layer above the ground, in Unit-s/m³, for each of the three time periods. The duration of the first time period is between 12 and 24 hours starting at release time. For a release before 12 UTC, it ends at 00 UTC; for a release after 12 UTC, it ends at 12 UTC the next day. The second time period is the 24 hours following the first time period. The third time period is the 24 hours following the second time period.

• total deposition (wet and dry) in units/m² from the release time to the end of the third time period.

The standard set of products was agreed upon at the First International Workshop on Users' Requirements for the Provision of Atmospheric Transport Model Products for Environmental Emergency Response. The CMC can also provide charts of air concentration estimates for the surface, 850, 700, 500, 300 and 250 hPa levels as well as total surface deposition estimates, at 3 or 6-hour intervals, if required. All the products can be transmitted by facsimile, in real time, during environmental emergencies. In addition, CMC is designed by the ICAO as a Volcanic Ash Advisory Centre (VAAC).

The Canadian Global model is used to prepare ozone and UV Index forecast at the 18-hour projection time based on 00 UTC data and at the 30-hour projection time based on 12 UTC data (Burrows et al., 1994). A Perfect Prog statistical method is used for forecasting total ozone, which is then supplemented with an error-feedback procedure. UV Index is calculated from the corrected ozone forecast. Charts of the total ozone forecast and of the UV Index forecast are prepared and transmitted to the Regional Offices. Bulletins giving the forecast UV Index at an ensemble of stations across Canada are also generated.

Sea-state forecasts of 48 hours over the Eastern Pacific and Western Atlantic are generated twice a day (00 UTC and 12 UTC) by the WAM (WAve Modeling) model. The Pacific version of the wave model uses the surface level winds from the global model while the Atlantic version uses the regional model wind outputs. Various parameters are plotted on the wave forecast chart (wave height, swell period, swell height and direction, etc.).

7.5    Extended range forecasts (10-30 days)

Fifteen-day temperature anomaly forecasts (Verret et al. 1998) are generated once a week using a perfect prog approach from the medium-range model described at section 7.2.2.

Monthly temperature forecasts based on numerical weather prediction techniques, are issued at the beginning and mid-month of every month. An ensemble of 5 runs, obtained from 24-hour time lag, is produced. The model used is very similar to the former operational spectral global model (Ritchie, 1991), except it has lower horizontal resolution (T63 L23) and has evolving geophysical forcing: the anomalies (analysis-climatology) of sea surface temperature (SST) and snow, observed during the previous 30 days, are added to the daily climatology during the integration. Direct model surface temperature outputs ensemble means are averaged over the 30-day period and subtracted from model climatology obtained from a 26-year hindcast period (see section 7.6). These temperature anomalies are then normalised by the model standard deviation multiplied by .43 (to get equiprobable classes) and categorised in above, below and normal classes. Charts are produced, showing above normal, and below normal temperatures. Monthly forecast products are on the Web (address http://www.cmc.ec.gc.ca/climate/htmletc/m_fcst-e.html).

Seasonal forecasts are issued 4 times a year (at the beginning of March, June, September and December). Seasonal products are distributed internationally and nationally through Internet (address http://www.cmc.ec.gc.ca/~cmcdev/saisons/seasons.html on the Web). They are also distributed nationally on the National Telecommunications System and to selected users by facsimile and made available on electronic bulletin boards. The charts are accompanied by a verification chart giving the performance of the forecast over the hindcast period. Also, verification charts, showing the previous season's prediction and a preliminary analysis of the observed anomaly, are provided.

Season 1 forecasts are produced using a numerical approach (Derome et al., 2000). Two ensembles of 6 runs, obtained from 24-hour time lag, are produced: 6 from the T63 L23 model described in section 7.5, 6 from a general circulation model (GCM) (McFarlane et al., 1992) (T32 L10). Both models use the same initial operational analyses. SST anomalies, that have been observed over the previous 30 days, are added to climatological values over the period; snow is relaxed towards climatology at the end of the first month, except for the GCM, where it is a prognostic variable. A simple statistical linear regression equations relates the 1000-500 hPa thickness anomalies (forecast minus model climatology) to surface temperature anomalies, using regression coefficients for 90-day forecasts. Maps are similar to monthly ones: 3 classes, separated using the .43 standard deviation of observed climatology.

The precipitation forecast is produced using a more direct approach: the two ensemble means of forecast precipitation are subtracted from their respective models’ climatologies, and normalised by models’ standard deviations. These normalised forecasts are then added, divided by two and used to produce a map, categorised in 3 classes, using the .43 value for separation.

Skill maps of temperature and precipitation, as obtained over the 26 years of historical runs, are shown for each of the 4 seasonal forecasts periods.

Seasonal forecasts at lead time of 3, 6 and 9 months are produced, using a Canonical Correlation Analysis technique (Shabbar and Barnston, 1996). The technique uses the SST anomalies observed over the last year to predict temperature and precipitation anomalies at Canadian stations (51 for temperatures; 69 for precipitation) for the following 3 seasons. Maps of above, normal and below temperature and precipitation are produced. These are accompanied by skill maps, as obtained from cross-validation over a 40-year period.

The objective verification of the operational numerical models is done on a continuing basis. S1 skill scores, bias and root mean square error are produced for the Canadian verification area. A monthly verification summary is produced and distributed to our clients.

A verification system following the WMO/CBS recommendations has been implemented in 1987. Results are regularly exchanged with the other participating centres. The table on the following page is a summary of the verification scores for 1998 according to the recommended format.
Since 1994, CMC began an electronic exchange of these verification scores with other NWP centres. This electronic exchange proved to be very useful and results from various operational NWP models can now be compared much more easily and thoroughly.

Improvements will be made to the 3D-Variational system: revision of error covariance statistics, analysis on model’s eta levels, assimilation of additional type of observations (ACARS/AMDARS, TOVS/ATOVS, high density SATOBS).

The regional model resolution will be increased to about 16 km, with improved surface and condensation schemes. The model code will be converted to MPI and will then run on multiple nodes on the SX-4 machine.

Ensemble prediction system will be used to issue probabilistic products.

Air quality forecasts will become fully operational with the implementation of an updated version of the statistical system and the implementation of a dynamically based ozone forecast model.

The UMOS system will become fully operational.

Instead of issuing 4 «seasonal» forecasts per year, 90 day forecasts are to be issued every month, blending dynamical and empirical techniques.

10. References

Brasnett, B. 1999: A global analysis of Snow Depth for Numerical Weather Prediction. J. Appl. Meteor., 38, 726-740.

Côté, J., S. Gravel, A. Méthot, A. Patoine, M. Roch and A. Staniforth, 1998a: The Operational CMC-MRB Global Environmental Multiscale (GEM) Model: Part I - Design Considerations and Formulation, Mon. Wea. Rev. 126, 1373-1395.

Côté, J., J.-G. Desmarais, S. Gravel, A. Méthot, A. Patoine, M. Roch and A. Staniforth, 1998b: The Operational CMC-MRB Global Environmental Multiscale (GEM) Model: Part II - Results, Mon. Wea. Rev. 126, 1397-1418.

Houtekamer, P. L., L. Lefaivre, J. Derome, H. Ritchie and H. L. Mitchell, 1996: A system simulation approach to ensemble prediction. Mon. Wea. Rev. 124, 1225-1242.