1 / 48

PP COLOBOC Status Report

PP COLOBOC Status Report. Jean-Marie Bettems / MeteoSwiss Jürgen Helmert / DWD COSMO GM Moscow, September 8 th , 2010. Goal (COSMO GM 2008). Co nsolidation of Lo wer Bo undary C onditions

presta
Download Presentation

PP COLOBOC Status Report

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. PP COLOBOCStatus Report Jean-Marie Bettems / MeteoSwissJürgen Helmert / DWD COSMO GM Moscow, September 8th, 2010

  2. Goal (COSMO GM 2008) • Consolidation of Lower Boundary Conditions • It is the main goal of this project to incorporate all activities related to the lower boundary conditions which have already reached an advanced state, and to consolidatethesedevelopments into well tested and documented software packages readily usable by the COSMO community.

  3. Actions (COSMO GM 2008) • Consolidate tools of general interest: • externalized TERRA module • software for generation of external parameters • Facilitate verification tasks: • facilitate access to and usage of soil/surface observations • Consolidate and extend external parametersdatabase • Find and validate an optimal configuration of TERRA with its associated external parameters and look-up tables • Revision of snow analysis and snow model • Deployment of urban module developed in Switzerland • Consolidate parameterization of land surface heterogeneity with the tile/mosaic approach

  4. Review – COLOBOC, task 0 Observation sets for SVAT model validation. Documentation Data poolCollect soil / surface / BL observations at selected sites, on behalf of the C-SRNWP Programme. Operational data with time lag from: Lindenberg (D), Payerne (CH), Toulouse (F), San Pietro (I), Sodankylaa (Finland), Cabauw (Netherland). Data of 6 participating sites for 2006-2008 soon available on-line on the COSMO web site. Year 2009 will follow afterwards. Common xls format. Special side meeting on this topic at the EMS conference in Zurich (Wed. Sept. 15th, 14h00 - 15h30). Demo

  5. Review – COLOBOC, task 0 Proposal for extension JM. Bettems/MCH, C.Heret, F.Beyrich/DWD, T.Andreadis/HNMS (0.3 FTE) Extend on-line documentation(e.g. effective data availability of participating sites) Check usability of collected data(Master student at Swiss Federal Institue of Technology) Collect 2009data Include additional observatories(Valday/Russia? Cardington/UK? More?) Make this effort sustainable

  6. Review – COLOBOC, task 0 Open questions Use SRNWP data pool for routine diagnostic of the soil, surface, and boundary layer of all COSMO models … interest from WG5 ?

  7. Review – COLOBOC, task 1 Consolidate externalized TERRA module. Updated package including tests and documentation soon available (integrated within COSMO source code, activated by compiler flag) Code remains fragile when used in a non tested configuration;code limitations will also be documented

  8. Review – COLOBOC, task 1 Task is terminated, but open questions remain Support after end of project, in particular with respect to the introduction of GRIB2 and the new physics interface Provision of COSMO temporal mean values for all driving fields (e.g. T_2M) Code re-write to achieve a production standard

  9. Review – COLOBOC, task 2 Consolidate software for generating external parameters. New code for the aggregation and interpolation of the raw data to the target grid is available (Fortran 90), technical documentation is available Raw data files are converted into NetCDF format,GRIB output will be available soon Reference system at DWD, accessible through a Web interface Demo

  10. Review – COLOBOC, task 3.1 • Document external parameters data set. • Documentation of datasets available on the COSMO web site(http://www.cosmo-model.org/content/model/modules/externalParams/default.htm)

  11. Review – COLOBOC, task 3.2 • Consolidate external parameters data set. • Consolidated raw data sets(GLOBE, GLC2000, DSMW) • 13 new external parameters available for any domain(aerosol, terra module, urban module, lake module) • External parameters for orographic radiationcorrection will soon be available (bug correction required)

  12. Description of External Parameter System GLOBE orography external parameters on target grid DSMW soil data GLC2000 land use

  13. ksh scripts Fortran programs (Fortran 2003/ Fortran 95) user defined grid Icon,GME COSMO ksh scripts Fortran programs global icosahedral grid ksh scripts Fortran programs limited area domain COSMO GME current status external parameter generation new library external parameter generation

  14. Current external parameter fields for the COSMO model, total 15 fields planned extensions for the COSMO model, total 30 fields NDVI, NDVI_MX, NDVI_RATIO, AER_BC, AER_DUST, AER_SO4, AER_SS, T2M_CL, URBAN, RSMIN, EMISS_RAD, FR_LAKE, DEPTH_LK, SLOPE_ASP, SLOPE_ANG, HORIZON, SKYVIEW HH, FR_LAND, SOILTYP, Z0, SSO_STDH, SSO_GAMMA, SSO_THETA, SSO_SIGMA, PLCOV_MN, PLCOV_MX, LAI_MN, LAI_MX, FOREST_D, FOREST_E, ROOTDP HH, FR_LAND, SOILTYP, Z0, SSO_STDH, SSO_GAMMA, SSO_THETA, SSO_SIGMA, PLCOV_MX, LAI_MX, FOREST_D, FOREST_E, ROOTDP Additional external parameters

  15. Example for new ext. parametersFraction lake and lake depth

  16. GME ni 256 based on GLCC GME ni 256 based on GLC2000 COSMO-DE based on GLC2000 COSMO-EU based on GLC2000 Example for uncertainties of ext. parametersEvergreen Forest: GLC2000 vs. GLCC

  17. Review – COLOBOC, task 3.2 Proposal for extension (DWD, MCH, CLM community) (1.0 – 1.5 FTE) Consolidate new MODIS derived solar albedo("background albedo" derived from MODIS data at MPI Hamburg) Consolidate alternative sets of soil types in Europe(Harmonized World Soil Database , European Soil DataBase) Add support for vertically dependent soil information(e.g. depth of water reservoir or inactive layer, texture) Consolidate root depth Alternative vegetation characteristics using MODIS calibrated phenology model(work by R.Stöckli) Allow correct representation of scale separation for z0 / SSO / resolved scales(filter, option for topo smoothing in EXPAR instead of INT2LM) Provide topography at higher resolution than GLOBE when possible(e.g. SRTM or ASTER GDEM)

  18. Review – COLOBOC, task 3.2 Open questions Use of alternative land use data "GLOBCOVER" (higher resolution than GLC2000) Evaluate usage of crop life cycle information (e.g. CLM AgroIBIS, J.L.Roujean/MeteoFrance product)

  19. Review – COLOBOC, task 4 • Revision of TERRA and the associated look-up tables. • Developments of TERRA have been integrated in official COSMO code • Experiments are running to calibrate the land-surface scheme (parametrizations, external parameters, look-up tables). • A consolidated and recommended configuration is expectedtill the end of 2010.

  20. Tested adaptions in int2lm and TERRAEXP 7224 • Aerosol climatology • Emissivity • Vegetation climatology (LAI, PLCOV) • Minimum stomatal resistance • Non-uniform root distribution • Ground water with upward diffusion • Soil moisture dependent heat conductivity

  21. Conclusions of EXP 7224 • Result: less evapotranspiration, dry and warm PBL • Reason : Tested combination of external parameters and TERRA adaptions could lead to increase in PRS/LAI ratio ? • SMA tries to compensate T2M-Bias, wet soil

  22. BATS Plant transpiration (Dickinson, 1984) rstom stomatal resistance rlc res. btw. leaves and canopy air rca res. btw. canopy air and atm. level E LAI rstom

  23. PRS_MIN/LAI_MX ROUTI 1999-2007 PRS_MIN = 150 s/m

  24. PRS_MIN/LAI JuneROUTI 2007- Now

  25. PRS_MIN/LAI JuneEXP Masson (2003) - Ecoclimap

  26. PRS_MIN/LAI JuneEXP Dorman and Sellers (1989)

  27. Latest TERRA developmentsGlobal Verification EXP 7601ECOCLIMAP 01.06.2009-01.07.2009 00 UTC i192f SK: -3.30 SK: -5.08 SK: 2.00 Europe SK: 0.88 SK: -9.72 SK: -6.31 East Sibiria SK: -3.80 SK: -8.23 North America SK: -5.43 SK: 2.73 SK: 1.22 Tropics Africa SK: 3.07 Tropics America SK: 4.49 SK: 3.86 SK: 0.09

  28. Latest TERRA developmentsGlobal Verification EXP 7602D&S89 01.06.2009-01.07.2009 00 UTC i192f SK: -1.88 SK: -4.00 SK: 1.83 Europe SK: -7.35 SK: -5.87 SK: -3.14 East Sibiria SK: -3.55 North America SK: -3.75 SK: -4.81 SK: 2.23 SK: 0.18 Tropics Africa SK: 3.23 Tropics America SK: 4.57 SK: 4.54 SK: 0.02

  29. Latest TERRA developmentsConclusions and Outlook • Result: less evapotranspiration, dry and warm PBL • Reason : Tested combination of external parameters and TERRA adaptions could lead to increase in PRS/LAI ratio ? • Global numerical experiments showed the impact of PRS/LAI ratio on evapotranspiration and screen level temperature • Results from global experiments will be used in COSMO experiments

  30. Review – COLOBOC, task 4 • Proposal for extension • (DWD, MCH, CLM community ) (1.0 FTE) • Revision of TERRA rainfall interception and surface water treatment • Implementation of an orography dependent surface runoff • Detailed comparison of COSMO/TERRA with COSMO/CLM in weather mode (Master student at ETHZ, in group of S.Seneviratne) • Tests in climate mode

  31. Review – COLOBOC, task 4 • Open questions • Introduction of a leaf temperature to replace T_2M(T_2M may be tuned to best match SYNOP observations, which is not necessarily the most adequate value for representing the temperature of the plant environment) • Optimal stretching factor of the soil mesh(a study by A.Will on the characteristics of the numerical solution seems to indicate that the current stretching factor is too large; this could explain the inexistent daily cycle of low cloudiness) • Integrate parameterization of mires (Alla Yurova)

  32. Review – COLOBOC, task 5.1 • Verify and consolidate the new multi-layer snow model • New multi-layers snow model, including • Snow compaction by metamorphism and gravity • Explicit description of radiation effects • Phase transition of liquid water within snow pack • Water percolation • Code is available in latest COSMO release • Tests have shown some improvements of snow depth with the new multi-layers snow model, specially during melting phase and in complex topography. • Some issues remain with the density of the (fresh) snow and with the representation of partial snow cover • Ongoing tests at DWD and MeteoSwiss

  33. Multi-layers snow modelTerra standalone experiments Station at 2450m on the southern slope of the Alps Model at2350m Dh = -100m September 2, 2007 to May 31, 2008

  34. Snow characteristics: SWE and snow depth (1.12.09-28.02.10) north 83 mm 91 mm 95 mm Medvegegorsk : SWE – measurements and 72h forecasts COLOBOC Workshop. Moscow, 6 September 2010

  35. Review – COLOBOC, task 5.1 • Proposal for extension • (E. Machulskaya/DWD, Roshydromet, CLM community) (1.0 – 1.5 FTE) • Correct fresh snow density and feedback effects on soil temperature • Improve representation of snow in forest covered area (albedo …) • Improve partial snow cover representation, in particular by using the tile/mosaic approach • Tests in climate mode

  36. Review – COLOBOC, task 5.2 • Improved snow analysis • MeteoSwiss has extended the DWD snow analysis by introducing a MSG derived snow mask and by tuning the Cressman analysis of snow height observations. Composite fractional snow cover

  37. Review – COLOBOC, task 5.2 • DWD and MeteoSwiss codes have been merged (sfcana 3.0),test are on-going • Altitudinal interpolation has still to be integrated • Temporal stability of analysis has to be improved

  38. Review – COLOBOC, task 5.2 • Proposal for extension • (JM. Bettems/MCH, M. Lange/DWD) • Adapt analysis for the new multi layers snow model

  39. Review – COLOBOC, task 6 Urban module. Code and documentation are ready. Code has already been used(S.Schubert PIK/Potsdam and B.Sändig IfT/Leipzig )

  40. Review – COLOBOC, task 7 Parameterization of land surface heterogeneity Modifications implemented in COSMO v4.8 Documentation has been adapted Still bugs in current code Not all branches/options of the COSMO code have been updated (e.g. multi-layers snow model, recent soil model developments)

  41. Review – COLOBOC, task 7 Proposal for extension (J.Helmert/DWD, A.Mazur/IMGW) (1.0 FTE) Implement tile approach (%) Full support of tile and mosaic in official COSMO code Impact studies (e.g. tile with nature/urban/lake/sea, mosaic for partial snow in complex topo)

  42. Objective Account for non-linear effects of sub-grid inhomegeneities at surface on the exchange of energy and moisture between atmosphere and surface (cf. Ament&Simmer, 2006) tile approach N dominant classes (e.g. water, snow, grass) mosaic approach surface divided in N subgrid cells (Figure taken from Ament&Simmer, 2006)

  43. Objective Coupling of coarse atmosphere and high resolution surface E.g. Latent Heat Flux for one patch: atmospheric variables Grid box average surface variables

  44. Options of TILE & MOSAIC MOSAIC versus TILE approach (Figure taken from Ament, 2006) preference for tile approach

  45. Implications for NWP • Implemention of tile approach requires: • development and implementation of corresponding extensions in external parameter software (i.e. landuse dependend parameters for a no. of dominant classes within each atmospheric grid cell) • code structure to support multiple ‚soil columns‘ within each grid cell (TERRA adaptions in COLOBOC) • physics interface routine or multi-layer soil model, which controls the computation over (flexible) number of classes within each cell and performs necessary aggregation (&disaggregation) • suitable diagnostics (within soil model) to allow proper validation of tile scheme • a computationally efficient and flexible implementation (vectorisation?)

  46. Implications for NWP • AROME (SURFEX) • 4 tiles: nature , town, sea, inland water • Nature: ISBA 3L (Boone et al 1999) • 1L snow scheme (Douville, 1995) • Town • Sea, inland water: constant T_s, Charnock formula UM (Jules) 9 tiles, 5 veg + 4 non-veg Broadleaf and needleleaf trees, temperate and tropical grasses, Shrubs, urban, inland water, bare soil, land ice. IFS (HTESSEL) 6 land-surface tiles High vegetation, low vegetation, interception reservoir, bare ground, snow on ground and low vegetation, Snow under high vegetation

  47. Summary Task 0 [EXTENDED]: Consolidate data exchange action (SRNWP) Task 1 [FINISHED, 09.2010]: Externalized TERRA module Task 2 [FINISHED, 09.2010]: Software for external parameters (incl. GRIB output, web access) Task 3 [FINISHED, 09.2010]: Consolidated raw data sets, new external parameters[EXTENDED]:Further improvements and extensions Task 4 [DELAYED, 12.2010]: Recommended configuration of TERRA and look-up tables[EXTENDED]: Additional improvements of TERRA (interception, surface runoff) Task 5.1 [DELAYED, 11.2010]: Tested multi-layers snow model [EXTENDED] : Fresh snow density, snow in forest, partial snow cover Task 5.2 [DELAYED, 11.2010]: Snow analysis (incl. altitudinal interpolation)[EXTENDED] : Adapt for multi-layers snow model Task 6 [FINISHED, 09.2010]: Urban module available Task 7 [DELAYED, 09.2011]: Parameterization of land surface heterogeneity

  48. Thank you for your attention! COSMO GM Offenbach | postprocessing tooljean-marie.bettems@meteoswiss.ch

More Related