Dynamic Habitat Indices

Composite RGB 3-1-2
Band 1 - cumulative
Band 2 - minimum
Band 3 - seasonality
Show political borders


All DHIs datasets are single composite RGB images stored in GeoTIFF format at 1 kilometer spatial resolution. Map projection is WGS84 (EPSG:4326).
NDVI16 (Normalized Difference Vegetation Index)
EVI16 (Enhanced Vegetation Index)
FPAR8 (Fraction absorbed Photosynthetically Active Radiation)
LAI8 (Leaf Area Index)
GPP8 (Gross Primary Productivity)
All datasets are freely available, no restrictions are placed on their use, and they can be shared and redistributed. For questions, please contact Volker Radeloff (radeloff@wisc.edu), and we love to learn about applications of the DHIs, and appreciate notifications when there are problems with the datasets. When referring to the DHIs please cite: Hobi, M.L., Dubinin, M., Graham, C.H., Coops, N.C., Clayton, M.K., Pidgeon, A.M., & Radeloff, V.C. (2017). A comparison of Dynamic Habitat Indices derived from different MODIS products as predictors of avian species richness. Remote Sensing of Environment, 195, 142-152.
The DHIs are designed for biodiversity assessments and to describe habitats of different species. Three individual indices comprise the DHIs (Fig. 1):
DHI cum – cumulative DHI, i.e., the area under the phenological curve of a year
DHI min – minimum DHI, i.e., the minimum value of the phenological curve of a year
DHI var – seasonality DHI, i.e., the coefficient of variation of the phenological curve of a year

Fig. 1: Calculation of the three DHIs using productivity (p) at different time periods (t) over the course of a year. Cumulative DHI (DHI Cum) where productivity values are summed up for all the time periods over a year, minimum DHI (DHIMin) where the minimum productivity value within a year is extracted and variation DHI (DHI Var) indicating the seasonality of the productivity by calculating the coefficient of variation using the standard deviation (σ) and themean (μ) over a year. Example of three median phenology curves calculated from 12 years of MODIS FPAR data. The number of time periods within a year is 46 for the 8-day products (FPAR, LAI and GPP) and 23 for the 16-day products (NDVI and EVI)


Downloading and mosaicking of data
We downloaded the yearly MODIS HDF files, and converted them to GeoTIFFs. Dependent on the temporal resolution of the data product this resulted in 46 data sets (one for every 8 days for FPAR, LAI and GPP) or 23 data sets (one for every 16 days in the case of NDVI and EVI) for a given year. For each of these time steps in each year, we mosaicked all tiles together to get a global coverage per time step.
Quality assessment rules
To reduce noise due to clouds or haze, we extracted the associated quality assessment (QA) metadata to identify those pixels that had to be excluded from the analysis. Two different rules had to be applied due to the different quality assessment data provided for the FPAR, LAI and GPP data on one hand, and the NDVI and EVI data on the other hand. For FPAR/LAI/GPP, we applied a threshold for good pixels with a QA < 83 and for the NDVI/EVI data we only used pixels classified as ‘land’ or ‘ocean coastlines and lake shorelines’ (Explicit QA rule: QA<5’411 or (QA>=18’433 & QA<=21’798) or (QA>=34’817 & QA<=38’378) or (QA>=51’201 & QA<=54’574)).
Handling of fill values
For the FPAR and LAI data, we applied an additional rule for the fill values, which we set to zero or no data dependent on their land cover type. Fill values assigned to perennial snow and ice areas (fill value 252) as well as rock, tundra or desert (fill value 253) were set to zero, because these areas are ecologically meaningful in the calculation of the DHI, and vegetation cover is very low. In contrast, we assigned fill values reflecting urban/built-up and populated areas (fill value 250), permanent wetlands/inundated marshlands (fill value 251), as well as perennial salt and inland fresh water (fill value 254) to no data. This mask of no data areas was then also applied to the NDVI, EVI and GPP data sets to make the five DHIs data sets comparable.
Missing data in northern latitudes
In all the data sets we corrected for missing data at the start and the end of the season in the northern latitudes. Because the baseline for these areas is unknown we set these time steps at the beginning and end of the season to zero, if there were records of vegetation productivity during mid-season.
Composite DHI product
The DHIs, especially minimum productivity, are relatively sensitive to occasional changes in data values that are typical for MODIS data, as well as climate events, such as droughts that can affect data for a given year. For some biodiversity analyses, such as analyses of species richness based on range maps, this interannual variation is not relevant. This is why we provide a single composite phenology curve from all MODIS data from 2003 to 2014 in addition to our DHIs for single years. The composite phenology curve represents the median value for each of the 12 observations that were available for each of the 23 or 46 time steps dependent on the MODIS input data. In addition, we only calculated the median if there were at least three valid pixels among the 12 possible years. We then used this composite phenology curve to calculate the DHIs.
Smoothed yearly DHI product
To eliminate noise from the raw data representing a yearly phenology curve, we applied a two-step filtering procedure. The steps included an iterative median and a Savitzky-Golay filter (see Fig. 2). The iterative median filter is based on Chen et al. (2004) and allows to reconstruct time-series. The main feature of the iterative median filter is to eliminate noise that reflects depressed raw values. The main difference of our procedure to the one of Chen et al. (2004) is that we use a median filter in the first step followed by the Savitzky-Golay.

Fig. 2: Two examples of the two filters within the smoothing process based on MODIS FPAR data. The first filter is an iterative median filter and second one is a Savitzky-Golay filter which is applied after the first filter.

Code for data processing
Preparation of MODIS raw data (Python, GDAL)
Smoothing of yearly phenology curves (Python, GRASS)
Unpacking of MODIS QA (ODS tables)


The DHIs are calculated globally based on vegetation productivity estimates derived from Terra and Aqua MODIS Collection 5 satellite data:

   Terra MOD13A2 – NDVI, EVI (16-day temporal resolution)
   Terra/Aqua MOD15A2 – FPAR, LAI (8-day temporal resolution)
   Terra MOD17A2 – GPP (8-day temporal resolution)