1. Introduction
Deforestation and forest degradation in tropical regions, particularly in the Amazon rainforest, represent critical environmental challenges with far-reaching implications. The Amazon, the world’s largest tropical rainforest, serves as a crucial carbon storage unit and a hotspot for biodiversity, playing an essential role in regulating the planet’s climate and providing habitat for numerous species [
1]. However, this vital ecosystem faces unprecedented threats driven by complex socio-economic factors and land use changes [
2,
3].
Over the past few decades, the Brazilian Amazon has experienced alarming rates of forest loss [
4]. This deforestation severely threatens the forest and is primarily driven by agricultural expansion [
5], cattle ranching [
6], infrastructure development [
7], and illegal logging [
8]. Forest degradation is equally serious because it frequently precedes total deforestation. Degraded forests are much more susceptible to further deforestation, with typical decreases in canopy height and aboveground biomass of up to 60% and 65%, respectively, according to recent research conducted in tropical moist forests [
9].
Understanding the dynamics of land use and land cover (LULC) changes is essential for developing effective conservation strategies and sustainable land management practices. The Brazilian Amazon has experienced varying rates of change across its vast expanse [
4,
10]. Thus, this heterogeneity presents both challenges and opportunities for researchers and policymakers seeking to understand and mitigate the impacts of land use changes. The advent of satellite imagery and advanced classification techniques has revolutionized our ability to monitor and analyze LULC changes at various spatial and temporal scales [
11].
Remote sensing has emerged as a crucial tool for monitoring deforestation and forest degradation in the Amazon. Satellite imagery provides consistent, large-scale observations that enable researchers to track changes in forest cover over time [
12]. The advent of medium and high-resolution satellite imagery has revolutionized our ability to monitor and analyze LULC changes at various spatial and temporal scales [
13]. While lower resolution imagery (e.g., 30 m Landsat, 250–1000 m MODIS) offers advantages in terms of broader coverage and long-term historical records, it may overlook fine-scale disturbances. Medium-resolution imagery, such as Sentinel-2 with its 10 m bands, shows a balance between coverage and resolution, allowing for the detection of smaller-scale disturbances that might be missed by coarser-resolution sensors [
14]. These advancements in spatial resolution allow for more accurate detection of subtle changes in forest structure and composition, which is critical for understanding degradation processes [
15].
To monitor LULC change and deforestation in the Amazon, pixel-based methods have been traditionally the primary approach [
16]. For instance, the Brazilian National Institute for Space Research (INPE) uses a pixel-based approach in its PRODES (Monitoring Deforestation in the Brazilian Amazon by Satellite Project) system, which has been the official method for monitoring deforestation in the Brazilian Amazon since 1988 [
17]. This monitoring system utilizes Landsat satellite imagery to provide annual assessments of deforestation, offering information with a minimum mapping unit of 6.25 hectares. However, recent advancements in remote sensing technology and analysis techniques have opened new possibilities for more accurate and detailed monitoring. Object-based image analysis (OBIA) has proven particularly advantageous for detecting LULC changes. Unlike pixel-based methods, OBIA considers the spatial context of image objects, allowing for a more nuanced classification of complex landscapes [
18]. When combined with satellite imagery from Sentinel-2 and machine learning algorithms like Random Forest, OBIA can achieve high classification accuracies and provide more detailed information about forest structure and composition.
Despite the potential advantages of OBIA combined with Sentinel-2 imagery and Random Forest classification, relatively few studies have applied this combination for monitoring deforestation and forest degradation in the Amazon. This approach offers the potential for more accurate detection of changes in forest cover, which is particularly important for monitoring forest degradation. For instance, Souza-Filho et al. used OBIA to investigate the influence of mining projects on LULC changes in the southeastern area of the Brazilian Amazon, resulting in overall accuracies of all the classified maps higher than 94% [
19]. Bueno et al. demonstrated the effectiveness of Random Forest with OBIA for land cover change detection in the Brazilian savanna (Cerrado) biome, achieving overall accuracies of around 88% [
20].
Regarding the forest areas, conservation units play a vital role in protecting the Amazon’s biodiversity and ecosystem services. Established according to the Brazilian National System of Nature Conservation Units (SNUCs), these areas are crucial for preserving the forest’s integrity and helping mitigate climate change [
21]. However, their effectiveness can vary depending on management strategies and local contexts. Recent studies have shown that while protected areas generally reduce deforestation rates, their impact can differ significantly across regions and management types [
22,
23]. Understanding these variations is essential for optimizing conservation strategies.
The heterogeneity of the Amazon presents both challenges and opportunities for researchers and policymakers seeking to understand and mitigate the impacts of land use changes [
24]. This study focuses on two distinct areas within the Brazilian Amazon: Manaus and Porto Velho. These areas represent different stages of development and deforestation pressures, offering valuable insights into the diverse challenges faced across the Amazon.
The objectives of this study are as follows: (1) to produce accurate LULC maps for Manaus and Porto Velho from 2018 to 2023 using Sentinel-2 imagery and OBIA with Random Forest classification; (2) to analyze LULC distribution in the study areas; (3) to compare LULC changes and conservation effectiveness between Manaus and Porto Velho, considering their different deforestation and degradation yearly rates; and (4) to quantify the total area of deforestation and forest degradation of conservation units with different management strategies and their surrounding areas to evaluate their effectiveness.
By addressing these objectives, this study aims to contribute valuable insights to the ongoing discourse on Amazon conservation and inform evidence-based policymaking for sustainable forest management in the region. The use of advanced remote sensing techniques and machine learning algorithms provides a detailed and accurate assessment of LULC changes, while the comparison between different types of conservation units offers crucial information for improving protected area management strategies.
2. Materials and Methods
2.1. Study Areas
The Brazilian Amazon rainforest (
Figure 1b) is located in the northern region of Brazil and covers the territory of 9 states (
Figure 1a). It experiences a tropical climate (A-Zone in the Köppen classification) with average temperature in the coldest month exceeding 20 °C in the center of the Amazon basin. Rainfall can reach over 3000 mm annually, and in some areas, it may even exceed 8000 mm [
25]. The forest is dominated by evergreen broadleaved trees and forest structure is typically described as having multiple layers, with structures that could reach heights of around 20–30 meters. Dominant tree species include
Eschweilera coriacea,
Euterpe precatoria, and
Protium altissimum [
26].
This study incorporates two areas within the Brazilian Amazon region: Manaus City and the district of Porto Velho (
Figure 1).
Manaus is the capital city of Amazonas state (
Figure 1c), and its area encompasses the urban area and its surrounding forest, covering approximately 1,142,000 ha. Manaus is situated in the heart of the Brazilian Amazon. It is a major urban center with a population exceeding two million people. The area’s location at the confluence of the Negro and Solimões rivers makes it a critical hub for transportation and commerce. Manaus has experienced rapid urban expansion, leading to increased pressure on surrounding forested areas [
27].
The district of Porto Velho (
Figure 1d), within the municipality of Porto Velho, capital of Rondônia state, is located on the east bank of the Madeira River in the southwestern part of the Brazilian Amazon. The area encompasses the urban area and its environs, spanning approximately 923,000 ha. Unlike Manaus, Porto Velho is more rural, with a strong reliance on agriculture and livestock farming. The construction of major infrastructure projects, such as highways and hydroelectric dams, has contributed to deforestation in the region [
27].
The district area of Porto Velho was chosen for this study, rather than the entire municipality, due to logistical and computational constraints. The district represents an important region of the Madeira River basin (Médio Madeira). Additionally, it houses the urban administrative hub of the municipality, and is comparable in size to Manaus, allowing for a more balanced comparison. The district contains 13 conservation units within its boundaries.
This study examines a total of 29 conservation units (
Table 1) across different management levels—sustainable use and full protection—and administrative scales—local, regional, and national—to assess their effectiveness in preventing deforestation and forest degradation. According to the SNUC legislation [
28], in full protection units, only indirect use of natural resources is allowed, and the rules and regulations are restrictive. Sustainable use units, on the other hand, reconcile nature conservation with the sustainable use of part of natural resources.
Buffer zones around these conservation units were also included in the analysis, as they play a crucial role in the management and protection of these areas [
29]. According to the SNUC, buffer zones are areas surrounding conservation units where human activities are subject to specific rules and restrictions to minimize negative impacts on the protected area. For this study, a three-kilometer buffer was considered based on the management plans of most units. However, it is important to note that the SNUC excludes Environmental Protection Areas (a subcategory of sustainable use units) from having buffer zones, as these areas often cover large territories with diverse land uses, including private properties and urban areas.
2.2. Data
2.2.1. Sentinel 2 Imagery Acquisition
This study utilized Sentinel-2 multispectral imagery for the period 2018–2023, chosen for its free availability, spatial resolution (10 m, 20 m, and 60 m depending on the band), and frequent revisit time, making it ideal for monitoring LULC changes. The Sentinel-2 MSI: Multi-Spectral Instrument Level-2A collection [
30] was accessed through the Google Earth Engine (GEE) platform, using the Cloud Score+ dataset [
31]. These products include radiometric and geometric corrections, providing highly accurate geolocated images without the need for further preprocessing, such as atmospheric correction.
The Cloud Score+ dataset functions as a quality assessment (QA) processor for optical satellite imagery, effectively removing clouds and cloud shadows while identifying relatively clear pixels. It employs two quality bands, cs, and cs_cdf, which rate each pixel’s usefulness about surface visibility on a continuous scale from 0 to 1, where 0 denotes “not clear” observations and 1 denotes “clear” observations.
To mitigate the impact of extensive cloud cover typically observed in tropical regions, median composite images were created from a group of images within a specific time interval. Sentinel-2 images were filtered from the collection over 4 months, from the beginning of May to the end of September, corresponding to the dry season in the Amazon region. This temporal selection was made to minimize cloud cover interference and ensure the highest quality data for analysis.
The resulting dataset consists of six annual median composites, corresponding to the years 2018 to 2023. These composites provide a consistent and cloud-free representation of the study area, enabling accurate detection and analysis of LULC changes over the study period. The use of median composites helps reduce the influence of outliers and temporary land cover changes, providing a more stable representation of the landscape for each year.
2.2.2. Training and Validation Samples
The study areas were classified into distinct land cover classes. Manaus was assigned eight classes: (1) Terra-Firme Forest, (2) Floodplain Forest, (3) Secondary Forest, (4) Agriculture/Pastureland, (5) Burned Areas, (6) Barren Land, (7) Development Areas, and (8) Water Bodies. Porto Velho included all these classes with the addition of (9) Savanna areas, which are found in small patches throughout the region. Detailed descriptions of these land use and land cover classes can be found in
Table 2.
Training and validation samples were collected through a combination of methods to ensure accuracy and representativeness. The primary method involved the manual visual assessment of high-resolution imagery from Google Earth Pro and Planet Labs, supplemented by spectral and vegetation indices analysis. This approach allowed for a comprehensive evaluation of land cover types across the study areas. Additionally, to enhance the ground-truth data, 23 sample points were collected in 2022 using a Garmin GPS (Garmin Ltd., Olathe, KS, USA) instrument in both Manaus and Porto Velho.
To mitigate the impact of spatial autocorrelation while still capturing the gradient of each land cover type, all samples were collected as points rather than polygons. Sample locations were chosen arbitrarily, with the constraint that they be at least 10 m apart from the nearest sample point. This strategy helps to ensure independence between samples and improves the robustness of the classification model.
A total of approximately 3000 sample points were collected for each year of the study period. This substantial dataset was then divided into two subsets: 70% of the samples were used to train the classification model, while the remaining 30% were reserved for validation and accuracy assessment.
2.3. Segmentation and Object-Based Image Classification
This study employed an object-based image analysis (OBIA) approach for LULC classification using Google Earth Engine (GEE). OBIA was chosen for its ability to reduce the “salt-and-pepper” effect common in pixel-based classifications and to incorporate spatial context into the classification process [
18,
32]. The OBIA method consists of two primary steps: image segmentation and object classification.
The segmentation process divides the entire image into multiple areas of varying sizes based on distinct spectral and textural properties. These segments are then merged to form larger objects corresponding to land cover classes. During this process, the geometric properties, topology, and adjacency relationships of the objects are exploited, with the anticipation that these objects will correspond more easily to land cover types than individual pixels would.
For the segmentation step, this study utilized the Simple Non-Iterative Clustering (SNIC) algorithm, which efficiently groups similar pixels and identifies potential individual objects [
33]. SNIC is initiated using a uniform grid of seeds, generated by the “Image.Segmentation.seedGrid” function, which requires specifying a superpixel seed location spacing in pixels to influence the size of the resulting clusters. The optimal spacing value can be determined through experimentation. The algorithm then identifies objects (clusters) based on input parameters and produces a multi-band raster output, which includes the clusters themselves and additional layers containing average values of the input features [
34]. SNIC requires several key parameters for setting:
Compactness factor: Influences cluster shape, with larger values producing more compact clusters.
Connectivity: Can be set to 4 or 8, determining whether to use Rook’s or Queen’s contiguity for merging adjacent clusters.
Neighborhood size: Helps avoid tile boundary artifacts.
After visual inspection of different combinations of parameter values and considering the characteristics of the LULC classes in the study areas, the following set of parameters was used in this study: compactness = 0, connectivity = 8, and neighborhood size = 256. Following segmentation, the OBIA approach combines spectral and spatial information with texture and contextual information from the image to perform the final object classification [
35].
2.4. Random Forest Algorithm
This study employed the Random Forest (RF) classifier, a robust and widely used machine learning algorithm, for the classification of segmented image objects. The Random Forest algorithm has been proven to be a reliable and accurate classification technique in numerous remote sensing applications [
36,
37,
38,
39,
40]. It was selected for this study due to its ability to handle high-dimensional data effectively and its capacity to assess feature importance.
Random Forest is an ensemble machine learning method designed to enhance the performance of Classification and Regression Trees (CARTs). It operates by constructing multiple decision trees and aggregating their outputs. In a classification analysis, each tree casts a vote for the most probable class of the input data, with the final classification determined by a majority vote. The RF algorithm employs two key techniques: bagging and random subspace selection. The process begins with the creation of multiple binary classification trees (
ntrees) using bootstrap samples drawn with replacements from the original dataset [
41]. For this study, the
ntrees parameter was set to 100 trees in classification mode, balancing computational efficiency with classification accuracy.
An important feature of the Random Forest algorithm is its use of out-of-bag (OOB) samples. These are observations not included in a particular bootstrap sample during the tree-building process. OOB samples serve a dual purpose in the RF model: they help estimate the misclassification error, provide an internal validation mechanism, and assist in assessing the importance of different variables in the model.
2.5. Auxiliary Features
To enhance classification accuracy and capture the complex landscape characteristics of the study areas, a range of spatial features were derived from spectral bands or acquired from various datasets. This study focused on the Blue, Green, Red, NIR, SWIR1, and SWIR2 bands of Sentinel-2 imagery, excluding other spectral channels to optimize computational efficiency while maintaining relevant spectral information. Band 1 (Coastal Aerosol) was omitted due to its primary application in studying coastal waters and atmospheric aerosol properties, which were not pertinent to the objectives of this study.
Topographic features, including elevation and slope, were extracted from the Shuttle Radar Topography Mission (SRTM) digital elevation dataset provided by NASA (Washington, DC, USA). These features contribute valuable information about the terrain characteristics that can influence land cover patterns. Principal Component Analysis (PCA) was applied to the spectral bands to reduce data dimensionality and enhance information content. PCA is widely used in classification and change detection analyses, particularly in unsupervised classification, as it eliminates redundant information and improves data representation. The first three PC levels were utilized in this study, providing a compact representation of the most significant spectral variations in the imagery.
Five spectral indices (
Table 3) were derived from the Sentinel-2 multispectral bands to enhance specific land cover characteristics [
42]. These indices were chosen due to their prominent use in LULC classification studies [
43,
44] and their ability to differentiate between various land features, particularly forest classes, from other types of land covers.
Textural information was incorporated using the Gray-Level Co-occurrence Matrix (GLCM), an effective method for extracting textural indices even from grayscale images [
50]. The texture attributes used in this study included Homogeneity (
‘idm’), Angular Second Moment (
‘asm’), Sum Average (
‘savg’), Entropy (
‘ent’), Contrast (
‘contrast’), and Correlation (
‘corr’).
Lastly, to account for seasonal variations in land cover, a multi-temporal seasonal image collection from the wet season was also included. Maximum and median values were extracted from this collection for each year, along with NDVI and NDWI values, to capture seasonal dynamics that might influence land cover classification.
2.6. Accuracy Assessment
The accuracy of the Random Forest (RF) classification was evaluated using a confusion matrix (CM), a widely used methodology in remote sensing for comparing classification outputs with ground-truth data. The predicted classes produced by the RF algorithm were compared to the labeled points in the testing dataset to assess the classification outcomes comprehensively. From the confusion matrix, several specific accuracy measures were derived [
51,
52]. Overall Accuracy (OA) represents the probability that a randomly chosen point on the map will be correctly classified. It is expressed as a percentage value, providing a general indication of the classification’s performance across all classes.
Producer’s Accuracy (PA), also known as “recall”, indicates the probability that a specific type of land cover on the ground is accurately classified on the map. It is computed by dividing the number of correctly identified pixels in each category by the total number of pixels in that category, as determined from the reference data. User’s Accuracy (UA), also known as “precision”, represents the probability that a pixel labeled as belonging to a particular category on the map truly reflects that category on the ground. It is calculated by dividing the number of correctly identified pixels within a category by the total number of pixels classified in that category.
The Kappa coefficient of agreement is a statistic that measures the agreement between classification and ground-truth values, taking into account the agreement that occurs by chance. It provides a more robust measure of classification accuracy, especially when comparing different classification results. These accuracy measures collectively provide a comprehensive assessment of classification performance, allowing for the evaluation of both overall classification accuracy and class-specific accuracies.
2.7. Change Detection Analysis
This study employed post-classification comparison, a widely used change detection method, implemented through QGIS 3.36 and ArcGIS Pro 2.8.0 post-processing tools. This approach has been shown to accurately represent land use changes in previous studies [
53,
54]. Change detection was carried out within the two areas, their consecutive conservation units, and buffer areas.
The analysis involved cross-tabulating classification results for consecutive years (2018–2019, 2019–2020, 2020–2021, 2021–2022, 2022–2023) and the overall study period (2018–2023). This process identified how each land cover class changed between consecutive years and between the initial and final dates of the study period. Transition maps were then generated to visualize these changes, providing a spatial representation of land cover dynamics. To calculate changes inside conservation units, polygons representing the boundaries of the units were used to extract LULC change information and summarize deforestation and forest degradation.
To assess the frequency of change, all six LULC maps were combined into a single raster in ArcGIS 2.8.0 Pro using the Combine tool. Each unique combination of LULC classes across all maps was assigned a unique value, with class pixel value information remaining for each year, which was organized in columns, or “fields”, in the raster’s attribute table. Using Field Calculation, new columns were created and changes were calculated between consecutive years. Pixels that underwent changes were assigned a value of 1, while pixels with no changes were assigned a value of 0. Subsequently, this information was aggregated, and a new raster was created using the Frequency tool to summarize the change frequency across the study areas. The tool sums all changes that occurred over the years, assigning each pixel a value between 0 (no changes) and 5 (changed every year), providing a comprehensive view of landscape dynamics throughout the study period.
2.8. Deforestation and Forest Degradation Detection
To detect deforestation and forest degradation in the study sites and conservation units, land cover classes were reclassified into broader categories: Primary Forest (Terra-Firme Forest and Floodplain Forest), Secondary Forest, and Non-Forest (Agriculture/Pastureland, Development Areas, Barren Land, and Burned Areas). Water Bodies and Savanna classes were excluded from the analysis to focus on forest dynamics. Using the change detection analysis, deforestation was then calculated from the conversion of Primary Forest classes to Non-Forest classes, while degradation was quantified as the transition from Primary Forest classes to the Secondary Forest class. Annual rates of deforestation and forest degradation were calculated for the whole study area. Total areas of deforestation and forest degradation within 5 years, using the initial and final land cover maps, were calculated for conservation units and buffer zones.
5. Conclusions
This study provided a comprehensive analysis of deforestation and forest degradation using LULC changes in Manaus and Porto Velho, Brazil, within and around conservation units from 2018 to 2023. By employing advanced remote sensing techniques and object-based classification methods, this study has achieved high-accuracy mapping of LULC dynamics in the Brazilian Amazon landscape.
Our findings reveal several key insights:
1. Despite the continued dominance of Terra-Firme Forest in both study areas, significant forest loss and degradation were observed, highlighting ongoing pressures on ecosystems in the Brazilian Amazon. This study revealed a gain of 22,362 ha in Secondary Forest areas in Manaus and 29,088 ha in Agriculture/Pastureland in Porto Velho within the study period.
2. The effectiveness of conservation units varied considerably across management types and administrative levels. Full protection units generally showed lower deforestation compared to sustainable use units, but both categories still face substantial challenges in preventing forest degradation.
3. Buffer zones around conservation units experienced larger deforestation and degradation compared to the conservation units, highlighting the effectiveness of these protected areas.
4. The conversion of forest to Agriculture/Pastureland emerged as the primary driver of land cover change and forest loss, reflecting the ongoing tension between conservation and agricultural expansion in the Amazon.
5. Forest degradation consistently surpassed deforestation in most areas.
These results have important implications for conservation policy and practice in the Brazilian Amazon. They suggest that while protected areas play a crucial role in forest conservation, current strategies may be insufficient to fully safeguard these ecosystems.
The high-accuracy LULC classification achieved in this study demonstrates the potential of advanced remote sensing and machine learning techniques for monitoring complex tropical forest landscapes. These methods provide valuable tools for tracking progress toward conservation goals and informing adaptive management strategies.
In conclusion, while conservation units in Manaus and Porto Velho play a vital role in forest protection, they face significant challenges from both internal and external pressures. By providing a detailed analysis of LULC dynamics in these areas, this study contributes to the evidence base needed to develop more effective and resilient conservation policies for these areas in the Brazilian Amazon.