Determination of Vegetation Thresholds for Assessing Land Use and Land Use Changes in Cambodia using the Google Earth Engine Cloud-Computing Platform

Leveraging Google Earth Engine for a More Effective Grassland Management: A Decision Support Application Perspective

Grasslands cover a substantial portion of the earth's surface and agricultural land and is crucial for human well-being and livestock farming. Ranchers and grassland management authorities face challenges in effectively controlling herders' grazing behavior and grassland utilization due to underdeveloped infrastructure and poor communication in pastoral areas. Cloud-based grazing management and decision support systems (DSS) are needed to address this issue, promote sustainable grassland use, and preserve their ecosystem services. These systems should enable rapid and large-scale grassland growth and utilization monitoring, providing a basis for decision-making in managing grazing and grassland areas. In this context, this study contributes to the objectives of the EU LIFE IMAGINE project, aiming to develop a Web-GIS app for conserving and monitoring Umbria's grasslands and promoting more informed decisions for more sustainable livestock management. The app, called "Praterie" and developed in Google Earth Engine, utilizes historical Sentinel-2 satellite data and harmonic modeling of the EVI (Enhanced Vegetation Index) to estimate vegetation growth curves and maturity periods for the forthcoming vegetation cycle. The app is updated in quasi-real time and enables users to visualize estimates for the upcoming vegetation cycle, including the maximum greenness, the days remaining to the subsequent maturity period, the accuracy of the harmonic models, and the grassland greenness status in the previous 10 days. Even though future additional developments can improve the informative value of the Praterie app, this platform can contribute to optimizing livestock management and biodiversity conservation by providing timely and accurate data about grassland status and growth curves.

Evaluation of urban green space per capita with new remote sensing and geographic information system techniques and the importance of urban green space during the COVID-19 pandemic

A recently conducted study by the Centers for Disease Control and Prevention encouraged access to urban green space for the public over the prevalence of COVID-19 in that exposure to urban green space can positively affect the physical and mental health, including the reduction rate of heart disease, obesity, stress, stroke, and depression. COVID-19 has foregrounded the inadequacy of green space in populated cities. It has also highlighted the extant inequities so as to unequal access to urban green space both quantitatively and qualitatively. In this regard, it seems that one of the problems related to Malatya is the uncoordinated distribution of green space in different parts of the city. Therefore, knowing the quantity and quality of these spaces in each region can play an effective role in urban planning. The aim of the present study has been to evaluate urban green space per capita and to investigate its distribution based on the population of the districts of Battalgazi county in Malatya city through developing an integrated methodology (remote sensing and geographic information system). Accordingly, in Google Earth Engine by images of Sentinel-1 and PlanetScope satellites, it was calculated different indexes (NDVI, EVI, PSSR, GNDVI, and NDWI). The data set was prepared and then by combining different data, classification was performed according to support vector machine algorithm. From the landscaping maps obtained, the map was selected with the highest accuracy (overall accuracy: 94.43; and kappa coefficient: 90.5). Finally, by the obtained last map, the distribution of urban green space per capita and their functions in Battalgazi county and its districts were evaluated. The results of the study showed that the existing urban green spaces in the Battalgazi/Malatya were not distributed evenly on the basis of the districts. The per capita of urban green space is twenty-four regions which is more than 9m² and in twenty-three ones is less than 9m². The recommendation of this study was that Türkiye city planners and landscape designers should replan and redesign the quality and equal distribution of urban green spaces, especially during and following COVID-19 pandemic. Additionally, drawing on the Google Earth Engine cloud system, which has revolutionized GIS and remote sensing, is recommended to be used in land use land cover modeling. It is straightforward to access information and analyze them quickly in Google Earth Engine. The published codes in this study makes it possible to conduct further relevant studies.

Assessment of Forest Cover Changes in Vavuniya District, Sri Lanka: Implications for the Establishment of Subnational Forest Reference Emission Level

