Restoration and Conservation of Priority Areas of Caatinga’s Semi-Arid Forest Remnants Can Support Connectivity within an Agricultural Landscape

Land-use and land-cover (LULC) changes are major drivers of biodiversity loss in semi-arid regions, such as the Caatinga biome located in the Northeast of Brazil. We investigated landscape dynamics and fragmentation in an area of the São Francisco Valley in the Brazilian Caatinga biome and measured the effect of these dynamics on ecological, functional and structural connectivity over a 33-year period (1985–2018). We calculated landscape connectivity indices based on graph theory to quantify the effect of further agricultural expansion on ecological connectivity at the landscape scale. We used a multicriteria decision analysis that integrates graph-based connectivity indices at the habitat patch scale, combined with an index of human disturbance to identify patches that, if conserved and restored, preserve the connectivity of the landscape most effectively. In the period studied, agriculture increased at a rate of 2104 ha/year, while native Caatinga vegetation decreased at a rate of 5203 ha/year. Both dense and open Caatinga became more fragmented, with the number of fragments increasing by 85.2% and 28.6%, respectively, whilst the average fragment size decreased by 84.8% and 6.1% for dense and open Caatinga, respectively. If agriculture patches were to expand by a 300 m buffer around each patch, the overall ecological connectivity could be reduced by 6–15%, depending on the species’ (smallto mid-size terrestrial vertebrates) mobility characteristics for which the connectivity indices were calculated. We provided explicit spatial connectivity and fragmentation information for the conservation and restoration of the Caatinga vegetation in the studied area. This information helps with conservation planning in this rapidly changing ecosystem.

Caatinga vegetation is under a lot of land-use pressure, and present high levels of chronic, human-caused disturbance [35].
The development of sustainable agriculture policies requires the increase in the understanding of the ecological impacts on biodiversity at the landscape level of farming systems [40]. In the Caatinga biome, conservation programs around agriculture land need to urgently stop species and natural ecosystem loss [34]. To fully understand the impact of agricultural development on Caatinga biodiversity and provide a quantification of recent changes and the current landscape situation, the aim of this paper is to analyze the changes in composition, configuration, and connectivity of the São Francisco Valley landscape over time. We anticipate that intensive agriculture is the main driver of land-use change within this region [33], and increases progressively over time, leading to reduced native vegetation cover (decreasing structural connectivity), greater fragmentation and a decline in ecological connectivity (increased barriers for terrestrial animal movement). In this paper, we combined Landsat images with a set of high-resolution land cover maps from the MapBiomas project collection 4.1 [41], supported by ground-truthing from field surveys, to address the following questions: (1) How have the cover and spatial configuration of Caatinga and agricultural habitat changed over the period of 1985-2018 in the São Francisco Valley landscape? (2) How would different patterns of future hypothetical continuing of agricultural expansion influence the connectivity of the landscape from the perspective of species with different levels of dispersal mobility? (3) If conserved or restored, which native vegetation areas preserve the connectivity of the landscape most effectively under agricultural expansion?

