^{1}

^{*}

^{2}

^{*}

^{1}

^{*}

Application of fractal theory on the evolution of nonlinear study of the hydrological system, which found its internal rules from the complex hydrologic system, could make us more fully understand the hydrodynamic characteristics of the complex motion of this system. Taking Weihe River as study area, this paper analyzes daily runoff series’ multi-fractal character and relative fluctuation feature by using the De-trended Fluctuation Analysis ( DFA) method. Result shows that the daily runoff series of main channel and branches of Weihe river all shows multi-fractal characteristics clearly, and the turns of multi-fractal intensity of daily runoff series in Weihe river are: Xianyang station (1.388) > Yingge station (0.697) > Linjiacun station (0.665) > Zhangjiashan station (0.662) > Zhuangtou station (0.635). Rainfall, evaporation, water income and human activity and other factors affect the fluctuation character and multi-fractal intensity of daily runoff series through these factors’ superimposition and pining down mutually. This study could provide a theoretical supply for obtaining the quantitative indicators on multi-fractal characteristics about eco-environment situation of watershed, and for runoff forecasting.

Runoff is one of key links of water cycle, its form and development affected by the comprehensive factors which include the climate, the terrain landform, soil and vegetation, so, describing it with multi-scales tool could be better, it may be used to well describe the fluctuation distributed condition of different partial conditions which influence to the entire series, and to analyze the function of influence factor to the small fluctuation. At present, in the outside and inside China, researches on runoff fractal mainly focus on the flood fractal characteristic. When Gupta and Waymire (1990) carried on the region analysis of flood, they successfully introduced the conception of scale invariable supposition, thus found a new way for the flood research [

The multiple fractal analysis (MF analysis) may calculate multiple fractal spectra (MF spectrum), including the sole fractal dimension. Therefore, in carrying on analysis to the runoff succession’s variation process, compared to sole fractal analysis, multiple fractal analysis may provide more information, and have the superiority in describing the non-homogenizing process. The correlation research indicated that daily runoff series can manifest the fluctuation degree of runoff process better. It is necessary to do the analysis on multiple fractal characteristics of daily runoff series of the branches and main channel of Weihe River, and further discuss the relationship between daily runoff process’ multiple fractal and its fluctuation, provides a solid theoretical foundation for the runoff research and the forecast.

Multiple fractal can also be called multiple scale fractal or multiple fractal measurement, multiple fractal may define initially as:

Supposes R d is a d dimension Euclidean space or Metric space, X is a d dimension subset of R d , and is also a measurement subset, or an invariant set of a dynamic system. To do appropriate division to X and entrusts with the invariable measure μ 0 , supposes α is a related parameter closely with division, after n steps’ division, an obtained subset could be defined as X n ( α ) , if X n = lim n → ∞ X n ( α ) is a fractal set, then it could be called ( X , μ ) fractal subset. Under this division, fractal subsets produced from ( X , μ ) may express as co-fractal subsets of many fractal subsets, and each fractal subset has different fractal dimension, then this fractal set is multiple fractals.

Non-tendency fluctuation analysis method (DFA) is a time series’ long period relevant method, Peng et al. (1994) proposed it when they research riptide behavior’s machine-made process [

Concrete analysis process (Yuan Ying et al. 2007; Xie Xianhong et al. 2008) of Multiple fractal non-tendency fluctuation analysis method (MF-DFA) is shown as follows [

a) Accumulated deviation computation of time series. To a given time series x t ( t = 1 , 2 , ⋯ , N ) , the first operation is accumulated and transforms this series into a new series Y ( i ) :

Y ( i ) = ∑ t = 1 i ( x ( t ) − x ¯ ) , i = 1 , 2 , ⋯ , N (1)

Here x ¯ is the average value of time series x ( t ) ;

b) Former transformed new series Y ( i ) has been divided into n new subinterval opening with same steps, here n = int ( N / s ) . As time series length N maybe not a multiple time larger than the time length s, in order to enable the rear part remainder data could be effectively used, uses the similar method to deal with the reverse order series Y ( i ) , therefore may obtain 2n equal length subintervals.

c) Through fitting, each sub-sector’s ( v = 1 , 2 , ⋯ , 2 n ) partial tendency could be obtained. Inner the subinterval, use least squares method to do k steps multinomial data fitting on s observation value, furthermore to look for the most superior fitting multinomial, then we can get partial tendency P v ( i ) .

P v ( i ) = a 0 + a 1 i + a 2 i 2 + ⋯ + a k i k (2)

Here: a_{i} is the multinomial fitting coefficient ( i = 1 , 2 , ⋯ , s ) ; k is the most higher number of times of multinomial fitting, when k is 1, 2, 3, … We call them separately DFA1, DFA2, DFA3, ….

d) To each sub-sector v, utilizes the multinomial come from fitting to eliminate the sub-partial tendency, then calculates each sub-sector’s mean-square deviation to eliminate the tendency:

{ F 2 ( v , s ) = 1 s ∑ i s { Y [ ( v − 1 ) s + i ] − P v ( i ) } 2 , ( v = 1 , 2 , ⋯ , n ) F 2 ( v , s ) = 1 s ∑ i s { Y [ N − ( v − n ) s + i ] − P v ( i ) } 2 , ( v = n + 1 , n + 2 , ⋯ , 2 n ) (3)

e) In order to enlarge or deduce the pantograph variance size to emphasized the fluctuation degree, to all same length sectors, look for their average q ( q ∈ R ) steps DFA waving function

{ F q s = 1 2 n ∑ v = 1 2 n { [ F 2 ( v , s ) ] q / 2 } 1 / q , q ≠ 0 F q s = exp { 1 4 n ∑ v = 1 2 n ln ( F 2 ( v , s ) ) } , q = 0 (4)

f) Determination of wave function scale coefficient. In each step q’s double logarithmic diagram ( s , F q s ) , inspect wave function’s scale behavior. If the long distance power law of primitive series { X t } is related, then here lies the power law relation between wave function F_{q}s and the criterion s, namely:

F q s ∼ s h ( q ) (5)

To each division length s, may extract corresponding fluctuation function value F_{q}s, take slope ratio h(q) that obtained from the linearity fitting value through least squares method, as q steps generalized Hurst coefficient. When h(q) is a constant and independent from q, the series shows as a sole fractal; When h(q) is the function on q, the series should be the multiple fractal.

Another parameter could be used to mark the multiple fractal is the strange spectrum. Through MF-DFA, generalized Hurst coefficient h(q) and Renyi coefficient τ(q) draw from multiple fractal formula system could be obtained. And following relation between the strange coefficient α and strange score f ( α ) exist here:

τ ( q ) = q h ( q ) − 1 α ( q ) = ∂ ∂ q τ ( q ) = h ( q ) + q h ′ ( q ) f ( α ) = q α − τ ( q ) = q [ α − h ( q ) ] + 1 (6)

Weihe River is the largest tributary of the Yellow River, China . It originated in the northern side of the Niaoshu Mountain, Gansu Province, the area of the watershed is 135,000 km^{2}, length of the main channel is 818 Km, of which Shaanxi Province accounted for 61.4%. The elevation of the watershed decreases from west to east, and changes from 329 m to 3495 m. There are four main types of landforms in the basin: loess hills, loess plateaus, river valleys basin and Qinling mountains. The tributaries of the Weihe River are dense, Jinghe River and the Beiluohe River are the first and second major tributaries of it, both lie in the north area of the watershed. Other north tributaries include: Xianhe River , Sandu River , Hulu River , Niutou River , Qianhe River , Qishuihe River , Shichuanhe River , et al. On the south side, the larger tributaries include, the Bangshahe River, Jihe River, Shitouhe River, Heihe River, Fenghe river, and Bahe River (

Selects Linjiacun station, Xianyang station and the Huaxian station of the main channel, and Zhangjiashan station of Jinghe River (largest branch of Weihe River), Zhuangtou station of Beiluohe river and Yingge station of Shitouhe river as typical stations, use their daily runoff data series and through MF-DFA, to analyse their multiple fractal characteristic, at the same time, to discusses the fractal relationship between the branches’ daily runoff process and main channel’s daily runoff process, thus provide the data support for the complexity research of Weihe River basin’s runoff process.

Among these stations, 58 year’s daily runoff data (1944-2001) of Linjiacun, Zhangjiashan and Zhuangtou are selected, 23 years (1986-2008) of Xianyang station, Huaxian station, 48 years (1961-2008), Yingge station of Shitouhe river is 34 years (1974-2007).

Using MF-DFA method, to do the multiple fractal analyses on the daily runoff series of Linjiacun, Xianyang and Huaxian station. According to the requirement of multiple fractal research, determine the value of s as 12, thus values of q are −12, −11, ∙∙∙, 11, 12. To the exponent number k, takes it as k = 2, 3, 4, 5, then calculate Hurst coefficient h ( q ) separately on three hydrologic station’s daily runoff series. Results are shown in

To MF-DFA of 2 to 5 steps, the daily runoff process’ Hurst generalized coefficient h ( q ) of Linjiacun, Xianyang and Huaxian station are all not the constant number, also all changes along with q’s changing (see

q | Lingjiacun station | Xianyang Station | Huaxian station | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|

k = 2 | k = 3 | k = 4 | k = 5 | k = 2 | k = 3 | k = 4 | k = 5 | k = 2 | k = 3 | k = 4 | k = 5 | |

−12 | 1.58 | 1.583 | 1.586 | 1.641 | 1.973 | 2.067 | 2.099 | 2.148 | 1.386 | 1.462 | 1.524 | 1.601 |

−11 | 1.571 | 1.575 | 1.578 | 1.634 | 1.965 | 2.059 | 2.091 | 2.14 | 1.379 | 1.455 | 1.518 | 1.595 |

−10 | 1.56 | 1.565 | 1.569 | 1.625 | 1.955 | 2.049 | 2.082 | 2.13 | 1.371 | 1.448 | 1.511 | 1.587 |

−9 | 1.546 | 1.553 | 1.558 | 1.615 | 1.943 | 2.037 | 2.07 | 2.118 | 1.362 | 1.439 | 1.502 | 1.578 |

−8 | 1.528 | 1.538 | 1.544 | 1.601 | 1.928 | 2.022 | 2.056 | 2.104 | 1.35 | 1.428 | 1.491 | 1.567 |

−7 | 1.506 | 1.519 | 1.526 | 1.584 | 1.908 | 2.003 | 2.038 | 2.086 | 1.335 | 1.414 | 1.478 | 1.553 |

−6 | 1.477 | 1.494 | 1.502 | 1.561 | 1.883 | 1.979 | 2.015 | 2.062 | 1.316 | 1.397 | 1.461 | 1.535 |

−5 | 1.437 | 1.458 | 1.469 | 1.529 | 1.848 | 1.945 | 1.982 | 2.03 | 1.291 | 1.374 | 1.439 | 1.511 |

−4 | 1.379 | 1.406 | 1.421 | 1.481 | 1.797 | 1.896 | 1.936 | 1.984 | 1.256 | 1.342 | 1.408 | 1.478 |

−3 | 1.294 | 1.329 | 1.349 | 1.408 | 1.717 | 1.82 | 1.865 | 1.914 | 1.205 | 1.296 | 1.362 | 1.429 |

−2 | 1.168 | 1.214 | 1.241 | 1.295 | 1.584 | 1.694 | 1.749 | 1.799 | 1.131 | 1.226 | 1.292 | 1.355 |

−1 | 1.011 | 1.064 | 1.096 | 1.143 | 1.374 | 1.491 | 1.557 | 1.608 | 1.026 | 1.119 | 1.18 | 1.237 |

0 | 0.857 | 0.903 | 0.935 | 0.978 | 1.131 | 1.23 | 1.292 | 1.333 | 0.897 | 0.968 | 1.021 | 1.068 |

1 | 0.739 | 0.774 | 0.804 | 0.845 | 0.934 | 1.016 | 1.07 | 1.104 | 0.784 | 0.839 | 0.886 | 0.928 |

2 | 0.666 | 0.693 | 0.723 | 0.764 | 0.806 | 0.887 | 0.943 | 0.975 | 0.706 | 0.755 | 0.802 | 0.844 |

3 | 0.62 | 0.642 | 0.671 | 0.712 | 0.718 | 0.804 | 0.864 | 0.896 | 0.648 | 0.697 | 0.746 | 0.788 |

4 | 0.589 | 0.607 | 0.635 | 0.676 | 0.655 | 0.747 | 0.81 | 0.843 | 0.605 | 0.655 | 0.705 | 0.748 |

5 | 0.567 | 0.581 | 0.609 | 0.65 | 0.611 | 0.706 | 0.773 | 0.806 | 0.572 | 0.624 | 0.674 | 0.718 |

6 | 0.552 | 0.562 | 0.589 | 0.629 | 0.578 | 0.677 | 0.746 | 0.779 | 0.546 | 0.6 | 0.651 | 0.695 |

7 | 0.54 | 0.547 | 0.573 | 0.614 | 0.554 | 0.655 | 0.725 | 0.759 | 0.526 | 0.581 | 0.633 | 0.678 |

8 | 0.53 | 0.535 | 0.56 | 0.601 | 0.535 | 0.638 | 0.709 | 0.744 | 0.51 | 0.566 | 0.619 | 0.664 |

9 | 0.522 | 0.525 | 0.55 | 0.591 | 0.52 | 0.625 | 0.697 | 0.731 | 0.496 | 0.554 | 0.608 | 0.653 |

10 | 0.516 | 0.517 | 0.542 | 0.582 | 0.508 | 0.614 | 0.687 | 0.721 | 0.485 | 0.544 | 0.599 | 0.644 |

11 | 0.51 | 0.51 | 0.535 | 0.575 | 0.498 | 0.605 | 0.678 | 0.713 | 0.476 | 0.536 | 0.591 | 0.636 |

12 | 0.505 | 0.505 | 0.528 | 0.569 | 0.49 | 0.598 | 0.671 | 0.706 | 0.468 | 0.529 | 0.585 | 0.63 |

Here: MF-DFA is called as k steps method.

similar non-linear relation, and h ( q ) decreased along with increasing of q, thus shown as a digression function on q, namely, we can determine that the daily runoff process of Linjiacun, Xianyang and Huaxian station of Weihe River’s main channel has multiple fractal characteristic.

Determination of multi-fractal of daily runoff process of Weihe River’s tributaries, From

Making the k step multinomial fitting to each sector’s s datum points of every runoff series of branches in Weihe River, aims of doing this lies on eliminating the k step tendency fluctuation from the accumulation deviation series y ( i ) , and eliminating the k-1 step tendency fluctuation in the primitive series. Result shows in

In general, through analysis on the main channel and branches of Weihe river, when using different steps MF-DFA to estimate runoff series’ generalized Hurst coefficient h ( q ) , its value changes in different sizes, but the general tendency is same.

Above research told us that, daily runoff series of main channel and branches of Weihe river all have multiple fractal feature, and all coefficients show the same changing rule. So, when carrying on analysis to multiple fractal intensity of runoff series of branches and main channel, it is a better choice to use daily runoff series’ k = 2 step multinomial (after eliminated tendency) to do contrast analysis.

As to daily runoff series’ multiple fractal intensities, we can generally use several

q | Zhang Jiashan Station | Zhuangtou Station | Yingge Station | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|

k = 2 | k = 3 | k = 4 | k = 5 | k = 2 | k = 3 | k = 4 | k = 5 | k = 2 | k = 3 | k = 4 | k = 5 | |

−12 | 0.843 | 0.898 | 0.975 | 1.063 | 0.83 | 0.892 | 0.95 | 1.047 | 1.002 | 1.114 | 1.233 | 1.356 |

−11 | 0.837 | 0.891 | 0.968 | 1.056 | 0.826 | 0.886 | 0.944 | 1.04 | 0.995 | 1.107 | 1.226 | 1.349 |

−10 | 0.829 | 0.883 | 0.96 | 1.047 | 0.82 | 0.879 | 0.936 | 1.032 | 0.987 | 1.099 | 1.219 | 1.341 |

−9 | 0.82 | 0.874 | 0.95 | 1.037 | 0.814 | 0.871 | 0.928 | 1.022 | 0.977 | 1.09 | 1.209 | 1.332 |

−8 | 0.81 | 0.862 | 0.937 | 1.024 | 0.807 | 0.862 | 0.917 | 1.01 | 0.965 | 1.078 | 1.198 | 1.32 |

−7 | 0.797 | 0.848 | 0.921 | 1.007 | 0.798 | 0.851 | 0.904 | 0.996 | 0.95 | 1.063 | 1.183 | 1.305 |

−6 | 0.78 | 0.829 | 0.9 | 0.984 | 0.788 | 0.837 | 0.889 | 0.977 | 0.929 | 1.044 | 1.164 | 1.287 |

−5 | 0.758 | 0.805 | 0.872 | 0.955 | 0.775 | 0.819 | 0.869 | 0.953 | 0.902 | 1.018 | 1.139 | 1.263 |

−4 | 0.729 | 0.774 | 0.837 | 0.918 | 0.759 | 0.797 | 0.844 | 0.923 | 0.864 | 0.984 | 1.106 | 1.23 |

−3 | 0.691 | 0.734 | 0.793 | 0.872 | 0.739 | 0.769 | 0.813 | 0.885 | 0.813 | 0.937 | 1.06 | 1.185 |

−2 | 0.645 | 0.687 | 0.743 | 0.816 | 0.712 | 0.736 | 0.776 | 0.84 | 0.748 | 0.878 | 1 | 1.124 |

−1 | 0.595 | 0.635 | 0.688 | 0.756 | 0.678 | 0.697 | 0.736 | 0.791 | 0.677 | 0.811 | 0.93 | 1.045 |

0 | 0.548 | 0.583 | 0.635 | 0.698 | 0.639 | 0.655 | 0.692 | 0.74 | 0.609 | 0.74 | 0.856 | 0.96 |

1 | 0.506 | 0.536 | 0.586 | 0.648 | 0.597 | 0.612 | 0.649 | 0.691 | 0.545 | 0.666 | 0.779 | 0.876 |

2 | 0.469 | 0.497 | 0.546 | 0.61 | 0.557 | 0.572 | 0.608 | 0.646 | 0.476 | 0.586 | 0.7 | 0.793 |

3 | 0.438 | 0.464 | 0.515 | 0.582 | 0.522 | 0.538 | 0.572 | 0.608 | 0.406 | 0.506 | 0.626 | 0.716 |

4 | 0.411 | 0.437 | 0.491 | 0.561 | 0.491 | 0.511 | 0.544 | 0.579 | 0.345 | 0.438 | 0.564 | 0.655 |

5 | 0.387 | 0.415 | 0.472 | 0.547 | 0.465 | 0.489 | 0.522 | 0.558 | 0.298 | 0.388 | 0.517 | 0.609 |

6 | 0.367 | 0.397 | 0.457 | 0.536 | 0.443 | 0.472 | 0.505 | 0.542 | 0.262 | 0.35 | 0.482 | 0.575 |

7 | 0.349 | 0.381 | 0.445 | 0.527 | 0.424 | 0.458 | 0.492 | 0.53 | 0.236 | 0.323 | 0.457 | 0.55 |

8 | 0.333 | 0.368 | 0.445 | 0.521 | 0.409 | 0.447 | 0.481 | 0.521 | 0.216 | 0.302 | 0.437 | 0.53 |

9 | 0.32 | 0.356 | 0.426 | 0.515 | 0.395 | 0.431 | 0.473 | 0.513 | 0.188 | 0.286 | 0.421 | 0.515 |

10 | 0.309 | 0.347 | 0.419 | 0.51 | 0.384 | 0.431 | 0.466 | 0.508 | 0.188 | 0.273 | 0.409 | 0.503 |

11 | 0.299 | 0.338 | 0.413 | 0.506 | 0.375 | 0.424 | 0.46 | 0.503 | 0.178 | 0.263 | 0.399 | 0.493 |

12 | 0.291 | 0.331 | 0.407 | 0.503 | 0.367 | 0.419 | 0.455 | 0.499 | 0.17 | 0.254 | 0.39 | 0.484 |

coefficients to determine, include change scope Δ h ( q ) of the generalized Hurst coefficient of h ( q ) , multiple fractal spectral widths Δ α , non-linear relation curve of scale coefficient τ ( q ) and q, as well as the changing relation curve of multiple fractal strange scores f(α) and the strange degree coefficient α. Concrete judgment standard is: the bigger of Δ h ( q ) and Δ α , the bigger of multiple fractal intensity of corresponding runoff series; and the stronger of non-linearity on the τ ( q ) and q’s changing relation curve , the bigger of corresponding runoff series’ multiple fractal intensity; And, in changing relation curve of multiple fractal strange scores f(α) and the strange degree coefficient α, the bigger scope of the small fluctuation means the bigger multiple fractal intensity of runoff series. Among these, Δ h ( q ) could be calculated out directly, and Δ α , τ ( q ) and f(α) may obtain according to formula 6, and calculated result of Linjiacun, Xianyang and Huaxian station of main channel of Weihe river are shown in

At the main channel of Weihe River, Δ h ( q ) and Δ α of Linjiacun station’s daily runoff series is smaller than Xianyang station obviously (see

From

To Weihe River’s branches, Δ h ( q ) and Δ α of daily runoff series of Yingge station at Shitouhe is bigger than Zhangjiashan station at Jinghe River obviously (

Station | Linjiacun | Xianyang | Huaxian county |
---|---|---|---|

Δh(q) | 1.07 | 1.48 | 0.92 |

Δα | 1.23 | 1.66 | 1.08 |

Stations | Zhangjiashan | Zhuangtou | Yingge |
---|---|---|---|

Δh(q) | 0.81 | 0.87 | 1.35 |

Δα | 0.96 | 1.01 | 1.52 |

multiple fractal intensity of daily runoff series of Shitouhe river is stronger than Jinghe river and Beiluohe river. Meanwhile from the changing relation curve of τ ( q ) along with q in

From changing curve of multiple fractal strange spectrum f(α) along with α(q) on daily runoff series of each Weihe River’s branch station (

curve and all turn towards right. This indicates that the small fluctuation’s influence stands at the first place to Weihe River branches’ daily runoff series, this also is one of the most important reasons of why Weihe River branches’ daily runoff series have long distance relevant respectively. At the same time, in

From the changing feature of strange spectral width Δα (

Size of runoff series’ multiple fractal intensity represented the runoff series’ adulatory property, namely, the stronger of runoff series multiple fractality, the bigger of runoff series’ adulatory property, and the greater of factor’s change which could affect runoff process’s fluctuation. This provided the research basis for studying effect degree of factors’ changing to runoff fluctuation.

Multiple fractal strange scores f(α) along with α(q)’s change value of daily runoff series of each station could be obtained through fitting (

Above analysis shows us that fluctuation degree of Xianyang station’s runoff is the biggest one, the reason lies in that runoff process between Linjiacun station and Xianyang station strongly affected by water income and human activity, thus leading to the very fierce change, on the runoff fluctuation. But Huaxian station’s runoff fluctuation is weaker than Linjiacun and Xianyang station, the reason is Huaxian station situate at downstream of main channel of Weihe River, many factors, such as rainfall, evaporation, water income from tributaries and human activity have an effect of interaction and adjustment each other, therefore cause the result of smaller runoff process fluctuation and tends to steadily.

Similarly, the runoff series’ fluctuation degree of Zhangjiashan station, Jinghe River and Zhuangtou station, Beiluohe river is almost same, but Shitouhe’s runoff fluctuation degree is far bigger than above two stations. The reason lies in these facts, Zhangjiashan station place at Jinghe River, the biggest branch of Weihe river, and Zhuangtou station locate at Beiluohe river, the second big branch of Weihe river, controlled drainage area of these two hydrologic stations is big (occupying proportion of Weihe River drainage area is 33.7% and 20% respectively), factors that can affect their runoff’s formation are many more, so adjustment ability of these two basins is relative strong. But Shitouhe river originates from Qinling mountain area, and its Yingge hydrologic station only has a rather small controlling drainage area, the runoff forms quickly, the mutual function’s counterbalance among effect factors is also weak, so the factor’s change shows a remarkable influence to runoff fluctuation.

q | Lingjiacun | Xianyang | Zhangjiashan | Zhuangtou | Yingge | |||||
---|---|---|---|---|---|---|---|---|---|---|

α | f(α) | α | f(α) | α | f(α) | α | f(α) | α | f(α) | |

−12 | 1.722 | 0.036 | 1.54 | 0 | 1.455 | 0.037 | 1.538 | 0.044 | 2.028 | 0.029 |

−11 | 1.72 | 0.053 | 1.54 | 0 | 1.453 | 0.054 | 1.536 | 0.063 | 2.027 | 0.044 |

−10 | 1.718 | 0.078 | 1.54 | 0.001 | 1.451 | 0.079 | 1.533 | 0.091 | 2.025 | 0.065 |

−9 | 1.714 | 0.113 | 1.54 | 0.002 | 1.447 | 0.114 | 1.529 | 0.13 | 2.022 | 0.097 |

−8 | 1.708 | 0.162 | 1.54 | 0.006 | 1.441 | 0.164 | 1.523 | 0.183 | 2.016 | 0.142 |

−7 | 1.699 | 0.23 | 1.539 | 0.013 | 1.432 | 0.232 | 1.513 | 0.254 | 2.008 | 0.206 |

−6 | 1.685 | 0.319 | 1.536 | 0.03 | 1.418 | 0.322 | 1.499 | 0.347 | 1.994 | 0.292 |

−5 | 1.664 | 0.434 | 1.529 | 0.068 | 1.397 | 0.436 | 1.478 | 0.462 | 1.974 | 0.406 |

-4 | 1.634 | 0.571 | 1.511 | 0.146 | 1.367 | 0.573 | 1.448 | 0.596 | 1.943 | 0.544 |

−3 | 1.591 | 0.72 | 1.467 | 0.298 | 1.324 | 0.721 | 1.407 | 0.739 | 1.898 | 0.699 |

−2 | 1.534 | 0.86 | 1.364 | 0.55 | 1.268 | 0.861 | 1.354 | 0.87 | 1.838 | 0.848 |

−1 | 1.465 | 0.962 | 1.157 | 0.851 | 1.199 | 0.962 | 1.29 | 0.965 | 1.763 | 0.959 |

0 | 1.389 | 1.000 | 0.846 | 1.000 | 1.123 | 1.000 | 1.220 | 1.000 | 1.680 | 1.000 |

1 | 1.312 | 0.962 | 0.536 | 0.851 | 1.048 | 0.962 | 1.15 | 0.965 | 1.596 | 0.959 |

2 | 1.244 | 0.86 | 0.329 | 0.55 | 0.979 | 0.861 | 1.086 | 0.87 | 1.522 | 0.848 |

3 | 1.187 | 0.72 | 0.225 | 0.298 | 0.923 | 0.721 | 1.033 | 0.739 | 1.462 | 0.699 |

4 | 1.144 | 0.571 | 0.181 | 0.146 | 0.88 | 0.573 | 0.992 | 0.596 | 1.417 | 0.544 |

5 | 1.114 | 0.434 | 0.163 | 0.068 | 0.85 | 0.436 | 0.962 | 0.462 | 1.386 | 0.406 |

6 | 1.093 | 0.319 | 0.156 | 0.03 | 0.829 | 0.322 | 0.941 | 0.347 | 1.366 | 0.292 |

7 | 1.079 | 0.23 | 0.154 | 0.013 | 0.815 | 0.232 | 0.927 | 0.254 | 1.352 | 0.206 |

8 | 1.07 | 0.162 | 0.153 | 0.006 | 0.806 | 0.164 | 0.917 | 0.183 | 1.344 | 0.142 |

9 | 1.064 | 0.113 | 0.152 | 0.002 | 0.8 | 0.114 | 0.911 | 0.13 | 1.338 | 0.097 |

10 | 1.06 | 0.078 | 0.152 | 0.001 | 0.796 | 0.079 | 0.907 | 0.091 | 1.335 | 0.065 |

11 | 1.058 | 0.053 | 0.152 | 0 | 0.794 | 0.054 | 0.904 | 0.063 | 1.333 | 0.044 |

12 | 1.056 | 0.036 | 0.152 | 0 | 0.792 | 0.037 | 0.903 | 0.044 | 1.332 | 0.029 |

To 2 - 5 step MF-DFA, generalized Hurst coefficient h ( q ) of main channel’s three stations and the north and south branch stations’ daily runoff process is all not a constant in Weihe River, and all changes along with the change of q, h ( q ) displays as a digression function of q. These disclosed that the daily runoff process of Weihe River’s main channel and braches all have characteristic of multiple fractals.

At the main channel of Weihe River, fluctuation of the daily runoff process of upstream Lingjiacun station is stronger than downstream station, Xianyang station. At the same time, daily runoff series’ fluctuation property of branches in north is stronger than that branch in south obviously.

Runoff series’ multiple fractal intensity of downstream Xianyang station is bigger than upstream Lingjiacun station, it means that Lingjiacun station’s daily runoff series meet more external factors’ influence. The reason lies in that, between Linjiacun and Xianyang station, the runoff process affected seriously by water income and human activity, causes a fierce runoff fluctuation. But, situation goes to another side at downstream Huaxian station, here runoff fluctuation shows a weak fluctuation compared with Linjiacun and Xianyang station, the reason is this station situated at the downstream of Weihe River’s main channel, in a larger catchment area, regional rainfall, more evaporation, more branches’ flow converge and human activity, these factors superimpose and pin down mutually, cause the runoff process to tend to be steadily, fluctuation decreases on the contrary.

For same reason, runoff series’ multiple fractal intensity of Zhangjiashan and Zhuangtou—two tributary stations is greater than that of Yingge station, just because two former stations’ catchment area is bigger, its runoff influence factors are more than Yingge station’s, but still in a limiting score; although we still can’t tell what it is, Yinngge station’s drainage area is small, its runoff effects factors have less opportunity to counterbalance mutually, so the runoff fluctuation is fierce.

The authors thank Zhou Xu for assistance with processing of data. The research was supported by Chinese National Natural Sciences Fund Project (51179160).

The authors declare no conflicts of interest regarding the publication of this paper.

Yan, B.W., Zhou, A.K. and Song, S.B. (2019) Multi-Fractal Characteristics of Daily Runoff Series of Weihe River Watershed. Journal of Water Resource and Protection, 11, 1146-1160. https://doi.org/10.4236/jwarp.2019.119067