Qingsong Sun

Early Spring Post-Fire Snow Albedo Dynamics in High Latitude Boreal Forests Using Landsat-8 OLI Data

Taking advantage of the improved radiometric resolution of Landsat-8 OLI which, unlike previous Landsat sensors, does not saturate over snow, the progress of fire recovery progress at the landscape scale (< 100m) is examined. High quality Landsat-8 albedo retrievals can now capture the true reflective and layered character of snow cover over a full range of land surface conditions and vegetation densities. This new capability particularly improves the assessment of post-fire vegetation dynamics across low- to high- burn severity gradients in Arctic and boreal regions in the early spring, when the albedos during recovery show the greatest variation. We use 30 m resolution Landsat-8 surface reflectances with concurrent coarser resolution (500m) MODIS high quality full inversion surface Bidirectional Reflectance Distribution Functions (BRDF) products to produce higher resolution values of surface albedo. The high resolution full expression shortwave blue sky albedo product performs well with an overall RMSE of 0.0267 between tower and satellite measures under both snow-free and snow-covered conditions. While the importance of post-fire albedo recovery can be discerned from the MODIS albedo product at regional and global scales, our study addresses the particular importance of early spring post-fire albedo recovery at the landscape scale by considering the significant spatial heterogeneity of burn severity, and the impact of snow on the early spring albedo of various vegetation recovery types. We found that variations in early spring albedo within a single MODIS gridded pixel can be larger than 0.6. Since the frequency and severity of wildfires in Arctic and boreal systems is expected to increase in the coming decades, the dynamics of albedo in response to these rapid surface changes will increasingly impact the energy balance and contribute to other climate processes and physical feedback mechanisms. Surface radiation products derived from Landsat-8 data will thus play an important role in characterizing the carbon cycle and ecosystem processes of high latitude systems.

Monitoring Land Surface Albedo and Vegetation Dynamics Using High Spatial and Temporal Resolution Synthetic Time Series from Landsat and the MODIS BRDF/NBAR/Albedo Product

Seasonal vegetation phenology can significantly alter surface albedo which in turn affects the global energy balance and the albedo warming/cooling feedbacks that impact climate change. To monitor and quantify the surface dynamics of heterogeneous landscapes, high temporal and spatial resolution synthetic time series of albedo and the enhanced vegetation index (EVI) were generated from the 500 m Moderate Resolution Imaging Spectroradiometer (MODIS) operational Collection V006 daily BRDF/NBAR/albedo products and 30 m Landsat 5 albedo and near-nadir reflectance data through the use of the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). The traditional Landsat Albedo (Shuai et al., 2011) makes use of the MODIS BRDF/Albedo products (MCD43) by assigning appropriate BRDFs from coincident MODIS products to each Landsat image to generate a 30 m Landsat albedo product for that acquisition date. The available cloud free Landsat 5 albedos (due to clouds, generated every 16 days at best) were used in conjunction with the daily MODIS albedos to determine the appropriate 30 m albedos for the intervening daily time steps in this study. These enhanced daily 30 m spatial resolution synthetic time series were then used to track albedo and vegetation phenology dynamics over three Ameriflux tower sites (Harvard Forest in 2007, Santa Rita in 2011 and Walker Branch in 2005). These Ameriflux sites were chosen as they are all quite nearby new towers coming on line for the National Ecological Observatory Network (NEON), and thus represent locations which will be served by spatially paired albedo measures in the near future. The availability of data from the NEON towers will greatly expand the sources of tower albedometer data available for evaluation of satellite products. At these three Ameriflux tower sites the synthetic time series of broadband shortwave albedos were evaluated using the tower albedo measurements with a Root Mean Square Error (RMSE) less than 0.013 and a bias within the range of ±0.006. These synthetic time series provide much greater spatial detail than the 500 m gridded MODIS data, especially over more heterogeneous surfaces, which improves the efforts to characterize and monitor the spatial variation across species and communities. The mean of the difference between maximum and minimum synthetic time series of albedo within the MODIS pixels over a subset of satellite data of Harvard Forest (16 km by 14 km) was as high as 0.2 during the snow-covered period and reduced to around 0.1 during the snow-free period. Similarly, we have used STARFM to also couple MODIS Nadir BRDF Adjusted Reflectances (NBAR) values with Landsat 5 reflectances to generate daily synthetic times series of NBAR and thus Enhanced Vegetation Index (NBAR-EVI) at a 30 m resolution. While normally STARFM is used with directional reflectances, the use of the view angle corrected daily MODIS NBAR values will provide more consistent time series. These synthetic times series of EVI are shown to capture seasonal vegetation dynamics with finer spatial and temporal details, especially over heterogeneous land surfaces.

