Remote sensing (RS) systems have been collecting massive volumes of datasets for decades, managing and analyzing of which are not practical using common software packages and desktop computing ...resources. In this regard, Google has developed a cloud computing platform, called Google Earth Engine (GEE), to effectively address the challenges of big data analysis. In particular, this platform facilitates processing big geo data over large areas and monitoring the environment for long periods of time. Although this platform was launched in 2010 and has proved its high potential for different applications, it has not been fully investigated and utilized for RS applications until recent years. Therefore, this study aims to comprehensively explore different aspects of the GEE platform, including its datasets, functions, advantages/limitations, and various applications. For this purpose, 450 journal articles published in 150 journals between January 2010 and May 2020 were studied. It was observed that Landsat and Sentinel datasets were extensively utilized by GEE users. Moreover, supervised machine learning algorithms, such as Random Forest, were more widely applied to image classification tasks. GEE has also been employed in a broad range of applications, such as Land Cover/land Use classification, hydrology, urban planning, natural disaster, climate analyses, and image processing. It was generally observed that the number of GEE publications have significantly increased during the past few years, and it is expected that GEE will be utilized by more users from different fields to resolve their big data processing challenges.
Vegetation is the main component of the terrestrial Earth, and it plays an imperative role in carbon cycle regulation and surface water/energy exchange/balance. The coupled effects of climate change ...and anthropogenic forcing have undoubtfully impacted the vegetation cover in linear/non-linear manners. Considering the essential benefits of vegetation to the environment, it is vital to investigate the vegetation dynamics through spatially and temporally consistent workflows. In this regard, remote sensing, especially Normalized Difference Vegetation Index (NDVI), has offered a reliable data source for vegetation monitoring and trend analysis. In this paper, two decades (2000 to 2020) of Moderate Resolution Imaging Spectroradiometer (MODIS) NDVI datasets (MOD13Q1) were used for vegetation trend analysis throughout Iran. First, the per-pixel annual NDVI dataset was prepared using the Google Earth Engine (GEE) by averaging all available NDVI values within the growing season and was then fed into the PolyTrend algorithm for linear/non-linear trend identification. In total, nearly 14 million pixels (44% of Iran) were subjected to trend analysis, and the results indicated a higher rate of greening than browning across the country. Regarding the trend types, linear was the dominant trend type with 14%, followed by concealed (11%), cubic (8%), and quadratic (2%), while 9% of the vegetation area remained stable (no trend). Both positive and negative directions were observed in all trend types, with the slope magnitudes ranging between −0.048 and 0.047 (NDVI units) per year. Later, precipitation and land cover datasets were employed to further investigate the vegetation dynamics. The correlation coefficient between precipitation and vegetation (NDVI) was 0.54 based on all corresponding observations (n = 1785). The comparison between vegetation and precipitation trends revealed matched trend directions in 60% of cases, suggesting the potential impact of precipitation dynamics on vegetation covers. Further incorporation of land cover data showed that grassland areas experienced significant dynamics with the highest proportion compared to other vegetation land cover types. Moreover, forest and cropland had the highest positive and negative trend direction proportions. Finally, independent (from trend analysis) sources were used to examine the vegetation dynamics (greening/browning) from other perspectives, confirming Iran’s greening process and agreeing with the trend analysis results. It is believed that the results could support achieving Sustainable Development Goals (SDGs) by serving as an initial stage study for establishing conservation and restoration practices.
Accurate and timely mapping of crop types and having reliable information about the cultivation pattern/area play a key role in various applications, including food security and sustainable ...agriculture management. Remote sensing (RS) has extensively been employed for crop type classification. However, accurate mapping of crop types and extents is still a challenge, especially using traditional machine learning methods. Therefore, in this study, a novel framework based on a deep convolutional neural network (CNN) and a dual attention module (DAM) and using Sentinel-2 time-series datasets was proposed to classify crops. A new DAM was implemented to extract informative deep features by taking advantage of both spectral and spatial characteristics of Sentinel-2 datasets. The spectral and spatial attention modules (AMs) were respectively applied to investigate the behavior of crops during the growing season and their neighborhood properties (e.g., textural characteristics and spatial relation to surrounding crops). The proposed network contained two streams: (1) convolution blocks for deep feature extraction and (2) several DAMs, which were employed after each convolution block. The first stream included three multi-scale residual convolution blocks, where the spectral attention blocks were mainly applied to extract deep spectral features. The second stream was built using four multi-scale convolution blocks with a spatial AM. In this study, over 200,000 samples from six different crop types (i.e., alfalfa, broad bean, wheat, barley, canola, and garden) and three non-crop classes (i.e., built-up, barren, and water) were collected to train and validate the proposed framework. The results demonstrated that the proposed method achieved high overall accuracy and a Kappa coefficient of 98.54% and 0.981, respectively. It also outperformed other state-of-the-art classification methods, including RF, XGBOOST, R-CNN, 2D-CNN, 3D-CNN, and CBAM, indicating its high potential to discriminate different crop types.
Iran is among the driest countries in the world, where many natural hazards, such as floods, frequently occur. This study introduces a straightforward flood hazard assessment approach using remote ...sensing datasets and Geographic Information Systems (GIS) environment in an area located in the western part of Iran. Multiple GIS and remote sensing datasets, including Digital Elevation Model (DEM), slope, rainfall, distance from the main rivers, Topographic Wetness Index (TWI), Land Use/Land Cover (LULC) maps, soil type map, Normalized Difference Vegetation Index (NDVI), and erosion rate were initially produced. Then, all datasets were converted into fuzzy values using a linear fuzzy membership function. Subsequently, the Analytical Hierarchy Process (AHP) technique was applied to determine the weight of each dataset, and the relevant weight values were then multiplied to fuzzy values. Finally, all the processed parameters were integrated using a fuzzy analysis to produce the flood hazard map with five classes of susceptible zones. The bi-temporal Sentinel-1 Synthetic Aperture Radar (SAR) images, acquired before and on the day of the flood event, were used to evaluate the accuracy of the produced flood hazard map. The results indicated that 95.16% of the actual flooded areas were classified as very high and high flood hazard classes, demonstrating the high potential of this approach for flood hazard mapping.
Mapping and monitoring the spatio-temporal variations of the Surface Urban Heat Island (SUHI) and thermal comfort of metropolitan areas are vital to obtaining the necessary information about the ...environmental conditions and promoting sustainable cities. As the most populated city of Iran, Tehran has experienced considerable population growth and Land Cover/Land Use (LULC) changes in the last decades, which resulted in several adverse environmental issues. In this study, 68 Landsat-5 and Landsat-8 images, collected from the Google Earth Engine (GEE), were employed to map and monitor the spatio-temporal variations of LULC, SUHI, and thermal comfort of Tehran between 1989 and 2019. In this regard, planar fitting and Gaussian Surface Model (GSM) approaches were employed to map SUHIs and derive the relevant statistical values. Likewise, the thermal comfort of the city was investigated by the Urban Thermal Field Variance Index (UTFVI). The results indicated that the SUHI intensities have generally increased throughout the city by an average value of about 2.02 °C in the past three decades. The most common reasons for this unfavorable increase were the loss of vegetation cover (i.e., 34.72%) and massive urban expansions (i.e., 53.33%). Additionally, the intra-annual investigations in 2019 revealed that summer and winter, with respectively 8.28 °C and 4.37 °C, had the highest and lowest SUHI magnitudes. Furthermore, the decadal UTFVI maps revealed notable thermal comfort degradation of Tehran, by which in 2019, approximately 52.35% of the city was identified as the region with the worst environmental condition, of which 59.94% was related to human residents. Additionally, the relationships between various air pollutants and SUHI intensities were appraised, suggesting positive relationships (i.e., ranging between 0.23 and 0.43) that can be used for establishing possible two-way mitigations strategies. This study provided analyses of spatio-temporal monitoring of SUHI and UTFVI throughout Tehran that urban managers and policymakers can consider for adaption and sustainable development.
Mangroves are among the most productive ecosystems in existence, with many ecological benefits. Therefore, generating accurate thematic maps from mangrove ecosystems is crucial for protecting, ...conserving, and reforestation planning for these valuable natural resources. In this paper, Sentinel-1 and Sentinel-2 satellite images were used in synergy to produce a detailed mangrove ecosystem map of the Hara protected area, Qeshm, Iran, at 10 m spatial resolution within the Google Earth Engine (GEE) cloud computing platform. In this regard, 86 Sentinel-1 and 41 Sentinel-2 data, acquired in 2019, were employed to generate seasonal optical and synthetic aperture radar (SAR) features. Afterward, seasonal features were inserted into a pixel-based random forest (RF) classifier, resulting in an accurate mangrove ecosystem map with average overall accuracy (OA) and Kappa coefficient (KC) of 93.23% and 0.92, respectively, wherein all classes (except aerial roots) achieved high producer and user accuracies of over 90%. Furthermore, comprehensive quantitative and qualitative assessments were performed to investigate the robustness of the proposed approach, and the accurate and stable results achieved through cross-validation and consistency checks confirmed its robustness and applicability. It was revealed that seasonal features and the integration of multi-source remote sensing data contributed towards obtaining a more reliable mangrove ecosystem map. The proposed approach relies on a straightforward yet effective workflow for mangrove ecosystem mapping, with a high rate of automation that can be easily implemented for frequent and precise mapping in other parts of the world. Overall, the proposed workflow can further improve the conservation and sustainable management of these valuable natural resources.
Precipitation, as an important component of the Earth’s water cycle, plays a determinant role in various socio-economic practices. Consequently, having access to high-quality and reliable ...precipitation datasets is highly demanded. Although Gridded Precipitation Products (GPPs) have been widely employed in different applications, the lack of quantitative assessment of GPPs is a critical concern that should be addressed. This is because the inherent errors in GPPs would propagate into any models in which precipitation values are incorporated, introducing uncertainties into the final results. This paper aims to quantify the capability of six well-known GPPs (TMPA, CHIRPS, PERSIANN, GSMaP, IMERG, and ERA5) at multiple time scales (daily, monthly, and yearly) using in situ observations (over 1.7 million) throughout Iran over the past two decades (2000–2020). Both continuous and categorical metrics were implemented for precipitation intensity and occurrence assessment based on the point-to-pixel comparison approach. Although all metrics did not support the superior performance of any specific GPP, taking all investigations into account, the findings suggested the better performance of the Global Satellite Mapping of Precipitation (GSMaP) in estimating daily precipitation (CC = 0.599, RMSE = 3.48 mm/day, and CSI = 0.454). Based on the obtained continuous metrics, all the GPPs had better performances in dry months, while this did not hold for the categorical metrics. The validation at the station level was also carried out to present the spatial characteristics of errors throughout Iran, indicating higher overestimation/underestimation in regions with higher precipitation rates. The validation analysis over the last two decades illustrated that the GPPs had stable performances, and no improvement was seen, except for the GSMaP, in which its bias error was significantly reduced. The comparisons on monthly and yearly time scales suggested the higher accuracy of monthly and yearly averaged precipitation values than accumulated values. Our study provides valuable guidance to the selection and application of GPPs in Iran and also offers beneficial feedback for further improving these products.
Mangrove ecosystems provide numerous ecological services and serve as vital habitats for a wide range of flora and fauna. Thus, accurate mapping and monitoring of relevant land covers in mangrove ...ecosystems are crucial for effective conservation and management efforts. In this study, we proposed a novel approach for mangrove ecosystem mapping using a Hybrid Selective Kernel-based Convolutional Neural Network (HSK-CNN) framework and multi-temporal Sentinel-2 imagery. A time series of the Normalized Difference Vegetation Index (NDVI) products derived from Sentinel-2 imagery was produced to capture the temporal behavior of land cover types in the dynamic ecosystem of the study area. The proposed algorithm integrated Selective Kernel-based feature extraction techniques to facilitate the effective learning and classification of multiple land cover types within the dynamic mangrove ecosystems. The model demonstrated a high Overall Accuracy (OA) of 94% in classifying eight land cover classes, including mangrove, tidal zone, water, mudflat, urban, and vegetation. The HSK-CNN demonstrated superior performance compared to other algorithms, including random forest (OA = 85%), XGBoost (OA = 87%), Three-Dimensional (3D)-DenseNet (OA = 90%), Two-Dimensional (2D)-CNN (OA = 91%), Multi-Layer Perceptron (MLP)-Mixer (OA = 92%), and Swin Transformer (OA = 93%). Additionally, it was observed that the structure of the network, such as the types of convolutional layers and patch sizes, affected the classification accuracy using the proposed model and, thus, the optimum scenarios and values of these parameters should be determined to obtain the highest possible classification accuracy. Overall, it was observed that the produced map could offer valuable insights into the distribution of different land cover types in the mangrove ecosystem, facilitating informed decision-making for conservation and sustainable management efforts.
When natural disasters occur, timely and accurate building damage assessment maps are vital for disaster management responders to organize their resources efficiently. Pairs of pre- and post-disaster ...remote sensing imagery have been recognized as invaluable data sources that provide useful information for building damage identification. Recently, deep learning-based semantic segmentation models have been widely and successfully applied to remote sensing imagery for building damage assessment tasks. In this study, a two-stage, dual-branch, UNet architecture, with shared weights between two branches, is proposed to address the inaccuracies in building footprint localization and per-building damage level classification. A newly introduced selective kernel module improves the performance of the model by enhancing the extracted features and applying adaptive receptive field variations. The xBD dataset is used to train, validate, and test the proposed model based on widely used evaluation metrics such as F1-score and Intersection over Union (IoU). Overall, the experiments and comparisons demonstrate the superior performance of the proposed model. In addition, the results are further confirmed by evaluating the geographical transferability of the proposed model on a completely unseen dataset from a new region (Bam city earthquake in 2003).
The Baltic Sea is one of the fastest-warming marginal seas globally, and its temperature rise has adversely affected its physical and biochemical characteristics. In this study, forty years ...(1982–2021) of sea surface temperature (SST) data from the advanced very high resolution radiometer (AVHRR) were used to investigate spatial and temporal SST variability of the Baltic Sea. To this end, annual maximum and minimum SST stacked series, i.e., time series of stacked layers of satellite data, were generated using high-quality observations acquired at night and were fed to an automatic algorithm to detect linear and non-linear trend patterns. The linear trend pattern was the dominant trend type in both stacked series, while more pixels with non-linear trend patterns were detected when using the annual minimum SST. However, both stacked series showed increases in SST across the Baltic Sea. Annual maximum SST increased by an average of 0.062 ± 0.041 °C per year between 1982 and 2021, while annual minimum SST increased by an average of 0.035 ± 0.017 °C per year over the same period. Averaging annual maximum and minimum trends produces a spatial average of 0.048 ± 0.022 °C rise in SST per year over the last four decades.