Study Area
The study area is located in the Brazilian Caatinga biome, in the northeast of the country. Within the Caatinga, we selected an area covering 846,400 ha that lies between the following geographic coordinates: 8 • 58 4 and 9 • 39 7 S; and 40 • 00 2 and 41 • 10 3 W. This area is located at the sub middle section of the São Francisco River basin, between the federal states of Bahia (BA) and Pernambuco (PE) (Figure 1). The study area includes part of the São Francisco River, urban areas of the cities of Casa Nova, Sobradinho, Petrolina, Juazeiro and Lagoa Grande, pastures, native vegetation and irrigated agriculture [42].
The native vegetation is a mosaic of arboreal or shrub forests, comprising mainly trees and low shrubs, many of which have thorns, microphylla, and some xerophytic features, as well as open areas [27,29,43]. The Caatinga's natural vegetation is characterized by leaf shedding in the dry season (August-November) [44] and is categorized as a dry tropical forest, formed by two types: a dense hyperxerophyllous forest and an open forest in different stages of development, featuring poaceae and open areas [44]. Caatinga can be broadly classified into two types and we followed the classification that has previously been used [41,42], which divided Caatinga into (1) the dense hyperxerophyllous forest, which has a woody cover of approximately 60-80% (dense Caatinga); and (2) the open forest, which is characterized by being in different stages of development, featuring grasses and abandoned areas, and has a woody cover of approximately 20-60%. Natural open areas are common around rock formations and in the areas too dry to sustain trees [32].
More recently, open areas that are characterized by low vegetation cover, density, and basal area, have been linked to historical clear-cuts for timber and high grazing pressures by livestock [34]; these can be characteristic of Caatinga degradation and are common in our study region [33,34]. In this area, the maximum temperature ranges from 29.7 • C (July) to 34.2 • C (November), with an annual mean of 32.3 • C; the minimum temperature ranges from 20.0 • C (July) to 23.5 • C (December) with an annual mean of 22.2 • C [32]. Mean monthly rainfall ranges from 1.4 mm (August) to 114.1 mm (March); most of it (70.4%) falls between January and April, with an annual mean of 482.6 mm [32].

Processing and Classification of Satellite Images
A 33-year spatiotemporal analysis (1985, 2000 and 2018) was conducted to analyze temporal changes in the landscape structure. Working on the Google Earth Engine (GEE) platform and Landsat images (sensors TM, ETM+ and OLI), we used the "best available pixel" method [45] to generate image composites with clear skies and the greatest possible absence of clouds [46] that were representative of the dry season, between August and November [47]. Surface reflectance products were subjected to topographic, radiometric, and atmospheric corrections [48,49].
To classify the Landsat images, we performed field surveys in 2019 to collect control points distributed in the agroecosystems and the native vegetation within the study area. These control points were complemented with LULC data obtained by photointerpretation of high-resolution images available from GEE and LULC maps from MapBiomas project collection 4.1 [41]. These maps cover the entire territory of Brazil at a spatial resolution of 30 m, from 1985 to 2018, with more than 80% accuracy [41]. MapBiomas LULC maps are classified into 19 classes and the native vegetation of the study area is classified as savanna [41]. A random forest (RF) supervised learning algorithm [50] was applied to classify the images in different LULC classes, obtaining one image for each year considered [47,51,52]. We classified the native vegetation into two classes (open and dense Caatinga) according to their structure.

Landscape Composition and Configuration Analysis
The resulting LULC maps were used to estimate changes in the spatial patterns during the study period . The quantification and comparison of the spatial layout of native vegetation fragments were carried out based on the following set of landscape metrics, which were chosen after reviewing native vegetation fragmentation studies [15,54,56]: (1) number of patches; (2) mean patch area (ha); (3) largest patch index (landscape percentage covered by the largest patch); and (4) the mean Euclidean nearest neighbor distance (m). These metrics were quantified in R version 4.0.4 [57], using the "landscape metrics" package [58].

