Development and evaluation of a cloud-gap-filled MODIS daily snow-cover product

The utility of the Moderate Resolution Imaging Spectroradiometer (MODIS) snow-cover products is limited by cloud cover which causes gaps in the daily snow-cover map products. We describe a cloud-gap-filled (CGF) daily snow-cover map using a simple algorithm to track cloud persistence, to account for the uncertainty created by the age of the snow observation. Developed from the 0.05° resolution climate-modeling grid daily snow-cover product, MOD10C1, each grid cell of the CGF map provides a cloud-persistence count (CPC) that tells whether the current or a prior day was used to make the snow decision. Percentage of grid cells “observable” is shown to increase dramatically when prior days are considered. The effectiveness of the CGF product is evaluated by conducting a suite of data assimilation experiments using the community Noah land surface model in the NASA Land Information System (LIS) framework. The Noah model forecasts of snow conditions, such as snow-water equivalent (SWE), are updated based on the observations of snow cover which are obtained either from the MOD10C1 standard product or the new CGF product. The assimilation integrations using the CGF maps provide domain-averaged bias improvement of ~11%, whereas such improvement using the standard MOD10C1 maps is ~3%. These improvements suggest that the Noah model underestimates SWE and snow depth fields, and that the assimilation integrations contribute to correcting this systematic error. We conclude that the gap-filling strategy is an effective approach for increasing cloud-free observations of snow cover.     

MODIS Collection 5 global land cover: Algorithm refinements and characterization of new datasets

Information related to land cover is immensely important to global change science. In the past decade, data sources and methodologies for creating global land cover maps from remote sensing have evolved rapidly. Here we describe the datasets and algorithms used to create the Collection 5 MODIS Global Land Cover Type product, which is substantially changed relative to Collection 4. In addition to using updated input data, the algorithm and ancillary datasets used to produce the product have been refined. Most importantly, the Collection 5 product is generated at 500-m spatial resolution, providing a four-fold increase in spatial resolution relative to the previous version. In addition, many components of the classification algorithm have been changed. The training site database has been revised, land surface temperature is now included as an input feature, and ancillary datasets used in post-processing of ensemble decision tree results have been updated. Further, methods used to correct classifier results for bias imposed by training data properties have been refined, techniques used to fuse ancillary data based on spatially varying prior probabilities have been revised, and a variety of methods have been developed to address limitations of the algorithm for the urban, wetland, and deciduous needleleaf classes. Finally, techniques used to stabilize classification results across years have been developed and implemented to reduce year-to-year variation in land cover labels not associated with land cover change. Results from a cross-validation analysis indicate that the overall accuracy of the product is about 75% correctly classified, but that the range in class-specific accuracies is large. Comparison of Collection 5 maps with Collection 4 results show substantial differences arising from increased spatial resolution and changes in the input data and classification algorithm.     

Use of high-resolution satellite imagery in an integrated model to predict the distribution of shade coffee tree hybrid zones

In New Caledonia (21°S, 165°E), shade-grown coffee plantations were abandoned for economic reasons in the middle of the 20th century. Coffee species (Coffea arabica, C. canephora and C. liberica) were introduced from Africa in the late 19th century, they survived in the wild and spontaneously cross-hybridized. Coffee species were originally planted in native forest in association with leguminous trees (mostly introduced species) to improve their growth. Thus the canopy cover over rustic shade coffee plantations is heterogeneous with a majority of large crowns, attributed to leguminous trees. The aim of this study was to identify suitable areas for coffee inter-specific hybridization in New Caledonia using field based environmental parameters and remotely sensed predictors. Due to the complex structure of tropical vegetation, remote sensing imagery needs to be spatially accurate and to have the appropriate bands for monitoring vegetation cover. Quickbird panchromatic (black and white) imagery at 0.6 to 0.7 m spatial resolutions and multispectral imagery at 2.4 m spatial resolution were pansharpened and used for this study. The two most suitable remotely sensed indicators, canopy heterogeneity and tree crown size, were acquired by the sequential use of tree crown detection (neural network), image processing (such as textural analysis) and classification. All models were supervised and trained on learning data determined by human expertise. The final model has two remotely sensed indicators and three physical parameters based on the Digital Elevation Model: elevation, slope and water flow accumulation. Using these five predictive variables as inputs, two modelling methods, a decision tree and a neural network, were implemented. The decision tree, which showed 96.9% accuracy on the test set, revealed the involvement of ecological parameters in the hybridization of Coffea species. We showed that hybrid zones could be characterized by combinations of modalities, underlining the complexity of the environment concerned. For instance, forest heterogeneity and large crown size, steep slopes (> 53.5%) and elevation between 194 and 429 m asl, are favourable factors for Coffea inter-specific hybridization. The application of the neural network on the whole area gave a predictive map that distinguished the most suitable areas by means of a nonlinear continuous indicator. The map provides a confidence level for each area. The most favourable areas were geographically localized, providing a clue for the detection and conservation of favourable areas for Coffea species neo-diversity.     

A scalable approach to mapping annual land cover at 250 m using MODIS time series data: A case study in the Dry Chaco ecoregion of South America

