Next Article in Journal
Rural Ecosystem Health Assessment and Spatial Divergence—A Case Study of Rural Areas around Qinling Mountain, Shaanxi Province, China
Next Article in Special Issue
Selecting the Most Sustainable Phosphorus Adsorbent for Lake Restoration: Effects on the Photosynthetic Activity of Chlorella sp.
Previous Article in Journal
Hydrogel Applications in Nitrogen and Phosphorus Compounds Recovery from Water and Wastewater: An Overview
Previous Article in Special Issue
Distribution Pattern of Species Richness of Endemic Genera in Mountainous Areas of Southwest China and Its Influencing Factors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characteristics and Influencing Factors of Landscape Pattern Gradient Transformation of Small-Scale Agroforestry Patches in Mountain Cities

by
Canhui Cheng
1,
Zhong Xing
1,*,
Lin Ye
1,
Junyue Yang
2 and
Zhuoming Xie
1
1
College of Architecture and Urban Planning, Chongqing University, Chongqing 400044, China
2
College of Architecture and Urban Planning, Guizhou University, Guiyang 550025, China
*
Author to whom correspondence should be addressed.
Sustainability 2024, 16(15), 6322; https://doi.org/10.3390/su16156322
Submission received: 18 May 2024 / Revised: 4 July 2024 / Accepted: 15 July 2024 / Published: 24 July 2024
(This article belongs to the Special Issue Biodiversity, Biologic Conservation and Ecological Sustainability)

Abstract

:
Small-scale agroforestry patches possess irreplaceable value compared to large-scale patches. In southwestern mountainous cities of China, the complex terrain and urbanization have led to the presence of numerous small, fragmented agroforestry patches around urban areas. These patches serve as crucial habitats for endemic species and provide essential space for wild food sources, thereby contributing to a range of ecosystem services. Consequently, their proper conservation and utilization planning are of paramount importance. This study investigates the transformation characteristics of landscape patterns of mountainous small-scale agroforestry patches and their constituent elements across urban–rural gradients, identifying the driving factors behind these transformations to support conservation and utilization planning. From an urban–rural gradient perspective, four directional transects were selected and divided into uniform sample grids. Using Fragstats 4.3, landscape indices of small-scale agroforestry patches were calculated, analyzing the transformation characteristics of these patches and their elements across different gradients. Spearman correlation coefficients in SPSS were employed to assess the influence of terrain and relevant anthropogenic factors on the transformation of agroforestry patches. The findings reveal the following: (1) Small-scale agroforestry patches and their elements exhibit similar patterns in terms of size, fragmentation, dispersion, and connectivity, showing an “increasing trend in size and connectivity, decreasing fragmentation, and fluctuating dispersion” from urban centers to natural areas, with slight variations in orchard patches. However, patch cohesion and shape complexity display nonlinear differentiated transformation characteristics. (2) Overall, small-scale agroforestry patches are significantly influenced by anthropogenic construction factors, with the landscape pattern of forest patches notably affected by terrain factors. (3) Across urban–rural gradient zones, the landscape patterns of small-scale agroforestry patches in urban centers, suburbs, and rural natural areas are more affected by terrain factors, whereas those in urban construction zones are significantly influenced by anthropogenic construction factors. The findings of this study provide a scientific basis for the conservation and planning of mountainous small-scale agroforestry patches.

1. Introduction

Small-scale agroforestry patches (SAPs) have been recognized as an important resource type in countries such as those in Europe since the early 21st century [1] and have been designated as key conservation targets [2]. These patches, also referred to as small-scale high nature value farmlands and home gardens, have seen a growing body of research since the 1980s. A search for the terms “High Value Farmland” and “home garden” in the Web of Science and Scopus databases reveals that the number of scientific papers published annually on this topic did not exceed nine until 2006. By 2015, this number had increased to 19 papers per year. France published the first peer-reviewed article on SAPs in 2006, while Italy and Portugal began publishing in this field in 2010 and 2012, respectively. A total of 25 European countries have published peer-reviewed articles focusing on themes such as biodiversity, environmental science, and agricultural conservation policies related to SAPs [3]. Furthermore, research has demonstrated that SAPs possess significant species diversity value [4,5,6], habitat ecological maintenance value [7,8,9], wild food production value [10,11], diversified food supply value [12], medicinal plant supply [13,14,15], landscape recreation [16,17], preservation of traditional horticultural practices [18], and green infrastructure value [19]. In addition, the mountainous cities of southwestern China, due to their complex topography and climatic conditions, boast high forest coverage, creating environments rich in biodiversity [20,21]. The landscape resource value and species diversity value of SAPs in these mountainous areas are particularly notable [22,23]. In conclusion, SAPs in mountainous cities hold significant conservation value and present substantial research potential. However, rapid urbanization in China has significantly changed the spatial structure of urban and rural areas. Construction land has expanded uncontrollably, encroaching on surrounding natural resources such as farmland, forest land, grassland, and wetlands [24,25,26]. This has caused a series of environmental problems [27]. Rural settlements and SAPs, already fragmented by mountainous terrain and human activities [28], have become even more fragmented. Urban construction has further divided them with roads, public service facilities, and buildings, leading to their gradual erosion [29].
In China’s land-use classification, SAPs refer to small, scattered agricultural and forestry lands not designated as protected areas. These patches are fragmented and dispersed due to complex mountainous terrain and human factors. According to national and local laws in China, farmland, orchards, and forest land under 3.3 hectares are not strictly regulated. Only forest land above 3.3 hectares and contiguous farmland and orchards are subject to clear regulatory requirements. In Chongqing, farmland and orchards above 3.3 hectares are designated as high-standard farmland, while smaller plots have no clear regulatory requirements and their occupation is not penalized. This makes them susceptible to encroachment by construction land in urban planning. Therefore, this study focuses on SAPs under 3.3 hectares that are not within protected areas. These patches are abundant in mountainous regions, and their total area is not less than that of protected agroforestry lands. Due to a lack of official recognition of their value, these patches lack protection and are gradually eroded and disappear. Research on these unique mountain agroforestry resources is crucial for urban development planning in mountainous cities.
To prevent further degradation due to urban expansion and effectively safeguard small agricultural and forest patches, understanding the transformation dynamics of small agricultural and forestry patches (SAPs) along the urban–rural gradient is crucial for devising scientifically sound protection and planning strategies. The concept of the urban-to-rural gradient, introduced by McDonnel in 2008, examines the ecological processes in complex urban–rural spatial environments and the impact of urbanization on such spaces. McDonnel extensively investigated urban–rural gradients from downtown New York to rural Litchfield County, Connecticut (along a 140 × 20 km research transect), unveiling the characteristics and causes of forest ecosystem structure, soil composition, climate, impermeability, land cover, and cultural composition along the gradient. This research aimed to develop a rational conservation strategy for natural resources along this trajectory [30]. Additionally, various studies have explored woodlands [31,32], cultivated land [33], insects [34], birds [35,36], changes in heat island intensity [37], water quality transformation [38], and landscape indices transformation [39] based on urban–rural gradients. Analyzing the changing characteristics of each factor along the urban–rural gradient provides a scientific foundation for planning and adaptive strategies.
Landscape pattern indices are metrics used to quantify and describe the spatial structure and distribution of landscapes. By analyzing the size, shape, number, and arrangement of various land use types, these indices reveal landscape complexity, heterogeneity, and ecological processes [40]. They can be used to assess the health of landscape ecosystems. For example, the patch area, shape index, and connectivity index can evaluate habitat quality and biodiversity [41]. They also indirectly represent the ecosystem service capacity of landscape elements [42]. Calculating landscape indices helps to understand the spatial distribution and interactions of ecological processes within landscapes [43]. Connectivity is crucial for species migration and gene flow, while heterogeneity affects ecosystem stability and productivity. Additionally, landscape pattern indices provide scientific support for land-use planning and ecological management. Understanding the current landscape structure allows for the development of more effective conservation and management strategies [44]. Similarly, calculating landscape indices can reflect the current spatial pattern characteristics of SAPs in different gradient zones. This information is useful for formulating targeted planning strategies. Currently, there are many types of landscape indices [45]. For analyzing the transformation characteristics of SAPs in urban–rural gradient zones, indices such as the patch area, patch density, edge density (dispersion index, landscape division), aggregation index (landscape connectivity index), fractal dimension, and largest patch index (LPI) are commonly used [46]. By calculating these indices, we can comprehensively analyze the current structural characteristics of the landscape.
The transformation characteristics of landscape Indices for SAPs and their elements (farmland, orchards, forest land) can represent the current landscape structure along the urban–rural gradient. The spatial distribution of SAPs and each agricultural and forestry element along this gradient determines the direction of conservation strategies. For example, by calculating patch density in different zones, we can understand the fragmentation of farmland patches. Based on the current state of fragmentation, we can predict further fragmentation trends and develop control measures to prevent the further fragmentation of farmland in specific zones. For already fragmented areas, we need to increase farmland connectivity [47]. Furthermore, understanding the specific factors driving the transformation of SAPs can inform the development of targeted management measures for different zones. Scholars generally believe that both natural and human factors influence these transformations [48]. Natural factors include topography, hurricanes, and geological disturbances [49]. The current landscape pattern of SAPs in mountainous areas is greatly influenced by complex topography. Different elevations, slopes, and terrain undulations lead to significant differences in the distribution of these patches [50]. Human factors affecting these patches include construction intensity, road size, and density, as well as urbanization and economic activities [51]. Understanding the impact of these driving factors on SAPs helps in formulating more rational and differentiated management rules for SAPs in various zones.
To sum up, this paper includes three research objectives:
(1) What are the transformation characteristics of landscape patterns for SAPs and their elements (farmland, orchards, forest land) along the urban–rural gradient in mountainous city planning areas?
(2) What are the overall driving factors for SAPs and their elements (farmland, orchards, forest land) along the urban–rural gradient within the planning area of mountainous cities?
(3) What are the driving factors for SAPs and their elements (farmland, orchards, forest land) in different zones along the urban–rural gradient within the planning area of mountainous cities?

