A Spatio-Temporal Enhancement Method for medium resolution LAI (STEM-LAI)

Satellite remote sensing has been used successfully to map leaf area index (LAI) across landscapes, but advances are still needed to exploit multi-scale data streams for producing LAI at both high spatial and temporal resolution. A multi-scale Spatio-Temporal Enhancement Method for medium resolution LAI (STEM-LAI) has been developed to generate 4-day time-series of Landsat-scale LAI from existing medium resolution LAI products. STEM-LAI has been designed to meet the demands of applications requiring frequent and spatially explicit information, such as effectively resolving rapidly evolving vegetation dynamics at sub-field (30 m) scales. In this study, STEM-LAI is applied to Moderate Resolution Imaging Spectroradiometer (MODIS) based LAI data and utilizes a reference-based regression tree approach for producing MODIS-consistent, but Landsat-based, LAI. The Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) is used to interpolate the downscaled LAI between Landsat acquisition dates, providing a high spatial and temporal resolution improvement over existing LAI products. STARFM predicts high resolution LAI by blending MODIS and Landsat based information from a common acquisition date, with MODIS data from a prediction date. To demonstrate its capacity to reproduce fine-scale spatial features observed in actual Landsat LAI, the STEM-LAI approach is tested over an agricultural region in Nebraska. The implementation of a 250 m resolution LAI product, derived from MODIS 1 km data and using a scale consistent approach based on the Normalized Difference Vegetation Index (NDVI), is found to significantly improve accuracies of spatial pattern prediction, with the coefficient of efficiency (E) ranging from 0.77–0.94 compared to 0.01–0.85 when using 1 km LAI inputs alone. Comparisons against an 11-year record of in-situ measured LAI over maize and soybean highlight the utility of STEM-LAI in reproducing observed LAI dynamics (both characterized by r² = 0.86) over a range of plant development stages. Overall, STEM-LAI represents an effective downscaling and temporal enhancement mechanism that predicts in-situ measured LAI better than estimates derived through linear interpolation between Landsat acquisitions. This is particularly true when the in-situ measurement date is greater than 10 days from the nearest Landsat acquisition, with prediction errors reduced by up to 50%. With a streamlined and completely automated processing interface, STEM-LAI represents a flexible tool for LAI disaggregation in space and time that is adaptable to different land cover types, landscape heterogeneities, and cloud cover conditions.     

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.     

Intercomparison and validation of MODIS and GLASS leaf area index (LAI) products over mountain areas: A case study in southwestern China

The validation study of leaf area index (LAI) products over rugged surfaces not only gives additional insights into data quality of LAI products, but deepens understanding of uncertainties regarding land surface process models depended on LAI data over complex terrain. This study evaluated the performance of MODIS and GLASS LAI products using the intercomparison and direct validation methods over southwestern China. The spatio-temporal consistencies, such as the spatial distributions of LAI products and their statistical relationship as a function of topographic indices, time, and vegetation types, respectively, were investigated through intercomparison between MODIS and GLASS products during the period 2011–2013. The accuracies and change ranges of these two products were evaluated against available LAI reference maps over 10 sampling regions which standed for typical vegetation types and topographic gradients in southwestern China.The results show that GLASS LAI exhibits higher percentage of good quality data (i.e. successful retrievals) and smoother temporal profiles than MODIS LAI. The percentage of successful retrievals for MODIS and GLASS is vulnerable to topographic indices, especially to relief amplitude. Besides, the two products do not capture seasonal dynamics of crop, especially in spring over heterogeneously hilly regions. The yearly mean LAI differences between MODIS and GLASS are within ±0.5 for 64.70% of the total retrieval pixels over southwestern China. The spatial distribution of mean differences and temporal profiles of these two products are inclined to be dominated by vegetation types other than topographic indices. The spatial and temporal consistency of these two products is good over most area of grasses/cereal crops; however, it is poor for evergreen broadleaf forest. MODIS presents more reliable change range of LAI than GLASS through comparison with fine resolution reference maps over most of sampling regions. The accuracies of direct validation are obtained for GLASS LAI (r = 0.35, RMSE = 1.72, mean bias = −0.71) and MODIS LAI (r = 0.49, RMSE = 1.75, mean bias = −0.67). GLASS performs similarly to MODIS, but may be marginally inferior to MODIS based on our direct validation results. The validation experience demonstrates the necessity and importance of topographic consideration for LAI estimation over mountain areas. Considerable attention will be paid to the improvements of surface reflectance, retrieval algorithm and land cover types so as to enhance the quality of LAI products in topographically complex terrain.     