Determining the Distance Threshold for Landscape Connectivity Analysis
To determine graph-based connectivity indices, it is necessary to choose or calculate a threshold of dispersal distances of the different species of interest [19,59,60]. We calculated two graph-based landscape connectivity indices-the "Integral Index of Connectivity" (IIC) and the "Class Coincidence Probability" (CCP)-that allowed us to define distance thresholds between fragments [20]. The IIC integrates into a single value of connectivity, the area of fragments of native vegetation (hereafter habitat patches), and the connectivity of the landscape [61]. The habitat patches of the focal species are represented as nodes, and the topological relations or connections between them are represented as links based on graph theory [62]. The values of IIC range from 0 to 1, increasing as the connectivity improves, thus a value of 1 would be reached in the hypothetical case of a landscape being totally occupied by the habitat. The IIC is calculated as follows [20]: where a (i or j) is the area of each patch; nl ij is the number of links in the shortest path (topological distance) between patches i and j. For patches that are not connected (belong to different components), the numerator in the sum of the equation equals zero (nl ij = ∞). When i = j then nl ij = 0 (no links needed to reach a certain patch from itself). n is the total number of patches in the landscape; A L is the total area of the landscape.
CCP is defined as the probability that two randomly selected points in a habitat belong to the same network [20]. CCP is calculated as follows: where NC is the total number of components in the landscape, where a component (connected region) is a set of nodes for which a path exists between every pair of nodes (an isolated patch makes up a component itself); c i is the total area of each component (i.e., the sum of the areas of the patches belonging to that component); and A C is the total area of patches in the landscape. CCP values range from 0 to 1. Higher values indicate higher connectivity.
To calculate the distance thresholds upon which to base the connectivity indices for the landscapes, IIC and CCP indices were calculated at 10 m intervals for distances between 10 and 1500 m. This was done for the dense Caatinga class because they provide better soil protection and have a greater ecological role than open Caatinga [63], and because these patches better represent the original native vegetation within the studied region, where open Caatinga is often degraded as a result of agricultural expansion, clear-cuts and overgrazing [33,34]. Based on these preliminary results that are in agreement with some species of mid-sized animals found in the area, the connectivity distances used to analyze the dense Caatinga were set at 100, 560, 1140 and 20,000 m (see Section 3.3). These (100-1140 m) cover a wide range of different species' responses to landscape patterning [64][65][66], including the following: Tropidurus hispidus, which is a reptile that can travel distances between patches nearly 100 m apart [67]; Hydrochoerus hydrochaeris, which is a mammal that has been reported to travel distances of approximately 500 m [68]; and Procyon cancrivorus, which is a mammal that can travel distances around 1900 m [69]. The 20,000 m distance was chosen to assess the landscape perception by Herpailurus yagouaroundi, a feline that travels long distances and is one of the most dispersive terrestrial animals found in the study area [70], inhabiting several areas of the Caatinga biome [71,72]; it has been reported in the fragmented forest in the state of São Paulo [64]. All these indices were calculated with the Graphab software version 2.4 [73].

Assessment of Hypothetical Agriculture Area Growth Projections
In order to assess the effect of the conversion of natural habitat into agriculture (irrigated croplands) on the functional connectivity to target taxa (fauna) in the study range, we established three hypothetical projections of agriculture class expansion based on LULC map classification. We created three new LULC maps in which patches of the 'agriculture' land class existing in 2018 were expanded by a buffer of new agricultural land, 100, 200 or 300 m in width. The first projection, in which all agricultural patches increased by a buffer of 100 m, is equivalent to an agricultural expansion of 5673 ha (6.2% more agricultural land than in 2018); in the second projection, agriculture expanded by a buffer of 200 m equivalent to 10,061 ha (+10.9%); in the third projection, agriculture expanded by a buffer of 300 m, equivalent to 13,352 ha (+14.5%). These projections are consistent with the observed increases in agriculture within the study area between 1985 and 2018 (see Table A5, Appendix B). We used the "Equivalent Connected Area index" (ECA) metric [14,74] to measure the landscape connectivity of the three projections, with and without increases in agricultural land. ECA is the size (area) of a single patch that can provide the same amount of connectivity as observed in the landscape of interest [74]. The formula is as follows: For a landscape with n habitat patches, a is the area extension of a habitat patch. When i = j, a i is patch area i and a j is patch area j. P is the probability of connectivity between the two patches. P ij = 1 when the two patches are connected with a link (i.e., they are closer together than the predefined distance thresholds), and P = 0 if otherwise. When i = j, i and j are the same patch, so P ij = 1. Therefore, ECA takes into account the connected area that exists within a habitat patch. ECA generally increases if (1) the size of a habitat patch increases, (2) a new habitat patch is created, or (3) a new forest corridor is created between habitat patches that were previously isolated [74].
ECA was used to determine to what extent only dense Caatinga patches contribute to overall landscape connectivity. When there is a change in the landscape from 2018, such as an increase or decrease in habitat patches-which most likely also implies a variation in the spatial layout of the habitat within the landscape-we can calculate the relative change in the amount of attainable habitat. This is determined by dECA = (finalECA − initialECA)/initialECA [74].
To analyze the impact of the loss of dense Caatinga patches near agricultural land, we calculated a dECA series comparing the landscape with no increase in the agricultural area (initialECA) to landscapes of the three agriculture expansion projections (finalECA). dECA was calculated for the four established distances described in 2.4 (100, 560, 1140 and 20,000 m). These indices were calculated with the Graphab software version 2.4 [73].