2. Materials and Methods

2.1. Study Area

For the present case study, we selected Yongchuan District in Chongqing City, located in Southwest China. Chongqing stands out as one of the most developed cities in China, serving as a core city within the Cheng-Yu economic circle. Characterized by diverse topography and challenging landforms, Chongqing is a representative mountainous city, posing difficulties in construction [52]. By visiting the districts of Chongqing, it is found that there are a lot of SAPs existing in the urban planning area, which form the special landscape features of mountainous areas and are the potential resource elements of mountainous areas, as a part of urban and rural green infrastructure, these can play social, economic, ecological and other composite functions (Figure 1). These fragmented agricultural and forestry lands show different morphological characteristics along the urban–rural gradients, especially in the urban planning area/affected area, and these SAPs will gradually disappear as the city expands, planning areas are therefore considered as study areas. This led to the selection of the Chongqing Yongchuan Urban Planning Area as a case study (Figure 2), an area of 175 square kilometers.

2.2. Research Framework

The study is divided into four steps: first, through the introduction part of the status of Agriculture and forestry in the Southwest Mountain region of policy analysis, we selected 3.3 hectares of unprotected agricultural and forestry land. Secondly, the samples were randomly selected from four different directions from rural to the urban center, and divided into an equal-distance grid of sample points. Thirdly, through the analysis of the data, the agricultural and forest patch and each element landscape indices were calculated by Fragstats 4.3. At the same time, the gradient change in the related factors was calculated; SPSS 19.0 Spearman correlation analysis was used to study the relationship between small-scale agroforestry patches and their influencing factors, and to determine which factors are the main driving forces to promote the segmentation, fragmentation, and even disappearance of agroforestry patches. Through the above research process, on the one hand, we can grasp the change characteristics of small-scale agricultural and forest patches and provide a scientific basis for the protection and utilization of small-scale agricultural and forest patches in urban planning areas; on the other hand, we know the factors that affect the characteristics of the patch pattern, which can provide the basis for the management and control of small-scale patches.

2.3. Study Object Selection

In the first chapter, the object of study has been clearly defined and determined in accordance with laws and regulations and the requirements for the protection of agricultural and forestry land, referring to small-scale agroforestry land with an area of less than 3.3 hectares. Therefore, this part of the land use should be statistically screened (Figure 3).
Through GIS, the space of agriculture and forestry under 3.3 ha is 8694.27 ha, accounting for 75.60% of the total area, and the space of agriculture and forestry over 3.3 ha is 2805.77 ha, accounting for 24.40% of the area. Overall, 5157.88 hectares (73.4%), 1268.29 hectares (91.1%), and 2268.13 hectares (73.5%) are under 3.3 hectares of arable land, garden land, and woodland, respectively (Table 1). And through the analysis of land-use characteristics, we attain an SAP spatial gradient of the transformation of the status analysis.

2.4. Sample Selection

To accurately depict the landscape pattern characteristics of SAPs across diverse gradients, a spatially informed approach was adopted. Commencing from the city center area, defined as the starting point of the city’s administrative office space, the investigation extended in four directions within the planning area. Samples were selected with dimensions of 22.4 km × 4 km along both horizontal and vertical axes. These samples were evenly divided into eight squares, each measuring 2.8 km × 4 km (Figure 4). According to the urban–rural transect zoning theory [53], the samples are divided into an Urban Center Zone, General Urban Zone (T2), Sub-urban Zone (T3), Rural Zone (T4), and Natural Zone (T5). The samples in the Urban Center Zone are 1-1, 2-1, 3-1, and 4-1. The samples in the Urban Zone are 1-2, 2-2, 3-2, and 4-2. The samples in the Sub-urban Zone are 1-3, 2-3, 3-3, and 4-3. The samples in the Rural Zone are 1-4, 2-4, and 4-4. The sample in the Natural Zone is 1–5.

2.5. Data Analysis

2.5.1. Calculation of SAP and Agroforestry Elements Landscape Indices

SAPs are classified according to Chinese land-use standards; agricultural and forestry land can be divided into arable land (dryland, paddy field, irrigated land), garden land (orchard, tea garden land, oak garden land, other garden land), and woodland (arbor woodland, bamboo woodland, shrub woodland, other woodland)) (Table 2). The change in SAP landscape indices and the change in various factors were calculated to reflect their transformation characteristics.
Landscape indices can not only simplify the landscape information, but also reflect the structure and spatial configuration of landscape information. They represent a simple and applicable quantitative index. In this study, Fragstats 4.3 software was used to analyze the spatial landscape pattern of agriculture and forestry in the study area. The study divided the agroforestry space into 3 Mm × 3 m grids and calculated the landscape pattern index of cultivated land, garden land, and woodland within each 2.8 km × 4 km plot. Three levels of analysis were included: patch, class, and landscape. For the purposes of calculation, the parameters of landscape type, depth effect distance of landscape edge, difference weight of adjacent landscape type, and similarity of adjacent landscape type were set. Within this context, the term “landscape type” encompasses arable land, garden land, and woodland. The “landscape edge depth effect distance” denotes the extent of the influence exerted by the edge effect of one patch on adjacent patches. The “difference weight of adjacent landscape type” signifies the extent of dissimilarity between two neighboring patches, while the “similarity of adjacent landscape type” quantifies the degree of resemblance between two adjacent patches.
Utilizing the tripartite data, Fragstats 4.3 was employed to compute landscape indices characteristics for agroforestry spaces, yielding the small agricultural and forestry patches (SAPs) landscape pattern indices for each segment. The indices encompass the patch area (PA), patch density (PD), edge density (ED), largest patch index (LPI), fractal dimension, and aggregation index (AI). These metrics collectively capture the scale, fragmentation, segmentation (spatial distribution), predominant types of agroforestry elements, shape complexity, and connectivity of agroforestry patches (Table 3).

2.5.2. Statistical Analysis of Natural and Artificial Factors

Analysis of Topographic and Geomorphological Factors

Research shows that natural succession and disturbances have a relatively small impact on the overall changes in agroforestry patches and take a long time to manifest. In mountainous environments, key factors include complex topography such as elevation, slope, and terrain ruggedness. These topographic factors determine the layout of agroforestry spaces [50,54]. The underlying reason is that topography affects soil nutrients, climate, sunlight, and humidity [55,56], which indirectly determine the current distribution of agroforestry patches. Therefore, this paper chooses topography as a natural factor for analysis.
Based on the DEM data of 12.5 m, the elevation, slope, and undulation (Figure 5) were analyzed and calculated by Arcgis 10.2. The slope calculation was obtained by 3D Analyst, relief was calculated by the Spatial Analyst Tool neighborhood analysis focus statistics. Then, the average values of elevation, slope, and fluctuation of each subarea were calculated by the Spatial Analyst Tool regional analysis regional statistics (Table 4).

Artificial Factor Analysis

Artificial factors include the impervious area, road density, building coverage, and so on. The impervious area is actually the superposition of the building’s surface coverage and the road area. In addition, other factors affecting SAP distribution are the construction of existing parks and the excavation of artificial ponds and lakes. First of all, the construction of the park will be combined with the existing distribution of agricultural and forestry resources, choosing to occupy these resources to build the park, so it will affect the transformation of the SAP status quo; at the same time, artificial pit pond lakes are usually used to irrigate cultivated land, garden land, and part of forestland, and their spatial distribution and scale are likely to affect the spatial pattern of SAPs [50]. Therefore, the selection of artificial factors includes building surface coverage, roads, parks, and the artificial excavation of the pit pond lake.
Based on the data of the third survey, the area scale of the road, park, building surface, and artificial pond and lake surface on the urban–rural gradient is calculated (Table 5).