Assessment of forest cover changes is required to establish the forest reference emission level (FREL) at any scale. Due to civil conflict, such assessments have not yet been undertaken in Sri Lanka, especially in the conflict zone. Here, we assessed the forest cover changes in Vavuniya District, Sri Lanka, from 2001 to 2020, using a combination of the Google Earth Engine (GEE) platform and the phenology-based threshold classification (PBTC) method. Landsat 5 TM data for 2001, 2006, and 2010, and Landsat 8 OLI data for 2016 and 2020 were used to classify forest cover by categories, and their related changes could be assessed by four categories, namely dry monsoon forest, open forest, other lands, and water bodies. With an overall average accuracy of 87% and an average kappa coefficient of 0.83, forest cover was estimated at 57.6% of the total land area in 2020. There was an increase of 0.46% per annum for the entire district between 2001 and 2010, but a drastic loss of 0.60% per year was observed between 2010 and 2020. Specifically, the dry monsoon forest lost 0.30%, but open forest gained 3.62% annually over the same period. Loss and gain of forest cover resulted in carbon emissions and removals of 165,306.6 MgCO2 and 24,064.5 MgCO2 annually, respectively, over the same period. Our findings could be used to set the baseline trend of deforestation, based on which, a subnational forest reference emission level can be established as an emission benchmark, against which comparisons of carbon emissions following the implementation of REDD+ activities can be made, and result-based payment can be claimed under the Paris Agreement.

AUnet: A Deep Learning Framework for Surface Water Channel Mapping Using Large-Coverage Remote Sensing Images and Sparse Scribble Annotations from OSM Data

Water is a vital component of life that exists in a variety of forms, including oceans, rivers, ponds, streams, and canals. The automated methods for detecting, segmenting, and mapping surface water have improved significantly with the advancements in satellite imagery and remote sensing. Many strategies and techniques to segment water resources have been presented in the past. However, due to the variant width and complex appearance, the segmentation of the water channel remains challenging. Moreover, traditional supervised deep learning frameworks have been restricted by the scarcity of water channel datasets that include precise water annotations. With this in mind, this research presents the following three main contributions. Firstly, we curated a new dataset for water channel mapping in the Pakistani region. Instead of employing pixel-level water channel annotations, we used a weakly trained method to extract water channels from VHR pictures, relying only on OpenStreetMap (OSM) waterways to create sparse scribbling annotations. Secondly, we benchmarked the dataset on state-of-the-art semantic segmentation frameworks. We also proposed AUnet, an atrous convolution inspired deep learning network for precise water channel segmentation. The experimental results demonstrate the superior performance of the proposed AUnet model for segmenting using weakly supervised labels, where it achieved a mean intersection over union score of 0.8791 and outperformed state-of-the-art approaches by 5.90% for the extraction of water channels.

Two-Stepwise Hierarchical Adaptive Threshold Method for Automatic Rapeseed Mapping over Jiangsu Using Harmonized Landsat/Sentinel-2

Rapeseed distribution mapping is a crucial issue for food and oil security, entertainment, and tourism development. Previous studies have used various remote sensing approaches to map rapeseed. However, the time-consuming and labor-intensive sample data used in these supervised classification methods greatly limit the development of large-scale mapping in rapeseed studies. Regarding threshold methods, some empirical thresholding methods still need sample data to select the optimal threshold value, and their accuracies decrease when a fixed threshold is applied in complex and diverse environments. This study first developed the Normalized Difference Rapeseed Index (NDRI), defined as the difference in green and short-wave infrared bands divided by their sum, to find a suitable feature to distinguish rapeseed from other types of crops. Next, a two-stepwise hierarchical adaptive thresholding (THAT) algorithm requiring no training data was used to automatically extract rapeseed in Xinghua. Finally, two adaptive thresholding methods of the standalone Otsu and Otsu with Canny Edge Detection (OCED) were used to extract rapeseed across Jiangsu province. The results show that (1) NDRI can separate rapeseed from other vegetation well; (2) the OCED-THAT method can accurately map rapeseed in Jiangsu with an overall accuracy (OA) of 0.9559 and a Kappa coefficient of 0.8569, and it performed better than the Otsu-THAT method; (3) the OCED-THAT method had a lower but acceptable accuracy than the Random Forest method (OA = 0.9806 and Kappa = 0.9391). This study indicates that the THAT model is a promising automatic method for mapping rapeseed.

Characterization of the main land processes occurring in Europe (2000-2018) through a MODIS NDVI seasonal parameter-based procedure

