1 / 37

François Bouyssel

Presentation of surface analyses (temperature and soil moiture over continental area, sea surface temperature, sea ice cover, snow). François Bouyssel with inputs from G. Balsamo, E. Bazile, K. Bergaoui, D. Giard, M. Harrouche, S. Ivatek-Sahdan, M. Jidane, F. Taillefer, L. Taseva,.

risa
Download Presentation

François Bouyssel

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Presentation of surface analyses (temperature and soil moiture over continental area, sea surface temperature, sea ice cover, snow) François Bouyssel with inputs from G. Balsamo, E. Bazile, K. Bergaoui, D. Giard, M. Harrouche, S. Ivatek-Sahdan, M. Jidane, F. Taillefer, L. Taseva, ... 12-14/11/2007 Budapest HIRLAM / AAA workshop on surface assimilation

  2. Plan • Introduction • Soil moisture and soil temperature analysis • Others analyses: snow, sea surface temperature, sea ice, ...

  3. Introduction

  4. Surface fluxes : key role in the evolution of meteorological fields near the ground, in the boundary layer and in the troposphere These fluxes depend strongly on surface variables which have strong variabilities in time and space (pronostic variables)  Necessity of same degree of sophistication between surface scheme, physiographic database, surface analysis Surface analyses are performed separately from upper air analysis Several surface analyses are used for different surface parameters (Soil temperature and Soil moisture, Snow, SST, Sea ice, ...) +6hfcst +6hfcst +6hfcst +6hfcst Introduction 00UTC 06UTC 12UTC 18UTC 00UTC

  5. For the time being surface analyses are performed separately from upper air analysis. In theory a single analysis would be better but it is much more difficult implement: 1) definition of B between upper air and surface variables, 2) time scale evolutions may be different, ... For the time being several surface analyses are used for simplicity and because very different surface parameters (Soil temperature and Soil moisture, Snow, SST, Sea ice, ...) +6hfcst +6hfcst +6hfcst +6hfcst Surface analyses and upper-air analysis Atmospheric analysis and several surface analyses are done separately and combined at the end to provide the final analysis for the forecast 00UTC 06UTC 12UTC 18UTC 00UTC (6h sequantial analysis is just an exemple)

  6. Soil moisture and soil temperature analysis

  7. Surface Parameterization scheme (ISBA) Operational version : Noilhan & Planton (1989), Noilhan & Mahfouf (1996), Bazile (1999), Giard & Bazile (2000) Water Energy Analysis surface temperature mean soil temperature superficial soil water content total soil water content ~1-2 h ~1-2 day ~6-12 h ~10 days Research versions :interactive vegetation module (Calvet et al. 1998), sub grid-scale runoff and sub-root layer (Boone et al 1999), explicit 3-layers snow scheme (Boone & Etchevers 2001), tiling, multi-layer soil scheme, urban scheme

  8. The link between soil moisture and atmosphere • The main interaction of soil moisture and atmosphere is due to evaporation and vegetation transpiration processes. E 0 <SWI< 1 Eg Etr Ws Wp Wpvegetation Ws bare ground

  9. stable surface conditions : Low surface fluxes. Influence of surface limited near the ground instable surface conditions : Strong surface fluxes. Influence on PBL evolution and sometimes more (trigger deep convection) Soil moisture is very important under strong solar radiation at the surface because it determines the repartition of incoming energy into sensible and latent heat fluxes. Importance of initialization: Wr << Ws << Wp according capacity and time scale evolution. Accumulation of model error may degrade significantly the forecast during long period. Soil temperature is important in case of stable conditions because it affects low level temperature. Importance of initialization: Ts << Tp Necessity of same degree of sophistication between surface scheme, physiographic database, surface analysis Importance of soil moisture and temperature analysis

  10. 10 days 5 days 48 h 24 h 12 h 6 h t t0 T2m forecast error But strong sensibilityof surface fluxes (sensible and latent) to the soil temperature and moisture To illustrate the memory effect, the impact of aprescribed initial error in the soil moisture fieldis shown for different forecast ranges (T2m, RH2m) RH2m forecast error Wp initial error

  11. Specificities of soil moisture and temperature analysis • Strong soil and vegetation spatial heterogeneities (mountains, coastal regions, forest, bare ground, various cultures, towns, lakes, ...) • Strong spatial variability of soil moisture (linked with surface and soil properties and precipitations) • Lack of direct observations (very expensive and problem of representativeness) • Large variety of time scales in soil processes (up to several weeks or months)

  12. Precipitations observations (rain gauges, radars) : + direct link with the variations of soil water content Satellite observations: + global coverage + infrared: clear sky, low vegetation, geostationnary satellites : high temporal and spatial resolutions (energy budget), strong sensitivity to low level wind, surface roughness + microwave: active and passive instruments measure directly the soil moisture in the first few centimeters (scatterometer (ERS,ASCAT), passive or active radiameters (SMOS, HYDROS): resolution ~20/40km, frequency ~0.3/1 per day 2m observations (temperature et humidity): + good global coverage of existing network + close links with the fields in the ground in some meteorological conditions Available observations for soil moisture analysis

  13. Climatological relaxation of deep soil parameters (uncertainties in these climatologies (GSWP), interannual variability not taken into account) Off line surface scheme driven by forecasted or analysed fields and fluxes (flux of precipitation, of radiation, fields near the surface T2m, HU2m, V10m, Ps) Exemple: SAFRAN-ISBA-MODCOU Few operational use of satellite data for temperature and soil moisture analysis (near future) Assimilation of 2m observations of temperature and relative humidity Operational initialization methods

  14. Off-line method (SIM exemple) Run operationally over France at 8 km : SAFRAN (upperair analysis:Ta, qa, U, SW, LW , RR, …) , ISBA, MODCOUHydrological model

  15. Off-line method (SIM exemple)

  16. SMOSREX SIM Off-line method (SIM exemple) • Validation: river flow & snow depth & measurement site water table (Seine) Soil Wetness Index(SMOSREX)

  17. Strengths: + with good precipitation, radiation and atmospheric forcings provides realistic soil moisture evolution even at high temporal evolution (useful for NWP, but also agriculture, water managment, ...) + cheap model (just the surface), work on PC, allows multi-years reanalysis + allows the use of complex surface model + high spatial resolution (RR analysis, MSG radiation fluxes) Limitations: + no analysis & perfect model hypothesis while surface processes are complex and physiographic database not perfect: model bias may exist on soil moisture and soil temperature and remain for a long period + restricted to some geographical areas (good obs coverage) Off-line method (pros & cons)

  18. Assimilation of 2m observations

  19. T2m t RH2m t Wp t 6-h 12-h 18-h 0-h Optimum Interpolation methodCoiffier 1987, Mahfouf 1991, Bouttier 1993, Giard and Bazile 2000 1) Optimum Interpolation of T2m and RH2m using SYNOP observations interpolated at the model grid-point (by a 2m analysis) D T2m = T2ma - T2mb D RH2m = RH2ma - RH2mb 2) Correction of surface parameters (Ts, Tp, Ws, Wp) using 2m increments between analysed and forecasted values Sequential analysis (every 6h) xa = xb + BHT(HBHT + R)-1(y - H(xb)) Tsa - Tsb = D T2m Tpa - Tpb = D T2m / 2p Wsa - Wsb = aWsT D T2m + aWsRHD RH2m Wpa - Wpb = aWpT D T2m + aWpRHD RH2m OI coefficients

  20. Optimum Interpolation coefficients Very strong dependency of these backgroung error statistics to physiographic properties and meteorological conditions MonteCarlo method under summer anticyclonic conditions to get the dependency to physiography (deriving analytical formulation of OI coefficients) + empirical additional dependency to meteorological conditions Long and difficult work (in principle should be redo with model or physiography evolutions!) aWp/sT/RH = f (t, veg, LAI/Rsmin, texture, atmospheric conditions)

  21. Soil moisture analysis From Giard and Bazile 2000

  22. Soil moisture analysis • March 98: - Operational implementation with ISBA (Giard and Bazile, 2001) • Characteristics: - OI coefficients according Bouttier et al. (1993) - Reduction of OI coefficients under specific meteorological conditions (see below) - Analysis is not allowed to make Wp/Ws jump outside the range : Veg.Wwilt < Wp < Wfc and 0 < Ws < Wfc (LIMVEG=T and LHUMID=T) - Temporal smoothing of total soil moisture increments (LISSEW=T) DWp = 0.25* { DWp(HH) + DWp(HH-6) + DWp(HH-12) + DWp(HH-18) } - Removing bias on DT2m analysis increments DT2m* = (1-SCOEFT).DT2m* + SCOEFT.DT2m with SCOEFT=0.5 DT2m = DT2m - DT2m* - Relaxation towards a climatology for Tp, Wp, Sn |Model Fields Threshold | Min solar time duration J_min 6 hMax wind velocity Vmax 10 ms-1Max precipitation P_max 0.3 mmMin surface evaporation E_min 0.001 mmMax soil ice W_imax 5.0 mm Presence of snow Sn_max 0.001 kg/m2

  23. Soil moisture analysis • October 99: - Factor 3 reduction of OI coefficients on Wp because the initial OI coeff have been computed with an NMC method using initial Wp values ranging uniformly between 0 < SWI < 1 but without any rescaling of sWp~0.25.(Wpfc-Wpwilt) - Continuous formulations for OI coefficients - Cloudiness is taken into account in OI coefficients New with factor 3 reduction OI Lonnberg-Holl. method

  24. Soil moisture analysis • May 03: - Spatial smoothing of Soil Wetness Index (SWI) - Improved 2m background error statistics (smaller scales) - Factor 2 reduction of OI coefficients on Wp - Zenith solar angle is taken into account - Remove temporal smooting of Wp analysis increments - No bias correction on T2m analysis increments • Improvments of SYNOP scores on T2m and H2m in summer • More realistic soil moisture T2m H2m OPER /NEW

  25. Illustration of problem with first implementation:42h ALADIN forecast for 17th June 2000 at 18h UTC Humide Sec

  26. Smoothing of Soil Wetness Index (SWI) - II SWI for 19th June 2002 00 UTC - summer example

  27. Soil wetness index (SWI) pour le 2 mai 2004 OPER NEW

  28. Soil Wetness Index in SIM (left) et in ARPEGE (right)11 July 2005

  29. Analysis increments (May-June 2006) Daily mean of absolute analysis increments |DT2m| Cumulated analysis increments on Wp (in mm)

  30. Comparison of statistical and dynamical OI • masking of low sensitivity grid-points (coherence of masked areas) • dependency from radiation rather than vegetation • evaluation of the overall correction of the OI A comparison with OI (Gain Matrix and OI coefficients) is useful to point out some properties of the variational approach OI DynOI Veg. cover (%) Radiation (W/m2)

  31. Strengths: + suitable in most area in the world, quite cheap analysis + work for soil moiture and soil temperature + take into account model errors (surface model, physiographic database) to provide suitable soil moisture for fitting 2m observations (if no model error sensible and latent heat fluxes are correct). Limitations: + OI coefficients are climatological ones (empirical adjustment to climatological conditions, should be recomputed when model changes)  dynamical OI or variational method + Instantaneous analysis (no assimilation of asynoptic observations) + Not suitable to analyse fast superficial soil moisture evolution  use of precipitation observations (raingauges, radars) Optimal interpolation with 2m obs (pros & cons)

  32. Others analyses: Sea surface temperature Sea ice Snow ...

  33. SST and Sea Ice cover analysis Optimal interpolation assimilating buoys and ships (~1300 obs by rXX) Relaxation towards SST NESDIS analysis 0.5°*0.5° (~5 days time scale) Use SSMI observations to determine Sea Ice (once a day). Temporal consistency in sea ice cover analysis. No lake temperature analysis Snow correction Snow analysis developed in CANARI, but never operational Research study to use either IFS or NESDIS snow cover analysis Snow melting in case of warm T2m observations Frozen soil correction Melting of frozen soil in case of warm T2m observations

  34. Assimilation of SST obs at higher resolution in CANARI Analyse Nesdis à 1/12° Utilisée opérationnellement par ARPEGE Analyse OISST produite par la NOAA Mise en place récente Analyse OSTIA Utilisation des données micro-ondes

  35. 1D simulation over Sodankyla (Finland) T2m juin 05 Eau du sol gelée OBS OPER OLD OPER OLD

  36. Perspectives • Surface analyses implementation in ALADIN (done at CHMI) and under development for AROME • Soil moisture / Soil temperture: - Assimilation of satelite obs (ASCAT, SMOS, SEVIRI, ...) - How to use analysed atmospheric forcings (radiation, precip) - New algorithms (2D-Var, Dyn-OI, ...) - Better consistency with upperair analysis • SST analysis at higher spatial and temporal resolution, Sea Ice Cover analysis, Develop a lake temperature analysis • Need a snow cover analysis (HIRLAM, NESDIS, SAF-Land, ...) • Development of analyses for LAI, VEG, albedo, ...

More Related