[go: up one dir, main page]

 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (5,585)

Search Parameters:
Keywords = landsat

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
20 pages, 20184 KiB  
Article
Snow Cover Extraction from Landsat 8 OLI Based on Deep Learning with Cross-Scale Edge-Aware and Attention Mechanism
by Zehao Yu, Hanying Gong, Shiqiang Zhang and Wei Wang
Remote Sens. 2024, 16(18), 3430; https://doi.org/10.3390/rs16183430 (registering DOI) - 15 Sep 2024
Viewed by 231
Abstract
Snow cover distribution is of great significance for climate change and water resource management. Current deep learning-based methods for extracting snow cover from remote sensing images face challenges such as insufficient local detail awareness and inadequate utilization of global semantic information. In this [...] Read more.
Snow cover distribution is of great significance for climate change and water resource management. Current deep learning-based methods for extracting snow cover from remote sensing images face challenges such as insufficient local detail awareness and inadequate utilization of global semantic information. In this study, a snow cover extraction algorithm integrating cross-scale edge perception and an attention mechanism on the U-net model architecture is proposed. The cross-scale edge perception module replaces the original jump connection of U-net, enhances the low-level image features by introducing edge detection on the shallow feature scale, and enhances the detail perception via branch separation and fusion features on the deep feature scale. Meanwhile, parallel channel and spatial attention mechanisms are introduced in the model encoding stage to adaptively enhance the model’s attention to key features and improve the efficiency of utilizing global semantic information. The method was evaluated on the publicly available CSWV_S6 optical remote sensing dataset, and the accuracy of 98.14% indicates that the method has significant advantages over existing methods. Snow extraction from Landsat 8 OLI images of the upper reaches of the Irtysh River was achieved with satisfactory accuracy rates of 95.57% (using two, three, and four bands) and 96.65% (using two, three, four, and six bands), indicating its strong potential for automated snow cover extraction over larger areas. Full article
Show Figures

Figure 1

Figure 1
<p>True color CSWV_S6 data synthesized from the red, green, and blue bands (the numbering in the figure corresponds to the original naming in the acquired files).</p>
Full article ">Figure 2
<p>(<b>a</b>) RGB composite of Landsat 8 imagery (red: band 4, green: band 3, blue: band 2). (<b>b</b>) Land cover types.</p>
Full article ">Figure 3
<p>CEFCSAU-net network model architecture. The input size was (512, 512, C), where C denotes the number of channels, and experiments in this paper utilized either 3 or 4; during the model’s operation on a GPU, intermediate feature maps were stored as tensors.</p>
Full article ">Figure 4
<p>Attention mechanism module for channel and space mixing. Here, (H, W, C) represent the height, width, and number of channels of the feature data, respectively, with values determined by input features at different stages. CF and CF’ denote feature maps from various intermediate operations within the channel attention mechanism. Cat Sf, Sf, Sf’ represent feature maps from different intermediate operations of the spatial attention mechanism. The SA feature denotes the feature map post-spatial attention mechanism, the CA feature represents those post-channel attention mechanisms, and the CSA feature illustrates feature maps following the CSA module.</p>
Full article ">Figure 5
<p>Cross-scale edge-aware feature fusion module. Sobelx F, Sobely F, and Laplacian F denote feature maps resulting from various edge detection operations. Shallow F refers to feature maps following shallow feature convolution. Fusion F illustrates feature maps resulting from the fusion of shallow and deep features. Deep F’ represents feature maps after a series of operations on deep features.</p>
Full article ">Figure 6
<p>Snow extraction results of CSWV_S6 data on different segmentation models (a set of two rows is the same image data, and rows two, four, and six are zoomed-in images of local details corresponding to rows one, three, and five. The blue area is snow, the white area is non-snow, and the red area is false detection).</p>
Full article ">Figure 7
<p>Snow extraction results from different deep learning models for Landsat 8 OLI imagery (blue areas are snow, white areas are non-snow, and red areas are false detections).</p>
Full article ">Figure 8
<p>(<b>a</b>) CSWV_S6 test set scores on different models for each type of metrics and (<b>b</b>) Landsat 8 OLI test set scores for various metrics on different models.</p>
Full article ">Figure 9
<p>Score results of the three CSWV_S6 test sets’ example data on the evaluation metrics on each model, with 0.08% of snow image elements in the first row of data, 0.95% of snow image elements in the second row of data, and 1.73% of data in the third row of data.</p>
Full article ">Figure 10
<p>Map of CEFCSAU-net model’s snow extraction in the cloud–snow confusion scenario of CSWV_S6 test set.</p>
Full article ">Figure 11
<p>Heat map comparing the mean values of ablation experiments on the test set with different data sets: (<b>a</b>) CSWV_6 dataset, (<b>b</b>) Landsat8 OLI dataset.</p>
Full article ">Figure 12
<p>(<b>a</b>) Input data image; (<b>b</b>) feature map of the first 8 channels before the intermediate feature data first pass through the CSA module; (<b>c</b>) feature map of the first 8 channels after the feature data first pass through the CSA module.</p>
Full article ">Figure 13
<p>This figure displays average feature maps following skip connections at various stages under different configurations of the CEFCSAU-net model. In this figure, (<b>a1</b>–<b>a4</b>) represent the model configuration without both the CSA and CEF modules; (<b>b1</b>–<b>b4</b>) indicate configurations without the CSA module yet including the CEF module; and (<b>c1</b>–<b>c4</b>) depict configurations featuring both CSA and CEF modules. The dimensions of the four columns of feature maps are sequentially 512 × 512, 256 × 256, 128 × 128, and 64 × 64.</p>
Full article ">
24 pages, 5994 KiB  
Article
Mapping Natural Populus euphratica Forests in the Mainstream of the Tarim River Using Spaceborne Imagery and Google Earth Engine
by Jiawei Zou, Hao Li, Chao Ding, Suhong Liu and Qingdong Shi
Remote Sens. 2024, 16(18), 3429; https://doi.org/10.3390/rs16183429 (registering DOI) - 15 Sep 2024
Viewed by 208
Abstract
Populus euphratica is a unique constructive tree species within riparian desert areas that is essential for maintaining oasis ecosystem stability. The Tarim River Basin contains the most densely distributed population of P. euphratica forests in the world, and obtaining accurate distribution data in [...] Read more.
Populus euphratica is a unique constructive tree species within riparian desert areas that is essential for maintaining oasis ecosystem stability. The Tarim River Basin contains the most densely distributed population of P. euphratica forests in the world, and obtaining accurate distribution data in the mainstream of the Tarim River would provide important support for its protection and restoration. We propose a new method for automatically extracting P. euphratica using Sentinel-1 and 2 and Landsat-8 images based on the Google Earth Engine cloud platform and the random forest algorithm. A mask of the potential distribution area of P. euphratica was created based on prior knowledge to save computational resources. The NDVI (Normalized Difference Vegetation Index) time series was then reconstructed using the preferred filtering method to obtain phenological parameter features, and the random forest model was input by combining the phenological parameter, spectral index, textural, and backscattering features. An active learning method was employed to optimize the model and obtain the best model for extracting P. euphratica. Finally, the map of natural P. euphratica forests with a resolution of 10 m in the mainstream of the Tarim River was obtained. The overall accuracy, producer’s accuracy, user’s accuracy, kappa coefficient, and F1-score of the map were 0.96, 0.98, 0.95, 0.93, and 0.96, respectively. The comparison experiments showed that simultaneously adding backscattering and textural features improved the P. euphratica extraction accuracy, while textural features alone resulted in a poor extraction effect. The method developed in this study fully considered the prior and posteriori information and determined the feature set suitable for the P. euphratica identification task, which can be used to quickly obtain accurate large-area distribution data of P. euphratica. The method can also provide a reference for identifying other typical desert vegetation. Full article
Show Figures

Figure 1