Land use and land cover (LULC) maps from remote sensing are vital for monitoring, understanding and predicting the effects of complex human-nature interactions that span local, regional and global scales. We present a method to map annual LULC at a regional spatial scale with source data and processing techniques that permit scaling to broader spatial and temporal scales, while maintaining a consistent classification scheme and accuracy. Using the Dry Chaco ecoregion in Argentina, Bolivia and Paraguay as a test site, we derived a suite of predictor variables from 2001 to 2007 from the MODIS 250 m vegetation index product (MOD13Q1). These variables included: annual statistics of red, near infrared, and enhanced vegetation index (EVI), phenological metrics derived from EVI time series data, and slope and elevation. For reference data, we visually interpreted percent cover of eight classes at locations with high-resolution QuickBird imagery in Google Earth. An adjustable majority cover threshold was used to assign samples to a dominant class. When compared to field data, we found this imagery to have georeferencing error < 5% the length of a MODIS pixel, while most class interpretation error was related to confusion between agriculture and herbaceous vegetation. We used the Random Forests classifier to identify the best sets of predictor variables and percent cover thresholds for discriminating our LULC classes. The best variable set included all predictor variables and a cover threshold of 80%. This optimal Random Forests was used to map LULC for each year between 2001 and 2007, followed by a per-pixel, 3-year temporal filter to remove disallowed LULC transitions. Our sequence of maps had an overall accuracy of 79.3%, producer accuracy from 51.4% (plantation) to 95.8% (woody vegetation), and user accuracy from 58.9% (herbaceous vegetation) to 100.0% (water). We attributed map class confusion to limited spectral information, sub-pixel spectral mixing, georeferencing error and human error in interpreting reference samples. We used our maps to assess woody vegetation change in the Dry Chaco from 2002 to 2006, which was characterized by rapid deforestation related to soybean and planted pasture expansion. This method can be easily applied to other regions or continents to produce spatially and temporally consistent information on annual LULC.     

Assessing the utility of airborne hyperspectral and LiDAR data for species distribution mapping in the coastal Pacific Northwest, Canada

To effectively manage forested ecosystems an accurate characterization of species distribution is required. In this study we assess the utility of hyperspectral Airborne Imaging Spectrometer for Applications (AISA) imagery and small footprint discrete return Light Detection and Ranging (LiDAR) data for mapping 11 tree species in and around the Gulf Islands National Park Reserve, in coastal South-western Canada. Using hyperspectral imagery yielded producer's and user's accuracies for most species ranging from > 52-95.4 and > 63-87.8%, respectively. For species dominated by definable growth stages, pixel-level fusion of hyperspectral imagery with LiDAR-derived height and volumetric canopy profile data increased both producer's (+ 5.1-11.6%) and user's (+ 8.4-18.8%) accuracies. McNemar's tests confirmed that improvements in overall accuracies associated with the inclusion of LiDAR-derived structural information were statistically significant (p < 0.05). This methodology establishes a specific framework for mapping key species with greater detail and accuracy then is possible using conventional approaches (i.e., aerial photograph interpretation), or either technology on its own. Furthermore, in the study area, acquisition and processing costs were lower than a conventional aerial photograph interpretation campaign, making hyperspectral/LiDAR fusion a viable replacement technology.     

Remote sensing of photosynthetic light-use efficiency across two forested biomes: Spatial scaling

Eddy covariance (EC) measurements have greatly advanced our knowledge of carbon exchange in terrestrial ecosystems. However, appropriate techniques are required to upscale these spatially discrete findings globally. Satellite remote sensing provides unique opportunities in this respect, but remote sensing of the photosynthetic light-use efficiency (ε), one of the key components of Gross Primary Production, is challenging. Some progress has been made in recent years using the photochemical reflectance index, a narrow waveband index centered at 531 and 570 nm. The high sensitivity of this index to various extraneous effects such as canopy structure, and the view observer geometry has so far prevented its use at landscape and global scales. One critical aspect of upscaling PRI is the development of generic algorithms to account for structural differences in vegetation. Building on previous work, this study compares the differences in the PRI: ? relationship between a coastal Douglas-fir forest located on Vancouver Island, British Columbia, and a mature Aspen stand located in central Saskatchewan, Canada. Using continuous, tower-based observations acquired from an automated multi-angular spectro-radiometer (AMSPEC II) installed at each site, we demonstrate that PRI can be used to measure ? throughout the vegetation season at the DF-49 stand (r² = 0.91, p < 0.00) as well as the deciduous site (r² = 0.88, p < 0.00). It is further shown that this PRI signal can be also observed from space at both sites using daily observations from the Moderate Resolution Imaging Spectro-radiometer (MODIS) and a multi-angular implementation of atmospheric correction (MAIAC) (r² = 0.54 DF-49; r² = 0.63 SOA; p < 0.00). By implementing a simple hillshade model derived from airborne light detection and ranging (LiDAR) to approximate canopy shadow fractions (α_s), it is further demonstrated that the differences observed in the relationship between PRI and ε at DF-49 and SOA can be attributed largely to differences in α_s. The findings of this study suggest that algorithms used to separate physiological from extraneous effects in PRI reflectance may be more broadly applicable and portable across these two climatically and structurally different biome types, when the differences in canopy structure are known.     