Retrieving understorey dynamics in the Australian tropical savannah from time series decomposition and linear unmixing of MODIS data

ABSTRACT Retrieval from remote sensing of separate temporal dynamics for the understorey layer in tropical savannahs would be beneficial for monitoring fuel loads, biomass for livestock, interrelationships between trees and grasses, and modelling of savannah systems. In this study, we combined unmixing of fractional cover with normalized difference vegetation index (NDVI) and the short wave infrared ratio (SWIR32) with time series decomposition of the NDVI to attempt to fully resolve the dynamics of the herbaceous understorey in the Australian tropical savannah based on the fractions of photosynthetic herbaceous vegetation (FPVH) and non-photosynthetic vegetation (FNPV), from the woody overstorey, represented by the fraction of photosynthetic vegetation in the tree canopy (FPVW). Evaluation of FPVH against field data gave moderate relationships between predicted and observed values (R2 between 0.5 and 0.6); since semivariogram metrics of representativeness indicated that field sites were relatively unrepresentative of variation at the Moderate Resolution Imaging Spectroradiometer MODIS) pixel scale. Both FPVW and FPVH produced strong linear relationships (root mean square error < 0.1 units) with high-resolution Orbview 3 cover fractions classified from tasselled cap transformations. However, FNPVH (non-photosynthetic herbaceous cover fraction) retrievals at southern arid locations produced an evaluation relationship with a greater deviation from the 1:1 line than for northern locations. This suggested that there may be limitations on the NDVI–SWIR32 unmixing approach in more sparsely vegetated savanna. Maps of average annual maximum FPVH, FNPVH, and total herbaceous cover fraction could be used as indicators of savannah productivity and landscape health. However, close examination of the limitations of the NDVI–SWIR32 response may be required for application of this method in other global savannahs.

Relationships between vegetation indices, fractional cover retrievals and the structure and composition of Brazilian Cerrado natural vegetation

ABSTRACT This study explores the use of the relationship between the normalized difference vegetation index (NDVI) and the shortwave infrared ratio (SWIR32) vegetation indices (VI) to retrieve fractional cover over the structurally complex natural vegetation of the Cerrado of Brazil using a time series of imagery from the Moderate Resolution Imaging Spectroradiometer (MODIS). Data from the EO-1 Hyperion sensor with 30 m pixel resolution is used to sample geographic and seasonal variation in NDVI, SWIR32, and the hyperspectral cellulose absorption index (CAI), and to derive end-member values for photosynthetic vegetation (PV), non-photosynthetic vegetation (NPV), and bare soil (BS) from a suite of protected and/or natural vegetation sites across the Cerrado. The end-members derived from relatively pure 30 m pixels are then applied to a 500 m pixel resolution MODIS time series using linear spectral unmixing to retrieve PV, NPV, and BS fractional cover (FPV, FNPV, and FBS). The two-way interaction response of MODIS-equivalent NDVI and SWIR32 was examined for regions of interest (ROI) collected within protected areas and nearby converted lands. The MODIS NDVI, SWIR32 and retrieved FPV, FNPV, and FBS are then compared to detailed cover and structural composition data from field sites, and the influence of the structural and compositional variation on the VIs and cover fractions is explored. The hyperion ROI analysis indicated that the two-way NDVI–SWIR32 response behaved as an effective surrogate for the two-way NDVI–CAI response for the campo limpo/grazed pasture to cerrado sensu stricto woody gradient. The SWIR32 sensitivity to the NPV and BS variation increased as the dry season progressed, but Cerrado savannah exhibited limited dynamic range in the NDVI–CAI and NDVI–SWIR32 two-way responses compared to the entire landscape, which also comprises fallow croplands and forests. Validation analysis of MODIS retrievals with Quickbird-2 images produced an RMSE value of 0.13 for FPV. However, the RMSE values of 0.16 and 0.18 for FBS and FNPV, respectively, were large relative to the seasonal and inter-annual variation. Analysis of site composition and structural data in relation to the MODIS-derived NDVI, SWIR32 and FPV, FNPV, and FBS, indicated that the VI signal and derived cover fractions were influenced by a complex mix of structure and cover but included a strong year-to-year seasonal effect. Therefore, although the MODIS NDVI–SWIR32 response could be used to retrieve cover fractions across all Cerrado land covers including bare cropland, pastures and forests, sensitivity may be limited within the natural Cerrado due to sub-pixel heterogeneity and limited BS and NPV sensitivity.

