Land cover in Beijing experienced a dramatic change due to intensive human activities, such as urbanization and afforestation. However, the spatial patterns of the dynamics are still unknown. The ...archived Landsat images provide an unprecedented opportunity to detect land cover changes over the past three decades. In this study, we used the Normalized Difference Vegetation Index (NDVI) trajectory to detect major land cover dynamics in Beijing. Then, we classified the land cover types in 2015 with the Google Earth Engine (GEE) cloud calculation. By overlaying the latest land cover types and the spatial distribution of land cover dynamics, we determined the main types where a land cover change occurred. The overall change detection accuracy for three types (vegetation loss associated with negative change in NDVI, vegetation gain associated with positive change in NDVI, and no changes) is 86.13%. We found that the GEE is a fast and powerful tool for land cover mapping, and we obtained a classification map with an overall accuracy of 86.61%. Over the past 30years, 1402.28km2 of land was with vegetation loss and 1090.38km2 of land was revegetated in Beijing. The spatial pattern of vegetation loss and vegetation gain shows significant differences in different zones from the center of the city. We also found that 1162.71km2 of land was converted to urban and built-up, whereas 918.36km2 of land was revegetated to cropland, shrub land, forest, and grassland. Moreover, 202.67km2 and 156.75km2 of the land was transformed to forest and shrub land in the plain of Beijing that were traditionally used for cropland and housing.
•We used annual NDVI time-series to detect the land cover dynamics in Beijing.•We used cloud computation in the GEE to map the most recent land cover in 2015.•We obtained a classification map with an overall accuracy of 86.61%.•We found vegetation loss and vegetation gain patterns over the past three decades.•1402.28km2 of land was with vegetation loss and 1090.38km2 was revegetated.
High quality training samples are essential for global land cover mapping. Traditionally, training samples are collected by field work or via manual interpretation based on high-resolution Google ...Earth images. Due to the difficulty of training sample collection, regular global land cover mapping is still a challenge. In this study, we developed an automatic training sample migration method based on the first all-season sample set in 2015 and all available archived Landsat 5 TM images in the Google Earth Engine cloud-based platform. By measuring the spectral similarity and spectral distance between the reference spectral and image spectral, we detected and identified the change state of training sample pixels in 2010, 2005, 2000, 1995, and 1990. Overall, 170,925 (66%), 118,586 (64%), 112,092 (67%), 154,931 (63%), and 147,267 (60%) respective training sample pixels were found with no changes over each five-year period. The detection (user's) accuracies of migrated training sample pixels as no change for the first four time periods were 99.25%, 97.65%, 95.03%, and 92.98%, respectively, by comparing with CCI-LC (Climate Change Initiative Land Cover) maps. Classification experiment showed that the migrated training samples can obtain a similar classification accuracy of 71.42% in 2010, when compared to the classification result in 2015 using the same number of training samples. Our study provides a potential solution to resolve the problem of lack of training samples for dynamic global land cover mapping efforts.
Urban boundaries, an essential property of cities, are widely used in many urban studies. However, extracting urban boundaries from satellite images is still a great challenge, especially at a global ...scale and a fine resolution. In this study, we developed an automatic delineation framework to generate a multi-temporal dataset of global urban boundaries (GUB) using 30 m global artificial impervious area (GAIA) data. First, we delineated an initial urban boundary by filling inner non-urban areas of each city. A kernel density estimation approach and cellular-automata based urban growth modeling were jointly used in this step. Second, we improved the initial urban boundaries around urban fringe areas, using a morphological approach by dilating and eroding the derived urban extent. We implemented this delineation on the Google Earth Engine platform and generated a 30 m resolution global urban boundary dataset in seven representative years (i.e. 1990, 1995, 2000, 2005, 2010, 2015, and 2018). Our extracted urban boundaries show a good agreement with results derived from nighttime light data and human interpretation, and they can well delineate the urban extent of cities when compared with high-resolution Google Earth images. The total area of 65 582 GUBs, each of which exceeds 1 km2, is 809 664 km2 in 2018. The impervious surface areas account for approximately 60% of the total. From 1990 to 2018, the proportion of impervious areas in delineated boundaries increased from 53% to 60%, suggesting a compact urban growth over the past decades. We found that the United States has the highest per capita urban area (i.e. more than 900 m2) among the top 10 most urbanized nations in 2018. This dataset provides a physical boundary of urban areas that can be used to study the impact of urbanization on food security, biodiversity, climate change, and urban health. The GUB dataset can be accessed from http://data.ess.tsinghua.edu.cn.
South Asia experienced a weakening of summer monsoon circulation in the past several decades, resulting in rainfall decline in wet regions. In comparison with other tropical ecosystems, quantitative ...assessments of the extent and triggers of vegetation change are lacking in assessing climate‐change impacts over South Asia dominated by crops. Here, we use satellite‐based Normalized Difference Vegetation Index (NDVI) to quantify spatial–temporal changes in vegetation greenness, and find a widespread annual greening trend that stands in contrast to the weakening of summer monsoon circulation particularly over the last decade. We further show that moisture supply is the primary factor limiting vegetation activity during dry season or in dry region, and cloud cover or temperature would become increasingly important in wet region. Enhanced moisture conditions over dry region, coinciding with the decline in monsoon, are mainly responsible for the widespread greening trend. This result thereby cautions the use of a unified monsoon index to predict South Asia's vegetation dynamics. Current climate–carbon models in general correctly reproduce the dominant control of moisture in the temporal characteristics of vegetation productivity. But the model ensemble cannot exactly reproduce the spatial pattern of satellite‐based vegetation change mainly because of biases in climate simulations. The moisture‐induced greening over South Asia, which is likely to persist into the wetter future, has significant implications for regional carbon cycling and maintaining food security.
Spatial distributions of GIMMS NDVI trends over the South Asia during the three periods: 1982–2014 (first row), 1982–2001 (second row), and 2002–2014 (third row). The trends are estimated on the annual (a, d, and g), wet season (b, e, and h), and dry season basis (c, f, and i), respectively. The inset panels show the pixels where NDVI trends are statistically significant at p < .05.
The spatiotemporal distribution of ecosystem service values (ESVs) and ecological risk are critical indicators to represent the regional ecological protection level and potential of sustainable ...development, which largely depend on land-use patterns. Aiming to contribute to global climate mitigation, China has proposed dual-carbon goals that would remarkably influence the land-use/cover change (LUCC) distribution. Based on the Landsat land cover data of 2000, 2010 and 2020 and multisource satellite products, several driving factors are integrated into the patch-generating land use simulation (PLUS) model to simulate future LUCC patterns for the Guangdong–Hong Kong–Macao Greater Bay Area (GBA) under rapid urbanization, cropland protection and carbon neutral (CN) scenarios from 2020 to 2050. Spatial–temporal ecosystem service and ESVs are allocated using INVEST and the equivalent factor method and thus ecological risks are evaluated using the entropy method. Results indicate that forest growth is the largest under the CN scenario, especially in the northwestern and northeastern GBA, exceeding 25,800 km2 in 2050, which results in both the highest habitat quality and carbon storage. The largest ESVs, reaching higher than 5210 yuan/pixel, are found in the CN scenario, particularly expanding toward the suburban area, leading to the lowest ecological risks. From 2020 to 2050, habitat quality, carbon storage and ESVs improve, while ecological risks decline in the CN scenario. This research provides implications for economic and ecological balanced development and gives references to the carbon-neutral pathway for the GBA.
As Arctic warming continues, its impact on vegetation greenness is complex, variable and inherently scale-dependent. Studies with multiple spatial resolution satellite observations, with 30 m ...resolution included, on tundra greenness have been implemented all over the North American tundra. However, finer resolution studies on the greenness trends in the Russian tundra have only been carried out at a limited local or regional scale and the spatial heterogeneity of the trend remains unclear. Here, we analyzed the fine spatial resolution dataset Landsat archive from 1984 to 2018 over the entire Russian tundra and produced pixel-by-pixel greenness trend maps with the support of Google Earth Engine (GEE). The entire Russian tundra was divided into six geographical regions based on World Wildlife Fund (WWF) ecoregions. A Theil–Sen regression (TSR) was used for the trend identification and the changed pixels with a significance level p < 0.05 were retained in the final results for a subsequent greening/browning trend analysis. Our results indicated that: (1) the number of valid Landsat observations was spatially varied. The Western and Eastern European Tundras (WET and EET) had denser observations than other regions, which enabled a trend analysis during the whole study period from 1984 to 2018; (2) the most significant greening occurred in the Yamal-Gydan tundra (WET), Bering tundra and Chukchi Peninsula tundra (CT) during 1984–2018. The EET had a greening trend of 2.3% and 6.6% and the WET of 3.4% and 18% during 1984–1999 and 2000–2018, respectively. The area of browning trend was relatively low when we first masked the surface water bodies out before the trend analysis; and (3) the Landsat-based greenness trend was broadly similar to the AVHRR-based trend over the entire region but AVHRR retrieved more browning areas due to spectral mixing adjacent effects. Higher resolution images and field measurement studies are strongly needed to understand the vegetation trend over the Russian tundra ecosystem.
The last two decades have witnessed increasing awareness of the potential of terrestrial laser scanning (TLS) in forest applications in both public and commercial sectors, along with tremendous ...research efforts and progress. It is time to inspect the achievements of and the remaining barriers to TLS-based forest investigations, so further research and application are clearly orientated in operational uses of TLS. In such context, the international TLS benchmarking project was launched in 2014 by the European Spatial Data Research Organization and coordinated by the Finnish Geospatial Research Institute. The main objectives of this benchmarking study are to evaluate the potential of applying TLS in characterizing forests, to clarify the strengths and the weaknesses of TLS as a measure of forest digitization, and to reveal the capability of recent algorithms for tree-attribute extraction. The project is designed to benchmark the TLS algorithms by processing identical TLS datasets for a standardized set of forest attribute criteria and by evaluating the results through a common procedure respecting reliable references. Benchmarking results reflect large variances in estimating accuracies, which were unveiled through the 18 compared algorithms and through the evaluation framework, i.e., forest complexity categories, TLS data acquisition approaches, tree attributes and evaluation procedures. The evaluation framework includes three new criteria proposed in this benchmarking and the algorithm performances are investigated through combining two or more criteria (e.g., the accuracy of the individual tree attributes are inspected in conjunction with plot-level completeness) in order to reveal algorithms’ overall performance. The results also reveal some best available forest attribute estimates at this time, which clarify the status quo of TLS-based forest investigations. Some results are well expected, while some are new, e.g., the variances of estimating accuracies between single-/multi-scan, the principle of the algorithm designs and the possibility of a computer outperforming human operation. With single-scan data, i.e., one hemispherical scan per plot, most of the recent algorithms are capable of achieving stem detection with approximately 75% completeness and 90% correctness in the easy forest stands (easy plots: 600 stems/ha, 20 cm mean DBH). The detection rate decreases when the stem density increases and the average DBH decreases, i.e., 60% completeness with 90% correctness (medium plots: 1000 stem/ha, 15 cm mean DBH) and 30% completeness with 90% correctness (difficult plots: 2000 stems/ha, 10 cm mean DBH). The application of the multi-scan approach, i.e., five scans per plot at the center and four quadrant angles, is more effective in complex stands, increasing the completeness to approximately 90% for medium plots and to approximately 70% for difficult plots, with almost 100% correctness. The results of this benchmarking also show that the TLS-based approaches can provide the estimates of the DBH and the stem curve at a 1–2 cm accuracy that are close to what is required in practical applications, e.g., national forest inventories (NFIs). In terms of algorithm development, a high level of automation is a commonly shared standard, but a bottleneck occurs at stem detection and tree height estimation, especially in multilayer and dense forest stands. The greatest challenge is that even with the multi-scan approach, it is still hard to completely and accurately record stems of all trees in a plot due to the occlusion effects of the trees and bushes in forests. Future development must address the redundant yet incomplete point clouds of forest sample plots and recognize trees more accurately and efficiently. It is worth noting that TLS currently provides the best quality terrestrial point clouds in comparison with all other technologies, meaning that all the benchmarks labeled in this paper can also serve as a reference for other terrestrial point clouds sources.
Abstract
Spatial data of urban green spaces (UGS) are critical for cities worldwide to evaluate their progress towards achieving the urban sustainable development goals on UGS. However, UGS maps at ...the global scale with acceptable accuracies are not readily available. In this study, we mapped UGS of all 1039 mid- and large-sized cities across the globe in 2015 with dense remote sensing data (i.e. 51 494 Landsat images) and Google Earth Engine (GEE) platform. Also, we quantified the spatial distribution and accessibility of UGS within the cities. By combining the greenest pixel compositing method and the percentile-based image compositing method, we were able to obtain the maximum extent of UGS in cities while better differentiating UGS from other vegetation such as croplands. The mean overall classification accuracy reached 89.26% (SD = 3.26%), which was higher than existing global land cover products. Our maps showed that the mean UGS coverage in 1039 cities was 38.46% (SD = 20.27%), while the mean UGS accessibility was 82.67% (SD = 22.89%). However, there was a distinctive spatial equity issue as cities in high-income countries had higher coverage and better accessibility than cities in low-income countries. Besides developing a protocol for large-scale UGS mapping, our study results provide key baseline information to support international endeavors to fulfill the relevant urban sustainable development goals.
Urban green spaces can yield considerable health benefits to urban residents. Assessing these health benefits is a key step for managing urban green spaces for human health and wellbeing in cities. ...In this study, we assessed the change of health benefits generated by urban green spaces in 28 megacities worldwide between 2005 and 2015 by using availability and accessibility as proxy indicators. We first mapped land covers of 28 megacities using 10,823 scenes of Landsat images and a random forest classifier running on Google Earth Engine. We then calculated the availability and accessibility of urban green spaces using the land cover maps and gridded population data. The results showed that the mean availability of urban green spaces in these megacities increased from 27.63% in 2005 to 31.74% in 2015. The mean accessibility of urban green spaces increased from 65.76% in 2005 to 72.86% in 2015. The increased availability and accessibility of urban green spaces in megacities have brought more health benefits to their residents.
•A framework is proposed to identify classification anomalies in large-scale land cover datasets without ground truth.•Multi-level randomness are built in to reduce the influence of inherent ...uncertainty in target dataset.•No prior knowledge of data distribution and original classifiers is required.
In the last decade, land cover products are produced at a global scale and updated with an unprecedent speed with the development of earth observation and mapping techniques. However, assessing large-scale land cover datasets is always a challenging task because of lack of ground truth and high dependency on manual inspection. To promote the efficiency of assessment, this study introduces a generic framework that identifies potential land cover classification errors without traditional reference data. Depending on rich features in remote sensing data, the overall procedure of can be regarded as a multi-class anomaly detection problem. To improve the performance in dealing with complex classification scheme, a pairing strategy is firstly proposed for pairwise analysis on each two classes, so that the multi-class problem is decomposed into multiple smaller binary-class problems. Secondly, the proximity matrix of each class pair is generated on the basis of extra trees model, which is applicable to mixed numerical and categorical data and provides robust proximity measurement by including multiple levels of randomness. Finally, an improved proximity-based anomaly detection algorithm is applied, and the detection results of each class pair are ensembled to obtain the final anomaly score of each land instance. Experiments show that the overall area under the receiver operating characteristic curve score reaches approximately 0.9 for synthetic datasets and 0.8 for real-life datasets. The proposed method is expected to provide users valuable data for a more efficient assessment in the absence of ground truth in practice.