Updating Landsat-based forest cover maps with MODIS images using multiscale spectral-spatial-temporal superresolution mapping

With the high deforestation rates of global forest covers during the past decades, there is an ever-increasing need to monitor forest covers at both fine spatial and temporal resolutions. Moderate Resolution Imaging Spectroradiometer (MODIS) and Landsat series images have been used commonly for satellite-derived forest cover mapping. However, the spatial resolution of MODIS images and the temporal resolution of Landsat images are too coarse to observe forest cover at both fine spatial and temporal resolutions. In this paper, a novel multiscale spectral-spatial-temporal superresolution mapping (MSSTSRM) approach is proposed to update Landsat-based forest maps by integrating current MODIS images with the previous forest maps generated from Landsat image. Both the 240 m MODIS bands and 480 m MODIS bands were used as inputs of the spectral energy function of the MSSTSRM model. The principle of maximal spatial dependence was used as the spatial energy function to make the updated forest map spatially smooth. The temporal energy function was based on a multiscale spatial-temporal dependence model, and considers the land cover changes between the previous and current time. The novel MSSTSRM model was able to update Landsat-based forest maps more accurately, in terms of both visual and quantitative evaluation, than traditional pixel-based classification and the latest sub-pixel based super-resolution mapping methods The results demonstrate the great efficiency and potential of MSSTSRM for updating fine temporal resolution Landsat-based forest maps using MODIS images.     

Phenology from Landsat when data is scarce: Using MODIS and Dynamic Time-Warping to combine multi-year Landsat imagery to derive annual phenology curves

Green-leaf phenology describes the development of vegetation throughout a growing season and greatly affects the interaction between climate and the biosphere. Remote sensing is a valuable tool to characterize phenology over large areas but doing at fine- to medium resolution (e.g., with Landsat data) is difficult because of low numbers of cloud-free images in a single year. One way to overcome data availability limitations is to merge multi-year imagery into one time series, but this requires accounting for phenological differences among years. Here we present a new approach that employed a time series of a MODIS vegetation index data to quantify interannual differences in phenology, and Dynamic Time Warping (DTW) to re-align multi-year Landsat images to a common phenology that eliminates year-to-year phenological differences. This allowed us to estimate annual phenology curves from Landsat between 2002 and 2012 from which we extracted key phenological dates in a Monte-Carlo simulation design, including green-up (GU), start-of-season (SoS), maturity (Mat), senescence (Sen), end-of-season (EoS) and dormancy (Dorm). We tested our approach in eight locations across the United States that represented forests of different types and without signs of recent forest disturbance. We compared Landsat-based phenological transition dates to those derived from MODIS and ground-based camera data from the PhenoCam-network. The Landsat and MODIS comparison showed strong agreement. Dates of green-up, start-of-season and maturity were highly correlated (r 0.86-0.95), as were senescence and end-of-season dates (r > 0.85) and dormancy (r > 0.75). Agreement between the Landsat and PhenoCam was generally lower, but correlation coefficients still exceeded 0.8 for all dates. In addition, because of the high data density in the new Landsat time series, the confidence intervals of the estimated keydates were substantially lower than in case of MODIS and PhenoCam. Our study thus suggests that by exploiting multi-year Landsat imagery and calibrating it with MODIS data it is possible to describe green-leaf phenology at much finer spatial resolution than previously possible, highlighting the potential for fine scale phenology maps using the rich Landsat data archive over large areas.     

Evaluating MODIS snow products for modelling snowmelt runoff: Case study of the Rio Grande headwaters

