Global, 30-m resolution continuous fields of tree cover: Landsat-based rescaling of MODIS vegetation continuous fields with lidar-based estimates of error

Abstract

We developed a global, 30-m resolution dataset of percent tree cover by rescaling the 250-m MOderate-resolution Imaging Spectroradiometer (MODIS) Vegetation Continuous Fields (VCF) Tree Cover layer using circa- 2000 and 2005 Landsat images, incorporating the MODIS Cropland Layer to improve accuracy in agricultural areas. Resulting Landsat-based estimates maintained consistency with the MODIS VCF in both epochs (RMSE =8.6% in 2000 and 11.9% in 2005), but showed improved accuracy in agricultural areas and increased discrimination of small forest patches. Against lidar measurements, the Landsat-based estimates exhibited accuracy slightly less than that of the MODIS VCF (RMSE=16.8% for MODIS-based vs. 17.4% for Landsat-based estimates), but RMSE of Landsat estimates was 3.3 percentage points lower than that of the MODIS data in an agricultural region. The Landsat data retained the saturation artifact of the MODIS VCF at greater than or equal to 80% tree cover but showed greater potential for removal of errors through calibration to lidar, with post-calibration RMSE of 9.4% compared to 13.5% in MODIS estimates. Provided for free download at the Global Land Cover Facility (GLCF) website (www.landcover.org), the 30-m resolution GLCF tree cover dataset is the highest-resolution multi-temporal depiction of Earth's tree cover available to the Earth science community.     

大尺度土地覆盖数据集在中国及周边区域的精度评价

大尺度土地覆盖数据是全球陆地表层过程研究、生态系统评估、环境建模等科学研究的重要基础,研究现有数据集的特点对数据使用者及生产新的数据集都具有指导意义。本研究以中国及周边区域为研究区,根据不同分类体系对地物的定义,研究不同分类体系中对应地物的相关系数,并将所有分类体系转换为IGBP分类体系;然后,从定性和定量两方面分析现有5种土地覆盖数据集(IGBP DISCover、UMD、GLC2000、MOD12Q1和GlobCover 2005)的空间一致性;并利用Google Earth高分影像选取两期验证样本评价5种土地覆盖数据集的精度。结果表明:同种地物在不同土地覆盖数据集之间的空间分布格局差异较大,且不同土地覆盖数据集之间的总体一致性系数较低;5种土地覆盖数据集中,GLC2000的总体精度和Kappa系数均最高,GlobCover 2005的总体精度和Kappa系数均最低。     

地表覆盖分类数据对区域森林叶面积指数反演的影响

以江西省吉安市为研究区,将5种全球地表覆盖分类数据(包括美国地质调查局(USGS)、马里兰大学(UMD)和波士顿大学(BU)生成的3套数据和欧洲生成的2套数据)以及由TM影像生成的区域地表覆盖分类数据,分别与MODIS1km反射率资料结合,利用基于4尺度几何光学模型的LAI反演方法生成研究区的LAI。在1km和4km两种尺度上将反演的LAI与TM资料生成的LAI进行比较,评价地表覆盖分类数据对LAI反演结果的影响。结果表明,TM和欧洲太空局的GLOBCOVER地表覆盖分类数据用于反演LAI的结果较好,在1km尺度上,反演的LAI与统计模型估算的TMLAI相关的R2分别为0.44和0.40,在4km尺度上的R2分别为0.57和0.54;其次为波士顿大学的MODIS地表覆盖分类数据,据其反演的LAI与TMLAI相关的R2在1km和4km尺度上分别为0.38和0.51;而马里兰大学的UMD和欧洲的GLC2000地表覆盖分类数据会导致反演的LAI存在较大误差,据其反演的LAI与TMLAI之间的一致性较差,在1km和4km两种尺度上平均偏低20%左右;LAI的反演结果对聚集度系数具有强的敏感性。该研究表明,为了提高区域/全球LAI反演精度,需要有高质量的地表覆盖分类数据。     

Production of global land cover data – GLCNMO

Abstract