Figure 1
<p>Geographical location of the study area and the distribution of sample points. (<b>a</b>): location of the study area in Xinjiang province in China; (<b>b</b>): training dataset distribution; (<b>c</b>): detailed sample area showing <span class="html-italic">P. euphratica</span> and non–<span class="html-italic">P. euphratica</span> in a Sentinel-2 false-color image.</p>
Full article ">Figure 2
<p>Distribution of validation dataset. The black solid line represents the range of the study area; the red and yellow points represent <span class="html-italic">P. euphratica</span> and non–<span class="html-italic">P. euphratica</span>, respectively.</p>
Full article ">Figure 3
<p>Workflow of the research.</p>
Full article ">Figure 4
<p>Threshold segmentation effect of MNDWI and NDVI. (<b>a</b>): false color image of Jieran Lik Reservoir in Xinjiang Province; (<b>b</b>): statistical result of the corresponding frequency distribution of MNDWI values of water and other ground objects in area (<b>a</b>); (<b>c</b>): false color image of Pazili Tamu in Xinjiang; (<b>d</b>): statistical result for the corresponding frequency distribution of NDVI values of desert bare land and other ground objects in region (<b>c</b>).</p>
Full article ">Figure 5
<p>Comparison of NDVI data before and after spatiotemporal fusion: (<b>a</b>) NDVI data derived from Sentinel-2 before fusion, (<b>b</b>) NDVI data after fusion.</p>
Full article ">Figure 6
<p>Comparison of the effects of different filter functions for: (<b>a</b>) <span class="html-italic">P. euphratica</span>; (<b>b</b>) <span class="html-italic">Tamarix</span>; (<b>c</b>) allee tree; (<b>d</b>) farmland; (<b>e</b>) wetland; (<b>f</b>) urban tree.</p>
Full article ">Figure 7
<p>Comparison between phenological curves of six typical vegetation species. Phenology parameters of (<b>a</b>) <span class="html-italic">P. euphratica</span>, (<b>b</b>) <span class="html-italic">Tamarix</span>, (<b>c</b>) allee tree, (<b>d</b>) farmland, (<b>e</b>) wetland, and (<b>f</b>) urban tree.</p>
Full article ">Figure 8
<p>Importance of different features in the RF classification.</p>
Full article ">Figure 9
<p>Natural <span class="html-italic">P. euphratica</span> forest maps extracted using four feature combinations: (<b>a</b>) PS, (<b>b</b>) PSB, (<b>c</b>) PST, and (<b>d</b>) PSBT.</p>
Full article ">Figure 10
<p>Comparison of <span class="html-italic">P. euphratica</span> extraction results using different feature combinations on Sentinel-2 standard false color images. Rows 1 to 4 show the identification of <span class="html-italic">P. euphratica</span> in desert areas, <span class="html-italic">P. euphratica</span>-dense areas, agricultural areas, and large river areas, respectively. The green area represents the classification result of <span class="html-italic">P. euphratica</span>. The yellow circle corresponding to each row is the area where the extraction results of different feature combinations are quite different.</p>
Full article ">Figure 11
<p>(<b>a</b>) Distribution of natural <span class="html-italic">P. euphratica</span> forest in the mainstream of the Tarim River. (<b>b</b>): UAV image of healthy <span class="html-italic">P. euphratica</span>, (<b>c</b>): classification result of healthy <span class="html-italic">P. euphratica</span>, (<b>d</b>): UAV image of unhealthy <span class="html-italic">P. euphratica</span>, (<b>e</b>): classification result of unhealthy <span class="html-italic">P. euphratica</span>, (<b>f</b>): UAV image of dense <span class="html-italic">P. euphratica</span>, (<b>g</b>): classification result of dense <span class="html-italic">P. euphratica</span>, (<b>h</b>): UAV image of sparse <span class="html-italic">P. euphratica</span>, (<b>i</b>): classification result of sparse <span class="html-italic">P. euphratica</span>. The green area represents the classification results of <span class="html-italic">P. euphratica</span>.</p>
Full article ">Figure 12
<p>Mixed pixel problems associated with <span class="html-italic">P. euphratica</span>: (<b>a</b>) <span class="html-italic">P. euphratica</span> occupying less than one pixel; (<b>b</b>) sandy soil interfering with the reflected signal of <span class="html-italic">P. euphratica</span>. The red box represents a pixel on the images for clearer observation. Basemaps of row 1-2 are UAV images while row 3 are Sentinel-2 standard false color images.</p>
Full article ">
17 pages, 13310 KiB  
Article
Spatiotemporal Dynamics and Drivers of Coastal Wetlands in Tianjin–Hebei over the Past 80 Years
by Feicui Wang, Fu Wang, Ke Zhu, Peng Yang, Tiejun Wang, Yunzhuang Hu and Lijuan Ye
Water 2024, 16(18), 2612; https://doi.org/10.3390/w16182612 (registering DOI) - 14 Sep 2024
Viewed by 349
Abstract
Coastal wetland ecosystems are critical due to their diverse ecological and economic benefits, yet they have been significantly affected by human activities over the past century. Understanding the spatiotemporal changes and underlying factors influencing these ecosystems is crucial for developing effective ecological protection [...] Read more.
Coastal wetland ecosystems are critical due to their diverse ecological and economic benefits, yet they have been significantly affected by human activities over the past century. Understanding the spatiotemporal changes and underlying factors influencing these ecosystems is crucial for developing effective ecological protection and restoration strategies. This study examines the Tianjin–Hebei coastal wetlands using topographic maps from the 1940s and Landsat satellite imagery from 1975, 2000, and 2020, supplemented by historical literature and field surveys. The aim is to analyze the distribution and classification of coastal wetlands across various temporal intervals. The findings indicate an expansion of the Tianjin–Hebei coastal wetlands from 7301.34 km2 in the 1940s to 8041.73 km2 in 2020. However, natural wetlands have declined by approximately 44.36 km2/year, while constructed wetlands have increased by around 53.61 km2/year. The wetlands have also become increasingly fragmented, with higher numbers of patches and densities. The analysis of driving factors points to human activities—such as urban construction, cultivated land reclamation, sea aquaculture, and land reclamation—as the primary contributors to these changes. Furthermore, the study addresses the ecological and environmental issues stemming from wetland changes and proposes strategies for wetland conservation. This research aims to enhance the understanding among researchers and policymakers of the dynamics and drivers of coastal wetland changes, as well as the major challenges in their protection, and to serve as a foundation for developing evidence-based conservation and restoration strategies. Full article
(This article belongs to the Special Issue Soil and Groundwater Quality and Resources Assessment)
Show Figures

Figure 1

Figure 1
<p>Location map of the study area.</p>
Full article ">Figure 2
<p>Distribution of wetlands in different periods.</p>
Full article ">Figure 3
<p>Nature and constructed wetland in different periods (Unit: km<sup>2</sup>).</p>
Full article ">Figure 4
<p>Spatial distribution of the main trajectory codes for wetland changes in the study area.</p>
Full article ">Figure 5
<p>Illustrates the wetland distribution around Tianjin Port.</p>
Full article ">Figure 6
<p>Schematic diagram of coastal wetland restoration locations and projects in the study area.</p>
Full article ">
20 pages, 7101 KiB  
Article
Tracking Evapotranspiration Patterns on the Yinchuan Plain with Multispectral Remote Sensing
by Junzhen Meng, Xiaoquan Yang, Zhiping Li, Guizhang Zhao, Peipei He, Yabing Xuan and Yunfei Wang
Sustainability 2024, 16(18), 8025; https://doi.org/10.3390/su16188025 - 13 Sep 2024
Viewed by 353
Abstract
Evapotranspiration (ET) is a critical component of the hydrological cycle, and it has a decisive impact on the ecosystem balance in arid and semi-arid regions. The Yinchuan Plain, located in the Gobi of Northwest China, has a strong surface ET, which has a [...] Read more.
Evapotranspiration (ET) is a critical component of the hydrological cycle, and it has a decisive impact on the ecosystem balance in arid and semi-arid regions. The Yinchuan Plain, located in the Gobi of Northwest China, has a strong surface ET, which has a significant impact on the regional water resource cycle. However, there is a current lack of high-resolution evapotranspiration datasets and a substantial amount of time is required for long-time series remote sensing evapotranspiration estimation. In order to assess the ET pattern in this region, we obtained the actual ET (ETa) of the Yinchuan Plain between 1987 and 2020 using the Google Earth Engine (GEE) platform. Specifically, we used Landsat TM+/OLI remote sensing imagery and the GEE Surface Energy Balance Model (geeSEBAL) to analyze the spatial distribution pattern of ET over different seasons. We then reproduced the interannual variation in ET from 1987 to 2020, and statistically analyzed the distribution patterns and contributions of ET with regard to different land use types. The results show that (1) the daily ETa of the Yinchuan Plain is the highest in the central lake wetland area in spring, with a maximum value of 4.32 mm day−1; in summer, it is concentrated around the croplands and water bodies, with a maximum value of 6.90 mm day−1; in autumn and winter, it is mainly concentrated around the water bodies and impervious areas, with maximum values of 3.93 and 1.56 mm day−1, respectively. (2) From 1987 to 2020, the ET of the Yinchuan Plain showed an obvious upward and downward trend in some areas with significant land use changes, but the overall ET of the region remained relatively stable without dramatic fluctuations. (3) The ETa values for different land use types in the Yinchuan Plain region are ranked as follows: water body > cultivated land > impervious > grassland > bare land. Our results showed that geeSEBAL is highly applicable in the Yinchuan Plain area. It allows for the accurate and detailed inversion of ET and has great potential for evaluating long-term ET in data-scarce areas due to its low meteorological sensitivity, which facilitates the study of the regional hydrological cycle and water governance. Full article
(This article belongs to the Section Sustainable Water Management)
Show Figures

Figure 1

