*Article* **Relationships between Landscape Patterns and Hydrological Processes in the Subtropical Monsoon Climate Zone of Southeastern China**

**Chong Wei 1,2, Xiaohua Dong 1,2,\*,†, Yaoming Ma 3,4,5,6,7,8,†, Menghui Leng 1,2, Wenyi Zhao 1,2, Chengyan Zhang 1,2, Dan Yu 1,2 and Bob Su <sup>9</sup>**

	- Islamabad 45320, Pakistan

**Abstract:** With rapid economic development, extensive human activity has changed landscape patterns (LPs) dramatically, which has further influenced hydrological processes. However, the effects of LPs changes on hydrological processes, especially for the streamflow–sediment relationship in the subtropical monsoon climate zone, have not been reported. In our study, 10 watersheds with different sizes in the subtropical monsoon climate zone of southeastern China were chosen as the study area, and the effect of the 14 most commonly used landscape metrics (LMs) on 4 typical hydrological indices (water yields (WY), the runoff coefficient (RC), the soil erosion modulus (SEM), and the suspended sediment concentration (SSC)) were analyzed based on land use maps and hydrological data from 1990 to 2019. The results reveal that the LP characteristics within the study area have changed significantly. The number of patches and landscape shape indices were significantly positively correlated with watershed size (*p* < 0.01). For most watersheds, the largest patch index was negatively correlated with WY, RC, and SEM, and the perimeter area fractal dimension was positively correlated with WY, RC, SEM, and SSC. The effects of several LMs on the hydrological indices had scale effects. WY/RC and the interspersion and juxtaposition index were negatively correlated in most larger watersheds but were positively correlated in most smaller watersheds. Similar results were found for Shannon's diversity/evenness index and SEM. In general, an increase in a small patch of landscape and in landscape diversity would increase WY, the fragmentation of LPs would result in more soil erosion, and LPs would affect the relationship between streamflow and sediment yield. As a result, a proper decrease in landscape fragmentation and physical connectivity in the subtropical monsoon climate zone of southeastern China would benefit soil erosion prevention. These results enhance the knowledge about the relationship between LPs and hydrological processes in the subtropical monsoon climate zone of southeastern China and benefit local water and soil conservation efforts.

**Keywords:** landscape pattern; runoff coefficient; soil erosion modulus; suspended sediment concentration; subtropical monsoon climate zone; southeastern China

**Citation:** Wei, C.; Dong, X.; Ma, Y.; Leng, M.; Zhao, W.; Zhang, C.; Yu, D.; Su, B. Relationships between Landscape Patterns and Hydrological Processes in the Subtropical Monsoon Climate Zone of Southeastern China. *Remote Sens.* **2023**, *15*, 2290. https:// doi.org/10.3390/rs15092290

Academic Editor: Konstantinos X. Soulis

Received: 6 March 2023 Revised: 18 April 2023 Accepted: 25 April 2023 Published: 26 April 2023

**Copyright:** © 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

### **1. Introduction**

With the rapid development of human society, land use has changed dramatically [1], which has altered ecosystem structures, functions, and services [2], further influencing eco-hydrological processes [3,4]. The characteristics of land use changes mainly contain land use quantity and landscape patterns (LPs). Scholars have investigated the differences in hydrological effects for various land uses [5–8] to better determine a proper quantity for each land use type for these study regions. Afterward, some scholars further explored the evolution of land use spatial distribution based on various land use prediction models such as the CLUE-S [9], CA [10], and FLUS [11] models. In addition, the spatial distribution and configuration of various land use types could also affect hydrological processes, even with the same quantity. In recent decades, more and more scholars have focused on the effect of LPs on hydrological processes such as runoff and soil erosion [12–14].

LPs are the spatial configuration characteristics of various landscape units with different sizes and shapes [15]. Various landscape metrics (LMs) have been developed to describe the characteristics of LPs [16–18]. Over recent decades, scholars have indicated that LPs have significant correlations with runoff [12], soil erosion [13,14], water quality [19,20], and organic matter decomposition [21,22]. Sadeghi et al. [23] investigated the relationships between LMs and hydrographic components within the Galazchai watershed in Iran. The results show that the number of disjunct core areas (NDCA) was positively correlated with flood volume, and patch density (PD) was positively correlated with peak discharge. Zhao and Huang [24] examined the effects of LPs on runoff within a small watershed in southern China. The results indicate that PD and landscape shape index (LSI) were negatively correlated with runoff, whereas Shannon's diversity index (SHDI) and the landscape division index (DIVISION) were positively correlated with runoff. In addition, PD, SHDI, and largest patch index (LPI) had the largest impacts on variations in runoff. Zhang et al. [25] attempted to reveal the effect of changes in LPs on soil erosion in a small watershed. The results suggest that soil erosion was positively correlated with the patch cohesion index (COHESION) but negatively correlated with LPI, modified Simpson's evenness index (MSIEI), and aggregation index (AI). Another study conducted on the island of Crete in Greece [26] found that there was a significant correlation between average soil erosion and PD, edge density (ED), LPI, and percentage of landscape (PLAND). The topic of the effects of LPs variations on various hydrological processes has been a research hotspot over the last decade.