Prioritization of Caatinga Patches to Conserve and Restore
Dense Caatinga patches have a higher conservation priority than open Caatinga patches within our study area, characterized by rapid expansion and intensification of agriculture, as they provide higher ecosystem services and protection against anthropogenic disturbance and aridity for functional diversity [75]. In fact, a recent study from Brazil found the loss of forests across agricultural landscapes to have a negative effect on agricultural productivity by reducing crop pollination [76]. As a result, within our study region, conservation priority was given to dense Caatinga patches, while open Caatinga patches were prioritized for restoration, as they have a reduced capacity for natural regeneration and contain grasses and abandoned areas [44]. This analysis considered the existing patches based on the 2018 land-cover map.
The importance of each patch in terms of its individual contribution to the maintenance of overall connectivity was measured, using the "Delta IIC" index (dIIC), which is usually represented by the removal of individual patches from the landscape and the quantification of the corresponding decrease in connectivity [59,77]. dIIC is calculated as follows: Several studies have highlighted the importance of patches with high network centrality in maintaining and enhancing connectivity [19,59,60]. We used the Betweenness Centrality (BC) index, which assigns a value to each patch and suggests that high-scoring patches constitute the "backbone of the landscape," as a disproportionately high number of the shortest paths in the entire landscape pass through these patches [77]. BC is calculated as follows [59]: BC quantifies the quantity of a patch and stands in between other patches in a network. It calculates the sum of all shortest paths between all pairs that go through the patch k divided by the sum of all the shortest paths between all pairs of patches. These indices were calculated with the Graphab software version 2.4 [73].
To determine the patches of dense Caatinga that should be conserved, the dIIC and BC indices were calculated with the dispersal distance of H. yagouaroundi (20 km; [70]), which was captured with camera traps in the study area (unpublished data). Then, the index values were normalized between 0 and 1 to facilitate their integration, where 0 is the lowest possible score and 1 is the highest possible one.
Vegetation patch conservation studies based on functional connectivity should include the intensity of human disturbance [21]. As a supplement to the aforementioned indices, we included the "Chronic Anthropogenic Disturbance Index" (CDI) [35]. This index measures the intensity of human disturbance and depends on several man-made environmental vectors, such as human population, infrastructure, grazing, timber extraction and wildfires; here, we consider the recent increases in overgrazing by livestock and the associated clear-cutting to promote herbaceous growth, extensive timber extraction and increased fire risk associated with human activities (see Introduction section). CDI values range from 0 to 1, where 0 is the least disturbed pixel and 1 is the most disturbed pixel [35]. The CDI is in raster format with a spatial resolution of 90 m and may be freely downloaded (https://bit.ly/2XEd7v1, accessed on 20 April 2021). In our analysis, higher CDI scores indicate a higher priority for conservation/restoration. We used dIIC, BC and CDI indices simultaneously to select conservation priority patches. To this end, we conducted a multi-criteria decision analysis (MCDA), which provides a systematic methodology for combining the values of these indices in order to define conservation priorities [78]. To evaluate the different results that we might obtain if we favored connectivity or the magnitude of human disturbance, we created two scenarios with different weights for each index and calculated the weighted linear sum of the indices in each scenario [78]. In the first scenario, the weights were 0.6, 0.2 and 0.2 for CDI, dIIC and BC, respectively, and in the second scenario, the weights were 0.2, 0.4 and 0.4 for CDI, dIIC and BC, respectively. In this sense, the first scenario considers species sensitive to human disturbances (more weight to CDI) and the second scenario considers species tolerant to the human matrix (less weight to CDI). The prioritization of open Caatinga patches for restoration followed the same methodology. Higher MCDA scores indicate a higher priority for conservation/restoration.

Spatial-Temporal Change in Land Use and Land Cover
Over the entire study period (1985-2018), agricultural land increased by 8.2%, at a rate of 2104 ha/year, and open and dense Caatinga forests, combined, decreased by 20.2%, at a rate of −5203 ha/year. The net change in area (ha) seen between 1985 and 2018 can be found in Appendix B (Table A4). Image classification had an overall accuracy of 87%, 88% and 90% for the years 1985, 2000 and 2018, respectively. The Kappa coefficient was 0.85, 0.86 and 0.88 for each year, respectively, indicating high accuracy according to this statistical index [79]. The confusion matrices for each period can be found in Appendix A. The results showed that the landscape of the study area has undergone continuous temporal and spatial changes that have transformed its structure and configuration ( Figure 2). During the first period (1985)(1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000), the most notable changes were a 6% (50,738

Determining the Connectivity Distance Threshold
For dense Caatinga, it was found that the distances showing the greatest change in IIC and CCP values were 100, 560 and 1140 m (Figure 5a,b). These distances were selected to conduct an ex-post analysis, as they cover a wide range of responses by different species to the landscape pattern [64][65][66].

Changes in Connectivity under Projections of Agricultural Expansion
In the hypothetical projection, where agricultural patches increased by a buffer of 100 m, a loss in connectivity of approximately 2% was observed for dispersal distances of 100, 560 and 1140 m ( Figure 6). In contrast, connectivity for the dispersal distance of 20,000 m was reduced by 5.8% ( Figure 6). When agricultural land increased by a 200 m buffer, the four smaller dispersal distances experienced a loss in connectivity of approximately 4% ( Figure 6) and the 20,000 m dispersal distance reduced its connectivity by 10.8% ( Figure 6). In the last projection (300 m buffer), the loss of connectivity ranged from approximately 5% to 6% for the 100, 560 and 1140 m dispersal distances (Figure 6), whilst the connectivity reduced by 15% for the 20,000 m dispersal distance ( Figure 6).

Identified Patches for Conservation and Restoration
The results of the prioritization to conserve patches of dense Caatinga for the selected distances showed that under both evaluated scenarios, the patches with the highest priority for conservation were concentrated in Juazeiro (9.4 • S, 40.1 • W) (Figure 7). In the first scenario, 3% (1849 ha) of the patches had conservation priority values greater than 0.6, compared to 11.9% (7362 ha) of patches in the second (Figure 1, Appendix C). The results of the prioritization to restore patches of open Caatinga found, under both scenarios, that the highest number of patches to be restored was located around Petrolina (9.1 • S, 40.7 • W) (Figure 7). In the first scenario, 3.4% (12,210 ha) of the patches had restoration priority values greater than 0.6, compared to 10.8% (38,471 ha) of the patches in the second scenario ( Figure 1, Appendix C).

Discussion
Our study showed a steady increase in agricultural land in the studied landscape, with Caatinga's native vegetation being the most negatively affected. This change is the result of landscape fragmentation, where areas of dense Caatinga have become fragmented (85.2% increase in the number of patches) and decreased in average size (by 84.8%, from 42.2 ha in 1985 to 6.4 ha in 2018). The LULC change patterns found in this study are similar to those reported by other authors for the same biome [7,[80][81][82] and similar to other seasonally dry areas, such as the Mediterranean landscapes in Turkey [83], semi-arid zones in Chile [56], and native forests in Mexico [84]. It is important to consider that open Caatinga patches within the studied area have large tracts of disturbed edges, which have liken arisen due to chronic anthropogenic disturbance through agricultural expansion, overgrazing or harvesting timber [33,34], and are likely to be at further increased risk of habitat conversion to pasture and agriculture [35].
By selecting a range of dispersal distances (from 100 m to 20,000 m), we showed that connectivity indices varied according to species' dispersal abilities [21,59,85]. The distances we calculated empirically, based on the behavior of graph-based connectivity indices for our landscape, were similar to the published dispersal distances of animals which are found in the study area. The species present in our study area have different habitat requirements, which could vary between the rainy and dry seasons, as well as with the distance to urban areas and with edge quality, etc. [86]. Even so, we believe that our calculated connectivity indices represent meaningful scales for the medium-sized mammals and reptiles that have been recorded in the study landscape.
The three projections of agricultural expansion (100,200 and 300 m increase in area around 2018 agriculture) had varying impacts on landscape connectivity and are likely to af-fect species differently, depending on their dispersal ability ( Figure 6). For H. yagouaroundi, the connectivity reduction was approximately 6-15% under our agricultural class expansion projections, as its long-term persistence and viability requires large, interconnected habitats [72]. A landscape fragmentation analysis conducted in the southern Amazonia showed that the minimum vegetation patch size required by H. yagouaroundi is about 790 ha [70]. Patches of this size are scarce in our study area, further decreasing the landscape connectivity for this feline, meaning that maintaining and increasing patch connectivity and conserving the suitable habitat will be key for its persistence.
Decreases in connectivity in response to agricultural expansion were assumed to be smaller for species with shorter dispersal distances, which were based on the published literature, and our findings found reductions in connectivity for these species to be between approximately 2% and 6% ( Figure 6). These decreases in connectivity may have negative consequences for T. hispidus due to its low dispersal capacity [67]. Defining priority native vegetation patches for conservation and restoration, using the dispersal distance of H. yagouaroundi as a reference, could potentially provide protection for other coexisting species, as we expect H. yagouaroundii's requirements to be the largest [87]. These results can be used to make spatial and territorial planning decisions [7], where, depending on the different scenarios, it can be decided whether to favor increasing connectivity, or minimizing potential anthropogenic threats, or both.
We used MCDA to prioritize habitat patches for the conservation of dense Caatinga or the restoration of open Caatinga, using two scenarios, which differed according to the relative weighting assigned to each index of connectivity or disturbance at the scale of habitat patches (Figure 7). In the first scenario, which gave a higher weighting to the CDI, patches with the highest prioritization values were immersed in the agricultural matrix and were very large, so they were the most exposed to agricultural management [64]. Vegetation patches close to the agricultural matrix may, in some cases, provide resources for certain species [88], such as H. yagouaroundi, which has been observed in agricultural areas [89] where there is a high abundance of prey [72]. In the second scenario, most of the patches had low prioritization values since the contribution to connectivity of most of the patches was low. Nevertheless, the first scenario suggested a higher number of patches as being important for restoration, indicating that restoration should happen mostly within open Caatinga patches, which are characterized by being degraded by agricultural activities, high grazing pressure and clear cuts within our studied area. These degraded open Caatinga patches around agricultural areas are large and widely distributed throughout the study area [34], thus, restoring them to promote the growth of native herbaceous species and higher vegetation coverage and density could prove beneficial for biodiversity [7], as large patches can conserve species with large populations [90], maintain a higher number of species [64] and provide a high-quality interior habitat [91]. Furthermore, as the mean Euclidean distance between open Caatinga patches was smaller than that of dense Caatinga patches (106.6 m vs. 148.9 m), restoring these patches could increase landscape connectivity [7]. In agreement with our results, recent studies showed that disturbed Caatinga vegetation fragments, such as the open Caatinga within our study area, still maintain a reasonable structural integrity but have been heavily exposed to anthropogenic pressure in recent decades [26,35], reinforcing the need for landscape planning for conservation and restoration [34].
Agricultural expansion made possible by irrigation, and increases in areas grazed as pastures, both of which are linked to increased urbanization, pose a constant threat to the Caatinga biome. Even though Brazil establishes a regulatory framework for agricultural expansion [18,92], the intensification of agriculture in both size and agrochemical use adversely impacts biodiversity [93]. In addition, sheep-based subsistence livestock farming in which animals are allowed to breed freely and feed on native vegetation are degrading the Caatinga and weakening ecosystem resilience [33,94]. These factors could account for the increase in pastures and the decrease in native vegetation observed in the study area in the São Francisco Valley. Landscape planning policies that include society and its local context are key to achieving sustainable development in the region [95].
Our results should be interpreted considering their limitations. One possible criticism of our analysis is that we assume that open Caatinga, within the study area, is open because of degradation primarily caused by contemporary overgrazing. An alternative perspective is that open Caatinga, maintained by herbivory, is a 'natural' state. This alternative perspective is supported by the observation that some Caatinga species assemblies contain many species, which are shared with other South American dry open forests and savannahs, such as the Argentinian Chaco [28,96]. However, despite such biogeographic considerations, we point to studies that relate contemporary human activities with Caatinga degradation in areas surrounded by agriculture [26,33] and expansion of invasive species associated with livestock pressures [93]. It is, therefore, important to note that there is no evidence that suggests that a 'natural', open Caatinga would biologically resemble the degraded conditions that are common around farmland and other human-disturbed areas. Nevertheless, it is important to caveat our conclusions that in some settings, and for some species groups, it may be preferable to restore degraded Caatinga to an open-structured condition, and that in those circumstances, herbivory may play an important role.
The selected studied area only encompasses one large area of Caatinga, which is characterized by strong agricultural pressure, thus our results may not be applicable to the whole of Caatinga and may have been different had we studied the whole biome. Our results should be taken as showing key patterns and trends, as the exact proportions of land changes may be associated with errors. For example, in the early 1980s, agricultural expansion rates were high, standing at 10,500 ha/year across the entire São Francisco Valley [97], while our results from the first study period show an increase of 3382 ha/year. Another potential limitation is that the land cover categories used in this study are based on geographic information systems (www.mapbiomas.org, accessed on 20 April 2021) and some categories, such as pastures, could include native herbs as part of its community, but the scale of the study does not allow to differentiate this level of detail. Thus, despite these limitations, we present trends that coincide with those from other studies [33,94] and that contribute to the knowledge of this area by considering both the historical and connectivity analysis for spatial and territorial planning. We believe that our results can inform the work of restoration ecologists and policymakers who can further build on these results to develop local and specific conservation goals for sustainable agriculture.
Our results show that vegetation fragments close to agricultural land play a fundamental role in landscape connectivity, acting as stepping-stones and increasing connectivity at the landscape level [64]. The restoration and conservation priorities proposed in this study have the long-term potential to be used as a conservation strategy in the Caatinga [34]. In this way, other studies have shown that the quality of the overall agricultural landscape can be enhanced for the species that live in the remnant habitat, patched via crop diversification and through increased connectivity, indicating a complementary strategy for the São Francisco Valley [17,18,98].
It has been suggested that sustainable agricultural landscapes in dry tropical regions, such as this study area, should have at least 50% natural vegetation to deliver an optimum mix of essential ecosystem services [99]. Although our study area still has a large amount of native vegetation (52.8%) remaining, we have identified historical trends in LULC change consistent with increases in livestock farming, which other studies have shown to have the potential to degrade even large remnant patches of the Caatinga [26,35,100]. If the native vegetation in the study area continues to decline, the degradation of this semi-arid ecosystem could become more severe [7,26,99].

Conclusions
The unique biodiversity of the Caatinga biome requires urgent actions to promote the conservation and restoration of the remnants of the degraded dry forest. We showed that the recent change in LULC between 1985 and 2018 caused habitat loss and fragmentation in the studied Caatinga area in the São Francisco Valley in the Northeast of Brazil. Irrigated agriculture, through indirect conversion (conversion from dense Caatinga to open and then to agricultural land) and direct conversion (Caatinga to agricultural land) are the main causes of the recent land-use change. Fragmentation and loss of connectivity impact species differently according to their dispersal distances and may affect ecosystem functionality. Our analyses showed that, despite the fragmentation of the studied landscape, patches can still be functional for connectivity purposes. We explicitly identified zones that support connectivity for conservation and restoration planning in the area. The results support conservation strategies that aim to increase landscape-scale connectivity for biodiversity and highlight the importance of the Caatinga remnants in the study area's expanding agricultural matrix. The recent agricultural development in the São Francisco Valley needs to achieve sustainable landscape management. Our study provides quantified information that should be used to inform long-term goal planning and to co-produce integrated farm interventions.

Conflicts of Interest:
The authors declare no conflict of interest.