Snow-covered area (SCA) is a key variable in the Snowmelt-Runoff Model (SRM) and in other models for simulating discharge from snowmelt. Landsat Thematic Mapper (TM), Enhanced Thematic Mapper (ETM + ) or Operational Land Imager (OLI) provide remotely sensed data at an appropriate spatial resolution for mapping SCA in small headwater basins, but the temporal resolution of the data is low and may not always provide sufficient cloud-free dates. The coarser spatial resolution Moderate Resolution Imaging Spectroradiometer (MODIS) offers better temporal resolution and in cloudy years, MODIS data offer the best alternative for mapping snow cover when finer spatial resolution data are unavailable. However, MODIS’ coarse spatial resolution (500 m) can obscure fine spatial patterning in snow cover and some MODIS products are not sensitive to end-of-season snow cover. In this study, we aimed to test MODIS snow products for use in simulating snowmelt runoff from smaller headwater basins by a) comparing maps of TM and MODIS-based SCA and b) determining how SRM streamflow simulations are changed by the different estimates of seasonal snow depletion. We compared gridded MODIS snow products (Collection 5 MOD10A1 fractional and binary SCA; SCA derived from Collection 6 MOD10A1 Normalised Difference Snow Index (NDSI) Snow Cover), and the MODIS Snow Covered-Area and Grain size retrieval (MODSCAG) canopy-corrected fractional SCA (SCA_MG), with reference SCA maps (SCA_REF) generated from binary classification of TM imagery. SCA_MG showed strong agreement with SCA_REF; excluding true negatives (where both methods agreed no snow was present) the median percent difference between SCA_REF and SCA_MG ranged between −2.4% and 4.7%. We simulated runoff for each of the four study years using SRM populated with and calibrated for snow depletion curves derived from SCA_REF. We then substituted in each of the MODIS-derived depletion curves. With efficiency coefficients ranging between 0.73 and 0.93, SRM simulation results from the SCA_MG runs yielded the best results of all the MODIS products and only slightly underestimated discharge volume (between 7 and 11% of measured annual discharge). SRM simulations that used SCA derived from Collection 6 NDSI Snow Cover also yielded promising results, with efficiency coefficients ranging between 0.73 and 0.91.In conclusion, we recommend that when simulating snowmelt runoff from small basins (<4000 km²) with SRM, we recommend that users select either canopy-corrected MODSCAG or create their own site-specific products from the Collection 6 MOD10A1 NDSI.     

Methodology for estimating availability of cloud-free image composites: A case study for southern Canada

Image composites are often used for earth surface phenomena studies at regional or national level. The compromise between residual clouds and the length of compositing period is a necessary corollary to the choice of satellite optical data for monitoring earth surface phenomena dynamics. This paper introduced a methodology for estimating availability of cloud-free image composites for optical sensors with various revisiting intervals, using MODIS MOD06 L2 cloud fraction product in the period of 2000–2008. The methodology starts with downscaling of the cloud fraction product to 1 km × 1 km cloud cover binary images. The binary images are then used for the exploration of spatial and temporal characteristics of cloud dynamics, and subsequently for the simulation of cloud-free composite availability with various revisiting intervals of optical sensors. Using Canada's southern provinces as an application case, the study explored several factors important for the design of environmental monitoring system using optical sensors of earth observation, in particular, cloud dynamics and its inter-annual variability, sensors’ revisiting intervals, and cloud-free threshold for targeting composites. While the cloud images used in the analysis are at 1 km × 1 km resolution, our analysis suggests that the simulated availabilities of cloud-free image composites may also provide reasonable estimates for optical sensors with higher than 1 km × 1 km resolution, though the closer to 1 km × 1 km resolution the optical sensor, the more pertinent the application. Also, the methodology can be parameterised to different temporal period and different spatial region, depending on applications.     

Spatially detailed retrievals of spring phenology from single-season high-resolution image time series

Vegetation indices derived from satellite image time series have been extensively used to estimate the timing of phenological events like season onset. Medium spatial resolution (≥250 m) satellite sensors with daily revisit capability are typically employed for this purpose. In recent years, phenology is being retrieved at higher resolution (≤30 m) in response to increasing availability of high-resolution satellite data. To overcome the reduced acquisition frequency of such data, previous attempts involved fusion between high- and medium-resolution data, or combinations of multi-year acquisitions in a single phenological reconstruction. The objectives of this study are to demonstrate that phenological parameters can now be retrieved from single-season high-resolution time series, and to compare these retrievals against those derived from multi-year high-resolution and single-season medium-resolution satellite data. The study focuses on the island of Schiermonnikoog, the Netherlands, which comprises a highly-dynamic saltmarsh, dune vegetation, and agricultural land. Combining NDVI series derived from atmospherically-corrected images from RapidEye (5 m-resolution) and the SPOT5 Take5 experiment (10m-resolution) acquired between March and August 2015, phenological parameters were estimated using a function fitting approach. We then compared results with phenology retrieved from four years of 30 m Landsat 8 OLI data, and single-year 100 m Proba-V and 250 m MODIS temporal composites of the same period. Retrieved phenological parameters from combined RapidEye/SPOT5 displayed spatially consistent results and a large spatial variability, providing complementary information to existing vegetation community maps. Retrievals that combined four years of Landsat observations into a single synthetic year were affected by the inclusion of years with warmer spring temperatures, whereas adjustment of the average phenology to 2015 observations was only feasible for a few pixels due to cloud cover around phenological transition dates. The Proba-V and MODIS phenology retrievals scaled poorly relative to their high-resolution equivalents, indicating that medium-resolution phenology retrievals need to be interpreted with care, particularly in landscapes with fine-scale land cover variability.     

