With automatic delineation, we verified that the eastern direct influence area of Iguaçu National Park in the Gonçalves Dias River sub-basin region has a total surface area of 49,522 ha, with 23,106 ha (46.66%) belonging to Iguaçu National Park and 26,417 ha (53.34%) corresponding to the agricultural matrix surface. The landscape of the Gonçalves Dias River sub-basin presents an abrupt transition between the forest cover of the eastern extremity of Iguaçu National Park and the mosaic that makes up the agricultural matrix of the eastern buffer zone.
We found that by 1986, the Gonçalves Dias river sub-basin was already in an advanced fragmentation process. As reported by Ribeiro et al. (2009) and Teixeira et al. (2009), classification based on TM and OLI sensor scenes enabled the identification of fragments with tree structure, present for at least 15 years in the agricultural matrix of the Gonçalves Dias river sub-basin.
In the direct influence area with agricultural matrix surface, supervised classification identified 3,378 ha (OA = 85.33%, K = 82) of native forest cover in the year 1986, representing 12.78% of the total. In the year 2001, we observed 3,401 ha (OA = 89.33%, K = 87%) representing 12.87%, and in the year 2016, 4,158 ha (OA = 90.28%, K = 88%) representing 15.73% of native forest cover (Fig. 2).
According to Roderjan et al. (2002), Paraná had 83% of its territory covered by native vegetation; however, the process of colonization and expansion of agricultural frontiers, more intense in the northern and western regions, resulted in the current scenario, where only 11.7% of its territory maintains remnants of the Atlantic Forest (Atlântica, 2016).
With the analysis of slope on the buffer zone surface, we observed a predominance of flat to gently undulating terrain (10,653.80 ha or 40.33%), followed by undulating terrain (9,147.79 ha or 34.63%), strongly undulating terrain (6,498.20 ha or 24.60%), and only 116.97 ha (0.44%) in mountainous terrain.
We observed the influence of physical restrictions directed by IPARDES (1995) for the implementation of agricultural and forestry activities, as native forest cover showed a larger area on strongly undulating terrain 2,030 ha (49%), followed by undulating terrain 1,321 ha (32%), flat to gently undulating terrain 730.58 ha (17%), and mountainous terrain with only 76.5 ha (2%) of the agricultural matrix in the buffer zone of Iguaçu National Park.
We found that between 0 to 20% slope, corresponding to 74.96% (19,801.59 ha) of the agricultural matrix surface, only 2,052.36 ha (10.36%) of native forest cover were identified. In these areas, a lower rate of native forest cover was expected, considering that the only restriction for land use is the implementation of conservation practices (Borges & Rezende, 2023).
In contrast, when measuring the area located on slopes above 20%, it was found that, together, they correspond to 25.04% of the agricultural matrix (6,615.17 ha), of which 2,106.24 ha (31.84%) maintain native forest cover. There is a total restriction for mechanized agriculture on slopes, and even non-mechanized agriculture could only be admitted under moderate to severe restrictions.
Only 116.97 ha of surface area was identified in mountainous terrain (above 45%), however, only 76.41 ha maintain native forest cover, exposing a deficit of 40.56 ha. This deficit exposes probable areas for the restoration of native vegetation cover, as according to the new Forest Code (BRASIL, 2012), slopes or areas where the gradient exceeds 45% are designated as permanent preservation areas (PPAs), rendering them ineligible for exploitation or certain land use activities under environmental regulations.
In addition to physical restrictions, these areas highlight generally shallow soils (IPARDES, 1995), which can enhance erosive processes and degradation throughout the entire watershed (Rodrigues et al., 2011), negatively altering the quantity, shape, and quality of water resources offered upstream (Tucci & Mendes, 2006).
After vectorization, the native vegetation identified in the agricultural matrix in the Landsat image of 2016 depicted 204 native forest fragments, totaling 4,158.60 ha (15.74% of the total sub-basin area). The indices generated in V-LATE for the size classes of the fragments are grouped in Table 1.
Table 1
Landscape ecology indices for fragment size classes in the agricultural matrix of the Gonçalves Dias river sub-basin. Abbreviations: CA (Total area of fragments); MPS (Mean patch size); PSSD (Standard deviation of patch size); PSCoV (Coefficient of variation of patch size); MSI (Mean shape index); MFRACT (Mean fractal dimension); MPI (Mean proximity index).
Size classes (ha) | CA area (ha) | Number of fragments | Density and size | Shape | MPI proximity |
| | | MPS (ha) | PSSD (ha) | PSCoV (%) | MSI | MFRACT | |
3 a 5 | 260.3 | 68 | 3.83 | 0.54 | 14.09 | 2.07 | 1.38 | 0.65 |
5 a 10 | 299.57 | 44 | 6.81 | 1.38 | 20.31 | 2.28 | 1.37 | 0.22 |
10 a 15 | 287.52 | 24 | 12.4 | 1.26 | 10.13 | 2.50 | 1.37 | 0.36 |
15 a 20 | 327.07 | 19 | 17.21 | 1.54 | 8.93 | 2.83 | 1.38 | 0.95 |
20 a 25 | 222.12 | 10 | 22.21 | 1.62 | 7.29 | 3.30 | 1.40 | 0.40 |
25 a 30 | 211.77 | 8 | 26.47 | 1.25 | 4.73 | 3.12 | 1.38 | 0.23 |
30 a 50 | 453.04 | 11 | 41.19 | 6.41 | 15.58 | 3.27 | 1.38 | 4.53 |
50 a 100 | 883.59 | 13 | 67.97 | 13.33 | 19.62 | 3.72 | 1.38 | 1.29 |
100 a 250 | 941.54 | 6 | 156.92 | 31.08 | 19.81 | 4.37 | 1.38 | 10.02 |
250 a 500 | 262.08 | 1 | 262.08 | 0.00 | 0.00 | 4.56 | 1.38 | 0.00 |
In the eastern direct influence area with agricultural matrix, fragments with an area below 50 ha (n = 186) contribute to 91% of the total number of fragments and 50% of the native forest cover (2,087.81 ha), while forest fragments with an area larger than 50 hectares (n = 18) contribute only 9% of the total fragments and correspond to 50% of the total area (2,087 ha) (Fig. 3).
The quantity obtained in this study corroborates with the results of Ribeiro et al. (2009), Juvanhol et al. (2011), Pirovani et al. (2014), Azevedo, Gomes, and Moraes (2016), Mello et al. (2016), which found a higher representation of fragments with an area smaller than 50 ha in studies in the Atlantic Forest.
The fragmentation process is closely linked to the expansion of agricultural frontiers, mechanization, and anthropization (Mello et al., 2016; Silva et al., 2016). As a consequence of changes in the matrix, alterations occur in the internal environments of the fragments and in the canopy composition.
In Atlantic Forest fragments, especially those with an area of less than 10 ha, the canopy is dominated by a small number of pioneer or opportunistic gap trees; this dominance increases with decreasing fragment size, resulting in the disappearance of the emergent layer in the smaller, older fragments (Tabarelli et al., 2002).
At the landscape scale, the mean shape index (MSI) was 2.59, the proximity index (MPI) was 1.05, and the mean fractal dimension index (MFRACT) was 1.38. Similar to findings by Mello et al. (2016), in the Gonçalves Dias river sub-basin, larger fragments tend to exhibit greater complexity in their shape, and the mean fractal dimension index (MFRACT) shows that the landscape of the sub-basin is composed of fragments of moderately complex forms.
Our value of the mean fractal dimension index (1.38) was similar to other studies conducted in the Atlantic Forest, such as those obtained by Pirovani et al. (2014) in the Itapemirim River basin, Espírito Santo (between 1.36 and 1.40), and by Azevedo et al. (2016), who obtained 1.32 in the Buranhém River basin, extreme south of Bahia.
Gomide and Lingnau (2008) obtained an average fractal dimension index of 1.10 for fragments located in the APA do Iraí, Curitiba, Paraná, the same value obtained by Juvanhol et al. (2011) for the ecological corridor between the Forno Grande and Pedra Azul State Parks, Espírito Santo.
The average degree of isolation between fragments (MPI) showed variation among size classes. We observed a higher degree of isolation among fragments belonging to the "3 to 5 ha" class (MPI = 0.65) than among the "5 to 10" (MPI = 0.22), "10 to 15" (MPI = 0.36), "20 to 25" (MPI = 0.40), and "25 to 30" (MPI = 0.23) classes.
As fragmentation progresses and promotes habitat loss, reducing the original area of the fragments, the remaining small fragments become more accessible to hunting and timber exploitation. Similarly, the percentage of forest under edge effects increases, leading to the establishment of ruderal species, further exacerbating habitat loss (Tabarelli et al., 2002).
Field studies in the Atlantic Forest have accumulated evidence of synergy among phenomena correlated with habitat loss, fragmentation, edge effects, forest product exploitation, elimination of disperser species, and the frequency of forest fires (Castro et al., 2021; Lacerda et al., 2022; Scussel et al., 2020).