Figure 1
<p>Location of the meteorological stations and Landsat images that were used to illustrate the land cover conditions in the Yinchuan Plain area.</p>
Full article ">Figure 2
<p>Comparison between meteorological stations with large ET and geeSEBAL ET. Diamond-shaped points represent outliers lying outside the 150% inter-quartile range.</p>
Full article ">Figure 3
<p>Comparison between ET<sub>p</sub> and geeSEBAL ET (<b>a</b>–<b>d</b>). Compared with large-scale evapotranspiration, R<sup>2</sup> has significantly improved, indicating that the model correlation is influenced by external factors.</p>
Full article ">Figure 4
<p>Comparison between small ET and geeSEBAL ET ((<b>a</b>–<b>f</b>) presents different meteorological stations in the Yinchuan Plain).</p>
Full article ">Figure 5
<p>Seasonal ET<sub>a</sub> changes in the Yinchuan Plain. (<b>a</b>) Spring ET distribution; (<b>b</b>) Summer ET distribution; (<b>c</b>) Autum ET distribution; (<b>d</b>) Winter ET distribution.</p>
Full article ">Figure 6
<p>Trends in ET<sub>a</sub> on the Yinchuan Plain from 1987 to 2020.</p>
Full article ">Figure 7
<p>Area of Yinchuan Plain land use types.</p>
Full article ">Figure 8
<p>ET<sub>a</sub> in different subsurface types.</p>
Full article ">Figure 9
<p>Comparison of remote sensing imagery (<b>a</b>), land use classification (<b>b</b>), and ET imagery (<b>c</b>) on the Yinchuan Plain.</p>
Full article ">Figure 10
<p>Impervious areas misclassified in some intersecting land types. water bodies are identified as impervious areas (red color).</p>
Full article ">Figure 11
<p>ET<sub>a</sub> contribution of different subsurface types.</p>
Full article ">Figure 12
<p>geeSEBAL with the batch image estimation mode.</p>
Full article ">
28 pages, 20281 KiB  
Article
Spatiotemporal Prediction of Conflict Fatality Risk Using Convolutional Neural Networks and Satellite Imagery
by Seth Goodman, Ariel BenYishay and Daniel Runfola
Remote Sens. 2024, 16(18), 3411; https://doi.org/10.3390/rs16183411 - 13 Sep 2024
Viewed by 309
Abstract
As both satellite imagery and image-based machine learning methods continue to improve and become more accessible, they are being utilized in an increasing number of sectors and applications. Recent applications using convolutional neural networks (CNNs) and satellite imagery include estimating socioeconomic and development [...] Read more.
As both satellite imagery and image-based machine learning methods continue to improve and become more accessible, they are being utilized in an increasing number of sectors and applications. Recent applications using convolutional neural networks (CNNs) and satellite imagery include estimating socioeconomic and development indicators such as poverty, road quality, and conflict. This article builds on existing work leveraging satellite imagery and machine learning for estimation or prediction, to explore the potential to extend these methods temporally. Using Landsat 8 imagery and data from the Armed Conflict Location & Event Data Project (ACLED) we produce subnational predictions of the risk of conflict fatalities in Nigeria during 2015, 2017, and 2019 using distinct models trained on both yearly and six-month windows of data from the preceding year. We find that predictions at conflict sites leveraging imagery from the preceding year for training can predict conflict fatalities in the following year with an area under the receiver operating characteristic curve (AUC) of over 75% on average. While models consistently outperform a baseline comparison, and performance in individual periods can be strong (AUC > 80%), changes based on ground conditions such as the geographic scope of conflict can degrade performance in subsequent periods. In addition, we find that training models using an entire year of data slightly outperform models using only six months of data. Overall, the findings suggest CNN-based methods are moderately effective at detecting features in Landsat satellite imagery associated with the risk of fatalities from conflict events across time periods. Full article
(This article belongs to the Special Issue Weakly Supervised Deep Learning in Exploiting Remote Sensing Big Data)
Show Figures

Figure 1

Figure 1
<p>Overview of the methodology applied for training and validating conflict fatality models at each discrete time period.</p>
Full article ">Figure 2
<p>Accuracy distribution across all tests based on time period (inset year range) and temporal window (x-axis). Boxplot outline colors reflect the same temporal window across the different time periods. The blue baseline is the average baseline accuracy across all time periods and temporal windows.</p>
Full article ">Figure 3
<p>ROC curve produced from CNNs trained on 2014 imagery (January–December) to predict the probability of a fatality if there is a conflict event in 2015 (January–December).</p>
Full article ">Figure 4
<p>Location of conflict events [<a href="#B3-remotesensing-16-03411" class="html-bibr">3</a>] and prediction results in 2015 (<b>left</b>), 2017 (<b>center</b>), 2019 (<b>right</b>). Imagery ©2024 TerraMetrics, Map data ©2024 Google.</p>
Full article ">Figure A1
<p>ROC curve produced from model train on 2014 (January–December) imagery to predict the probability of a fatality if there is a conflict event in 2015 (January–December).</p>
Full article ">Figure A2
<p>ROC curve produced from model train on 2015 (January–December) imagery to predict the probability of a fatality if there is a conflict event in 2016 (January–December).</p>
Full article ">Figure A3
<p>ROC curve produced from model train on 2018 (January–December) imagery to predict the probability of a fatality if there is a conflict event in 2019 (January–December).</p>
Full article ">Figure A4
<p>ROC curve produced from model train on 2014 h1 (January–June) imagery to predict the probability of a fatality if there is a conflict event in 2015 h1 (January–June).</p>
Full article ">Figure A5
<p>ROC curve produced from model train on 2014 h2 (July–December) imagery to predict the probability of a fatality if there is a conflict event in 2015 h2 (July–December).</p>
Full article ">Figure A6
<p>ROC curve produced from model train on 2016 h1 (January–June) imagery to predict the probability of a fatality if there is a conflict event in 2017 h1 (January–June).</p>
Full article ">Figure A7
<p>ROC curve produced from model train on 2016 h2 (July–December) imagery to predict the probability of a fatality if there is a conflict event in 2017 h2 (July–December).</p>
Full article ">Figure A8
<p>ROC curve produced from model train on 2018 h1 (January–June) imagery to predict the probability of a fatality if there is a conflict event in 2019 h1 (January–June).</p>
Full article ">Figure A9
<p>ROC curve produced from model train on 2018 h2 (July–December) imagery to predict the probability of a fatality if there is a conflict event in 2019 h2 (July–December).</p>
Full article ">
14 pages, 2411 KiB  
Article
Enhanced Water Quality Inversion in the Ningxia Yellow River Basin Using a Hybrid PCWA-ResCNN Model: Insights from Landsat-8 Data
by Qi Li, Zhonghua Guo, Jialong Li, Xiaojun Li and Bo Ban
Appl. Sci. 2024, 14(18), 8264; https://doi.org/10.3390/app14188264 - 13 Sep 2024
Viewed by 367
Abstract
The real-time monitoring and evaluation of water quality provides a scientific basis for water resource management and promotes regional sustainable development. This study established a database using Landsat-8 satellite data and water quality data from the Ningxia Yellow River basin in China, spanning [...] Read more.
The real-time monitoring and evaluation of water quality provides a scientific basis for water resource management and promotes regional sustainable development. This study established a database using Landsat-8 satellite data and water quality data from the Ningxia Yellow River basin in China, spanning 2021 to 2023, and this paper proposes a custom residual convolutional neural network model with a hybrid attention mechanism, referred to as PCWA-ResCNN. The accuracy of the model in predicting turbidity, permanganate, ammonia nitrogen, and dissolved oxygen concentration was more than 95%. Compared to convolutional neural networks and long short-term memory models, this model performed better in predicting water quality parameters with significantly improved prediction performance. In terms of spatial distribution, the pollution degree in the middle reaches of the basin is relatively serious. However, the overall water quality is good, being mainly Class I and Class II water quality. The hybrid model established in this paper can better capture the complex nonlinear relationship between the observed values and the surface water reflectance, showing strong robustness. This model can be used for the water quality monitoring of complex inland rivers and lakes, and it can also provide effective support for relevant government departments to formulate scientific and reasonable water quality management policies. Full article
(This article belongs to the Section Ecology Science and Engineering)
Show Figures

Figure 1

Figure 1
<p>Map of sampling sites in the Ningxia Yellow River basin, China.</p>
Full article ">Figure 2
<p>PCWA-ResCNN model structure.</p>
Full article ">Figure 3
<p>Pearson correlation between water quality parameters and each band.</p>
Full article ">Figure 4
<p>Scatter diagram of prediction using the water quality inversion model.</p>
Full article ">Figure 5
<p>Spatial variation characteristics of water quality parameters.</p>
Full article ">
15 pages, 4826 KiB  
Article
Assessing Evapotranspiration Changes in Response to Cropland Expansion in Tropical Climates
by Leonardo Laipelt, Julia Brusso Rossi, Bruno Comini de Andrade, Morris Scherer-Warren and Anderson Ruhoff
Remote Sens. 2024, 16(18), 3404; https://doi.org/10.3390/rs16183404 - 13 Sep 2024
Viewed by 204
Abstract
The expansion of cropland in tropical regions has significantly accelerated in recent decades, triggering an escalation in water demand and changing the total water loss to the atmosphere (evapotranspiration). Additionally, the increase in areas dedicated to agriculture in tropical climates coincides with an [...] Read more.
The expansion of cropland in tropical regions has significantly accelerated in recent decades, triggering an escalation in water demand and changing the total water loss to the atmosphere (evapotranspiration). Additionally, the increase in areas dedicated to agriculture in tropical climates coincides with an increased frequency of drought events, leading to a series of conflicts among water users. However, detailed studies on the impacts of changes in water use due to agriculture expansion, including irrigation, are still lacking. Furthermore, the higher presence of clouds in tropical environments poses challenges for the availability of high-resolution data for vegetation monitoring via satellite images. This study aims to analyze 37 years of agricultural expansion using the Landsat collection and a satellite-based model (geeSEBAL) to assess changes in evapotranspiration resulting from cropland expansion in tropical climates, focusing on the São Marcos River Basin in Brazil. It also used a methodology for estimating daily evapotranspiration on days without satellite images. The results showed a 34% increase in evapotranspiration from rainfed areas, mainly driven by soybean cultivation. In addition, irrigated areas increased their water use, despite not significantly changing water use at the basin scale. Conversely, natural vegetation areas decreased their evapotranspiration rates by 22%, suggesting possible further implications with advancing changes in land use and land cover. Thus, this study underscores the importance of using satellite-based evapotranspiration estimates to enhance our understanding of water use across different land use types and scales, thereby improving water management strategies on a large scale. Full article
Show Figures

Figure 1