The most popular methods used to investigate the relationships between LPs and hydrological processes are various correlation analysis (the correlation coefficient method [27], stepwise regression analysis [28], multiple linear regression [25], etc.) methods. In addition, hydrological series have mainly been simulated using multiple hydrological or soil erosion models such as the SWAT [29], InVEST [30], IUH [23], RUSLE [26], WaTEM/SEDEM [31] models with various land use maps as inputs, which are further used for analysis with LMs calculated from the corresponding land use maps. The number of land use maps ranges from 3 to 8 in the previously mentioned works, which is insufficient and possibly brought uncertainty into the results. In addition, most of them used simulated hydrological series from multiple eco-hydrological models that ignored the effect of different LPs on hydrological processes. For instance, the SWAT model separates the watershed into several sub-watersheds and further divides them into more hydrological response units (HRUs) based on different combinations of land use, soil, and slope belt [32]. The hydrological processes were simulated at an HRU scale and then aggregated into a sub-watershed scale. The SWAT model is used to calculate the runoff and sediment generation based on specific combinations of land use, soil, and slope belt, but it does not consider the effects of LPs on the streamflow convergence process and the sediment transportation process. As reported by previous studies, the impacts of LPs on hydrological processes are mainly concentrated in the streamflow converge process and sediment transportation process [33]. Hence, simulated hydrological processes based on the SWAT model ignore the major effects of LPs on these hydrological processes. The RUSLE model simulates soil erosion using an experience

equation containing six factors [26], which also does not take LPs into consideration. In addition, the simulation results also contain errors compared to the observed data to some extent.

The existing relevant studies mainly focus on specific watersheds all over the world, and the areas of these watersheds range from 16.6 to 320,000 km2 [17,23,24,26]. As scale is a fundamental concept, analyses at various scales may have resulted in different results or even reverse results [34]. In addition, differences in soil type, topography, and meteorology within various regions could also have been a possible driver for the different results. For example, SHDI was positively correlated with runoff in some regions [17,35] and negatively correlated with runoff in other regions [36,37]; PD was positively correlated with the sediment yields load (SYL) in some regions [34,38] and negatively correlated with SYL in other regions [35,36]. Similar reverse results were also observed for COHESION, LPI, SHEI, AI, CONTAG, LSI, etc. in various regions [25,31,33,39]. These phenomena may suggest that it is insufficient to reveal the relationships between LPs and hydrological processes using relevant data from only one watershed. To investigate the impacts of LPs on water quality across Taiwan, Chiang et al. [40] revealed relationships between 12 water quality parameters and 12 LMs based on water quality records for 10 watersheds. The results show that temperature, PH, NO3-N, and TN are significantly correlated with most LMs in these watersheds. A similar relevant study has not been reported to have been conducted for runoff/SYL and LMs. In addition, the relationship between LMs and suspended sediment concentration (SSC) has not been reported until now. As SSC is an important hydrological index that could reflect the relationship between streamflow and suspended sediment in a watershed, understanding the driving factors behind SSC would benefit not only soil erosion prevention, but also the prevention and control of river reservoir siltation. As the existing studies report, LPs influence runoff and SYL, mainly during the transportation process. It is necessary to conduct a relevant study to investigate whether LPs changes affect SSC to better understand the driving factors behind SSC.

In previous studies, scholars investigated the relationships between LPs and hydrological processes in specific watersheds, but relevant studies focused on climate zones are still rare. Soil erosion brought by rainfall and runoff is strongly controlled by climate [41], and it is important to investigate the effects of LPs on hydrological processes in various climate zones. The subtropical monsoon climate zone is characterized by high annual rainfall and concentrated summer rainfall, which would result in severe soil erosion [42]. The subtropical monsoon climate zone covers ~29% of China's total area, is mainly distributed in southeastern China, and has experienced rapid economic development and land use changes. In addition, southeastern China contains a variety of landforms and shows a general conversion from mountains and hills to plains from south to north, which are more likely to suffer soil erosion. As pronounced by the Resource and Environment Science and Data Center, Chinese Academy of Sciences (RESDC), ~56.39% of Chinese land suffers from water erosion, with this percentage accounting for almost all of southeastern China [43]. However, how LPs changes affect hydrological processes within the subtropical monsoon climate zone of southeastern China is still unclear. Under the background of global climate change and land use change resulting from extensive human activity, it is significant to reveal the relationships between LPs and representative hydrological processes for better land use management for water and soil conservation. To settle this issue, 10 watersheds with area ranging from 1700 to 80,900 km<sup>2</sup> were chosen as the study area, and water yields (WY), runoff coefficient (RC), soil erosion modulus (SEM), and SSC were chosen as the representative hydrological indices to determine the correlations between the chosen hydrological indices and the most studied LMs from 1990 to 2019 in the subtropical monsoon climate zone of southeastern China. The results enhance the knowledge about the effects of LPs on hydrological processes in the subtropical monsoon climate zone of southeastern China for better land use management with an objective of soil and water conservation.

### **2. Study Area**

A total of 10 watersheds (Zhuji (ZJ), Dufengkeng (DFK), Hushan (HS), Lijiadu (LJD), Lanxi (LX), Boluo (BL), Chaoan (CA), Shijiao (SJ), Zhuqi (ZQ), and Waizhou (WZ)) of different sizes in southeastern China, which were dominated by the subtropical monsoon climate, were chosen as the study area (Figure 1). The study area covered a region between 112.1◦~120.7◦E, and 22.5◦~30.0◦N, which is mainly dominated by hills and mountainous terrain, with its highest elevation being 2191 masl and its lowest elevation being −136 masl. The main climatological characteristics in the subtropical monsoon climate zone are temperature and precipitation, which are higher in the summer and lower in the winter, and annual precipitation, which is high but unevenly distributed throughout the year. The annual average temperature and precipitation in southeastern China range from 15 to 22 ◦C and from 1500 to 2200 mm, respectively. The region is rich with river networks and has a developed water system, which is mainly replenished by precipitation. The main geographical and hydrological characteristics of the chosen watersheds are shown in Table 1. The dominant land use of the study area is forest and agricultural land, and LPs have changed rapidly in the last 30 years due to intensive human activity. The main LP changes were an expansion of urban land and water areas and a reduction in forest and agricultural land from 1990 to 2019.