基于HJ-1高光谱数据的植被覆盖度估测方法研究

植被覆盖度是衡量地表植被状况的一个重要参数,在水文、生态等方面有重要意义,同时,也是影响土壤侵蚀与水土流失的主要因子,是评价土地荒漠化最有效的指标。以环境一号(HJ-1)小卫星上搭载的新型传感器HSI获取的高光谱数据为数据源,通过选择合适的植被指数建立了植被覆盖度反演模型——像元二分模型。然后运用该模型提取了新疆石河子地区的植被覆盖度信息。通过与地面样方数据进行交互比较,对HJ-1/HSI数据反演植被覆盖度的精度进行了评价。研究结果表明,HJ-1/HSI数据能够得到较高精度的植被覆盖度反演结果,在植被动态及全球变化研究领域具有潜在应用价值。     

Tangential cover for thick digital curves

The recognition of digital shapes is a deeply studied problem. The arithmetical framework, initiated by Reveillès [Géométrie discrète, calcul en nombres entiers et algorithmique, Thèse d’Etat, 1991], provides a powerful theoretical basis, as well as many algorithms to deal with digital objects. The tangential cover, first presented in Feschet and Tougne [Optimal time computation of the tangent of a discrete curve: application to the curvature, in: G. Bertrand, M. Couprie, L. Perroton (Eds.), 8th Discrete Geometry for Computer Imagery, Lecture Notes in Computer Science, vol. 1568, Springer, Berlin, 1999, pp. 31-40] and Feschet [Canonical representations of discrete curves, Pattern Anal. Appl. 8(1-2) (2005) 84-94] is a useful tool for representing geometric digital primitives. It computes the set of all maximal segments of a digital curve and permits either to obtain minimal length polygonalization or asymptotic convergence of tangents estimations. Nevertheless, the arithmetical approach does not tolerate the introduction of irregularities, which are however inherent to the acquisition of digital shapes. The present paper is an extension of Faure and Feschet [Tangential cover for thick digital curves, in: D. Coeurjolly, I. Sivignon, L. Tougne, F. Dupont (Eds.), DGCI 2008, Lecture Notes in Computer Science, vol. 4992, Springer, Berlin, 2008, pp. 358-369], in which we propose a new definition for a class of the so-called “thick digital curves” that applies well to a large class of digital object boundaries. We then propose an extension of the tangential cover to thick digital curves and provide an algorithm with an O(nlogn) time complexity, where n denotes the number of points of specific subparts of the thick digital curve. In order to keep up with this low complexity, some critical points must be taken into account. We describe all required implementation details in this paper.     

Bryophyte cover estimation in a boreal black spruce forest using airborne lidar and multispectral sensors

Bryophytes are the dominant ground cover vegetation layer in many boreal forests and in some of these forests the net primary production of bryophytes exceeds the overstory. Therefore it is necessary to quantify their spatial coverage and species composition in boreal forests to improve boreal forest carbon budget estimates. We present results from a small exploratory test using airborne lidar and multispectral remote sensing data to estimate the percentage of ground cover for mosses in a boreal black spruce forest in Manitoba, Canada. Multiple linear regression was used to fit models that combined spectral reflectance data from CASI and indices computed from the SLICER canopy height profile. Three models explained 63-79% of the measured variation of feathermoss cover while three models explained 69-92% of the measured variation of sphagnum cover. Root mean square errors ranged from 3-15% when predicting feathermoss, sphagnum, and total moss ground cover. The results from this case study warrant further testing for a wider range of boreal forest types and geographic regions.     

Integration of MODIS-derived metrics to assess interannual variability in snowpack, lake ice, and NDVI in southwest Alaska

Impacts of global climate change are expected to result in greater variation in the seasonality of snowpack, lake ice, and vegetation dynamics in southwest Alaska. All have wide-reaching physical and biological ecosystem effects in the region. We used Moderate Resolution Imaging Spectroradiometer (MODIS) calibrated radiance, snow cover extent, and vegetation index products for interpreting interannual variation in the duration and extent of snowpack, lake ice, and vegetation dynamics for southwest Alaska. The approach integrates multiple seasonal metrics across large ecological regions.Throughout the observation period (2001-2007), snow cover duration was stable within ecoregions, with variable start and end dates. The start of the lake ice season lagged the snow season by 2 to 3 months. Within a given lake, freeze-up dates varied in timing and duration, while break-up dates were more consistent. Vegetation phenology varied less than snow and ice metrics, with start-of-season dates comparatively consistent across years. The start of growing season and snow melt were related to one another as they are both temperature dependent. Higher than average temperatures during the El Niño winter of 2002-2003 were expressed in anomalous ice and snow season patterns. We are developing a consistent, MODIS-based dataset that will be used to monitor temporal trends of each of these seasonal metrics and to map areas of change for the study area.