Figure 1
<p>São Marcos River Basin: location in Brazil (<b>a</b>), climate zones according to Köppen–Geiger classification and irrigation pivots (<b>b</b>), and Normalized Difference Vegetation Index (NDVI) values computed using average composition of Landsat 8 for 2021 (<b>c</b>).</p>
Full article ">Figure 2
<p>Daily average of <span class="html-italic">ET</span> illustrated as boxplot for each land cover and use (<b>a</b>) for Landsat scenes between 1986 and 2022. We also illustrated the seasonal monthly average of <span class="html-italic">ET</span> (<b>b</b>), and trends of annual average <span class="html-italic">ET</span> for different land types, with natural vegetation (forest and savanna) demonstrating positive trends over the years, as well as irrigated areas (<b>c</b>).</p>
Full article ">Figure 3
<p>Changes in the <span class="html-italic">ET</span> spatial patterns for the São Marcos River Basin from 1986 to 2022 (<b>a</b>). The contribution of the water usage for each land cover and use between 1986 and 2022 is shown in (<b>b</b>), whereas (<b>c</b>) illustrates changes in land cover and use.</p>
Full article ">Figure 4
<p>Annual composition ET (mm day<sup>−1</sup>) in the São Marcus River Basin between 1986 (<b>a</b>) and 2021 (<b>b</b>). Highlighted plots showed the expressive number of pivot irrigation systems over the basin for specific locations.</p>
Full article ">Figure 5
<p>Monthly <span class="html-italic">ET</span> in the São Marcus River basin was analyzed for each month of one water year (2019 and 2020). During the dry season (May to September), precipitation is limited and radiation availability is high, being a water-limited environment. Consequently, lower <span class="html-italic">ET</span> values are observed during the dry season, while the wet season increases <span class="html-italic">ET</span> rates due to higher precipitation availability.</p>
Full article ">Figure 6
<p>Seasonal differences in daily <span class="html-italic">ET</span> for irrigated and rainfed croplands in the São Marcus River Basin (<b>a</b>), and the difference between both estimations (<b>b</b>). We used a simplified method to fill the gap between Landsat scenes by interpolating <span class="html-italic">EF</span> over time and multiplying with the respective reference <span class="html-italic">ET</span>.</p>
Full article ">
35 pages, 6364 KiB  
Article
Mapping the Influence of Olympic Games’ Urban Planning on the Land Surface Temperatures: An Estimation Using Landsat Series and Google Earth Engine
by Joan-Cristian Padró, Valerio Della Sala, Marc Castelló-Bueno and Rafael Vicente-Salar
Remote Sens. 2024, 16(18), 3405; https://doi.org/10.3390/rs16183405 - 13 Sep 2024
Viewed by 397
Abstract
The Olympic Games are a sporting event and a catalyst for urban development in their host city. In this study, we utilized remote sensing and GIS techniques to examine the impact of the Olympic infrastructure on the surface temperature of urban areas. Using [...] Read more.
The Olympic Games are a sporting event and a catalyst for urban development in their host city. In this study, we utilized remote sensing and GIS techniques to examine the impact of the Olympic infrastructure on the surface temperature of urban areas. Using Landsat Series Collection 2 Tier 1 Level 2 data and cloud computing provided by Google Earth Engine (GEE), this study examines the effects of various forms of Olympic Games facility urban planning in different historical moments and location typologies, as follows: monocentric, polycentric, peripheric and clustered Olympic ring. The GEE code applies to the Olympic Games that occurred from Paris 2024 to Montreal 1976. However, this paper focuses specifically on the representative cases of Paris 2024, Tokyo 2020, Rio 2016, Beijing 2008, Sydney 2000, Barcelona 1992, Seoul 1988, and Montreal 1976. The study is not only concerned with obtaining absolute land surface temperatures (LST), but rather the relative influence of mega-event infrastructures on mitigating or increasing the urban heat. As such, the locally normalized land surface temperature (NLST) was utilized for this purpose. In some cities (Paris, Tokyo, Beijing, and Barcelona), it has been determined that Olympic planning has resulted in the development of green spaces, creating “green spots” that contribute to lower-than-average temperatures. However, it should be noted that there is a significant variation in temperature within intensely built-up areas, such as Olympic villages and the surrounding areas of the Olympic stadium, which can become “hotspots.” Therefore, it is important to acknowledge that different planning typologies of Olympic infrastructure can have varying impacts on city heat islands, with the polycentric and clustered Olympic ring typologies displaying a mitigating effect. This research contributes to a cloud computing method that can be updated for future Olympic Games or adapted for other mega-events and utilizes a widely available remote sensing data source to study a specific urban planning context. Full article
(This article belongs to the Special Issue Urban Planning Supported by Remote Sensing Technology II)
Show Figures

Figure 1