The identification and recognition of the land processes are of vital importance for a proper management of the ecosystem functions and services. However, on-ground land uses/land covers (LULC) characterization is a time-consuming task, often limited to small land areas, which can be solved using remote sensing technologies. The objective of this work is to investigate how the different MODIS NDVI seasonal parameters responded to the main land processes observed in Europe in the 2000-2018 period; characterizing their temporal trend; and evaluating which one reflected better each specific land process. NDVI time-series were evaluated using TIMESAT software, which extracted eight seasonality parameters: amplitude, base value, length of season, maximum value, left and right derivative values and small and large integrated values. These parameters were correlated with the LULC changes derived from COoRdination of INformation on the Environment Land Cover (CLC) for assessing which parameter better characterized each land process. The temporal evolution of the maximum seasonal NDVI was the parameter that better characterized the occurrence of most of the land processes evaluated (afforestation, agriculturalization, degradation, land abandonment, land restoration, urbanization; R² from 0.67-0.97). Large integrated value also presented significant relationships but they were restricted to two of the three evaluated periods. On the contrary, land processes involving CLC categories with similar NDVI patterns were not well captured with the proposed methodology. These results evidenced that this methodology could be combined with other classification methods for improving LULC identification accuracy or for identifying LULC processes in locations where no LULC maps are available. Such information can be used by policy-makers to draw LULC management actions associated with sustainable development goals. This is especially relevant for areas where food security is at stake and where terrestrial ecosystems are threatened by severe biodiversity loss.

Datasets of drought and flood impact on croplands in Southeast Asia from 1980 to 2019

Data on droughts and floods and their impacts on croplands and production are important for policy makers and the scientific community. This dataset was developed to provide data of the impacts of droughts and floods on agriculture in the Monsoon Climate Region and Equatorial Climate Region of Southeast Asia during the crop growing seasons over a 40-year period between 1980 and 2019. The data were generated using the TerraClimate global high-resolution gridded Palmer Drought Severity Index (PDSI) datasets in Google Earth Engine along with a set of algorithms. Datasets were available on 47,192 grid points of a 10 × 10 km resolution containing PDSI, their latitude longitude between 1980 and 2019 with five years interval, monthly temporal PDSI data, cropland drought and flood intensity data between 1980 and 2019.

Spatiotemporal Analysis of Land Cover and the Effects on Ecosystem Service Values in Rupandehi, Nepal from 2005 to 2020

Land cover (LC) is a crucial parameter for studying environmental phenomena. Cutting-edge technology such as remote sensing (RS) and cloud computing have made LC change mapping efficient. In this study, the LC of Rupandehi District of Nepal were mapped using Landsat imagery and Random Forest (RF) classifier from 2005 to 2020 using Google Earth Engine (GEE) platform. GEE eases the way in extracting, analyzing, and performing different operations for the earth’s observed data. Land cover classification, Centre of gravity (CoG), and their trajectories for all LC classes: agriculture, built-up, water, forest, and barren area were extracted with five-year intervals, along with their Ecosystem service values (ESV) to understand the load on the ecosystem. We also discussed the aspects and problems of the spatiotemporal analysis of developing regions. It was observed that the built-up areas had been increasing over the years and more centered in between the two major cities. Other agriculture, water, and forest classes had been subjected to fluctuations with barren land in the decreasing trend. This alteration in the area of the LC classes also resulted in varying ESVs for individual land cover and total values for the years. The accuracy for the RF classifier was under substantial agreement for such fragmented LCs. Using LC, CoG, and ESV, the paper discusses the need for spatiotemporal analysis studies in Nepal to overcome the current limitations and later expansion to other regions. Studies such as these help in implementing proper plans and strategies by district administration offices and local governmental bodies to stop the exploitation of resources.

Mangrove Forest Cover and Phenology with Landsat Dense Time Series in Central Queensland, Australia

Wetlands are one of the most biologically productive ecosystems. Wetland ecosystem services, ranging from provision of food security to climate change mitigation, are enormous, far outweighing those of dryland ecosystems per hectare. However, land use change and water regulation infrastructure have reduced connectivity in many river systems and with floodplain and estuarine wetlands. Mangrove forests are critical communities for carbon uptake and storage, pollution control and detoxification, and regulation of natural hazards. Although the clearing of mangroves in Australia is strictly regulated, Great Barrier Reef catchments have suffered landscape modifications and hydrological alterations that can kill mangroves. We used remote sensing datasets to investigate land cover change and both intra- and inter-annual seasonality in mangrove forests in a large estuarine region of Central Queensland, Australia, which encompasses a national park and Ramsar Wetland, and is adjacent to the Great Barrier Reef World Heritage site. We built a time series using spectral, auxiliary, and phenology variables with Landsat surface reflectance products, accessed in Google Earth Engine. Two land cover classes were generated (mangrove versus non-mangrove) in a Random Forest classification. Mangroves decreased by 1480 hectares (−2.31%) from 2009 to 2019. The overall classification accuracies and Kappa coefficient for 2008–2010 and 2018–2020 land cover maps were 95% and 95%, respectively. Using an NDVI-based time series we examined intra- and inter-annual seasonality with linear and harmonic regression models, and second with TIMESAT metrics of mangrove forests in three sections of our study region. Our findings suggest a relationship between mangrove growth phenology along with precipitation anomalies and severe tropical cyclone occurrence over the time series. The detection of responses to extreme events is important to improve understanding of the connections between climate, extreme weather events, and biodiversity in estuarine and mangrove ecosystems.