Detection and mapping the spatial distribution of bracken fern weeds using the Landsat 8 OLI new generation sensor

Bracken fern is an invasive plant that presents serious environmental, ecological and economic problems around the world. An understanding of the spatial distribution of bracken fern weeds is therefore essential for providing appropriate management strategies at both local and regional scales. The aim of this study was to assess the utility of the freely available medium resolution Landsat 8 OLI sensor in the detection and mapping of bracken fern at the Cathedral Peak, South Africa. To achieve this objective, the results obtained from Landsat 8 OLI were compared with those derived using the costly, high spatial resolution WorldView-2 imagery. Since previous studies have already successfully mapped bracken fern using high spatial resolution WorldView-2 image, the comparison was done to investigate the magnitude of difference in accuracy between the two sensors in relation to their acquisition costs. To evaluate the performance of Landsat 8 OLI in discriminating bracken fern compared to that of Worldview-2, we tested the utility of (i) spectral bands; (ii) derived vegetation indices as well as (iii) the combination of spectral bands and vegetation indices based on discriminant analysis classification algorithm. After resampling the training and testing data and reclassifying several times (n = 100) based on the combined data sets, the overall accuracies for both Landsat 8 and WorldView-2 were tested for significant differences based on Mann-Whitney U test. The results showed that the integration of the spectral bands and derived vegetation indices yielded the best overall classification accuracy (80.08% and 87.80% for Landsat 8 OLI and WorldView-2 respectively). Additionally, the use of derived vegetation indices as a standalone data set produced the weakest overall accuracy results of 62.14% and 82.11% for both the Landsat 8 OLI and WorldView-2 images. There were significant differences {U (100) = 569.5, z = −10.8242, p < 0.01} between the classification accuracies derived based on Landsat OLI 8 and those derived using WorldView-2 sensor. Although there were significant differences between Landsat and WorldView-2 accuracies, the magnitude of variation (9%) between the two sensors was within an acceptable range. Therefore, the findings of this study demonstrated that the recently launched Landsat 8 OLI multispectral sensor provides valuable information that could aid in the long term continuous monitoring and formulation of effective bracken fern management with acceptable accuracies that are comparable to those obtained from the high resolution WorldView-2 commercial sensor.     

Evaluation of seasonal water body extents in Central Asia over the past 27 years derived from medium-resolution remote sensing data

In this study medium resolution remote sensing data of the AVHRR and MODIS sensors were used for derivation of inland water bodies extents over a period from 1986 till 2012 for the region of Central Asia. Daily near-infrared (NIR) spectra from the AVHRR sensor with 1.1 km spatial resolution and 8-day NIR composites from the MODIS sensor with 250 m spatial resolution for the months April, July and September were used as input data. The methodological approach uses temporal dynamic thresholds for individual data sets, which allows detection of water pixel independent from differing conditions or sensor differences. The individual results are summed up and combined to monthly composites of areal extent of water bodies. The presented water masks for the months April, July, and September were chosen to detect seasonal patterns as well as inter-annual dynamics and show diverse behaviour of static, decreasing, or dynamic water bodies in the study region. The size of the Southern Aral Sea, as the most popular example for an ecologic catastrophe, is decreasing significantly throughout all seasons (R² 0.96 for April; 0.97 for July; 0.96 for September). Same is true for shallow natural lakes in the northern Kazakhstan, exemplary the Tengiz-Korgalzhyn lake system, which have been shrinking in the last two decades due to drier conditions (R² 0.91 for July; 0.90 for September). On the contrary, water reservoirs show high seasonality and are very dynamic within one year in their areal extent with maximum before growing season and minimum after growing season. Furthermore, there are water bodies such as Alakol-Sasykol lake system and natural mountainous lakes which have been stable in their areal extent throughout the entire time period. Validation was performed based on several Landsat images with 30 m resolution and reveals an overall accuracy of 83% for AVHRR and 91% for MODIS monthly water masks. The results should assist for climatological and ecological studies, land and water management, and as input data for different modelling applications.     

