The Normalized Difference Vegetation Index (NDVI) and the Enhanced Vegetation Index (EVI) have gained considerable attention in ecological research and management as proxies for landscape-scale vegetation quantity and quality. In the Greater Yellowstone Ecosystem (GYE), these indices are especially important for mapping spatiotemporal variation in the forage available to migratory elk (Cervus elaphus). Here, we examined how the accuracy of using MODIS-derived NDVI and EVI as proxies for forage biomass and quality differed across elevation-related phenology and land use gradients, determined if polynomial NDVI/EVI, site, and season effects improved these models, and then mapped spatiotemporal variation in the abundance of high quality forage available to elk across the Upper Yellowstone River Basin (UYRB) of the GYE. Models with a polynomial NDVI effect explained 19%–55% more variation in biomass than the linear NDVI and EVI models. Models with linear season effect explained 14%–20% more variation in chlorophyll, 37%–69% more variation in crude protein, and 26%–50% more variation in in vitro dry matter digestibility (IVDMD) than the linear NDVI and EVI models. Linear NDVI models explained more variation in biomass and quality across the UYRB than the linear EVI models. The accuracy of these models was lowest in grasslands with late onset of growth, in irrigated agriculture, and after the peak in biomass. Forage biomass and quality varied across the elevation-related phenology and land use gradients in the UYRB throughout the season. At their seasonal peak, the abundance of high quality forage for elk was 50% greater in grasslands with late onset of growth and 200% greater in irrigated agriculture than in all other grasslands, suggesting that these grasslands play an especially important role in the movement and fitness of migratory elk. These results provide novel information on the utility of NDVI and EVI for mapping spatiotemporal patterns of forage biomass and quality.

There is a pressing need to map changes in forest structure from the earliest time period possible given forest management policies and accelerated disturbances from climate change. The availability of Landsat data from over four decades helps researchers study an ecologically meaningful length of time. Forest structure is most often mapped utilizing lidar data, however these data are prohibitively expensive and cover a narrow temporal window relative to the Landsat archive. Here we describe a technique to use the entire length of the Landsat archive from Multispectral Scanner to Operational Land Imager (M2O) to produce three novel outcomes: (1) we used the M2O dataset and standard change vector analysis methods to classify annual forest structure in northwestern Montana from 1972 to 2015, (2) we improved the accuracy of each yearly forest structure classification by applying temporal continuity rules to the whole time series, with final accuracies ranging from 97% to 68% respectively for two and six-category classifications, and (3) we demonstrated the importance of pre-1984 Landsat data for long-term change studies. As the Landsat program continues to acquire Earth imagery into the foreseeable future, time series analyses that aid in classifying forest structure accurately will be key to the success of any land management changes in the future.

Moderate resolution satellite imagery traditionally has been thought to be inadequate for mapping vegetation at the species level. This has made comprehensive mapping of regional distributions of sensitive species, such as whitebark pine, either impractical or extremely time consuming. We sought to determine whether using a combination of moderate resolution satellite imagery (Landsat Enhanced Thematic Mapper Plus), extensive stand data collected by land management agencies for other purposes, and modern statistical classification techniques (boosted classification trees) could result in successful mapping of whitebark pine. Overall classification accuracies exceeded 90%, with similar individual class accuracies. Accuracies on a localized basis varied based on elevation. Accuracies also varied among administrative units, although we were not able to determine whether these differences related to inherent spatial variations or differences in the quality of available reference data.