Multi-Temporal Arable Land Monitoring in Arid Region of Northwest China Using a New Extraction Index

Development of a high-accuracy method to extract arable land using effective data sources is crucial to detect and monitor arable land dynamics, servicing land protection and sustainable development. In this study, a new arable land extraction index (ALEI) based on spectral analysis was proposed, examined by ground truth data, and then applied to the Hexi Corridor in northwest China. The arable land and its change patterns during 1990–2020 were extracted and identified using 40 Landsat TM/OLI images acquired in 1990, 2000, 2010, and 2020. The results demonstrated that the proposed method can distinguish arable land areas accurately, with the User’s (Producer’s) accuracy and overall accuracy (kappa coefficient) exceeding 0.90 (0.88) and 0.89 (0.87), respectively. The mean relative error calculated using field survey data obtained in 2012 and 2020 was 0.169 and 0.191, respectively, indicating the feasibility of the ALEI method in arable land extracting. The study found that arable land area in the Hexi Corridor was 13217.58 km2 in 2020, significantly increased by 25.33% compared to that in 1990. At 10-year intervals, the arable land experienced different change patterns. The study results indicate that ALEI index is a promising tool used to effectively extract arable land in the arid area.

A Classification of Tidal Flat Wetland Vegetation Combining Phenological Features with Google Earth Engine

The composition and distribution of wetland vegetation is critical for ecosystem diversity and sustainable development. However, tidal flat wetland environments are complex, and obtaining effective satellite imagery is challenging due to the high cloud coverage. Moreover, it is difficult to acquire phenological feature data and extract species-level wetland vegetation information by using only spectral data or individual images. To solve these limitations, statistical features, temporal features, and phenological features of multiple Landsat 8 time-series images obtained via the Google Earth Engine (GEE) platform were compared to extract species-level wetland vegetation information from Chongming Island, China. The results indicated that (1) a harmonic model obtained the phenological characteristics of wetland vegetation better than the raw vegetation index (VI) and the Savitzky–Golay (SG) smoothing method; (2) classification based on the combination of the three features provided the highest overall accuracy (85.54%), and the phenological features (represented by the amplitude and phase of the harmonic model) had the greatest impact on the classification; and (3) the classification result from the senescence period was more accurate than that from the green period, but the annual mapping result on all seasons was the most accurate. The method described in this study can be applied to overcome the impacts of the complex environment in tidal flat wetlands and to effectively classify wetland vegetation species using GEE. This study could be used as a reference for the analysis of the phenological features of other areas or vegetation types.

Mask R-CNN and OBIA Fusion Improves the Segmentation of Scattered Vegetation in Very High-Resolution Optical Sensors

Vegetation generally appears scattered in drylands. Its structure, composition and spatial patterns are key controls of biotic interactions, water, and nutrient cycles. Applying segmentation methods to very high-resolution images for monitoring changes in vegetation cover can provide relevant information for dryland conservation ecology. For this reason, improving segmentation methods and understanding the effect of spatial resolution on segmentation results is key to improve dryland vegetation monitoring. We explored and analyzed the accuracy of Object-Based Image Analysis (OBIA) and Mask Region-based Convolutional Neural Networks (Mask R-CNN) and the fusion of both methods in the segmentation of scattered vegetation in a dryland ecosystem. As a case study, we mapped Ziziphus lotus, the dominant shrub of a habitat of conservation priority in one of the driest areas of Europe. Our results show for the first time that the fusion of the results from OBIA and Mask R-CNN increases the accuracy of the segmentation of scattered shrubs up to 25% compared to both methods separately. Hence, by fusing OBIA and Mask R-CNNs on very high-resolution images, the improved segmentation accuracy of vegetation mapping would lead to more precise and sensitive monitoring of changes in biodiversity and ecosystem services in drylands.