Estimating winter wheat biomass by assimilating leaf area index derived from fusion of Landsat-8 and MODIS data

A sufficient number of satellite acquisitions in a growing season are essential for deriving agronomic indicators, such as green leaf area index (GLAI), to be assimilated into crop models for crop productivity estimation. However, for most high resolution orbital optical satellites, it is often difficult to obtain images frequently due to their long revisit cycles and unfavorable weather conditions. Data fusion algorithms, such as the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) and the Enhanced STARFM (ESTARFM), have been developed to generate synthetic data with high spatial and temporal resolution to address this issue. In this study, we evaluated the approach of assimilating GLAI into the Simple Algorithm for Yield Estimation model (SAFY) for winter wheat biomass estimation. GLAI was estimated using the two-band Enhanced Vegetation Index (EVI2) derived from data acquired by the Operational Land Imager (OLI) onboard the Landsat-8 and a fusion dataset generated by blending the Moderate-Resolution Imaging Spectroradiometer (MODIS) data and the OLI data using the STARFM and ESTARFM models. The fusion dataset had the temporal resolution of the MODIS data and the spatial resolution of the OLI data. Key parameters of the SAFY model were optimised through assimilation of the estimated GLAI into the crop model using the Shuffled Complex Evolution-University of Arizona (SCE-UA) algorithm. A good agreement was achieved between the estimated and field measured biomass by assimilating the GLAI derived from the OLI data (GLAI_L) alone (R² = 0.77 and RMSE = 231 g m⁻²). Assimilation of GLAI derived from the fusion dataset (GLAI_F) resulted in a R² of 0.71 and RMSE of 193 g m⁻² while assimilating the combination of GLAI_L and GLAI_F led to further improvements (R² = 0.76 and RMSE = 176 g m⁻²). Our results demonstrated the potential of using the fusion algorithms to improve crop growth monitoring and crop productivity estimation when the number of high resolution remote sensing data acquisitions is limited.     

Satellite-based investigation of flood-affected rice cultivation areas in Chao Phraya River Delta,Thailand

The occurrence of catastrophic floods in Thailand in 2011 caused significant damage to rice agriculture. This study investigated flood-affected rice cultivation areas in the Chao Phraya River Delta (CRD) rice bowl, Thailand using time-series moderate resolution imaging spectroradiometer (MODIS) data. The data were processed for 2008 (normal flood year) and 2011, comprising four main steps: (1) data pre-processing to construct time-series MODIS vegetation indices (VIs), to filter noise from the time-series VIs by the empirical mode decomposition (EMD), and to mask out non-agricultural areas in respect to water-related cropping areas; (2) flood-affected area classification using the unsupervised linear mixture model (ULMM); (3) rice crop classification using the support vector machines (SVM); and (4) accuracy assessment of flood and rice crop mapping results. The comparisons between the flood mapping results and the ground reference data indicated an overall accuracy of 97.9% and Kappa coefficient of 0.62 achieved for 2008, and 95.7% and 0.77 for 2011, respectively. These results were reaffirmed by close agreement (R² > 0.8) between comparisons of the two datasets at the provincial level. The crop mapping results compared with the ground reference data revealed that the overall accuracies and Kappa coefficients obtained for 2008 were 88.5% and 0.82, and for 2011 were 84.1% and 0.76, respectively. A strong correlation was also found between MODIS-derived rice area and rice area statistics at the provincial level (R² > 0.7). Rice crop maps overlaid on the flood-affected area maps showed that approximately 16.8% of the rice cultivation area was affected by floods in 2011 compared to 4.9% in 2008. A majority of the flood-expanded area was observed for the double-cropped rice (10.5%), probably due to flood-induced effects to the autumn–summer and rainy season crops. Information achieved from this study could be useful for agricultural planners to mitigate possible impacts of floods on rice production.     