Figure 1
<p>Location of the Olympic Game cities from 1972 to 2024, which are included in the Google Earth Engine code. The cities used as examples in this paper, representing four Olympic urban planning patterns, are highlighted in yellow. Source: Author’s own elaboration based on data from Open Street Map (@OpenStreetMap contributors) and International Olympic Committee (IOC) information [<a href="#B59-remotesensing-16-03405" class="html-bibr">59</a>].</p>
Full article ">Figure 2
<p>Area of interest (city AOI) of the eight cities analysed (red outline), and its corresponding area of interest (Olympic AOI) of the Olympic facilities (yellow outline). (<b>a</b>) In the Paris case, the city AOI is defined by the Ille de France administrative boundaries. (<b>b</b>) In the Tokyo case, the city AOI is defined by some municipalities of the Tokyo Metropolitan Area administrative boundaries. (<b>c</b>) In the Rio case, the city AOI is defined by the Rio de Janeiro Municipality administrative boundaries. (<b>d</b>) In the Beijing case, the city AOI is defined by Beijing’s central urban area. (<b>e</b>) In the Sydney case, the city AOI is defined by some municipalities of the North South Wales administrative boundaries. (<b>f</b>) In the Barcelona case, the city AOI is defined by the administrative boundaries of Barcelonès. (<b>g</b>) In the Seoul case, the city AOI is defined by Keijo Teukbyeolsi administrative boundaries. (<b>h</b>) In the Montreal case, the city AOI is defined by the Champlain, Communauté Urbaine de Montréal and Laval administrative boundaries.</p>
Full article ">Figure 3
<p>Overall methodology and processing chain.</p>
Full article ">Figure 4
<p>Area of interest (city AOI) of the eight cities analysed (red outline), and the corresponding areas of interest (Olympic AOI) of the Olympic facilities (yellow outline). (<b>a</b>) In the Paris case, the city AOI is defined by the Ille de France administrative boundaries and the Olympic AOI is clustered. (<b>b</b>) In the Tokyo case, the city AOI is defined by some municipalities of the Tokyo Metropolitan Area administrative boundaries and the Olympic AOI is polycentric. (<b>c</b>) In the Rio case, the city AOI is defined by the Rio de Janeiro Municipality administrative boundaries and the Olympic AOI is peripheric. (<b>d</b>) In the Beijing case, the city AOI is defined by Beijing’s central urban area and the Olympic AOI is polycentric. (<b>e</b>) In the Sydney case, the city AOI is defined by some municipalities of North South Wales administrative boundaries and the Olympic AOI is peripheric. (<b>f</b>) In the Barcelona case, the city AOI is defined by Barcelonès administrative boundaries and the Olympic AOI is clustered. (<b>g</b>) In the Seoul case, the city AOI is defined by Keijo Teukbyeolsi administrative boundaries and the Olympic AOI is monocentric. (<b>h</b>) In the Montreal case, the city AOI is defined by the Champlain, Communauté Urbaine de Montréal and Laval administrative boundaries and the Olympic AOI is monocentric.</p>
Full article ">Figure 5
<p>Normalized difference vegetation index (NDVI) maps created using the median synthetic image over a 5-year period, for each of the 8 cities analysed in this study. NDVI was calculated using the NIR and red bands of the synthetic image (see Equation (3)). Additionally, there is a focus on the main Olympic facilities. The real data range is [−1 to 1] but this was stretched to [−0.25 to 0.25] in all of the maps for better understanding and visualization. (<b>a</b>) In the Paris case, higher NDVI levels are in the periphery, but some Olympic facilities (i.e. Champs-de-Mars) take advantage of inner green areas. (<b>b</b>) In the Tokyo case, there are sparse but important green spaces in the central area. (<b>c</b>) In the Rio case, there are elevated NDVI levels for the entire urban area, but not in the main Olympic facilities. (<b>d</b>) In the Beijing case, the central urban area has sparse green spaces, with overall moderate NDVI levels, where Olympic facilities where placed. (<b>e</b>) In the Sydney case, elevated NDVI levels suggest that their urban area has a considerable amount of green space, including some parts pf the Olympic Park. (<b>f</b>) In the Barcelona case, sparse green spaces can be found, highlighting Montjuïc Olympic ring area. (<b>g</b>) In the Seoul case, higher NDVI levels are located in the north and the south periphery, not where the Olympic Park was placed. (<b>h</b>) In the Montreal case, the overall urban area presents high NDVI levels, and some of the Olympic Park area was also located in a green area.</p>
Full article ">Figure 6
<p>Normalized difference built-up index (NDBI) maps created using the median synthetic image over a 5-year period, for each of the 8 cities analysed in this study. NDBI was calculated using the NIR and SWIR1 bands of the synthetic image (see Equation (4)). Additionally, there is a focus on the main Olympic facilities. The real data range is [−1 to 1] but this was stretched to [−0.25 to 0.25] in all the maps for better understanding and visualization. (<b>a</b>) In the Paris case, radial urban configuration shows a dense urbanised centre with high NDBI levels, such as the Stade de France. (<b>b</b>) In the Tokyo case, the area is densely urbanised, such as the Tokyo Dome complex, but with interstitial green spaces. (<b>c</b>) In the Rio case, there are several densely urbanised focuses with high NDBI levels, such as the Barra Olímpica complex, limited by densely vegetated areas. (<b>d</b>) In the Beijing case, the concentric pattern leads to a densely urbanised city with high NDBI levels, such as the National Stadium. (<b>e</b>) In the Sydney case, the extensive urbanization sprawl is combined with green spaces, such as the Olympic Park. (<b>f</b>) In the Barcelona case, the gridded configuration shows urban continuity and density with very high NDBI levels, such as the Olympic Village. (<b>g</b>) In the Seoul case, the urbanisation is dense around the Han River, with generalized high NDBI levels, such as the Jasmil Sports Complex. (<b>h</b>) In the Montreal case, the gridded pattern presents dense build-up areas combined with inner green spaces, such as the Olympic Park and the adjacent Botanical Garden.</p>
Full article ">Figure 7
<p>Normalized land surface temperature (NLST) maps created using the median synthetic image over a 5-year period, for each of the 8 cities analysed in this study. NLST was calculated using the thermal band of the synthetic image and the minimum and maximum LST values in each city (see Equation (4)). Additionally, there was a focus on the main Olympic facilities. The transect used to calculate the thermal profile of the NLST in each city is included. (<b>a</b>) In the Paris case, the relative high temperature focuses are in the central, north and south areas, also some Olympic facilities such as the Stade de France. (<b>b</b>) In the Tokyo case, the relative high temperature focuses are on the port and around the centre of the SUHI, while Olympic facilities are relative lower temperature zones. (<b>c</b>) In the Rio case, the relative low temperature focuses of the forested areas can be seen in the centre of the AOI, and Olympic venues are acting as relative hotspots. (<b>d</b>) In the Beijing case, the relative high temperature focuses are on the south-west, south and east areas, locating the Olympic facilities in relative lower temperature areas. (<b>e</b>) In the Sydney case, the Olympic venues act as hotspots in relation with the surroundings. (<b>f</b>) In the Barcelona case, the Olympic ring is a relative green spot in comparison with the urban area. (<b>g</b>) In the Seoul case, the Olympic facilities are acting as relative high temperature areas in the SUHI. (<b>h</b>) In the Montreal case, the location of the Stade Olympique is part of the relative higher temperature areas.</p>
Full article ">Figure 8
<p>Transects created using the NSLT, the NDVI and the NDBI images. A segment was digitized that crossed the city and the Olympic facilities to obtain the thermal profile from the NLST. Additionally, to compare results, the NDVI and the NDBI profiles were added. This was undertaken with the Profile Tool v.4.2.6 QGIS plugin, which essentially intersects the segment with the target raster, and extracts the value of the overlapped pixels. The result is a table with values that can be plotted in the GIS v.3.32 software or exported to another software to edit the graph. (<b>a</b>) In the Paris case, there is a peak in the Stade de France NLST transect graph, indicating a hotspot in this location in relation to the Paris UHI. (<b>b</b>) In the Tokyo case, the thermal peak is located over the Stadium and the Dome. (<b>c</b>) In the Rio case, the thermal peak is located over Barra Olímpica and a secondary peak is found over Maracaná. (<b>d</b>) In the Beijing case, the hottest location is the Beijing National Stadium. (<b>e</b>) In the Sydney case, the extensive and low-density neighbourhoods, with many green spaces, contrasts with the Olympic Stadium and the central and dense downtown, where the thermal peaks are located. (<b>f</b>) In the Barcelona case, the highest surface temperature is in the industrial area, and the Olympic Ring has low relative temperatures due to its vegetated park areas. (<b>g</b>) In the Seoul case, the Han River presents the lowest relative surface temperatures, with the higher temperatures located on dense urban areas and over the Olympic Stadium. (<b>h</b>) In the Montreal case, the higher relative surface temperatures are found on the dense residential areas, and there is observed a peak just over the Olympic Park.</p>
Full article ">Figure 9
<p>Boxplots relating the normalized land surface temperature (NLST) within each urban area and within its Olympic facilities. (<b>a</b>) In the Paris case, the boxplot indicates that the Olympic facilities contribute to a slight increase in the relative LST in Paris’s urban area. (<b>b</b>) In the Tokyo case, the boxplot indicates a strong contribution of the Olympic facilities to reducing the overall LST in the Tokyo urban area. (<b>c</b>) In the Rio case, boxplot indicates that the Olympic facilities contribute to increase the overall SUHI LST. (<b>d</b>) In the Beijing case, the boxplot shows a median and average LST lower in the Olympic facilities, thus a strong contribution to the reduced overall LST in the Beijing urban area. (<b>e</b>) In the Sydney case, the NLST median and average values within the Olympic area are much higher than in the overall urban area; thus, the Olympic facilities contribute to the overall increase in LST of the resulting urban area after the games. (<b>f</b>) In the Barcelona case, the median and average LST is lower in the Olympic facilities, which significantly contributes to the overall reduction of LST in the urban area of Barcelona. (<b>g</b>) In the Seoul case, the boxplot suggest that the Olympic facilities have led to a relative rise of LST in the Seoul urban area. (<b>h</b>) In the Montreal case, the average and median values, as well as the higher position of the 1st and 3rd quartile, suggest that the Olympic venues have contributed to an overall relative increase of the LST in the resulting Montreal urban area after the games.</p>
Full article ">Figure 10
<p>Linear simple regressions relating the LST and the NLST pixels overlapped by the thermal transect defined in each city. The LST(K) and the NLST are expected to perfectly correlate in a simple linear regression because they are simply converted using the scaling method in all the cases (<b>a</b>–<b>h</b>) (see <a href="#sec2dot3dot1-remotesensing-16-03405" class="html-sec">Section 2.3.1</a>).</p>
Full article ">
19 pages, 6418 KiB  
Article
Evaluating Sugarcane Yield Estimation in Thailand Using Multi-Temporal Sentinel-2 and Landsat Data Together with Machine-Learning Algorithms
by Jaturong Som-ard, Savittri Ratanopad Suwanlee, Dusadee Pinasu, Surasak Keawsomsee, Kemin Kasa, Nattawut Seesanhao, Sarawut Ninsawat, Enrico Borgogno-Mondino and Filippo Sarvia
Land 2024, 13(9), 1481; https://doi.org/10.3390/land13091481 - 13 Sep 2024
Viewed by 433
Abstract
Updated and accurate crop yield maps play a key role in the agricultural environment. Their application enables the support for sustainable agricultural practices and the formulation of effective strategies to mitigate the impacts of climate change. Farmers can apply the maps to gain [...] Read more.
Updated and accurate crop yield maps play a key role in the agricultural environment. Their application enables the support for sustainable agricultural practices and the formulation of effective strategies to mitigate the impacts of climate change. Farmers can apply the maps to gain an overview of the yield variability, improving farm management practices and optimizing inputs to increase productivity and sustainability such as fertilizers. Earth observation (EO) data make it possible to map crop yield estimations over large areas, although this will remain challenging for specific crops such as sugarcane. Yield data collection is an expensive and time-consuming practice that often limits the number of samples collected. In this study, the sugarcane yield estimation based on a small number of training datasets within smallholder crop systems in the Tha Khan Tho District, Thailand for the year 2022 was assessed. Specifically, multi-temporal satellite datasets from multiple sensors, including Sentinel-2 and Landsat 8/9, were involved. Moreover, in order to generate the sugarcane yield estimation maps, only 75 sampling plots were selected and surveyed to provide training and validation data for several powerful machine-learning algorithms, including multiple linear regression (MLR), stepwise multiple regression (SMR), partial least squares regression (PLS), random forest regression (RFR), and support vector regression (SVR). Among these algorithms, the RFR model demonstrated outstanding performance, yielding an excellent result compared to existing techniques, achieving an R-squared (R2) value of 0.79 and a root mean square error (RMSE) of 3.93 t/ha (per 10 m × 10 m pixel). Furthermore, the mapped yields across the region closely aligned with the official statistical data from the Office of the Cane and Sugar Board (with a range value of 36,000 ton). Finally, the sugarcane yield estimation model was applied to over 2100 sugarcane fields in order to provide an overview of the current state of the yield and total production in the area. In this work, the different yield rates at the field level were highlighted, providing a powerful workflow for mapping sugarcane yields across large regions, supporting sugarcane crop management and facilitating decision-making processes. Full article
Show Figures

Figure 1