Combination of Landsat 8 OLI and Sentinel-1 SAR Time-Series Data for Mapping Paddy Fields in Parts of West and Central Java Provinces, Indonesia

The rise of Google Earth Engine, a cloud computing platform for spatial data, has unlocked seamless integration for multi-sensor and multi-temporal analysis, which is useful for the identification of land-cover classes based on their temporal characteristics. Our study aims to employ temporal patterns from monthly-median Sentinel-1 (S1) C-band synthetic aperture radar data and cloud-filled monthly spectral indices, i.e., Normalized Difference Vegetation Index (NDVI), Modified Normalized Difference Water Index (MNDWI), and Normalized Difference Built-up Index (NDBI), from Landsat 8 (L8) OLI for mapping rice cropland areas in the northern part of Central Java Province, Indonesia. The harmonic function was used to fill the cloud and cloud-masked values in the spectral indices from Landsat 8 data, and smile Random Forests (RF) and Classification And Regression Trees (CART) algorithms were used to map rice cropland areas using a combination of monthly S1 and monthly harmonic L8 spectral indices. An additional terrain variable, Terrain Roughness Index (TRI) from the SRTM dataset, was also included in the analysis. Our results demonstrated that RF models with 50 (RF50) and 80 (RF80) trees yielded better accuracy for mapping the extent of paddy fields, with user accuracies of 85.65% (RF50) and 85.75% (RF80), and producer accuracies of 91.63% (RF80) and 93.48% (RF50) (overall accuracies of 92.10% (RF80) and 92.47% (RF50)), respectively, while CART yielded a user accuracy of only 84.83% and a producer accuracy of 80.86%. The model variable importance in both RF50 and RF80 models showed that vertical transmit and horizontal receive (VH) polarization and harmonic-fitted NDVI were identified as the top five important variables, and the variables representing February, April, June, and December contributed more to the RF model. The detection of VH and NDVI as the top variables which contributed up to 51% of the Random Forest model indicated the importance of the multi-sensor combination for the identification of paddy fields.

Mapping the Natural Distribution of Bamboo and Related Carbon Stocks in the Tropics Using Google Earth Engine, Phenological Behavior, Landsat 8, and Sentinel-2

Although vegetation phenology thresholds have been developed for a wide range of mapping applications, their use for assessing the distribution of natural bamboo and the related carbon stocks is still limited, especially in Southeast Asia. Here, we used Google Earth Engine (GEE) to collect time-series of Landsat 8 Operational Land Imager (OLI) and Sentinel-2 images and employed a phenology-based threshold classification method (PBTC) to map the natural bamboo distribution and estimate carbon stocks in Siem Reap Province, Cambodia. We processed 337 collections of Landsat 8 OLI for phenological assessment and generated 121 phenological profiles of the average vegetation index for three vegetation land cover categories from 2015 to 2018. After determining the minimum and maximum threshold values for bamboo during the leaf-shedding phenology stage, the PBTC method was applied to produce a seasonal composite enhanced vegetation index (EVI) for Landsat collections and assess the bamboo distributions in 2015 and 2018. Bamboo distributions in 2019 were then mapped by applying the EVI phenological threshold values for 10 m resolution Sentinel-2 satellite imagery by accessing 442 tiles. The overall Landsat 8 OLI bamboo maps for 2015 and 2018 had user’s accuracies (UAs) of 86.6% and 87.9% and producer’s accuracies (PAs) of 95.7% and 97.8%, respectively, and a UA of 86.5% and PA of 91.7% were obtained from Sentinel-2 imagery for 2019. Accordingly, carbon stocks of natural bamboo by district in Siem Reap at the province level were estimated. Emission reductions from the protection of natural bamboo can be used to offset 6% of the carbon emissions from tourists who visit this tourism-destination province. It is concluded that a combination of GEE and PBTC and the increasing availability of remote sensing data make it possible to map the natural distribution of bamboo and carbon stocks.

Applications of the Google Earth Engine and Phenology-Based Threshold Classification Method for Mapping Forest Cover and Carbon Stock Changes in Siem Reap Province, Cambodia