Semi-automated mapping of burned areas in semi-arid ecosystems using MODIS time-series imagery

Understanding spatial and temporal patterns of burned areas at regional scales, provides a long-term perspective of fire processes and its effects on ecosystems and vegetation recovery patterns, and it is a key factor to design prevention and post-fire restoration plans and strategies. Remote sensing has become the most widely used tool to detect fire affected areas over large tracts of land (e.g., ecosystem, regional and global levels). Standard satellite burned area and active fire products derived from the 500-m Moderate Resolution Imaging Spectroradiometer (MODIS) and the Satellite Pour l’Observation de la Terre (SPOT) are available to this end. However, prior research caution on the use of these global-scale products for regional and sub-regional applications. Consequently, we propose a novel semi-automated algorithm for identification and mapping of burned areas at regional scale. The semi-arid Monte shrublands, a biome covering 240,000 km² in the western part of Argentina, and exposed to seasonal bushfires was selected as the test area. The algorithm uses a set of the normalized burned ratio index products derived from MODIS time series; using a two-phased cycle, it firstly detects potentially burned pixels while keeping a low commission error (false detection of burned areas), and subsequently labels them as seed patches. Region growing image segmentation algorithms are applied to the seed patches in the second-phase, to define the perimeter of fire affected areas while decreasing omission errors (missing real burned areas). Independently-derived Landsat ETM+ burned-area reference data was used for validation purposes. Additionally, the performance of the adaptive algorithm was assessed against standard global fire products derived from MODIS Aqua and Terra satellites, total burned area (MCD45A1), the active fire algorithm (MOD14); and the L3JRC SPOT VEGETATION 1 km GLOBCARBON products. The correlation between the size of burned areas detected by the global fire products and independently-derived Landsat reference data ranged from R² = 0.01–0.28, while our algorithm performed showed a stronger correlation coefficient (R² = 0.96). Our findings confirm prior research calling for caution when using the global fire products locally or regionally.     

Comparative assessment of thematic accuracy of GLC maps for specific applications using existing reference data

Inputs to various applications and models, current global land cover (GLC) maps are based on different data sources and methods. Therefore, comparing GLC maps is challenging. Statistical comparison of GLC maps is further complicated by the lack of a reference dataset that is suitable for validating multiple maps. This study utilizes the existing Globcover-2005 reference dataset to compare thematic accuracies of three GLC maps for the year 2005 (Globcover, LC-CCI and MODIS). We translated and reinterpreted the LCCS (land cover classification system) classifier information of the reference dataset into the different map legends. The three maps were evaluated for a variety of applications, i.e., general circulation models, dynamic global vegetation models, agriculture assessments, carbon estimation and biodiversity assessments, using weighted accuracy assessment. Based on the impact of land cover confusions on the overall weighted accuracy of the GLC maps, we identified map improvement priorities. Overall accuracies were 70.8 ± 1.4%, 71.4 ± 1.3%, and 61.3 ± 1.5% for LC-CCI, MODIS, and Globcover, respectively. Weighted accuracy assessments produced increased overall accuracies (80–93%) since not all class confusion errors are important for specific applications. As a common denominator for all applications, the classes mixed trees, shrubs, grasses, and cropland were identified as improvement priorities. The results demonstrate the necessity of accounting for dissimilarities in the importance of map classification errors for different user application. To determine the fitness of use of GLC maps, accuracy of GLC maps should be assessed per application; there is no single-figure accuracy estimate expressing map fitness for all purposes.     

Probabilistic mapping of flood-induced backscatter changes in SAR time series

The information content of flood extent maps can be increased considerably by including information on the uncertainty of the flood area delineation. This additional information can be of benefit in flood forecasting and monitoring. Furthermore, flood probability maps can be converted to binary maps showing flooded and non-flooded areas by applying a threshold probability value p_F = 0.5. In this study, a probabilistic change detection approach for flood mapping based on synthetic aperture radar (SAR) time series is proposed. For this purpose, conditional probability density functions (PDFs) for land and open water surfaces were estimated from ENVISAT ASAR Wide Swath (WS) time series containing >600 images using a reference mask of permanent water bodies. A pixel-wise harmonic model was used to account for seasonality in backscatter from land areas caused by soil moisture and vegetation dynamics. The approach was evaluated for a large-scale flood event along the River Severn, United Kingdom. The retrieved flood probability maps were compared to a reference flood mask derived from high-resolution aerial imagery by means of reliability diagrams. The obtained performance measures indicate both high reliability and confidence although there was a slight under-estimation of the flood extent, which may in part be attributed to topographically induced radar shadows along the edges of the floodplain. Furthermore, the results highlight the importance of local incidence angle for the separability between flooded and non-flooded areas as specular reflection properties of open water surfaces increase with a more oblique viewing geometry.     

