Using the Rescaled Range Anal Ysis for the Stu Dy Of
Total Page:16
File Type:pdf, Size:1020Kb
Homage to Ramon Marga/ef; or Why there is such p/e asure in studyin g nature 21 (J. D. Ros & N. Prat, eds.). 1991 . Oec% gia aquatica. 10: 21-33 USING THE RESCALED RANGE ANALYSIS FOR THE STUDY OF HYDROLOGICAL RECORDS: THE RIVER TER AS AN EXAMPLE 1 1 2 JOAN ARMENGOL , SERGI SABATER , ANTONI VIDAL & FRANCESC 1 SABATER Departament d'Ecologia. Facultat de Biologia. Universitat de Barcelona. Av. Diagonal, 4 Barcelona. 1 6 5. 08028 Spain 2 Junta d' Aigües. Abastament d' Aigua del Ter. Cardedeu. 08440 Barcelona. Spain Received: November 1990 SUMMARY Water discharges into the Ter river have a strong influence on the water chemistry and on the dynamics of benthic cornmunities. For this reason, the temporal pattems followed by monthly runoff from 1954 to 1988 have been studied. From a statistical point of view, discharges show a very regular seasonal cycle, but in practice, climatic events introduce a source of variability on runoffs. To determine whether the series of 420 monthly discharges fluctuates with sorne regularity or at random, a rescaled range analilsis was made. The persistence (H) of the series of discharges was measured by the equation R / S = (T / 2) . In fact, H is a measure of the existence of clear trends or periodicities in the records of persistent stocastic processes. An additional measure of the persistence was obtained by means of the relationship D = 2 - H, where D is the local fractal dimensiono For comparative reasons, a random series of the same length, mean and standard deviation was generated and analysed in the same way as the Ter discharges were. The results were H = 0.68 and D = 1 .38 and H = 0.44 and D = 1.52, for the Ter data and for the random series respectively. These results show that the simulated series is a case of ordinary Brownian motion while the Ter discharges series has an intermediate value of persistence. Because of the value of H in the Ter series, there are periods whose values were estimated by the autocorrelation function and the periodogram of the series. The results show the existence of short fluctuations of 3, 6 and 12 months deterrnining the seasonality of the annual cycle, and large fluctuations with periods of 5.5, 8.6, 10.1 and 11.7 months which can be considered an expression of cycles of circa 7 and 1 1 years. KEYWORD S: River discharge, fractal dimension, Hurst exponent, fractional Brownian motion, persistence . In the case of the landscape, rain fa lling at random and moving from each point of the surface along the steepest profi le, through eros ion and convergence of runoffs, generates a secondary "hydrologic" relief Will the fractal quality be preserved (MARGALEF, ? 1988) INTRODUCTION seasonal c1imatic changes. Nevertheless, this pattern is usuall y modified by irregular In Mediterranen ecosystems , annual climatic fluctuations that cause sudden cyc1es are always described in terms of changes in their dynamics. Temporal ponds regular fluctuations associated with and rivers are among the most affected 22 ARMENGOL, SABATER, VI DAL & SABATER ecosystems because they are mainly taking place (SCHINDLER, 1988; controlled by the periodicity and intensity ARMENGOL, 1988). of the rains. The tumover time of the water In rivers, strong rains increase the acts as a key factor in most of the processes intensity of the advective transport of solid 250 A 200 (") 150 S :::c a 100 50 3 2 -2 -3 1955 1960 1965 1970 1975 1980 1985 3 I FIGURE 1. River Ter monthly discharges for the 1954 - 1988 periodo A) original series of data in Hm month- . B) normalized series by Y = . C) the standardized sequence. XO.22 RESCALED RANGE ANALYS IS OF THE TER RIVER 23 and dissolved substances. Because of this, (backward steps) (MARGALEF, 1977). the spates act as huge inputs of external In the Ter river, water flow is the most energy that simplify the structure of the important factor determining the annual river. Most of the organisms are carried variability in the environmental conditions down and the chemical composition of (SABATER & ARMENGOL, 1986) and in water beco mes more similar along the ri ver. the dynamics of phytobenthic (SABATER The annual cyc1e of the cornmunities can et al., 1988) and macroinvertebrate (PUIG be seen as a sequence of populations with et al., 1987) communities. These relation short life span, well adapted to seasonal ships are well established for short periods changes but able to reconstruct the of time, such as the annual cyc1e, but for cornmunity after each event (MARGALEF, longer periods it is assumed that c1imatic 1980). If longer periods of time are fluctuations have a similar effect. If we considered, changes in the communities can accept these relationships for the Ter river, be seen as a sequence of periods of the long-term study of a record of organisation (secondary successions) hydrological events as large as possible can followed by quick disruptive episodes provide us with sorne insight into how it A 3 2 3 B 2 o -1 -2 -3 1955 1960 1965 1970 1975 1980 1985 FIGURE 2. Simulated sequence of 420 values (35 years x 12 months) nonnally distributed with the same mean and standard deviation as the data of figure lB . B) the standardized series. 24 ARMENGOL, SABATER, VIDAL & SABATER functions from an ecological point of view. (MANDELB ROT & WALLIS, 1969); for For longer periods it is well established that this reason the existence of temporal the water discharges affect the chemical "cycles" has been studied here by mean s of pollution and the eutrophication levels in autocorrelation and spectral analysis. the reservoirs situated in the middle stretch of the river (VIDAL, 1977, 1987, 1989; ARMENGOL et al., 1986). MATERIAL AND METHODS Long-term records of hydrological events show erratic behaviour at different The record of hydrological data used in time scales (FEDER, 1988). These records this study is composed by the monthly can be studied by means of the rescaled discharges at Roda de Ter station from range analysis (HURST et al., 1965) and 1954 to 1988 (Fig. l). The data were characterized by an exponent H, called the obtained from the Dirección General de Hurst exponent, which is a measure of the Obras Hidráulicas gauging yearbooks persistence or "memory" of one series (D.G.O.H., 1972, 1973, 1974 and 1975; (MANDELBROT & VAN NESS, 1968). D.G.O.H., 1976) for the period 1965- 1977, MANDELBROT (1985) has shown that and were provided by the Servei Hidrologic this record is a self-affine fractal curve with de la Generalitat de Cata]unya for the a dimension D = 2 - H. For a practical periods 1954- 1 965 and 1978-1 988. perspective, the H exponent measures how Roda de Ter is situated in the middle far a temporal series of data departs from a stretch of the Ter river, 88.6 Km from the random series. source (river length 200,3 Km), where it 2 We have used a record of 35 years of has a drainage area of 1,476.8 Km (48. 11 monthly average discharges to analyse the % of the total). At this point the Ter is a persistence in the data sequence. The aim fourth order river. This station was selected of this study is to determine the temporal because it is situated aboye the reservoirs, trend followed by the series and how this thus excluding their regulatory effect, and key factor carried weight on the river at the same time the drainage area is large considered as an ecosystem. It is well enough to avoid the effect of sudden known that wet and dry years appear in discharge due to a local rainfall. clusters, the "Joseph effect" The hydrological characteristics of the 100 B A 75 75 ..c::..., s:: o 8 50 "'.......... 8 :r: 25 25 JF MAMJ JASo ND FIG URE Average annual cycle showing the mean monthly discharges with confidence interval: A) rol' 3. 95 % Ter discharges (Fig. lB) and B) for random series (Fig. 2 A). RESCALED RANGE ANALYSIS OF THE TER RIVER 25 3 Ter are typical of the Mediterranean region, for the period studied was 4l.3 3.74 Hm ±3 but the proximity of the Pyrenees also has with ranges between 23 l.4 Hm and 0.70 3 an influence. The annual cyc1e is Hm , while, if the data are grouped by 3 characterized by two periods of maximum years, a mean of 554.6 35.7 Hm is ± 3 water flow; the first, in May (28.6 6.0 obtained with a maximum of 1180 Hm in 3 1 ± 3 Hm month- , p <0.05 for all the data of 1972 and a rninimum of 258 Hm in the this paper), is the highest due to the effect year 1957. of the spring rains plus the thaw of the The Ter hydrological record does not snow accumulated in the Pyrenees, while follow a Gaussian distribution (the kurtosis the second, in November is less important is 4.9 12 and the skewness is 2.001). For ¡ 3 - (17.6 ± 5.7 Hm month ) and exc1usively this reason, it is not possible to use many of due to the autumn rains. the parametrical statistical methods. In Two minima are situated in between, the order to transform the original data into a 3 I lower in August (1 1.4 2.6 Hm month- new series with a near normality ± d and the other in February (14.2 3.9 Hm distribution, a power transformation has 1 ± month- ). The average monthly water flow been sought by means of the Box & Cox A ..c:: s:= 100 O E ...........