2.5.3. Correlation Analysis

The paper uses Spearman correlation analysis to determine the relationship between SAPs and topographic and human factors. Spearman is a non-parametric test method that assesses the strength of the relationship between two random variables based on the size of the coefficient [57]. By calculating the rank correlation coefficient rs(Xi, Yi) and p-value between SAPs and various factors, we can determine the strength of their association. This helps to identify which factors drive the changes in the landscape indices of SAPs and elements across urban–rural gradient zones. The mathematical principle is as follows: Let {(Xi, Yi)}, denote n pairs of data pairs with an independent and identical distribution whose parent body is a binary continuous distribution. To arrange Xi from small to large, we form a new set of data pairs X (1) < X (2) < X (3) … < X (N). Where), where the Yi corresponding to the order statistic of X is called an adjoint of Xi. The formula is as follows:
r s X i , Y i = 1 6 1 n ( P i Q i ) 2 n n 2 1
In the formula, P i means Xi is the K position in the sequence {( Xi)}; K is the rank of Xi; similarly, Q i is the rank of Yi. The positive (negative) correlation coefficient indicated that there was an influence relationship among the factors. Rs < 0.7 was a low correlation, 0.7–0.9 was a moderate correlation, and > 0.9 was a high correlation. If p > 0.05, the two functions were independent. The effect of each factor on patch transformation was estimated by the Spearman coefficient. The paper uses SPSS 19.0 to calculate the Spearman coefficients (rs and p-value) between the landscape indices of SAPs and various topographic and human factors, which helps to identify which factors drive the transformation of the landscape patterns of SAPs.

3. Results and Discussion

3.1. The Transformation Characteristics of Agricultural and Forestry Landscape Pattern Gradients

3.1.1. Overall Transformation Characteristics of Agroforestry Patches

Through the calculation of four sample landscape indices, it was found that various landscape indices of agroforestry patches in different directions exhibit non-linear changes (Figure 6). Firstly, the area of SAPs, maximum patch index, and aggregation index all show an overall increasing trend from urban centers to natural areas, with the patch area being largest in rural areas and slightly decreasing thereafter, while patch density decreases linearly. These findings align with natural spatial distribution patterns, where patches tend to become larger, more aggregated, more connected, and less fragmented. The largest patch areas in rural areas are attributed to the concentration of rural residents, who have the highest demand for SAPs.
Secondly, the boundary density of SAPs shows a stepwise change pattern, with similar boundary densities observed between urban centers and rural areas, indicating similar dispersion levels of these patches in these two zones.
Lastly, the fractal dimension of SAPs also displays a stepwise change pattern, reaching the highest values in suburban and rural areas, indicating the highest complexity of patch shapes in these regions. In contrast, the fractal dimension is smallest in urban development zones, indicating more regular shapes of agroforestry patches. Similar results are observed in the fractal dimension of SAPs between urban centers and natural areas.

3.1.2. Transformation Characteristics of Agricultural and Forestry Elements

Transformation Characteristics of Landscape Indices of Woodland

The transformation characteristics of small-scale forest patches’ landscape indices share similarities and differences with those of overall agroforestry patches. Firstly, based on the research findings (Figure 7), the area of forest patches shows a linear increase from urban centers to natural areas, while patch density exhibits the opposite trend. Compared to overall SAPs, the forest patch area reaches its maximum in natural areas, consistent with the pattern of natural ecological spatial changes. The density and edge density transformation characteristics of small-scale forest patches are similar to those of overall agroforestry patches.
Secondly, the maximum patch index, fractal dimension, and aggregation index of small-scale forest patches exhibit their own transformation characteristics. The maximum patch index is highest in urban centers and natural areas, and lowest in rural areas and urban development zones, indicating that small-scale forest patches have a high connectivity in areas with both high and low levels of urbanization. The fractal dimension is largest in urban centers and suburbs, and smallest in natural areas, reflecting the varying complexity of small-scale forest patch shapes in these regions. The aggregation index is lowest in urban development zones and highest in natural areas, indicating the poorer connectivity of small-scale forest patches in urban development areas.

Transformation Characteristics of Patch Landscape Indices of Cultivated Land

The transformation characteristics of small-scale cropland patches share both similarities and differences with overall agroforestry patches (Figure 8). Firstly, the transformation pattern of small-scale cropland patch area is identical to that of overall agroforestry patches, adhering to normal conversion rules. Additionally, the edge density of small-scale cropland patches mirrors that of overall agroforestry and forest patches. The density of small-scale cropland patches also shows a general decline, with similar densities observed in suburban and rural areas. The maximum patch index exhibits an increasing trend, indicating the higher connectivity of cropland patches. In contrast, the fractal dimension of cropland patches displays entirely different transformation characteristics from agroforestry and forest patches, showing a consistent upward trend, reflecting increasingly complex shapes. Similarly, the aggregation index trends upward in alignment with agroforestry and forest patches.

Landscape Indices Transformation Characteristics of Garden Land

The garden land patches differ from other agroforestry elements and the overall agroforestry landscape indices in all aspects except patch area (Figure 9). Overall, the garden land patch area shows a continuous increase in size. The density, edge density, maximum patch index, and aggregation index of garden land patches exhibit a similar piecewise transformation pattern. Specifically, the fragmentation, dispersion, patch size, and connectivity are highest in urban center areas, sharply decrease in urban construction zones, show moderate growth in suburban areas, reach nearly their lowest in rural regions, and then gradually increase again in natural areas. This transformation aligns closely with urbanization patterns and natural features. Additionally, the fractal dimension of garden land patches shows a transformation pattern of initial increase followed by decrease, with the lowest fractal dimension observed in urban center areas and the highest in urban construction zones, gradually declining thereafter.

3.2. Drivers of Gradient Transformation of Agroforestry Landscape Patterns

3.2.1. The Driving Factors of the Overall Transformation of Agroforestry Landscape Pattern

Based on the correlation analysis of landscape indices and natural or anthropogenic factors, it is evident that for SAPs across the urban–rural gradient, terrain factors have a relatively weak influence on agroforestry landscape indices (Figure 10). Firstly, the area of SAPs shows a strong negative correlation with road use, parkland, and building surfaces, indicating a predominant influence of anthropogenic factors on patch size. Patch density correlates negatively with terrain undulation but positively with road use and building surfaces, suggesting that flatter terrains are more susceptible to fragmentation due to construction activities. This is supported by the negative correlation between the aggregation index and the road and building land area. Other landscape indices (edge density, LPI, and fractal dimension) show no significant correlations with terrain or anthropogenic factors.
On the other hand, in the urban center, there is a significant positive correlation between the patch area and the artificially excavated pond area. Additionally, the edge density of SAPs in the urban center is significantly influenced by terrain features, showing a positive correlation. This may be due to the fact that most SAPs in the center are located in terraced land difficult for construction, characterized by a high elevation, slope, and undulation. Moreover, the maximum patch index shows a positive correlation with the road area, indicating that roads may serve as corridors connecting SAPs in the urban center. Other landscape indices in this area show no correlations with influencing factors. In urban construction zones, patch area shows a negative correlation with road and building coverage, indicating the erosion of small agroforestry patches by construction activities. Other landscape indices show no significant correlations with influencing factors. In suburban areas, there is a significant positive correlation between the maximum patch index and slope, indicating that steeper areas tend to have more aggregated and larger SAPs. Outside urban areas, in rural and natural regions, patch density shows a significant negative correlation with terrain features and a significant positive correlation with road use area, suggesting that flatter areas are more susceptible to fragmentation due to infrastructure such as roads (Figure 11).
Therefore, from the correlation analysis of landscape indices and natural and anthropogenic factors of SAPs, it is evident that anthropogenic construction activities predominantly drive the transformation of patch area. In flat areas, the transformation drivers for patch density (dispersion) and the aggregation index (connectivity) are primarily influenced by construction activities, whereas complex terrain environments are less affected. Additionally, there are significant differences in the driving factors of landscape indices across different urban–rural gradient zones. SAPs in urban centers exhibit a certain pattern stability, with less overall impact from anthropogenic construction activities. In contrast, patches in urban construction zones are dynamically changing due to ongoing construction activities. SAPs in steep areas of suburban zones are relatively well-preserved. In urban–rural and natural areas, SAPs in flat areas are fragmented due to road construction activities.

3.2.2. Drivers of Transformation of Each Element Landscape Pattern

Analysis of Driving Factors of Woodland Conversion