Global land cover is one of the fundamental contents of Digital Earth. The Global Mapping project coordinated by the International Steering Committee for Global Mapping has produced a 1-km global land cover dataset – Global Land Cover by National Mapping Organizations. It has 20 land cover classes defined using the Land Cover Classification System. Of them, 14 classes were derived using supervised classification. The remaining six were classified independently: urban, tree open, mangrove, wetland, snow/ice, and water. Primary source data of this land cover mapping were eight periods of 16-day composite 7-band 1-km MODIS data of 2003. Training data for supervised classification were collected using Landsat images, MODIS NDVI seasonal change patterns, Google Earth, Virtual Earth, existing regional maps, and expert's comments. The overall accuracy is 76.5% and the overall accuracy with the weight of the mapped area coverage is 81.2%. The data are available from the Global Mapping project website (http://www.iscgm.org/). The MODIS data used, land cover training data, and a list of existing regional maps are also available from the CEReS website. This mapping attempt demonstrates that training/validation data accumulation from different mapping projects must be promoted to support future global land cover mapping.     

Evaluation of global land cover maps for cropland area estimation in the conterminous United States

Global land cover data could provide continuously updated cropland acreage and distribution information, which is essential to a wide range of applications over large geographical regions. Cropland area estimates were evaluated in the conterminous USA from four recent global land cover products: MODIS land cover (MODISLC) at 500-m resolution in 2010, GlobCover at 300-m resolution in 2009, FROM-GLC and FROM-GLC-agg at 30-m resolution based on Landsat imagery circa 2010 against the US Department of Agriculture survey data. Ratio estimators derived from the 30-m resolution Cropland Data Layer were applied to MODIS and GlobCover land cover products, which greatly improved the estimation accuracy of MODISLC by enhancing the correlation and decreasing mean deviation (MDev) and RMSE, but were less effective on GlobCover product. We found that, in the USA, the CDL adjusted MODISLC was more suitable for applications that concern about the aggregated county cropland acreage, while FROM-GLC-agg gave the least deviation from the survey at the state level. Correlation between land cover map estimates and survey estimates is significant, but stronger at the state level than at the county level. In regions where most mismatches happen at the county level, MODIS tends to underestimate, whereas MERIS and Landsat images incline to overestimate. Those uncertainties should be taken into consideration in relevant applications. Excluding interannual and seasonal effects, R² of the FROM-GLC regression model increased from 0.1 to 0.4, and the slope is much closer to one. Our analysis shows that images acquired in growing season are most suitable for Landsat-based cropland mapping in the conterminous USA.     

全球30m地表覆盖遥感制图的总体技术

本文针对全球30米分辨率地表覆盖遥感制图这一世界性难题,提出了以多源影像最优化处理、参考资料服务化整合、覆盖类型精细化提取、产品质量多元化检核为主线的总体研究思路,研发了影像几何与辐射重建、异质异构服务化集成、对象化分层分类、知识化检核处理等主体技术方法;用于制定了相应数据产品规范、生产技术规范,研发了多项生产型软件,用于研制了2000和2010两个基准年的全球）30米地表覆盖数据产品,将空间分辨率提高了1个数量级。     

Primitives as building blocks for constructing land cover maps

Land cover maps play an integral role in environmental management. However, countries and institutes encounter many challenges with producing timely, efficient, and temporally harmonized updates to their land cover maps. To address these issues we present a modular Regional Land Cover Monitoring System (RLCMS) architecture that is easily customized to create land cover products using primitive map layers. Primitive map layers are a suite of biophysical and end member maps, with land cover primitives representing the raw information needed to make decisions in a dichotomous key for land cover classification. We present best practices to create and assemble primitives from optical satellite using computing technologies, decision tree logic and Monte Carlo simulations to integrate their uncertainties. The concept is presented in the context of a regional land cover map based on a shared regional typology with 18 land cover classes agreed on by stakeholders from Cambodia, Laos PDR, Myanmar, Thailand, and Vietnam. We created annual map and uncertainty layers for the period 2000–2017. We found an overall accuracy of 94% when taking uncertainties into account. RLCMS produces consistent time series products using free long term historical Landsat and MODIS data. The customizable architecture can include a variety of sensors and machine learning algorithms to create primitives and the best suited smoothing can be applied on a primitive level. The system is transferable to all regions around the globe because of its use of publicly available global data (Landsat and MODIS) and easily adaptable architecture that allows for the incorporation of a customizable assembly logic to map different land cover typologies based on the user's landscape monitoring objectives     

