|
Progress in research for ALADIN along the first 10 months of 2002
(contributions from thematic coordinators)
(postscript version)
Menu:
Coupling
· Spectral coupling :
Coding is achieved and basic tests with spectral coupling
have been done. Technically it works, but meaningful case studies are still
missing for tuning an optimal combination of Davies and spectral coupling. Now
experiments with a 1d shallow-water model are being performed.
· Surface pressure tendency coupling :
Situation did not change much. The code has been ported to
the last library, AL25T2, but further experimentation is required.
· Time interpolation :
Intensive work has been done on this subject. No satisfactory
solution has been found that would be competetive with increasing time
frequency of coupling files. Spectral analysis of the time evolution at the
Christmas' storm case indicated a good potential in finding diagnostics that
can be computed in the coupling model to determine the time-dependent coupling
frequency required by the situation.
· The work on
two-way-nesting was delayed. Instead an
investigation of the problem of coupling in ALADIN following a similar
approach as Aidan Mc Donald used for HIRLAM was launched.
Dynamics
· Stabilization of NH dynamics
A new prognostic variable for pseudo-divergence ( d4
) has been implemented in the main code of ALADIN. Stability analyses and
2d experiments showed it could improve stability in the presence of orography.
Besides a new time-discretisation called "predictor/corrector" was
also coded, to allow the use of a two-time-level semi-Lagrangian scheme, but
only for the old prognostic variables as a first step. The merging of the two
ingredients will be completed in the next weeks. These two ingredients are
expected to offer a robust enough scheme to be used with 2 time-levels
discretisations. Preliminary 2D tests are very encouraging.
Besides an analysis of the anelastic approach in ALADIN was
done. It was found that the anelastic approximation is not practically
applicable using the Laprise coordinate in a spectral semi-implicit model.
· Evaluation of dynamics at high resolution
The properties of ALADIN dynamics (accuracy, stability,
efficiency, consistency) were evaluated for increasing horizontal and vertical
resolutions, using the "quasi-academic" ALPIA framework. It consists
in 4 embedded models covering the French Alps, with mesh-sizes decreasing from
10 km to 1.25 km and the number of vertical levels increasing simultaneously
from 30 to 85.
The first part, dedicated to the comparison between Eulerian
and semi-Lagrangian stability criteria, was completed. The experiments were
done with three-time-level Eulerian and semi-Lagrangian non-hydrostatic
schemes. As outcome, there is now a robust result that the three-time-level
semi-Lagrangian scheme can be used still with at least double of the Eulerian
time-step without meeting Lipschitz criterion.
Then the experiments were resumed to investigate in details
the properties of ALADIN dynamics at very high resolution : hydrostatic versus
non-hydrostatic, explicit versus semi-implicit, Eulerian versus
semi-Lagrangian, ...
· Lower boundary condition for NH
An alternative, more consistent, implementation of the lower
boundary condition has been introduced in the main code of ALADIN. The
principle is to use vertical velocity (w) as a gridpoint prognostic
variable instead of the vertical divergence (e.g. d4), in the advection
process. This ingredient improves the response over mountains for the
semi-Lagrangian scheme in non-hydrostatic dynamics. Merging with the two other
main modifications is in progress.
· Radiative upper boundary condition
The work is slowly restarting, following a new approach based on recursive filtering.
· Gridpoint semi-Lagrangian diffusion
The work was resumed at the beginning of August, by a new
student, and has to start with an update of the corresponding code. No
significant progress is expected within the end of the year.
· Improved filtering of orography
A new method was proposed to smooth orography at high
resolution, which was proved necessary (to avoid spurious small-scale noise,
especially when using "linear" spectral truncations). An additional
spectral term is added to the cost function used to compute and optimize the
spectral orography, which acts at damping the smallest scales. This should
provide smoother spectra than the present method (importing the
"quadratic" spectral orography when using a "linear"
truncation), which abruptly sets to zero the contribution from the smallest
scales (the last third of the spectrum), and will be far more flexible.
However the first tests interfered with problems in the minimization library,
and the present formulation is to be reconsidered.
· Relaxation of "thin
layer" hypothesis
The corresponding modifications are now implemented both in
ARPEGE and ALADIN. Experimentation is in progress, to evaluate the practical
impact in terms of scores or forecast accuracy, mainly in tropical regions. An
extension to non-hydrostatic dynamics with the Laprise coordinate is under study.
Physics : precipitations (convection, microphysics, ...)
A development work plan on moist processes was adopted at the
meeting in Medulin.
Two main approaches of the cloud-precipitation scheme are
investigated, namely the Functional Boxes approach focusing on
"resolved" cloud and precipitation (and using the operational
diagnostic convective scheme as another box), and a more integrated approach
focusing on convection and using a microphysical scheme derived from the one
developed by Philippe Lopez.
· Integrated approach
(L. Gerard)
A version of the physical package under AL15 with a
microphysical scheme derived from Ph. Lopez's, but separating ice and liquid
cloud condensate and considering a pseudo-historical precipitation content, is
now available.
The prognostic convection scheme developed last year by L.
Gerard is now running under AL15. At this occasion, we observed an undesirable
stiffness of the tunings with the re-introduction of the modulation of
moisture convergence by the mesh size (GCOMOD=1). \\ Validation tests of the
scheme are going on in Vienna( S. Greilberger).
A version of the code with both modifications mentioned above
is also ready.
Work is going on on the new integrated scheme, implying a
separation of the downdraught and the updraught; the updraught routine under
test implies an explicit treatment of the suspended cloud condensate, full 3D
prognostic variables for the updraught velocity \SE{and} the updraught mesh
fraction, the production of detrained condensate instead of precipitation.
· Functional boxes approach
A communication of E. Bazile in HIRLAM Newsletter 41 (2002),
pp79-83, synthesizes the current status.
The problems of the ice phase seem now solved -- except maybe
a small difficulty around the triple point.
The functional boxes seem very sensitive to time and vertical
instabilities of the parameterisations, and can often reveal hidden weaknesses
in them.
·
Entrainment in clouds
The sensitivity of the deep convection scheme to the humidity
of the middle atmosphere has been tested within the EUROCS project frame. The
present entrainment formulations confer a pretty good sensitivity to the
scheme as compared to Cloud Resolving Models.
· Stratocumulus clouds
The underestimation of stratocumulus clouds has been
investigated on the FIRE I case. A diagnostic approach using a new cloud
contents and cloudiness has been proposed. A prognostic approach, base on the
same "Philippe Lopez" as mentioned above, is also under tests and progress.
· Other aspects
It was decided in June to start the development of a
shallow-convection scheme, based on a mass flux approach (Th. Haiden, J.M. Piriou).
Extended collaborations participations were planned in
Medulin, but the communication about this stays weak.
Physics : turbulent processes, boundary layer
· Andre Simon
: Study of the relationship between turbulent fluxes
in stable PBL and cyclogenesis (ALATNET PhD).
- Main goals are :
- explain the sensitivity of the forecast of rapid cyclogenesis on
parameterization of turbulent heat exchange (tuned above all by parameter USURID).
- find new methods to adjust the current parameterization of critical
Richardson number while avoiding the erosion of inversion
Results :
- The heat flux parameterization used in cycle 24 was tested on case study
of explosive cyclogenesis on 20 December 1998. Main sensitivity was found in
area eastward from Newfoundland. The role of increased heat exchange seems to
be in a decrease of static stability in the mentioned area, that allows bigger
vertical velocities and slower drop of pressure during the decaying period of
the cyclone. Thus the influence is indirect, because it acts before the rapid
reinforcement of the storm (42 hours later).
- The used value of USURID parameter (0.042 - operational in the beginning
of 2002) gives unrealistic forecasts of static stability in the PBL in areas
crucial for development of the cyclone. Even very high values of USURID (0.14)
give less static stability in comparison to model analysis.
- The Richardson number limitation is more important by values between 1
and 100 as by very high values (relative differences in K-profiles for current
parameterization reach their maximum by Ri=4.5 for various USURID and by
constant wind shear).
Suggestions :
- make the parameterization of heat exchange in the PBL dependent on
surface stability. This is not a final solution for the cyclogenesis problem
but keep more stability over cold surfaces (land) and allow instability over
warm lowest layers (sea). Similarly, differentiation should be made for
momentum exchange as well.
- make the height of K- decrease situation dependent (decrease or increase
the parameter UHDIFV for instable/stable PBL layers)
- vertical discretization of Rid gives very small differences, while
tuning it via parameter USURIDE. Thus new kind of vertical discretization of
ZUSURID should be used to allow "more freedom" for tuning.
· Janko Merse
is finishing an article for Meteorology and Atmospheric Physics
(79, 195-213, 2002) with the title: "Turbulent dissipation of the
cold-air pool in a basin: comparison of observed and simulated
development", where the MM5 model has been used to simulate the process
and it's been shown that MY 2.5 parameterization of the PBL (taken from the
ETA model) seems to describe the process quite well.
· Martin Gera
: Improved representation of boundary layer, main
goals for year 2002 :
- Analysis of anisotropy and inhomogeneity of turbulent processes.
Influence of stratification and surface to the behaviour of turbulence.
- Creation of new integral stratification dependence functions for
friction velocity and friction potential temperature, which express their
variability nearest the surface. The bases for the derivation are Businger`s
empirical functions and the static equation for TKE.
- Anomalous dissipation of energy close to the surface is expressed by
blocked kinetic energy. Some comparisons are done with classical Prandtl
mixing length expression near the surface. This approach modifies current
ALADIN mixing length expression and looks better when using the prognostic TKE
equation (it's in progress now).
- Implementation of TKE computation for stationary regime, it was the
first step for solving prognostic equation. We have concentrated on the
behaviour near the surface.
- Implementation of prognostic TKE equation and consecutive computation of
exchange coefficients. They have a direct dependency on TKE.
- The scheme is still unstable, we are assessing methods for solving the
nonlinear TKE equation.
· Alexander Kann and
Harald Seidl :
Strange oscillations of humidity and consequently cloudiness
at certain levels were encountered within the 1d version of ALADIN (cy22t1op8)
which did not appear in the 3D runs (al12 +cycora-bis). Concerning our stratus
diagnosis it seems that we have to lower the threshold values within the code
in order to gain more realistic Lifted Fog Forecasts. This seems to be
necessary as the initial state of vertical profiles already deviates from
observed soundings.
· Martina Tudor
performed a few case studies connected with Bura wind on the Adriatic for
the Croatian ministry for traffic (because the road is being built there). She
has used the dynamical adaptation of wind field to 2km resolution for that.
Physics : radiation
· Introduction of ozone
climatologic profiles varying with the month and the location, instead of a
unique profile. The impact has not yet been tested.
· Comparison of the
physics from the Climate model to the one from the NWP model. This test lead
to a closer assessment of the operational radiation scheme, and to testing
more other schemes in single column model (Morcrette, Graals and modified
versions of the current operational scheme).
· It appears that the
radiation has a strong impact on dynamics much sooner than was thought up to now.
Physics : surface
· The use of new
databases is presently in standby. Some tests of O. Latinne showed a slight
deterioration of the scores during the summer, and also dubious values of the
Leaf Area Index over boreal forests. This lead to light modifications of the
databases, and we are awaiting some feedback from the Climate Group which had
to test it. The main problem is the want for a person who would redo a whole
set of tests over one summer in assimilation mode.
· Phasing is on the way
for the improved snow parameterisation, on cycle CY25 which should be running
in double chain before the winter. This scheme and its tunings are described
in the proceedings of SRNWP/HIRLAM meeting of 22-22 October 2001.
Physics : orography
· From the side of
orographic forcing, the main event was the specialised mini workshop in
Toulouse that identified and prioritised three needs:
- the one of finding a new way to use data
bases, in particular regarding the current "quadratic cut-off" for
linear grid applications that should be replaced by a smoother transition,
this action being possibly twined with a revival of the concept of
"envelope at constant volume";
- the one to have again available the kind
of momentum budget diagnostic tools once developed by G. Gregorcic, in order
to better understand how the several components of our orographic forcing
(resolved mean, resolved envelope, gravity wave drag, form drag, lift -if
activated-) do interact at fine scale to deliver an as consistent as possible
large scale balance of angular momentum;
- that of a rather independent search for
the causes of the vexing problem of ill-placed (slopes rather than summit)
orographically induced precipitation.
· In addition, discussions taking place
at the recent informal RC-LACE training course, raised the following
hypothesis: are not we unable to modify the current unsatisfactory situation
where we need envelope despite form-drag, where lift cannot be incorporated
without deterioration of the results and where our form-drag-coefficient is
much too large, simply because we try to evolve incrementally (while a
simultaneous jump to more logical choices might be successful at once) ?
· Concerning the
optimization of the orography for operational models, nothing was done.
Physics : simplified physics
· Cornel Soci proposed
and tested successfully modifications of the parametrization of large-scale
precipitations, for use at a resolution of 10 km.His coming stay (early 2003)
in Toulouse will be the occasion to transfer his knowledge about specific LAM
constraints to C. Loo.
· In ARPEGE, some
modifications and tunings in vertical diffusion have been done in order that
the simplified version would be in better agreement with the complete one and
in order to avoid some instabilities appeared in linear models from increase
of resolution. A validation of the simplified direct physical package was also done.
Physics : validation
· Methods
Debugging of the Horizontal Domains Diagnostic package in the LAM is on the way.
Discussion in June in Medulin also planned :
- to set up a ``Gaussian'' control of the tendencies,
- to enter the code for extracting profiles from the 3D model for the
Single Column Model, into the main cycle,
but this is not yet realized.
· Experimentation
One is referred to J.F. Geleyn's document of "25 items".
It must be recognised that ALADIN physics is just emerging
from an incredible "black series" (not completely finished though)
as if the past five years of extreme stability and of concerted evolution via
big CYCORA-type packages had to be "paid" (probably because of
insufficient consolidation work at the time). The itemisation of the causes
(repeated below) is also very revealing of the diversity of sources for the
different problems that makes any strategical effort to prevent a future
repeat difficult - but for a strong increase in post-coding non-numerical
validation (something everyone hates to do!).
Statistically the 25 items are distributed as follows :
- Severe bugs : 2
- Small bugs : 2
- Design errors prior to or at the CYCORA-bis level : 4
- Design errors following CYCORA-bis : 2
- Forced scientific changes : 6
- Preventive actions : 1
- Scientific innovations : 5
- Tuning : 3
Last, a "nutshell" description of the physics operational in
ALADIN-France was produced.
· The planned "validation
chart" is not ready.
Physics : interface to dynamics
· Study of the interactions between
non-hydrostatic features and physical parameterisations
The problem of the introduction of diabatic forcing within
the predictor-corrector approach was addressed during the dedicated
mini-workshop. Safety locks to preserve the presently used alternatives will
be introduced. A more thorough investigation is delayed, waiting for a stable
dynamical core and more impact studies.
Besides a number of empirically derived approximations to the
treatment of the diabatic heating are being tested. There are labelled
quasi-anelastic and quasi-hydrostatic approximations to the diabatic heating
because they produce results analogous to those obtained with the
correspondingly approximated dynamical models, when using short time-steps.
· Use of new prognostic variables for physics
The interface for introducing new prognostic variables in
convection has been developed and should enter the next library (end 2002),
together with the parallel changes for using a TKE-based turbulence parameterization.
· Sensitivity of physics to the
length of time-step and the vertical discretization, fibrillations
The large vertical oscillations exhibited by Martina Tudor
last autumn while trying to diagnose the PBL height, and problems noticed when
increasing the operational vertical resolution, led to reconsider the
anti-fibrillation scheme. First, it is now obvious that a retuning is required
for each change in resolution. Second, the vertical variations of the
coefficient controlling the anti-fibrillation scheme are limited, with a
maximum equal to the value at the lowest level, in order to avoid the vertical
slicing. This prevented problems in operations, waiting for a more thorough
investigation of the behaviour of physical parameterizations as longer
time-steps are used.
Eric Bazile started, using the single-column model and
addressing the parameterizations of vertical diffusion and convection. The
response of the vertical diffusion scheme iss very regular, but for convection
very steep jumps are observed. This is due to the fact that the triggering of
convection depends on the tendencies issued from vertical diffusion at two
adjacent vertical levels, and the scheme is very noisy on the vertical. A
correction was proposed and the first 1d tests performed.
Martina Tudor afterwards undertook a thorough study of the
ARPEGE/ALADIN-physics' susceptibility to stiffness and nonlinear instability.
This was done by halving the time-step used to compute fluxes (and not the one
to apply their effect) for detecting stiffness and by applying the mirror
trick of the previous one in space (thicker levels only at the time of the
flux computation) for tracking nonlinear instability.
In the first case, the vertical diffusion problem was
identified again and the new anti-fibrillation cure positively reassessed.
More surprisingly a potential source of stiffness was detected in the
condensation/rain-reevaporation cycle but the maintaining mechanism is
"dry" (either in vertical advection and/or in vertical diffusion)
which forbids any practical action in the current state of the physics.
In the second case the test allowed to find a small bug in
the protection against nonlinear instability for the pseudo-advection by the
convective mass-flux. This discovery indirectly proves the power of the method
since no other problem was spotted once the bug being corrected.
Interestingly, the remaining open question corresponds to one known problem in
the TL/AD regularised physics (together with the one of fibrillations). Hence
the whole picture in this area of research is now more consistent, pending
further study of the huge amount of accumulated results.
· Smoothing of the shallow-convection parameterization
The smoothing proposed last year by Martin Bellus on a
suggestion of Eric Bazile was cleanly introduced into the code and revalidated
in order to prepare a possible operational implementation. Its main aim is to
suppress a coupled time/space oscillation (the second aspect being mainly but
not only in the vertical) created by the on/off character of the
shallow-convection parameterisation, which cannot be handled by the
anti-fibrillation scheme.
Most problems are damped, but some oscillations still resist.
The ad-hoc character of the time-smoothing scheme underlines the need to code
a new scheme.
· Impact of time-stepping
Ilian Gospodinov started to analyse various techniques for
introducing physical tendencies within the semi-implicit approach, with
hydrostatic dynamics.
Data assimilation : methods
· Local pre-operational
3D-Var assimilation suites using raw ATOVS or aircraft data (including local
tunings of statistics, observation use, cycling and initialization) :
This work is still under progress, with significative
advances in every center involved :
- Budapest :
Several solutions for a local assimilation cycle based on
3D-Var are under evaluation, with emphasis on coupling issues (choice of
coupling data, simple versus double nesting) and injection of local raw
radiance data. In parallel, experience on the objective evaluation of
assimilation cycles is gained. A pre-operational, stabilized 3D-Var
assimilation cycle is planned for completion for ???
- Casablanca :
Preliminary assimilation suites, based on the principle of a
combined use of DFI-blending and 3D-Var analysis, have been installed in
research mode for the ALADIN/MOROCCO model. A new assimilation cycle, again
blending+ 3D-Var, is under construction for the ALADIN/NORTH-AFRICA domain.
This "blendvar" cycle should be run using ATOVS raw radiances for
pre-operational, then operational, use early 2003.
- Prague :
The evaluation of the "blendvar" type of
assimilation cycles has been continued, in research mode only, in order to
study both individual cases (where pathological behaviours such as when very
localized storms can appear) and perform objective verification (scores). The
future plans for operations will be finalized after migration to the new
computer in 2003, with probably tests of DFI-blending+SurfCANARI followed by
DFI-blending+ 3D-Var+SurfCANARI.
- Toulouse :
An important work has been devoted to the stabilization and
validation of the main algorithms under the new technical environment (ODB,
raw radiances, CY25). A renovated version of the assimilation scripts is in
progress, which should eventually lead to a stable, user-friendly tool for
"switch-on/switch-off" ALADIN assimilation cycles.
· Local implementation
and tuning of blending :
3 centers have implemented blending cycles :
- Casablanca :
in research mode for ALADIN/MORROCCO and under preparation
for ALADIN/NORTH-AFRICA; more evaluation to be done in late 2002 and 2003
- Prague :
maintenance of the LACE operational blending suite
- Toulouse :
A renovated blending suite is in preparation, for the
ALADIN/FRANCE linear grid (41 vertical levels). This suite will be evaluated
in the 4th quarter of 2002, for operational purpose.
· Prototype version of
3d-FGAT :
This work has started in September 2002, mainly as part of
Vincent Guidard's PhD developments. The technical work has progressed
(4D-screening is ok), and a first scientific evaluation will be performed
using the well-known MAP IOP14 case. This work will be continued by a more
in-depth re-evaluation of the problem of coupling between LAM and global data
assimilation suites : evaluation of LAM versus LBC error growth, evaluation of
different strategies for the injection of the global model solution into the
LAM (2003-2004).
· Investigation of the
problems related to double-nesting :
This work has started under the frame of Steluta Alexandru's
ALATNET PhD. Several strategies for the nesting have been defined, with
coupling data chosen either from ARPEGE or from ALADIN/LACE (simple versus
double nesting) and either space- or time-consistent lateral boundary data.
Comparisons with the operational ALADIN/HUNGARY dynamical adaptation are also
performed. So far, the results show a significative fit towards the
observations with 3D-Var, but 6 hour forecast scores only marginally benefit
from the analyses. The potential imbalances of the 3D-Var increments have also
been checked, and it has been shown that dynamical imbalances are not
generated significantly by the 3D-Var analyses.
· A posteriori
validation tools and tuning of Jb / Jo; completion of the PhD work of Wafaa
Sadiki :
This work has been re-started this summer in Toulouse,
keeping a quite simple formulation of the a-posteriori functions. For a
reminder, the main diagnostics are : the ratio Jmin/P, the estimated global
background error variance rescaling "Sb" (=Jbmin/Tr(KH)) and the
estimated global observation error variance rescaling "So" (=Jomin/Tr(Ip-HK)).
The main results so far:
- There is a clear need to compare the a-posteriori results with
independent diagnostics (such as explicit statistics on the innovation
vector), in order to have a more in-depth understanding of the possible
misfits between modelled and real error statistics. Indeed, it appears that
too many ingredients (possible misfits) enter these diagnostics so that extra
knowledge is needed for the understanding. For example, Wafaa is now comparing
her results with the output of a Lönnberg-Hollingsworth method for
auto-covariances, for comparison of both variances and correlation length-scales.
- There remains a fundamental difficulty in assessing the amount of
scale-selective error variances, since the innovation vector which is involved
in the a-posteriori validation contains all the scales. A scale-selective
approach could be a possible future direction of research, for a potential candidate.
- The diagnostics exhibit clearly case-dependent behaviours, which
indicates that the error statistics in the LAM are sensitive to weather
systems (eg. synoptic storm inside the ALADIN domain, anticyclonic state, ...)
- A posteriori validation apparently can be used for re-tuning of
statistical parameters, although the work with real data obliges one to depart
from the strict theoretical framework and imagine tractable, pragmatical
applications. For example, the perturbation method to compute the
"Trace" operators demands to aggregate daily datasets together, and
to verify afterwards that the (neglected) time-correlations do not spoil the
results too much. Also, one needs to assess the stationnarity of the
diagnostics over the period of aggregation, which in practice means that the
diagnostics are more like monthly climatology.
Wafaa Sadiki will write her PhD manuscript in 2002/03, for a
scheduled PhD defense in early summer 2003. The work on a posteriori
validation will be one major contribution to her PhD. A continuation of this
work is planned in 2003 in Budapest.
· Improved Jb (new
methods, geographical modulation, part of humidity, ...) :
The work on the geographical dependency of the background
error statistics has been continued. One-dimensional variations of the
standard-deviations and of the length-scales with respect to latitude were
calculated from the bi-Fourier coefficients, and they were successfully
validated by comparison with the results of Mohamed Raouindi based on
mono-Fourier coefficients. This diagnostic work will be continued by
investigating the two-dimensional variations.
A formalization of the different variants of the NMC method
was proposed. This gives some indications on the formal contributions to
forecast differences, which can come from differences in the initial
conditions, in the lateral boundary conditions, and in the models (ARPEGE
versus ALADIN). For instance the emphasis on meso-scale features that is
provided by the so-called "lagged NMC method" is consistent with the
use of initial differences related to ARPEGE and ALADIN model differences, and
with the absence of lateral boundary differences.
An innovative work on the use of wavelet functions for the B
matrix representation has also started. Standard-deviations and correlation
functions implied by a diagonal wavelet covariance matrix were examined in
particular. The results appeared to be intermediate between those implied by a
diagonal approach in spectral and gridpoint spaces respectively.
Background error standard deviations have also been evaluated
for the pre-processing of IASI sounding data (Malgorzata Szczech in Toulouse,
ALATNET PhD).
In the frame of ARPEGE data, the use of ensemble analyses for
the computation of the background errror statistics is under test (Margarida
Belo Pereira, in her ALATNET PhD). The production of such ensembles has been
performed, and the evaluation of the physical properties of the derived
statistics, with a comparison with the NMC-derived data, is under way.
· Sensitivity studies;
completion of the PhD work of Cornel Soci :
The ALATNET PhD work has been completed in the first quarter
of 2002. The overall conclusions concerning the evaluation of the adjoint
model solutions in ALADIN, and the sensitivity of precipitation forecasts to
changes in the initial conditions, are as followin g:
- The perturbation of the lateral boundary conditions by
projected elements from the adjoint model is probably leading to a dead-end
issue, in the frame of the classical Davies' relaxation.
- The use of the ARPEGE simplified physics exhibits the
development of sometimes very strong instabilities, but that are the
high-resolution equivalents of the already known flaws of the large-scale
linearized physics. At this stage, it remains an open question whether one
should devote a significative amount of work to the linear physics for high
resolution 4D-Var or rely on very simple formulations, as was done a few years
ago in the global models. Especially, a very selective application of the
linear physics, and possibly more in-depth reformulation of the simplified
parametrizations, cannot be excluded. It is believed that this issue will
re-appear in a few years, for high-resolution adjoint computations.
- Cornel's test cases have shown a wide range of behaviours,
with respect to the sensitivity to initial conditions, for cases of
misforecast precipitations : no impact because the adjoint solution moves
outside the limited area of the model, no impact because model shortcomings
(in the parametrizations) are at stake, existing impact either due to wrong
vertical wind profiles or wrong low-level conditions (controlling the onset of convection).
Cornel Soci is now writing his PhD manuscript for the
University of Bucarest, and plans his defense for early 2003.
Data assimilation : observations
·
ODB : training, design of new tools for ALADIN + local implementations
- training in Toulouse for ALADINers.
-
design of new tools for ALADIN : geographical selection of observations for
C+I zone.
-
ODB is running with cycle 24T1 in Budapest, Prague and Casablanca, in research mode.
- Morocco
has established a stable, local observational database and can now handle a
local, autonomous database collection (the next step will be to add local
radiance data in it).
·
Screening : problem of observations close to lateral boundaries, local tunings
No
significant development was performed for screening, except on the thinning of
AIREP data. It is planned to look in more details at the ATOVS raw radiance
data sampling.
·
Using aircraft data in local ALADIN 3D-Var suites
No
major progress.
· Using observations of Hu2m for
upperair assimilation
The purpose was to study the possible use of screen-level
humidity observations in 3D-Var analysis within ALADIN-LACE. The results with
lagged statistics compared to standard ones lead to decrease the values of the
increments for specific humidity and temperature both in vertical and
horizontal directions.
·
Impact of the drift of radiosondes
This
work has been accomplished in the frame of LACE Center activities. It was
based on the monitoring of the raw sounding data from 2 stations, with
explicit recovering of the (X,Y) positions from each individual measurement.
The study has revealed a weak positive impact.
·
Use of raw ATOVS data in ALADIN 3D-Var
Promissing
progresses are reported from Budapest but there are still problems in the
reference library (CY24T1) in Toulouse for assimilating the ATOVS raw
radiances in ALADIN 3D-Var. An
adaptation on ALADIN of bias correction procedure, originally designed for
ARPEGE, is under development, but already well advanced in Budapest and Toulouse.
Sensitivity
studies on sampling (when extracting the data from the database) and thinning
(when assimilating the data) with respect to Jb structure functions (standard
Jb / lagged Jb) need to be performed (in 2003).
This
topic will receive again much attention in 2003, as it is interesting several
Centers ready to participate actively (Budapest, Casablanca, Toulouse).
·
Use of pseudo-observations (from MAP) in 3D-Var
One
case study of using humidity pseudo-observations from the MAP POI-14 provided
by GMME (Météo-France) has shown slightly positive impact on
precipitation forecast in ALADIN-FRANCE.
More
case studies should be performed in 2003 with improved pseudo-profiles on the
North Africa domain within a coordination/collaboration with GMME
·
Use of AIRS/
IASI data
over land : continuation
Normal
progress in the framework of an Alatnet PhD (Malgorzata Szczech) :
-
Specification of the background error covariance matrix is done for surface
temperature and will be done for emissivity.
-
Homogeneous bands have been selected in the infrared domain, based on mean
transmittances on a profile database. For each of these 18 infrared spectral
bands, emissivity maps will be produced from a database used by the MODIS team
for 18 types of surfaces in the model. Validation will be performed with MODIS emissivities.
-
Channel selection and retrieval of surface parameters is planned for 2003.
Data assimilation : surface
·
New background error statistics for 2m analyses
The Lönnberg-Hollingsworth's method has been used to
compute the background error statistics for 2m temperature and relative
humidity. The obtained 2m background error correlations decrease more rapidly
with distance than the operational ones, in agreement with the statistics used
in CANARI-Diagpack for mesoscale screen-level analyses. A new function
representing the background error correlations and rescalings have been proposed.
·
Spatial smoothing of soil moisture
Because of the small scales described by the model
(precipitation, cloudiness, surface characteristics), the threshold conditions
used in the surface assimilation, and the long time-scale of the evolution of
the soil moisture, the soil water content in ARPEGE and ALADIN models present
some unrealistic high spatial heterogeneities. This feature is not controlled
by the surface analysis since its length-scale is much smaller than the
surface observation network density, and is responsible for spatially noisy 2m
forecasts in case on strong solar radiation.
A spatial smoothing has been developed for ARPEGE and ALADIN
within CANARI analysis, controlled by two tuning parameters : the spatial
characteristic length of the smoothing and the number of iterations of the
process. The optimal values have been studied for ARPEGE and ALADIN. The
impact of this smoothing has been evaluated on a 2 months summer period with a
small positive impact on 2m forecasts and without any long term drift.
·
Combination of upperair and surface assimilations
The potential combinations between O.I. surface analysis and
upperair assimilation were examined in two frameworks : in ARPEGE with 4D-Var
and in ALADIN-LACE with spectral dfi-blending.
The upperair fields analysed with 4D-Var are currently used
as input for CANARI surface analysis. This choice aimed at reducing the 2m
errors related to the lowest model errors, but this is not fully in agreement
with the optimal interpolation theory which considers two uncorrelated sources
of informations : the observations and a first guess. Moreover since the
replacement of 3D-Var by 4D-Var the upperair analysed fields are in
equilibrium with surface fields at initial time, which differ from the
background surface fields (via a 3h forecast with different atmospheric
fields). This inconsistency may have locally some negative impacts. There is
also one technical argument in favour of using upperair guess as input to
surface assimilation : it allows more flexibility on the operational suite
because the upperair and surface analyses could be run in parallel. For these
reasons the impact of using the upperair guess has been evaluated. The
differences between the two methods on 2m analysed fields are generally small.
It appears very difficult to prove the superiority of one choice compared to
the other.
The
objective was to replace the surface blending by the CANARI surface analysis
in ALADIN. The experiments have shown that the use of CANARI surface analysis
in the blending cycle improves the low level forecasts on temperature and
humidity. The experiment where the surface analysis is done first before the
upperair blending is the most balanced one.
· Variational surface analysis
A variational surface analysis is being developed in the
framework of the European ALATNET project. The validity of the linear
hypothesis for the evaluation of the observation operator and the decoupling
hypothesis between gridpoints have been verified. Different strategies for the
evaluation of the observation operator were studied. The configuration with
two chess-type perturbations was found to be the best compromise between an
accurate evaluation and a reasonable computational cost. No empirical
condition is necessary to control the soil moisture corrections, which is an
improvement compared to the current operational optimal interpolation surface
analysis. An article is submitted to QJRMS and G. Balsamo is writing his PhD
manuscript. This work will continue in 2003 in the framework of the ELDAS
(European Land Data Assimilation System) project.
· Snow cover analysis
No
progress.
· SST analysis
Some work was done in 2002 to improve the use of a land/sea
mask in the guess estimation at the observation point for the SST
measurements. We noticed during the first months of 2002 some spots in the SST
analysis, especially near the coasts, in the North Sea and in the Baltic Sea.
Investigations have shown that these spots were due to a wrong use of SST
measurements from oil platforms and coastal stations. So, some modifications
have been made in the code in order to consider platforms as SHIPs on one hand
and to produce two first-guess values for surface temperature on the other
hand : one valid for land measurements (T2m and Hu2m observations) and the
other one for sea measurements (SST), for the use of coastal observations. A
procedure to avoid using observations at isolated land or water points has
been designed too.
The first steps in the use of satellite data to better
describe sea-ice limits are ready (interface to observations).
· Diagpack
An interface to ODB and new tools were developed. The work on
the use of aircraft data is progressing, but no significant impact was noticed
despite of retunings.
The CANARI documentation has been updated.
Predictability
In the medium-term research plan it was the first time that
the predictability topic appeared. However for the time being effective work
was carried out only at Météo-France, where ensembles were
computed taking into account 10 predictions (and generating own
perturbations). Meanwhile the interest is growing at the Member States and on
top of that the first LAM-EPS workshop was organised in the SRNWP framework in
Madrid at the beginning of October. Based on the discussions and
recommendations of this workshop and internal ALADIN discussions a working
plan is planned to be compiled within the end of this year, where the main
objectives, steps, milestones will be identified, and the effective work can
start at the beginning of 2003. For the time being there is a clear interest
for the topic in Hungary and Morocco, however all the Partners will be asked
to identify their interest and their level of contribution to this very new
and important project of ALADIN.
|