The results indicate that, overall, terrain features and anthropogenic factors similarly influence forest landscape indices (Figure 12). Forest patch area shows significant correlations with various factors, being positively correlated with terrain features such as elevation, slope, and undulation. This suggests that small-scale forest patches are predominantly located in areas with poorer construction conditions. Conversely, forest patch area shows a significant negative correlation with road, park, and building areas, indicating substantial anthropogenic erosion. Patch density correlates negatively with elevation and undulation but positively with the road and building land area, highlighting increased fragmentation driven by construction in flatter environments. The aggregation index of forest patches shows a significant positive correlation with terrain features and a negative correlation with building surfaces and pond/lake areas, suggesting increased connectivity with more complex terrain while being hindered by construction activities.
Furthermore, the maximum patch index exhibits significant positive correlations with slope and undulation, and negative correlations with pond/lake areas, indicating a higher patch connectivity in steeper areas, yet it is inhibited by pond/lake construction. In urban centers, forest patch area shows a significant positive correlation with terrain features and a negative correlation with building area, implying larger forest patches in complex terrains where construction is less suitable, although still subject to erosion. Additionally, patch density and the maximum patch index show significant negative correlations with pond/lake areas, which are primarily artificial irrigation ponds, and positive correlations with arable land patch area (as mentioned below), suggesting space occupied by ponds or arable land contributes to forest fragmentation.
In urban construction zones, the forest patch area shows significant positive correlations with elevation, indicating more forest preservation in higher elevation areas. Edge density shows a significant positive correlation with slope and a negative correlation with pond/lake area, indicating natural dispersion driven by slopes hindered by pond/lake construction. In suburban areas, the forest patch area shows positive correlations with slope, edge density with undulation, and the maximum patch index with undulation, suggesting forest preservation in steep areas where the land is less flat; thus, forest patches are more dispersed due to human interference in flat land.
In rural and natural areas, the forest patch area shows significant positive correlations with terrain features and the aggregation index, and negative correlations with road-use area. Patch density shows negative correlations with the pond/lake area, and the maximum patch index shows positive correlations with the pond/lake area, indicating better forest preservation and higher connectivity in complex terrains while roads hinder connectivity. Additionally, pond/lake construction reduces forest fragmentation (Figure 13).

Driver Factor Analysis of Cultivated Land Conversion

From the analysis results, it is evident that overall, arable patch area is significantly influenced by anthropogenic factors (Figure 14). It shows a significant negative correlation with road size, parkland, and built surfaces, while positively correlating with pond/lake areas. Arable patch density correlates negatively with slope and undulation, and positively with road and building land area, similar to forest patches. This indicates that areas with a greater slope and undulation have a lower patch density but a higher connectivity of arable patches, while larger construction areas result in a higher fragmentation of arable land. Additionally, the arable patch edge density and aggregation index show positive correlations with pond/lake areas, suggesting that larger ponds/lakes lead to a more dispersed distribution of arable land without affecting its connectivity.
In urban centers, the distribution of arable patch area and density is significantly influenced by ponds/lakes, with the former showing a positive correlation and the latter a negative correlation with pond/lake size, mirroring the overall influences on arable patch area and density. In urban construction zones, the shape complexity of arable patches correlates negatively with slope and positively with pond/lake area, indicating that higher slopes lead to less complex shapes of arable patches, suggesting more regular field patterns, while ponds/lakes promote greater shape complexity. In suburban areas, arable patch density correlates negatively with undulation, shape complexity positively correlates with undulation, and the aggregation index negatively correlates with slope. Flat areas in suburban regions exhibit a higher fragmentation and higher connectivity of arable patches but a lower complexity of shapes. In rural and natural areas, the arable patch area shows a positive correlation with building area, indicating that larger building areas correspond to larger arable patch areas, reflecting the development of residential land and subsequent expansion of surrounding small arable patches (Figure 15).

Analysis of Driving Factors of Garden Land Conversion

The overall area of garden patches is significantly influenced by human factors, showing a negative correlation with road and park areas, and a positive correlation with pond/lake areas. This reflects that the construction of roads and parks typically occupies existing garden spaces, while agricultural land is not only taken up by roads and parks but also eroded by building construction. Both of these factors are closely associated with ponds/lakes, which require appropriate irrigation, thereby increasing the size of garden and agricultural areas. Additionally, garden patch density correlates negatively with slope and undulation, and positively with the road and building area, similar to patterns observed in forest patches. This suggests that areas with a flat terrain exhibit a higher fragmentation but a higher connectivity in complex terrains, where road and building construction are significant contributors to fragmentation, as evident from the maximum patch index. The fractal dimension also shows a moderate positive correlation with artificial pond/lake surfaces, indicating that these are excavated for irrigation purposes, influencing the complexity of garden patch shapes, with more pond/lake surfaces increasing the likelihood of garden cultivation (Figure 16).
In urban centers, the area, density, and edge density of garden patches correlate positively with the topographical environment but negatively with the building area, while the aggregation index shows a negative correlation with topography and a positive correlation with building area. This indicates that the size, fragmentation, and dispersion of garden patches in central areas are more significantly influenced by natural factors, with less disturbance in complex terrains. However, high fragmentation and dispersion in flat areas do not imply independence among patches but rather a higher connectivity. The results also show that larger fractal dimensions correspond to smaller road areas, highlighting roads as primary drivers of regular garden patch shapes. The situation in urban construction zones regarding arable patch area and density mirrors that of urban centers. In suburban areas, the garden patch fractal dimension correlates negatively with ponds/lakes, suggesting that the pond/lake layout determines the complexity of garden patch shapes. In rural and natural areas, the garden patch fractal dimension correlates positively with the building area, akin to agricultural land in rural natural areas (Figure 17).

4. Discussion

4.1. Landscape Indices Transformation Characteristics and Driving Factors of Small-Scale Agroforestry Patches in Mountainous Urban Areas

Based on the research findings, the landscape indices’ transformations of SAPs in mountainous urban areas exhibit significant differentiation. Overall, patch area shows an increasing trend across urban–rural gradients, while patch density demonstrates a decreasing trend. This transformation pattern reflects the conversion dynamics between agricultural and forestry land uses and natural ecological spaces [58]. However, small-scale garden patches exhibit variation, with the highest patch densities in central areas and the lowest in rural regions. Edge density shows an overall fluctuating trend, indicating varying degrees of patch dispersion or isolation, generally highest in central and rural areas and lowest in natural areas. Small-scale cropland and garden patches also exhibit slight differences. Additionally, the maximum patch index, fractal dimension, and aggregation index of SAPs and their elements demonstrate distinct transformation characteristics. These transformation patterns indicate that SAPs in mountainous urban areas differ from those in plain urban areas [59]. Mountainous SAPs consist of a semi-natural mix of various agroforestry elements, with different types exhibiting varied landscape patterns due to natural environments and urbanization processes. Furthermore, the area of SAPs is significantly influenced by human activities, leading to decreased patch sizes and increased fragmentation, consistent with existing research conclusions [29,60]. Notably, the landscape patterns of small-scale woodland patches are more significantly influenced by natural terrain factors. These findings underscore the need for tailored strategies in planning the conservation and utilization of SAPs in mountainous areas, based on their distinct typological differences.
Based on the analysis of driving factors across urban–rural gradients, it is evident that SAPs and their elements in different zones are significantly influenced by varying factors. In urban centers, suburban areas, and rural natural zones, these patches are primarily influenced by natural terrain factors, which contrasts with general perceptions where urban centers and suburbs are typically more affected by construction processes. SAPs in urban development areas are notably disrupted by construction activities. Moreover, specific influencing factors vary across different zones. Therefore, differentiated zoning in planning and management may be crucial for guiding the conservation and utilization of SAPs in mountainous regions.

4.2. Planning and Management Methods for the Conservation and Utilization of Small-Scale Agroforestry Patches in Mountainous Urban Areas Based on Landscape Pattern Characteristics

4.2.1. Assessment of the Current Status of Small-Scale Agroforestry Patches in Mountainous Areas Should Employ a Differentiated Zoning Assessment Approach

The gradient changes in landscape indices reflect the diverse current structural characteristics of SAPs in mountainous areas across different zones. Moreover, the influencing factors on landscape patterns vary between these zones, suggesting the need for differentiated assessment systems. Research indicates significant variations in the current status, structure, and functions of agroforestry land along urban–rural gradients. This variation arises from differences in residents’ demands and the current state of agroforestry land across different gradient zones, leading to distinct functional values [61]. Studies by Penghui et al [59] further highlight that socio-economic and ecological differences contribute to the spatial and functional differentiation of agroforestry areas along urban–rural gradients.
These findings collectively suggest that a uniform approach is not suitable for assessing small-scale agroforestry in environmentally distinct mountainous regions. For instance, applying criteria from basic farmland protection zones in plain areas to assess SAPs often fails to identify high-value agroforestry land accurately. This study’s results also demonstrate that using indicators such as contiguity and slope criteria from basic farmland delineation standards [62] in assessing SAPs would likely exclude many valuable patches from protection zones. Given their fragmented nature and challenging geographic locations to meet regulatory requirements, many fragmented SAPs with unique species values near urban development areas and high slopes would struggle to qualify under minimum area standards, exacerbating their exclusion from protection zones.