Comparison of MODIS vegetation continuous field — based forest density maps with IRS-LISS III derived maps

The study compared forest cover maps derived using coarse resolution vegetation continuous fields (MODIS VCF; 500m resolution) with the maps derived from medium resolution (24m; IRS LISS-III) data. The comparison of VCF, per cent tree cover product, for the years 2000 to 2004 with LISS III forest density class maps of 2001 and 2003 was carried out for two sites representing hilly (Uttarakhand) and undulating terrains (Madhya Pradesh). Slicing VCF to corresponding forest crown cover, i.e., 0–10%, 10–40%, 40–70% and >70% produced considerable difference in forest area estimates when compared to original LISS III derived crown cover area. The corresponding value range in VCF for 0–10% of actual forest cover were 0–31% and 0–25% in 2 sites respectively, and the respective limit was consistent at 1–20% when VCF range were sliced with respect to upscaled LISS III at 500m resolution. Similarly, all other class limits were also found through iterative process. These limits were similar, within a site, across five years. Spatial Kappa match between these two data indicated higher match in 40–70% class, and also in undulating site. When compared at same resolution, similar forest area cover estimated with weighted area upscaling gave closest match. The study is useful in knowing the usability and limits of VCF product, and utility of spatial Kappa.     

Land cover mapping applications with MODIS: a literature review

Abstract

Land use/land cover monitoring and mapping is crucial to efficient management of the land and its resources. Since the late 1980s increased attention has been paid to the use of coarse resolution optical data. The Moderate Resolution Imaging Spectroradiometer (MODIS) has features, which make it particularly suitable to earth characterization purposes. MODIS has 10 products dedicated mainly to land cover characterization and provides three kinds of data: angular, spectral and temporal. MODIS data also includes information about the data quality through the ‘Quality Assessment’ product. In this paper, we review how MODIS data are used to map land cover including the preferred MODIS products, the preprocessing and classification approaches, the accuracy assessment, and the results obtained.     

Major forest changes and land cover transitions based on plant functional types derived from the ESA CCI Land Cover product

Land use and land cover change are of prime concern due to their impacts on CO₂ emissions, climate change and ecological services. New global land cover products at 300 m resolution from the European Space Agency (ESA) Climate Change Initiative Land Cover (CCI LC) project for epochs centered around 2000, 2005 and 2010 were analyzed to investigate forest area change and land cover transitions. Plant functional types (PFTs) fractions were derived from these land cover products according to a conversion table. The gross global forest loss between 2000 and 2010 is 172,171 km², accounting for 0.6% of the global forest area in year 2000. The forest changes are mainly distributed in tropical areas such as Brazil and Indonesia. Forest gains were only observed between 2005 and 2010 with a global area of 9844 km², mostly from crops in Southeast Asia and South America. The predominant PFT transition is deforestation from forest to crop, accounting for four-fifths of the total increase of cropland area between 2000 and 2010. The transitions from forest to bare soil, shrub, and grass also contributed strongly to the total areal change in PFTs. Different PFT transition matrices and composition patterns were found in different regions. The highest fractions of forest to bare soil transitions were found in the United States and Canada, reflecting forest management practices. Most of the degradation from grassland and shrubland to bare soil occurred in boreal regions. The areal percentage of forest loss and land cover transitions generally decreased from 2000–2005 to 2005–2010. Different data sources and uncertainty in the conversion factors (converting from original LC classes to PFTs) contribute to the discrepancy in the values of change in absolute forest area.     

Open land cover from OpenStreetMap and remote sensing