Digital and scalable technologies are increasingly important for rapid and large-scale assessment and monitoring of land cover change. Until recently, little research has existed on how these technologies can be specifically applied to the monitoring of Reducing Emissions from Deforestation and Forest Degradation (REDD+) activities. Using the Google Earth Engine (GEE) cloud computing platform, we applied the recently developed phenology-based threshold classification method (PBTC) for detecting and mapping forest cover and carbon stock changes in Siem Reap province, Cambodia, between 1990 and 2018. The obtained PBTC maps were validated using Google Earth high resolution historical imagery and reference land cover maps by creating 3771 systematic 5 × 5 km spatial accuracy points. The overall cumulative accuracy of this study was 92.1% and its cumulative Kappa was 0.9, which are sufficiently high to apply the PBTC method to detect forest land cover change. Accordingly, we estimated the carbon stock changes over a 28-year period in accordance with the Good Practice Guidelines of the Intergovernmental Panel on Climate Change. We found that 322,694 ha of forest cover was lost in Siem Reap, representing an annual deforestation rate of 1.3% between 1990 and 2018. This loss of forest cover was responsible for carbon emissions of 143,729,440 MgCO2 over the same period. If REDD+ activities are implemented during the implementation period of the Paris Climate Agreement between 2020 and 2030, about 8,256,746 MgCO2 of carbon emissions could be reduced, equivalent to about USD 6-115 million annually depending on chosen carbon prices. Our case study demonstrates that the GEE and PBTC method can be used to detect and monitor forest cover change and carbon stock changes in the tropics with high accuracy.

Synergy of Active and Passive Remote Sensing Data for Effective Mapping of Oil Palm Plantation in Malaysia

Oil palm is recognized as a golden crop, as it produces the highest oil yield among oil seed crops. Malaysia is the world’s second largest producer of palm oil; 16% of its land is planted with oil palm. To cope with the ever-increasing global demand on edible oil, additional areas of oil palm are forecast to increase globally by 12 to 19 Mha by 2050. Multisensor remote sensing plays an important role in providing relevant, timely, and accurate information that can be developed into a plantation monitoring system to optimize production and sustainability. The aim of this study was to simultaneously exploit the synthetic aperture radar ALOS PALSAR 2, a form of microwave remote sensing, in combination with visible (red) data from Landsat Thematic Mapper to obtain a holistic view of a plantation. A manipulation of the horizontal–horizontal (HH) and horizontal–vertical (HV) polarizations of ALOS PALSAR data detected oil palm trees and water bodies, while the red spectra L-band from Landsat data (optical) could effectively identify built up areas and vertical–horizontal (VH) polarization from Sentinel C-band data detected bare land. These techniques produced an oil palm area classification with overall accuracies of 98.36% and 0.78 kappa coefficient for Peninsular Malaysia. The total oil palm area in Peninsular Malaysia was estimated to be about 3.48% higher than the value reported by the Malaysian Palm Oil Board. The over estimation may be due the MPOB’s statistics that do not include unregistered small holder oil palm plantations. In this study, we were able to discriminate most of the rubber areas.

Recording Urban Land Dynamic and Its Effects during 2000–2019 at 15-m Resolution by Cloud Computing with Landsat Series

Cities, the core of the global climate change and economic development, are high impact land cover land use change (LCLUC) hotspots. Comprehensive records of land cover land use dynamics in urban regions are essential for strategic climate change adaption and mitigation and sustainable urban development. This study aims to develop a Google Earth Engine (GEE) application for high-resolution (15-m) urban LCLUC mapping with a novel classification scheme using pan-sharpened Landsat images. With this approach, we quantified the annual LCLUC in Changchun, China, from 2000 to 2019, and detected the abrupt changes (turning points of LCLUC). Ancillary data on social-economic status were used to provide insights on potential drivers of LCLUC by examining their correlation with change rate. We also examined the impacts of LCLUC on environment, specifically air pollution. Using this approach, we can classify annual LCLUC in Changchun with high accuracy (all above 0.91). The change detection based on the high-resolution wall-to-wall maps show intensive urban expansion with the compromise of cropland from 2000 to 2019. We also found the growth of green space in urban regions as the result of green space development and management in recent years. The changing rate of different land types were the largest in the early years of the observation period. Turning points of land types were primarily observed in 2009 and 2010. Further analysis showed that economic and industry development and population migration collectively drove the urban expansion in Changchun. Increasing built-up areas could slow wind velocity and air exchange, and ultimately led to the accumulation of PM2.5. Our implement of pan-sharpened Landsat images facilitates the wall-to-wall mapping of temporal land dynamics at high spatial resolution. The primary use of GEE for mapping urban land makes it replicable and transferable by other users. This approach is a first crucial step towards understanding the drivers of change and supporting better decision-making for sustainable urban development and climate change mitigation.