**Figure 1.** The distribution of the study regions. "ZJ", "DFK", "HS", "LJD", "LX", "BL", "CA", "SJ", "ZQ", and "QZ" represent the "Zhuji", "Dufengkeng", "Hushan", "Lijiadu", "Lanxi", "Boluo", "Chaoan", "Shijiao", "Zhuqi", and "Waizhou" watersheds, respectively.


**Table 1.** Hydrological and geographical characteristics of the chosen watersheds from 1990 to 2021.

Note: DA, AE, and PRE represent drainage area, average elevation, and precipitation, respectively.

### **3. Data and Methods**

### *3.1. Data Description*

The main data used in this study consist of land use maps, precipitation records, and annual runoff and sediment yield records. The land use maps from 1990 to 2019 were obtained from Wuhan University [44] with a spatial resolution of 30 m. The land use maps were generated using the random forest classifier and the visual interpretation method, and they were validated with three sources of test samples. They divided the land use type into 9 categories: cropland, forest, shrub, grassland, water, snow and ice, barren, imperious, and wetland. They can be obtained at https://zenodo.org/record/4417810 (accessed on 28 December 2022), and they have been employed in many studies [19,45,46]. The precipitation records were obtained from China Scientific Data with a temporal resolution of monthly (1960–2020) and a spatial resolution of 1 km [47]. The precipitation records were interpreted with ANUSPLIN 4.4 based on ~2400 precipitation monitoring stations within China and were validated by comparing them with the precipitation records from precipitation gauges and the Chinese hydrological yearbook. The precipitation records that were used can be downloaded at http://doi.org/10.11922/sciencedb.01607 (accessed on 25 December 2022). The measured annual runoff and sediment yield records from 1990 to 2019 of the DFK, LJD, LX, BL, CA, SJ, ZQ, and WZ watersheds and those from 2002 to 2019 of the ZJ and HS watersheds were obtained from the Bulletin of River Sediment in China and were downloaded at http://www.mwr.gov.cn/sj/tjgb/zghlnsgb/ (accessed on 30 December 2022). The bulletins were compiled jointly by the Department of Hydrology of the Ministry of Water Resources, the Monitoring and Forecasting Center for Hydrology and Water Resources of the Ministry of Water Resources, the International Sediment Research and Training Center, and the hydrology bureaus of the river basin institutions, which are very credible. The annual runoff and sediment yield records were used to represent the historical hydrology conditions within the chosen watersheds.

### *3.2. LMs Selection*

LMs were used to represent the spatial distribution characteristics of the various land use patches, and they have been widely used to investigate the relationships between LPs and eco-hydrological processes [21,25,39,48]. The LMs were divided into three categories based on different levels: patch-level, class-level, and landscape-level metrics [39]. The LMs on a landscape level integrated all patch types or classes over the entire study region, and they were separated into four categories reflecting different types of characteristics of landscapes: edge area, shape, aggregation, and diversity metrics [49]. The main purpose of this study was to reveal the relationships between landscape patterns and representative hydrological indices in the subtropical monsoon climate zone. Hence, 14 widely used LMs, which were found to have different relationships with runoff/sediment in various regions, were chosen for this study [17,25,27,31,33,35,36,38,39,50,51]. The selected LMs were LPI, ED, perimeter area fractal dimension (PAFRAC), number of patches (NP), DIVI- SION, AI, interspersion and juxtaposition index (IJI), contiguity index (CONTAG), PD, LSI, COHESION, SHDI, Shannon's evenness index (SHEI), and MSIEI. Their definitions and relevant literature are shown in Table A1, and these LMs were calculated using Fragstats v4.2.1 software for all land use maps from 1990 to 2019 for all chosen watersheds. Before calculation, the land use types were reclassified into seven categories: agricultural land, forest land, shrubland, grassland, water area, urban land, and bare land.

### *3.3. Hydrological Indices*

To better evaluate the relationships between hydrological processes and LPs, five fundamental hydrological indices (WY (m3); RC; SYL (t); SEM (t/km2); SSC (mg/L)) relevant to runoff, sediment yields, and streamflow–sediment relationships were chosen for the study and obtained on an annual scale. These hydrological indices have been employed in previous studies to investigate the hydrological characteristics of various regions [29,52–55]. In this paper, WY, SYL, and SSC values were collected from the Bulletin of River Sediment in China. RC and SEM were calculated using Equations (1) and (2).

$$\text{RC} = \frac{\text{WY}}{1000 \cdot \text{DA} \cdot \text{PRE}} \tag{1}$$

$$\text{SEM} = \frac{\text{SYL}}{\text{DA}} \tag{2}$$

where WY represents annual water yield (m3); DA represents drainage area (km2); PRE represents annual precipitation (mm); SYL represents the annual sediment yield load (t).

### *3.4. Analysis Methods*

Linear regression [56] was used to analyze the change trends of various hydrological indices and of LMs from 1990 to 2019 for the chosen watersheds. This type of analysis has been widely employed in previous relevant works and has shown reasonable results [43,57,58]. In addition, the F-test [59] was adopted to test the significance of the trend for each series at confidence levels of 0.05 and 0.01. The equation form of the linear regression is as follows:

