CERFACS CERFACS Scientific Activity Report Jan. 2010 â Dec. 2011
CERFACS CERFACS Scientific Activity Report Jan. 2010 â Dec. 2011
CERFACS CERFACS Scientific Activity Report Jan. 2010 â Dec. 2011
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
DATA ASSIMILATION FOR ATMOSPHERIC CHEMISTRY<br />
address the monitoring of air quality with MOCAGE-Valentina. The MOCAGE model has the advantage of<br />
offering the possibility to work on four nested domains, from global to local scale. Working with the regional<br />
domain is essential for the modelling of air quality in order to have a sufficient horizontal resolution. For<br />
computational reasons it is not possible to increase the resolution of the global domain up to scales useful<br />
for air quality studies. The Valentina system has thus been further developed to produce chemical analyses<br />
on the limited area covered by the regional version of MOCAGE.<br />
The first requirement of the MACC project was to produce analyses of the concentration at the surface<br />
of several key atmospheric components. The analyses had to use the surface measurements from groundbased<br />
stations available each day for the day before. As Valentina is able to assimilate only one chemical<br />
component, we selected the ozone molecule. This choice was driven by the fact that ozone is one of the most<br />
observed species and because ozone is actively involved in air pollution. We decided to use the full chemical<br />
scheme of MOCAGE to obtain values for all the key atmospheric components of MACC. Because even if<br />
only ozone is assimilated, modifying the ozone concentration produces modifications in the concentration<br />
of other species due to the chemical links between them.<br />
The version of MOCAGE-Valentina that assimilates the surface ozone measurements was delivered to<br />
Météo-France/CNRM in mid-<strong>2010</strong> for a pre-operational daily production of the ozone surface analysis<br />
of the day before. We intensively used this pre-operational production to evaluate our analyses. In particular<br />
we have developed a monitoring suite that produces each day and each week some plots that are displayed<br />
on a web site to easily follow the evolution of the products. This monitoring allowed us to detect that we<br />
had to periodically update the error statistics both for the forecast and for the observation. We have then<br />
developed tools based on a posteriori diagnostics to automatically update the forecast and the observation<br />
error covariances matrices. The monitoring was also useful to point out the stations where the measurements<br />
were repeatedly not coherent with the model. This has lead to the creation of a blacklist for such stations,<br />
the blacklist being also automatically updated.<br />
(a) West-East length-scales<br />
(b) North-South length-scales<br />
FIG. 3.1: Diagnosed length-scales (in km) of the horizontal correlation of the forecast error, averaged on<br />
the period between 1 st and 10 th July <strong>2010</strong> : West-East (left) and North-South (right) length-scales.<br />
Among the work to ensure a good quality of the daily analysis of the surface ozone, we had begun a thorough<br />
investigation on the background error covariance matrix (BECM). In the operation suite, the background<br />
error variances comes from a posteriori diagnostics. The background error correlations are fixed in time and<br />
had never changed. We first diagnosed the BECM using an ensemble method with a methodology similar<br />
to the one we used for the global scale (see section 3). We have found that there are strong space and time<br />
variations in the background error covariances (Fig. 3.1), the use of constant values for the correlations as<br />
in the pre-operational suite is therefore not a good option [4]. We then derived several configurations of the<br />
78 <strong>Jan</strong>. <strong>2010</strong> – <strong>Dec</strong>. <strong>2011</strong>