Mapping seagrass coverage and spatial patterns with high spatial resolution IKONOS imagery

Seagrass habitats in subtidal coastal waters provide a variety of ecosystem functions and services and there is an increasing need to acquire information on spatial and temporal dynamics of this resource. Here, we explored the capability of IKONOS (IKO) data of high resolution (4 m) for mapping seagrass cover [submerged aquatic vegetation (%SAV) cover] along the mid-western coast of Florida, USA. We also compared seagrass maps produced with IKO data with that obtained using the Landsat TM sensor with lower resolution (30 m). Both IKO and TM data, collected in October 2009, were preprocessed to calculate water depth invariant bands to normalize the effect of varying depth on bottom spectra recorded by the two satellite sensors and further the textural information was extracted from IKO data. Our results demonstrate that the high resolution IKO sensor produced a higher accuracy than the TM sensor in a three-class % SAV cover classification. Of note is that the OA of %SAV cover mapping at our study area created with IKO data was 5–20% higher than that from other studies published. We also examined the spatial distribution of seagrass over a spatial range of 4–240 m using the Ripley’s K function [L(d)] and IKO data that represented four different grain sizes [4 m (one IKO pixel), 8 m (2 × 2 IKO pixels), 12 m (3 × 3 IKO pixels), and 16 m (4 × 4 IKO pixels)] from moderate-dense seagrass cover along a set of six transects. The Ripley’s K metric repeatedly indicated that seagrass cover representing 4 m × 4 m pixels displayed a dispersed (or slightly dispersed) pattern over distances of <4–8 m, and a random or slightly clustered pattern of cover over 9–240 m. The spatial pattern of seagrass cover created with the three additional grain sizes (i.e., 2 × 24 m IKO pixels, 3 × 34 m IKO pixels, and 4 × 4 m IKO pixels) show a dispersed (or slightly dispersed) pattern across 4–32 m and a random or slightly clustered pattern across 33–240 m. Given the first report on using satellite observations to quantify seagrass spatial patterns at a spatial scale from 4 m to 240 m, our novel analyses of moderate-dense SAV cover utilizing Ripley’s K function illustrate how data obtained from the IKO sensor revealed seagrass spatial information that would be undetected by the TM sensor with a 30 m pixel size. Use of the seagrass classification scheme here, along with data from the IKO sensor with enhanced resolution, offers an opportunity to synoptically record seagrass cover dynamics at both small and large spatial scales.     

Vegetation chlorophyll estimates in the Amazon from multi-angle MODIS observations and canopy reflectance model

As a preparatory study for future hyperspectral missions that can measure canopy chemistry, we introduce a novel approach to investigate whether multi-angle Moderate Resolution Imaging Spectroradiometer (MODIS) data can be used to generate a preliminary database with long-term estimates of chlorophyll. MODIS monthly chlorophyll estimates between 2000 and 2015, derived from a fully coupled canopy reflectance model (ProSAIL), were inspected for consistency with eddy covariance fluxes, tower-based hyperspectral images and chlorophyll measurements. MODIS chlorophyll estimates from the inverse model showed strong seasonal variations across two flux-tower sites in central and eastern Amazon. Marked increases in chlorophyll concentrations were observed during the early dry season. Remotely sensed chlorophyll concentrations were correlated to field measurements (r² = 0.73 and r² = 0.98) but the data deviated from the 1:1 line with root mean square errors (RMSE) ranging from 0.355 μg cm⁻² (Tapajós tower) to 0.470 μg cm⁻² (Manaus tower). The chlorophyll estimates were consistent with flux tower measurements of photosynthetically active radiation (PAR) and net ecosystem productivity (NEP). We also applied ProSAIL to mono-angle hyperspectral observations from a camera installed on a tower to scale modeled chlorophyll pigments to MODIS observations (r² = 0.73). Chlorophyll pigment concentrations (Chl_A+B) were correlated to changes in the amount of young and mature leaf area per month (0.59 ≤  r² ≤  0.64). Increases in MODIS observed Chl_A+B were preceded by increased PAR during the dry season (0.61 ≤ r² ≤  0.62) and followed by changes in net carbon uptake. We conclude that, at these two sites, changes in LAI, coupled with changes in leaf chlorophyll, are comparable with seasonality of plant productivity. Our results allowed the preliminary development of a 15-year time series of chlorophyll estimates over the Amazon to support canopy chemistry studies using future hyperspectral sensors.     