$$\mathbf{y} = \mathbf{a}\mathbf{x} + \mathbf{b} \tag{3}$$

where y represents the dependent variable (WY, RC, SYL, SEM, and SSC); a represents the change trend per year; x represents year; b represents the constant.

To investigate the exchanges between different land use types from 1990 to 2019, the period was separated into two stages: 1990–2005 and 2005–2019. Land use conversion analysis was conducted for land use maps in 1990, 2005, and 2019 in the ArcGIS 10.2 platform. To better show the exchanges between different land use types, land with no land use type exchanges was ignored when conducting relevant analyses. With this method, the exchanges between each land use type from 1990 to 2019 were analyzed for the chosen watersheds.

The effects of LPs changes that propagate to hydrological processes have lag times that vary according to differences in watershed size, climate, landscape, topography, soil type, etc. [60–62]. In previous studies, lag times of 0–4 years were needed before significant responses of hydrological processes to landscape changes could be observed [63,64]. Hence, a slip correlation analysis [43] was conducted with a lag time of 0–4 years to investigate the relationships between LPs and hydrological indices. A description of the slip correlation analysis is as follows:


(3) The correlation coefficient (CCi) between X and Y was calculated, and the result was the lag correlation between X and Y with a lag time of i years. The highest value of CC*<sup>i</sup>* with a significance of 0.05 was regarded as the lag time between X and Y.

As SEM was calculated by dividing SYL by DA, the SEM series was linearly correlated with SYL and removed the effect of watershed size on soil erosion. In this paper, a slip correlation analysis was conducted between the 14 LMs and 4 hydrological indices (WY, RC, SEM, and SSC) in the chosen watersheds.

### **4. Results**

### *4.1. Variations in SY, RC, SYL, SEM, and SSC*

Table 1 shows the annual average values of various hydrological indices, and Figure 2 shows variations in WY, RC, SYL, SEM, and SSC. The annual average WY, SYL, RC, SEM, and SSC values for the chosen watersheds ranged from 12.17 × <sup>10</sup><sup>8</sup> to 713.52 × <sup>10</sup><sup>8</sup> <sup>m</sup>3, from 8.9 × <sup>10</sup><sup>4</sup> to 463.7 × <sup>10</sup><sup>4</sup> t, from 0.39 to 0.55, from 47.3 to 162.75 t/km2, and from 0.05 to 0.13 mg/L, respectively. The variations in WY were higher within ZQ and WZ, followed by LJD, LX, BL, CA, and SJ, and the variations in SYL were obviously higher within CA, SJ, ZQ, and WZ compared to the other watersheds. In addition, the variations in SEM and SSC were higher within DFK, HS, LX, and CA compared to the other watersheds, which indicates that soil erosion in these watersheds was much more serious compared to the other watersheds.

**Figure 2.** Variations in hydrological indices from 1990 to 2019. (**a**–**e**) represent WY, RC, SYL, SEM, and SSC, respectively.

The change trends of the hydrological indices are shown in Table 2, and their temporal variations are shown in Figure A1. WY showed a decreased trend in the DFK, BL, CA, SJ, and WZ watersheds while exhibiting an increased trend in the ZJ, HS, LJD, LX, and ZQ watersheds. RC increased in all watersheds except for CA and SJ. In particular, DFK, BL, and WZ suffered decreases in WY, whereas their RC values increased. As a whole, the change trends of WY and RC were not significant for all watersheds. Differing from WY and RC, SYL and SEM decreased in all watersheds except for the ZJ, HS, and LX watersheds. SYL and SEM values increased significantly (*p* < 0.01) in HS and decreased significantly (*p* < 0.01) in CA and WZ. Interestingly, LJD and ZQ increased in WY and RC but decreased in SYL and SEM. This phenomenon may have been caused by LP changes and human activity in these regions. As for SSC, it significantly increased in HS (*p* < 0.01) and

significantly (*p* < 0.05) decreased in LJD, ZQ, BL, CA, and WZ. In general, soil erosion was mitigated in most watersheds, whereas water resources decreased in the larger watersheds and increased in the smaller watersheds.

**Table 2.** Change trend statistics of hydrological indices from 1990 to 2019.


Note: The units of WY, SYL, SEM, and SSC are 108 m3, 104 t, t/km2, and mg/l, respectively, and RC is a dimensionless quantity. The change trends of these indices are shown as the average variations for each year. \* and \*\* represent significant correlations at confidence levels of 95% and 99%, respectively.

### *4.2. Temporal Variations in LPs in the Chosen Watersheds*

4.2.1. Land Use Changes from 1990 to 2019

As shown in Figure 3, the dominant land use type of the chosen watersheds was forest land, followed by agricultural land. The ratio of agricultural land was higher in ZJ, LJD, LX, BL, and WZ compared to the other watersheds, which were also more intensively distributed. In addition, from 1990 to 2019, urban land expanded drastically in all watersheds, especially in ZJ, LJD, LX, and WZ. The regions in which the land use types for different watersheds were altered occupied about 3.01–8.46% (average of 6.03%), 4.22–9.78% (average of 6.81%), and 5.11–11.47% (average of 8.83%) of the total area from 1990 to 2005, from 2005 to 2019, and from 1990 to 2019, respectively. Interestingly, the average area of the regions whose land use types had been altered from 1990 to 2019 was much lower than the sum of these regions from 1990 to 2005 and from 2005 to 2019. This may suggest that some unreasonable land use conversions from 1990 to 2005 had been corrected from 2005 to 2019. From 1990 to 2019, urban expansion was the only long-lasting land use conversion process.

