Monitoring post-burn recovery of chaparral vegetation in southern California using multi-temporal satellite data

Monitoring the post-burn recovery condition of chaparral vegetation in southern California is important for managers to determine the appropriate time to conduct controlled burns. Due to the difficulty of monitoring post-fire recovery over large areas and the absence of detailed fire records in many areas, we examined the possibility of using satellite observations to establish the postfire recovery stage of chamise chaparral stands in this region. SPOT XS data collected on three dates between 1986 and 1992 were analysed to determine if temporal changes in a spectral vegetation index tracked the expected post-fire recovery trajectory of the above-ground biomass of chamise chaparral stands of varying post-fire ages. Results of the study indicated that neither the normalized difference vegetation index nor the soil adjusted vegetation index followed the expected post-fire recovery patterns in these stands. These findings are explained by interannual variations in precipitation having a larger than expected effect on the growth of this drought-resistant evergreen community, with changes in green leaf area dominating the temporal variations in the spectral vegetation indices.     

Parallel adjustments in vegetation greenness and ecosystem CO2 exchange in response to drought in a Southern California chaparral ecosystem

Some form of the light use efficiency (LUE) model is used in most models of ecosystem carbon exchange based on remote sensing. The strong relationship between the normalized difference vegetation index (NDVI) and light absorbed by green vegetation make models based on LUE attractive in the remote sensing context. However, estimation of LUE has proven problematic since it varies with vegetation type and environmental conditions. Here we propose that LUE may in fact be correlated with vegetation greenness (measured either as NDVI at constant solar elevation angle, or a red edge chlorophyll index), making separate estimates of LUE unnecessary, at least for some vegetation types. To test this, we installed an automated tram system for measurement of spectral reflectance in the footprint of an eddy covariance flux system in the Southern California chaparral. This allowed us to match the spatial and temporal scales of the reflectance and flux measurements and thus to make direct comparisons over time scales ranging from minutes to years. The 3-year period of this study included both “normal” precipitation years and an extreme drought in 2002. In this sparse chaparral vegetation, diurnal and seasonal changes in solar angle resulted in large variation in NDVI independent of the actual quantity of green vegetation. In fact, one would come to entirely different conclusions about seasonal changes in vegetation greenness depending on whether NDVI at noon or NDVI at constant solar elevation angle were used. Although chaparral vegetation is generally considered “evergreen”, we found that the majority of the shrubs were actually semi-deciduous, leading to large seasonal changes in NDVI at constant solar elevation angle. LUE was correlated with both greenness indices at the seasonal timescale across all years. In contrast, the relationship between LUE and PRI was inconsistent. PRI was well correlated with LUE during the “normal” years but this relationship changed dramatically during the extreme drought. Contrary to expectations, none of the spectral reflectance indices showed consistent relationships with CO₂ flux or LUE over the diurnal time-course, possibly because of confounding effects of sun angle and stand structure on reflectance. These results suggest that greenness indices can be used to directly estimate CO₂ exchange at weekly timescales in this chaparral ecosystem, even in the face of changes in LUE. Greenness indices are unlikely to be as good predictors of CO₂ exchange in dense evergreen vegetation as they were in the sparse, semi-deciduous chaparral. However, since relatively few ecosystems are entirely evergreen at large spatial scales or over long time spans due to disturbance, these relationships need to be examined across a wider range of vegetation types.     

Post-fire recovery of leaf area index in California chaparral: A remote sensing-chronosequence approach

Fire is a major driver of land surface transformation in California Mediterranean-type shrublands (i.e. chaparral). The re-growth of leaves following fire impacts a wide variety of ecosystem processes and information on the post-fire recovery of leaf area index (LAI) is often required in eco-hydrologic modelling studies. A few studies have reported LAI values for chaparral, but none have tracked LAI dynamics over the entire post-fire recovery sequence. In this study we used a chronosequence approach with satellite imagery to determine the post-fire development sequence of LAI for chaparral shrublands in central California. Moreover, we explored how LAI varied with differences in annual antecedent precipitation conditions (APC) and physical site factors. LAI recovery following fire was most rapid during the first 15 years, after which it remained relatively constant with increasing stand age. For a given stand age, LAI varied nonlinearly with annual APC, while spatial variations in LAI were associated with differences in topographic aspect and landscape wetness potential. However, a better understanding of the nature and interaction of these controls on LAI is needed if realistic post-fire LAI trajectories (for historic, present and future periods) for eco-hydrological modelling studies in chaparral catchments are to be developed in the future.     