Dry-season vegetation mass and cover fraction from SWIR1.6 and SWIR2.1 band ratio: Ground-radiometer and MODIS data in the Sahel

The potential of the short-wave infrared (SWIR) bands to detect dry-season vegetation mass and cover fraction is investigated with ground radiometry and MODIS data, confronted to vegetation data collected in rangeland and cropland sites in the Sahel (Senegal, Niger, Mali). The ratio of the 1.6 and 2.1 μm bands (called STI) acquired with a ground radiometer proved well suited for grassland mass estimation up to 2500 kg/ha with a linear relation (r² = 0.89). A curvilinear regression is accurate for masses ranging up to 3500 kg/ha. STI proved also well suited to retrieve vegetation cover fraction in crop fields, fallows and rangelands. Such dry-season monitoring, with either ground or satellite data, has important applications for forage, erosion risk and fire risk assessment in semi-arid areas.     

Resolution vs. image quality in pre-tsunami imagery used for tsunami impact models in Aceh,Indonesia

Land cover roughness coefficients (LCRs) have been used in multivariate spatial models to test the mitigation potential of coastal vegetation to reduce impacts of the 2004 tsunami in Aceh, Indonesia. Previously, a Landsat 2002 satellite imagery was employed to derive land cover maps, which were then combined with vegetation characteristics, i.e., stand height, stem diameter and planting density to obtain LCRs. The present study tested LCRs extracted from 2003 and 2004 Landsat (30 m) images as well as a combination of 2003 and 2004 higher spatial resolution SPOT (10 m) imagery, while keeping the previous vegetation characteristics. Transects along the coast were used to extract land cover, whenever availability and visibility allowed. These new LCRs applied in previously developed tsunami impact models on wave outreach, casualties and damages confirmed previous findings regarding distance to the shoreline as a main factor reducing tsunami impacts. Nevertheless, the models using the new LCRs did not perform better than the original one. Particularly casualties models using 2002 LCRs performed better (δAIC > 2) than the more recent Landsat and SPOT counterparts. Cloud cover at image acquisition for Landsat and low area coverage for SPOT images decreased statistical predictive power (fewer observations). Due to the large spatial heterogeneity of tsunami characteristics as well as topographic and land-use features, it was more important to cover a larger area. Nevertheless, if more land cover classes would be referenced and high resolution imagery with low cloud cover would be available, the full benefits of higher spatial resolution imagery used to extract more precise land use roughness coefficients could be exploited.     

Modelling the standing timber volume of Baden-Württemberg—A large-scale approach using a fusion of Landsat,airborne LiDAR and National Forest Inventory data

Remote sensing-based timber volume estimation is key for modelling the regional potential, accessibility and price of lignocellulosic raw material for an emerging bioeconomy. We used a unique wall-to-wall airborne LiDAR dataset and Landsat 7 satellite images in combination with terrestrial inventory data derived from the National Forest Inventory (NFI), and applied generalized additive models (GAM) to estimate spatially explicit timber distribution and volume in forested areas. Since the NFI data showed an underlying structure regarding size and ownership, we additionally constructed a socio-economic predictor to enhance the accuracy of the analysis. Furthermore, we balanced the training dataset with a bootstrap method to achieve unbiased regression weights for interpolating timber volume. Finally, we compared and discussed the model performance of the original approach (r² = 0.56, NRMSE = 9.65%), the approach with balanced training data (r² = 0.69, NRMSE = 12.43%) and the final approach with balanced training data and the additional socio-economic predictor (r² = 0.72, NRMSE = 12.17%). The results demonstrate the usefulness of remote sensing techniques for mapping timber volume for a future lignocellulose-based bioeconomy.