The files for soil depth, wilting capacity, and field capacity are collectively known as the 'soil wetness index' files, and contain all the relevant information to calculate the index from the volumetric soil moisture dataset model output. They are provided separately for each of the three soil layers represented in LISFLOOD. They are defined as follows: - Soil depth: The soil depth for each of the three soil layers at each grid point. The value is relative from the top of the land surface to the bottom of each layer respectively.
- Wilting point: The minimal amount of water in the soil that the plant requires not to wilt. If the soil water content decreases to this or any lower point a plant wilts and can no longer recover its turgidity when placed in a saturated atmosphere for 12 hours.
- Field capacity: The amount of soil moisture or water content held in the soil after excess water has drained away and the rate of downward movement has decreased.
More detail on the datasets and how they should be used to derive the soil moisture index is provided here. | Code Block |
|---|
| language | bash |
|---|
| title | Soil depth NetCDF metadata (for soil level 3, and v4.0) |
|---|
| collapse | true |
|---|
| dimensions:
y = 950 ;
x = 1000 ;
variables:
double y(y) ;
y:_FillValue = NaN ;
y:units = "Meter" ;
y:long_name = "Y coordinate of projection" ;
y:standard_name = "projection_Y_coordinate" ;
double x(x) ;
x:_FillValue = NaN ;
x:units = "Meter" ;
x:long_name = "x coordinate of projection" ;
x:standard_name = "projection_x_coordinate" ;
double soilLayer ;
soilLayer:_FillValue = NaN ;
soilLayer:long_name = "original GRIB coordinate for key: level(soilLayer)" ;
soilLayer:units = "1" ;
float sod(y, x) ;
sod:_FillValue = NaNf ;
sod:GRIB_paramId = 260367LL ;
sod:GRIB_dataType = "cf" ;
sod:GRIB_numberOfPoints = 950000LL ;
sod:GRIB_typeOfLevel = "soilLayer" ;
sod:GRIB_stepUnits = 1LL ;
sod:GRIB_stepType = "instant" ;
sod:GRIB_gridType = "lambert_azimuthal_equal_area" ;
sod:GRIB_NV = 0LL ;
sod:GRIB_cfName = "unknown" ;
sod:GRIB_cfVarName = "sod" ;
sod:GRIB_gridDefinitionDescription = "Lambert azimuthal equal area projection " ;
sod:GRIB_missingValue = 9999LL ;
sod:GRIB_name = "Soil depth" ;
sod:GRIB_shortName = "sod" ;
sod:GRIB_totalNumber = 51LL ;
sod:GRIB_units = "m" ;
sod:long_name = "Soil depth" ;
sod:units = "m" ;
sod:standard_name = "unknown" ;
sod:grid_mapping = "lambert_azimuthal_equal_area" ;
sod:coordinates = "soilLayer" ;
int lambert_azimuthal_equal_area ;
lambert_azimuthal_equal_area:semi_major_axis = 6378137. ;
lambert_azimuthal_equal_area:EPSG_code = "EPSG:3035" ;
lambert_azimuthal_equal_area:latitude_of_projection_origin = 52. ;
lambert_azimuthal_equal_area:inverse_flattening = 298.257223563 ;
lambert_azimuthal_equal_area:longitude_of_projection_origin = 10. ;
lambert_azimuthal_equal_area:proj4_params = "+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +units=m +no_defs" ;
lambert_azimuthal_equal_area:false_northing = 3210000. ;
lambert_azimuthal_equal_area:grid_mapping_name = "lambert_azimuthal_equal_area" ;
lambert_azimuthal_equal_area:false_easting = 4321000. ;
lambert_azimuthal_equal_area:coordinates = "soilLayer" ;
byte land_binary_mask(y, x) ;
land_binary_mask:standard_name = "land_binary_mask" ;
land_binary_mask:long_name = "masking_area" ;
land_binary_mask:units = "1" ;
land_binary_mask:grid_mapping = "lambert_azimuthal_equal_area" ;
land_binary_mask:esri_pe_string = "PROJCS[\"ETRS_1989_LAEA\",GEOGCS[\"GCS_ETRS_1989\",DATUM[\"D_ETRS_1989\",SPHEROID[\"GRS_1980\",6378137.0,298.257222101]],PRIMEM[\"Greenwich\",0.0],UNIT[\"Degree\",0.0174532925199433]],PROJECTION[\"Lambert_Azimuthal_Equal_Area\"],PARAMETER[\"false_easting\",4321000.0],PARAMETER[\"false_northing\",3210000.0],PARAMETER[\"central_meridian\",10.0],PARAMETER[\"latitude_of_origin\",52.0],UNIT[\"Meter\",1.0]]" ;
land_binary_mask:coordinates = "soilLayer" ;
// global attributes:
:GRIB_edition = 2LL ;
:GRIB_centre = "ecmf" ;
:GRIB_centreDescription = "European Centre for Medium-Range Weather Forecasts" ;
:GRIB_subCentre = 0LL ;
:Conventions = "CF-1.7" ;
:institution = "European Centre for Medium-Range Weather Forecasts" ;
:history = "2021-11-10T10:31 GRIB to CDM+CF via cfgrib-0.9.9.1/ecCodes-2.23.0 with {\"source\": \"/cache/tmp/afa118b9-8b07-482c-aed7-26976ba41bb7-adaptor.mars.external-1636540185.293234-5781-3-tmp.grib\", \"filter_by_keys\": {}, \"encode_cf\": [\"parameter\", \"time\", \"geography\", \"vertical\"]}" ;
} |
| Code Block |
|---|
| language | bash |
|---|
| title | Wilting point NetCDF metadata (for soil level 3, and v4.0) |
|---|
| collapse | true |
|---|
| dimensions:
x = 1000 ;
y = 950 ;
variables:
double x(x) ;
x:long_name = "x coordinate of projection" ;
x:units = "Meter" ;
x:axis = "X" ;
double y(y) ;
y:long_name = "y coordinate of projection" ;
y:units = "Meter" ;
y:axis = "Y" ;
float wiltingpoint(y, x) ;
wiltingpoint:standard_name = "thetar3" ;
wiltingpoint:long_name = "residual_moisture3" ;
wiltingpoint:units = "mm" ;
wiltingpoint:esri_pe_string = "PROJCS[\"ETRS_1989_LAEA\",GEOGCS[\"GCS_ETRS_1989\",DATUM[\"D_ETRS_1989\",SPHEROID[\"GRS_1980\",6378137.0,298.257222101]],PRIMEM[\"Greenwich\",0.0],UNIT[\"Degree\",0.0174532925199433]],PROJECTION[\"Lambert_Azimuthal_Equal_Area\"],PARAMETER[\"false_easting\",4321000.0],PARAMETER[\"false_northing\",3210000.0],PARAMETER[\"central_meridian\",10.0],PARAMETER[\"latitude_of_origin\",52.0],UNIT[\"Meter\",1.0]]" ;
// global attributes:
:CDI = "Climate Data Interface version 1.8.2 (http://mpimet.mpg.de/cdi)" ;
:Conventions = "CF-1.6" ;
:history = "Tue Dec 11 15:47:51 2018: cdo setname,wiltingpoint test6.nc thmin3.nc\n",
"Tue Dec 11 15:47:49 2018: cdo div test5.nc SDmap3.nc test6.nc\n",
"Tue Dec 11 15:47:49 2018: cdo add test2.nc test4.nc test5.nc\n",
"Tue Dec 11 15:47:49 2018: cdo mul test1.nc fracnoforest.nc test2.nc\n",
"Tue Dec 11 15:47:48 2018: cdo mul thetar3.nc soildepth3_o.nc test1.nc\n",
"Created Tue Feb 03 18:06:55 2015" ;
:source = "Lisflood European maps - pb2015" ;
:conventions = "CF-1.6" ;
:source_software = "Python netCDF3_Classic" ;
:title = "Lisflood maps for European setting Januar 2015" ;
:keywords = "Lisflood, Europe" ;
:institition = "JRC H01" ;
:CDO = "Climate Data Operators version 1.8.2 (http://mpimet.mpg.de/cdo)" ;
} |
| Code Block |
|---|
| language | bash |
|---|
| title | Field capacity NetCDF metadata (for soil level 3, and v4.0) |
|---|
| collapse | true |
|---|
| dimensions:
x = 1000 ;
y = 950 ;
variables:
double x(x) ;
x:long_name = "x coordinate of projection" ;
x:units = "Meter" ;
x:axis = "X" ;
double y(y) ;
y:long_name = "y coordinate of projection" ;
y:units = "Meter" ;
y:axis = "Y" ;
float fieldcapacity(y, x) ;
fieldcapacity:standard_name = "thetas3" ;
fieldcapacity:long_name = "saturated_moisture2" ;
fieldcapacity:units = "mm" ;
fieldcapacity:esri_pe_string = "PROJCS[\"ETRS_1989_LAEA\",GEOGCS[\"GCS_ETRS_1989\",DATUM[\"D_ETRS_1989\",SPHEROID[\"GRS_1980\",6378137.0,298.257222101]],PRIMEM[\"Greenwich\",0.0],UNIT[\"Degree\",0.0174532925199433]],PROJECTION[\"Lambert_Azimuthal_Equal_Area\"],PARAMETER[\"false_easting\",4321000.0],PARAMETER[\"false_northing\",3210000.0],PARAMETER[\"central_meridian\",10.0],PARAMETER[\"latitude_of_origin\",52.0],UNIT[\"Meter\",1.0]]" ;
// global attributes:
:CDI = "Climate Data Interface version 1.8.2 (http://mpimet.mpg.de/cdi)" ;
:Conventions = "CF-1.6" ;
:history = "Tue Dec 11 15:53:48 2018: cdo setname,fieldcapacity thetas3.nc thmax3.nc\n",
"Created Tue Feb 03 18:06:56 2015" ;
:source = "Lisflood European maps - pb2015" ;
:conventions = "CF-1.6" ;
:source_software = "Python netCDF3_Classic" ;
:title = "Lisflood maps for European setting Januar 2015" ;
:keywords = "Lisflood, Europe" ;
:institition = "JRC H01" ;
:CDO = "Climate Data Operators version 1.8.2 (http://mpimet.mpg.de/cdo)" ;
} |
|