Characterizing post-fire recovery of fynbos vegetation in the Western Cape Region of South Africa using MODIS data

Post-fire recovery trajectories of five fynbos vegetation stands in the Western Cape Region of South Africa were characterized using moderate-resolution imaging spectroradiometer (MODIS) normalized difference vegetation index (NDVI) 250 m data. Indices of NDVI recovery relative to pre-fire values or values from unburnt control plots indicated full recovery within 7 years and particularly rapid recovery in the first two post-fire years. Intra-stand variability of pixel NDVIs generally increased after fires and also exhibited a rapid recovery to pre-fire conditions. While stand age was the dominant determinant of NDVI recovery, drought interrupted the recovery pathways and this effect was amplified on drier, equator-facing slopes. Post-fire recovery characteristics of fynbos NDVI were found to be similar to those documented for chaparral vegetation in California despite contrasting rainfall and soil nutrient conditions in the two regions.     

The effects of phenology on indirect measures of aboveground biomass in annual grasses

Remote sensing is increasingly being used to quantify vegetation biomass across large areas, often with algorithms based on calibrated relationships between biomass and indices such as the normalized difference vegetation index (NDVI). To improve capacity to evaluate grassland dynamics over time, we examined the influence of phenological changes on NDVI–biomass relationships in annual grasses. Our findings support the use of NDVI throughout early growth and the beginning stages of canopy maturation, but suggest caution for later stages. In contrast, measurements of fractional photosynthetically active radiation (fAPAR) absorbed by the canopy and leaf area index (LAI) served as good season‐long surrogates for canopy biomass. Canopies reached maximum biomass approximately 40 days after maximum greenness, with biomass increasing by approximately 20% during senescence. For multi‐year studies of management impact (i) avoid using seasonal comparisons from dates much after the point of maximum greenness or (ii) consider non‐NDVI‐based approaches.     

Stability,normalization and accuracy of MODIS‐derived estimates of live fuel moisture for southern California chaparral

Time series of spectral vegetation indices (SVIs) derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) and live fuel moisture (LFM) data for chaparral vegetation of southern California were extended to include one of the driest (2004) and one of the wettest (2005) years on record. Independent, spatially varying field‐based estimates of LFM enabled accuracy of MODIS‐derived estimates to be quantified. Pixel‐based scaling of SVI values based on maximum and minimum values of the time series reduced effects of varying vegetation cover and substantially reduced root mean square errors for two of the three SVIs tested.     

Retrieval of leaf area index in different vegetation types using high resolution satellite data

With the successful launch of the IKONOS satellite, very high geometric resolution imagery is within reach of civilian users. In the 1-m spatial resolution images acquired by the IKONOS satellite, details of buildings, individual trees, and vegetation structural variations are detectable. The visibility of such details opens up many new applications, which require the use of geometrical information contained in the images. This paper presents an application in which spectral and textural information is used for mapping the leaf area index (LAI) of different vegetation types. This study includes the estimation of LAI by different spectral vegetation indices (SVIs) combined with image textural information and geostatistical parameters derived from high resolution satellite data. It is shown that the relationships between spectral vegetation indices and biophysical parameters should be developed separately for each vegetation type, and that the combination of the texture indices and vegetation indices results in an improved fit of the regression equation for most vegetation types when compared with one derived from SVIs alone. High within-field spatial variability was found in LAI, suggesting that high resolution mapping of LAI may be relevant to the introduction of precision farming techniques in the agricultural management strategies of the investigated area.     

MODIS-informed greenness responses to daytime land surface temperature fluctuations and wildfire disturbances in the Alaskan Yukon River Basin

Pronounced climate warming and increased wildfire disturbances are known to modify forest composition and control the evolution of the boreal ecosystem over the Yukon River Basin (YRB) in interior Alaska. In this study, we evaluate the post-fire green-up rate using the normalized difference vegetation index (NDVI) derived from 250 m 7 day eMODIS (an alternative and application-ready type of Moderate Resolution Imaging Spectroradiometer (MODIS) data) acquired between 2000 and 2009. Our analyses indicate measureable effects on NDVI values from vegetation type, burn severity, post-fire time, and climatic variables. The NDVI observations from both fire scars and unburned areas across the Alaskan YRB showed a tendency of an earlier start to the growing season (GS); the annual variations in NDVI were significantly correlated to daytime land surface temperature (LST) fluctuations; and the rate of post-fire green-up depended mainly on burn severity and the time of post-fire succession. The higher average NDVI values for the study period in the fire scars than in the unburned areas between 1950 and 2000 suggest that wildfires enhance post-fire greenness due to an increase in post-fire evergreen and deciduous species components.     

