A methodological approach to water quality ... - Springer Link

6 downloads 979 Views 301KB Size Report
Apr 15, 2008 - concentration was lower than during maximum flow and organic carbon ..... Oxford. Gordon, B. N. D., McMahon, T. A., & Finlayson, B. L. (1994).
GeoJournal (2007) 70:281–288 DOI 10.1007/s10708-008-9134-z

A methodological approach to water quality assessment in an ungauged basin, Buenos Aires, Argentina Silvana Arreghini Æ Laura de Cabo Æ Rafael Seoane Æ Nicola´s Tomazin Æ Roberto Serafini Æ Alicia Fabrizio de Iorio

Published online: 15 April 2008 Ó Springer Science+Business Media B.V. 2008

Abstract The Reconquista River is one of the most polluted rivers in Latin America. This paper aims at identifying the dynamics of water quality in an area with low or ‘‘background’’ concentrations of pollutants within the Reconquista River system in order to better define levels of pollution in the main system. In S. Arreghini (&)  R. Serafini  A. F. de Iorio Departamento de Recursos Naturales y Ambiente, Facultad de Agronomı´a, Universidad de Buenos Aires, Av. San Martı´n, 4453, 1417 Buenos Aires, Argentina e-mail: [email protected] A. F. de Iorio e-mail: [email protected] S. Arreghini  L. de Cabo  R. Serafini Museo Argentino de Ciencias Naturales ‘‘Bernardino Rivadavia’’, Av. Angel Gallardo 470, 1405 Buenos Aires, Argentina L. de Cabo e-mail: [email protected] R. Seoane  N. Tomazin Instituto Nacional del Agua, AU Ezeiza-Can˜uelas, Tramo J. Newbery km 1.620, 1804 Buenos Aires, Argentina R. Seoane Departamento de Hidra´ulica, Facultad de Ingenierı´a, Universidad de Buenos Aires, Av. Las Heras 2214, 1127 Buenos Aires, Argentina R. Seoane National Research Council (CONICET), Buenos Aires, Argentina

order to describe the dynamics of water quality in the background area, we propose a methodology based on flow estimation with the instantaneous unit hydrograph model and on measurements of physical and chemical water variables under different hydrological conditions. Because of high dissolved oxygen and low ammonium and o-phosphate concentrations, the Arroyo Durazno, a tributary stream of the Reconquista River, is defined as a background area. When a storm event begins, the concentration of nitrates and the electrical conductivity diminish. An increase in dissolved organic carbon suggests an important input of carbon from hillslope runoff. The proportion of fulvic and humic acids also increases. On the receding limb of the hydrograph, nitrate concentration was lower than during maximum flow and organic carbon concentration remained high. This behavior, known as the ‘‘flushing effect’’, suggests that the soluble material accumulated in the drainage area during dry periods is transported to the stream by leaching or ‘‘lixiviation’’ and surface runoff, thus raising solute concentrations during the first few hours of the storm. Water quality changes rapidly, even in background areas, due to its dependence on the flow. The methodology followed in this paper can also be applied to other basins with similar characteristics. Due to the difficulty in defining baseline areas for surface waters, a knowledge of background water quality and its dynamics is essential for understanding pollution trends and anthropogenic impacts on rivers.

123

282

Keywords Background area  Background water quality  Unit hydrograph  Temperate river  Hydrological model  Ecohydrology  Pollution Introduction Hydrologic cycle components, land use, soil properties, topography, and the activity of organisms in the drainage basin have an impact on the chemical composition of a water body. In recent years, attempts have been made to integrate knowledge of hydrological and biological processes at basin scale (Zalewski 2002). Identification and quantification of natural and anthropogenic impacts on the chemical composition of river water should be an important part of land and water management in a basin (Petts and Calow 1996). Water quality is a major issue in sustainable development, especially in humid areas of Latin America, where water problems relate more to quality preservation than to shortages. In Argentina, direct estimation of rainfall/flow data is difficult due to a lack of sufficient data from the basins. Therefore, the hydrologic model parameters needed to determine the flows that will be used in chemical data interpretation must be estimated indirectly with, for example, the instantaneous unit Fig. 1 Map of the Reconquista river basin

123

GeoJournal (2007) 70:281–288

hydrograph theory (IUH) (Clark 1945). A basin’s IUH makes it possible to obtain direct flow estimates with a certain level of associated uncertainty, such as the use of a constant value of the Curve Number obtained from the tables as a function of the basin’s soil type and vegetation cover (Soil Conservation Service 1972). The Reconquista River is a lowland river of the Pampa Ondulada in the great Plata Basin. Its main streambed is about 55 km long and drains an area of some 1670 km2 (Fig. 1). Forty-three percent of the land is used as natural pastureland for cattle grazing and for agriculture. The other 57% is urbanized (http://www.indec.mecon.ar). Almost 8% of the country’s population, more than 3 million people, lives in the basin. It comprises a critical area with respect to the conservation status of fresh water ecoregions in Latin America (Olson et al. 1998). Though many studies have assessed the water quality of this river (Castan˜e´ et al. 1998, de Cabo et al. 2000), there is no detailed information on the water quality of a background area, considering as background the natural condition of a system that may include some anthropic components and that is used as the starting point for monitoring activities (Edmunds et al. 2003). This paper seeks to identify a background area for a temperate lowland basin with no regular

GeoJournal (2007) 70:281–288

hydrometeorological records, the Reconquista River, and to describe the dynamics of water quality in that background area using a method of hydrograph estimation.

283

Statistical analysis The Kruskal–Wallis test was used to compare all the data sets, and two-sample comparisons were made using Kolmogorov–Smirnov test when deemed appropriate. A 5% level of significance was used.

Materials and methods Flow estimation Ten water samplings were undertaken in the Reconquista River basin between April 1996 and March 1998 at five sites along the main streambed—CAS, REY, GOR, SMT, BAN (Fig. 1)––and at two sites in the tributaries––DUR and CHO. Downstream from the Roggero Reservoir, where the tributaries meet, industrial and urbanization activities have inexorably led to an inflow of pollutants from point and non-point sources. In addition, 19 water samples were taken at DUR under different hydrological conditions, including a heavy storm (March 2001). Physical–chemical analysis The following measurements were taken in situ: dissolved oxygen (DO) using a YSI 51B oxymeter; pH using an Orion pH meter 250 A; electrical conductivity (EC) with a Luftman conductimeter; temperature; transparency (with a Secchi disk); and depth of the streambed. Triplicate surface water samples were filtered through Whatman GF/C filters, and sent to the laboratory at 4°C for analysis. Concentrations were determined using the following methods: soluble reactive phosphorus (SRP), nitrates (N–NO3 ), and nitrites (N–NO2 ) according to Strickland and Parsons (1972); ammonium (N–NH+4 ) and bicarbonates (HCO3 ) according to Mackereth et al. (1989); calcium (Ca2+), magnesium (Mg2+), sodium (Na+), potassium (K+), sulphates (SO24 ) and chlorides (Cl-) according to the method of APHA (1992); suspended solids (SS), particulate organic carbon (POC) and total organic carbon (TOC) according to Golterman et al. (1978). The content of dissolved inorganic nitrogen (DIN) was calculated as the sum of ammonium, nitrates and nitrites. Dissolved organic carbon was determined as the difference between TOC and POC. Absorbance values at 250 and 365 nm were recorded for the calculation of E2/E3 and for estimating the humic and fulvic acid ratio in water (Peuravuori and Pihlaja 1997).

Due to a lack of flow data on the Arroyo Durazno, the direct flow corresponding to each rainfall event was estimated with the Instantaneous Unit Hydrograph (IUH) model. The method for estimating the hydrograph consists in calculating effective rainfall for each precipitation event by means of an infiltration method that uses the Curve Number and rainfall value (Soil Conservation Service 1972), and then using the IUH model (Clark 1945) for estimating the flow hydrograph for the stream. The associated uncertainty of the model is due to factors such as intrinsic parameter estimation, lack of information on spatial rainfall distribution, and lack of flow measurements to calibrate it. However, when only rainfall data are available, the IUH model is one of the most commonly used to estimate direct runoff from small basins. The software used was HEC1-HMS 2.0 (Hoggan 1989). The parameters needed to estimate rainfall-flow values are: basin concentration time tc (hours) and the storage constant K (hours), both estimated from the basin’s topographic characteristics. The parameters were calculated according to Wilson and Brown (1992) for use in basins with minimal hydrological data. Basin concentration time was estimated following Kirpich’s (1940) formula. In each sampling, flow (Q) was estimated in situ from channel cross-sectional area and stream velocity. To this end, the time it takes a floating object to travel a given distance was measured (Gordon et al. 1994). These results were compared with those estimated by the IUH model.

Results The highest DO median concentrations (p \ 0.01) were found at DUR and CAS; values lower than 4 mg/ l––the permissible limit for aquatic life protection (USEPA 1986)––were never recorded. Concentrations of SS, nitrates, nitrites, POC and TOC, and EC at both sites did not differ significantly (Table 1). Lower

123

284

GeoJournal (2007) 70:281–288

Table 1 Chemical and physical variables measured at different sites on the Reconquista river Site

DO (mg l-1)

DUR Median

6.9

EC (lS cm-1) 1891

SS (mg l-1) 33

N–NH+4 (lg l-1)

N–NO3 (lg l-1)

403

365

N–NO2 (lg l-1) 70

SRP (lg l-1)

POC (mg l-1)

366

1.3

TOC (mg l-1) 23.2

Min

4.8

352

8

43

117

46

67

0.4

3.8

Max

12.5

2600

111

1219

1038

108

684

3.1

51.8

n

10

CHO Median Min Max n CAS Median

10

10

10

10

10

10

1190

20

4716 *

320

114 *

917 *

0.4

362

8

1170

52

11

324

4.4

1533

70

10262

867

680

3606

9

10

10

10

10

10

1.8 *

10

10 2.10 * 0.8 5.4 10

10 24.1 8.8 36.9 10

7.4

962

54

311 *

378

30

600 *

1.5

Min

4.6

371

17

181

82

7

195

0.9

2.2

Max

9.8

2000

86

1192

1556

112

1161

2.5

29.0

n REY Median

10 4.5 *

9

10

10

10

10

10

1073

35

6323 *

775 *

322 *

1168 *

10

22.0

10

2.08 *

27.0

Min

2.4

973

19

2813

342

180

659

1.2

21.4

Max n

6.8 6

1357 5

610 6

14424 6

2218 6

443 6

1367 6

3.6 6

43.7 6

1.0 *

2.29 *

17.0

GOR Median

1217

31

8264 *

1287 *

390 *

2091 *

Min

0.4

526

17

2565

59

3

854

0.7

2.2

Max

7.8

2350

74

11047

3588

726

4164

4.6

33.7

n SMT Median Min Max n BAN Median

10

9

9

10

10

10

10

1550

72 *

10235 *

150

27

1905 *

0.0

890

44

7250

23

\4

732

0.4

2125

156

22459

530

1164

3477

9

10

10

10

10

10

0.2 *

10

10

6.56 *

52.68 *

4.3

24.4

22.3

80.8

10

10

1545

53 *

11067 *

47 *

20

1826 *

7.63 *

34.36 *

Min

0.0

889

29

6798

12

\4

858

4.3

23.8

Max

0.4

2118

119

14820

220

31

3735

12.4

50.8

9

10

10

10

10

10

10

10

n

0.3 *

10

10

* indicates significant differences between DUR and the other sites

median ammonium concentrations were detected at CAS, followed by those at DUR (Table 1). The highest concentrations were found at SMT and BAN. The low nitrate levels detected at CHO, SMT and BAN coincide with the lower DO levels that promote denitrification. The lowest concentration of o-phosphate was found at DUR (p \ 0.01). Suspended solids were significantly higher at SMT and BAN (p \ 0.01). According to these results, it seems that the best water quality of all sampled sites was that of DUR and CAS, but, since CAS is downstream from the Roggero Reservoir and land use there is mainly urban, we consider that DUR is the site that best represents background conditions for the Reconquista River.

123

The 360 km2-sub-basin of the Arroyo Durazno is 32 km long, and has the lowest population density in the basin. Sixty percent of the land is used as natural pastureland for cattle grazing, and 20% for agriculture (INDEC 1988). The main topographic and geo-morphological characteristics are: area ratio (RA), 3.4; bifurcation ratio (RB), 3.2; length ratio (RL), 1.7; length of the highestorder stream (LX), 26 km; and area of the watershed (AX), 257 km2. Relief consists of very flat plains and low hills (about 20 m a.s.l.) with a 0.05% slope and low-permeability, poorly drained soils. The section under study belongs to secondorder streams, according to the Strahler classification (Gordon et al. 1994).

GeoJournal (2007) 70:281–288

285

16 CN=50-81 14

CN=50-75 CN=60-81 CN=60-75 CN=55- 75

10

3

Direct flow (m /s)

12

8

The storm event between March 1st and 5th, 2001, produced 124 mm of rain. According to the basin’s characteristics, surface flow prevailed during the first hours of the event. Nitrate concentration and EC diminished, SRP and POC concentrations did not change, and DOC and the E2/E3 rate increased before peak flow (Fig. 3). On the receding limb of the hydrograph, nitrate concentration was lower than during peak flow. Concentrations of POC and DOC remained high. However, the E2/E3 rate did not change.

6

4

Discussion

2

Because baseline areas along many surface watercourses are difficult to define, knowledge of background water quality is essential for understanding pollution trends and human impact on rivers (Edmunds et al. 2003). In Argentina, runoff estimates cannot be made due to lack of sufficient rainfall and discharge data from the basins, some of which are slightly to seriously polluted. This is why we propose a methodological approach to water quality assessment in basins that satisfy the hypotheses of the IUH theory (Clark 1945). This methodology involves an analysis of land use and population density, physical and chemical characterization of river water, identification of a background area for the basin, rainfall/ flow modeling with the IUH, relationships between water quality and flow at the background site (water quality dynamics), and water quality monitoring. We compared different sets of data with the median instead of the mean. Several authors (such as Edmunds et al. 2003) recommend using the median because it is more robust and less affected by outliers. When compared to other sites in the Reconquista River basin, water quality at the Arroyo Durazno is most similar to that of natural environments. Rendina et al. (2001) determined that heavy metal mobility in the Arroyo Durazno is the lowest in the basin. As land is used mainly for agriculture and cattle-raising and population density is low, we identified DUR as the background site for the Reconquista River. With an indirect estimation technique, a good approximation was found between flow values determined by the IUH and values estimated at site for the rainfall/flow event. However, direct flow estimation was influenced by the CN variation. McCuen (2002)

53

45

49

41

37

33

29

25

21

17

13

9

5

1

0 Time (h)

Fig. 2 Hydrographs for the rainfall-runoff event observed in March 2001 using different CN values

Parameters used for the IUH of the Arroyo Durazno sub-basin were: tc= 27.8 h and K = 28.9 h. The impact of the Curve Number (CN) variation on direct flow estimation for the Arroyo Durazno was analyzed by means of five tests of the IUH model for the rainfall/flow event of March 2001, using different CN values, each associated with different previous humidity conditions (Fig. 2). The CN values used fall within the confidence intervals suggested by McCuen (2002) for basins of similar characteristics. The analysis of Fig. 2 indicates sensitivity of the form of the hydrographs estimated for the CN range under consideration. The CN values selected for the March 2001 event were 55 and 75, 76 h after the storm event began (Fig. 2). At the background site (DUR), the pH was slightly alkaline (about 8). Sodium was the prevailing cation (Na [[ Ca = Mg [ K) and bicarbonates the prevailing anion in the water (HCO3 [[ Cl [ SO4) (Table 2). Dissolved organic carbon was the main TOC fraction (75–95%). Under high flow conditions, concentrations of TOC and DOC increased, Secchi/ depth values decreased (Table 2), and the water was darker. Concentrations of nitrates, ammonium and major ions, and EC were significantly lower under high flow conditions.

123

286

GeoJournal (2007) 70:281–288

Table 2 Chemical and physical variables measured under different hydrological conditions at DUR Low flow Median DO (mg l-1) Temperature (°C) EC (lS cm-1) Secchi:Depth (%) SS (mg l-1) pH

High flow Min

Max

n

Median

Min

Max

n

6.8

4.5

12.5

17

5.9

3.9

8.1

12

18.0

6.0

29.1

17

20.2

7.4

24.4

12

1650

850

2600

17

349

229

785

12

100

60

100

17

27

15

53

12

26

6

111

17

32

19

86

12

8.0

7.5

8.4

17

7.5

7.1

8.2

12

N–NH+4 (lg l-1)

88

19

1219

17

49

15

1205

12

-1 N–NO3 (lg l )

324

35

3121

17

116

11

402

12

41

3

108

17

25

0

108

12

235

70

684

17

211

118

561

N–NO2

-1

(lg l )

SRP (lg l-1) POC (mg l-1)

1.0

TOC (mg l-1)

17.2

Ca+2 (mg l-1)

33

0.1 3.8 12

2.7

17

1.2

37.9

17

34.6

45

17

6

0.4

3.1

12 12

6.9

67.4

12

3

13

12

Mg+2 (mg l-1)

21

11

33

17

4

2

6

12

Na+ (mg l-1) K+ (mg l-1)

259 15

198 12

476 18

17 17

68 10

22 8

193 14

12 12

-1 HCO3 (mg l )

539

411

745

17

178

118

528

12

81

46

183

17

14

1

70

12

174

95

317

17

22

12

44

-1 SO24 (mg l )

Cl- (mg l-1) Q (m3 s-1)

0.07

0.04

0.70

pointed out that although the CN is treated as a constant, it varies for each observed rainfall event, which means that the confidence intervals are inadequate to assess such variation. At the background site and under high- and lowflow conditions, the ionic composition of the riverwater is consistent with soil characteristics in the sub-basin, where almost 50% have a natric horizon. Data on chlorides are considered as good tracers for surface water source identification (Ribolzi et al. 2000). During dry periods, chloride concentrations and EC were similar to those determined by EASNE (1973) for groundwater, thus showing that under these hydrologic conditions the water table is the predominant source of stream water. The increase in TOC and DOC concentrations during high-flow conditions suggests a significant input of carbon from runoff. Meybeck (1982) reported high concentrations of TOC in other rivers with poorly drained areas. The main DOC fraction in water from the soil is represented by humic substances, as evidenced by the dark color of water during high-flow

123

17

3.85

2.67

8.60

12 12

conditions and the low Secchi/depth values. Also, the increase in the E2/E3 ratio shows an increase in the fulvic and humic acid proportion during the high-flow period. In the presence of fulvic acids, heavy metal solubility, mobility, environmental prevalence, and toxicity increase (Varshal et al. 2000), thus increasing pollution risk during high-flow conditions. Besides, nitrate, ammonium and major ion concentrations and EC decreased due to the dilution effect of the massive inflow of water into the stream. On the receding limb of the hydrograph, nitrate concentration was lower than during maximum flow and organic carbon concentration remained high. This behavior, known as the ‘‘flushing effect’’ (Mitchell et al. 1997), suggests that the soluble material accumulated in the drainage area during dry periods is transported to the stream by lixiviation and surface runoff, thus raising solute concentrations during the first storm hours. Water quality in rivers of low to moderate flow changes rapidly, even in background areas due to its dependence on the flow.

GeoJournal (2007) 70:281–288 0

Rain (mm)

Fig. 3 Hydrographs and physico-chemical variables corresponding to the rainfall-runoff event of March 2001

287

20

40

60 1,4

10

9

9

8

8

7

7

1

6

0,8

5 4 3

5 4

0,6

3

0,4

2

2

1

1

0

0 50

100

150

200

0,2

0

250

Discharge in situ

POC

10

40

1600

9

35

8

1400

8

1200

6

1000

5 800

4

600

3

d Discharge (m3 s-1)

7

6

25

5

20

4

15

3

2

2

1

200

1

0

0

0

50

100

150

200

30

7

400

250

10 5 0

50

EC

8

600

7 6

500

5

400

4

300

3

200

2

f

2,5

5

2

4

1,5

3

0

At the same time, a lack of rainfall/flow records in many parts of the world complicates the integration of hydrologic and geochemical observations. In view of this situation, our proposal here is a contribution towards an improved description of processes for proper surface water quality management based on eco-hydrological processes.

1

2 1

Rainfall event (hours) Discharge (IUH model) N-NO3

3

6

0 250

3,5

7

100 200

4

9

0 150

DOC

10

1 100

0 250

200

8 Discharge (m3 s-1)

700 N-NO3- and SRP (µg l -1)

800

9

50

150

Discharge (IUH model)

10

0

100

Rainfall event (hours)

Discharge (IUH model)

3 -1

0 250

1800

Rainfall event (hours)

Discharge (m s )

200

9

0

e

150

10

EC (µS cm -1)

Discharge (m3 s-1)

c

100

Rainfall event (hours) Discharge (IUH model)

Rainfall event (hours) Discharge (IUH model)

50

E2/E3

0

POC (mg l -1)

6

1,2

DOC (mg l-1)

b

10

Discharge (m3 s-1)

Discharge (m3 s-1)

a

0,5 0

50

100

150

200

0 250

Rainfall event (hours) SRP

Discharge (IUH model)

E2/E3

Acknowledgments The authors wish to acknowledge the logistic support of the Prefectura Naval Argentina during sampling operations.

References APHA. (1992). Me´todos normalizados para el ana´lisis de aguas potables y residuales. Madrid.

123

288 Castan˜e´, P. M., Loez, C., Olguı´n H., Puig, A., Rovedatti, M. G., Topalia´n, M., & Salibia´n, A. (1998). Caracterizacio´n y variacio´n espacial de para´metros fisicoquı´micos y del plancton en un rı´o contaminado (Rı´o Reconquista, Argentina). Revista Internacional de la Contaminacio´n Ambiental, 14(2), 69–77. Clark, C. O. (1945). Storage and the unit hydrograph. Transactions of American Society of Civil Engineers, 110, 1419–1488. de Cabo, L., Arreghini, S., Fabrizio de Iorio, A., Rendina, A., Bargiela, M., Vella, R., & Bonetto, C. (2000). Impact of the Moro´n stream on water quality of the Reconquista River (Buenos Aires, Argentina). Revista del Museo Argentino Ciencias Naturales, n.s. 2(2), 123–130. EASNE (Estudios de Aguas Subterra´neas del Nor-Este de la Provincia de Buenos Aires). (1973). Contribucio´n al estudio geohidrolo´gico del noreste de la Provincia de Buenos Aires. Consejo de Estudios de Aguas Subterra´neas Nor-Este. Serie Te´cnica N°24. Buenos Aires. 144 pp. Edmunds, W. M., Shand, P., Hart, P., & Ward, R. S. (2003). The natural (baseline) quality of groundwater: A UK pilot study. Science of the Total Environment, 310, 25–35. Golterman, H., Clymo, R., & Ohndtad, M. (1978). Methods for the physical and chemical examination of freshwaters. Oxford. Gordon, B. N. D., McMahon, T. A., & Finlayson, B. L. (1994). Stream hydrology. An introduction for ecologists. New York. Hoggan, D. (1989). Computer assisted floodplain hydrology and hydraulics. McGraw-Hill Publishing Company, 518 pp. INDEC (Instituto Nacional de Estadı´stica y Censos). (1988). Censo Nacional Agropecuario. Resultados GeneralesProvincia de Buenos Aires. Argentina: Buenos Aires. Kirpich, Z. P. (1940). Time of concentration of small agricultural watersheds. Civil Engineering American Society, 10, 362. Mackereth, F., Heron, J., & Talling, J. (1989). Water analysis: Some revised methods for limnologists, 2nd ed. UK: Cumbria. McCuen, R. H. (2002). Approach to confidence interval estimation for curve numbers. Journal of Hydrologic Engineering, 7(1), 43–48. Meybeck, M. (1982). Carbon, nitrogen, and phosphorus transport by world rivers. American Journal of Science, 282, 401–450.

123

GeoJournal (2007) 70:281–288 Mitchell, A. W., Bramley, R. G. V., & Johnson, A. K. L. (1997). Export of nutrients and suspended sediment during a cyclone-mediated flood event in the Herbert River catchment, Australia. Marine and Freshwater Research, 48, 79–88. Olson, D. E., Dinerstein, E., Canevari, P., Davidson, I., Castro, G., Morisset, V., Abell, R., & Toledo, E. (1998). Freshwater biodiversity of Latin America and the Caribbean: A conservation assessment. Washington DC. Petts, G., & Calow, P. (Eds). (1996). River restoration. Oxford. Peuravuori, J., & Pihlaja, K. (1997). Molecular size distribution and spectroscopic properties of aquatic humic substances. Analytica Chimica Acta, 337, 133–149. Rendina, A., Iorio, A.F. de, Bargiela, M., de Cabo, L., & Arreghini, S. (2001). Geochemical distribution and mobility factors of Zn and Cu in sediments of the Reconquista river (Argentina). Revista Internacional de la Contaminacio´n Ambiental, 17(4), 187–192. Ribolzi, O., Andrieux, P., Valles, V., Bouzigues, R., Bariac, T., & Voltz, M. (2000). Contribution of groundwater and overland flows to storm flow generation in a cultivated Mediterranean catchment. Quantification by natural chemical tracing. Journal of Hydrology, 233, 241–257. Soil Conservation Service. (1972). National engineering handbook. Washington, DC. Strickland, J., & Parsons, T. (1972). A practical handbook of seawater analysis. 2nd ed. Bull. No. 167, Ottawa, Canada: Fisheries Research Board. USEPA. (1986). Quality criteria for water. Washington, DC: United States Environmental Protection Agency, 440/586-001. Varshal, G., Velyukhanova, T., Khushvakhtova, S., Koshcheeva, I., & Tyutyunnik, O. (2000). The role of humus acids in the processes of dispersion and accumulation of heavy metals in natural waters, soils and carbonaceous rocks. International Conference on heavy metals in the Environment, 2000, Ann Arbor, MI. Wilson, B., & Brown, W. (1992). Development and evaluation of a dimensionless unit hydrograph. Water Resources Bulletin, 28(2), 397–408. Zalewski, M. (2002). Ecohydrology-the use of ecological and hydrological processes for sustainable management of water resources. Hydrological Sciences Journal, 47(5), 823–832.