Figure 4 shows the details of land use conversions among different land use types from 1990 to 2005 and from 2005 to 2019. We found that land use conversions mainly happened on agricultural land, forest land, and urban land. From 1990 to 2005, the main land use conversions were forest land–agricultural land, agricultural land–urban land, and agricultural land–forest land; part of the agricultural land had also been converted into water areas within ZJ, HS, LJD, LX, BL, and SJ. From 2005 to 2019, the main land use conversions were agricultural land–urban land, agricultural land–forest land, and forest land–agricultural land; some urban land was converted from forest land within DFK, HS, LX, CA, ZQ, and WZ. These results reveal that the urban expansion in southeastern China was mainly facilitated by occupying agricultural land, especially from 1990 to 2005. During recent decades, urban expansion also happened in some forest regions for most watersheds. This was reasonable, as urban development had extended to the surrounding rural areas first in most regions, which had large areas of agricultural land and forest land [65–67]. In addition, conversions between agricultural land and forest land were relatively greater. These results may have been caused by land use exchanges. For better soil erosion prevention, some sloped agricultural land was abandoned and converted into forest land, whereas some flat forest land was converted into agricultural land to preserve the amount of agricultural land and promote guaranteed crop production.

**Figure 3.** Land use maps of ZJ, DFK, HS, LJD, LX, BL, CA, SJ, ZQ, and WZ in 1990, 2005, and 2019.

4.2.2. LM Changes from 1990 to 2019

As shown in Figure 5, the variations in all LMs of ZJ were the largest among all watersheds. In addition, several LMs of LX, ZQ, and WZ also showed larger variations compared to the other watersheds. The average values of LPI and PAFRAC had similar trends with increases in watershed size. The analogical character was also shown for DIVISION, IJI, PD, SHDI, SHEI, and MSIEI. It is notable that the LPs within ZJ changed intensively in 1998, 2001, 2007, and 2017, with obvious variations in LPI, DIVISION, and COHESION. The LPs within LX changed intensively in 1993, 2007, and 2016 (Figure A2). Different from the other LMs, COHESION changed little in all watersheds except ZJ compared to the other LMs. This may suggest that COHESION was not as sensitive as other LMs to LPs changes.

**Figure 4.** Land use conversion statistics between 1990, 2005, and 2019 for the chosen watersheds. Land whose land use type did not change is not included in this figure. (**a**–**j**) represent the ZJ, DFK, HS, LJD, LX, BO, CA, SJ, ZQ, and WZ watersheds, respectively.

The change trend statistic results of the chosen LMs from 1990 to 2019 are shown in Table 3. For most watersheds, most LMs changed significantly (*p* < 0.05), which demonstrates that LPs characteristics changed significantly during this period. Among all watersheds, the LPs within ZQ was the most intensively changed, as all LMs of ZQ changed significantly (*p* < 0.01). In detail, COHESION (ED, PAFRAC, AI, IJI, CONTAG, LSI, and SHEI) and SHDI (LPI and DIVISION) significantly (*p* < 0.05) changed in nine (individually, eight and seven, respectively) watersheds; NP and PD (MSIEI) significantly (*p* < 0.01) changed in six (four) watersheds. Among these significant trends, ED, DIVISION, IJI, LSI, SHDI, SHEI, and MSIEI (AI, CONTAG, COHESION, LPI, and PAFRAC) all significantly (*p* < 0.05) increased (or decreased) in most watersheds. These uniformities in the change trends of the LMs may suggest that the development directions in land use planning were similar across the regions of southeastern China. In addition, COHESION significantly decreased in most watersheds, which did not show obvious variations. The change trends of SHDI and SHEI were similar in all watersheds, and the change trend of SHDI was more obvious than that of SHEI (Figure A2). In general, the LPs within the chosen watersheds significantly changed from a landscape perspective.

**Figure 5.** Variations of LMs from 1990 to 2019. (**a**–**n**) represent LPI, ED, PAFRAC, NP, DIVISION, AI, IJI, CONTAG, PD, LSI, COHESION, SHDI, SHEI, and MSIEI, respectively.


**Table 3.** Change trend statistics of LMs within the chosen watersheds from 1990 to 2019.

Note: The units of LPI, ED, DIVISION, AI, IJI, CONTAG, and PD are %, meters/hectare, proportion, %, %, %, and number/100 hectares, respectively, and PAFRAC, NP, LSI, COHESION, SHDI, SHEI, and MSIEI are dimensionless quantities. The change trends of these indicators are shown as the average variations each year. \* and \*\* represent significant correlations at confidence levels of 95% and 99%, respectively.

### *4.3. Relationships between Hydrological Indices and LPs*

The CC between the hydrological indices (WY, RC, SEM, and SSC) and LMs (LPI, ED, PAFRAC, NP, DIVISION, AI, IJI, CONTAG, PD, LSI, COHESION, SHDI, SHEI, and MSIEI) of the chosen watersheds are shown in Tables 4 and 5. There were discrepancies in the lag times of the correlations between the hydrological indices and LMs in different watersheds in the subtropical monsoon climate zone of southeastern China. In addition, the lag times between the hydrological indices and LMs did not show an obvious relationship with watershed size.

**Table 4.** CC values between the landscape metrics and the WY and RC.


Note: bold font indicates a correlation was significant at a confidence level of 95%, bold and underlined font indicates a correlation was significant at a confidence level of 99%, and a superscript number indicates the lag time of the correlation between the LMs and the WY or RC.