Dynamics of the relationship between NDVI and SWIR32 vegetation indices in southern Africa: implications for retrieval of fractional cover from MODIS data

ABSTRACT Fractional cover of photosynthetic vegetation (FPV), non-photosynthetic vegetation (FNPV), and bare soil (FBS) has been retrieved for Australian tropical savannah based on linear unmixing of the two-dimensional response envelope of the normalized difference vegetation index (NDVI) and short wave infrared ratio (SWIR)32 vegetation indices (VI) derived from Moderate Resolution Imaging Spectroradiometer (MODIS) reflectance data. The approach assumes that cover fractions are made up of a simple mixture of green leaves, senescent leaves, and bare soil. In this study, we examine retrieval of fractional cover using this approach for a study area in southern Africa with a more complex vegetation structure. Region-specific end-members were defined using Hyperion images from different locations and times of the season. These end-members were applied to a 10-year time series of MODIS-derived NDVI and SWIR32 (from 2002 to 2011) to unmix FPV, FNPV, and FBS. Results of validation with classified high-resolution imagery indicated major bias in estimation of FNPV and FBS, with regression coefficients for predicted versus observed data substantially less than 1.0 and relatively large intercept values. Examination with Hyperion images of the inverse relationship between the MODIS-equivalent SWIR32 index and the Hyperion-derived cellulose absorption index (CAI) to which it nominally approximates revealed: (1) non-compliant positive regression coefficients for certain vegetation types; and (2) shifts in slope and intercept of compliant regression curves related to day of year and geographical location. The results suggest that the NDVI–SWIR32 response cannot be used to approximate the NDVI–CAI response in complex savannah systems like southern Africa that cannot be described as simple mixtures of green leaves, dry herbaceous material high in cellulose, and bare soil. Methods that use a complete set of multispectral channels at higher spatial resolution may be needed for accurate retrieval of fractional cover in Africa.

A Multisensor Approach to Global Retrievals of Land Surface Albedo

Satellite-based retrievals offer the most cost-effective way to comprehensively map the surface albedo of the Earth, a key variable for understanding the dynamics of radiative energy interactions in the atmosphere-surface system. Surface albedo retrievals have commonly been designed separately for each different spaceborne optical imager. Here, we introduce a novel type of processing framework that combines the data from two polar-orbiting optical imager families, the Advanced Very High-Resolution Radiometer (AVHRR) and Moderate Resolution Imaging Spectroradiometer (MODIS). The goal of the paper is to demonstrate that multisensor albedo retrievals can provide a significant reduction in the sampling time required for a robust and comprehensive surface albedo retrieval, without a major degradation in retrieval accuracy, as compared to state-of-the-art single-sensor retrievals. We evaluated the multisensor retrievals against reference in situ albedo measurements and compare them with existing datasets. The results show that global land surface albedo retrievals with a sampling period of 10 days can offer near-complete spatial coverage, with a retrieval bias mostly comparable to existing single sensor datasets, except for bright surfaces (deserts and snow) where the retrieval framework shows degraded performance because of atmospheric correction design compromises. A level difference is found between the single sensor datasets and the demonstrator developed here, pointing towards a need for further work in the atmospheric correction, particularly over bright surfaces, and inter-sensor radiance homogenization. The introduced framework is expandable to include other sensors in the future.

Evaluation of VIIIRS daily BRDF, Albedo, and NBAR product using the MODIS Collection V006 product and in situ measurements

Bidirectional Reflectance Distribution Function (BRDF), Albedo, and Nadir BRDF Adjusted Reflectance (NBAR) products are being produced for the Visible Infrared Imaging Radiometer Suite (VIIRS) onboard the Suomi-National Polar-orbiting Partnership (NPP) satellite in order to extend the MODerate resolution Imaging Spectroradiometer (MODIS) record for research and operational users. The VIIRS product is evaluated by comparison with the MODIS Collection V006 BRDF, Albedo, and NBAR products and in situ albedo collected at spatially representative sites. Preliminary results show that VIIRS can provide comparable BRDF, Albedo, NBAR products as with MODIS. Furthermore, the VIIRS, MODIS and in situ albedos agree well at spatially representative evaluation sites. The accuracy of both products therefore meet the requirements for climate and biosphere models and long term monitoring studies.