Detecting post-fire burn severity and vegetation recovery using multitemporal remote sensing spectral indices and field-collected composite burn index data in a ponderosa pine forest

It is challenging to detect burn severity and vegetation recovery because of the relatively long time period required to capture the ecosystem characteristics. Multitemporal remote sensing data can provide multitemporal observations before, during and after a wildfire, and can improve the change detection accuracy. The goal of this study is to examine the correlations between multitemporal spectral indices and field-observed burn severity, and to provide a practical method to estimate burn severity and vegetation recovery. The study site is the Jasper Fire area in the Black Hills National Forest, South Dakota, that burned during August and September 2000. Six multitemporal Landsat images acquired from 2000 (pre-fire), 2001 (post-fire), 2002, 2003, 2005 and 2007 were used to assess burn severity. The normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), normalized burn ratio (NBR), integrated forest index (IFI) and the differences of these indices between the pre-fire and post-fire years were computed and analysed with 66 field-based composite burn index (CBI) plots collected in 2002. Results showed that differences of NDVI and differences of EVI between the pre-fire year and the first two years post-fire were highly correlated with the CBI scores. The correlations were low beyond the second year post-fire. Differences of NBR had good correlation with CBI scores in all study years. Differences of IFI had low correlation with CBI in the first year post-fire and had good correlation in later years. A CBI map of the burnt area was produced using regression tree models and the multitemporal images. The dynamics of four spectral indices from 2000 to 2007 indicated that both NBR and IFI are valuable for monitoring long-term vegetation recovery. The high burn severity areas had a much slower recovery than the moderate and low burn areas.     

Post-fire vegetation response as a proxy to quantify the magnitude of burn severity in tropical peatland

In recent years, fires in tropical forests in Southeast Asia have become more frequent and widespread, resulting in an increased need to evaluate fire impacts at a landscape scale. We examine whether post-fire vegetation regrowth can be used as a proxy to evaluate burn severity in a peatland landscape in Central Kalimantan, Indonesian Borneo, that has been subject to frequent fires. Several single- and bi-temporal indices as well as spectral fraction endmembers derived from either a post-fire image or a combination of pre- and post-fire images obtained by the Landsat sensor were examined. Spectral data were correlated with vegetation variables obtained from in situ measurements collected 4 years after the last fire. Of the tested spectral data, the bi-temporal and single normalized burn ratio (dNBR and NBR) showed the strongest correlations with the sets of vegetation variables (i.e. total woody aboveground biomass, tree density, and number of trees <10 cm diameter at breast height (DBH)). The results of an analysis of variance (ANOVA) and Tukey's multiple comparison of means test confirmed that NBR, dNBR, and the normalized difference water index could delineate four regrowth classes, thus confirming their utility in separating areas subjected to a single fire from those affected by multiple fires (MFs) as well as for discrimination between fires of differing severity. The results (a) provide evidence of the long-lasting impact that MFs have on forest recovery in this ecosystem and (b) confirm that vegetation response can be used as a proxy to quantify burn severity in locations affected by MFs.     

Comparison of vegetation phenology in the western USA determined from reflected GPS microwave signals and NDVI

The normalized microwave reflection index (NMRI) is a measure of multipath scattering calculated daily from continuously operating GPS sites. GPS satellites transmit L-band microwave signals, and thus NMRI is sensitive to the amount of water in vegetation, not plant greenness or dry biomass. The sensing footprint is approximately 1000 m², although reflections from a distance of hundreds of metres are important at some sites. NMRI exhibits clear seasonal variations that are linked to the changes in vegetation water content that accompany plant growth and senescence. In this paper, NMRI and the normalized difference vegetation index (NDVI) are compared for the period 2008–2012. NMRI data are derived from 184 GPS sites in the western USA. NDVI data are from the 250 m, 16-day pixel containing each GPS station. Amplitude of the annual growth cycle and correlation between NMRI and NDVI are estimated, with and without lags. Phenology metrics are calculated from both indices (i.e. the start of the growing season, timing of peak growth, and season length).