4.2.2. The Conservation of Small-Scale Agroforestry Patches in Mountainous Areas Should Adopt a Zoning-Based Targeted Planning and Management Approach

Landscape indices of SAPs in mountainous areas are notably influenced by both human-made and natural factors, with significant differences observed among the constituent elements of these patches. Specifically, only forest patches are significantly affected by natural factors, whereas overall, human-made construction factors exert a more pronounced influence, consistent with drivers observed in agroforestry patches in Chinese plain cities [63]. However, SAPs and their elements across urban–rural gradients in mountainous cities do not universally share the same driving factors. For instance, in urban center areas, the area of SAPs correlates closely with the size of remaining ponds, while in construction zones, significant impacts are attributed to road and building construction.
Differences in characteristics among zoned agroforestry patches necessitate targeted planning and control strategies based on differential assessments. Specific planning interventions should be selected based on the specific driving factors of agroforestry patch landscape indices in each zone [64]. For example, based on current research, SAPs near urban areas should focus on artificial ecological restoration, utilizing landscape corridors to connect fragmented agroforestry patches into a green infrastructure network. Conversely, SAPs located further from urban centers should be enhanced through protective measures, delineating conservation boundaries and implementing specific control regulations [65].
Furthermore, within the same zone, driving factors influencing the landscape indices transformation of various small-scale agroforestry elements also differ. For example, in urban center areas, the area of forest and garden patches is closely related to terrain and environmental factors, while the area of arable land patches is strongly correlated with pond areas. In urban construction zones, forest patches are significantly influenced by elevation, while garden patches are notably affected by road and building construction disturbances. Therefore, planning and control methods should not only be formulated based on zone differences but also classified according to the types of small-scale agroforestry elements.

4.3. Study Limitations

There are three limitations in this paper: First, although the characteristics of the landscape pattern transformation of SAPs and its elements have been described, there are still some problems in the study, such as the small number of samples and the short gradient distance between urban and rural areas; however,, larger samples and longer urban-rural gradient distances could yield more accurate study results.. Secondly, the selection of influencing factors of SAPs can be more diversified, and the changes in agroforestry landscape pattern may be influenced by other factors; examples include urban development policies and future directions [66], typhoons and earthquakes [49], and artificial and natural fires [67]. Therefore, all kinds of relevant factors should be considered in order to make SAP protection decisions more reasonable. Finally, this paper only explores the transformation characteristics of the SAPs landscape pattern from space. In fact, it is also affected by the time change, as the landscape indices of the same plot will be different at different times [68]. In order to understand the characteristics of SAP transformation, we must study it from the two dimensions of time and space, which is the content to be further improved in the future.

5. Conclusions

This study used Yongchuan, a mountainous urban area, as a case study to investigate the transformation characteristics and influencing factors of SAPs across urban–rural gradient zones. The research findings are as follows:
(1) The scale, fragmentation, dispersion, and connectivity of SAPs and their elements generally exhibit similar patterns. Specifically, from urban centers to natural areas, there is a trend of “fluctuating increases in scale and connectivity, decreasing fragmentation, and fluctuating changes in dispersion”, with garden patches showing slight variations. However, patch cohesion and shape complexity display nonlinear differentiated transformation characteristics.
(2) In general, SAP are significantly influenced by human-made construction factors, with the landscape pattern of forest patches particularly affected by terrain and topographical factors.
(3) Within urban–rural gradient areas, the landscape configurations of SAPs in urban centers, suburbs, and rural natural settings are more significantly influenced by terrain factors, whereas those within urban development zones are distinctly affected by human-made construction factors.
The analysis of transformation characteristics and influencing factors of SAPs in mountainous areas suggests that while most landscape indices exhibit nearly identical transformation characteristics across urban–rural gradient zones, conservation and utilization strategies should be tailored to the differences in zone and element landscape patterns. It is essential to develop targeted strategies rather than applying uniform control methods like those used for basic farmland protection areas.

Author Contributions

Conceptualization, Z.X. (Zhong Xing); methodology, C.C.; software, Z.X. (Zhuoming Xie); validation, L.Y., J.Y. and C.C.; formal analysis, C.C.; investigation, C.C.; resources, Z.X. (Zhong Xing); data curation, C.C.; writing—original draft preparation, C.C.; writing—review and editing, L.Y.; visualization, Z.X. (Zhuoming Xie); supervision, J.Y.; project administration, Z.X. (Zhong Xing); funding acquisition, Z.X. (Zhong Xing). All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by The National Natural Science Foundation of China (grant number: 52178032).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