Figure 1
<p>Flow chart of the implemented methodology for mapping sugarcane yield in 2022 at Tha Khan Tho District, Thailand using multi-temporal Sentinel-2 (S2) and Landsat 8/9 (L8/9) dataset together with the several machine-learning methods.</p>
Full article ">Figure 2
<p>Study area (<b>a</b>): background shows Sentinel-2 (S2) imagery (image composites: during November 2022) with false color (Red = band 8: Green = band 4: Blue = band 3). The 60 yellow sampling plots are used for training datasets and remaining 15 blue plots were used for validating the mapped results. (<b>b</b>) is a location of the Tha Khan Tho District, Kalasin Province, Thailand (study region), (<b>c</b>) shows sampling plot with size of 10 m × 10 m.</p>
Full article ">Figure 3
<p>The sugarcane field dataset (2364 fields) was visually interpreted using very high-resolution imagery as Planet imagery during November 2022.</p>
Full article ">Figure 4
<p>The ranking importance of the features using the random forest (RF) method for the year 2022 with sampling plots and the multi-temporal Sentinel-2 (S2) and Landsat 8/Landsat 9 (L8/9) data.</p>
Full article ">Figure 5
<p>Zoom box of the estimated yield maps with yield value from 59 to 108 based on productive models: multiple linear regression (MLR) (<b>a</b>); stepwise multiple regression (SMR) (<b>b</b>); partial least squares regression (PLS) (<b>c</b>); random forest regression (RFR) (<b>d</b>); and support vector regression (SVR) (<b>e</b>). The sugarcane fields the entire study area (<b>f</b>) are shown.</p>
Full article ">Figure 6
<p>The scatter plots of sugarcane yield estimation results using five productive models together with multi-temporal Sentinel-2 (S2) and Landsat 8/9 (L8/9) data: multiple linear regression (MLR) (<b>a</b>); stepwise multiple regression (SMR) (<b>b</b>); partial least squares regression (PLS) (<b>c</b>); random forest regression (RFR) (<b>d</b>); and support vector regression (SVR) (<b>e</b>).</p>
Full article ">Figure 7
<p>A comparison of observed yield (t/ha) and estimated yield (t/ha) of 15 sampling fields across the study area based on the best random forest regression (RFR)-predictive model together with Sentinel-2 (S2) and Landsat data.</p>
Full article ">Figure 8
<p>The spatial distribution of estimated yield in 2022 using the best random forest regression (RFR) together with Sentinel-2 (S2) and Landsat data.</p>
Full article ">Figure 9
<p>Histogram of the frequency distribution of estimated yield (t/ha) across the Tha Khan Tho District, Thailand, from 10 m × 10 m Sentinel-2 (S2). The red dotted line is the mean value of the estimated yield in this region.</p>
Full article ">
20 pages, 11776 KiB  
Article
Leveraging Machine Learning and Remote Sensing for Water Quality Analysis in Lake Ranco, Southern Chile
by Lien Rodríguez-López, Lisandra Bravo Alvarez, Iongel Duran-Llacer, David E. Ruíz-Guirola, Samuel Montejo-Sánchez, Rebeca Martínez-Retureta, Ernesto López-Morales, Luc Bourrel, Frédéric Frappart and Roberto Urrutia
Remote Sens. 2024, 16(18), 3401; https://doi.org/10.3390/rs16183401 - 13 Sep 2024
Viewed by 490
Abstract
This study examines the dynamics of limnological parameters of a South American lake located in southern Chile with the objective of predicting chlorophyll-a levels, which are a key indicator of algal biomass and water quality, by integrating combined remote sensing and machine learning [...] Read more.
This study examines the dynamics of limnological parameters of a South American lake located in southern Chile with the objective of predicting chlorophyll-a levels, which are a key indicator of algal biomass and water quality, by integrating combined remote sensing and machine learning techniques. Employing four advanced machine learning models (recurrent neural network (RNNs), long short-term memory (LSTM), recurrent gate unit (GRU), and temporal convolutional network (TCNs)), the research focuses on the estimation of chlorophyll-a concentrations at three sampling stations within Lake Ranco. The data span from 1987 to 2020 and are used in three different cases: using only in situ data (Case 1), using in situ and meteorological data (Case 2), using in situ, and meteorological and satellite data from Landsat and Sentinel missions (Case 3). In all cases, each machine learning model shows robust performance, with promising results in predicting chlorophyll-a concentrations. Among these models, LSTM stands out as the most effective, with the best metrics in the estimation, the best performance was Case 1, with R2 = 0.89, an RSME of 0.32 µg/L, an MAE 1.25 µg/L and an MSE 0.25 (µg/L)2, consistently outperforming the others according to the static metrics used for validation. This finding underscores the effectiveness of LSTM in capturing the complex temporal relationships inherent in the dataset. However, increasing the dataset in Case 3 shows a better performance of TCNs (R2 = 0.96; MSE = 0.33 (µg/L)2; RMSE = 0.13 µg/L; and MAE = 0.06 µg/L). The successful application of machine learning algorithms emphasizes their potential to elucidate the dynamics of algal biomass in Lake Ranco, located in the southern region of Chile. These results not only contribute to a deeper understanding of the lake ecosystem but also highlight the utility of advanced computational techniques in environmental research and management. Full article
(This article belongs to the Special Issue Remote Sensing Band Ratios for the Assessment of Water Quality)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Location of Chile in South America; (<b>b</b>) XIV Región de los Ríos; (<b>c</b>) Lake Ranco. Monitoring stations shown in yellow circles; meteorological stations shown in orange triangles and bathymetry.</p>
Full article ">Figure 2
<p>Land use, road network, and populated areas within a 10 km buffer zone around Lake Ranco. (<b>a</b>) Town of Futrono; (<b>b</b>) town of Lake Ranco; and (<b>c</b>) town of Llifén.</p>
Full article ">Figure 3
<p>Recurrent neural network architecture. Adapted from [<a href="#B58-remotesensing-16-03401" class="html-bibr">58</a>].</p>
Full article ">Figure 4
<p>Long short-term memory architecture, adapted from [<a href="#B61-remotesensing-16-03401" class="html-bibr">61</a>].</p>
Full article ">Figure 5
<p>Gated recurrent unit architecture adapted from [<a href="#B63-remotesensing-16-03401" class="html-bibr">63</a>].</p>
Full article ">Figure 6
<p>TCN architecture adapted from [<a href="#B65-remotesensing-16-03401" class="html-bibr">65</a>].</p>
Full article ">Figure 7
<p>Meteorological conditions in Lake Ranco over 1989–2022.</p>
Full article ">Figure 8
<p>Chl-a estimation values Case 1 in each lake station.</p>
Full article ">Figure 9
<p>Chl-a estimation values Case 2 in each lake station.</p>
Full article ">Figure 10
<p>Chl-a estimation values for Case 3 at each lake station.</p>
Full article ">
28 pages, 15371 KiB  
Article
Research on the Spatial-Temporal Evolution of Changsha’s Surface Urban Heat Island from the Perspective of Local Climate Zones
by Yanfen Xiang, Bohong Zheng, Jiren Wang, Jiajun Gong and Jian Zheng
Land 2024, 13(9), 1479; https://doi.org/10.3390/land13091479 - 12 Sep 2024
Viewed by 326
Abstract
Optimizing urban spatial morphology is one of the most effective methods for improving the urban thermal environment. Some studies have used the local climate zones (LCZ) classification system to examine the relationship between urban spatial morphology and Surface Urban Heat Islands (SUHIs). However, [...] Read more.
Optimizing urban spatial morphology is one of the most effective methods for improving the urban thermal environment. Some studies have used the local climate zones (LCZ) classification system to examine the relationship between urban spatial morphology and Surface Urban Heat Islands (SUHIs). However, these studies often rely on single-time-point data, failing to consider the changes in urban space and the time-series LCZ mapping relationships. This study utilized remote sensing data from Landsat 5, 7, and 8–9 to retrieve land surface temperatures in Changsha from 2005 to 2020 using the Mono-Window Algorithm. The spatial-temporal evolution of the LCZ and the Surface Urban Heat Island Intensity (SUHII) was then examined and analyzed. This study aims to (1) propose a localized, long-time LCZ mapping method, (2) investigate the spatial-temporal relationship between the LCZ and the SUHII, and (3) develop a more convenient SUHI assessment method for urban planning and design. The results showed that the spatial-temporal evolution of the LCZ reflects the sequence of urban expansion. In terms of quantity, the number of built-type LCZs maintaining their original types is low, with each undergoing at least one type change. The open LCZs increased the most, followed by the sparse and the composite LCZs. Spatially, the LCZs experience reverse transitions due to urban expansion and quality improvements in central urban areas. Seasonal changes in the LCZ types and the SUHI vary, with differences not only among the LCZ types but also in building heights within the same type. The relative importance of the LCZ parameters also differs between seasons. The SUHI model constructed using Boosted Regression Trees (BRT) demonstrated high predictive accuracy, with R2 values of 0.911 for summer and 0.777 for winter. In practical case validation, the model explained 97.86% of the data for summer and 96.77% for winter. This study provides evidence-based planning recommendations to mitigate urban heat and create a comfortable built environment. Full article
Show Figures

Figure 1