Estimating Proportion of Vegetation Cover at the Vicinity of Archaeological Sites Using Sentinel-1 and -2 Data, Supplemented by Crowdsourced OpenStreetMap Geodata

Monitoring vegetation cover is an essential parameter for assessing various natural and anthropogenic hazards that occur at the vicinity of archaeological sites and landscapes. In this study, we used free and open access to Copernicus Earth Observation datasets. In particular, the proportion of vegetation cover is estimated from the analysis of Sentinel-1 radar and Sentinel-2 optical images, upon their radiometric and geometric corrections. Here, the proportion of vegetation based on the Radar Vegetation Index and the Normalized Difference Vegetation Index is estimated. Due to the medium resolution of these datasets (10 m resolution), the crowdsourced OpenStreetMap service was used to identify fully and non-vegetated pixels. The case study is focused on the western part of Cyprus, whereas various open-air archaeological sites exist, such as the archaeological site of “Nea Paphos” and the “Tombs of the Kings”. A cross-comparison of the results between the optical and the radar images is presented, as well as a comparison with ready products derived from the Sentinel Hub service such as the Sentinel-1 Synthetic Aperture Radar Urban and Sentinel-2 Scene classification data. Moreover, the proportion of vegetation cover was evaluated with Google Earth red-green-blue free high-resolution optical images, indicating that a good correlation between the RVI and NDVI can be generated only over vegetated areas. The overall findings indicate that Sentinel-1 and -2 indices can provide a similar pattern only over vegetated areas, which can be further elaborated to estimate temporal changes using integrated optical and radar Sentinel data. This study can support future investigations related to hazard analysis based on the combined use of optical and radar sensors, especially in areas with high cloud-coverage.

Mapping Annual Land Use and Land Cover Changes in the Yangtze Estuary Region Using an Object-Based Classification Framework and Landsat Time Series Data

A system understanding of the patterns, causes, and trends of long-term land use and land cover (LULC) change at the regional scale is essential for policy makers to address the growing challenges of local sustainability and global climate change. However, it still remains a challenge for estuarine and coastal regions due to the lack of appropriate approaches to consistently generate accurate and long-term LULC maps. In this work, an object-based classification framework was designed to mapping annual LULC changes in the Yangtze River estuary region from 1985–2016 using Landsat time series data. Characteristics of the inter-annual changes of LULC was then analyzed. The results showed that the object-based classification framework could accurately produce annual time series of LULC maps with overall accuracies over 86% for all single-year classifications. Results also indicated that the annual LULC maps enabled the clear depiction of the long-term variability of LULC and could be used to monitor the gradual changes that would not be observed using bi-temporal or sparse time series maps. Specifically, the impervious area rapidly increased from 6.42% to 22.55% of the total land area from 1985 to 2016, whereas the cropland area dramatically decreased from 80.61% to 55.44%. In contrast to the area of forest and grassland, which almost tripled, the area of inland water remained consistent from 1985 to 2008 and slightly increased from 2008 to 2016. However, the area of coastal marshes and barren tidal flats varied with large fluctuations.

Assessing the Importance of Tree Cover Threshold for Forest Cover Mapping Derived from Global Forest Cover in Myanmar

Comprehensive forest cover mapping is essential for making policy and management decisions. However, creating a forest cover map from raw remote sensing data is a barrier for many users. Here, we investigated the effects of different tree cover thresholds on the accuracy of forest cover maps derived from the Global Forest Change Dataset (GFCD) across different ecological zones in a country-scale evaluation of Myanmar. To understand the effect of different thresholds on map accuracy, nine forest cover maps having thresholds ranging from 10% to 90% were created from the GFCD. The accuracy of the forest cover maps within each ecological zone and at the national scale was assessed. The overall accuracies of ecological zones other than tropical rainforest were highest when the threshold for tree cover was less than 50%. The appropriate threshold for tropical rainforests was 80%. Therefore, different optimal tree cover thresholds were required to achieve the highest overall accuracy depending on ecological zones. However, in the unique case of Myanmar, we were able to determine the threshold across the whole country. We concluded that the threshold for tree cover for creating a forest cover map should be determined according to the areal ratio of ecological zones determined from large-scale monitoring. Our results are applicable to tropical regions having similar ecological zones.