**Table 5.** CC values between LMs and the SEM and SSC.

Note: bold font indicates a correlation was significant at a confidence level of 95%, bold and underlined font indicates a correlation was significant at a confidence level of 99%, and a superscript number indicates the lag time of the correlation between the LMs and the WY or RC.

In most watersheds, WY and RC were negatively correlated with LPI, AI, and COHE-SION and positively correlated with ED, PAFRAC, NP, DIVISION, PD, LSI, SHDI, SHEI, and MSIEI; SEM was negatively correlated with LPI and IJI and positively correlated with PAFRAC, NP, PD, COHESION, and MSIEI; SSC was negatively correlated with IJI and positively correlated with PAFRAC, NP, PD, COHESION, and MSIEI. In addition, the correlation between WY/RC and IJI varied in watersheds with different sizes. In general, runoff decreases when a landscape becomes more fragmented or its shape gains complexity [51], and sediment transportation is disturbed when a landscape patch becomes larger or fragmented, further reducing sediment yields [35]. In addition, the results show that the effect of LPs on SSC and SEM were similar, which indicates that the LPs affect SSC by affecting SEM. It can be seen in Table 4 that most watersheds with larger areas had negative correlations between WY/RC and IJI, whereas most watersheds with smaller areas had positive correlations between WY/RC and IJI. A similar phenomenon was found between SHDI/SHEI and SEM. Most watersheds with larger areas had negative correlations between SHDI/SHEI and SEM, whereas most watersheds with smaller areas had positive correlations between SHDI/SHEI and SEM.

The effects of LMs on different hydrological indices varied, even in the same watersheds. For instance, for DKF, no significant (*p* < 0.05) correlations were found between WY and any LMs, but all LMs except for PAFRAC, IJI, and COHESION were significantly (*p* < 0.05) correlated with RC; for WZ, all LMs except for LPI and DIVISION did not show significant (*p* < 0.05) correlations with WY, whereas all LMs except for PAFRAC, NP, PD, and COHESION were significantly (*p* < 0.05) correlated with RC. In general, the effects of the LMs on SEM and SSC were more significant than their effects on WY and RC. Among the chosen watersheds, only four watersheds (LJD, BL, SJ, and ZQ) did not show significant (*p* < 0.05) correlations between SEM and the LMs. In particular, all LMs were significantly (*p* < 0.01) correlated with the SEM in HS; seven, ten, and eleven LMs were significantly (*p* < 0.05) correlated

with the SEMs in ZJ, LX, and CA, respectively. Based on the CC between SSC and all LMs, it seems that LPs changes could closely affect the relationship between streamflow and sediment yields. SSC was found to be significantly (*p* < 0.05) correlated with at least one LM in all of the chosen watersheds. In particular, all LMs were significantly (*p* < 0.01) correlated with the SSC in HS, and eight, eleven, eight, and twelve LMs were significantly (*p* < 0.05) correlated with the SSCs in ZJ, LX, BL, and CA, respectively.

### **5. Discussion**

### *5.1. Correlations between LMs and Watershed Size*

Different LMs represent different LPs characteristics, which may be different in different regions due to the varied landscape. In our study, the average LMs were correlated with watershed size (Figure 6). In detail, ED, DIVISION, PD, SHDI, SHEI, and MSIEI had similar correlations to watershed size. These LMs were negatively correlated with watershed size when it ranged from 20,000 to 60,000 km2 and positively correlated with watershed size when it was in other ranges. LPI, AI, and CONTAG had similar correlations with watershed size. These LMs showed reverse correlations with watershed size compared to ED, DIVISION, PD, SHDI, SHEI, and MSIEI. Generally, ED, DIVISION, PD, SHDI, SHEI, and MSIEI represent the fragmentation of the landscape, whereas LPI, AI, and CONTAG represent the aggregation of the landscape. This result indicates that the LMs that represent similar landscape characteristics had similar correlations with watershed size, whereas the LMs with opposite landscape characteristics had reverse correlations with watershed size. The thresholds of 20,000 and 60,000 km2 may suggest that fragmentation of the landscape would be decreased in southeastern China where the area ranges from 20,000 to 60,000 km2. PAFRAC and COHESION showed similar correlations with watersheds, with slow increasing trends where the watershed size ranged from 0 to 30,000 km2. Interestingly, NP and LSI were significantly (*p* < 0.01) positively correlated with watershed size. The increasing trend was smaller where the watershed size ranged from 25,000 to 55,000 km<sup>2</sup> compared to other ranges. NP represents the number of patches in the watershed, and LSI represents the ratio of edge and area of landscape. Generally, NP and LSI would increase with the increase in watershed size, but the slight increasing trend of NP and LSI where the watershed size ranges from 25,000 to 55,000 km<sup>2</sup> may suggest that the number of patches and the ratio of the edge to the area of a landscape would not be changed obviously where the area of watershed ranges from 25,000 to 55,000 km2. When the watershed area fell in this scope, the increased area mainly had less fragmentation. These correlations between LMs and watershed size may suggest that the effects of LMs on hydrological processes also have scale effects, as suggested by Xiao, Cao, Liu, and Lu [34]. More LMs should be taken into consideration within more watersheds with wider size ranges to further verify this interesting hypothesis.

### *5.2. Effects of Various LMs on Hydrological Processes*