OpenStreetMap (OSM) tags were used to produce a global Open Land Cover (OLC) product with fractional data gaps available at osmlanduse.org. Data gaps in the global OLC map were filled for a case study in Heidelberg, Germany using free remote sensing data, which resulted in a land cover (LC) prototype with complete coverage in this area. Sixty tags in the OSM were used to allocate a Corine Land Cover (CLC) level 2 land use classification to 91.8% of the study area, and the remaining gaps were filled with remote sensing data. For this case study, complete are coverage OLC overall accuracy was estimated 87%, which performed better than the CLC product (81% overall accuracy) of 2012. Spatial thematic overlap for the two products was 84%. OLC was in large parts found to be more detailed than CLC, particularly when LC patterns were heterogeneous, and outperformed CLC in the classification of 12 of the 14 classes. Our OLC product represented data created in different periods; 53% of the area was 2011–2016, and 46% of the area was representative of 2016–2017.     

基于直方变差图的土地覆盖数据集多尺度精度评价

针对目前精度评价尺度单一的问题,提出基于直方变差图的多尺度精度评价方法,分别在像元尺度和亚像元尺度进行土地覆盖数据集精度评价。在像元尺度利用驻点作为采样工具直接评价数据集精度;亚像元尺度上,则利用非严格定义的驻点和驻点直方变差图对不同面积和空间结构的优势类进行精度评价。并以浙江北部典型区域为实验区,Landsat TM/ETM+为参考数据,对UMD、IGBP DISCover、MOD12Q1-2001、GLC2000、GlobCover2009等5种大尺度土地覆盖数据集进行多尺度精度评价实验。结果表明,多尺度精度评价方法能够全面地评价土地覆盖数据集的精度,提供更加丰富的多尺度精度信息。像元尺度精度评价可在一定程度上消除由于参考数据与数据集间的空间匹配造成的误差,评价结果更加客观;亚像元尺度精度评价能有效反映亚像元尺度优势地物面积及空间结构与精度的关系。     

Collaborative validation of GlobeLand30:Methodology and practices

30-m Global Land Cover(GLC)data products permit the detection of land cover changes at the scale of most human land activities,and are therefore used as fundamental information for sustainable development,environmental change studies,and many other societal benefit areas.In the past few years,increasing efforts have been devoted to the accuracy assessment of GlobeLand30 and other finer-resolution GLC data products.However,most of them were conducted either within a limited percentage of map sheets selected from a global scale or in some individual countries(areas),and there are still many areas where the uncertainty of 30-m resolution GLC data products remains to be validated and documented.In order to promote a comprehensive and collaborative validation of 30-m GLC data products,the GEO Global Land Cover Community Activity had organized a project from 2015 to 2017,to examine and explore its major problems,including the lack of international agreed validation guidelines and on-line tools for facilitating collaborative validation activities.With the joint effort of experts and users from 30 GEO member countries or participating organizations,a technical specification for 30-m GLC validation was developed based on the findings and experiences.An on-line validation tool,GLCVal,was developed by integrating land cover validation procedures with the service computing technologies.About 20 countries(regions)have completed the accuracy assess-ment of GlobeLand30 for their territories with the guidance of the technical specification and the support of GLCVal.     

Comparison and relative quality assessment of the GLC2000, GLOBCOVER,MODIS and ECOCLIMAP land cover data sets at the African continental scale

Land cover dynamics at the African continental scale is of great importance for global change studies. Actually, four satellite-derived land cover maps of Africa now available, e.g. ECOCLIMAP, GLC2000, MODIS and GLOBCOVER, are based on images acquired in the 2000s. This study aims at stressing the compliances and the discrepancies between these four land cover classifications systems. Each of them used different mapping initiatives and relies on different mapping standards, which supports the present investigation. In order to do a relative comparison of the four maps, a preamble was to reconcile their thematic legends into more aggregated categories after a projection into the same spatial resolution. Results show that the agreement between the four land cover products is between 56 and 69%. While all these land cover datasets show a reasonable agreement in terms of surface types and spatial distribution patterns, mapping of heterogeneous landscapes in the four products is not very successful. Land cover products based on remote sensing imagery can indeed significantly be improved by using smarter algorithms, better timing of image acquisition, improved class definitions. Either will help to improve the accuracy of future land cover maps at the African continental scale. Data producers may use the areas of spatial agreement for training area selection while users might need to verify the information in the areas of disagreement using additional data sources.     

Linear unmixing of MODIS albedo composites to infer subpixel land cover type albedos

