ECOCLIMAP Second Generation is the latest version of ECOCLIMAP, produced at 300m-resolution and following a new philosophy.
First, the land cover map is directly composed of the vegetation and urban types used in ISBA and TEB: each grid point of the map represents a pure type, either inland water bodies, or sea and ocean, or one vegetation type or one urban type. In other words, the notion of "cover" or ecosystem of homogeneous land cover type containing several fractions of vegetation types is abandoned.
This novelty leads to simplify the interpretation of the map inside the SURFEX code.
Secondly, the vegetation primary parameters (LAI, ground depths, height of trees, visible and near infrared soil and vegetation albedos), that are defined and averaged by cover in classic ECOCLIMAP, are now given to the PGD step throught the namelist NAM_DATA_ISBA. Consequently, they can be uniform values or input maps coming from satellite data, for example.
Some default maps or values of these parameters are provided with the current version of ECOCLIMAP-SG. However, each user can use his own maps of parameters according to his specific case, he can also add new maps for parameters that are not externalized by the namelist in the default version.
To reduce the amount of memory needed to store these maps, and also to speed up their reading in SURFEX PGD, they were manually compressed following a made in-house method.
Legend of the land cover map and technical documentation¶
The land cover types represented in ECOCLIMAP-SG land cover map are:
1. sea and oceans
4. bare land
5. bare rock
6. permanent snow
7. boreal broadleaf deciduous
8. temperate broadleaf deciduous
9. tropical broadleaf deciduous
10. temperate broadleaf evergreen
11. tropical broadleaf evergreen
12. boreal needleleaf evergreen
13. temperate needleleaf evergreen
14. boreal needleleaf deciduous
16. boreal grassland
17. temperate grassland
18. tropical grassland
19. winter C3 crops
20. summer C3 crops
21. C4 crops
22. flooded trees
23. flooded grassland
24. LCZ1: compact high-rise
25. LCZ2: compact midrise
26. LCZ3: compact low-rise
27. LCZ4: open high-rise
28. LCZ5: open midrise
29: LCZ6: open low-rise
30: LCZ7: lightweight low-rise
31: LCZ8: large low-rise
32: LCZ9: sparsely built
33: LCZ10: heavy industry
A special directory is dedicated to the realization of the maps of parameters: save_scripts_param.tar.gz.
In this directory, the file make_all_fields_param.sh allows to build the maps of parameters without the Kalman filter; the txt file process_kalman.txt explains how to reproduce the creation of the maps with the Kalman filter.
How to download ECOCLIMAP-SG¶
Given the amount of externalized data needed to run ECOCLIMAP-SG, they were dropped off on an external ftp site: ftp.umr-cnrm.fr
To get the data, you need to connect to this ftp site as:
There, you will find several directories:
- COVER: contains ecosg_final_map.dir.gz and ecosg_final_map.hdr.gz, the ECOCLIMAP-SG land cover map input files for SURFEX
- HT: contains new_ht_c.dir.gz and new_ht_c.hdr.gz, the compressed map for the heights of trees.
- LAI/300M: contains the 36 couples of compressed files LAI_MMJJ_c.dir.gz and LAI_MMJJ_c.hdr.gz, one couple for each 10-day period. These LAI data come from the 2014-2016 copernicus satellite LAI data at 300m-resolution (still in "demonstration" stage). Warning: A large fraction of the data in northern latitudes is missing during winter. These missing data have been reconstructed, but are of poor quality.
- LAI/1KM: contains the 36 couples of compressed files LAI_A_MMJJ_c.dir.hz and LAI_MMJJ_c.hdr.gz, one couple for each 10-day period. These LAI data come from the 1999-2016 copernicus satellite LAI data at 1km-resolution, brought to 300m resolution following the method explained in Carrer et al., 2014, that uses a Kalman filter to separate the contributions of each vegetation type present in the pixel to the signal at 1km.
-> If you run a PGD with ECOCLIMAP-SG, you need to choose to use one or the other of these two LAI datasets.
- ALBNIR_SAT & ALBVIS_SAT: contain the 36 couples of compressed files AL-BH-NI_MMJJ_c.dir.gz and AL-BH-NI_MMJJ_c.hdr.gz, one couple for each 10-day period. These albedo data come from the 1999-2016 copernicus satellite albedo data at 1km-resolution.
- ANS, ANV, AVS, AVV: contain the 36 couples of compressed files ANS/ANV/AVS/AVV_DESAG_MMJJ.dir.gz and ANS/ANV/AVS/AVV_DESAG_MMJJ.hdr.gz, one couple for each 10-dau period. These albedo data come from the 1999-2016 copernicus satellite albedo data at 1km resolution, with a separation of the contributions from the soil and the vegetation, following the method explained in Carrer et al., 2014, that uses a Kalman filter to perform this separation.
ANS means "albedo of the soil in the near infrared", ANV "albedo of the vegetation in the near infrared", AVS "albedo of the soil in the visible", AVV "albedo of the soil in the near infrared".
-> If you run a PGD with ECOCLIMAP-SG, you need to choose to use one or the other of these two albedo datasets. If you choose the first one, you will take the same albedo for the soil and for the vegetation.
NB: All the files need to be unzipped (gunzip [file]) to be used as input in SURFEX.
Currently, we use uniform values for the root depth and the ground depth, because we don't have maps for these parameters.
How to install and to use ECOCLIMAP-SG in SURFEX V81¶
If you want to run a PGD with ECOCLIMAP-SG in SURFEX V81, you will first need to activate a new namelist key:
LECOSG = T
Then, you need to link all the input files you want to use in the run directory.
After that, you will need to refer to the correct name in:
YCOVER = 'ecosg_final_map'
YCOVERFILETYPE = 'DIRECT'
And finally, you will need to update your namelist NAM_DATA_ISBA with the names of the files you use for the primary parameters.
Here an example for NAM_DATA_ISBA
In this you will find the values we used for the soil depths during the tests. You also can notice that the type for the compressed maps of parameters is not "DIRECT" but "DIRTYP".
How to visualize an input map (land cover map of parameter map)¶
In the script, you will have to replace LAI_0105_c.dir by the name of the file or a group of file to want to uncompress.
The uncompressed file from a file named [file] will be named [file]_2.