The relationships between four hydrological indices (WY, RC, SEM, and SSC) and fourteen LMs (LPI, ED, PAFRAC, NP, DIVISION, AI, IJI, CONTAG, PD, LSI, COHESION, SHDI, SHEI, and MSIEI) within ten representative watersheds with various area sizes in the subtropical monsoon climate zone of southeastern China were evaluated using the slip correlation analysis method. In our study, for most watersheds, the land use maps and hydrological series for 30 consecutive years were used to investigate the relationships between landscapes and hydrological processes. Compared to previous studies, the amount of land use maps was much greater, and the hydrological series were measured rather than simulated with various models, which could greatly reduce the uncertainty brought by various models. In addition, the measured hydrological series under different landscapes could directly reflect the effect of landscape on hydrological processes. Taking the SWAT model as an example, as mentioned in the Introduction section, it simulates hydrological processes based on different combinations of land use, soil types, and slope belt. The SWAT model considers the land use amount but ignores the effects of land use spatial variations

on hydrological processes. For example, an HRU in a subwatershed would distribute in several places within the subwatershed, which makes it hard to represent the effect of land use spatial distributions on hydrological processes. Hence, it is more reliable to investigate the relationships between landscapes and hydrological processes using measured data compared to using simulated data from various models.

**Figure 6.** Statistics of the watershed area and average LMs for the chosen watersheds from 1990 to 2019. (**a**–**n**) represent LPI, ED, PAFRAC, NP, DIVISION, AI, IJI, CONTAG, PD, LSI, COHESION, SHDI, SHEI, and MSIEI, respectively.

Similar to previous studies, the relationships between the different hydrological indices and LMs in various watersheds had discrepancies [35]. In our study, some hydrological indices and LMs had the same correlations in most watersheds, and a few hydrological indices and some LMs even had the same correlation in all watersheds. For most watersheds (≥7), WY was negatively correlated with LPI, AI, CONTAG, and COHESION and positively correlated with ED, PAFRAC, DIVISION, LSI, SHDI, SHEI, and MSIEI; RC was negatively correlated with LPI, CONTAG, and COHESION and positively correlated with PAFRAC, DIVISION, SHDI, SHEI, and MSIEI; SEM was negatively correlated with LPI and IJI; and SEM and SSC were positively correlated with PAFRAC, NP, PD, COHESION, and MSIEI. Interestingly, several LMs had the same relationships with WY and RC in all watersheds. As for SEM and SSC, there were no consistent correlations with any LMs for any watersheds. This discrepancy may have been caused by the soil erosion process being easier to disturb with human activity than the runoff process [68,69]. Hence, the relationships between LMs and SEM/SSC have discrepancies in various watersheds with differences in human activity intensity and mode.

LPI and ED were both related to the edge-area characteristic of LPs, and they had reverse effects on the runoff process. As mentioned before, LPI was negatively correlated with WY and RC, whereas ED was reversed. This result suggests that runoff increases in the subtropical monsoon climate zone when a watershed is dominated by a small patch of landscape. Similar results have also been reported in other regions [17], and a reverse result was also obtained [38]. In addition, SEM was negatively correlated with LPI in most watersheds, whereas it had no consistent correlation with ED. The dominant land use was forest and agricultural land, which may suggest that the dispersal of forest and agricultural land would increase soil erosion. The negative correlation between IJI and SEM also proved this result. A similar result was obtained by Zhang, Fan, Li and Yi [25]. In our study, seven aggregation LMs were chosen to analyze their correlation with four hydrological indices. For most watersheds, NP, PD, LSI, and DIVISION were positively correlated with WY and RC, whereas NP, PD, and COHESION were positively correlated with SEM and SSC. These results suggest that a proper increase in the fragmentation and physical connectivity of LPs was beneficial for water and soil conservation. A similar result was reported within the Upper Du River Basin in the middle of the Yangtze River [51], which is also within the subtropical monsoon climate zone. In addition, reverse results have been found in northern China, which belongs to the semi-arid continental monsoon climate zone [36,39] and may be caused by differences in rainfall patterns in different climate zones. The average annual precipitation in the study regions ranges from 1686 to 2162 mm (Table 1), whereas the average annual precipitation within the semi-arid continental monsoon climate zone ranges from 300 to 500 mm [70]. Furthermore, Ouyang, Skidmore, Hao, and Wang [33] pointed out that properly improving the fragile landscape status could prevent soil erosion, which may suggest that there is a threshold for landscape fragility. The changes in the fragmentation in the upper or lower fragmentation thresholds may be different, which may be another reason for the different correlations between relevant LMs and soil erosion. WY and RC were positively correlated with SHDI, SHEI, and MSIEI in most watersheds, which suggests that an increase in landscape diversity would increase runoff generation in the subtropical monsoon climate zone. As for SEM and SSC, the effect of landscape diversity varied in different regions [36,39,71], which may have resulted from discrepancies in land use type conversions. The main land use conversions from 1990 to 2019 within these watersheds were forest land–agricultural land, agricultural land–urban land, and agricultural land–forest land (Figure 4). Generally, small agricultural patches increased in a large forest patch, small urban patches increased in a large agricultural patch, and small forest patches increased in a large agricultural patch, which increased landscape diversity but may have affected soil erosion differently. A further investigation should be conducted to verify this hypothesis.