We appreciate the assistance of the Natural Resources Bureau of Yongchuan District, Chongqing, for providing the base data.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Plieninger, T.; Bieling, C. Resilience-Based Perspectives to Guiding High-Nature-Value Farmland through Socioeconomic Change. Ecol. Soc. 2013, 18, 20. [Google Scholar] [CrossRef]
  2. Ratcliffe, D.A.; Smith, J.E.; Clapham, A.R. Nature conservation: Aims, methods and achievements. Proc. R. Soc. Lond. Ser. B Biol. Sci. 1997, 197, 11–29. [Google Scholar] [CrossRef]
  3. Benedetti, Y. Trends in High Nature Value farmland studies: A systematic review. Eur. J. Ecol. 2017, 3, 19–32. [Google Scholar] [CrossRef]
  4. Álvarez Díaz, J.E.; Santa Regina, M.D.C.; Santa Regina, I. Can the stepping stone enhance the establishment, competition and distribution of sown grassland species during recovery on exarable lands? Open J. Ecol. 2016, 6, 579–597. [Google Scholar] [CrossRef]
  5. Haggar, J.; Pons, D.; Saenz, L.; Vides, M. Contribution of agroforestry systems to sustaining biodiversity in fragmented forest landscapes. Agric. Ecosyst. Environ. 2019, 283, 106567. [Google Scholar] [CrossRef]
  6. Hodel, U.; Gessler, M. In Situ Conservation of Plant Genetic Resources in Home Gardens of Southern Vietnam: A Report of Home Garden Surveys in Southern Vietnam, December 1996–May 1997; (Working Paper); International Plant Genetic Resources Institute: Italy, Rome, 1999; n. 106; ISBN 978-92-9043-419-1/ 92-9043-419-8. [Google Scholar]
  7. Han, L.; Wang, Z.; Wei, M.; Wang, M.; Shi, H.; Ruckstuhl, K.; Yang, W.; Alves, J. Small patches play a critical role in the con-nectivity of the Western Tianshan landscape, Xinjiang, China. Ecol. Indic. 2022, 144, 109542. [Google Scholar] [CrossRef]
  8. Brothers, T.S.; Spingarn, A. Forest Fragmentation and Alien Plant Invasion of Central Indiana Old-Growth Forests. Conserv. Biol. 1992, 6, 91–100. [Google Scholar] [CrossRef]
  9. Körner, C. Mountain Biodiversity, Its Causes and Function. AMBIO J. Hum. Environ. 2004, 33, 11–17. [Google Scholar] [CrossRef]
  10. Fu, Y.; Guo, H.; Chen, A.; Cui, J.; Padoch, C. Relocating plants from swidden fallows to gardens in Southwestern China. Econ. Bot. 2003, 57, 389–402. [Google Scholar] [CrossRef]
  11. Eyzaguirre, P.B.; Linares, O.F. Homegardens and Agro Biodiversity; Smithsonian Press: Washinton, DC, USA, 2004; 296p, ISBN 158834-112-7. [Google Scholar]
  12. Ivanova, T.; Bosseva, Y.; Chervenkov, M.; Dimitrova, D. Enough to Feed Ourselves!—Food Plants in Bulgarian Rural Home Gardens. Plants 2021, 10, 2520. [Google Scholar] [CrossRef]
  13. Abdullah, A.; Khan, S.M.; Pieroni, A.; Haq, A.; Haq, Z.U.; Ahmad, Z.; Sakhi, S.; Hashem, A.; Al-Arjani, A.-B.F.; Alqarawi, A.A.; et al. A Comprehensive Appraisal of the Wild Food Plants and Food System of Tribal Cultures in the Hindu Kush Mountain Range; a Way Forward for Balancing Human Nutrition and Food Security. Sustainability 2021, 13, 5258. [Google Scholar] [CrossRef]
  14. Farina, A. The Cultural Landscape as a Model for the Integration of Ecology and Economics. BioScience 2000, 50, 313–320. [Google Scholar] [CrossRef]
  15. Halada, L.; Evans, D.; Romão, C.; Petersen, J.-E. Which habitats of European importance depend on agricultural practices? Biodivers. Conserv. 2011, 20, 2365–2378. [Google Scholar] [CrossRef]
  16. Byrne, J.; Sipe, N.; Searle, G. Green around the gills? The challenge of density for urban greenspace planning in SEQ. Aust. Plan. 2010, 47, 162–177. [Google Scholar] [CrossRef]
  17. O’Rourke, E.; Kramm, N. Changes in the Management of the Irish Uplands: A Case-Study from the Iveragh Peninsula. Eur. Countrys. 2009, 1, 53–66. [Google Scholar] [CrossRef]
  18. Xu, J.; Ma, E.T.; Tashi, D.; Fu, Y.; Lu, Z.; Melick, D. Integrating Sacred Knowledge for Conservation: Cultures and Landscapes in Southwest China. Ecol. Soc. 2005, 10, 7. [Google Scholar] [CrossRef]
  19. Sanesi, G.; Colangelo, G.; Lafortezza, R.; Calvo, E.; Davies, C. Urban green infrastructure and urban forests: A case study of the Metropolitan Area of Milan. Landsc. Res. 2017, 42, 164–175. [Google Scholar] [CrossRef]
  20. Qian, S.; Qin, D.; Wu, X.; Hu, S.; Hu, L.; Lin, D.; Zhao, L.; Shang, K.; Song, K.; Yang, Y. Urban growth and topographical fac-tors shape patterns of spontaneous plant community diversity in a mountainous city in southwest China. Urban For. Urban Green. 2020, 55, 126814. [Google Scholar] [CrossRef]
  21. Zhang, Y. Advantages and rational utilization of agricultural climatic resources in mountainous areas of China. Mt. Res. 1992, 1, 11–18. [Google Scholar]
  22. Xing, Z.; Tang, X.; Gu, Y.; He, Y.; Lei, Q. Urban-rural shared “Nature-Agriculture Park”: Research on the protection planning method for high-value small-scale agricultural and forestry land in urban-rural fringe areas. Urban Plan. 2021, 45, 54–64. [Google Scholar]
  23. Tang, X.; Xing, Z. Management of small-scale agricultural and forestry land systems in urban fringe areas in the context of spa-tial planning. Planners 2020, 36, 50–57. [Google Scholar]
  24. Del Mar López, T.; Aide, T.M.; Thomlinson, J.R. Urban Expansion and the Loss of Prime Agricultural Lands in Puerto Rico. AMBIO J. Hum. Environ. 2001, 30, 49–54. [Google Scholar] [CrossRef]
  25. Mao, D.; Wang, Z.; Wu, J.; Wu, B.; Zeng, Y.; Song, K.; Yi, K.; Luo, L. China’s wetlands loss to urban expansion. Land Degrad. Dev. 2018, 29, 2644–2657. [Google Scholar] [CrossRef]
  26. van Vliet, J. Direct and indirect loss of natural area from urban expansion. Nat. Sustain. 2019, 2, 755–763. [Google Scholar] [CrossRef]
  27. Bhatta, B. Causes and Consequences of Urban Growth and Sprawl. In Analysis of Urban Growth and Sprawl from Remote Sensing Data; Advances in Geographic Information Science; Bhatta, B., Ed.; Springer: Berlin/Heidelberg, Germany, 2010; pp. 17–36. [Google Scholar] [CrossRef]
  28. Tan, S.; Heerink, N.; Qu, F. Land fragmentation and its driving forces in China. Land Use Policy 2006, 23, 272–285. [Google Scholar] [CrossRef]
  29. Su, S.; Hu, Y.; Luo, F.; Mai, G.; Wang, Y. Farmland fragmentation due to anthropogenic activity in rapidly developing region. Agric. Syst. 2014, 131, 87–93. [Google Scholar] [CrossRef]
  30. McDonnell, M.; Pickett, S.T.A.; Groffman, P.; Bohlen, P.; Pouyat, R.; Zipperer, W.; Parmelee, R.; Carreiro, M.; Medley, K. Eco-system Processes Along an Urban-to-Rural Gradient. Urban Ecosyst. 2008, 1, 299–313. [Google Scholar] [CrossRef]
  31. Burton, M.L.; Samuelson, L.J.; Pan, S. Riparian woody plant diversity and forest structure along an urban-rural gradient. Urban Ecosyst. 2005, 8, 93–106. [Google Scholar] [CrossRef]
  32. Carreiro, M.M.; Tripler, C.E. Forest Remnants Along Urban-Rural Gradients: Examining Their Potential for Global Change Research. Ecosystems 2005, 8, 568–582. [Google Scholar] [CrossRef]
  33. Li, Y.; Li, Y.; Westlund, H.; Liu, Y. Urban–rural transformation in relation to cultivated land conversion in China: Implications for optimizing land use and balanced regional development. Land Use Policy 2015, 47, 218–224. [Google Scholar] [CrossRef]
  34. Niemelä, J.; Kotze, D.J.; Venn, S.; Penev, L.; Stoyanov, I.; Spence, J.; Hartley, D.; de Oca, E.M. Carabid beetle assemblages (Col-eoptera, Carabidae) across urban-rural gradients: An international comparison. Landsc. Ecol. 2002, 17, 387–401. [Google Scholar] [CrossRef]
  35. Clergeau, P.; Savard, J.-P.L.; Mennechez, G.; Falardeau, G. Bird Abundance and Diversity along an Urban-Rural Gradient: A Comparative Study between Two Cities on Different Continents. Condor 1998, 100, 413–425. [Google Scholar] [CrossRef]
  36. Garaffa, P.I.; Filloy, J.; Bellocq, M.I. Bird community responses along urban–rural gradients: Does the size of the urbanized area matter? Landsc. Urban Plan. 2009, 90, 33–41. [Google Scholar] [CrossRef]
  37. Yang, Y.; Guangrong, S.; Chen, Z.; Hao, S.; Zhouyiling, Z.; Shan, Y. Quantitative analysis and prediction of urban heat island intensity on urban-rural gradient: A case study of Shanghai. Sci. Total Environ. 2022, 829, 154264. [Google Scholar] [CrossRef]
  38. Wear, D.N.; Turner, M.G.; Naiman, R.J. Land Cover Along an Urban–Rural Gradient: Implications for Water Quality. Ecol. Appl. 1998, 8, 619–630. [Google Scholar] [CrossRef]
  39. Kaminski, A.; Bauer, D.M.; Bell, K.P.; Loftin, C.S.; Nelson, E.J. Using landscape metrics to characterize towns along an ur-ban-rural gradient. Landsc. Ecol. 2021, 36, 2937–2956. [Google Scholar] [CrossRef]
  40. Tlapáková, L.; Stejskalová, D.; Karásek, P.; Podhrázská, J. Landscape Metrics as a Tool for Evaluation Landscape Struc-ture—Case Study Hustopeče. Eur. Countrys. 2013, 5, 52–70. [Google Scholar] [CrossRef]
  41. Dale, V.H.; Offerman, H.; Frohn, R.; Gardner, R.H. Landscape Characterization and Biodiversity Research (No. CONF-9408220-1); Oak Ridge National Lab. (ORNL): Oak Ridge, TN, USA, 1995. [Google Scholar]
  42. Zhang, Y.; Su, T.; Ma, Y.; Wang, Y.; Wang, W.; Zha, N.; Shao, M. Forest ecosystem service functions and their associations with landscape patterns in Renqiu City. PLoS ONE 2022, 17, e0265015. [Google Scholar] [CrossRef]
  43. Tischendorf, L. Can landscape indices predict ecological processes consistently? Landsc. Ecol. 2001, 16, 235–254. [Google Scholar] [CrossRef]
  44. Uuemaa, E.; Antrop, M.; Roosaare, J.; Marja, R.; Mander, Ü. Landscape metrics and indices: An overview of their use in land-scape research. Living Rev. Landsc. Res. 2009, 3, 1–28. [Google Scholar]
  45. Cardille, J.A.; Turner, M.G. Understanding Landscape Metrics. In Learning Landscape Ecology: A Practical Guide to Concepts and Techniques; Gergel, S.E., Turner, M.G., Eds.; Springer: New York, NY, USA, 2017; pp. 45–63. [Google Scholar] [CrossRef]
  46. Li, H.; Reynolds, J.F. A Simulation Experiment to Quantify Spatial Heterogeneity in Categorical Maps. Ecology 1994, 75, 2446–2455. [Google Scholar] [CrossRef]
  47. Liu, J.; Jin, X.; Xu, W.; Zhou, Y. Evolution of cultivated land fragmentation and its driving mechanism in rural development: A case study of Jiangsu Province. J. Rural Stud. 2022, 91, 58–72. [Google Scholar] [CrossRef]
  48. Forman, R.T.T. Land Mosaics: The Ecology of Landscapes and Regions; Cambridge University Press: Cambridge, UK, 1995. [Google Scholar]
  49. Lin, Y.; Chang, T.; Wu, C.; Chiang, T.; Lin, S. Assessing Impacts of Typhoons and the Chi-Chi Earthquake on Chenyulan Wa-tershed Landscape Pattern in Central Taiwan Using Landscape Metrics. Environ. Manag. 2006, 38, 108–125. [Google Scholar] [CrossRef]
  50. Xiao, L.; Yu, M.; Wang, H. Hill, pond, field, forest, and residence: Analysis of rural settlement landscapes in the hilly areas of Sichuan and Chongqing. Landsc. Archit. 2023, 40, 71–77. [Google Scholar]
  51. Yang, X.; Liu, Z. Quantifying landscape pattern and its change in an estuarine watershed using satellite imagery and landscape metrics. Int. J. Remote Sens. 2005, 26, 5297–5323. [Google Scholar] [CrossRef]
  52. Huang, G. Ecological considerations on the spatial structure of mountainous cities. Urban Plan. 2005, 01, 57–63. [Google Scholar]
  53. Duany, A. Introduction to the Special Issue: The Transect. J. Urban Des. 2002, 7, 251–260. [Google Scholar] [CrossRef]
  54. Ren, H.; Zhao, Y.; Ge, Y. Spatial correlation between farmland fragmentation and landform types in karst mountainous areas: A case study of Xiuwen County, Guizhou Province. J. Guizhou Norm. Univ. Nat. Sci. Ed. 2020, 38, 1–9. [Google Scholar] [CrossRef]
  55. Xu, X. Ecological strategy research on green city design based on bioclimatic conditions. Ph.D. Thesis, Southeast University, Nanjing, China, 2007. [Google Scholar]
  56. Deng, O.; Zhou, X.; Huang, P.; Deng, L. Study on the correlation between soil nutrient spatial differentiation and topographic factors in the purple soil area of central Sichuan. Resour. Sci. 2013, 35, 2434–2443. [Google Scholar]
  57. Spearman, C. The proof and measurement of association between two things. Int. J. Epidemiol. 2010, 39, 1137–1150. [Google Scholar] [CrossRef]
  58. Wadduwage, S.; Millington, A.; Crossman, N.D.; Sandhu, H. Agricultural Land Fragmentation at Urban Fringes: An Applica-tion of Urban-To-Rural Gradient Analysis in Adelaide. Land 2017, 6, 28. [Google Scholar] [CrossRef]
  59. Penghui, J.; Dengshuai, C.; Manchun, L. Farmland landscape fragmentation evolution and its driving mechanism from rural to urban: A case study of Changzhou City. J. Rural Stud. 2021, 82, 1–18. [Google Scholar] [CrossRef]
  60. Jiang, P.; Li, M.; Lv, J. The causes of farmland landscape structural changes in different geographical environments. Sci. Total Environ. 2019, 685, 667–680. [Google Scholar] [CrossRef]
  61. Gao, X.; Song, Z.; Li, C.; Charlie, L.; Liang, S.; Tang, H. Spatial differentiation characteristics of multifunctional value of culti-vated land under urban-rural gradients. Trans. Chin. Soc. Agric. Eng. 2021, 37, 251–259. [Google Scholar]
  62. Yang, J.; Zhao, L.; Xu, F.; Yue, Y.; Du, Z.; Zhu, D. High-standard basic farmland construction zoning based on the connectivity of cultivated land. Trans. Chin. Soc. Agric. Mach. 2017, 48, 142–148. [Google Scholar]
  63. Lambin, E.F.; Gibbs, H.K.; Ferreira, L.; Grau, R.; Mayaux, P.; Meyfroidt, P.; Morton, D.C.; Rudel, T.K.; Gasparri, I.; Munger, J. Estimating the world’s potentially available cropland using a bottom-up approach. Glob. Environ. Chang. 2013, 23, 892–901. [Google Scholar] [CrossRef]
  64. Tang, X. “Agriculture-Nature Park” planning. Ph.D. Thesis, Chongqing University, Chongqing, China, 2018. [Google Scholar]
  65. Xing, Z.; Tang, X.; Zhou, Q.; Gu, Y.; Chen, Z. Planning for a green infrastructure network in urban fringe areas: Ensuring public welfare output. Urban Plan. 2020, 44, 57–69. [Google Scholar]
  66. Moreau, M. Transect urbanism: Readings in human ecology. Urban Res. Pract. 2021, 14, 483–484. [Google Scholar] [CrossRef]
  67. Kashian, D.M.; Tinker, D.B.; Turner, M.G.; Scarpace, F.L. Spatial heterogeneity of lodgepole pine sapling densities following the 1988 fires in Yellowstone National Park, Wyoming, USA. Can. J. For. Res. 2004, 34, 2263–2276. [Google Scholar] [CrossRef]
  68. Southworth, J.; Nagendra, H.; Tucker, C. Fragmentation of a Landscape: Incorporating landscape metrics into satellite analyses of land-cover change. Landsc. Res. 2002, 27, 253–269. [Google Scholar] [CrossRef]