We use a linear unmixing approach to test how land use and forestry maps, in combination with the MODIS BRDF/albedo product, can be used to estimate land cover type albedos in boreal regions. Operational land use maps from three test areas in Finland and Canada were used to test the method. The resulting endmember albedo estimates had low standard errors of the mean and were realistic for the main land cover types. The estimated albedos were fairly consistent with albedo measurements conducted with a telescope mast and pure pixel albedos. Problems with the method are the possible errors in the land cover maps, lack of good quality winter MODIS albedo composites and the mismatch between the MODIS pixels and the true observation area. The results emphasize the role of tree species as determinant of forest albedo. Comprehensive spatial and temporal measurements of land cover albedo are usually not possible with in situ mast measurements, and the spatial resolution of MODIS albedo product is often too low to allow direct comparison of pixel albedos and land cover types in areas with heterogeneous vegetation. Hence, and since local forestry maps exist for most temperate and boreal regions, we believe that the proposed method will be useful in estimating average regional land cover type albedos as well as in tracking changes in them.     

Land cover classification of the North China Plain using MODIS_EVI time series

Moderate Resolution Imaging Spectroradiometer (MODIS) data have played an important role in global environmental and resource research. However, its low spatial resolution has been an impediment to researchers pursuing more accurate classification results. In this research, the high temporal resolution of MODIS was employed to improve the accuracy of land cover classification of the North China Plain using MODIS_EVI time series from 2003. Harmonic Analysis of Time Series (HANTS) was performed on the MODIS_EVI image time series to reduce cloud and other noise effects. The improved MODIS_EVI time series was then classified into 100 clusters by the Iterative Self-Organizing Data Analysis Technique (ISODATA). To distinguish ambiguous land cover classes, a decision tree was built on five phenological features derived from EVI profiles, Land Surface Temperature (LST) and topographic slope. The overall accuracy of the final land cover map was 75.5%, indicating the promise of using MODIS EVI time series and decision trees for broad area land cover classification.     

Accuracy assessment of the global forest watch tree cover 2000 in China

Inaccurate information on forest resources could hamper forest conservation, reforestation and sustainable management. Remote-sensing products have emerged as key tools in forest cover monitoring. The Global Forest Watch (GFW) dataset as an interactive remote sensing product, is now applied by more than 2 million users including researchers, conservationists and local communities for analyzing forest cover changes. The quality of this product varies spatially, and local validations are recommended before using the data for inventory and management tasks. Our study evaluated the accuracy and suitability of the GFW dataset for analyzing China’s forest cover. We conducted a validation based on a streamlined visual interpretation procedure using high-resolution optical imagery on Google Earth to map the uncertainties and inaccuracies of GFW Tree Cover 2000 in China. We then estimated China’s forest area after considering the data uncertainty, made a comparison with the data reported by the National Forest Inventory of China (CNFI) to understand where and how the land-based inventory differs from the presence/absence-based remote sensing data. The results showed that the overall accuracy of the GFW Tree Cover 2000 data reached 94.5 %. The user’s and producer’s accuracy of forest classification was 89.26 % and 82.13 %. The sample-based area estimation using GFW showed a larger forest area than the figure reported by CNFI in mainland China, while data discrepancy varied at provincial levels. The study provides a detailed performance assessment of GFW in terms of accuracy of defining forest, and we advise the consideration of data uncertainty in forest cover estimates for future forest management.     

Evaluating land cover changes in Eastern and Southern Africa from 2000 to 2010 using validated Landsat and MODIS data

In this study, we assessed land cover land use (LCLU) changes and their potential environmental drivers (i.e., precipitation, temperature) in five countries in Eastern & Southern (E&S) Africa (Rwanda, Botswana, Tanzania, Malawi and Namibia) between 2000 and 2010. Landsat-derived LCLU products developed by the Regional Centre for Mapping of Resources for Development (RCMRD) through the SERVIR (Spanish for “to serve”) program, a joint initiative of NASA and USAID, and NASA’s Moderate Resolution Imaging Spectroradiometer (MODIS) data were used to evaluate and quantify the LCLU changes in these five countries. Given that the original development of the MODIS land cover type standard products included limited training sites in Africa, we performed a two-level verification/validation of the MODIS land cover product in these five countries. Precipitation data from CHIRPS dataset were used to evaluate and quantify the precipitation changes in these countries and see if it was a significant driver behind some of these LCLU changes. MODIS Land Surface Temperature (LST) data were also used to see if temperature was a main driver too.Our validation analysis revealed that the overall accuracies of the regional MODIS LCLU product for this African region alone were lower than that of the global MODIS LCLU product overall accuracy (63–66% vs. 75%). However, for countries with uniform or homogenous land cover, the overall accuracy was much higher than the global accuracy and as high as 87% and 78% for Botswana and Namibia, respectively. In addition, the wetland and grassland classes had the highest user’s accuracies in most of the countries (89%–99%), which are the ones with the highest number of MODIS land cover classification algorithm training sites.Our LCLU change analysis revealed that Botswana’s most significant changes were the net reforestation, net grass loss and net wetland expansion. For Rwanda, although there have been significant forest, grass and crop expansions in some areas, there also have been significant forest, grass and crop loss in other areas that resulted in very minimal net changes. As for Tanzania, its most significant changes were the net deforestation and net crop expansion. Malawi’s most significant changes were the net deforestation, net crop expansion, net grass expansion and net wetland loss. Finally, Namibia’s most significant changes were the net deforestation and net grass expansion.The only noticeable environmental driver was in Malawi, which had a significant net wetland loss and could be due to the fact that it was the only country that had a reduction in total precipitation between the periods when the LCLU maps were developed. Not only that, but Malawi also happened to have a slight increase in temperature, which would cause more evaporation and net decrease in wetlands if the precipitation didn’t increase as was the case in that country. In addition, within our studied countries, forestland expansion and loss as well as crop expansion and loss were happening in the same country almost equally in some cases. All of that implies that non-environmental factors, such as socioeconomics and governmental policies, could have been the main drivers of these LCLU changes in many of these countries in E&S Africa. It will be important to further study in the future the detailed effects of such drivers on these LCLU changes in this part of the world.     

Accuracy assessment of four cloud-free snow cover products over the Qinghai-Tibetan Plateau

Four up-to-date daily cloud-free snow products – IMS (Interactive Multisensor Snow products), MOD-SSM/I (combination of the MODIS and SSM/I snow products), MOD-B (Blending method basing on the MODIS snow cover products) and TAI (Terra–Aqua–IMS) – with high-resolutions over the Qinghai-Tibetan Plateau (QTP) were comprehensively assessed. Comparisons of the IMS, MOD-SSM/I, MOD-B and TAI cloud-free snow products against meteorological stations observations over 10 snow seasons (2004–2013) over the QTP indicated overall accuracies of 76.0%, 89.3%, 92.0% and 92.0%, respectively. The K_hat values of the IMS, MOD-SSM/I, MOD-B and TAI products were 0.084, 0.463, 0.428 and 0.526, respectively. The TAI products appear to have the best cloud-removal ability among the four snow products over the QTP. Based on the assessment, an I-TAI (Improvement of Terra–Aqua–IMS) snow product was proposed, which can improve the accuracy to some extent. However, the algorithms of the MODIS series products show instability when identifying wet snow and snow under forest cover over the QTP. The snow misclassification is an important limitation of MODIS snow cover products and requires additional improvements.     

Quantifying urban land cover change between 2001 and 2006 in the Gulf of Mexico region

We estimated urbanization rates (2001–2006) in the Gulf of Mexico region using the National Land Cover Database (NLCD) 2001 and 2006 impervious surface products. An improved method was used to update the NLCD impervious surface product in 2006 and associated land cover transition between 2001 and 2006. Our estimation reveals that impervious surface increased 416 km² with a growth rate of 5.8% between 2001 and 2006. Approximately 1110.1 km² of non-urban lands were converted into urban land, resulting in a 3.2% increase in the region. Hay/pasture, woody wetland, and evergreen forest represented the three most common land cover classes that transitioned to urban. Among these land cover transitions, more than 50% of the urbanization occurred within 50 km of the coast. Our analysis shows that the close-to-coast land cover transition trend, especially within 10 km off the coast, potentially imposes substantial long-term impacts on regional landscape and ecological conditions.