Study on the correlation between river network patterns and topography in the Haihe River basin

In recent decades, the river network patterns (RNPs) in China’s Haihe River basin have changed dramatically, and the topology of the river network has become increasingly complex. It is important to quantitatively study the correlation between river network patterns and topography (CRNPT) and the changes in the correlation. In this paper, the Haihe River basin was spatially gridded (4 km × 4 km), and different geomorphological areas were extracted for a multiarea study. We selected topographic and river network indicators and proposed new indicators to characterize regional topographic ‘stressfulness’ and then used redundancy analysis for correlation studies. The results showed that the variance of RNP explained by topography was 53.39%. The combined contribution of the topographic wetness index (TWI) and topographic wetness stress index (TSI) ranged from 35.66% to 78.29% in multiple areas, and the TSI showed stronger explanatory power. The regional effect of the CRNPT was significant, with mountains and transition areas having higher effects than plain areas. Compared to the natural river network, the CRNPT of the current river network was significantly lower. Among the RNP indicators, the artificial channel proportion ( P a c ) had the highest proportion of variance, and the CRNPT was strongly influenced by artificial channels. Artificial channels changed the consistency of topography with the RNP and reduced the topographic interpretation of the RNP, which may weaken the stability and hydrological connectivity of the river network. The variation in interpretation was related to the distribution of artificial channels, which showed a logarithmic function relationship between them.


Introduction
Rivers are multiscale systems that interact with climate, tectonics, erosion, sediment transfer patterns, and human activity (Carbonneau et al., 2012;Scorpio et al., 2018;Sarker et al., 2019;Walley et al., 2020), and their network branching models have a complicated topology (Barabási and Albert, 1999;Albert and Barabási, 2002). The river network is a fundamental feature of the basin and is a channel for horizontal effects (Heasley et al., 2019), influencing hydrological processes throughout the basin and regulating the biological population diversity and the spatial distribution of channel and riparian physical attributes (Benda et al., 2004). The river network is not a random topology (Perron et al., 2012), and its spatial pattern is determined by the physical properties of the catchment (Zanardo et al., 2013). It is critical to study the structural, geometric, and topological properties of the river network as well as to quantify the physical factors that shape differential river network patterns (RNPs) to better understand the dynamics of the river network in changing environments (Biswal and Marani, 2010;Shelef and Hilley, 2014;Abed Elmdoust et al., 2016).
In the landscape pattern of a catchment, the evolution of the topography and changes in the river network are strongly linked (Willett et al., 2014). Topography determines the path and residence time of water flow (Curie et al., 2007), which greatly influences the spatial pattern of the river network. Therefore, exploring the intrinsic links between river networks and topography is important for predicting the topological development of river networks and their stability, as well as for watershed planning. Chang. (1979) divided rivers into straight, meandering, and braided patterns and found that the channel pattern is straight when the channel slope is equal to the valley slope. However, when multiple channel slopes exist on a uniform valley slope, the river is sinuous. Alabyan and Chalov. (1998) analysed the stream power defined by river slope and discharge and noted that it is one of the most important factors in river development. The greater the stream power is, the stronger the branching tendency. Niemann et al. (2001) confirmed the importance of elevation in river network changes using river network growth models. Changes in elevation cause significant changes in the physical properties of the river network. In areas with gentle topography, the boundaries and water flow are less constrained, resulting in rivers that meander and twist (Piégay et al., 2005). Chen et al. (2019) subdivided China's Yangtze and Yellow River basins into 60 subbasins. The study found that slope had a significant positive relationship with the river network density, bifurcation ratio, and side branch ratio. Based on a study of 100 small watersheds in the United States, Sangireddy et al. (2016) found a strong positive correlation between river network density and topographic relief in wet environments and a weaker correlation in arid and semiarid environments. Zhao et al. (2016) discovered that tectonic uplift controls gully (channel) density, which is correlated with overall slope and vegetation cover, and that longitudinal slope and elevation integrals of the river network were correlated with gully density on China's Loess Plateau. Li et al. (2022) used correlation analysis, partial correlation analysis and multiple linear regression to analyse the quantitative relationships between topographic and river network parameters in the 18 subbasins of the Haihe River basin. The results showed significant correlations but strong spatial heterogeneity. The influence of topography on the formation and development of the river network varies from region to region and has obvious regional effects (Finnegan et al., 2005;Piégay et al., 2005;Tong et al., 2020).
These studies are mostly qualitative or small-scale, with no clear quantification of the correlation between river network patterns and topography (CRNPT). Furthermore, hydrological and geomorphological studies on river networks have long been based on networks with specific assumptions, such as that water always flows in the direction of the steepest downslope and that there is no ring structure in the network (Ni J R, 1998;Chen et al., 2019;Sarker et al., 2019). Most of these river networks are created using DEMs and are similar to natural dendritic river networks in terms of morphological features and topology (Vörösmarty et al., 2000;Cohen et al., 2018). They follow Horton's law (Horton, 1945) and the law of optimal channel networks (Rodríguez Iturbe et al., 1992;Abed Elmdoust et al., 2016) with the lowest total energy. However, there are significant differences between generated river networks and the actual "natural + artificial" river networks, particularly in basins with extensive floodplains, where the actual river network is a complex network structure subjected to multiple influences rather than a natural dendritic structure. Thus, the quantitative relationship between large watershed topography and actual network structure is still difficult to describe clearly.
The Haihe River basin is one of China's seven major river basins, and it has a significant role in politics, economy, and agriculture and a highly complicated river network. In recent years, China proposed and included in the 14th Five-Year Plan the "National Water Network Major Project," with the goal of speeding the building of the main skeleton and artery of the national river network. Watershed management based on river networks necessitates a greater scientific understanding of river morphology and function at multiple scales, as well as enhanced interdisciplinary studies of watershed topography, geomorphology, and hydrology (Heasley et al., 2019). This interdisciplinary approach is critical for understanding the complex relationships between factors in river systems (Dollar et al., 2007). In this paper, we gridded the Haihe River basin and extracted several topographic and geomorphic areas. Then, we used redundancy analysis to investigate the relationship between the actual "natural + artificial" river network and topography and to answer three important questions based on a thorough understanding of the RNP: 1) To what extent does topography explain the spatial distribution of the actual river network patterns? Do different areas of the basin behave in the same way? 2) Which topographic features have the most influence on river network patterns, and how can we understand the coupling effect between them? 3) Is the relationship between RNP and topography stable, and does the addition of artificial channels affect their consistency (i.e., how does the CRNPT change)? The findings of this paper can provide a basis for the study of the consistency of river networks and topography in large watersheds around the world.

Study flow
The aim of this paper was to quantify the correlation between RNP and topography and to compare the difference in the CRNPT between the current river network and the natural river network. Figure 1 shows the flow chart of the research in this paper. This paper used the Haihe River basin in China as the study area, and its main flow is as follows: 1) divide the spatial grid as a spatial sample; 2) determine the analysis area for comparative study; 3) select and calculate characteristic indicators, including topographical indicators and river network indicators; 4) conduct collinearity analysis on the characteristic indicators and remove those that are collinear; 5) use topographical indicators as explanatory variables and river network indicators as response variables for redundancy analysis; and 6) analyse the CRNPT of the current river network and the natural river network and discuss the differences in the results of the current river network and their reasons combined with the findings of previous studies.

Study area
The Haihe River basin, with a basin area of 320,600 km 2 , is in the eastern part of China's Bohai Sea region. It has a dense population and many large and medium-sized cities, with the Taihang Mountains and Yanshan Mountains in the upper reaches and the North China Plain, an important grain-producing region, in the lower reaches. The Haihe River basin is fan-shaped, with distinct mountain and plain boundaries and a trend of contraction from upstream to downstream. The upstream mountain area is rugged, while the downstream plain is gentle, and the topography and geomorphology are distributed spatially in a remarkable way. The Haihe River basin belongs to the temperate East Asian monsoon climate zone, with an average annual temperature of 0-14.5°C and an average annual relative humidity of 50%-70%. The average annual precipitation in the basin is 535 mm, which is representative of a semi-humid and semiarid zone. The average annual temperature and average annual precipitation basically decrease from southeast to northwest, but there are two extreme centres of annual precipitation in front of the windward slopes in the western and northeastern mountains. The density of the river network is 0.51, the natural length is 85,638.02 km, and the artificial length is 75,267.49 km, which indicates the complexity of the river network structure. Artificial channels are primarily constructed to serve flood control, irrigation, and water supply functions and have become an important component of the river network in the basin, causing drastic changes in the topology and spatial characteristics of the original natural river network, which has a dense mesh structure. The river network is very well connected, particularly in the vast downstream plain areas. Figure 2 provides a schematic diagram of the Haihe River basin's river network.

Topography and river data
The development of topographic data and remote sensing image data has aided in the quantification of structural and morphological characteristics of river networks (Benstead and Leigh, 2012), forming the foundation for the CRNPT study. The river network data in this paper were obtained from the National Catalogue Service for Geographic Information (https://www.webmap.cn/main.do? method=index). These data were formed by actual surveying and mapping and describe the actual river network in the Haihe River basin. The DEM and slope data were from the Geospatial Data Cloud (https://www.gscloud.cn/), and the climate and vegetation data were from the Resource and Environment Science and Data Center (https://www.resdc.cn/Default.aspx).

Topography and river network indicators
Many metrics, such as the Horton ratio (Horton, 1945), river network density (Horton, 1932;Abrahams, 1984), branching angle (Devauchelle et al., 2012), width function (Leopold, 1971), or ranking the degree of river branching by hierarchical methods such as Horton, Strahler, and others (Horton, 1945;Strahler, 1952;Tokunaga, 1984), can be used to characterize river network geometry and topology. However, in the Haihe River basin, where there are many artificial channels, the river network has very different characteristics than the natural network. Because the natural dendritic structure is not satisfied, the flow sequence in the reaches cannot be determined, and calculating the Horton ratio, branching angle, and other indicators characterizing the natural network is difficult. Finally, indicators such as network density (D r ), fractal dimension (D n ), junctions (N j ), and artificial channel proportion (P ac ) were chosen to characterize the river network. The analysis of RNP not only focuses on the network structure but also includes the network composition, which can more comprehensively reflect the characteristics of the river network and be integrated with basin management practice.
River network density is expressed as the ratio of total river length to total watershed area. River network density is a characteristic of the spatial distribution of rivers and can reflect Frontiers in Earth Science frontiersin.org 03 the density of regional water flow channels. It is calculated by the following formula (Horton, 1932): where D r is the network density (km/km 2 ); L i is the length of the river within grid i (km); and A is the area of the grid (km 2 ). Box dimension is the most basic and widely used characterization indicator of fractal dimension, and it can be used to study the scalar properties of complex structures describing river morphology and processes (Mandelbrot, 1982). Additionally, the box dimension of a river system can reflect the developmental stage of river erosion as follows (Hong, 1996): 1) a river is in the early (young) stage if the box dimension is less than 1.6; 2) a river is in the middle (mature) stage if the box dimension value is between 1.6 and 1.89; and 3) a river is in the late (old) stage if the box dimension value is between 1.89 and 2.0.

lgr
( 2) where D n is the box dimension; N (r) is the number of boxes with rivers; and r is the box side length (m). The box size should reflect the planform of the river, and its side length is determined by the degree of curvature of the river in the study area, taking multiple values from small to large (Mandelbrot, 1982). The calculation workload (total number of boxes) should also be taken into account. Finally, we used r values of 400 m, 500 m, 800 m, 1,000 and 2000 m based on the density and curvature of the rivers in the Haihe River basin.
Junctions are the number of intersections of the regional river network that can be obtained through spatial analysis and statistical calculation. Junctions characterize the degree of physical connectivity of regional rivers. The artificial channel proportion is the ratio of the length of the artificial channel to the total length of the river, and it indicates the composition structure of the network. The proportion of artificial channels can characterize the extent to which the river network has been modified and rebuilt by humans. The higher the human activity in the watershed is, the higher the artificial channel proportion. The calculation formula is as follows: where P ac is the artificial channel proportion; L i is the total length of the river in grid i (km); and L r,i is the length of the artificial channel in grid i (km). Not all topographic and geomorphological indicators were correlated with the RNP, so we chose some representative indicators based on the characteristics of the Haihe River basin, including elevation, slope, relief, topographic wetness index (TWI), elevation stress index (ESI), slope stress index (SSI), relief stress index (RSI), topographic wetness stress index (TSI), and farmland area (Farmland). Furthermore, the precipitation and moisture index Frontiers in Earth Science frontiersin.org (IM) were chosen to analyse the impact of climatic factors on the RNP. The climatic factor was not the focus of this paper, and it was included to assist in validating the results of the topography analysis (i.e., whether the results of climate and topography were consistent). Therefore, only two climatic indicators were selected. These indicators were chosen because they had a significant impact on the topology of the river network or the hydrological processes in the basin, as evidenced by many studies (Seibert and McGlynn, 2006;Zanardo et al., 2013;Kleidon, 2016;Zhao et al., 2016;Chen et al., 2019;Tong et al., 2020). Among them, elevation (m), slope (°), relief (m), Farmland (km 2 ), and precipitation (mm) were in standard units, and the TWI, ESI, SSI, RSI, TSI, and IM were dimensionless indicators. In this paper, the topographic indicators were extracted and calculated using ArcGIS 10.3. The topographic wetness index is a function of the regional catchment area and slope. It can characterize the effect of regional topography on runoff flow and storage and reflect spatially produced flow characteristics. The higher the value of the TWI is, the larger the flow-producing area and storage capacity are. It is calculated as follows (Beven and Kirkby, 1979): where α s is the upstream confluence area per unit contour length (km 2 ) and β is the slope (rad). The elevation stress index is an index proposed based on existing research. Riparian features and topographic constraints influence the range of channel stability (Piégay et al., 2005), and Tong et al. (2020) demonstrated that topography influences the spatial variations in channel width and slope, and nearby topographic constraints on river boundaries and water flow can determine the degree of channel curvature, which is well reflected by free-flowing meandering and braided channels on floodplains. The topography of the hillsides on both sides of the river therefore had a large limiting effect on the morphology and structure of the river, which could be reflected by the difference between the hillside elevation and the riverbed elevation. To reflect the limiting effect of topography, we defined the stress (or limiting) index by calculating the ratio of the terrain index without river pixels to the geomorphology index with river pixels based on the raster pixels. The higher the stress index is, the greater the limiting effects of the slope on the morphology and structure of the river are. The SSI, RSI, and TSI are all calculated in the same way: where SI is the stress index; T n is the topographic index value without river pixels; and T h is the topographic index value of the river pixels. Explanatory factors may have serious collinearity issues, thus increasing parameter variance estimation and causing large errors in the analysis results (Greene, 2003). To assess the degree of collinearity, the variance inflation factor (VIF) was used (O Brien, 2007). When VIF>10, it is widely assumed that there is a serious collinearity problem (Hair et al., 1992;Kennedy, 2008), and the indicators were selected again following this standard (Table 1). Figure 3 shows the spatial distribution of the final selected network and topographic indicators.

Spatial grids and areas
The Haihe River basin is a large area with a dense network of rivers. According to the topographic characteristics and RNP of the basin, and considering the calculation efficiency, it can be divided into 19,826 square grids with a spatial resolution of 4 km × 4 km. To investigate the local CRNPT as well as the differences, the Haihe River basin was divided into mountain, transition, and plain areas based on topographic distribution, and the windward and leeward sides of the Taihang Mountains and Yanshan Mountains were extracted as different spatial areas (Figure 4). Table 2 displays the statistical data for each spatial area.

Redundancy analysis
Multivariate statistical analysis is a powerful tool for investigating the interactions between hydrological elements and natural and human factors (Wang et al., 2017;Rakotondrabe et al., 2018). Redundancy analysis (RDA) (Rao, 1964) is a widely used form of asymmetric canonical analysis (Legendre et al., 2011). Redundancy analysis (RDA) examines the causes of variation in the response variable by explicitly modelling it as a function of the explanatory variables (Israels, 1984). Its goal is to determine the rate at which a given set of explanatory variables explains the variation (or dispersion) in the response variable. The greatest advantage of RDA is that this statistical method can independently maintain the contribution of each explanatory variable for each dependent variable without performing a simple analysis for the explanatory variable vector and by converting some of the variables into virtual complex variables. RDA has already been successfully implemented and is commonly used in geomorphic studies (Zhao et al., 2015), ecological studies (Wagner, 2004), hydrobiological studies (Obolewski and Strzelczak, 2009), and many others. The main steps of RDA are as follows (Makarenkov and Legendre, 1999;Mazouz et al., 2013): I) For each response variable Y, a multiple regression on all explanatory variables X is performed.
II) Calculate the covariance matrix S of the fitted value matrixŶ.
III) Eigenvalue and eigenvector decomposition is performed on the covariance matrix S.
where λ k is the regular eigenvalue and u k is the corresponding regular feature vector.
IV) The ranking of the samples (rows of Y) in the dependent variable space is obtained directly from the central matrix Y. The principal components are calculated using the standard equations.
Similarly, the ranking of the samples in the space of explanatory variables is obtained from the central matrix X.
Frontiers in Earth Science frontiersin.org

FIGURE 3
Spatial distribution of network and topographic indicators in the Haihe River basin.
Frontiers in Earth Science frontiersin.org 06 V) The contributions of the explanatory variables to the various canonical axes are estimated by calculating the linear correlations between the variables X and the ranking axes using Eq. 9 in Y-space and Eq. 10 in X-space.
In this paper, river network indicators were used as response variables, and topographic features were used as explanatory variables for redundancy analysis. Furthermore, the precipitation and moisture index were used as explanatory variables to compare the correlation between the RNP and climatic factors.
3 Results and discussion 3.1 Spatial distribution of river network and topography Figure 3 shows the spatial distribution of the RNP, topography, and climate indicators, while Figure 5 depicts the box plots for each indicator in various spatial areas. For the mountain-transition-plain areas, the mean values of D r were 0.42, 0.47, and 0.73, those of D n were 1.19, 1.22, and 1.35, those of N j were 0.49, 0.91, and those of 2.84 and P ac are 0.00, 0.45, and 0.72, respectively. The network indicator values increased gradually from the mountains to the plains. Under the influence of human activity, the river network in the plain area had a high density and many artificial channels. The topology of the river network in the plain area was complex, with a high degree of river development and physical connectivity. The greatest differences in spatial distribution were found for ) interval was small in the mountain area, as were the spatial differences in the river network. The spatial differences in river network topology were greater in the plains than in the mountains, and the regional river network varied more dramatically. Comparing the windward and leeward sides of Taihang and Yanshan, the mean values of D r were 0.41, 0.46, 0.40, and 0.40, those of D n were 1.18, 1.22, 1.17, and 1.17, those of N j were 0.49, 0.90, 0.44, and 0.49, and those of P ac were 0.0, 0.02, 0.0, and 0.0; additionally, the mean variances of the indicators were 0.028, 0.037, 0.026, and 0.029, respectively. The spatial differences in the river network leeward of Taihang were significant, while the differences in the other three areas were minor.
Along the mountain-transition-plain gradient, the mean elevation values were 1,016.10, 91.57, and 17.63, respectively. The ESI values were 1.11, 1.07, and 1.01; the SSI values were 1.41, 1.10, and 0.93; and the RSI values were 1.44, 1.19, and 0.97. The upstream mountain area was steep, and the topography had a greater restrictive effect on the river. The topography in the plain area was gentle, with the lowest restrictive effect. The topographic characteristics of the transition area were between those of the two areas. The TWIs were 7.39, 9.60, and 10.60, respectively, with the plain

Correlation analysis
The river network indicators had a negative correlation with the elevation, ESI, SSI, and RSI, a positive correlation with the TWI, TSI, and farmland, and an unstable correlation with the precipitation and IM ( Figure 6). The correlation coefficients of D r , D n , and N j with elevation, TWI, ESI, SSI, RSI, TSI, and Farmland were between 0.1 and 0.39, indicating weak correlations. The coefficients of P ac with the TWI and TSI were greater than 0.7, indicating strong correlations. P ac had a moderate correlation with elevation and farmland and a weak correlation with the ESI, SSI, RSI, and IM. The mean coefficient between P ac and environmental factors (including topography and climate) was 0.42, which was significantly higher than D r (0.25), D n (0.20), and N j (0.21). The intensity of human activity was influenced by topography, with more artificial channels in plain areas with gentle topography and fewer in mountainous areas. The distribution of artificial channels reflected the spatial nature of the topography. The TWI and TSI had correlation coefficients of 0.49 and 0.45 with network indicators, which were greater than those of elevation (−0.39), ESI (−0.22), SSI (−0.17), RSI (−0.18), and Farmland (0.38), indicating that regional runoff flow and storage capacity had a greater influence on the RNP.
Correlations were weaker in the local area than in the whole basin. The D r , D n , N j , and P ac correlation coefficients with environmental factors in mountain areas were 0.23, 0.17, 0.18, and 0.24, respectively. D r , D n , and N j were not correlated with environmental factors in the transition area, and the correlation coefficient between P ac and environmental factors was 0.31. The correlation between the river network and environmental factors was relatively weak in the plains. In the mountain area, topography

Frontiers in Earth Science
frontiersin.org played a greater role in influencing the RNP, and the correlation was higher than it was downstream. The correlation coefficients between network and environmental factors in Taihang were 0.17, 0.12, 0.13, and 0.19 for the windward side and 0.31, 0.24, 0.25, and 0.31 for the leeward side, with the leeward being higher than the windward. Climate had a greater influence on the leeward side of the Taihang than on the windward side. There were small differences in the correlation coefficients between network and environmental factors Heatmap of Pearson's correlation coefficient (* means significance test p < 0.05)

FIGURE 5
Box plots of network indicators, topographic indicators, and climate indicators in various areas.
Frontiers in Earth Science frontiersin.org 09  Table 3 shows the results of the redundancy analysis, where the p values (>F) for the first two significant axes, RDA1 and RDA2, were both 0.001, indicating a significant correlation between the river network reflected by the two ranking axes and the topography. The variance explained by the topography on the network indicators was 53.39%. The elevation, TWI, and TSI contributed 21.46%, 25.45%, and 30.07%, respectively. The coefficients of the network ranking axis determine the contribution, and the explanatory variables in the redundancy analysis are based on the variance of the response variables. P ac had the highest variance (0.169), and the coefficient's absolute value in the first ranking axis RDA1 (variance explained 52.68%) was 4.75 (Table 4), so the TSI, which had the highest correlation with P ac (0.73), was the indicator with the largest contribution. Artificial channels changed the natural network pattern, and P ac became the most interpretative network indicator, which also overshadowed the interpretation of other indicators. There was a negative correlation between elevation and the river network. The RNP in the plain area was severely affected by artificial channels because human activities were primarily concentrated in the gentle terrain area (Figure 2). The TWI and TSI were most strongly correlated with the RNP; the TWI reflected the regional stream power, and the TSI represented the corresponding physical restrictiveness. This is an important factor in determining the RNP (Chang, 1979;Alabyan and Chalov, 1998). The efficiency of kinetic power dissipation in the  channel is reduced in areas with large catchment areas and small slopes, and the more energetic flow during a flood may remove its channel in locations where the restriction effect is weak, promoting the lateral development of the river network and the increasing density of the network. The human renovation of the river network for flood protection is also based on this phenomenon. The ESI, SSI, and RSI are indicators that characterize only topography compared to the TSI, which contains the definition of stream power. The SSI describes the limiting effect of slope but cannot define the water flow, and differences in values did not cause a strong response in the RNP (weak correlation). As a result, the impact of the ESI, SSI, and RSI on the river network was minor. People have been excavating diversion irrigation channels in areas with a high concentration of Farmland and are constantly changing the river network structure, which causes Farmland to be correlated with RNP. Table 5 The variance in the river network explained by topography was 28.25% in the mountain area, 30.13% in the transition area, and 9.59% in the plain area. The contribution rates of the TWI, TSI, and Farmland in mountain areas were 31.82%, 32.96%, and 18.05%, respectively. P ac had the largest coefficient in the ranking axis, so the TSI had the largest contribution rate. The river network and TSI distribution in the mountain area are depicted in Figure 7A. The artificial channels were spatially consistent with the distribution of high TSI values, and the alignment of artificial channels was always along the narrow area of high TSI values. In the transition area, D r , D n , and N j were not correlated with topography, and the variance explanation of the RNP mainly relied on P ac . The correlation coefficient between the TSI and P ac was 0.64, and the contribution of the TSI was 69.20%, making the explanation rate of the transition area greater than that of the mountain and plain areas. Figure 7B shows the river network and TSI distribution in the transition area. The stream power constraint was large upstream of the transition area, and the river network was dominated by natural channels. As the constraint decreased downstream, the artificial channels increased, and the composition structure of the river network changed, which was consistent with the distribution of the TSI. In the plain area, the numerical interval of topography varied very little compared to the river network ( Figure 5). The correlation coefficients between network indicators and topography were low, indicating that the consistency between the regional actual network and topography in the region was poor.

Redundancy analysis
The variances in the river network explained by topography were 21.18% and 33.48% for the windward and leeward sides of Taihang, respectively. The CRNPT was greater leeward than windward. The variances in the river network explained by topography were 16.79% and 13.21% for the windward and leeward sides of Yanshan, respectively. The CRNPT was greater windward than leeward, but the difference between the two was small. The river network and topography on the windward and leeward sides did not follow a consistent pattern, and the differences in correlation were dependent on the respective regional characteristics. Taihang is located inland, trends northeast to southwest and is influenced by the southeastern monsoon and geomorphological evolution, with significant topographic and climatic differences between the windward and leeward sides. The numerical interval of the river network and topographic indicators on the leeward side is large, and the influence of topography increases the spatial variation in the river network.
On the windward side, the internal differentiation is reduced, resulting in a larger difference between the CRNPT on the windward and leeward sides. Yanshan has a smaller spatial extent and elevation than Taihang, as well as a different geographic location and mountain range orientation. The Yanshan is oriented east-west and is close to the Bohai Sea, with similar overall topography and climate, so the CRNPT differences were minor.
The network variance explained by climatic factors in all areas did not exceed 15%, and the value was 10.62% for the whole basin, with the correlation coefficients for climatic and river network indicators being low. Existing research in other regions found that climate has a large influence on river network geometry and topology (Abrahams, 1984;Zanardo et al., 2013;Seybold et al., 2017;Ranjbar et al., 2018;Chen et al., 2019) and that there was a correlation between them, but no significant effect was found in the Haihe basin. The inconsistency is mainly due to the following factors. First, the choice of response variables and statistical analysis results were dependent on the variables analysed, and the correlation between different indicators varied significantly. This view is demonstrated in the study of Zanardo et al. (2013), where the indicators characterizing the network branching may be more strongly correlated with climate. Second, it may be related to the analysis area, where the climatic differences within the Haihe River basin are small, resulting in a weak correlation between the RNP and climatic factors. The difference in precipitation was only 22.04 mm, and the difference in IM was 5.77 on the mountain-plain gradient with a significant distribution of the RNP, whereas studies by Seybold et al. (2017), Ranjbar et al. (2018), and others were based on the entire U.S. or larger areas, with significant differences in climate. Finally, the distribution of artificial channels diminished the consistency between climate and river networks.
The above results show that there is a significant regional effect of topography on the interpretation of the RNP. On the one hand, there were differences in linear correlations between indicators in different areas; on the other hand, there were differences in the interpretation of each river network indicator. When multiple indicators characterized the RNP, those with large variance and more correlation with topography had a stronger role in determining the explanatory rate, as is the case with P ac in multiple areas. The variance of P ac was 0.169 (80% of the total variance), and the correlation coefficient was 0.49 with topography in the whole basin. The variation in the indicators decreased in the local area, and the correlations between river network indicators and topography were all weaker than those in the whole basin. The P ac variance and proportion in the mountain, transition, and plain areas were 0.056 (0.58), 0.167 (0.73), and 0.103 (0.66), respectively. Therefore, the explanatory rate was greater in the whole basin than in local areas. The P ac variance and its windward and leeward proportions in Taihang were 0.045 (0.39) and 0.086 (0.58), respectively, and the correlation coefficients were 0.21 and 0.32, respectively. The P ac variance and proportion of windward and leeward Yanshan were 0.02 (0.19) and 0.018 (0.15), respectively, and the correlation coefficients were 0.15 and 0.15, respectively. As a result, the explanation rate of the leeward side of Taihang was significantly higher than that of the Frontiers in Earth Science frontiersin.org windward side, and the difference between the windward and leeward sides of Yanshan was very small. Although there was a regional effect within the basin and a low explanation rate in the local areas, the combined contributions of the TWI and TSI in each area were 55.52%, 64.78%, 78.29%, 35.66%, 59.68%, 62.12%, 52.12%, and 51.48%, respectively. In fact, the power dissipation of water flow, as described by the TWI, is also an important component of channel erosion forces (Thorne et al., 1997), and it has been demonstrated that a dense branching network is continuously and significantly influenced by erosion forces (Perron et al., 2012). Both the TSI and the riverbanks are important factors that interfere with channel flow movement, which can limit or facilitate the cutting of side slopes and encroachment on nearby channels, affecting the lateral stability of the channels. The TWI is also an important geomorphological indicator for characterizing hydrological properties, as it reflects the spatial distribution of soil moisture, surface saturation, and the runoff generation process (Beven and Kirkby, 1979). Hydrological similarity exists between two points of equal TWI, indicating an inherent relationship between the RNP and hydrology. Slope, as an intrinsic factor of the TWI and TSI, has a strong influence on the RNP and hydrological processes, and the river network develops more rapidly in areas with low slope and good hydrological conditions. The TWI and TSI are the most strongly correlated indicators with the river network and are consistent with the river network pattern in different areas.

Correlation between river network patterns and topography influence factor analysis
The RNP in the Haihe River basin has complex spatial characteristics, and the regional effect of the CRNPT was

FIGURE 8
Correlation of the natural river network with topographic and climatic factors (* means significance test p < 0.05) Frontiers in Earth Science frontiersin.org significant. However, the correlation between the topography and actual RNP was not very strong in multiple areas, which was inconsistent with some previous studies. Several studies have revealed the potential relationship between the RNP and topography, involving river networks of different characteristics in Asia, Africa, and the Americas. As mentioned in the introduction, significant relationships between numerous topographic factors (slope, relief, tectonic movement, etc.) and river network indicators (river network density, branching ratio, confluence angle, etc.) have been demonstrated by statistical, geomorphological, and remote sensing methods (Perron et al., 2012;Molin and Corti, 2015;Sangireddy et al., 2016;Chen et al., 2019;Beeson and McCoy, 2022;Li et al., 2022). Furthermore, results from small-scale studies in the laboratory have shown a correlation between water flow channels and topography (Darboux et al., 2002).
The inconsistency in the conclusions may be due to differences in the scale of the study (grid size) and especially the type of river network data (actual river network was used in this study, natural river network was used in previous studies). There are numerous artificial channels in the Haihe River basin, whose combined length is close to that of natural rivers (0.88:1). The variance proportion of P ac in the whole basin, mountain area, transition area, plain area, and Taihang leeward area was greater than 50%, which was the most interpretative indicator of the river network. The correlation between the topography and actual RNP was strongly influenced by the distribution of artificial channels, but we did not have a clear understanding of their extent. Therefore, this paper continues to explore the response of the CRNPT to changes in the river network and clarifies the impact of artificial channels on the quantitative natural correlation between the RNP and topography. First, to compare the effect of artificial channels on the spatial distribution of the river network, Moran's I was used to assess the spatial distribution of the RNP (Fortin, 1999;Haining and Haining, 2003;Wagner and Fortin, 2005). Moran's I is a spatial autocorrelation measure that generally ranges from −1 (discrete) to 1 (clustering), with values close to 0 indicating spatially random data (Qi and Wu, 1996;Fortin, 1999). To investigate the impact of artificial channels, the indicator P ac was removed, and only three indicators, D r , D n , and N j , were used to characterize the RNP. Table 6 compares the Moran's I of the natural and actual networks in each area, with the mean values of the indices as follows: whole basin (0.3, 0.5), mountain area (0.19, 0.33), transition area (0.32, 0.32), plain area (0.28, 0.55), Taihang windward (0.15, 0.21), Taihang leeward (0.28, 0.35), Yanshan windward (0.06, 0.07), and Yanshan leeward (0.14, 0.16); these results indicate a positive spatial autocorrelation of the RNP. Due to the needs of flood control, irrigation, etc., numerous artificial channels are concentrated in areas with gentle topography, which increases the autocorrelation of the RNP in each area, and the spatial distribution is more aggregated, reducing the spatial heterogeneity of the RNP. There are fewer artificial channels in areas with rugged topography and little change in spatial distribution. The correlation results showed that artificial channel construction changed the original intrinsic relationship between the natural river network and the topography. The natural river network is positively and weakly correlated with the elevation (0.29), ESI (0.34), SSI (0.17), and RSI (0.22) and negatively and weakly correlated with the TWI (−0.28) (Figure 8). After combining the artificial channels, the actual river network is negatively and weakly correlated with the elevation (−0.31), ESI (−0.19), SSI (−0.14), and RSI (−0.15) and positively and weakly correlated with the TWI (0.42). The trend is consistent across all areas, and the artificial channels cause a positive or negative shift in the weak correlation between the RNP and the topography, but the difference in correlation degree is not significant. This shift in

FIGURE 9
Comparison of variance explained by the RNP; the river network indicators are D r , D n , and N j Frontiers in Earth Science frontiersin.org correlation is progressively promoted with the construction of artificial channels. The RNP is constantly changing, and its correlation with topography is also changing towards the same trend.
The artificial channels altered the correlation between the natural river network and the topography, and Figure 9 depicts the comparative results of the redundancy analysis. The explained rates of natural and actual river network in the whole basin were 73.69% and 17.96%, 32.05% and 15.94% in mountain area, 76.51% and 1.14% in transition area, 82.27% and 14.32% in plain area, 23. 65% and 11.19% in Taihang windward,52.71% and 13.91% in Taihang leeward,19.44% and 16.13% in Yanshan windward, and 13.05% and 13.96% in Yanshan Leeward. The construction of artificial channels significantly reduces the correlation between the original natural network and topography, and the length of the artificial channels increases the variations in the explanation rate. Additionally, the relationship between the two showed a logarithmic function (R 2 =0.87). Therefore, the spatial stability of the CRNPT was poor and influenced by the disturbance effect of human activities, and the construction of artificial channels will cause the CRNPT to change continuously. The consistency between the natural river network and topography is reduced or even reversed, making its topography significantly less capable of explaining it. The natural network is generally dense in the upstream higher terrain areas and sparse in the downstream plain area (Chen et al., 2019). The Haihe River basin has completely changed the natural RNP due to the concentrated construction of quantities of artificial channels in local areas, raising the plain area D r higher than the mountain area and effectively increasing the spatial autocorrelation, but changing the original natural CRNPT, and the correlation continues to decrease. This is an unavoidable trend in the basins' downstream area to ensure flood control, irrigation, and water supply.
The variation in the CRNPT reflected changes in consistency. Topography was the determining factor of hydrological processes, and the river network pattern matched the hydraulic gradient to ensure the healthy operation of the network. The consistency of the river network with the topography in its natural state is conducive to ecosystem health and biodiversity, while changes in consistency can affect the stability of the river network ecosystem. The river network has evolved from a natural topographic feature to a natural-artificial dual feature under human influence (Sun et al., 2021), which can increase the network connectivity by providing routes around natural barriers and other water conveyance channels (Cote et al., 2009). However, the addition of a large number of artificial channels has severely fragmented the flow of the river network, leaving the lower river network without sufficient flow to maintain a healthy ecosystem (and the disappearance of large floods), ultimately severely disrupting the hydrological connectivity and habitat health of the network (Xingyuan et al., 2023). Moreover, multiple modifications and reconstructions within the river network may cumulatively impair ecologically relevant flow characteristics (Merenlender et al., 2008;Deitch et al., 2009). The serious problems facing the river network in the lower reaches of the Hai River basin also validate this conclusion. The lower reaches of the Hai River basin have serious ecological problems, with the river network often in a state of disconnection and little water flowing out of the estuary, resulting in the severe shrinkage of many estuaries (Hu et al., 2001).
The water quality of the lower river network is poor, and sediments contain high levels of pollutants (Liu et al., 2014). Moreover, the Haihe River basin faces threats such as habitat degradation and biodiversity reduction (Chen et al., 2012). The changing river network and the reduced consistency are continuously driving these problems. This inconsistency will be a long-term state, as human activity in the basin greatly limits the dynamics of the river network derived from natural factors such as topography and climate.
The renovation and reconstruction of large river networks is always based on single purposes and lacks systematic planning and assessment, leading to an evolutionary imbalance in the structure and function of the river network. We revealed the disruption of the correlation between the RNP and topography by artificial channels and its consequences. Basin management should pay attention to changes in the physical characteristics of the river network and runoff processes caused by the instability between the RNP and topography and the impact of spatially concentrated construction of artificial channels on the RNP. As a large catchment of 320,000 km 2 , we explored the relationship between the RNP and topography based on a gridding approach. In China, similar changes in river networks are very common in large watersheds (Yangtze River, Taihu Lake, etc.) (Deng et al., 2018;Song et al., 2019). The human influence on river networks in watersheds is so great that large-scale analysis can reveal the intrinsic relationships at the watershed level as well as the changes in such relationships, and this approach can be better generalized for use in other watersheds.

Conclusion
This paper examined the potential relationship between the current river network and topography, taking the Haihe River basin as an example, and revealed the change law of the CRNPT of the basin, as well as the influence of artificial channels. This research can serve as a foundation for watershed river network planning and management. The main conclusions are as follows.
The mountainous terrain is rugged, and the river network is sparse, with natural rivers dominating. The plain area has gentle topography and a dense river network, and it has 72% of the artificial channels. The variance explanation of topography on the RNP was 53.39%. With variance explanations of 28.25%, 30.13%, and 9.59% for mountain, transition, and plain areas, respectively, the regional effect of the CRNPT was significant, indicating that the more complex the topography is, the more the RNP is affected. The RNP and topography in the windward and leeward directions did not follow a consistent pattern.
The elevation, TWI, and TSI had the highest correlation with the RNP, and their contributions to variance explanation were 21.46%, 25.45%, and 30.07%, respectively, which had the greatest influence on the RNP. The combined contribution of the TWI and TSI in each area ranged from 35.66% to 78.29%. The TWI characterized the spatial water flow power dissipation and soil hydrology, and the TSI characterized the limiting effect of regional topography on hydrology and runoff processes, which demonstrated a stable correlation with the evolution of the RNP.
Artificial channels played the most important role in the CRNPT. The concentrated construction of artificial channels