Figure 1
<p>The location of the study area.</p>
Full article ">Figure 2
<p>Distribution of LCZ Parameters from 2005 to 2020.</p>
Full article ">Figure 2 Cont.
<p>Distribution of LCZ Parameters from 2005 to 2020.</p>
Full article ">Figure 3
<p>The semivariogram model of building height.</p>
Full article ">Figure 4
<p>Various schematic diagrams of local climate zones in Changsha City.</p>
Full article ">Figure 5
<p>The LCZ maps in the years 2005, 2010, 2015, and 2020.</p>
Full article ">Figure 6
<p>The spatial variation of the LCZ types from 2005 to 2020.</p>
Full article ">Figure 7
<p>The urban structural development directions from 2005 to 2020.</p>
Full article ">Figure 8
<p>Spatiotemporal distribution of the LST in Changsha in summer and winter in 2005, 2010, 2016 and 2020: (<b>a</b>) 2005, (<b>b</b>) 2010, (<b>c</b>) 2016, and (<b>d</b>) 2020 in summer; (<b>e</b>) 2005, (<b>f</b>) 2010, (<b>g</b>) 2016, and (<b>h</b>) 2020 in winter; A: Lugu High-Tech Industrial Park, B: Changsha Economic and Technological Development Zone, C: Changsha Tianxin Economic Development Zone.</p>
Full article ">Figure 9
<p>Changes of the SUHII in the LCZ in summer and winter in 2005, 2010, 2016, and 2020: (<b>a</b>) 2005, (<b>b</b>) 2010, (<b>c</b>) 2016, and (<b>d</b>) 2020 in summer; (<b>e</b>) 2005, (<b>f</b>) 2010, (<b>g</b>) 2016 and (<b>h</b>) 2020 in winter. The boxplots represent the variation of SUHII values for each LCZ type, while the strip plots indicate the mean SUHII value for each LCZ type.</p>
Full article ">Figure 10
<p>Relative influences of the LCZ parameters in the two seasons.</p>
Full article ">Figure 11
<p>BRT model’s prediction results.</p>
Full article ">Figure 12
<p>The location and LST of Wangcheng District.</p>
Full article ">
25 pages, 9415 KiB  
Article
Spatial and Seasonal Variation and the Driving Mechanism of the Thermal Effects of Urban Park Green Spaces in Zhengzhou, China
by Yuan Feng, Kaihua Zhang, Ang Li, Yangyang Zhang, Kun Wang, Nan Guo, Ho Yi Wan, Xiaoyang Tan, Nalin Dong, Xin Xu, Ruizhen He, Bing Wang, Long Fan, Shidong Ge and Peihao Song
Land 2024, 13(9), 1474; https://doi.org/10.3390/land13091474 - 11 Sep 2024
Viewed by 382
Abstract
Greenscaping, a key sustainable practice, helps cities combat rising temperatures and climate change. Urban parks, a pivotal greenscaping element, mitigate the urban heat island (UHI) effect. In this study, we utilized high-resolution remote sensing imagery (GF-2 and Landsat 8, 9) and in situ [...] Read more.
Greenscaping, a key sustainable practice, helps cities combat rising temperatures and climate change. Urban parks, a pivotal greenscaping element, mitigate the urban heat island (UHI) effect. In this study, we utilized high-resolution remote sensing imagery (GF-2 and Landsat 8, 9) and in situ measurements to analyze the seasonal thermal regulation of different park types in Zhengzhou, China. We calculated vegetation characteristic indices (VCIs) and landscape patterns (LMs) and employed boosted regression tree models to explore their relative contributions to land surface temperature (LST) across different seasons. Our findings revealed that urban parks lowered temperatures by 0.65 °C, 1.41 °C, and 2.84 °C in spring, summer, and autumn, respectively, but raised them by 1.92 °C in winter. Amusement parks, comprehensive parks, large parks, and water-themed parks had significantly lower LSTs. The VCI significantly influenced LST in autumn, with trees having a stronger cooling effect than shrubs. LMs showed a more prominent effect than VCIs on LST during spring, summer, and winter. Parks with longer perimeters, larger and more dispersed green patches, higher plant species richness, higher vegetation heights, and larger canopies were associated with more efficient thermal reduction in an urban setting. The novelty of this study lies in its detailed analysis of the seasonal thermal regulation effects of different types of urban parks, providing new insights for more effective urban greenspace planning and management. Our findings assist urban managers in mitigating the urban surface heat effect through more effective urban greenspace planning, vegetation community design, and maintenance, thereby enhancing cities’ potential resilience to climate change. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Location of Zhengzhou city, Henan, China. (<b>a</b>) Distribution of 123 selected parks in the study area versus 805 sampling points, (<b>b</b>–<b>f</b>) Distribution of sample sites in selected parks, (<b>b’</b>–<b>f’</b>) Selected parkland classification results.</p>
Full article ">Figure 2
<p>Flowchart of this study.</p>
Full article ">Figure 3
<p>(<b>a</b>–<b>d</b>) LST in different seasons (Units: °C); (<b>e</b>–<b>h</b>) Land surface temperature of parks in different seasons.</p>
Full article ">Figure 4
<p>Differences in LST categories among parks in different seasons. The letters a, b, and c denote significant disparities identified via Fisher’s least significant difference test (<span class="html-italic">p</span> &lt; 0.05) across various park types during different seasons. (<b>a</b>–<b>c</b>), (<b>d</b>–<b>f</b>), (<b>g</b>–<b>i</b>), and (<b>j</b>–<b>l</b>) represent the LST distribution of parks classified by different standards in spring, summer, autumn, and winter, respectively.</p>
Full article ">Figure 5
<p>Analysis of the differences in driving factors among different park types. The letters a, b, and c denote significant disparities identified via Fisher’s least significant difference test (<span class="html-italic">p</span> &lt; 0.05) among different park types. (<b>a</b>–<b>o</b>) are vegetation characteristic indices, and (<b>p</b>–<b>y</b>) are landscape pattern indices.</p>
Full article ">Figure 6
<p>Spearman correlation coefficients of influencing factors with LST across seasons. (**. Correlation is significant at the 0.01 level; *. Correlation is significant at the 0.05 level).</p>
Full article ">Figure 7
<p>Relative contribution of each influencing factor to surface temperature under different seasons.</p>
Full article ">Figure 8
<p>Relative importance of different park types for surface temperature in different seasons (%).</p>
Full article ">Figure 9
<p>Partial dependence plot of the driving factors’ impact on LST.</p>
Full article ">Figure A1
<p>The correlation between each driving factor and LST in different park types and seasons.</p>
Full article ">
27 pages, 10360 KiB  
Article
Soil Moisture-Derived SWDI at 30 m Based on Multiple Satellite Datasets for Agricultural Drought Monitoring
by Jing Ning, Yunjun Yao, Joshua B. Fisher, Yufu Li, Xiaotong Zhang, Bo Jiang, Jia Xu, Ruiyang Yu, Lu Liu, Xueyi Zhang, Zijing Xie, Jiahui Fan and Luna Zhang
Remote Sens. 2024, 16(18), 3372; https://doi.org/10.3390/rs16183372 - 11 Sep 2024
Viewed by 478
Abstract
As a major agricultural hazard, drought frequently occurs due to a reduction in precipitation resulting in a continuously propagating soil moisture (SM) deficit. Assessment of the high spatial-resolution SM-derived drought index is crucial for monitoring agricultural drought. In this study, we generated a [...] Read more.
As a major agricultural hazard, drought frequently occurs due to a reduction in precipitation resulting in a continuously propagating soil moisture (SM) deficit. Assessment of the high spatial-resolution SM-derived drought index is crucial for monitoring agricultural drought. In this study, we generated a downscaled random forest SM dataset (RF-SM) and calculated the soil water deficit index (RF-SM-SWDI) at 30 m for agricultural drought monitoring. The results showed that the RF-SM dataset exhibited better consistency with in situ SM observations in the detection of extremes than did the SM products, including SMAP, SMOS, NCA-LDAS, and ESA CCI, for different land cover types in the U.S. and yielded a satisfactory performance, with the lowest root mean square error (RMSE, below 0.055 m3/m3) and the highest coefficient of determination (R2, above 0.8) for most observation networks, based on the number of sites. A vegetation health index (VHI), derived from a Landsat 8 optical remote sensing dataset, was also generated for comparison. The results illustrated that the RF-SM-SWDI and VHI exhibited high correlations (R ≥ 0.5) at approximately 70% of the stations. Furthermore, we mapped spatiotemporal drought monitoring indices in California. The RF-SM-SWDI provided drought conditions with more detailed spatial information than did the short-term drought blend (STDB) released by the U.S. Drought Monitor, which demonstrated the expected response of seasonal drought trends, while differences from the VHI were observed mainly in forest areas. Therefore, downscaled SM and SWDI, with a spatial resolution of 30 m, are promising for monitoring agricultural field drought within different contexts, and additional reliable factors could be incorporated to better guide agricultural management practices. Full article
(This article belongs to the Special Issue Quantitative Remote Sensing of Vegetation and Its Applications)
Show Figures

Figure 1

Figure 1
<p>The spatial distribution of training and testing stations used in the downscaling framework. The map of land cover types of the substudy area and the locations of the in situ observation stations appear at the <b>top left</b> and <b>bottom</b>, respectively.</p>
Full article ">Figure 2
<p>Downscaling framework for the surface SM at 30 m through the integration of multiple datasets.</p>
Full article ">Figure 3
<p>Scatterplots of the comparison for the RF-SM data and SM derived from in situ observations at (<b>a</b>) 170 training stations and (<b>b</b>) 72 independent validation stations. The color indicates the density of the samples distributed in the area.</p>
Full article ">Figure 4
<p>Permutation importance of RF-SM. The features (i.e., input variables) include the SM products (SMAP, SMOS, ESA CCI, and NCA-LDAS), the soil properties (clay, sand, and silt), and the reflectance at visible and near-infrared bands (from SR_b4 to SR_b7), as well as the surface temperature (ST_b10) derived from Landsat 8.</p>
Full article ">Figure 5
<p>Boxplots of the in situ SM, RF-SM data, and the four SM products (SMAP, SMOS, NCA-LDAS, and KGE) for different land cover types. In the single boxplots, the red cross-dots denote outliers; the lowest and highest lines denote minimum and maximum results, respectively, except for extreme values (outliers); and the lower bound of the box, red line in the box, and upper bound of the box represent the lower quartile (25%), the median, and upper quartile (75%), respectively.</p>
Full article ">Figure 6
<p>Diagrams of the statistics (R<sup>2</sup>, RMSE, Bias, and KGE) for the comparison between the RF-SM dataset and the four SM products (SMAP, SMOS, NCA-LDAS, KGE) for the different observation networks.</p>
Full article ">Figure 7
<p>Temporal variations in precipitation (P) and surface SM derived from RF-SM and the four products at the representative stations in the substudy area during 2016.</p>
Full article ">Figure 8
<p>Spatial distributions of the RF-SM in the substudy area during 2016.</p>
Full article ">Figure 9
<p>Comparison between the RF-SM-SWDI and VHI based on the Pearson correlation coefficient (R) from 242 in situ stations in 2016.</p>
Full article ">Figure 10
<p>Temporal variations in SM-SWDI, RF-SM-SWDI, VHI, and precipitation (P) anomalies at the representative stations in the substudy area in 2016.</p>
Full article ">Figure 11
<p>Spatial distributions of the RF-SM-SWDI and VHI in the substudy area in 2016.</p>
Full article ">Figure 12
<p>Comparison between the RF-SM-SWDI and VHI based on the Pearson correlation coefficient (R) in the substudy area in 2016.</p>
Full article ">Figure 13
<p>Comparison between the RF-SM-SWDI and two VHI components: (<b>a</b>) VCI and (<b>b</b>) TCI, based on the Pearson correlation coefficient (R) in the substudy area in 2016.</p>
Full article ">Figure 14
<p>Spatial distributions of the RF-SM-SWDI, RF-SM-SWDI after resampling, and the short-term drought blend (STDB) in the substudy area in 2016.</p>
Full article ">
28 pages, 2705 KiB  
Article
Estimating Evapotranspiration of Rainfed Winegrapes Combining Remote Sensing and the SIMDualKc Soil Water Balance Model
by Wilk S. Almeida, Paula Paredes, José Basto, Isabel Pôças, Carlos A. Pacheco and Teresa A. Paço
Water 2024, 16(18), 2567; https://doi.org/10.3390/w16182567 - 10 Sep 2024
Viewed by 398
Abstract
Soil water balance (SWB) in woody crops is sometimes difficult to estimate with one-dimensional models because these crops do not completely cover the soil and usually have a deep root system, particularly when cropped under rainfed conditions in a Mediterranean climate. In this [...] Read more.
Soil water balance (SWB) in woody crops is sometimes difficult to estimate with one-dimensional models because these crops do not completely cover the soil and usually have a deep root system, particularly when cropped under rainfed conditions in a Mediterranean climate. In this study, the actual crop evapotranspiration (ETc act) is estimated with the soil water balance model SIMDualKc which uses the dual-Kc approach (relating the fraction of soil cover with the crop coefficients) to improve the estimation of the water requirements of a rainfed vineyard, using data from a deep soil profile. The actual basal crop coefficient (Kcb act) obtained using the SIMDualKc model was compared with the Kcb act estimated using the A&P approach, which is a simplified approach based on measurements of the fraction of ground cover and crop height. Spectral vegetation indices (VIs) derived from Landsat-5 satellite data were used to determine the fraction of ground cover (fc VI) and thus the density coefficient (Kd). The SIMDualKc model was calibrated using available soil water (ASW) measurements down to a depth of 1.85 m, which significantly improved the conditions for using an SWB estimation model. The test of the model was performed using a different ASW dataset. A good agreement between simulated and field-measured ASW was observed for both data sets along the crop season, with RMSE < 12.0 mm and NRMSE < 13%. The calibrated Kcb values were 0.15, 0.60, and 0.52 for the initial, mid-season, and end season, respectively. The ratio between ETc act and crop evapotranspiration (ETc) was quite low between veraison and maturity (mid-season), corresponding to 36%, indicating that the rainfall was not sufficient to satisfy the vineyard’s water requirements. VIs used to compute fc VI were unable to fully track the plants’ conditions during water stress. However, ingestion of data from remote sensing (RS) showed promising results that could be used to support decision making in irrigation scheduling. Further studies on the use of the A&P approach using RS data are required. Full article
(This article belongs to the Section Soil and Water)
Show Figures

Figure 1

Figure 1
<p>Study area location in Santarém, Portugal. (Vineyard approximate boundaries in black and sample collection area (SCA) in red). 1D (▲) and 2A (▄) are the locations of soil water content measurements used for calibrating and testing the SIMDualKc model, respectively. SAVI is the Soil Adjusted Vegetation Index.</p>
Full article ">Figure 2
<p>The annual cycle of the vine and the crop growth stages. The crop growth stages are delimited according to the FAO segmented curve [<a href="#B27-water-16-02567" class="html-bibr">27</a>].</p>
Full article ">Figure 3
<p>Average monthly reference evapotranspiration (ET<sub>o</sub>) and precipitation (P) for the period 1982–1987 and for the study year (1987) in Santarém, Portugal.</p>
Full article ">Figure 4
<p>Schematic partial top view (<b>a</b>) of the experimental layout and the access tubes (AT, blue circles) (<b>b</b>) locations at the experimental field, Santarém, Portugal.</p>
Full article ">Figure 5
<p>Available soil water dynamics for the calibration (<b>a</b>) and test (<b>b</b>) of the SIMDualKc model. Dots represent observations while the curve represents simulations of the available soil water (ASW). TAW represents the total available water, and RAW denotes the rapidly available water.</p>
Full article ">Figure 6
<p>Severity of water stress in a rainfed vineyard, according to the predawn leaf water potential (ψ<sub>p</sub>) limits proposed by [<a href="#B26-water-16-02567" class="html-bibr">26</a>], grown in Santarém, Central Portugal.</p>
Full article ">Figure 7
<p>Standard and actual basal crop coefficients (K<sub>cb</sub>, K<sub>cb act</sub>), soil evaporation coefficient (K<sub>e</sub>), and actual crop coefficient (K<sub>c act</sub> = K<sub>cb act</sub> + K<sub>e</sub>) and precipitation (P) computed by the SIMDualKc model in a rainfed vineyard in Santarém, Portugal.</p>
Full article ">Figure 8
<p>Simulated soil water balance components: precipitation, evapotranspiration (ET), soil water content variation (Δ ASW), runoff (RO), and capillary rise (CR) (all variables in mm) after accurate SIMDualKc model calibration.</p>
Full article ">
20 pages, 2679 KiB  
Article
Spatio-Temporal Analysis of Green Infrastructure along the Urban-Rural Gradient of the Cities of Bujumbura, Kinshasa and Lubumbashi
by Henri Kabanyegeye, Nadège Cizungu Cirezi, Héritier Khoji Muteya, Didier Mbarushimana, Léa Mukubu Pika, Waselin Salomon, Yannick Useni Sikuzani, Kouagou Raoul Sambieni, Tatien Masharabu and Jan Bogaert
Land 2024, 13(9), 1467; https://doi.org/10.3390/land13091467 - 10 Sep 2024
Viewed by 273
Abstract
This study analyses the dynamics of green infrastructure (GI) in the cities of Bujumbura, Kinshasa, and Lubumbashi. A remote sensing approach, combined with landscape ecology metrics, characterized this analysis, which was based on three Landsat images acquired in 2000, 2013, and 2022 for [...] Read more.
This study analyses the dynamics of green infrastructure (GI) in the cities of Bujumbura, Kinshasa, and Lubumbashi. A remote sensing approach, combined with landscape ecology metrics, characterized this analysis, which was based on three Landsat images acquired in 2000, 2013, and 2022 for each city. Spatial pattern indices reveal that GI was suppressed in Bujumbura and Kinshasa, in contrast to Lubumbashi, which exhibited fragmentation. Furthermore, the values of stability, aggregation, and fractal dimension metrics suggest that Bujumbura experienced rather intense dynamics and a reduction in the continuity of its GI, while Kinshasa showed weaker dynamics and tendencies towards patch aggregation during the study period. In contrast, Lubumbashi exhibited strong dynamics and aggregation of its GI within a context of significant anthropization. The evolution of the Normalized Difference Vegetation Index demonstrates a sawtooth pattern in the evolution of tall vegetation patches in Bujumbura, compared to a gradual decrease in Kinshasa and Lubumbashi. It is recommended that urban growth in these cities should be carefully planned to ensure the integration of sufficient GI. Full article
Show Figures

Figure 1

Figure 1
<p>Bujumbura is located in Burundi and Kinshasa, and Lubumbashi is in the Democratic Republic of the Congo.</p>
Full article ">Figure 2
<p>Land use maps of Bujumbura (Burundi), Kinshasa, and Lubumbashi (DRC) from supervised classification of Landsat images from 2000, 2013, and 2022 based on the Random Forest algorithm.</p>
Full article ">Figure 3
<p>Trends in number of patches, total area, average patch area, and vegetation class dominance for the cities of Bujumbura (Burundi), Kinshasa, and Lubumbashi (DRC) for the years 2000, 2013, and 2022.</p>
Full article ">Figure 4
<p>Normalized difference vegetation maps of the cities of Bujumbura (Burundi) (<b>A</b>), Kinshasa (<b>B</b>), and Lubumbashi (DRC) (<b>C</b>) for the years 2020, 2013, and 2022.</p>
Full article ">Figure 5
<p>Proportions of GI areas over NDVI intervals were found for the cities of Bujumbura, Kinshasa, and Lubumbashi for the years 2000, 2013, and 2022.</p>
Full article ">
Back to TopTop