Figure 1. Diagram of the current status of SAPs in the urban built-up area and urban fringe within the urban planning zone.
Figure 1. Diagram of the current status of SAPs in the urban built-up area and urban fringe within the urban planning zone.
Sustainability 16 06322 g001
Figure 2. The location of the study area. The upper right figure illustrates the geographic location of Chongqing Municipality within China, while the lower right figure indicates the position of Yongchuan District within Chongqing Municipality. The left figure depicts the location of the research area within Yongchuan District, including the extent of urban built-up areas within the research area.
Figure 2. The location of the study area. The upper right figure illustrates the geographic location of Chongqing Municipality within China, while the lower right figure indicates the position of Yongchuan District within Chongqing Municipality. The left figure depicts the location of the research area within Yongchuan District, including the extent of urban built-up areas within the research area.
Sustainability 16 06322 g002
Figure 3. Identification of small-scale agricultural and forestry land. Using ArcGIS 10.2, small-scale agroforestry land patches of under 3.3 hectares were selected, which are intricately intertwined with agroforestry protection zones, spanning across both urban construction areas and the peripheries of cities.
Figure 3. Identification of small-scale agricultural and forestry land. Using ArcGIS 10.2, small-scale agroforestry land patches of under 3.3 hectares were selected, which are intricately intertwined with agroforestry protection zones, spanning across both urban construction areas and the peripheries of cities.
Sustainability 16 06322 g003
Figure 4. Sample partitioning. Spatial transects of urban–rural areas were delineated from the city center to the boundaries of natural conservation areas, partitioned evenly in 2.8 km × 4 km intervals along the cardinal directions of southeast, southwest, northeast, and northwest.
Figure 4. Sample partitioning. Spatial transects of urban–rural areas were delineated from the city center to the boundaries of natural conservation areas, partitioned evenly in 2.8 km × 4 km intervals along the cardinal directions of southeast, southwest, northeast, and northwest.
Sustainability 16 06322 g004
Figure 5. Analysis of topographic and geomorphological. (a) represents the elevation within the spline space, (b) denotes the slope within the spline space, and (c) refers to the undulation within the spline space.
Figure 5. Analysis of topographic and geomorphological. (a) represents the elevation within the spline space, (b) denotes the slope within the spline space, and (c) refers to the undulation within the spline space.
Sustainability 16 06322 g005
Figure 6. Landscape indices transformation of SAPs. The figure illustrates the trend of spline space landscape indices transformation for SAPs in different urban–rural zones.
Figure 6. Landscape indices transformation of SAPs. The figure illustrates the trend of spline space landscape indices transformation for SAPs in different urban–rural zones.
Sustainability 16 06322 g006aSustainability 16 06322 g006b
Figure 7. Representation of woodland landscape indices transformation. The figure illustrates the trend of spline space landscape indices transformation for small-scale woodland in different urban–rural zones.
Figure 7. Representation of woodland landscape indices transformation. The figure illustrates the trend of spline space landscape indices transformation for small-scale woodland in different urban–rural zones.
Sustainability 16 06322 g007
Figure 8. Representation of cultivated land landscape indices transformation. The figure illustrates the trend of spline space landscape indices transformation for small-scale cultivated land in different urban–rural zones.
Figure 8. Representation of cultivated land landscape indices transformation. The figure illustrates the trend of spline space landscape indices transformation for small-scale cultivated land in different urban–rural zones.
Sustainability 16 06322 g008
Figure 9. Representation of garden land landscape indices transformation. The figure illustrates the trend of spline space landscape indices transformation for small-scale garden land in different urban–rural zones.
Figure 9. Representation of garden land landscape indices transformation. The figure illustrates the trend of spline space landscape indices transformation for small-scale garden land in different urban–rural zones.
Sustainability 16 06322 g009
Figure 10. Spearman correlation analysis of landscape indices for SAPs. Note: ** indicates p < 0.01; * indicates p < 0.05; RS < 0.7 is a low correlation, 0.7–0.9 is a moderate correlation, and >0.9 is a high correlation.
Figure 10. Spearman correlation analysis of landscape indices for SAPs. Note: ** indicates p < 0.01; * indicates p < 0.05; RS < 0.7 is a low correlation, 0.7–0.9 is a moderate correlation, and >0.9 is a high correlation.
Sustainability 16 06322 g010
Figure 11. Spearman’s correlation analysis of landscape indices of SAPs across the urban–rural gradient in each urban–rural zone.
Figure 11. Spearman’s correlation analysis of landscape indices of SAPs across the urban–rural gradient in each urban–rural zone.
Sustainability 16 06322 g011
Figure 12. Spearman correlation coefficients for landscape indices of small-scale woodland patches. Note: ** indicates p < 0.01; * indicates p < 0.05; RS < 0.7 is a low correlation, 0.7–0.9 is a moderate correlation, and >0.9 is a high correlation.
Figure 12. Spearman correlation coefficients for landscape indices of small-scale woodland patches. Note: ** indicates p < 0.01; * indicates p < 0.05; RS < 0.7 is a low correlation, 0.7–0.9 is a moderate correlation, and >0.9 is a high correlation.
Sustainability 16 06322 g012
Figure 13. Spearman’s correlation analysis of landscape indices of small-scale woodland patches across the urban–rural gradient in each sub-district.
Figure 13. Spearman’s correlation analysis of landscape indices of small-scale woodland patches across the urban–rural gradient in each sub-district.
Sustainability 16 06322 g013aSustainability 16 06322 g013b
Figure 14. Spearman correlation coefficients for landscape indices of small-scale cultivated land patches. Note: ** indicates p < 0.01; * indicates p < 0.05; RS < 0.7 is a low correlation, 0.7–0.9 is a moderate correlation, and >0.9 is a high correlation.
Figure 14. Spearman correlation coefficients for landscape indices of small-scale cultivated land patches. Note: ** indicates p < 0.01; * indicates p < 0.05; RS < 0.7 is a low correlation, 0.7–0.9 is a moderate correlation, and >0.9 is a high correlation.
Sustainability 16 06322 g014
Figure 15. Spearman’s correlation analysis of landscape indices of small-scale cultivated land patches across the urban–rural gradient in each urban–rural zone.
Figure 15. Spearman’s correlation analysis of landscape indices of small-scale cultivated land patches across the urban–rural gradient in each urban–rural zone.
Sustainability 16 06322 g015
Figure 16. Spearman correlation coefficients for landscape indices of small-scale garden land patches. Note: ** indicates p < 0.01; * indicates p < 0.05; RS < 0.7 is a low correlation, 0.7–0.9 is a moderate correlation, and >0.9 is a high correlation.
Figure 16. Spearman correlation coefficients for landscape indices of small-scale garden land patches. Note: ** indicates p < 0.01; * indicates p < 0.05; RS < 0.7 is a low correlation, 0.7–0.9 is a moderate correlation, and >0.9 is a high correlation.
Sustainability 16 06322 g016
Figure 17. Spearman’s correlation analysis of landscape indices of small-scale garden land patches across the urban–rural gradient in each urban–rural zone.
Figure 17. Spearman’s correlation analysis of landscape indices of small-scale garden land patches across the urban–rural gradient in each urban–rural zone.
Sustainability 16 06322 g017
Table 1. Yongchuan agricultural and forestry space statistics of various types of land area.
Table 1. Yongchuan agricultural and forestry space statistics of various types of land area.
TypeWoodlandCultivated LandGarden Land
Area (hm2) Percentage (%) Area (hm2) Percentage (%) Area (hm2) Percentage (%)
Over 3.3 ha 815.826.451865.426.56124.548.94
Under 3.3 ha2268.1373.555157.8873.441268.2991.06
Table 2. Classification of land uses for territorial space survey, planning, and use control.
Table 2. Classification of land uses for territorial space survey, planning, and use control.
Cultivated landPaddy fields
Irrigated land
Drylands
Garden landOrchards
Tea garden land
Oak park
Other garden land plots
Woodland Arbor woodland
Bamboo woodland
Shrub woodland
Other woodland
Table 3. Description of various types of landscape indices.
Table 3. Description of various types of landscape indices.
Landscape IndicesFormulaExplainExplain
Patch area (PA) P A = j = 1 n a i j 1 10000 a i j is the area of patch ijPA can reflect the change in the agroforestry scale on the urban–rural gradient.
Plaque density (PD) P D = n i A 10000 / 100 n i is the number of patches contained by patch type I, and A is the area of the whole landscape, including the background within the landscape. Units are 1/100 hm2.The patch density, defined as the number of patches per unit area, serves as an indicator of the fragmentation level in agricultural and forestry patches.
Edge density (ED) Ed = E/AE is the total perimeter of the plaque, and A is the total area of the patch.The higher the boundary density, the more the land-use type is divided, and the more dispersed the layout.
Largest plaque index (LPI) L P I = m a x ( a i j ) A a i j represents the area of patch ij; and A denotes the total landscape area, encompassing both the landscape interior and background.Reflecting dominant land-use patch types and landscape dominance in different partitions.
Fractal dimension (FD) F D = 2 l n ( 0.25 p i j ) l n ( a i j ) p i j is the perimeter of patch ij; a i j is the area of patch ij.To quantify the complexity of patch shapes, values closer to 1 indicate simpler shapes, while values closer to 2 signify more intricate shapes, indicative of heightened artificial influence.
Aggregation index (AI) G i = g i i k = 1 m g i k m i n e i
A I = G i p i p i   ( G i < p i   a n d   p i < 0.5 ) G i p i 1 p i
G i is the similar adjacency ratio, p i is the area proportion of the patch type in the landscape.Indicates the level of connectivity between patches; higher values denote increased connectivity.
Table 4. Statistics of gradient transformation of surface terrain data.
Table 4. Statistics of gradient transformation of surface terrain data.
SamplesAverage ElevationAverage SlopeAverage Rise and Fall
1-1276.415.2127.89
1-2279.586.1735.56
1-3283.265.1625.11
1-4293.688.1343.12
1-5345.3412.3284.50
2-1281.755.4728.65
2-2351.1712.2788.75
2-3627.0820.11172.52
2-4412.8415.65127.00
3-1304.287.7644.24
3-2336.145.3434.18
3-3285.564.9530.21
4-1286.976.2233.77
4-2281.835.4229.66
4-3266.734.8627.84
4-4243.975.0727.05
Table 5. The gradient transformation of artificial elements.
Table 5. The gradient transformation of artificial elements.
SamplesRoad Sites (ha) Park Site (ha) Building Surface (ha) The Surface of Hang Tong Lake (ha)
1-1148.7663.02750.424.10
1-2148.7172.86662.544.55
1-367.65.85307.0746.91
1-442.680.62109.4452.50
1-525.39067.6225.28
2-1226.6884.76651.9616.62
2-254.0670.1183.2732.86
2-319.090.0727.427.48
2-423.35059.229.82
3-1129.864.78408.9835.41
3-233.230132.55111.58
3-300167.0392.53
4-1110.4459.44566.6934.74
4-2109.79110.6402.1895.40
4-3120.114.33536.5348.72
4-456.35095.1997.29
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cheng, C.; Xing, Z.; Ye, L.; Yang, J.; Xie, Z. Characteristics and Influencing Factors of Landscape Pattern Gradient Transformation of Small-Scale Agroforestry Patches in Mountain Cities. Sustainability 2024, 16, 6322. https://doi.org/10.3390/su16156322

AMA Style

Cheng C, Xing Z, Ye L, Yang J, Xie Z. Characteristics and Influencing Factors of Landscape Pattern Gradient Transformation of Small-Scale Agroforestry Patches in Mountain Cities. Sustainability. 2024; 16(15):6322. https://doi.org/10.3390/su16156322

Chicago/Turabian Style

Cheng, Canhui, Zhong Xing, Lin Ye, Junyue Yang, and Zhuoming Xie. 2024. "Characteristics and Influencing Factors of Landscape Pattern Gradient Transformation of Small-Scale Agroforestry Patches in Mountain Cities" Sustainability 16, no. 15: 6322. https://doi.org/10.3390/su16156322

APA Style

Cheng, C., Xing, Z., Ye, L., Yang, J., & Xie, Z. (2024). Characteristics and Influencing Factors of Landscape Pattern Gradient Transformation of Small-Scale Agroforestry Patches in Mountain Cities. Sustainability, 16(15), 6322. https://doi.org/10.3390/su16156322

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop