Mathematical and Numerical Modeling of Flow and Transport 2013View this Special Issue
Research Article | Open Access
Ni Zhihui, Wu Lichun, Wang Ming-hui, Yi Jing, Zeng Qiang, "The Fractal Dimension of River Length Based on the Observed Data", Journal of Applied Mathematics, vol. 2013, Article ID 327297, 9 pages, 2013. https://doi.org/10.1155/2013/327297
The Fractal Dimension of River Length Based on the Observed Data
Although the phenomenon that strictly meets the constant dimension fractal form in the nature does not exist, fractal theory provides a new way and means for the study of complex natural phenomena. Therefore, we use some variable dimension fractal analysis methods to study river flow discharge. On the basis of the flood flow corresponding to the waterline length, the river of the overall and partial dimensions are calculated and the relationships between the overall and partial dimensions are discussed. The law of the length in section of Chongqing city of Yangtze River is calibrated by using of variable fractal dimension. The results conclude that it does express a second-order accumulated variable-dimensional fractal phenomenon, and the dimension can reflect the degree of the river; the greater dimension, the more the river bend. It has different dimensions at a different location in the same river. In the same river, the larger dimension, the worse flow discharge capacity of the river and the more obvious of the flood will be on the performance.
A fractal dimension is a ratio providing a statistical index of complexity comparing how detail in a pattern changes with the scale at which it is measured. It has also been characterized as a measure of the space-filling capacity of a pattern that tells how a fractal scales differently than the space it is embedded in; a fractal dimension does not have to be an integer.
Fractals have been introduced in order to quantify the self-similarity observed in nature while at the same time to make the study of nondifferentiable processes possible. Given that this self-similar behavior has often a “local” character, the theory of fractals was generalized to multifractals, enabling the description of more complex phenomena with varying fractal properties. Examples of processes that have been thus treated are the energy dissipation in turbulence and the price increments in finance. In the field of geophysics and atmospheric physics, fractal and multifractal analyses have been extensively applied , since self-similarity is present in a wide variety of such phenomena, from the distribution of earthquake epicenters [2–4] and hypocenters  to climate change  and atmospheric turbulence [7, 8]. In the same field of research, fractal and multifractal methods have been used both to characterize the long-term behavior of related signals and to indicate possible precursors in experimental time series of their records, yielding very promising results [9–12].
Fractal processes have become one of the most widely used modeling tools in science and engineering, with diverse applications in finance, physics, network traffic, and recently geography sciences. In the study of geometric properties of dynamical systems or fractal measures, one is often interested in the asymptotic behaviour of local quantities associated with the underlying dynamical or geometric structure. For example, one is often interested in the ergodic average of a continuous function, the local entropy or the local Lyapunov exponent, or the local dimension of a measure. These quantities provide a description of various aspects of measures or dynamical systems, for example, chaoticity, sensitive dependence, and so forth. All these quantities provide important information about the underlying geometric or dynamical structure . The mixing length based on fractal theory has been calculated and analyzed . Jou et al.  by assuming a self-similar structure for the Kelvin waves along vortex loops with successive smaller scale features model the fractal dimension of a super fluid vortex tangle in the zero temperature limits. Their model assumes that at each step the total energy of the vortices is conserved but the total length can change. They obtain a relation between the fractal dimension and the exponent describing how the vortex energy per unit length changes with the length scale. In addition, many scholars have also concerned about relationship between the fractal theory and river systems [16–22].
The analysis of river flows has a long history; nevertheless some important issues have been lost. Many scholars use some fractal analysis methods to study river flow fluctuations. Sadegh Movahed and Hermanis  have studied one component of the climate system, the river flux, by using the novel approach in the fractal analysis like detrended fluctuation analysis, fourier-detrended fluctuation analysis and scaled windowed variance analysis methods. The statistical and fractal analysis of river flows should be an important issue in the geophysics and hydrological systems to recognize the influence of environmental conditions and to detect the relative effects. A set of most important results which can be given by using statistical tools are as follows: a concept of scale self-similarity for the topography of Earth’s surface , the hydraulic-geometric similarity of river system and floods forced by the heavy rain , and so forth. Already more than half a century ago the engineer Hurst found that runoff records from various rivers exhibit “long-range statistical dependencies.” Later, such long-term correlated fluctuation behavior has also been reported for many other geophysical records including precipitation data. These original approaches exclusively focused on the absolute values or the variances of the full distribution of the fluctuations, which can be regarded as the first and second moments of detrended fluctuation analysis [24, 26]. In the last decade it has been realized that a multifractal description is required for a full characterization of the runoff records . This multifractal description of the records can be regarded as a “fingerprint” for each station or river, which, among other things, can serve as an efficient nontrivial test bed for the state-of-the-art precipitation-runoff models.
In this work, the fractal dimension of the river (hereinafter referred to as dimension) is generated by studying characteristics of fractal structure. The dimension of the river is divided into river length and river network, to explore the fractal dimension from the view of the entire river length, to be called the unitary dimensions of river length. However, it has different dimensions at a different location in the same river. Meandering is different, so the dimension of rivers length (hereinafter referred to as part dimension) is, even vary considerably.
Chongqing is an oversize industrial city and the water-land transport hub that developed relying on the Yangtze River and Jialing River. The main section of Chongqing city is from Daduko to Tongluoxia; tributary section is from Jingkou of Jialing River to Chaotianmen. The total length of it is about 60 km. The section of Chongqing city of Yangtze River is located in the southeastern edge of Sichuan basin; the section of river is the lotus root shape. Because the river is affected by geological structure and its lithology changes, the longitudinal profile along the higher ups and downs changes greatly, and the river boundary conditions are very complicated as the section of Yangtze River is located in fluctuating backwater area of Three Gorges Project. Based on statistic water surface profile when flow discharge is 58000 m3/s of Cuntan station, the average river width of Jiulongpo and Caiyuanba is more than 900 m, the total length of the two sections is about 6 km, and the length of narrow reach which is below 600 m is 1.7 km. The reach in the section of Chongqing city performance is of continuous and irregular curve shape, which has 6 continuous bends. These curves are slowly bending 150° and 90° elbow. Chongqing section of the Yangtze River in 2007 is shown in Figure 1.
At present, the research on the fractal characteristics of river length in section of Chongqing city of Yangtze River is few, and it is still in the theoretical exploration and analysis phase. Chongqing is located in the upper reaches of the Yangtze River, is an important strategic position, and is the largest port city in the upper reaches of the Yangtze River. Therefore, measurement on fractal dimension of river length in section of Chongqing city of Yangtze River can give full play to the effect of golden waterway in Yangtze River. It will play a leading role in water transport of Chongqing, and provide better services for the upper reaches of Yangtze River and the western region. Thus, the study of the fractal dimension of river length in section of Chongqing city of Yangtze River is of great significance. In this paper, in the view of river length, the fractal dimension and its relationship with floods will to be studied.
2.1. Building of 1D Model Program
The former USSR, North America, Western Europe, and China have carried out research and application of hydrodynamic, sediment transport model since the 1950s. A one-dimensional mathematical model of water flow has been more mature after several decades’ development and application.
This research is studied through one-dimensional flow mathematical model. In view of the constant flow calculation which has been more mature, a one-dimensional constant flow mathematical model is taken to the whole river.
2.1.1. Basic Equations
Here is a 1D water movement and continuity equation: where is cross-section area, is flow discharge, is section average flow velocity, is water level, is river width, is horizontal direction, is time, is water density, and is shear stress of river bottom.
Simplified equations is used in the calculation process, and water movement is changed to where is downstream section level, is upstream section level, is downstream section average velocity, is upstream section average velocity, is downstream section kinetic energy correction factor, is upstream section kinetic energy correction factor, is frictional head loss, and is local head loss.
2.1.2. Boundary Condition
Control conditions of this model are the flow discharge of upstream and water level of downstream. Stage-discharge data of Cuntan station is the boundary condition of downstream. According to the hydrological data in 2011 and 2012 of Cuntan station, stage-discharge relation of Cuntan station is shown in Figure 2 and Table 1.
2.1.3. Model Verification
Roughness is vital to the calculation of water surface profile. Roughness must adjust and revise repeatedly until to the deviation between calculated value and measured value. Finally, the ideal roughness is 0.036.
There is water level date from May 1, 2009, to October 17 of the section of Chongqing city. When water level of downstream is 161.66 m, flow discharge of Cuntan station is 6320 m3/s. The water level of calculation and observation is shown in Table 2.
From Table 2 the following can be seen. The calculation water level of the section of Chongqing city is similar to the observe water level. The result of calculation is reasonable.
In this work, 1D model program mentioned above is used to calculate the length of waterline () corresponding to 25 class flow discharge ().
2.2. Building of Fractal Dimension
According to the definition of Mandelbrot, fractal refers to the body that the part is similar to the whole in some way. Mandelbrot (1967) put forward the formula of a statistical fractal dimension estimated for the self-similar fractal case: where is Euclidean length (waterline), is the measured size (flow discharge), is the fractal dimension, and is a proportion constant.
Take the natural logarithm to (3), then get it as follows:
And then paint and on the coordinates of -axis and -axis, respectively, and last, use the least square method to fit the straight line and its slope is . We can derive its fractal dimension.
If the fractal dimension is calculated as a constant, it is simple fractal dimension, and if not, it needs to be described as variable fractal dimension . In fact, the phenomenon of strictly meet the simple fractal dimension form does not exist in nature, a large number of complex phenomena need to use variable fractal dimension to describe.
In this thesis, fractal theory of cumulative sum sequence is used for calculation.
This specific method steps are as follows.(1)Determine the raw data , where orders from small to large, . There are some data points in log-log coordinates. And then calculate the slope of two adjacent points by the use of (5); the sub variable fractal dimension is . In general, the fractal dimensions change much and there is no law: (2)Construct the accumulated sum of a total order, and is the basic sequence. Then it is constructed by following the rules: where are the accumulated sum of first order, second order, and third order, (3)Establish variable fractal dimension model of the accumulated sum of a total order, taking the first order as an example, and the variable fractal dimension is the opposite slope of data points calculated by (6) in the log-log coordinates.According to the data of , fractal dimensions of are obtained, known as the fractal dimension sequence. is variable fractal dimension sequence of the accumulated sum of a total order, .(4)Determine the better order of the accumulated sum, and identify the corresponding fractal dimension.
3. Data and Results
Based on the field observations of topographic map in the Chongqing section of the Yangtze River in 2007 (Figure 1), a valuable data and analysis results by (4) are given. In this paper, the results are shown in Figure 3.
(a) The left bank (median flood)
(b) The left bank (flood)
(c) The right bank (median flood)
(d) The right bank (flood)
From the double logarithmic coordinates it can be seen that (Figure 3) the data point is clearly not a straight line. From the analysis of it, the relationship that the river length of Chongqing section of the Yangtze River was able to meet the variable fractal dimension should be applicable to subdimensional variable fractal model. Therefore, we have adopted variable fractal model to calculate the fractal dimension. The calculation result of the fractal dimension during median flood period on the left bank D2 is −0.3321, the correlation coefficient is 0.9855; the fractal dimension on the right bank D2 is −0.3323, the correlation coefficient is 0.9854. The calculation result of the fractal dimension during flood period on the left bank D2 is −1.6844, and the correlation coefficient is 0.9834; the fractal dimension on the right bank D2 is −1.6859, and the correlation coefficient is 0.9834 (Tables 3 and 4).
From Figure 4 it can be seen that, after the transformation of second-order accumulated and variable dimensional fractal, the data points fit better to a straight line, which means that the river length in section of Chongqing city of Yangtze River has the characteristics of second-order accumulated variable dimensional fractal. Thus, the river length in the section of Chongqing city of Yangtze River has characteristics of second-order fractal dimension.
(a) The left bank (median flood)
(b) The left bank (flood)
(c) The right bank (median flood)
(d) The right bank (flood)
4. Analysis and Discussion
The fractal dimension of the river length reflects the degree of bending of the river. It is shown that the greater the fractal dimension of the river length, the more tortuous of the river. On the contrary, the river is straighter. It has different dimensions of the river length and the degree of bending at a different location in the same river. In terms of the flood, the possibility and intensity of flooding in a different reach are different. Thus, calculating the fractal dimension values of the whole river length has little significance. The river segmentation being carried out, which calculated the value of the fractal dimension in different sections, found out the correlation between the fractal dimension of each reach and flood.
From the qualitative analysis of the possibility of flood and its fractal dimension on various river reaches, it is shown that there are some relationships between them. The greater fractal dimensions of the river length, the more tortuous of the river. The worse flood carrying capacity of the rivers, the more obvious the flood will be on the performance. However, from quantitative analysis, what kind of relationship exists between the fractal dimension of river length and the flood? Based on Chongqing city segments of Yangtze River, the relationship between them is explored by using quantitative calculation and measured data.
4.1. Calculation of Local Fractal Dimension
The section of Chongqing city of Yangtze River is divided into six sections. So we can calculate the fractal dimension of each reach and the correlation coefficient by the left and right sides separately, and the concrete results were shown in Tables 5 and 6.
Where , .|
Where , .|
From Tables 5 and 6 it can be seen that the correlation coefficient on the fractal dimension of the length for each reach in section of Chongqing city of the Yangtze River is more than 98%. It is shown that the length has good characteristics of fractal dimension, and it can reflect the characteristics of it.
In Table 5, the fractal dimensions during flood period on the left bank of river length for the three reaches of Lijiatuo Bridge to Egongyan Bridge, Egongyan Bridge to Caiyuanba Bridge, and Caiyuanba Bridge to Shibanpo Bridge are −1.7203, −1.6844, and −1.6776, respectively. In Table 6, the fractal dimensions on the right bank of river length are −1.7071, −1.6805, and −1.6736, respectively. And the dimension can reflect the degree of bending of the river; the greater dimension, the more tortuous of the river. So we can get that the bending degree of the four reaches of Caiyuanba Bridge to Shibanpo Bridge, Egongyan Bridge to Caiyuanba Bridge, and Lijiatuo Bridge to Egongyan Bridge is more and more big.
4.2. The Local Fractal Dimension and Overall Fractal Dimension
The arithmetic average of river length does not mean the overall fractal dimension values (see Table 7). In order to further reveal the existence of the law, the situation of the level with the left bank in section of Chongqing city of the Yangtze River is calculated. The fractal dimensions with the upper reaches of the Yangtze River to Lijiatuo bridge and Lijiatuo bridge to Egongyan bridge are−0.3318, −0.3327, respectively. From the upper reaches of the Yangtse River to Egongyan bridge, the arithmetic mean is −0.3322, and it is not equal to fractal dimension −0.3325 (Table 5). The results are also consistent with other reaches of the river and the right bank. So, it can be found that the fractal dimension of length in section of Chongqing city of the Yangtze River is not equal to its part fractal dimension of the arithmetic mean.
4.3. The Relationship between Fractal Dimension and the Flow Discharge
In general, from Tables 5 and 6 and Figure 5, the dimension during median flood period is smaller than the dimension during flood period in the same observation station. And the absolute value of dimension during median flood period is inversely proportional to the flow discharge.
(a) Median flood
Some theory can be derived from the relationship between stage and discharge at the customs of Chongqing Cuntan hydrological station (Table 1); the larger flow, the higher level in the same observation station. Conversely, the greater dimension, the higher flow and the more obvious the flood will be on the performance. This just confirms the relationship between the possibility of flood and fractal dimension: the greater dimension, the more the river bend, and the larger dimension, the worse flow discharge capacity of the river and the more obvious the flood will be on the performance.
In this paper, taking the measured data in section of Chongqing city of Yangtze River as an example explored the fractal characteristics from the perspective of fractal scale. Through analysis, comparison, and discussion in this paper, it draws the following conclusions.(1)The phenomenon of variable dimension fractal, with second-order fractal dimension, exists on the main reaches in section of Chongqing city of the Yangtze River. The fractal dimension value during median flood period of the left bank is −0.3321, the right bank is −0.3323, the fractal dimension value during flood period of the left bank is −1.6884, and the right bank is −1.6859.(2)The dimension can reflect the degree of bending of the river; the greater dimension, the more tortuous of the river. It can be got that the bending degree of the three reaches of Caiyuanba Bridge to Shibanpo Bridge, Egongyan Bridge to Caiyuanba Bridge, and Lijiatuo Bridge to Egongyan Bridge is more and more big.(3)The fractal dimension of length in section of Chongqing city of the Yangtze River is not equal to its part fractal dimension of the arithmetic mean.(4)In the same river, the larger dimension, the more obvious the flood will be on the performance. Therefore, the fractal dimension of the river can be used as a quantitative indicator of flood forecasting. The larger fractal dimensions, the worse capacity of flood carrying. However, due to the impact of floods produced by many factors, such as water level and sediment, the fractal dimension of the river can only be one of the indicators as forecasting floods. Considered, we should identify more predictors of faster, more accurate prediction of flood.
The study was supported by the Natural Science Foundation Project of CQ CSTC (Grant no. cstc2012jjA30002), the Science and Technology Project of Chongqing Education Committee (Grants no. KJ110409, no. KJ111501), the National Engineering Research Center for Inland Waterway Regulation Program (Grant no. SLK2012A02), and the National Key Technology R&D Program (Grant no. 2012BAB05B03).
- S. Lovejoy and D. Schertzer, “Scale invariance, symmetries, fractals and stochastic simulations of atmospheric phenomena,” Bulletin of the American Meteorological Society, vol. 67, no. 1, pp. 21–32, 1986.
- Y. Y. Kagan and L. Knopoff, “Stochastic synthesis of earthquake catalogs,” Journal of Geophysical Research, vol. 86, no. B4, pp. 2853–2862, 1981.
- D. Marsan, C. J. Bean, S. Steacy, and J. McCloskey, “Spatio-temporal analysis of stress diffusion in a mining-induced seismicity system,” Geophysical Research Letters, vol. 26, no. 24, pp. 3697–3700, 1999.
- D. Kiyashchenko, N. Smirnova, V. Troyan, and F. Vallianatos, “Dynamics of multifractal and correlation characteristics of the spatio-temporal distribution of regional seismicity before the strong earthquakes,” Natural Hazards and Earth System Science, vol. 3, no. 3-4, pp. 285–298, 2003.
- M. C. Robertson, C. G. Sammis, M. Sahimi, and A. J. Martin, “Fractal analysis of three-dimensional spatial distributions of earthquakes with a percolation interpretation,” Journal of Geophysical Research, vol. 100, no. B1, pp. 609–620, 1995.
- Y. Ashkenazy, D. R. Baker, H. Gildor, and S. Havlin, “Nonlinearity and multifractality of climate change in the past 420,000 years,” Geophysical Research Letters, vol. 30, no. 22, 2003.
- D. Schertzer, S. Lovejoy, F. Schmitt, Y. Chigirinskaya, and D. Marsan, “Multifractal cascade dynamics and turbulent intermittency,” Fractals, vol. 5, no. 3, pp. 427–471, 1997.
- P. Hubert, Y. Tessier, S. Lovejoy et al., “Multifractals and extreme rainfall events,” Geophysical Research Letters, vol. 20, no. 10, pp. 931–934, 1993.
- P. A. Varotsos, N. V. Sarlis, E. S. Skordas, and M. S. Lazaridou, “Fluctuations, under time reversal, of the natural time and the entropy distinguish similar looking electric signals of different dynamics,” Journal of Applied Physics, vol. 103, no. 1, Article ID 014906, 2008.
- N. Scafetta and B. J. West, “Phenomenological reconstructions of the solar signature in the Northern Hemisphere surface temperature records since 1600,” Journal of Geophysical Research, vol. 112, no. D24, 2007.
- Y. Ida, M. Hayakawa, A. Adalev, and K. Gotoh, “Multifractal analysis for the ULF geomagnetic data during the 1993 Guam earthquake,” Nonlinear Processes in Geophysics, vol. 12, no. 2, pp. 157–162, 2005.
- P. A. Varotsos, N. V. Sarlis, E. S. Skordas, H. K. Tanaka, and M. S. Lazaridou, “Attempt to distinguish long-range temporal correlations from the statistics of the increments by natural time analysis,” Physical Review E, vol. 74, no. 2, Article ID 021123, 2006.
- L. Olsen, “Multifractal analysis of divergence points of deformed measure theoretical Birkhoff averages. IV. Divergence points and packing dimension,” Bulletin des Sciences Mathématiques, vol. 132, no. 8, pp. 650–678, 2008.
- Z. H. Ni, X. J. Zhang, and R. S. Xu, “Fractal study on the vertical concentration distribution of sediment flow in Yangtze River and Yellow River,” Yangtze River, vol. 42, no. 19, pp. 73–76, 2011 (Chinese).
- D. Jou, M. S. Mongiov, M. Sciacca, and C. F. Barenghi, “Vortex length, vortex energy and fractal dimension of superfluid turbulence at very low temperature,” Journal of Physics A, vol. 43, no. 20, Article ID 205501, pp. 1–10, 2010.
- V. I. Nikora, “Fractal structures of river plan forms,” Water Resources Research, vol. 27, no. 6, pp. 1327–1333, 1991.
- P. La Barbera and R. Rosso, “On the fractal dimension of stream networks,” Water Resources Research, vol. 25, no. 4, pp. 735–741, 1989.
- V. K. Gupta and E. Waymire, “Statistical self-similarity in river networks parameterized by elevation,” Water Resources Research, vol. 25, no. 3, pp. 463–476, 1989.
- A. Robert and A. G. Roy, “On the fractal interpretation of the mainstream length-drainage area relationship,” Water Resources Research, vol. 26, no. 5, pp. 839–842, 1990.
- D. G. Tarboton, R. L. Bras, and I. Rodriguez-Iturbe, “A physical basis for drainage density,” Geomorphology, vol. 5, no. 1-2, pp. 59–76, 1992.
- Z. H. Ni, Z. Y. Song, X. J. Zhang, L. C. Wu, and J. Yi, “A modification to vertical distribution of tidal flow Reynolds stress in shallow sea,” China Ocean Engineering, vol. 26, no. 3, pp. 431–442, 2012.
- W. Kinsner, “A unified approach to fractal dimensions,” International Journal of Cognitive Informatics and Natural Intelligence, vol. 1, no. 4, pp. 26–46, 2007.
- M. Sadegh Movahed and E. Hermanis, “Fractal analysis of river flow fluctuations,” Physica A, vol. 387, no. 4, pp. 915–932, 2008.
- B. B. Mandelbrot and J. R. Wallis, “Some long-run properties of geophysical records,” Water Resources Research, vol. 5, no. 2, pp. 321–340, 1969.
- P. Burlando and R. Rosso, “Scaling and multiscaling models of depth-duration-frequency curves for storm precipitation,” Journal of Hydrology, vol. 187, no. 1-2, pp. 45–64, 1996.
- C. Matsoukas, S. Islam, and I. Rodriguez-Iturbe, “Detrended fluctuation analysis of rainfall and streamflow time series,” Journal of Geophysical Research, vol. 105, no. D23, pp. 29165–29172, 2000.
- G. Pandey, S. Lovejoy, and D. Schertzer, “Multifractal analysis of daily river flows including extremes for basins of five to two million square kilometres, one day to 75 years,” Journal of Hydrology, vol. 208, no. 1-2, pp. 62–81, 1998.
- Y. H. Fu, “Fractal dimension and fractals in ocean engineering,” China Ocean Engineering, vol. 8, no. 3, pp. 285–292, 1994.
Copyright © 2013 Ni Zhihui et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.