Short time augmented extended Kalman filter for soil analysis: a feasibility study
|
|
- Joan Lambert
- 5 years ago
- Views:
Transcription
1 ATMOSPHERIC SCIENCE LETTERS Atmos. Sci. Let. 13: (212 Published online 13 July 212 in Wiley Online Library (wileyonlinelibrary.com DOI: 1.12/asl.394 Short time augmented extended Kalman filter for soil analysis: a feasibility study Alberto Carrassi, 1,2 * Rafiq Hamdi, 1 Piet Termonia 1,3 and Stéphane Vannitsem 1 1 Institut Royal Météorologique de Belgique, Bruxelles, Belgique 2 Institut Catalá deciéncies del Clima (IC3, Barcelona, Spain 3 Department of Physics and Astronomy, Ghent University, Ghent, Belgium *Correspondence to: A. Carrassi, Institut Catalá de Ciéncies del Clima, Carrer Doctor Trueta 23, 85 Barcelona, Spain. acarrassi@ic3.cat Received: 4 August 211 Revised: 17 May 212 Accepted: 24 May 212 Abstract This paper presents a soil analysis scheme based on an extended Kalman filter (EKF, the short time augmented extended Kalman filter (STAEKF, where the model parameters are estimated along with the system state. We use an off-line version of the interaction soilbiosphere-atmosphere model for the assimilation of screen-level temperature and relative humidity. The leaf area index, the albedo, and the minimum stomatal resistance are estimated with the STAEKF. Results show that the STAEKF gives encouraging results and indicate the benefit of the simultaneous estimation of the state and parameters in soil data assimilation. Copyright 212 Royal Meteorological Society Keywords: data assimilation; parameter estimation; soil analysis 1. Introduction Current soil moisture analysis schemes in most numerical weather prediction (NWP centers are based on analyzed or observed screen-level variables. An optimal interpolation technique (OI developed by Mahfouf (1991 on the basis of the short-range forecasts of 2 m temperature (T 2m and relative humidity (RH 2m is operational in a number of centers: Météo-France (Giard and Bazile, 2, Environment Canada (Bélair et al., 23, and the European Centre for Medium-Range Weather Forecasts (Douville et al., 2 among others. However, some fundamental limitations are inherent with the use of the OI (Balsamo et al., 24; Mahfouf et al., 29, among these the fact that it does not allow for the use of nonlinear observation operator (note that the physical link between near-surface observations and soil variables is nonlinear and that the OI coefficients are usually derived using very simple assumptions and/or model configurations. A particularly promising approach to address the above shortcomings is the use of an extended Kalman filter (EKF. The EKF and a simplified version (SEKF, in which a static forecast error covariance matrix is adopted, have been successfully applied to the assimilation of direct and/or indirect soil observations (Draper et al., 29; Drusch et al., 29; Mahfouf et al., 29; Albergel et al., 21; Mahfouf and Bliznak, 211. Recently, Carrassi and Vannitsem (211a introduced an alternative formulation of the EKF where the uncertain model parameters are estimated along with the system state variables. The algorithm, referred to as short time augmented extended Kalman filter (STAEKF, uses a deterministic formulation for the model error dynamics (Nicolis, 23; an approach already used for the treatment of the error arising from the unresolved scales (Carrassi and Vannitsem, 211b and in the context of variational assimilation (Carrassi and Vannitsem, 21. With the aim of contributing to the current discussion on soil data assimilation strategies, we undertake here a set of numerical twin experiments with observations designed to test the STAEKF for the assimilation of screen-level observations using an offline version of the land surface model, Interactions between surface, biosphere, and atmosphere (ISBA (Noilhan and Mahfouf, 1996 where errors in three land surface parameters, the leaf area index (LAI, the albedo, and the minimum stomatal resistance (RS min are introduced. 2. The land surface model The two-layer version of the land surface model ISBA describes the evolution of soil temperature and moisture contents based on the force-restore method. The model is available within a surface externalized platform (SLDAS; Mahfouf, 27. The equations can be formally written as a dynamical system, dx dt = g(x, λ. The state vector, x = (T s, T 2, w g, w 2, contains the surface and deep soil temperatures T s and T 2 and the corresponding water contents w g and w 2. The vector λ is taken to represent the set of Copyright 212 Royal Meteorological Society
2 STAEKF for Soil Analysis 269 model parameters. A detailed description of ISBA can be found in the study by Noilhan and Mahfouf ( Surface assimilation schemes 3.1. Extended Kalman filter (EKF Let us write the model equations as a mapping from time t k to t k+1, x f k+1 = M(xa k, where xf and x a are the forecast and analysis state respectively, of dimension I,andMis the nonlinear model forward operator. A set of M observations of the true system, stored as the components of an M -dimensional observation vector y o, is supposed available at the regularly spaced discrete times t k = t + kτ, k = 1, 2..., with τ being the assimilation interval, y o k = H(y k + ɛ k, with ɛ being the observation error, assumed to be Gaussian with known covariance matrix R and uncorrelated in time, and y is taken to represent the unknown true state; H is the (possibly nonlinear observation operator mapping model variables to observation space. The EKF analysis state update equation reads x a = [I KH]x f + Ky o. The analysis error covariance, P a, is updated through: P a = [I KH]P f.thei M gain matrix K is given by K = P f H T [HP f H T + R] 1 with P f being the I I forecast error covariance matrix and H the linearized observation operator. The matrix P f is obtained by linearizing the model around its trajectory between two successive analysis times. Model error is assumed to be a random uncorrelated noise whose effect is modeled by adding an error covariance matrix, P m, at the forecast step so that, P f = MP a M T + P m, with M being the tangent linear model. Note that while P f is propagated in time and is therefore flow dependent, P m is defined once for all and it is then kept constant Short time augmented extended Kalman filter (STAEKF The forecast model, is augmented with P model parameters: z f = [ x f λ f ] = F(z a = [ M(x a ] F λ (λ a, (1 z = (x, λ is the augmented state vector. The augmented dynamical system F includes the dynamical model for the system s state, M, and a dynamical model for the parameters F λ. In the absence of additional information, a persistence model for F λ is often assumed so that F λ = I and λ f t k+1 = λ a t k ; the same choice has been adopted here. The forecast/analysis error covariance matrix, P f,a z, for the augmented system reads P f,a z = ( P f,a x P f,a T xλ P f,a xλ P f,a λ (2 where the I I matrix P f,a x is the error covariance of the state estimate x f,a, P f,a λ is the P P parametric error covariance and P f,a xλ the I P error correlation matrix between the state vector, x, and the vector of parameters λ. These correlations are essential for the estimation of the parameters. In general, one does not have access to a direct measurement of the parameters, and information are only obtained through observations of the system s state. In the controlled experiments given in the next sections the knowledge of the true parameters is used to diagnose the algorithms performance. The forecast error propagation in the STAEKF is given by P f z = CPa z CT, with C being the STAEKF forward operator defined as ( g C = M λ λ aτ (3 I P where I P is the P P identity matrix. Equation (3 embeds the key feature of the STAEKF; the presence of the term g λ λaτ allows for accounting for the contribution of the parametric error to the forecast error as well as to the error correlation between model state and parameters. An augmented observation operator is introduced, H z = [H ] with H as for the standard EKF. Its linearization, H z is now a M (I + P matrix in which the last P columns contain zeros. The augmented state and covariance update complete the algorithm and are equivalent to those of the EKF except that they refer now to the augmented system, and the gain matrix has dimension (I + P M (see Carrassi and Vannitsem, 211a, for more details. 4. Experimental setup We follow here a standard procedure in data assimilation studies, known as observation system simulation experiments (OSSE, in which a trajectory, solution of the model equations, is taken to represent the true unknown evolution which is then sampled at discrete times to obtain simulated observations. The estimated model trajectory is given again by a model solution but using a different set of parameters and initial conditions. The forcing data are the same for the truth and the assimilation solutions. They consist of 1 hourly air temperature, specific humidity, atmospheric pressure, incoming global radiation, incoming long-wave radiation, precipitation rate, and wind speed relative to the ten summers in the decade extract from ECMWF Re-analysis ERA4 and then dynamically downscaled to 1 km horizontal resolution over Belgium (Hamdi et al., 211. The fields are then temporally interpolated to get data consistent with the time resolution of the integration scheme of ISBA (3 s. In this study, ISBA is run in one off-line single column Copyright 212 Royal Meteorological Society Atmos. Sci. Let. 13: (212
3 27 A. Carrassi et al. T s T 2.5 EKF STAEKF W g W [m 2 /m 2 ].3 Leaf Area Index.25.2 [m 2 /m 2 ] 2 1 P a LAI.5 Figure 1. RMS estimation error in the four state variables for the EKF (red and the STAEKF (blue. The RMS error in the estimate of LAI obtained with the STAEKF is shown in the bottom-left panel, while the corresponding estimated error variance, P a LAI, in the bottom-right panel. mode for a 9 period, and the forcing parameters are those relative to the grid point closest to Brussels. An one-point soil model has been also used by Orescanin et al. (29, for parameter estimation using an ensemble-based assimilation algorithm. The simulated observations are T 2m and RH 2m, interpolated between the forcing level ( 2 m and the surface with the Geleyn s interpolation scheme (1988, at, 6, 12, and 18 UTC. The assimilation interval is τ = 6 h, while the observational noise is drawn from a Gaussian, N(, R, with zero-mean and covariance given by the diagonal matrix R with elements: diag(r = (σ 2 T 2m,σ 2 w 2m = (1K 2, 1 2. As explained in Mahfouf et al. (29, the observation operator H, relating the state vector to the observation includes the model integration. The initial P a and P m required by the EKF are set as diagonal with elements diag(p a = (σ 2 T s,σ 2 T 2,σ 2 w g,σ 2 w 2 = (1 K 2, 1 K 2, 1 2, and diag(p m = (σ 2 T s, σ 2 T 2,σ 2 w g,σ 2 w 2 = ( K 2, K 2, 4 1 4, (Mahfouf, 27. Parametric errors are introduced by perturbing either alternatively or simultaneously the LAI, the albedo and the RS min. These parameters strongly influence the surface energy balance budget and partitioning, which in turn regulate the circulation patterns and modify the hydrological processes. According to the STAEKF setup, we compute (not shown the functional dependence of the model equations on the 1 The values of σ wg and σ w2 are expressed as soil wetness index SWI = (w w wilt /(w fc w wilt where w fc is the volumetric field capacity and w wilt is the wilting point. uncertain parameters, the term g λ λ, required in Equation (3, relative to the ISBA model for the LAI, the albedo, and the RS min. For each summer in the period , a reference trajectory is generated by integrating the model with LAI = 1 m 2 /m 2, albedo =.2 and RS min = 94 s/m. Around each of these trajectories, Gaussian samples of 1 initial conditions and uncertain parameters are used to initialize the assimilation cycles. The initial conditions are sampled from a distribution with standard deviation (σ Ts,σ T2,σ wg,σ w2 = (5 K, 5 K, 1, 1, whereas LAI, the albedo, and the RS min, are sampled with standard deviations, σ LAI =.5 m 2 /m 2, σ albedo =.5 and σ RSmin = 5 s/m respectively (Ghilain et al., 211. The initial P a λ in the STAEKF read P a LAI = 1 (m2 /m 2 2, P a albedo = 1 4 and P a RS min = 5 (s/m 2, P a x is taken as in the EKF while P a x,λ is initially set to zero. 5. Results In the first set of experiments, only LAI is uncertain. Figure 1 displays the time series of the RMS estimation error in the four state variables for the EKF (red and the STAEKF (blue. The root-mean-square (RMS in the estimate of LAI obtained with the STAEKF is shown in the bottom-left panel, while the corresponding estimated error variance, P a LAI, in the bottom-right panel. The error is averaged over the ensemble of 1 experiments and over the ten summers. Note that in view of the rapid fluctuations of T s and T 2, to improve Copyright 212 Royal Meteorological Society Atmos. Sci. Let. 13: (212
4 STAEKF for Soil Analysis 271 Table I. RMS estimation error in the four state variables averaged over the 9 s of experiment for the Open Loop, EKF and STAEKF. T s (K T 2 (K w g (m 3 /m 3 w 2 (m 3 /m 3 LAI Open Loop EKF STAEKF Albedo Open Loop EKF STAEKF RS min Open Loop EKF STAEKF LAI and Albedo Open Loop EKF STAEKF visualization the corresponding error curves are displayed every 15 time-steps (equivalent to 12.5 h in all the figures. For surface and deep temperature, T s and T 2,the EKF and STAEKF give similar results for the first month with a slight advantage for the STAEKF at the error peaks, while a larger improvement is visible afterward. For the water contents, w g and w s, the STAEKF takes a longer time to exhibit a benefit over the EKF, but after around 2 s the STAEKF estimation error attains a lower level. After a sufficiently long time the estimation error essentially comes from the parametric error and the STAEKF becomes advantageous. The estimation of LAI with the STAEKF is also satisfactory as depicted by the error reduction observed in the bottom-left panel of Figure 1. The reduction of the actual error is tracked by the associated estimated variance, as it can be observed in the bottom-right panel. A complete perspective on the relative performance of the assimilation algorithms is given in the upper-most part of Table I which shows the RMS estimation error averaged in time (over the 9 s of experiment; results relative to an experiment without assimilation (referred to as open loop are also shown for completeness. By looking at the error level for the open loop experiment, the benefit of implementing a data assimilation technique is evident; then the slight advantage of the STAEKF over the EKF comes from its additional parameter estimation feature. After around 2 s the component of the prediction error coming from the growth of the previous analysis error, becomes progressively smaller and the error mainly originates from the misspecification of the parameters. Figure 2 is the same as Figure 1, but for the experiments with parametric error in the albedo. Similar to the previous case, the STAEKF shows better results after around 2 s, but now the improvements in the analyzed temperatures are a bit more rapid. A smooth error reduction is observed in the estimate of the albedo and is accompanied by a corresponding reduction in the associated error variance. In the second block of Table I the total RMS error relative to the experiment with error in the albedo is provided. Comparing the skill of the open loop with that obtained when LAI was uncertain, it is possible to 1.5 EKF T s T 2 STAEKF W g W Albedo x P a Albedo Figure 2. The same as in Figure 1, but for an error in the albedo. Copyright 212 Royal Meteorological Society Atmos. Sci. Let. 13: (212
5 272 A. Carrassi et al. deduce the strong impact of the albedo on the estimates of temperatures. This is also reflected by the performance of the EKF and STAEKF: the estimation errors are now larger for the temperatures and lower for the water contents, compared with those obtained in the previous case. A further indication of the different roles played by these two parameters in soil modeling is given by observing their relative performances. When LAI is uncertain the STAEKF improves over the EKF of approximately 14 and 9% in T s and T 2 but the improvement raises to 25 and 24% for w g and w 2 respectively. The situation is reversed when the albedo is uncertain: the improvement is approximately 27 and 24% in T s and T 2 in contrast to 14 and 9% for w g and w 2. In practice, reducing the parametric error leads to improving the estimate of those variables which are more sensitive to the uncertain parameters. In analogy to Figures 1 and 2, Figure 3 refers to experiments with error in RS min. In this case the overall performance of the STAEKF is less satisfactory. The estimation of RS min appears more difficult (see bottom-left panel and although a trend of error reduction is visible, it is now lower than in the two previous cases, with a parametric error reduction of about 4 and 35% in the LAI and albedo, respectively, and a reduction of 14% for RS min now. The error is also subject to significant fluctuations. As a consequence, the STAEKF is not able to reach the skill of the EKF, although it improves over the open loop (Table I. The different behavior of the STAEKF when dealing with different parameters was already observed in the study by Carrassi and Vannitsem (211a in the context of a low-order atmospheric model, and it was found connected to the sensitivity of the model dynamics to the uncertain parameters. In the STAEKF, the effect of parametric error is accounted for on the basis of a first-order approximation, a truncation which turns out to be poorly accurate in certain cases, such as the RS min here, affected by very large errors. Finally, Figure 4 describes the results relative to the experiments with the simultaneous presence of error in the LAI and the albedo. Note that now the bottom panels display the parametric error in LAI and albedo respectively. Although similar conclusions to the previous cases can be drawn, some differences are worth mentioning: the STAEKF shows now a higher skill along the entire duration of the experiment and the parametric error in LAI and albedo is successfully reduced. As it can be further deduced by looking at Table I, the simultaneous presence of errors in two parameters enhances the benefit of the parameter estimation feature of the STAEKF. Overall, a main concern revealed by these results is the slow rate of error convergence of the STAEKF. This rate is notably related to the initial parametric error variance, P a λ, which in turn affects the correlation between observables and parameters as well as the effect of parametric uncertainty on the state estimation error. After a sufficient long time the STAEKF (as all Kalman filter-type algorithms loses memory about the initialization, and the error covariances are modulated by the dynamics. In this study, P a λ has been chosen as proportional to the actual initial parametric error variance; the value used in the experiment is then the outcome of a short tuning. A more extensive tuning would have led to an acceleration of the rate.5 T s T.5 2 EKF STAEKF.1.5 W g W Minimum Stomatal Resistance (RSmin 5 P a RSmin [s/m] 25 2 [s 2 /m 2 ] Figure 3. The same as in Figure 1, but for an error in the RS min. Copyright 212 Royal Meteorological Society Atmos. Sci. Let. 13: (212
6 STAEKF for Soil Analysis 273 T s T 2.8 EKF STAEKF W g.15 W [m 2 /m 2 ].32 Leaf Area Index Albedo Figure 4. The same as in Figure 1, but for error in both the LAI and the albedo. of error convergence, but is out of the scope of this study. 6. Conclusions This paper describes a first application of the STAEKF to soil analysis, for the assimilation of screen-level observations into an off-line version of the land surface model, ISBA. A twin experiment approach has been adopted for the summers for the decade Model error is simulated by perturbing either independently or simultaneously the LAI, the albedo, and the minimum stomatal resistance, three key parameters in soil modeling, usually affected by significant errors. The effect of model error is implicitly accounted for in the EKF forecast error covariance, while in the STAEKF the uncertain parameters are estimated at the analysis time along with the system s state. Other schemes such as an EKF including a bias correction scheme also exist in the literature (Dee and da Silva, Results indicate that the STAEKF is able to retrieve the true value of the uncertain parameter and progressively reduces the associated estimation error. The accuracy of these estimates is inherently related to the type of parameter to be estimated, given the model sensitivity to the specific parameter and the accuracy of the short-time approximation at the core of the STAEKF. These conditions have to be verified case by case. The present results suggest that the on-line estimation of LAI and/or the Albedo helps in improving the soil analysis skill to a larger extent than the estimate of RS min. The rate of error convergence in the STAEKF is related to the initial parametric error variance. Such a rate could be accelerated, to a certain extent, by an optimal tuning of the initial error covariance, a case-dependent procedure which has not been implemented in this feasibility study. By reducing the parametric error a better guess for the system state can be obtained and this in turn improves the analysis field and again the accuracy of the parameter estimate. Moreover, given that this feature is incorporated using the short-time formulation (Carrassi and Vannitsem, 211a, the additional computational cost with respect to the standard EKF is almost negligible. This proof-of-concept study is a first step in the validation of the STAEKF scheme and encourages its use in a more realistic model and observational scenario, possibly in the context of a full threedimensional NWP model. This research is currently on-going. Acknowledgements The authors are thankful to J.-F. Mahfouf (Météo-France for providing the SLDAS platform and to N. Ghilain for the useful discussion on the estimation of the leaf area index and the albedo. This work has been partially supported by the Belgian Federal Science Policy Program under contract MO/34/17 and MO/34/2. References Albergel C, Calvet J-C, Mahfouf J-F, Rudiger C, Barbu AL, Lafont S, Roujean J-L, Walker JP, Crapeau M, Wigneron J-P. 21. Copyright 212 Royal Meteorological Society Atmos. Sci. Let. 13: (212
7 274 A. Carrassi et al. Monitoring of water and carbon fluxes using a land data assimilation system: a case study for southwestern France. Hydrology and Earth System Sciences. Discussion 7: Balsamo G, Bouyssel F, Noilhan J. 24. A simplified bi-dimensional variational analysis of soil moisture from screen-level observations in a mesoscale numerical weather-prediction model. Quarterly Journal of the Royal Meteorological Society 13: Bélair S, Crevier L-P, Mailhot J, Bilodeau B, Delage Y. 23. Operational implementation of the ISBA land surface scheme in the Canadian regional weather forecast model. Part I: Warm season results. Journal of Hydrometeorology 4: Carrassi A, Vannitsem S. 21. Accounting for model error in variational data assimilation. A deterministic formulation. Monthly Weather Review 138: Carrassi A, Vannitsem S. 211a. State and parameter estimation with the extended Kalman filter: An alternative formulation of the model error dynamics. Quarterly Journal of the Royal Meteorological Society 137: Carrassi A, Vannitsem S. 211b. Treatment of the error due to the unresolved scales in sequential data assimilation. International Journal of Bifurcation and Chaos 21: Dee D, da Silva A Data assimilation in the presence of forecast bias. Quarterly Journal of the Royal Meteorological Society 117: Douville H, Viterbo P, Mahfouf J-F, Beljaars ACM. 2. Evaluation of optimal interpolation and nudging techniques for soil moisture analysis using FIFE data. Monthly Weather Review 128: Draper C, Mahfouf J-F, Walker J. 29. An EKF assimilation of AMSR-E soil moisture into the ISBA land surface scheme. Journal of the Geophysical Research 114: D214, DOI: 1.129/28JD1165. Drusch M, Scipal K, de Rosnay P, Balsamo G, Andersson E, Bougeault P, Viterbo P. 29. Towards a Kalman filter based soil moisture analysis system for the operational ECMWF Integrated Forecast System. Geophysical Research Letters 36: L141. DOI: 1.129/29GL Ghilain N, Arboleda A, Sepulcre-Cantò G, Batelaan O, Ardö J, Gellens-Meulenberghs F Improving evapotranspiration in land surface models by using biophysical parameters derived from MSG/SEVIRI satellite. Hydrology and Earth System Sciences. Discussion 8: Giard D, Bazile E. 2. Implementation of a new assimilation scheme for soil and surface variables in a global NWM model. Monthly Weather Review 128: Hamdi R, Van de Vyver H, Termonia P New cloud and micro-physics parameterization for use in high-resolution dynamical downscaling: Application for summer extreme temperature over Belgium. International Journal of Climatology DOI: 1.12/joc.249. Mahfouf J-F Analysis of soil moisture from near-surface parameters: a feasibility study. Journal of Applied Meteorology 3: Mahfouf J-F. 27. Soil analysis at Météo-France, Part I: Evaluation and perspectives at local scale (in French, Note de Cent. CNRM/GMME, 84, 58 pp., Groupe de Mtorol. Moyenne Echelle, Center National de la Recherche Meteorological Service, Météo- France, Toulouse, France. Mahfouf J-F, Bliznak V Combined assimilation of screenlevel observations and radar-derived precipitation for soil moisture analysis. Quarterly Journal of the Royal Meteorological Society 137: Mahfouf J-F, Bergaoui K, Draper C, Bouyssel F, Taillefer F, Taseva L. 29. A comparison of two off-line soil analysis schemes for assimilation of screen level observations. Journal of the Geophysical Research 114: D815, DOI: 1.129/28JD1177. Nicolis C. 23. Dynamics of model error: Some generic features. Journal of the Atmospheric Sciences 6: Noilhan J, Mahfouf J-F The ISBA land-surface parameterization scheme. Global and Planetary Change 13: Orescanin B, Rajkovic B, Zupanski M, Zupanski D. 29. Soil model parameter estimation with ensemble data assimilation. Atmospheric Science Letters 1: Copyright 212 Royal Meteorological Society Atmos. Sci. Let. 13: (212
Towards a Kalman Filter based soil moisture analysis system for the operational ECMWF Integrated Forecast System
GEOPHYSICAL RESEARCH LETTERS, VOL. 36, L10401, doi:10.1029/2009gl037716, 2009 Towards a Kalman Filter based soil moisture analysis system for the operational ECMWF Integrated Forecast System M. Drusch,
More informationComparing the SEKF with the DEnKF on a land surface model
Comparing the SEKF with the DEnKF on a land surface model David Fairbairn, Alina Barbu, Emiliano Gelati, Jean-Francois Mahfouf and Jean-Christophe Caret CNRM - Meteo France Partly funded by European Union
More informationAdvances in land data assimilation at ECMWF
Advances in land data assimilation at ECMWF Patricia de Rosnay 1, Matthias Drusch 1, Gianpaolo Balsamo 1, Anton Beljaars 1, Lars Isaksen 1, Drasko Vasiljevic 1, Clément Albergel 2, Klaus Scipal 1 1 ECMWF,
More informationSummary of activities with SURFEX data assimilation at Météo-France. Jean-François MAHFOUF CNRM/GMAP/OBS
Summary of activities with SURFEX data assimilation at Météo-France Jean-François MAHFOUF CNRM/GMAP/OBS Outline Status at Météo-France in 2008 Developments undertaken during 2009-2011 : Extended Kalman
More informationThe Canadian Land Data Assimilation System (CaLDAS)
The Canadian Land Data Assimilation System (CaLDAS) Marco L. Carrera, Stéphane Bélair, Bernard Bilodeau and Sheena Solomon Meteorological Research Division, Environment Canada Dorval, QC, Canada 2 nd Workshop
More informationEvaluation of remotely sensed and modelled soil moisture products using global ground-based in situ observations
Evaluation of remotely sensed and modelled soil moisture products using global ground-based in situ observations C. Albergel (1), P. de Rosnay (1), G. Balsamo (1),J. Muñoz-Sabater(1 ), C. Gruhier (2),
More informationExtended Kalman Filter soil-moisture analysis in the IFS
from Newsletter Number 127 Spring 211 METEOROLOGY Extended Kalman Filter soil-moisture analysis in the IFS doi:1.21957/ik7co53s This article appeared in the Meteorology section of ECMWF Newsletter No.
More informationRoot zone soil moisture from the assimilation of screen level variables and remotely sensed soil moisture
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 116,, doi:10.1029/2010jd013829, 2011 Root zone soil moisture from the assimilation of screen level variables and remotely sensed soil moisture C. S. Draper, 1,2 J.
More informationEUMETSAT LSA-SAF EVAPOTRANSPIRATION PRODUCTS STATUS AND PERSPECTIVES
EUMETSAT LSA-SAF EVAPOTRANSPIRATION PRODUCTS STATUS AND PERSPECTIVES Arboleda, N. Ghilain, F. Gellens-Meulenberghs Royal Meteorological Institute, Avenue Circulaire, 3, B-1180 Bruxelles, BELGIUM Corresponding
More informationAdvances in land data assimilation at Météo-France
Advances in land data assimilation at Météo-France Jean-François Mahfouf Météo-France/CNRS - CNRM/GAME 42, Avenue Coriolis, 31057 Toulouse Cedex 01, France jean-francois.mahfouf@meteo.fr ABSTRACT This
More informationECMWF. ECMWF Land Surface Analysis: Current status and developments. P. de Rosnay M. Drusch, K. Scipal, D. Vasiljevic G. Balsamo, J.
Land Surface Analysis: Current status and developments P. de Rosnay M. Drusch, K. Scipal, D. Vasiljevic G. Balsamo, J. Muñoz Sabater 2 nd Workshop on Remote Sensing and Modeling of Surface Properties,
More information1. Current atmospheric DA systems 2. Coupling surface/atmospheric DA 3. Trends & ideas
1 Current issues in atmospheric data assimilation and its relationship with surfaces François Bouttier GAME/CNRM Météo-France 2nd workshop on remote sensing and modeling of surface properties, Toulouse,
More informationAn EKF assimilation of AMSR-E soil moisture into the ISBA land surface scheme
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 114,, doi:10.1029/2008jd011650, 2009 An EKF assimilation of AMSR-E soil moisture into the ISBA land surface scheme C. S. Draper, 1 J.-F. Mahfouf, 2 and J. P. Walker
More informationCombined assimilation of screen-level observations and radar-derived precipitation for soil moisture analysis
Quarterly Journal of the Royal Meteorological Society Q. J. R. Meteorol. Soc. 137: 709 722, April 2011 A Combined assimilation of screen-level observations and radar-derived precipitation for soil moisture
More informationECMWF soil moisture data assimilation
ECMWF soil moisture data assimilation Janneke Ettema and Pedro Viterbo ECMWF, Shinfield Park, Reading RG2 9AX, United Kingdom J.Ettema@ECMWF.int, P.Viterbo@ECMWF.int 1. Introduction Soil moisture is a
More informationECMWF. ECMWF Land Surface modelling and land surface analysis. P. de Rosnay G. Balsamo S. Boussetta, J. Munoz Sabater D.
Land Surface modelling and land surface analysis P. de Rosnay G. Balsamo S. Boussetta, J. Munoz Sabater D. Vasiljevic M. Drusch, K. Scipal SRNWP 12 June 2009 Slide 1 Surface modelling (G. Balsamo) HTESSEL,
More informationDevelopment of the Canadian Precipitation Analysis (CaPA) and the Canadian Land Data Assimilation System (CaLDAS)
Development of the Canadian Precipitation Analysis (CaPA) and the Canadian Land Data Assimilation System (CaLDAS) Marco L. Carrera, Vincent Fortin and Stéphane Bélair Meteorological Research Division Environment
More informationLand surface data assimilation for Numerical Weather Prediction
Sixth WMO Symposium on Data Assimilation, University of Maryland, 7-11 October 2013 Land surface data assimilation for Numerical Weather Prediction P. de Rosnay, J. Muñoz Sabater, C. Albergel, G. Balsamo,
More informationSoil Moisture Prediction and Assimilation
Soil Moisture Prediction and Assimilation Analysis and Prediction in Agricultural Landscapes Saskatoon, June 19-20, 2007 STEPHANE BELAIR Meteorological Research Division Prediction and Assimilation Atmospheric
More informationIntroduction to Data Assimilation
Introduction to Data Assimilation Alan O Neill Data Assimilation Research Centre University of Reading What is data assimilation? Data assimilation is the technique whereby observational data are combined
More informationMultiscale performance of the ALARO-0 model for simulating extreme summer precipitation climatology in Belgium
Multiscale performance of the ALARO-0 model for simulating extreme summer precipitation climatology in Belgium Rozemien De Troch 1,2, Rafiq Hamdi 1, Hans Van de Vyver 1, Jean-François Geleyn 2,3, Piet
More informationDevelopment of a land data assimilation system at NILU
Development of a land data assimilation system at NILU W.A. Lahoz, wal@nilu.no 10 th June 2009, 2 nd Workshop on Remote Sensing and Modeling of Surface Properties, Toulouse, France Thanks to S.-E. Walker
More informationAssimilation of land surface satellite data for Numerical Weather Prediction at ECMWF
4 th workshop on Remote Sensing and Modelling of Surface properties Saint Martin d Hères, 14-16 March 2016 Assimilation of land surface satellite data for Numerical Weather Prediction at ECMWF P. de Rosnay,
More informationModelling and Data Assimilation Needs for improving the representation of Cold Processes at ECMWF
Modelling and Data Assimilation Needs for improving the representation of Cold Processes at ECMWF presented by Gianpaolo Balsamo with contributions from Patricia de Rosnay, Richard Forbes, Anton Beljaars,
More informationMeteo-France operational land surface analysis for NWP: current status and perspectives
Meteo-France operational land surface analysis for NWP: current status and perspectives Camille Birman, Jean-François Mahfouf, Yann Seity, Eric Bazile, Françoise Taillefer, Zied Sassi, Nadia Fourrié, Vincent
More informationSnow Analysis for Numerical Weather prediction at ECMWF
Snow Analysis for Numerical Weather prediction at ECMWF Patricia de Rosnay, Gianpaolo Balsamo, Lars Isaksen European Centre for Medium-Range Weather Forecasts IGARSS 2011, Vancouver, 25-29 July 2011 Slide
More informationSoil moisture analysis combining screen-level parameters and microwave brightness temperature: A test with field data
403 Soil moisture analysis combining screen-level parameters and microwave brightness temperature: A test with field data G. Seuffert, H.Wilker 1, P. Viterbo, J.-F. Mahfouf 2, M. Drusch, J.-C. Calvet 3
More informationAssimilation of satellite derived soil moisture for weather forecasting
Assimilation of satellite derived soil moisture for weather forecasting www.cawcr.gov.au Imtiaz Dharssi and Peter Steinle February 2011 SMOS/SMAP workshop, Monash University Summary In preparation of the
More informationHigh resolution land reanalysis
Regional Reanalysis Workshop 19-20 May 2016, Reading High resolution land reanalysis P. de Rosnay, G. Balsamo, J. Muñoz Sabater, E. Dutra, C. Albergel, N. Rodríguez-Fernández, H. Hersbach Introduction:
More informationLand Data Assimilation for operational weather forecasting
Land Data Assimilation for operational weather forecasting Brett Candy Richard Renshaw, JuHyoung Lee & Imtiaz Dharssi * *Centre Australian Weather and Climate Research Contents An overview of the Current
More informationEuropean High-Resolution Soil Moisture Analysis (EHRSOMA)
European High-Resolution Soil Moisture Analysis (EHRSOMA) Jasmin Vural EUMETSAT Fellow Day, 05.03.2018 European High-Resolution Soil Moisture Analysis (EHRSOMA) Jasmin Vural EUMETSAT Fellow Day, 05.03.2018
More informationRadar data assimilation using a modular programming approach with the Ensemble Kalman Filter: preliminary results
Radar data assimilation using a modular programming approach with the Ensemble Kalman Filter: preliminary results I. Maiello 1, L. Delle Monache 2, G. Romine 2, E. Picciotti 3, F.S. Marzano 4, R. Ferretti
More informationESA CONTRACT REPORT. Joaquín Muñoz Sabater, Patricia de Rosnay, Clément Albergel, Lars Isaksen
ESA CONTRACT REPORT Contract Report to the European Space Agency TNPII WP3401 & WP3402 SMOS Report on background and observation error scenarios Joaquín Muñoz Sabater, Patricia de Rosnay, Clément Albergel,
More informationERA-5 driven land surface reanalysis : LDAS-Monde applied to the Continental US
ERA-5 driven land surface reanalysis : LDAS-Monde applied to the Continental US Clement Albergel 1, Emanuel Dutra 2, Simon Munier 1, Jean-Christophe Calvet 1, Joaquin Munoz-Sabater 3, Patricia de Rosnay
More informationA simplified Extended Kalman Filter for the global operational soil moisture analysis at ECMWF
Quarterly Journal of the Royal Meteorological Society Q. J. R. Meteorol. Soc. 139: 1199 1213, July 2013 A A simplified Extended Kalman Filter for the global operational soil moisture analysis at ECMWF
More informationCurrent Limited Area Applications
Current Limited Area Applications Nils Gustafsson SMHI Norrköping, Sweden nils.gustafsson@smhi.se Outline of talk (contributions from many HIRLAM staff members) Specific problems of Limited Area Model
More informationStochastic methods for representing atmospheric model uncertainties in ECMWF's IFS model
Stochastic methods for representing atmospheric model uncertainties in ECMWF's IFS model Sarah-Jane Lock Model Uncertainty, Research Department, ECMWF With thanks to Martin Leutbecher, Simon Lang, Pirkka
More informationAssimilation of ASCAT soil wetness
EWGLAM, October 2010 Assimilation of ASCAT soil wetness Bruce Macpherson, on behalf of Imtiaz Dharssi, Keir Bovis and Clive Jones Contents This presentation covers the following areas ASCAT soil wetness
More informationProspects for radar and lidar cloud assimilation
Prospects for radar and lidar cloud assimilation Marta Janisková, ECMWF Thanks to: S. Di Michele, E. Martins, A. Beljaars, S. English, P. Lopez, P. Bauer ECMWF Seminar on the Use of Satellite Observations
More informationThe Canadian approach to ensemble prediction
The Canadian approach to ensemble prediction ECMWF 2017 Annual seminar: Ensemble prediction : past, present and future. Pieter Houtekamer Montreal, Canada Overview. The Canadian approach. What are the
More informationInvestigating the urban climate characteristics of two Hungarian cities with SURFEX/TEB land surface model
Investigating the urban climate characteristics of two Hungarian cities with SURFEX/TEB land surface model Gabriella Zsebeházi Gabriella Zsebeházi and Gabriella Szépszó Hungarian Meteorological Service,
More informationThe hybrid ETKF- Variational data assimilation scheme in HIRLAM
The hybrid ETKF- Variational data assimilation scheme in HIRLAM (current status, problems and further developments) The Hungarian Meteorological Service, Budapest, 24.01.2011 Nils Gustafsson, Jelena Bojarova
More informationTowards a probabilistic hydrological forecasting and data assimilation system. Henrik Madsen DHI, Denmark
Towards a probabilistic hydrological forecasting and data assimilation system Henrik Madsen DHI, Denmark Outline Hydrological forecasting Data assimilation framework Data assimilation experiments Concluding
More informationEnKF Review. P.L. Houtekamer 7th EnKF workshop Introduction to the EnKF. Challenges. The ultimate global EnKF algorithm
Overview 1 2 3 Review of the Ensemble Kalman Filter for Atmospheric Data Assimilation 6th EnKF Purpose EnKF equations localization After the 6th EnKF (2014), I decided with Prof. Zhang to summarize progress
More informationOutline. Part I (Monday) Part II (Tuesday) ECMWF. Introduction Snow analysis Screen level parameters analysis
Outline Part I (Monday) Introduction Snow analysis Screen level parameters analysis Part II (Tuesday) Soil moisture analysis OI and EKF analyses Use of satellite data: ASCAT and SMOS Summary and future
More informationImpact of improved soil moisture on the ECMWF precipitation forecast in West Africa
GEOPHYSICAL RESEARCH LETTERS, VOL. 37,, doi:10.1029/2010gl044748, 2010 Impact of improved soil moisture on the ECMWF precipitation forecast in West Africa A. Agustí Panareda, 1 G. Balsamo, 1 and A. Beljaars
More information4DEnVar. Four-Dimensional Ensemble-Variational Data Assimilation. Colloque National sur l'assimilation de données
Four-Dimensional Ensemble-Variational Data Assimilation 4DEnVar Colloque National sur l'assimilation de données Andrew Lorenc, Toulouse France. 1-3 décembre 2014 Crown copyright Met Office 4DEnVar: Topics
More informationThe Usage of Screen-Level Parameters and Microwave Brightness Temperature for Soil Moisture Analysis
516 JOURNAL OF HYDROMETEOROLOGY VOLUME 5 The Usage of Screen-Level Parameters and Microwave Brightness Temperature for Soil Moisture Analysis G. SEUFFERT European Centre for Medium-Range Weather Forecasts,
More informationA SATELLITE LAND DATA ASSIMILATION SYTEM COUPLED WITH A MESOSCALE MODEL: TOWARDS IMPROVING NUMERICAL WEATHER PREDICTION
A SATELLITE LAND DATA ASSIMILATION SYTEM COUPLED WITH A MESOSCALE MODEL: TOWARDS IMPROVING NUMERICAL WEATHER PREDICTION Mohamed Rasmy*, Toshio Koike*, Souhail Bousseta**, Xin Li*** Dept. of Civil Engineering,
More informationEnsemble of Data Assimilations methods for the initialization of EPS
Ensemble of Data Assimilations methods for the initialization of EPS Laure RAYNAUD Météo-France ECMWF Annual Seminar Reading, 12 September 2017 Introduction Estimating the uncertainty in the initial conditions
More informationA simplified Extended Kalman Filter for the global operational soil moisture analysis at ECMWF
662 A simplified Extended Kalman Filter for the global operational soil moisture analysis at ECMWF P. de Rosnay, M. Drusch, D. Vasiljevic, G. Balsamo, C. Albergel and L. Isaksen Research Department Submitted
More informationLand Surface Data Assimilation
Data Assimilation Training course Land Surface Data Assimilation Patricia de Rosnay Training course - Surface analysis Part I 10 March 2014 Slide 1 Outline Part I (Monday 10 March) Introduction Snow analysis
More informationConvective scheme and resolution impacts on seasonal precipitation forecasts
GEOPHYSICAL RESEARCH LETTERS, VOL. 30, NO. 20, 2078, doi:10.1029/2003gl018297, 2003 Convective scheme and resolution impacts on seasonal precipitation forecasts D. W. Shin, T. E. LaRow, and S. Cocke Center
More informationThe Planetary Boundary Layer and Uncertainty in Lower Boundary Conditions
The Planetary Boundary Layer and Uncertainty in Lower Boundary Conditions Joshua Hacker National Center for Atmospheric Research hacker@ucar.edu Topics The closure problem and physical parameterizations
More informationLinearized Physics (LP): Progress and Issues
Linearized Physics (LP): Progress and Issues Philippe Lopez ECMWF, Reading, UK with thanks to Marta Janisková, Alan Geer and Peter Bauer ECMWF-JCSDA Workshop on Assimilating Satellite Observations of Clouds
More informationRelative Merits of 4D-Var and Ensemble Kalman Filter
Relative Merits of 4D-Var and Ensemble Kalman Filter Andrew Lorenc Met Office, Exeter International summer school on Atmospheric and Oceanic Sciences (ISSAOS) "Atmospheric Data Assimilation". August 29
More informationSurface issues for High Resolution NWP
Surface issues for High Resolution NWP François Bouyssel - Jean-François Mahfouf with many contributors High Resolution NWP Workshop, 17-20 May 2010, Brac, Croatia Outline Introduction (specificities of
More informationEXPERIMENTAL ASSIMILATION OF SPACE-BORNE CLOUD RADAR AND LIDAR OBSERVATIONS AT ECMWF
EXPERIMENTAL ASSIMILATION OF SPACE-BORNE CLOUD RADAR AND LIDAR OBSERVATIONS AT ECMWF Marta Janisková, Sabatino Di Michele, Edouard Martins ECMWF, Shinfield Park, Reading, U.K. Abstract Space-borne active
More informationSENSITIVITY STUDY FOR SZEGED, HUNGARY USING THE SURFEX/TEB SCHEME COUPLED TO ALARO
SENSITIVITY STUDY FOR SZEGED, HUNGARY USING THE SURFEX/TEB SCHEME COUPLED TO ALARO Report from the Flat-Rate Stay at the Royal Meteorological Institute, Brussels, Belgium 11.03.2015 01.04.2015 Gabriella
More informationVariational data assimilation of lightning with WRFDA system using nonlinear observation operators
Variational data assimilation of lightning with WRFDA system using nonlinear observation operators Virginia Tech, Blacksburg, Virginia Florida State University, Tallahassee, Florida rstefane@vt.edu, inavon@fsu.edu
More informationThe role of soil moisture in influencing climate and terrestrial ecosystem processes
1of 18 The role of soil moisture in influencing climate and terrestrial ecosystem processes Vivek Arora Canadian Centre for Climate Modelling and Analysis Meteorological Service of Canada Outline 2of 18
More informationTangent-linear and adjoint models in data assimilation
Tangent-linear and adjoint models in data assimilation Marta Janisková and Philippe Lopez ECMWF Thanks to: F. Váňa, M.Fielding 2018 Annual Seminar: Earth system assimilation 10-13 September 2018 Tangent-linear
More informationMesoscale meteorological models. Claire L. Vincent, Caroline Draxl and Joakim R. Nielsen
Mesoscale meteorological models Claire L. Vincent, Caroline Draxl and Joakim R. Nielsen Outline Mesoscale and synoptic scale meteorology Meteorological models Dynamics Parametrizations and interactions
More informationRemote sensing derived evapotranspiration: comparisons to observations and models and results over the full MSG disk and selected basins
4 th LSA-SAF WORKSHOP Toulouse,France, 15 th to 17 th November 2010 Remote sensing derived evapotranspiration: comparisons to observations and models and results over the full MSG disk and selected basins
More informationCross-evaluation of modelled and remotely sensed surface soil moisture with in situ data in southwestern France
doi:10.5194/hess-14-2177-2010 Author(s) 2010. CC Attribution 3.0 License. Hydrology and Earth System Sciences Cross-evaluation of modelled and remotely sensed surface soil moisture with in situ data in
More informationA comparison of two off-line soil analysis schemes for assimilation of screen level observations
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 114,, doi:10.1029/2008jd011077, 2009 A comparison of two off-line soil analysis schemes for assimilation of screen level oservations J.-F. Mahfouf, 1 K. Bergaoui,
More informationGaussian Process Approximations of Stochastic Differential Equations
Gaussian Process Approximations of Stochastic Differential Equations Cédric Archambeau Dan Cawford Manfred Opper John Shawe-Taylor May, 2006 1 Introduction Some of the most complex models routinely run
More informationDrought Monitoring with Hydrological Modelling
st Joint EARS/JRC International Drought Workshop, Ljubljana,.-5. September 009 Drought Monitoring with Hydrological Modelling Stefan Niemeyer IES - Institute for Environment and Sustainability Ispra -
More informationTHE OFFLINE VERSION OF SURFEX COUPLED TO THE
SURFEX@RMI THE OFFLINE VERSION OF SURFEX COUPLED TO THE OPERATIONAL ALADIN FORECAST OVER BELGIUM: A TOOL TO IMPROVE WINTER SCREEN TEMPERATURE Royal Meteorological Institute, Brussels, Belgium Corresponding
More informationLATE REQUEST FOR A SPECIAL PROJECT
LATE REQUEST FOR A SPECIAL PROJECT 2016 2018 MEMBER STATE: Italy Principal Investigator 1 : Affiliation: Address: E-mail: Other researchers: Project Title: Valerio Capecchi LaMMA Consortium - Environmental
More informationData assimilation in mesoscale modeling and numerical weather prediction Nils Gustafsson
Data assimilation in mesoscale modeling and numerical weather prediction Nils Gustafsson Croatian USA Workshop on Mesometeorology June 2012 Perspective: What are the important issues for development of
More informationForecasting and data assimilation
Supported by the National Science Foundation DMS Forecasting and data assimilation Outline Numerical models Kalman Filter Ensembles Douglas Nychka, Thomas Bengtsson, Chris Snyder Geophysical Statistics
More informationNew soil physical properties implemented in the Unified Model
New soil physical properties implemented in the Unified Model Imtiaz Dharssi 1, Pier Luigi Vidale 3, Anne Verhoef 3, Bruce Macpherson 1, Clive Jones 1 and Martin Best 2 1 Met Office (Exeter, UK) 2 Met
More informationProf. Stephen G. Penny University of Maryland NOAA/NCEP, RIKEN AICS, ECMWF US CLIVAR Summit, 9 August 2017
COUPLED DATA ASSIMILATION: What we need from observations and modellers to make coupled data assimilation the new standard for prediction and reanalysis. Prof. Stephen G. Penny University of Maryland NOAA/NCEP,
More informationDynamic Inference of Background Error Correlation between Surface Skin and Air Temperature
Dynamic Inference of Background Error Correlation between Surface Skin and Air Temperature Louis Garand, Mark Buehner, and Nicolas Wagneur Meteorological Service of Canada, Dorval, P. Quebec, Canada Abstract
More informationRecent Data Assimilation Activities at Environment Canada
Recent Data Assimilation Activities at Environment Canada Major upgrade to global and regional deterministic prediction systems (now in parallel run) Sea ice data assimilation Mark Buehner Data Assimilation
More informationLAND SURFACE ALBEDO AND DOWNWELLING SHORTWAVE RADIATION FROM MSG GEOSTATIONARY SATELLITE (METHOD FOR RETRIEVAL, VALIDATION, AND APPLICATION)
LAND SURFACE ALBEDO AND DOWNWELLING SHORTWAVE RADIATION FROM MSG GEOSTATIONARY SATELLITE (METHOD FOR RETRIEVAL, VALIDATION, AND APPLICATION) Dominique Carrer, Jean-Louis Roujean, Olivier Hautecoeur, Jean-Christophe
More informationSensitivity of NWP model skill to the obliquity of the GPS radio occultation soundings
ATMOSPHERIC SCIENCE LETTERS Atmos. Sci. Let. 13: 55 60 (2012) Published online 1 November 2011 in Wiley Online Library (wileyonlinelibrary.com) DOI: 10.1002/asl.363 Sensitivity of NWP model skill to the
More informationLand Surface: Snow Emanuel Dutra
Land Surface: Snow Emanuel Dutra emanuel.dutra@ecmwf.int Slide 1 Parameterizations training course 2015, Land-surface: Snow ECMWF Outline Snow in the climate system, an overview: Observations; Modeling;
More informationEvapotranspiration monitoring with Meteosat Second Generation satellites: method, products and utility in drought detection.
Evapotranspiration monitoring with Meteosat Second Generation satellites: method, products and utility in drought detection. Nicolas Ghilain Royal Meteorological Institute Belgium EUMeTrain Event week
More informationMétéo-France seasonal forecast system 5 versus system 4
Météo-France seasonal forecast system 5 versus system 4 Robust scores June 2015 Page 1 Table of contents 1. Introduction...2 2. ENSO Scores...2 3. Anomaly correlations...3 4. Circulation indices...4 5.
More informationERA-CLIM: Developing reanalyses of the coupled climate system
ERA-CLIM: Developing reanalyses of the coupled climate system Dick Dee Acknowledgements: Reanalysis team and many others at ECMWF, ERA-CLIM project partners at Met Office, Météo France, EUMETSAT, Un. Bern,
More informationRobust Ensemble Filtering With Improved Storm Surge Forecasting
Robust Ensemble Filtering With Improved Storm Surge Forecasting U. Altaf, T. Buttler, X. Luo, C. Dawson, T. Mao, I.Hoteit Meteo France, Toulouse, Nov 13, 2012 Project Ensemble data assimilation for storm
More informationRegional Climate Simulations with WRF Model
WDS'3 Proceedings of Contributed Papers, Part III, 8 84, 23. ISBN 978-8-737852-8 MATFYZPRESS Regional Climate Simulations with WRF Model J. Karlický Charles University in Prague, Faculty of Mathematics
More informationTHORPEX Data Assimilation and Observing Strategies Working Group: scientific objectives
THORPEX Data Assimilation and Observing Strategies Working Group: scientific objectives Pierre Gauthier Co-chair of the DAOS-WG Atmospheric Science and Technology Branch Environment Canada Data assimilation
More informationImpact of Surface Parameter Uncertainties within the Canadian Regional Ensemble Prediction System
1506 M O N T H L Y W E A T H E R R E V I E W VOLUME 141 Impact of Surface Parameter Uncertainties within the Canadian Regional Ensemble Prediction System CHRISTOPHE LAVAYSSE Atmospheric and Oceanic Sciences,
More informationREVISION OF THE STATEMENT OF GUIDANCE FOR GLOBAL NUMERICAL WEATHER PREDICTION. (Submitted by Dr. J. Eyre)
WORLD METEOROLOGICAL ORGANIZATION Distr.: RESTRICTED CBS/OPAG-IOS (ODRRGOS-5)/Doc.5, Add.5 (11.VI.2002) COMMISSION FOR BASIC SYSTEMS OPEN PROGRAMME AREA GROUP ON INTEGRATED OBSERVING SYSTEMS ITEM: 4 EXPERT
More informationForest Processes and Land Surface Tiling in the ECMWF model
Forest Processes and Land Surface Tiling in the ECMWF model presented by Gianpaolo Balsamo Outline: Introduction on HTESSEL and its evolution Tiling and forest/snow contrasts Tiling and forest/lakes contrasts
More informationway and atmospheric models
Scale-consistent consistent two-way way coupling of land-surface and atmospheric models COSMO-User-Seminar 9-11 March 2009 Annika Schomburg, Victor Venema, Felix Ament, Clemens Simmer TR / SFB 32 Objective
More informationAn Investigation of Reforecasting Applications for NGGPS Aviation Weather Prediction: An Initial Study of Ceiling and Visibility Prediction
An Investigation of Reforecasting Applications for NGGPS Aviation Weather Prediction: An Initial Study of Ceiling and Visibility Prediction Kathryn L. Verlinden, Oregon State University David Bright, WFO
More informationREQUEST FOR A SPECIAL PROJECT
REQUEST FOR A SPECIAL PROJECT 2017 2019 MEMBER STATE: Sweden.... 1 Principal InvestigatorP0F P: Wilhelm May... Affiliation: Address: Centre for Environmental and Climate Research, Lund University Sölvegatan
More informationGuo-Yuan Lien*, Eugenia Kalnay, and Takemasa Miyoshi University of Maryland, College Park, Maryland 2. METHODOLOGY
9.2 EFFECTIVE ASSIMILATION OF GLOBAL PRECIPITATION: SIMULATION EXPERIMENTS Guo-Yuan Lien*, Eugenia Kalnay, and Takemasa Miyoshi University of Maryland, College Park, Maryland 1. INTRODUCTION * Precipitation
More informationApplication and verification of ECMWF products 2009
Application and verification of ECMWF products 2009 Hungarian Meteorological Service 1. Summary of major highlights The objective verification of ECMWF forecasts have been continued on all the time ranges
More informationThe Impact of Background Error on Incomplete Observations for 4D-Var Data Assimilation with the FSU GSM
The Impact of Background Error on Incomplete Observations for 4D-Var Data Assimilation with the FSU GSM I. Michael Navon 1, Dacian N. Daescu 2, and Zhuo Liu 1 1 School of Computational Science and Information
More informationData Assimilation: Finding the Initial Conditions in Large Dynamical Systems. Eric Kostelich Data Mining Seminar, Feb. 6, 2006
Data Assimilation: Finding the Initial Conditions in Large Dynamical Systems Eric Kostelich Data Mining Seminar, Feb. 6, 2006 kostelich@asu.edu Co-Workers Istvan Szunyogh, Gyorgyi Gyarmati, Ed Ott, Brian
More informationConvective-scale NWP for Singapore
Convective-scale NWP for Singapore Hans Huang and the weather modelling and prediction section MSS, Singapore Dale Barker and the SINGV team Met Office, Exeter, UK ECMWF Symposium on Dynamical Meteorology
More informationSnow-atmosphere interactions at Dome C, Antarctica
Snow-atmosphere interactions at Dome C, Antarctica Éric Brun, Vincent Vionnet CNRM/GAME (Météo-France and CNRS) Christophe Genthon, Delphine Six, Ghislain Picard LGGE (CNRS and UJF)... and many colleagues
More informationComparison results: time series Margherita Grossi
Comparison results: time series Margherita Grossi GOME Evolution Climate Product v2.01 vs. ECMWF ERAInterim GOME Evolution Climate Product v2.01 vs. SSM/I HOAPS4 In order to assess the quality and stability
More informationLocal Ensemble Transform Kalman Filter
Local Ensemble Transform Kalman Filter Brian Hunt 11 June 2013 Review of Notation Forecast model: a known function M on a vector space of model states. Truth: an unknown sequence {x n } of model states
More informationClimate Models and Snow: Projections and Predictions, Decades to Days
Climate Models and Snow: Projections and Predictions, Decades to Days Outline Three Snow Lectures: 1. Why you should care about snow 2. How we measure snow 3. Snow and climate modeling The observational
More information