The correlations between LMs and RC were more obvious than those between LMs and WY (Table 4). There are no glaciers in the subtropical monsoon climate zone, and the main source of water yields was precipitation recharge [72]. Precipitation variation plays a more important role in interannual WY changes compared to LPs variations [73]. RC removed the effect of precipitation variations and watershed size to some extent compared to WY, and it was mainly affected by LULC, soil property, and slope [74]. These results suggest that relevant studies should choose RC to represent runoff generation character rather than WY. In addition, the relationships between LMs and SEM/SSC were not consistent in different watersheds. This may have been caused by the erosion process, which was also influenced by precipitation patterns, soil properties, slope, and topography [73,75–77]. WY/RC and SEM/SSC were both positively correlated with PAFRAC, NP, PD, and MSIEI in most watersheds. This shows that the soil erosion in these regions is mainly caused by rainfall erosion and surface runoff transportation [78]. In addition, the extensive construction projects brought about by rapid urban expansion also could have led to a lot of soil erosion [79], further influencing these results.

### *5.3. Recommendations for LM Selection in Future Relevant Studies*

Our study selected 10 watersheds located in the subtropical monsoon climate zone in southeastern China. For the 10 chosen watersheds, the temporal variations and change trends of the 14 chosen LMs from 1990 to 2019 are shown in Figure A2 and Table 3. The results showed that the change trends of LPI and ED were opposite in all watersheds except for ZJ; the change trend of AI was opposite to LSI in all watersheds and similar to CONTAG and COHESION in all watersheds except for CA; and finally, SHDI and SHEI had similar change trends in all watersheds. In addition, the temporal variation in CONTAG was more obvious than those of AI and COHESION, but the change trend of COHESION was more significant compared to those of AI and CONTAG. SHDI was less disturbed by landscape changes compared to SHEI, as SHEI changed abruptly in some years, whereas similar changes did not appear for SHDI (Figure A2). In addition, the correlations with WY/RC for LPI and ED were opposite; the correlations with WY/RC for AI, CONTAG, and COHESION were similar, and COHESION had a positive correlation with SEM/SSC in most watersheds; the correlations with WY/RC for SHDI and SHEI were similar. Similar results were also reported in previous studies [17,35,37,80]. Hence, when analyzing the correlations between LPs and hydrological processes, it would be better to choose LPI, ED, PAFRAC, AI, COHESION, and SHDI as the representative LMs for fundamental LPs characteristics that include the edge area, shape, aggregation, and diversity aspects.

### *5.4. Implications, Limitations, and Prospects*

Understanding the relationships between LPs and hydrological processes has a benefit for water and environmental management [27]. The results suggest that a proper decrease in landscape fragmentation and connectivity in the subtropical monsoon climate zone of southeastern China would benefit soil erosion prevention. In addition, meaningless land use conversions (forest land–agricultural land and agricultural land–forest land) should be decreased during local economic development, which would benefit decreasing landscape fragmentation and connectivity. In addition, the effects of LPs on SSC were also revealed in our study. Soil erosion and sediment transportation have a similar response to LPs changes, which indicates that it is possible to adjust LPs to mitigate river and reservoir siltation in some special locations. Overall, these results could benefit land use management for better soil and water conservation under the rapid development of southeastern China.

The effects of human activity on the results were not taken into consideration in our study. As mentioned before, construction engineering results in serious soil erosion events [79], and agricultural activity and industrial development consume many water resources [81,82]. A water use dataset should be collected to be used in relevant investigations in the future to reduce the influence of these human activities. In addition, it has been reported that the effects of various LPs on hydrological processes may have discrepancies in different seasons during a year [40,48]. The variations in hydrological processes characteristic in different seasons were not taken into consideration in our study due to the limitations of the hydrological data. A further study should be conducted to clarify these differences. Finally, though the watershed size in our study ranges from 1700 to 80,900 km2, the number of watersheds used in this study was relatively few. More watersheds with various size should be taken into consideration to obtain more stable correlations between various LMs and hydrological processes.

### **6. Conclusions**

In our study, the relationships between 14 LPs and 4 hydrological indices were investigated for 10 watersheds with various areas in the subtropical monsoon climate zone of southeastern China. The main conclusions are as follows:


More studies should be conducted in various climate zones to enhance the understanding of the relationships between LPs and various hydrological indices, and human activity and seasonal discrepancies should be taken into consideration.

**Author Contributions:** Conceptualization, C.W.; data curation, Y.M. and W.Z.; formal analysis, X.D., M.L., D.Y. and C.Z.; funding acquisition, X.D., Y.M. and D.Y.; investigation, C.W., W.Z. and M.L.; methodology, C.W.; software, D.Y.; supervision, X.D. and Y.M.; visualization, C.W.; writing—original draft, C.W.; writing—review and editing, X.D., Y.M. and B.S. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research has been supported by the Second Tibetan Plateau Scientific Expedition and Research Program (STEP, grant number 2019QZKK0103), the European Space Agency (ESA) and National Remote Sensing Center of China (grant number 58516), the Power Construction Corporation of China (grant number DJ-ZDZX-2016-02-09), the Research Fund for Excellent Dissertation of China Three Gorges University (grant number No. 2021BSPY002), and the National Natural Science Foundation of China (grant number No. 52109058).

**Data Availability Statement:** Not applicable.

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

### **Glossary**


### **Appendix A**

**Table A1.** The descriptions of the selected landscape metrics and relevant references.


**Figure A1.** The temporal variations of WY, RC, SYL, SYM, and SSC in the ZJ, DFK, HS, LJD, LX, BL, CA, SJ, ZQ, and WZ watersheds.

**Figure A2.** *Cont*.

**Figure A2.** The temporal variations of LPI, ED, PAFRAC, NP, DIVISION, AI, IJI, CONTAG, PD, LSI, COHESION, SHDI, SHEI, and MSIEI in the ZJ, DFK, HS, LJD, LX, BL, CA, SJ, ZQ, and WZ watersheds.

### **References**


**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.