NMRI and NDVI are correlated at most sites, but the degree of correlation varies regionally. Correlation is lowest in California and coastal regions (R = 0.25), where NDVI increases earlier in the spring than NMRI. It is highest for mountain and prairie sites (R = 0.66 and 0.73, respectively). Allowing for a lag between NMRI and NDVI greatly increases the correlation. The lag that yields the greatest correlation is nearly 30 days for the California sites (R = 0.71 with lag), but only 10 days for mountain and prairie sites (R = 0.78 and 0.85 with lag). There are clear differences in phenology metrics extracted from NMRI and NDVI that are consistent with the correlation-lag analysis. Using NMRI, there is a later start to the season, later peak day of the year, and shorter season length. The greatest differences are in California where NDVI start of the season is nearly 60 days earlier than that calculated from NMRI. These data suggest that green-up precedes increases in vegetation water content, with the duration of offset varying regionally. This study is the first to compare GPS-derived microwave reflectance data with NDVI at multiple sites, and highlights both opportunities and limitations offered by NMRI data.     

Seasonal characterization and discrimination of savannah physiognomies in Brazil using hyperspectral metrics from Hyperion/EO-1

The seasonal characterization and discrimination of savannahs in Brazil are still challenging due to the high spatial variability of the vegetation cover and the spectral similarity between some physiognomies. As a preparatory study for future hyperspectral missions that will operate with large swath width and better signal-to-noise ratio than the current orbital sensors, we evaluated six Hyperion images acquired over the Estação Ecológica de Águas Emendadas, a protected area in central Brazil. We studied the seasonal variations in spectral response of the savannah physiognomies and tested their discrimination in the rainy and dry seasons using distinct sets of hyperspectral metrics. Floristic and structural attributes were inventoried in the field. We considered three sets of metrics in the data analysis: the reflectance of 146 Hyperion bands, 22 narrowband vegetation indices (VIs), and 24 absorption band parameters. The VIs were selected to represent vegetation structure, biochemistry, and physiology. The depth, area, width, and asymmetry of the major absorption bands centred at 680 nm (chlorophyll), 980, and 1200 nm (leaf water) and 1700, 2100, and 2300 nm (lignin-cellulose) were calculated on a per-pixel basis using the continuum removal method. Using feature selection and multiple discriminant analysis (MDA), we tested the discriminatory capability of these metrics and of their combined use for vegetation discrimination in the rainy and dry seasons. The results showed that the spectral modifications with seasonality were stronger with the savannah woodland-grassland gradient represented by decreasing tree height, basal area, tree density and biomass and by increasing canopy openness. We observed a reflectance increase in the red, red edge, and shortwave (SWIR) intervals towards the dry season. In the near-infrared, the reflectance differences between the physiognomies were smaller in the dry season than in the rainy season. From the 22 VIs, the visible atmospherically resistant index (VARI), visible green index (VIg), and normalized difference infrared index (NDII) were the most sensitive indices to water stress and vegetation cover, presenting the largest rates of changes between the rainy (March) and dry (August) seasons in shrub and grassland areas. Absorption band parameters associated with the lignin-cellulose spectral features in the SWIR increased towards the dry season with great amounts of non-photosynthetic vegetation (NPV) in the herbaceous stratum. The opposite was observed for the 680 nm chlorophyll absorption band and the 980 and 1200 nm leaf water features. In general, the number of selected metrics necessary for vegetation discrimination was lower in the dry season than in the rainy season. The best MDA-classification accuracy was obtained in the dry season using nine VIs (79.5%). The combination of different hyperspectral metrics increased the classification accuracy to 81.4% in the rainy season and to 84.2% in the dry season. This combination added a gain higher than 10% for the classification of shrub savannah, open woodland savannah and wooded savannah.     

Predicting slow‐drying fire weather index fuel moisture codes with NOAA‐AVHRR images in Canada's northern boreal forests

Fire danger predicted by the Canadian Fire Weather Index, a system based on point‐source weather records, is limited spatially. NOAA‐AVHRR images were used to model two slow‐drying fuel moisture codes, the duff moisture code and the drought code of the fire weather index, in boreal forests of a 250,000 km² portion of northern Alberta and the southern Northwest Territories, Canada. Temporal and spatial factors affecting both codes and spectral variables (normalized difference vegetation index, surface temperature, relative greenness, and the ratio between normalized difference vegetation index and surface temperature) were identified. Models were developed on a yearly and seasonal basis. They were strongest in spring, but had a tendency to saturate. Drought code was best modelled (R ² = 0.34–0.75) in the spring of 1995 when data were categorized spatially by broad forest cover types. These models showed improved spatial resolution by mapping drought code at the pixel level compared to broadly interpolated weather station‐based estimates. Limitations and possible improvements of the study are also discussed.     

Fusing Landsat and MODIS data to retrieve multispectral information from fire-affected areas over tropical savannah environments in the Brazilian Amazon

ABSTRACT

In this study, the combination of surface reflectance products from Terra- Moderate Resolution Imaging Spectroradiometer and Landsat-Enhanced Thematic Mapper Plus sensors are explored through the Flexible Spatiotemporal DAta Fusion (FSDAF) algorithm within the framework of forest fire studies over tropical savannah environments. Thus, 60 fusion-derived images were generated from four spectral bands [red, near-infrared, shortwave infrared (SWIR₁ and SWIR₂)] and six spectral indices [normalized difference vegetation index, normalized difference moisture index, global environment monitoring index, soil-adjusted vegetation index, normalized burn ratio (NBR), and differenced normalized burn ratio (dNBR)] over two selected study sites. For all fusion processes performed, the actual Landsat images for the corresponding dates are available, which supports validation of the blended images. Additionally, integration of blended spectral indices in the immediate post-fire evaluation and the generation of fire severity were analysed. The blended bands presented correlation and Structure Similarity Index Measure (SSIM) values that were consistently higher than 0.819 and root mean square error values of less than 0.027, which confirms good accuracy levels obtained from the model. Similar correlation and SSIM accuracy levels were observed in the blended indices assessment for both study sites, which enables its values to be well-integrated for an analysis of the immediately post-fire date. However, the fire severity mapping from fused images needs to be carefully implemented since the dNBR index is generally less accurate than other blended indices. FSDAF fusion proved to be a useful alternative to retrieving multispectral information from savannah environments affected by fires.     

Understanding burn severity sensing in Arctic tundra: exploring vegetation indices,suboptimal assessment timing and the impact of increasing pixel size

Little is known about how satellite imagery can be used to describe burn severity in tundra landscapes. The Anaktuvuk River Fire (ARF) in 2007 burned over 1000 km² of tundra on the North Slope of Alaska, creating a mosaic of small (1 m²) to large (>100 m²) patches that differed in burn severity. The ARF scar provided us with an ideal landscape to determine if a single-date spectral vegetation index can be used once vegetation recovery began and to independently determine how pixel size influences burn severity assessment. We determine and explore the sensitivity of several commonly used vegetation indices to variation in burn severity across the ARF scar and the influence of pixel size on the assessment and classification of tundra burn severity. We conducted field surveys of spectral reflectance at the peak of the first growing season post-fire (extended assessment period) at 18 field sites that ranged from high to low burn severity. In comparing single-date indices, we found that the two-band enhanced vegetation index (EVI2) was highly correlated with normalized burn ratio (NBR) and better distinguished among three burn severity classes than both the NBR and the normalized difference vegetation index (NDVI). We also show clear evidence that shortwave infrared (SWIR) reflectivity does not vary as a function of burn severity. By comparing a Quickbird scene (2.4 m pixels) to simulated 30 and 250 m pixel scenes, we are able to confirm that while the moderate spatial resolution of the Landsat Thematic Mapper (TM) sensor (30 m) is sufficient for mapping tundra burn severity, the coarser resolution of the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor (250 m) is not well matched to the fine scale of spatial heterogeneity in the ARF burn scar.     

Time-series analysis of NDVI from AVHRR data over the Hindu Kush–Himalayan region for the period 1982–2006

The Hindu Kush–Himalayan (HKH) region with its surrounding mountains in central Asia is a region that has been warming at an alarming rate and is sensitive to climate change due to its heterogeneous terrain and high altitude. In a region where research is limited due to the paucity of field-based biophysical observations, analysis of remotely sensed data such as the normalized difference vegetation index (NDVI) can provide invaluable information on spatio-temporal patterns in linkages among land use, climate and vegetative phenological cycles, and trends in vegetative cover. In this study, NDVI data with 8 km spatial resolution for each 15 day composite period from 1982 to 2006 were analysed using a seasonal trend analysis technique, where the first step determines the annual mean and seasonal NDVI patterns across the HKH region. The second step analyses the non-parametric trends in magnitude and timing of the annual mean and seasonal NDVI cycle. The seasonal vegetation cycles were compared for the first and last ten years of the time series and were also analysed across areas undergoing significant change. Results indicated an overall greening trend in NDVI magnitude in most areas, particularly over open shrubland, grassland and cropland. Trends in the annual seasonal timing of NDVI indicated an earlier green-up for most parts of this region. Results also confirmed deforestation trends observed in a few states in northeastern India and Myanmar (Shan state) within the HKH region.     

Hyperspectral remote sensing for estimating aboveground biomass and for exploring species richness patterns of grassland habitats

Dry grassland sites are amongst the most species-rich habitats of central Europe and it is necessary to design effective management schemes for monitoring of their biomass production. This study explored the potential of hyperspectral remote sensing for mapping aboveground biomass in grassland habitats along a dry-mesic gradient, independent of a specific type or phenological period. Statistical models were developed between biomass samples and spectral reflectance collected with a field spectroradiometer, and it was further investigated to what degree the calibrated biomass models could be scaled to Hyperion data. Furthermore, biomass prediction was used as a surrogate for productivity for grassland habitats and the relationship between biomass and plant species richness was explored. Grassland samples were collected at four time steps during the growing season to capture normally occurring variation due to canopy growth stage and management factors. The relationships were investigated between biomass and (1) existing broad- and narrowband vegetation indices, (2) narrowband normalized difference vegetation index (NDVI) type indices, and (3) multiple linear regression (MLR) with individual spectral bands. Best models were obtained from the MLR and narrowband NDVI-type indices. Spectral regions related to plant water content were identified as the best estimators of biomass. Models calibrated with narrowband NDVI indices were best for up-scaling the field-developed models to the Hyperion scene. Furthermore, promising results were obtained from linking biomass estimations from the Hyperion scene with plant species richness of grassland habitats. Overall, it is concluded that ratio-based NDVI-type indices are less prone to scaling errors and thus offer higher potential for mapping grassland biomass using hyperspectral data from space-borne sensors.     

Normalized sensitivity measures for leaf area index estimation using three-band spectral vegetation indices

A great number of spectral vegetation indices (SVIs) have been developed to estimate key biophysical parameters such as leaf area index (LAI). Considerable interest is often given to the local optimization, performance analysis and sensitivity of each spectral band and SVI for LAI estimation given that several confounding factors are present. In this regard, inclusion of shortwave infrared (SWIR) reflectance in traditionally near-infrared (NIR)-red (R)-based SVIs has played a great role for local optimization and increased sensitivity of SVIs to LAI. This study presents the enhanced and normalized sensitivity functions for evaluating (1) the sensitivity of each spectral band and SVI to LAI and (2) the generic performance analysis of empirical model to estimate LAI based on the SVIs. Several alternatives for three-band (NIR-R-SWIR) SVI modifications have been recommended and proven to be simplistic and unbiased way of local optimization.     

Mapping understory vegetation using phenological characteristics derived from remotely sensed data

Understory vegetation is an important component in forest ecosystems not only because of its contributions to forest structure, function and species composition, but also due to its essential role in supporting wildlife species and ecosystem services. Therefore, understanding the spatio-temporal dynamics of understory vegetation is essential for management and conservation. Nevertheless, detailed information on the distribution of understory vegetation across large spatial extents is usually unavailable, due to the interference of overstory canopy on the remote detection of understory vegetation. While many efforts have been made to overcome this challenge, mapping understory vegetation across large spatial extents is still limited due to a lack of generality of the developed methods and limited availability of required remotely sensed data. In this study, we used understory bamboo in Wolong Nature Reserve, China as a case study to develop and test an effective and practical remote sensing approach for mapping understory vegetation. Using phenology metrics generated from a time series of Moderate Resolution Imaging Spectroradiometer data, we characterized the phenological features of forests with understory bamboo. Using maximum entropy modeling together with these phenology metrics, we successfully mapped the spatial distribution of understory bamboo (kappa: 0.59; AUC: 0.85). In addition, by incorporating elevation information we further mapped the distribution of two individual bamboo species, Bashania faberi and Fargesia robusta (kappa: 0.68 and 0.70; AUC: 0.91 and 0.92, respectively). Due to its generality, flexibility and extensibility, this approach constitutes an improvement to the remote detection of understory vegetation, making it suitable for mapping different understory species in different geographic settings. Both biodiversity conservation and wildlife habitat management may benefit from the detailed information on understory vegetation across large areas through the applications of this approach.