**Current Challenges in Experimental Watershed Hydrology**

Wei-Zu Gu, Jiu-Fu Liu, Jia-Ju Lu and Jay Frentress

Additional information is available at the end of the chapter

http://dx.doi.org/10.5772/55087

## **1. Introduction**

The river basin, watershed or catchment is central to many of concepts in hydrology [1] including contaminant hydrology. In fact, two different studies in the Seine River basin during the end of seventeenth century, independently made by Pierre Perraut and Edme Mariotte have been identified by A. K. Biswas as the beginning of quantitative hydrology [2], Along with Perraut and Mariotte, Biswas also suggested Edmond Halley as co-founder of experi‐ mental hydrology. It was later accepted that scientific hydrology was founded on these two basin studies [3].

However, basin studies developed slowly until the end of nineteenth century when public demands accelerated. The first modern basin studies commenced in Emmental in Switzer‐ land during the 1890s and were focused on hydrological differences between two small, ca 60 ha, basins [1]. A multitude of basin studies have appeared since the early twentieth century from many parts of the world. These include: (1) Wagon Wheel Gap experiment of USA begun in 1910 in two forested basins [4]; (2) Valday Hydrological Laboratory of USSR begun in 1933 and focused on field investigations of multiple hydrological parame‐ ters in watersheds with different scales [5]; (3) Coweeta Hydrologic Laboratory of USA set up since 1934 for forest hydrology and, ecological research in two main basins with 32 subwatersheds under various treatments [6]; (4) Harz Mountains experiment of Germany with a pair of catchments, Wintertal and Large Bramke, begun in 1948 and focusing on hydro‐ logical effects of land use [7]; (5) Alrance experiment of France started in 1950 and focused on streamflow patterns[8]; and (6) Bluebrook Runoff Experiment of China established on 1953 and focused on drainage problems of the vast agricultural plain of the North Huai He River [9].

© 2013 Gu et al.; licensee InTech. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. © 2013 Gu et al.; licensee InTech. This is a paper distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

A period of rapid worldwide development of hydrological basin studies resulted from the International Hydrology Decade (IHD) Representative and Experimental Basin Programme 1965- 1974, with an estimated 3000 basin studies conducted during the Decade [1]. "Repre‐ sentative basins (RBs), which are selected as representative of a region of presumed hydro‐ logical similarity, are used for intensive investigations of specific problems of the hydrological cycle (or part thereof) under relatively stable, natural conditions" [10]. Experimental basins (EBs) are relatively homogeneous in soil and vegetation, have uniform physical characteristics, and are deliberately modified for study [10]. However, many EBs are just well-instrumented basins of small area [11]. Sometimes both RB and EB are included, as in the fruitful Plynlimon experiment in the UK [12].

**2. Challenges**

using experimental watershed studies.

**2.1. Demands of advancing hydrology science**

Numerous challenges face researchers as they move forward to advance hydrology science

Current Challenges in Experimental Watershed Hydrology

http://dx.doi.org/10.5772/55087

301

More than two decades ago, V. Klemes proclaimed that "hydrology is as yet lacking a solid scientific foundation needed for its development as a natural science" [17] and this sentiment has been reiterated by others. Bras and Eagleson said "in the modern science establishment, this niche is vacant" in their paper titled "Hydrology, the fogotten earth science" [18]. Jim Dooge asked "Is hydrology now an established science? Is hydrologic practice now firmly based on scientific principles?"[19]. McDonnell et al concluded that "the critical ideas and positive vision presented in that (Dooge) paper remain just as fresh, relevant and, unfortu‐ nately, very much unfulfilled" [13]. Sivapalan also suggested several steps toward a new unified hydrologic theory [14]. Substantial progress in hydrologic science "ultimately depends on new experimental work, new field observations, and new data collection networks" as summarized by Kirchner [15]. But, what shape will these new experiments and networks take?

*Coupling processes.* A watershed is an inherently dynamic ecological system composed of a variety of biotic and abiotic processes. For such systems, it is important not only to focus on hydrological process but also on linked ecological, biogeochemical, pedological and geological processes. Misunderstandings of the interrelatedness of these processes are commonplace (e.g., summarizing the net effect of all biological processes as evapotranspiration). The vegetation root system forms a dynamically complex net for preferential flowpaths. It changes runoff generation mechanisms not only for surface runoff but the subsurface flows, and the root system itself is effected via feedback from the various, changing flow patterns. Vegetation growth also adapts according to climate, infiltration, soil water and even evapotranspiration while feedback mechanisms reinforce the interconnectedness of the watershed system [20].

*Innovative measurements.* More than twenty years ago, JE Nash suggested that discovery in hydrology sciences has been limited by "a deficiency in our empiricism" and, "our tolerance of poor methods of observations" [21]. Since then, methods of observations have improved significantly. However, our observation techniques still inhibit hydrology and a fundamental change to systematic measurement programs is needed. The current observation programs are aimed mostly at classic watershed hydrological process. Instead we should focus on a holistic description of heterogeneity of not only the watershed hydrological process but all related

*Supporting and testing hydrological model.* Physically based hydrological models are one of the most promising directions for advancing hydrologic science. Grayson et al concluded that "the models are enabling us to ask more questions, many of which are fundamental to our under‐ standing of the natural systems that are the subject of our models" [22]. However, vital weaknesses in current watershed models exist. First, they are generally based on "well known small-scale theories such as Darcy's law and the Richards equation for coupled balance of mass

processes as described above, as well as landscape properties and climate inputs.

If we designate the first phase of basin study (until ca the middle twentieth century) as foundational and the second phase (during/after the IHD) as developmental, it would appear that experimental watershed hydrology is inevitably going into a third phase of transition and innovation. Most field experiments in watershed science to date remain largely descriptive, with results that are difficult to generalize [13]. In effect, "… catchment hydrology is trapped in a dead-end track, a theoretical impasse" [14]. Experimental watershed studies, the core of watershed hydrology, are now confronting tremendous and even more complicated challeng‐ es, which are raised mainly from changing environment conditions and partly from the weakness of current watershed studies. This chapter provides an overview of the main challenges facing experimental watershed hydrology. Addressing these challenges will hopefully lead to substantial innovation in the field.

An important and instructive paper addressing a new vision of watershed hydrology [13] asked "What's wrong with the status quo?" If we posed this same question to the experimental hydrologist, the limitations of most existing RBs/EBs might be characterized as mainly twofold. First, study basin designs have been limited by the black box concept and many misconceptions (e.g., the linearity, non-heterogeneity, additivity of hydrologic systems etc.). Second, operation has been substantially bounded by the hydraulic conception of these watersheds as isolated hydrological systems. All of these watershed studies monitor only total runoff at the streamoutlet and the subsurface responses of the watershed are only estimated by hydrograph separation, etc. These characteristics undermine the formulation of a unified theory of watershed hydrology [14] and the development of watershed models [13,15].

There is a clear need to move beyond the status quo and expand from this narrow hydrological perspective to generate hypotheses governing general behavior across places and scales [13], with the ultimate aim "to advance the science of hydrology" [15]. For the third phase, the stage of transition and innovation, instead of classic RBs as described above, another kind of experimental watershed, the Critical Zone Experimental Block (CZEB) is suggested, with its concept and infrastructures very different from that of current EB. Following the concept of Critical Zone defined by the National Research Council [16], the CZEB will have practical boundaries for different geomorphologic regions.

An ongoing trial application of the CZEB approach is presented in this chapter.

## **2. Challenges**

A period of rapid worldwide development of hydrological basin studies resulted from the International Hydrology Decade (IHD) Representative and Experimental Basin Programme 1965- 1974, with an estimated 3000 basin studies conducted during the Decade [1]. "Repre‐ sentative basins (RBs), which are selected as representative of a region of presumed hydro‐ logical similarity, are used for intensive investigations of specific problems of the hydrological cycle (or part thereof) under relatively stable, natural conditions" [10]. Experimental basins (EBs) are relatively homogeneous in soil and vegetation, have uniform physical characteristics, and are deliberately modified for study [10]. However, many EBs are just well-instrumented basins of small area [11]. Sometimes both RB and EB are included, as in the fruitful Plynlimon

300 Current Perspectives in Contaminant Hydrology and Water Resources Sustainability

If we designate the first phase of basin study (until ca the middle twentieth century) as foundational and the second phase (during/after the IHD) as developmental, it would appear that experimental watershed hydrology is inevitably going into a third phase of transition and innovation. Most field experiments in watershed science to date remain largely descriptive, with results that are difficult to generalize [13]. In effect, "… catchment hydrology is trapped in a dead-end track, a theoretical impasse" [14]. Experimental watershed studies, the core of watershed hydrology, are now confronting tremendous and even more complicated challeng‐ es, which are raised mainly from changing environment conditions and partly from the weakness of current watershed studies. This chapter provides an overview of the main challenges facing experimental watershed hydrology. Addressing these challenges will

An important and instructive paper addressing a new vision of watershed hydrology [13] asked "What's wrong with the status quo?" If we posed this same question to the experimental hydrologist, the limitations of most existing RBs/EBs might be characterized as mainly twofold. First, study basin designs have been limited by the black box concept and many misconceptions (e.g., the linearity, non-heterogeneity, additivity of hydrologic systems etc.). Second, operation has been substantially bounded by the hydraulic conception of these watersheds as isolated hydrological systems. All of these watershed studies monitor only total runoff at the streamoutlet and the subsurface responses of the watershed are only estimated by hydrograph separation, etc. These characteristics undermine the formulation of a unified theory of

There is a clear need to move beyond the status quo and expand from this narrow hydrological perspective to generate hypotheses governing general behavior across places and scales [13], with the ultimate aim "to advance the science of hydrology" [15]. For the third phase, the stage of transition and innovation, instead of classic RBs as described above, another kind of experimental watershed, the Critical Zone Experimental Block (CZEB) is suggested, with its concept and infrastructures very different from that of current EB. Following the concept of Critical Zone defined by the National Research Council [16], the CZEB will have practical

watershed hydrology [14] and the development of watershed models [13,15].

An ongoing trial application of the CZEB approach is presented in this chapter.

experiment in the UK [12].

hopefully lead to substantial innovation in the field.

boundaries for different geomorphologic regions.

Numerous challenges face researchers as they move forward to advance hydrology science using experimental watershed studies.

## **2.1. Demands of advancing hydrology science**

More than two decades ago, V. Klemes proclaimed that "hydrology is as yet lacking a solid scientific foundation needed for its development as a natural science" [17] and this sentiment has been reiterated by others. Bras and Eagleson said "in the modern science establishment, this niche is vacant" in their paper titled "Hydrology, the fogotten earth science" [18]. Jim Dooge asked "Is hydrology now an established science? Is hydrologic practice now firmly based on scientific principles?"[19]. McDonnell et al concluded that "the critical ideas and positive vision presented in that (Dooge) paper remain just as fresh, relevant and, unfortu‐ nately, very much unfulfilled" [13]. Sivapalan also suggested several steps toward a new unified hydrologic theory [14]. Substantial progress in hydrologic science "ultimately depends on new experimental work, new field observations, and new data collection networks" as summarized by Kirchner [15]. But, what shape will these new experiments and networks take?

*Coupling processes.* A watershed is an inherently dynamic ecological system composed of a variety of biotic and abiotic processes. For such systems, it is important not only to focus on hydrological process but also on linked ecological, biogeochemical, pedological and geological processes. Misunderstandings of the interrelatedness of these processes are commonplace (e.g., summarizing the net effect of all biological processes as evapotranspiration). The vegetation root system forms a dynamically complex net for preferential flowpaths. It changes runoff generation mechanisms not only for surface runoff but the subsurface flows, and the root system itself is effected via feedback from the various, changing flow patterns. Vegetation growth also adapts according to climate, infiltration, soil water and even evapotranspiration while feedback mechanisms reinforce the interconnectedness of the watershed system [20].

*Innovative measurements.* More than twenty years ago, JE Nash suggested that discovery in hydrology sciences has been limited by "a deficiency in our empiricism" and, "our tolerance of poor methods of observations" [21]. Since then, methods of observations have improved significantly. However, our observation techniques still inhibit hydrology and a fundamental change to systematic measurement programs is needed. The current observation programs are aimed mostly at classic watershed hydrological process. Instead we should focus on a holistic description of heterogeneity of not only the watershed hydrological process but all related processes as described above, as well as landscape properties and climate inputs.

*Supporting and testing hydrological model.* Physically based hydrological models are one of the most promising directions for advancing hydrologic science. Grayson et al concluded that "the models are enabling us to ask more questions, many of which are fundamental to our under‐ standing of the natural systems that are the subject of our models" [22]. However, vital weaknesses in current watershed models exist. First, they are generally based on "well known small-scale theories such as Darcy's law and the Richards equation for coupled balance of mass and momentum" [13] and do not accurately capture the spatial heterogeneity, inherent nonlinearity, and non-additive properties of natural watershed systems. Second, a large number of models are heavily over-parameterized, leading to equi-finality, wherein multiple combi‐ nations of parameters can yield equivalent results [13]. Grayson et al showed that comparable fits to a hydrograph could be achieved using a saturated overland flow model or Hortonian flow assumptions [22]. Kirchner argued that "parameter-rich models may succeed as mathe‐ matical marionettes dancing to match the calibration data even if their underlying premises are unrealistic" [15]. Kirchner concluded that if "the present trend away from physical processes and toward mathematistry ('blackboard hydrology') continues in hydrologic education and practice, hydrology will end up in a dead end as a science and become useless for applications" [17].

*Consequences of replumbed hydrology.* Monitoring and unraveling the hydrological and ecolog‐ ical responses and the feedbacks due to the replumbed hydrological cycle from water con‐ servancy projects, especially those in arid and semi-arid regions, is needed. Dams are perhaps the most dramatic examples of human capacity to transform nature in the name of develop‐ ment [27], but development programs and projects create both winners and losers [28]. Dams in arid and semi-arid basins, especially in endorheic basins, break the natural cycle, looting downstream groundwater and accelerating desertification [29]. Long distance water diversion may reasonably be viewed as an ecological "planning disaster"[30]. There are concerns that development is outpacing scientific understanding, a circumstance that may prove disastrous

Current Challenges in Experimental Watershed Hydrology

http://dx.doi.org/10.5772/55087

303

*Fundamental studies.* Various hydrological and biogeochemical fluxes through catchments are important to understanding the effects of hydrologic replumbing, but many of these are poorly understood and rarely measured. For example, various subsurface ecological and hydrological components cannot be measured directly and must instead be estimated from a small number of point measurements. Particular emphasis has previously been placed on the downward flux of catchment processes components (i.e., rainfall and discharge) rather than upward flux- land evaporation, transpiration, and upward recharge of groundwater flux, etc – even though "these fluxes serve as important regulators of the dynamics of the cycle"[23]. There are various unknown or poorly understood mechanisms regarding natural catchment behavior that are affected by the hydrologic replumbing occurring at local, regional and global scales. These include: (1) the old water paradox [31]; (2) network-like preferential flow[13], natural symmetry between water, landscape, soil and vegetation and the underlying organizing principles between soil, vegetation and other biotic elements [14]; (3) the combined mecha‐ nisms of the geological water-cycle and hydrologic water cycle in the shallow hydrosphere, as well as the deep circulation of groundwater, which go far beyond current hydrogeological boundaries [32]; (4) the puzzle of the missing sink for the remaining 40% of watershed carbon balance, which after years effort a consensus for it still has yet to emerge [33]; (5) the mecha‐ nisms involved in low flow hydrology and its relation to shallow and deep groundwater, regolith, land use, water quality and the biodiversity of aquatic ecosystems [34]; and (6) the coupling of basin geomorphologic processes with hydrologic and ecologic processes, includ‐ ing the poorly understood processes and mechanisms of environmental release, transport, and biological transformation of various contaminants in the unsaturated zone, shallow aquifers

*Historical data on water replumbling.* Historical data from sources other than hydrometric stage records will be very helpful to understand large-scale environmental changes, serving as a reference for the natural variation and anthropogenic impact. Figure 1 illustrates desertifica‐ tion in an endorheic basin as the result of human impacts. Four stages of man-earth relationship are illustrated, from the 'stage of natural harmony' until the stage of the environmental

*Design systematic experimental facilities.* The dynamic systems currently facing experimental hydrologists exhibit a wide range of heterogeneity and process complexity across large spatiotemporal scales. There are multiple ways to overcome such scaling problems [13,14], but

disasters, the punishment of nature (see Figure 1 caption).

for the environment for the generations to come.

and deep aquifers.

To achieve substantive improvements, models should be developed "in conjunction with, and as an integral part of, carefully planned and executed field studies established for the purpose of advancing our knowledge of the natural system"[22]. Key experimental basin requirements for model-generated hypothesis testing include: use of innovative basin designs and inclusion of sufficient variation in physiographic settings, hydroclimates, and watershed scales.

## **2.2. Hydrologic replumbing and natural climate oscillations**

The National Research Council of the US identified anthropogenic perturbation and replumb‐ ing of the hydrologic cycle [23] as a fundamental challenge [16] in watershed science. Human influences dominate "the natural cycle of freshwater causing environmental changes that have been argued to move the planet to a new geologic era termed the 'Anthropocene'" [23,24].

Natural climate oscillations with observed large, abrupt events, the widespread millennial scale climate changes of the last glaciation (consisting of Dansgaard-Oeschger oscillations, Heinrich events etc.) are hypothesized to have been forced by North Atlantic atmosphereocean-ice interactions [25]. It has been argued that these oscillations are ultimately driven by variations in eccentricity, axial tilt and precession of the Earth's orbit as described in the Milankovitch theory. The decadal and multidecadal climate variations are generally linked with recognized dynamics, such as the El Niño Southern Oscillation (ENSO) [23].

Anthropogenic modification of the water cycle could "push the climate into new regimes" because "climate change has taken the climate system out of the repeated cycle of glacialinterglacial episodes" [23]. These alterations may accelerate the arrival of millennial scale events and trigger a "tipping point" transition resulting in "major climatic perturbations on time scales of decades to centuries"[26]. How might experimentalists, a small proportion of total hydrologists, address these challenges with progressive basin studies?

*Long term monitoring.* Well-established representative basins and benchmark experimental basins distributed across different physiographic and hydroclimatic conditions are critical for assessing hydrologic changes due to replumbing. Basins instrumented during or before the 1950's may also be used to determine patterns of rainfall redistribution, resultant stream flow and related biological issues.

*Consequences of replumbed hydrology.* Monitoring and unraveling the hydrological and ecolog‐ ical responses and the feedbacks due to the replumbed hydrological cycle from water con‐ servancy projects, especially those in arid and semi-arid regions, is needed. Dams are perhaps the most dramatic examples of human capacity to transform nature in the name of develop‐ ment [27], but development programs and projects create both winners and losers [28]. Dams in arid and semi-arid basins, especially in endorheic basins, break the natural cycle, looting downstream groundwater and accelerating desertification [29]. Long distance water diversion may reasonably be viewed as an ecological "planning disaster"[30]. There are concerns that development is outpacing scientific understanding, a circumstance that may prove disastrous for the environment for the generations to come.

and momentum" [13] and do not accurately capture the spatial heterogeneity, inherent nonlinearity, and non-additive properties of natural watershed systems. Second, a large number of models are heavily over-parameterized, leading to equi-finality, wherein multiple combi‐ nations of parameters can yield equivalent results [13]. Grayson et al showed that comparable fits to a hydrograph could be achieved using a saturated overland flow model or Hortonian flow assumptions [22]. Kirchner argued that "parameter-rich models may succeed as mathe‐ matical marionettes dancing to match the calibration data even if their underlying premises are unrealistic" [15]. Kirchner concluded that if "the present trend away from physical processes and toward mathematistry ('blackboard hydrology') continues in hydrologic education and practice, hydrology will end up in a dead end as a science and become useless

To achieve substantive improvements, models should be developed "in conjunction with, and as an integral part of, carefully planned and executed field studies established for the purpose of advancing our knowledge of the natural system"[22]. Key experimental basin requirements for model-generated hypothesis testing include: use of innovative basin designs and inclusion of sufficient variation in physiographic settings, hydroclimates, and watershed scales.

The National Research Council of the US identified anthropogenic perturbation and replumb‐ ing of the hydrologic cycle [23] as a fundamental challenge [16] in watershed science. Human influences dominate "the natural cycle of freshwater causing environmental changes that have been argued to move the planet to a new geologic era termed the 'Anthropocene'" [23,24].

Natural climate oscillations with observed large, abrupt events, the widespread millennial scale climate changes of the last glaciation (consisting of Dansgaard-Oeschger oscillations, Heinrich events etc.) are hypothesized to have been forced by North Atlantic atmosphereocean-ice interactions [25]. It has been argued that these oscillations are ultimately driven by variations in eccentricity, axial tilt and precession of the Earth's orbit as described in the Milankovitch theory. The decadal and multidecadal climate variations are generally linked

Anthropogenic modification of the water cycle could "push the climate into new regimes" because "climate change has taken the climate system out of the repeated cycle of glacialinterglacial episodes" [23]. These alterations may accelerate the arrival of millennial scale events and trigger a "tipping point" transition resulting in "major climatic perturbations on time scales of decades to centuries"[26]. How might experimentalists, a small proportion of

*Long term monitoring.* Well-established representative basins and benchmark experimental basins distributed across different physiographic and hydroclimatic conditions are critical for assessing hydrologic changes due to replumbing. Basins instrumented during or before the 1950's may also be used to determine patterns of rainfall redistribution, resultant stream flow

with recognized dynamics, such as the El Niño Southern Oscillation (ENSO) [23].

total hydrologists, address these challenges with progressive basin studies?

**2.2. Hydrologic replumbing and natural climate oscillations**

302 Current Perspectives in Contaminant Hydrology and Water Resources Sustainability

for applications" [17].

and related biological issues.

*Fundamental studies.* Various hydrological and biogeochemical fluxes through catchments are important to understanding the effects of hydrologic replumbing, but many of these are poorly understood and rarely measured. For example, various subsurface ecological and hydrological components cannot be measured directly and must instead be estimated from a small number of point measurements. Particular emphasis has previously been placed on the downward flux of catchment processes components (i.e., rainfall and discharge) rather than upward flux- land evaporation, transpiration, and upward recharge of groundwater flux, etc – even though "these fluxes serve as important regulators of the dynamics of the cycle"[23]. There are various unknown or poorly understood mechanisms regarding natural catchment behavior that are affected by the hydrologic replumbing occurring at local, regional and global scales. These include: (1) the old water paradox [31]; (2) network-like preferential flow[13], natural symmetry between water, landscape, soil and vegetation and the underlying organizing principles between soil, vegetation and other biotic elements [14]; (3) the combined mecha‐ nisms of the geological water-cycle and hydrologic water cycle in the shallow hydrosphere, as well as the deep circulation of groundwater, which go far beyond current hydrogeological boundaries [32]; (4) the puzzle of the missing sink for the remaining 40% of watershed carbon balance, which after years effort a consensus for it still has yet to emerge [33]; (5) the mecha‐ nisms involved in low flow hydrology and its relation to shallow and deep groundwater, regolith, land use, water quality and the biodiversity of aquatic ecosystems [34]; and (6) the coupling of basin geomorphologic processes with hydrologic and ecologic processes, includ‐ ing the poorly understood processes and mechanisms of environmental release, transport, and biological transformation of various contaminants in the unsaturated zone, shallow aquifers and deep aquifers.

*Historical data on water replumbling.* Historical data from sources other than hydrometric stage records will be very helpful to understand large-scale environmental changes, serving as a reference for the natural variation and anthropogenic impact. Figure 1 illustrates desertifica‐ tion in an endorheic basin as the result of human impacts. Four stages of man-earth relationship are illustrated, from the 'stage of natural harmony' until the stage of the environmental disasters, the punishment of nature (see Figure 1 caption).

*Design systematic experimental facilities.* The dynamic systems currently facing experimental hydrologists exhibit a wide range of heterogeneity and process complexity across large spatiotemporal scales. There are multiple ways to overcome such scaling problems [13,14], but hydrometric measurement systems in natural basins still face challenges (e.g., the ability of current rain gauges with conventional deployment techniques to obtain precise precipitation inputs for natural basins remains controversial). It is necessary to have systematic experimental facilities consisting of laboratory physical models, hillslope "catchment", small-size natural watershed with uniformly dominant vegetation, small-size natural watershed with more complex vegetation, unchannelized catchments, nested sub-watersheds, and watersheds of large size. Such experimental facilities need to be situated in zones with different physio‐ graphic and hydroclimatic conditions. Obviously, this would require collaborative efforts at an international scale.

**3.1. Deficiencies in the design and operation of current experimental basins**

*Watershed surface laterality.* Ignoring the inherent interconnectedness of surface and subsurface watershed components reflects the fact that watershed study is driven by practical rather than scientific requirements. Hydrologists still lack a comprehensive understanding of surface responses to nonvisible and difficult to measure subsurface constraints. A physiological analogy might be studying the skin while ignoring the functions beneath. Focusing on surface responses using surface monitoring will quickly trap the observer with misconceptions. For instance, hydrologists delineate watersheds with surface topography, though actual water‐ shed boundaries may cross topographically determined boundaries depending on subsurface features. Even our most fundamental measurements of the watershed – its shape and size –

Current Challenges in Experimental Watershed Hydrology

http://dx.doi.org/10.5772/55087

305

are strongly influenced by our inability to peer beneath the surface of the catchment.

while only looking at blood vessels.

dominant terms in water balance [37].

controversial results.

*Hydrologic process laterality.* The hydrologic cycle is tightly coupled with other cycles [36]. The hydrologic cycle is also dynamic, not static, in nature. Studying the hydrologic cycle while ignoring the interconnectedness of other cycles is equivalent to studying blood circulation

*Downward components laterality.* There is a tendency to emphasize measurements of downward components (e.g. rainfall, discharge and infiltration) rather than upward components (e.g. evaporation and transpiration aforementioned), even though upward components can be the

**3.2. Lessons learned from the Chinese experimental basin studies over the last fifty years** The Chinese experimental watershed studies experienced a saddle-backed fluctuation with two peaks [9] during their fifty-year history. The lessons learned were summarized in

**1.** The research facilities need not only the RB and EB of natural conditions, but also those with controlled boundary conditions. All EBs of recent decades referred only to the natural boundary of the surface watershed, ignoring that of the bedrock. This frequently led to

**3.** The runoff response of an EB, including surface and subsurface components, should be monitored hydrometrically. For all EBs in past decades, these components were moni‐ tored incompletely or incorrectly, insufficient even for 'black box' evaluations. As a

**4.** It is necessary to use multiple tracers to look inside the physical processes of the hydro‐

**5.** Comparative EBs should be located in different climate and morphological regions and

It was also suggested that hydrological experimentation should not address "only to the hydrosphere but the interactions between atmosphere, lithosphere, biosphere and the intelligence-sphere." To this end, there is a need to "develop approaches using physical,

cooperation with J.J. McDonnell, C. Kendall and N.E. Peters, as following [9]:

**2.** The EB should be treated as an integrated system, from surface to bedrock.

consequence, most of the results can't be explained physically.

should have comparable monitoring strategies.

logical cycle and the anthropogenic impacts at the catchment scale.

**Figure 1.** Desertification of the Rocksheep River (Shiyang River) Basin [35], an endorheic basin with total drainage area of 40,690 km2. (a) prior to 121 BC, the "stage of natural harmony", the area of its terminal lake around 540 km2; (b) 420- 589 AD, the "stage of relying"; (c) ca 1900AD, the "stage of impact", the cultivated area during Qing Dynasty was developed quickly, which reached to an area about 4/5 of that of modern time. (d) ca 1990's, the "stage of dam‐ age", there are 4 dams since 1959, 22 dams until 1972 [35].

## **3. The Critical Zone Experimental Block (CZEB)**

For the future of hydrologic basin study, another kind of experimental basin is suggested: the Critical Zone Experimental *Block* (CZEB) customarily known as Critical Zone Experimental *Basin*.

## **3.1. Deficiencies in the design and operation of current experimental basins**

hydrometric measurement systems in natural basins still face challenges (e.g., the ability of current rain gauges with conventional deployment techniques to obtain precise precipitation inputs for natural basins remains controversial). It is necessary to have systematic experimental facilities consisting of laboratory physical models, hillslope "catchment", small-size natural watershed with uniformly dominant vegetation, small-size natural watershed with more complex vegetation, unchannelized catchments, nested sub-watersheds, and watersheds of large size. Such experimental facilities need to be situated in zones with different physio‐ graphic and hydroclimatic conditions. Obviously, this would require collaborative efforts at

304 Current Perspectives in Contaminant Hydrology and Water Resources Sustainability

**Figure 1.** Desertification of the Rocksheep River (Shiyang River) Basin [35], an endorheic basin with total drainage area of 40,690 km2. (a) prior to 121 BC, the "stage of natural harmony", the area of its terminal lake around 540 km2; (b) 420- 589 AD, the "stage of relying"; (c) ca 1900AD, the "stage of impact", the cultivated area during Qing Dynasty was developed quickly, which reached to an area about 4/5 of that of modern time. (d) ca 1990's, the "stage of dam‐

For the future of hydrologic basin study, another kind of experimental basin is suggested: the Critical Zone Experimental *Block* (CZEB) customarily known as Critical Zone Experimental

age", there are 4 dams since 1959, 22 dams until 1972 [35].

*Basin*.

**3. The Critical Zone Experimental Block (CZEB)**

an international scale.

*Watershed surface laterality.* Ignoring the inherent interconnectedness of surface and subsurface watershed components reflects the fact that watershed study is driven by practical rather than scientific requirements. Hydrologists still lack a comprehensive understanding of surface responses to nonvisible and difficult to measure subsurface constraints. A physiological analogy might be studying the skin while ignoring the functions beneath. Focusing on surface responses using surface monitoring will quickly trap the observer with misconceptions. For instance, hydrologists delineate watersheds with surface topography, though actual water‐ shed boundaries may cross topographically determined boundaries depending on subsurface features. Even our most fundamental measurements of the watershed – its shape and size – are strongly influenced by our inability to peer beneath the surface of the catchment.

*Hydrologic process laterality.* The hydrologic cycle is tightly coupled with other cycles [36]. The hydrologic cycle is also dynamic, not static, in nature. Studying the hydrologic cycle while ignoring the interconnectedness of other cycles is equivalent to studying blood circulation while only looking at blood vessels.

*Downward components laterality.* There is a tendency to emphasize measurements of downward components (e.g. rainfall, discharge and infiltration) rather than upward components (e.g. evaporation and transpiration aforementioned), even though upward components can be the dominant terms in water balance [37].

## **3.2. Lessons learned from the Chinese experimental basin studies over the last fifty years**

The Chinese experimental watershed studies experienced a saddle-backed fluctuation with two peaks [9] during their fifty-year history. The lessons learned were summarized in cooperation with J.J. McDonnell, C. Kendall and N.E. Peters, as following [9]:


It was also suggested that hydrological experimentation should not address "only to the hydrosphere but the interactions between atmosphere, lithosphere, biosphere and the intelligence-sphere." To this end, there is a need to "develop approaches using physical, chemical, isotopic, biogeochemical and hydrometeorological methods for basin studies as a 'hybrid' basin research"[38].

## **3.3. What is CZEB?**

The term "critical zone" has been used in many publications to refer to wide-ranging things in areas of geosciences and mineral, etc (e.g., the geological formation, the rhizosphere, the transitionalzones inalluvial coastalplainrivers, etc)[36].In2001,theNationalResearchCouncil of US recommended the integrated study of the "Critical Zone" (CZ) [36], defined as "the heterogeneous, near surface environment in which complex interactions involving rock, soil, water, air and living organisms regulate the natural habitat and determine availability of life sustaining resources" [16].TheCZis one ofthemost compellingresearchareas inEarthsciences in the 21st century [36]. This zone was further defined as extending from "the vegetation canopy to the zone of groundwater" in 2005 [39], or "the bedrock to the atmosphere boundary layer" [40] and, "top of the vegetation down to the bottom of the aquifers" [36] by 2010.

The Critical Zone Experimental Block (CZEB) geologically is a monolith-block with its surface, the watershed, bounded by topographical water divides, which define the surface boundary, and which lacks defined subsurface boundaries. A drainage basin conceptually is only its surface part, the visible face (Figure 2). The CZEB is actually an Experimental "Block" within the Critical Zone with a surface drainage basin (the watershed). It is a dynamic ecosystem coupled with various supporting systems but using hydrological processes as the unifying theme. It is a living, breathing, evolving boundary layer where rock, soil, water, air and living organisms interact [41]. The CZEB is the experimental hydrological watershed study in the CZ framework. It follows that the components within the CZEB (including monolith, hillslope, sub-watershed, etc) will be monitored according to the "downward and upward approaches to theory development in catchment hydrology" [14, 42]. "The Critical Zone Observatory network is the only type to integrate biological and geological sciences so tightly"[40], and following this, the CZEB network will perhaps be the only network to integrate the hydro‐ logical, biological and geological sciences so tightly as well.

Case III: There are stratified alluvia, potentially with multiple aquifers with aquicludes and aquitards, common in flat areas with thick deposit and deep bedrock, up to hundreds of meters or more. Groundwater recharge to the river can be separated into sensitive, active and passive zones. In this case, the bottom boundary of CZEB is defined as the bottom of the active recharge

Current Challenges in Experimental Watershed Hydrology

http://dx.doi.org/10.5772/55087

307

zone, the plane of TNL (Figure 3c).

**Figure 2.** A conceptualized Critical Zone Experimental Block (CZEB)

**Figure 3.** Boundaries of CZEB

#### **3.4. Boundaries of CZEB**

*Top.* The evaporation surface of the canopy in general is not the mean surface of the canopy but slightly lower. If the mean evaporation surface of the canopy is *h (m)* above the mean ground surface, then the top boundary of a CZEB is defined as *H(m)*, where *H*=(1.5 to 2.0)×*h* with coefficients of *h* varying according to the vegetation. This is mainly for the purpose of energy budget and eddy covariance flux observations. *H* is just the lower part of the atmos‐ phere boundary layer.

*Bottom.* There are three cases:

Case I: Bedrock is situated at a relatively shallow depth from the ground surface while the regolith is shallow, too. The bottom boundary is defined as the geological boundary (Figure 3a).

Case II: Bedrock is deep. The bottom boundary is defined as the plane where the tritium content of groundwater approaches zero or the detection limit of ±0.7 TU (the "tritium naught line" (TNL)[43]), which is same as that of Case III (Figure 3b).

**Figure 2.** A conceptualized Critical Zone Experimental Block (CZEB)

chemical, isotopic, biogeochemical and hydrometeorological methods for basin studies as a

306 Current Perspectives in Contaminant Hydrology and Water Resources Sustainability

The term "critical zone" has been used in many publications to refer to wide-ranging things in areas of geosciences and mineral, etc (e.g., the geological formation, the rhizosphere, the transitionalzones inalluvial coastalplainrivers, etc)[36].In2001,theNationalResearchCouncil of US recommended the integrated study of the "Critical Zone" (CZ) [36], defined as "the heterogeneous, near surface environment in which complex interactions involving rock, soil, water, air and living organisms regulate the natural habitat and determine availability of life sustaining resources" [16].TheCZis one ofthemost compellingresearchareas inEarthsciences in the 21st century [36]. This zone was further defined as extending from "the vegetation canopy to the zone of groundwater" in 2005 [39], or "the bedrock to the atmosphere boundary layer"

The Critical Zone Experimental Block (CZEB) geologically is a monolith-block with its surface, the watershed, bounded by topographical water divides, which define the surface boundary, and which lacks defined subsurface boundaries. A drainage basin conceptually is only its surface part, the visible face (Figure 2). The CZEB is actually an Experimental "Block" within the Critical Zone with a surface drainage basin (the watershed). It is a dynamic ecosystem coupled with various supporting systems but using hydrological processes as the unifying theme. It is a living, breathing, evolving boundary layer where rock, soil, water, air and living organisms interact [41]. The CZEB is the experimental hydrological watershed study in the CZ framework. It follows that the components within the CZEB (including monolith, hillslope, sub-watershed, etc) will be monitored according to the "downward and upward approaches to theory development in catchment hydrology" [14, 42]. "The Critical Zone Observatory network is the only type to integrate biological and geological sciences so tightly"[40], and following this, the CZEB network will perhaps be the only network to integrate the hydro‐

*Top.* The evaporation surface of the canopy in general is not the mean surface of the canopy but slightly lower. If the mean evaporation surface of the canopy is *h (m)* above the mean ground surface, then the top boundary of a CZEB is defined as *H(m)*, where *H*=(1.5 to 2.0)×*h* with coefficients of *h* varying according to the vegetation. This is mainly for the purpose of energy budget and eddy covariance flux observations. *H* is just the lower part of the atmos‐

Case I: Bedrock is situated at a relatively shallow depth from the ground surface while the regolith is shallow, too. The bottom boundary is defined as the geological boundary (Figure 3a). Case II: Bedrock is deep. The bottom boundary is defined as the plane where the tritium content of groundwater approaches zero or the detection limit of ±0.7 TU (the "tritium naught line"

[40] and, "top of the vegetation down to the bottom of the aquifers" [36] by 2010.

logical, biological and geological sciences so tightly as well.

(TNL)[43]), which is same as that of Case III (Figure 3b).

**3.4. Boundaries of CZEB**

phere boundary layer.

*Bottom.* There are three cases:

'hybrid' basin research"[38].

**3.3. What is CZEB?**

Case III: There are stratified alluvia, potentially with multiple aquifers with aquicludes and aquitards, common in flat areas with thick deposit and deep bedrock, up to hundreds of meters or more. Groundwater recharge to the river can be separated into sensitive, active and passive zones. In this case, the bottom boundary of CZEB is defined as the bottom of the active recharge zone, the plane of TNL (Figure 3c).

**Figure 3.** Boundaries of CZEB

#### *Lateral sides*

Part I: Above the ground surface up to height *H* as described above, delineated according to the surface topographic watershed boundary (Figure 3).

46]. Conservation of energy appears in the CZEB at all scales, with observable 'behaviors' of organization, symmetry, genes, etc. Lin suggested including "information" as one of the four general factors (i.e., conservation of energy and mass and, the accumulation of entropy and information, which dictates the evolutionary outcome and the functioning of the CZ ) [36].

Current Challenges in Experimental Watershed Hydrology

http://dx.doi.org/10.5772/55087

309

The CZEB can be separated into three interdependent zones as mentioned above with each zone consisting of the same base materials, the "*five-xing*" rock, soil, water, air and living organisms as well as non-material bases of energy, force and information. This coincides well with the '*Bagua*' ('*eight*-*gua'*) of the earliest Chinese philosophical work "*I –Jing"* ("The Book of Chang"). According to I-Jing, the origin of the universe was just a singularity of chaos ('*Taiji*'), consisting of '*Yang*' and '*Yin*', the unity of opposites, which produced '*eight-gua,*' the general principles governing development in the material world. Each "*gua*" is a combination of three whole lines ('*Yang*') and broken lines ('*Yin*'). As in the four general characteristics of CZ elaborated by Lin [36]; in the CZEB, the material base refers to the five '*gua*'s (the "*fivexing*") and the non-material base refers to another three '*gua*'s of energy, force and information, which together form a dynamic '*Bagua*'. Each zone has the same five material and three non-

**Figure 4.** Schematic 'reactors' of the CZEB. The dynamic *'Bagua*' in each 'reactor' operates with its own rate; 1- 5 : the material base (1- air, 2- water, 3- soil, 4- rock, 5- organisms); 6-8: the immaterial base (6- energy, 7- force, 8- informa‐

This is equally applicable to CZEB.

material components with different rates, Figure 4.

**3.6. The reactors of CZEB**

tion).

Part II: Below the ground surface, it is in general defined arbitrarily, except in the case of existing geological boundaries.

## **3.5. Functioning of CZEB**

CZEB is an ecological dynamic and evolving system. It also is a natural open system, exchang‐ ing mass and energy with its surrounding environment across 'arbitrary' boundaries. It is a dissipative, complex system, however, with some degree of self-organization [36].

*Interfaces.* CZEB encompasses "the near-surface biosphere and atmosphere, the entire pedo‐ sphere, and the surface and near-surface portion of the hydrosphere and lithosphere" [36]. Within the CZEB, various processes including hydrologic, atmospheric, lithospheric, geomor‐ phic and geochemical processes are coupled and dynamically interrelated. To simplify the organization of various interfacing processes throughout the CZEB, compartment zones can be broadly separated as follows:


This layering has a general trend of increasing density with depth, has a dampening effect on state variables with depth, and an increase in distance to energy input at the soil surface [36]. There is also "an overall trend of increasing response time" [36, 44].

*Mass – the material aspect of CZEB.* The material base of the CZEB is fixed and limited, involving rock, soil, water, air and living organisms [16]. On this point it happens to coincide with the ancient Chinese philosophy, the so-called "*Five-xing*", which holds that five fundamental elements form the universe and the Earth: "Jin"(metal), corresponding with the term "rock"; "*Mu*" (wood), with that of "living organisms"; "Shui"(water) with that of "water"; "Huo"(fire) with that of "air"; and "Tu"(soil).

*Energy and force – the driven aspect of CZEB.* Continuous energy fluxes and inherent forces drive the CZEB system. External solar energy is certainly the vital component of external energy source. Various inherent forces include: gravity, surface tension, intermolecular forces, capillary force, etc which control its dynamic situation.

*Organization and entropy – the philosophical aspect of CZEB.* This open system of structural dissipative processes and irreversible evolution tends to increase its entropy spontaneously and go towards disorder [36]. Even the exchange of entropy fluxes across its boundaries is continuous. The second law of thermodynamics, from which the concept of entropy was derived, is one of the fundamental natural laws. However, the role of feedbacks of this nonlinear system will promote "self-organization" as energy dissipates, providing opportu‐ nities for dissipating energy to act again within the system towards the direction of order. "Conservation without evolution is death. Evolution without conservation is madness" [45, 46]. Conservation of energy appears in the CZEB at all scales, with observable 'behaviors' of organization, symmetry, genes, etc. Lin suggested including "information" as one of the four general factors (i.e., conservation of energy and mass and, the accumulation of entropy and information, which dictates the evolutionary outcome and the functioning of the CZ ) [36]. This is equally applicable to CZEB.

#### **3.6. The reactors of CZEB**

*Lateral sides*

existing geological boundaries.

be broadly separated as follows:

with that of "air"; and "Tu"(soil).

**3.5. Functioning of CZEB**

Part I: Above the ground surface up to height *H* as described above, delineated according to

Part II: Below the ground surface, it is in general defined arbitrarily, except in the case of

CZEB is an ecological dynamic and evolving system. It also is a natural open system, exchang‐ ing mass and energy with its surrounding environment across 'arbitrary' boundaries. It is a

*Interfaces.* CZEB encompasses "the near-surface biosphere and atmosphere, the entire pedo‐ sphere, and the surface and near-surface portion of the hydrosphere and lithosphere" [36]. Within the CZEB, various processes including hydrologic, atmospheric, lithospheric, geomor‐ phic and geochemical processes are coupled and dynamically interrelated. To simplify the organization of various interfacing processes throughout the CZEB, compartment zones can

**2.** The unsaturated zone, "the belowground root zone and the deeper vadose zone "[36].

This layering has a general trend of increasing density with depth, has a dampening effect on state variables with depth, and an increase in distance to energy input at the soil surface [36].

*Mass – the material aspect of CZEB.* The material base of the CZEB is fixed and limited, involving rock, soil, water, air and living organisms [16]. On this point it happens to coincide with the ancient Chinese philosophy, the so-called "*Five-xing*", which holds that five fundamental elements form the universe and the Earth: "Jin"(metal), corresponding with the term "rock"; "*Mu*" (wood), with that of "living organisms"; "Shui"(water) with that of "water"; "Huo"(fire)

*Energy and force – the driven aspect of CZEB.* Continuous energy fluxes and inherent forces drive the CZEB system. External solar energy is certainly the vital component of external energy source. Various inherent forces include: gravity, surface tension, intermolecular forces,

*Organization and entropy – the philosophical aspect of CZEB.* This open system of structural dissipative processes and irreversible evolution tends to increase its entropy spontaneously and go towards disorder [36]. Even the exchange of entropy fluxes across its boundaries is continuous. The second law of thermodynamics, from which the concept of entropy was derived, is one of the fundamental natural laws. However, the role of feedbacks of this nonlinear system will promote "self-organization" as energy dissipates, providing opportu‐ nities for dissipating energy to act again within the system towards the direction of order. "Conservation without evolution is death. Evolution without conservation is madness" [45,

dissipative, complex system, however, with some degree of self-organization [36].

**1.** The zone above ground surface, the "aboveground vegetation zone" [36].

**3.** The saturated zone, the "saturated aquifer zone" [36].

capillary force, etc which control its dynamic situation.

There is also "an overall trend of increasing response time" [36, 44].

the surface topographic watershed boundary (Figure 3).

308 Current Perspectives in Contaminant Hydrology and Water Resources Sustainability

The CZEB can be separated into three interdependent zones as mentioned above with each zone consisting of the same base materials, the "*five-xing*" rock, soil, water, air and living organisms as well as non-material bases of energy, force and information. This coincides well with the '*Bagua*' ('*eight*-*gua'*) of the earliest Chinese philosophical work "*I –Jing"* ("The Book of Chang"). According to I-Jing, the origin of the universe was just a singularity of chaos ('*Taiji*'), consisting of '*Yang*' and '*Yin*', the unity of opposites, which produced '*eight-gua,*' the general principles governing development in the material world. Each "*gua*" is a combination of three whole lines ('*Yang*') and broken lines ('*Yin*'). As in the four general characteristics of CZ elaborated by Lin [36]; in the CZEB, the material base refers to the five '*gua*'s (the "*fivexing*") and the non-material base refers to another three '*gua*'s of energy, force and information, which together form a dynamic '*Bagua*'. Each zone has the same five material and three nonmaterial components with different rates, Figure 4.

**Figure 4.** Schematic 'reactors' of the CZEB. The dynamic *'Bagua*' in each 'reactor' operates with its own rate; 1- 5 : the material base (1- air, 2- water, 3- soil, 4- rock, 5- organisms); 6-8: the immaterial base (6- energy, 7- force, 8- informa‐ tion).

Each zone will behave as a 'feed-through reactor' according to Anderson et al [47]. It follows then that there are three feed-through reactors coupled together in the CZEB (Figure 4), with probably different rates and residence time of materials.

(e.g., "*isolation*") [9, 38]. Using detailed observation on ever-decreasing scales allows the investigator to focus and isolate processes from the naturally occurring heterogeneity (i.e., watershed, sub-watershed, unchanneled watershed, hillslope, monolith, a tree) (Figure 5). To some extent it is similar to dissection: splitting "problems into their smallest possible compo‐

Current Challenges in Experimental Watershed Hydrology

http://dx.doi.org/10.5772/55087

311

There are limitations in the field, however. Using a single tree as an example, the tree itself can be well instrumented, as well as the subsurface with pits up to meters deep, but still the fluxes from its bottom boundary would be unknown. In this sense, *isolation* refers to a system rather than the physical isolation of a natural unit, which can only be partially isolated in field. Strictly speaking, the field result from such natural systems can be explained only empirically, from a statistic or a stochastic view. Understanding causality mechanisms from observational studies alone is problematic. Thus, McDonnell et al (2007) noted that "most field experiments and observations in watershed science to date remain largely descriptive" [13]. They also discussed a shared 'philosophical path' which claims that "if we characterize enough hillslopes and watersheds around the world through detailed experimentations, some new understand‐ ing is bound to emerge eventually", the authors concluded that "what this approach to experimental design has succeeded in doing is to help characterize the idiosyncrasies of more and more watersheds, in different places and at different scales, but with little progress toward realizing the Dooge vision" of "hydrologic laws"[13]. This does not mean, however, that field studies have no place in hydrology. On the contrary, the problem becomes how to "put the pieces back together again" [48]. In our case it is how to put together hydrologic mechanisms

for a monolith, a slope, a sub-watershed and, a watershed.

**Figure 5.** Concepts of the two-way multi-scale approach: natural and artificial

*Control and synthesis in the artificial system.* A series of artificial systems are suggested and defined for "control and synthesis" [9, 38]. They encompass the pedon, polypedon, monolith,

nents"[48].

*Reactor I.* The first zone is the above ground surface zone, which also contains above ground vegetation up to its interface with atmosphere (Figure 4a). Evaporation and plant transpiration may account for 50% or more of the total local precipitation and use up to 50% of the total solar energy.

*Reactor II.* The second zone is the unsaturated zone, extending from the soil surface to the upper surface of the groundwater table. It consists of the whole soil profile: The O horizon (humus), A horizon (topsoil), B horizon (subsoil) and C, and potentially D, horizons (Figure 4a). The unsaturated soil zone has been recognized as "the most complicated biomaterials on the planet" [36].

*Reactor III.* The third zone is the saturated zone, extending from the groundwater table down to bedrock, including the capillary fringe (Figure 4a). There are two general cases for CZEB, including phreatic groundwater and confined aquifers.

## *Functioning of reactors*


## **4. A trial of the CZEB: The Chuzhou CZEB experimental system**

#### **4.1. The strategy: A two–way multi-scale approach**

*Isolation of the natural system.* In order to investigate the natural watershed system, which is filled with complexity and heterogeneity, simplification is incorporated as much as possible (e.g., "*isolation*") [9, 38]. Using detailed observation on ever-decreasing scales allows the investigator to focus and isolate processes from the naturally occurring heterogeneity (i.e., watershed, sub-watershed, unchanneled watershed, hillslope, monolith, a tree) (Figure 5). To some extent it is similar to dissection: splitting "problems into their smallest possible compo‐ nents"[48].

Each zone will behave as a 'feed-through reactor' according to Anderson et al [47]. It follows then that there are three feed-through reactors coupled together in the CZEB (Figure 4), with

*Reactor I.* The first zone is the above ground surface zone, which also contains above ground vegetation up to its interface with atmosphere (Figure 4a). Evaporation and plant transpiration may account for 50% or more of the total local precipitation and use up to 50% of the total solar

*Reactor II.* The second zone is the unsaturated zone, extending from the soil surface to the upper surface of the groundwater table. It consists of the whole soil profile: The O horizon (humus), A horizon (topsoil), B horizon (subsoil) and C, and potentially D, horizons (Figure 4a). The unsaturated soil zone has been recognized as "the most complicated biomaterials on the

*Reactor III.* The third zone is the saturated zone, extending from the groundwater table down to bedrock, including the capillary fringe (Figure 4a). There are two general cases for CZEB,

**1.** Reactor I includes hydrometeorological and ecohydrological processes, Reactor II includes hydropedological and hydroecological processes, while the Reactor III is more

**3.** Each Reactor has its own lateral flux exchanges *Lex* via the CZEB lateral boundaries (Figure 4b). The current hydrometric runoff data is only that from channel. Reactor I has vertical fluxes exchange *Vex* with atmosphere via the upper boundary of CZEB while Reactor III has *Vex* with deep aquifers and/or deep circulation via the bottom boundary of CZEB

**4.** Reactors are closely coupled within the boundaries of CZEB. There are flux exchanges *Wex* (Figure 5) between Reactor I and II and, and between II and III within the CZEB. There are flux channels *MC* (Figure 4) through I to III for materials, and flux channels *IMC*

**5.** The *MC* appears similar to blood vessel system and *IMC* to meridians and collaterals of the human body from the Chinese art of acupuncture. The three "reactors" are similar to the so-called human "*three Jiaos*" (Figure 4), a central concept of Chinese medicine

*Isolation of the natural system.* In order to investigate the natural watershed system, which is filled with complexity and heterogeneity, simplification is incorporated as much as possible

**2.** There is an overall trend of decreasing operation rates from Reactors I to III.

(Figure 4b). Fluxes include all material and immaterial components.

**4. A trial of the CZEB: The Chuzhou CZEB experimental system**

probably different rates and residence time of materials.

310 Current Perspectives in Contaminant Hydrology and Water Resources Sustainability

including phreatic groundwater and confined aquifers.

exclusively for hydrogeological processes.

(Figure 4) through I to III for non-materials.

**4.1. The strategy: A two–way multi-scale approach**

energy.

planet" [36].

*Functioning of reactors*

philosophy.

There are limitations in the field, however. Using a single tree as an example, the tree itself can be well instrumented, as well as the subsurface with pits up to meters deep, but still the fluxes from its bottom boundary would be unknown. In this sense, *isolation* refers to a system rather than the physical isolation of a natural unit, which can only be partially isolated in field. Strictly speaking, the field result from such natural systems can be explained only empirically, from a statistic or a stochastic view. Understanding causality mechanisms from observational studies alone is problematic. Thus, McDonnell et al (2007) noted that "most field experiments and observations in watershed science to date remain largely descriptive" [13]. They also discussed a shared 'philosophical path' which claims that "if we characterize enough hillslopes and watersheds around the world through detailed experimentations, some new understand‐ ing is bound to emerge eventually", the authors concluded that "what this approach to experimental design has succeeded in doing is to help characterize the idiosyncrasies of more and more watersheds, in different places and at different scales, but with little progress toward realizing the Dooge vision" of "hydrologic laws"[13]. This does not mean, however, that field studies have no place in hydrology. On the contrary, the problem becomes how to "put the pieces back together again" [48]. In our case it is how to put together hydrologic mechanisms for a monolith, a slope, a sub-watershed and, a watershed.

**Figure 5.** Concepts of the two-way multi-scale approach: natural and artificial

*Control and synthesis in the artificial system.* A series of artificial systems are suggested and defined for "control and synthesis" [9, 38]. They encompass the pedon, polypedon, monolith, hillslope, catchment, until subwatershed (Figure 5). Artificial systems are necessary in order to move beyond problems of boundary and parameter control found in natural systems. It is in fact inspired by the aforementioned empirical impasse of natural systems; artificial systems overcome heterogeneity by controlling all variables and leaving only test variables to fluctuate. **4.2. The trial practice of the CZEB: The Chuzhou CZEB experimental system**

34%, forest 5%, ponds 1.2%), 17.5 km2

monoliths (L1, L2) of 32 m2

*The natural CZEB, Nandadish*

(100%, 100%, 0% *idem*), 2.0 km2

watershed and the artificial catchment.

and bottom area (horizontal projection) of 490m2

km2

The Chuzhou CZEB can be attributed to the artistic conception of the *West Brook,* described in a poem of Tang Dynasty where intellectual exiles of Song Dynasty found their Arcadia ca 900 years ago. On the same brook in 1962, a hydrological experimental field station with three watersheds was founded. About fifteen years later, construction of the experimental watershed hydrology had barely survived from the cultural revolution of the 'Road to Perdition'. The Chuzhou CZEB was an adventurous plan for a two-way experimental watershed system including natural catchments and artificial catchments with controlled boundaries. The experiment was, unfortunately, interrupted again during the 'market-economy tides' of China and only partially completed. A turning point emerged a few years ago, with more efforts towards an innovative CZEB system. In the meantime, the components of the natural and artificial systems are: (1) the natural system with main infrastructure of hydrological moni‐ toring and isotope biomonitoring encompasses watersheds with drainage area of 82.1 km2 (hill

(56%, 10%, 0.6% *idem*), 4.5 km2

(5%, 2%, 30% *idem*), 0.06 km2

(Morning Glory Catchment, 100% weathered debris), and, a mini CZEB with directly observ‐

with measurable surface and subsurface responses encompasses a catchment with surficial

lysimeters. Here, the focus is on results from the two experimental watersheds - the mini CZEB

Nandadish has a surficial drainage area of 7897 m2 and rests on the consolidated bedrock of the concordant body of andesitic and tuffaceous facies with a thin weathered layer. The boundary condition of this mini CZEB belongs to the type shown in Figure 3a. The Quaternary regolith overlies the bedrock, its depth ranges from of 1 to 7 m with an average of 2.46 m; its bed rock topography was surveyed via 69 drillings (Figure 8). The regolith consists of bruni‐ solic soil of heavy loam, medium and clay loams; saprolite with prismatic and block structures. Horizontal and vertical fissures and cracks developed in the upper regolith. This watershed

each, a saturated monolith (LS) of 1 m2

able surface and subsurface responses with a drainage area of 7897m2

**Figure 7.** Reconciliation of two-way natural and artificial, *'Yin'* and *'Yang'* [14, 42].

(100%, 2%, 0% *idem*), 3.3

; (2) the artificial system

, and several weighing

, two grass

(80% paddy field), 4573 m2

, a crop monolith of 150m2

Current Challenges in Experimental Watershed Hydrology

http://dx.doi.org/10.5772/55087

313

*The laboratory system.* The two-way multi-scale approach can only be applied at the macroscopic level of watersheds or, maybe the mesoscopic level (e.g., catena monolith). In fact many problems (e.g., the old water paradox) can be trace back to mechanisms occurring on a microscopic level (e.g., the molecular and pore related mechanisms as that of microbial processes of some contaminants transformation). The laboratory system connects with and supplements the artificial system (Figure 5).

*Could this be a new vision for watershed hydrology and a unified theory of hydrology?* Lin summarized three systems of different levels of complexity and a vast gap between the two extremes [36] as shown in Figure 6. The coupled natural and artificial systems are analogous to '*Yin*' and '*Yang*', which is the successful philosophy of '*Taiji*' (Figure 7). Klmes suggested "a rational search for meaningful conceptualization in hydrology can proceed along two routes: upwards and downwards "[42], or, the "upward or bottom-up approach" and the "downward or topdown approach" [14]. The two-way multi-scale approach perhaps could also be upwards and downwards. Hopefully, it could be a way for a new vision for watershed hydrology [13] and, for the unified theory of hydrology [14].

**Figure 6.** The gap between extremes of complex catchment systems. (a) Three types of complex catchment systems related to different treatment methods [14,19,36]; (b) The supposed field and artificial two-way multi-scale frame‐ work and the gap depending on the level of control. The gap becomes the bridge.

## **4.2. The trial practice of the CZEB: The Chuzhou CZEB experimental system**

The Chuzhou CZEB can be attributed to the artistic conception of the *West Brook,* described in a poem of Tang Dynasty where intellectual exiles of Song Dynasty found their Arcadia ca 900 years ago. On the same brook in 1962, a hydrological experimental field station with three watersheds was founded. About fifteen years later, construction of the experimental watershed hydrology had barely survived from the cultural revolution of the 'Road to Perdition'. The Chuzhou CZEB was an adventurous plan for a two-way experimental watershed system including natural catchments and artificial catchments with controlled boundaries. The experiment was, unfortunately, interrupted again during the 'market-economy tides' of China and only partially completed. A turning point emerged a few years ago, with more efforts towards an innovative CZEB system. In the meantime, the components of the natural and artificial systems are: (1) the natural system with main infrastructure of hydrological moni‐ toring and isotope biomonitoring encompasses watersheds with drainage area of 82.1 km2 (hill 34%, forest 5%, ponds 1.2%), 17.5 km2 (56%, 10%, 0.6% *idem*), 4.5 km2 (100%, 2%, 0% *idem*), 3.3 km2 (100%, 100%, 0% *idem*), 2.0 km2 (5%, 2%, 30% *idem*), 0.06 km2 (80% paddy field), 4573 m2 (Morning Glory Catchment, 100% weathered debris), and, a mini CZEB with directly observ‐ able surface and subsurface responses with a drainage area of 7897m2 ; (2) the artificial system with measurable surface and subsurface responses encompasses a catchment with surficial and bottom area (horizontal projection) of 490m2 , a crop monolith of 150m2 , two grass monoliths (L1, L2) of 32 m2 each, a saturated monolith (LS) of 1 m2 , and several weighing lysimeters. Here, the focus is on results from the two experimental watersheds - the mini CZEB watershed and the artificial catchment.

**Figure 7.** Reconciliation of two-way natural and artificial, *'Yin'* and *'Yang'* [14, 42].

## *The natural CZEB, Nandadish*

hillslope, catchment, until subwatershed (Figure 5). Artificial systems are necessary in order to move beyond problems of boundary and parameter control found in natural systems. It is in fact inspired by the aforementioned empirical impasse of natural systems; artificial systems overcome heterogeneity by controlling all variables and leaving only test variables to fluctuate.

*The laboratory system.* The two-way multi-scale approach can only be applied at the macroscopic level of watersheds or, maybe the mesoscopic level (e.g., catena monolith). In fact many problems (e.g., the old water paradox) can be trace back to mechanisms occurring on a microscopic level (e.g., the molecular and pore related mechanisms as that of microbial processes of some contaminants transformation). The laboratory system connects with and

*Could this be a new vision for watershed hydrology and a unified theory of hydrology?* Lin summarized three systems of different levels of complexity and a vast gap between the two extremes [36] as shown in Figure 6. The coupled natural and artificial systems are analogous to '*Yin*' and '*Yang*', which is the successful philosophy of '*Taiji*' (Figure 7). Klmes suggested "a rational search for meaningful conceptualization in hydrology can proceed along two routes: upwards and downwards "[42], or, the "upward or bottom-up approach" and the "downward or topdown approach" [14]. The two-way multi-scale approach perhaps could also be upwards and downwards. Hopefully, it could be a way for a new vision for watershed hydrology [13] and,

**Figure 6.** The gap between extremes of complex catchment systems. (a) Three types of complex catchment systems related to different treatment methods [14,19,36]; (b) The supposed field and artificial two-way multi-scale frame‐

work and the gap depending on the level of control. The gap becomes the bridge.

supplements the artificial system (Figure 5).

312 Current Perspectives in Contaminant Hydrology and Water Resources Sustainability

for the unified theory of hydrology [14].

Nandadish has a surficial drainage area of 7897 m2 and rests on the consolidated bedrock of the concordant body of andesitic and tuffaceous facies with a thin weathered layer. The boundary condition of this mini CZEB belongs to the type shown in Figure 3a. The Quaternary regolith overlies the bedrock, its depth ranges from of 1 to 7 m with an average of 2.46 m; its bed rock topography was surveyed via 69 drillings (Figure 8). The regolith consists of bruni‐ solic soil of heavy loam, medium and clay loams; saprolite with prismatic and block structures. Horizontal and vertical fissures and cracks developed in the upper regolith. This watershed has an altitude difference of 12.9 m with a surface slope of 6.7% to 17.1% at different directions. The brook gradient is 6.7%. The coverage during the watershed's construction in 1979 was natural grasses with small shrubs and a few Masson pines aged 5 to 6 years. Since then, coverage has shifted to a dense forest with canopy height ca 8 m.

*The artificial catchment, Hydrohill*

surveying) - 512 m2

at the weirs.

The Hydrohill catchment with a drainage area of 490 m2

hillslope was first excavated to bedrock to create a bare catchment of ca 4700 m2

aquiclude was created above it and consisted of two intersecting slopes dipping towards each other at 10°with an overall downslope gradients of 16.9°(Figure 10a). Impermeable concrete walls enclose the catchment on all sides to prevent any flow of water across the catchment boundaries. Silt-loam soil was removed, layer by layer, from an agricultural field near the Hydrohill site and installed on the artificial aquiclude, layer-by-layer again. Also the bulk density of the soil was adjusted during filling to approximate the natural soil profile of the original agricultural field. Hence, the final 1-m soil 'profile' was identical, at least in compo‐ sition, to the natural soil profile. Grass was then planted over the surface and soil allowed to settle for 3 years (Figure 10b), after which time a central drainage trench was then excavated at the intersection of the two slopes and the water-sampling instrumentation was installed.

Five fiberglass troughs, each 40 cm wide and 40 m in length, were installed longitudinally in the trench. These troughs were stacked on top of each other to create a set of long zero-tension lysimeters (Figure 11a). Each trough has a 20 cm aluminium lip that extends horizontally into the soil layer to prevent leakage between layers. Water collected in troughs is routed through V-notch based logarithm weirs located in a gauging room (Figure 11b) under the hill where discharge is continuously monitored. Water samples are collected manually above the ponding

As illustrated in Figure 11a, the uppermost trough collects rain; the next lower trough collects surface runoff. The next three troughs collect subsurface flow from soil layers spanning the depths of 0-30, 30-60, and 60-100 cm. These troughs will be referred to as the 30 cm, 60 cm, and 100 cm troughs. The source of the water in these troughs (i.e., whether the water is derived from interflow or saturated flow) varies locally and during storms. The lowermost trough collects either saturated flow or interflow, depending on the height of the water table. When the water table is high, saturated flow may be collected in both of the lower two troughs.

A network of 21 aluminium alloy access tubes for neutron moisture gauges [51], a tensiometer scanner, and 22 wells for water-table measurements and groundwater sampling were installed (Figure 12). The wells were drilled to the aquiclude and are slotted along the lowermost 20 cm. After installation, the spaces around the slotted lengths were packed with sands to allow movement of groundwater to the well, space above the slotted lengths were packed with clay to prevent vertical drainage along the pipes (Figure 11c). The neutron probe access tubes (Figure 11c) and soil water potential tensiometers were positioned adjacent to the wells for soil water monitoring. Catchment evaporation was monitored by methods of water balance, energy budget and, soil water variations above zero flux potential. The plan view of the surface topography, the locations of wells and access tubes, and, the central stacked lysimeter troughs

and that of energy budget set are showing in Figure 12.

(horizontal projection by plane

http://dx.doi.org/10.5772/55087

Current Challenges in Experimental Watershed Hydrology

. A concrete

315

for its inclined surface [49, 50] - is sited on a small andesitic hill. The entire

The main infrastructure includes: (1) measurement for energy flux, water flux, geochemical flux, gas flux; (2) changes of watershed storages of energy, water, gas and that of geochemical ions; (3) CZ tree experiment; and (4) variations of precipitation isotope fingerprints in flux and in storage. The surface and subsurface runoff are monitored directly via a trench, which extends upslope to capture subsurface and surface flow. This CZEB block has deep soils near the divide but only 1-meter depths near the outlet, making the block easy to close via a concrete wall installed to the bedrock at the outlet point. In this way, all surface and subsurface flow drain into discharge measuring structures (Figure 9).

**Figure 8.** The mini CZEB Nandadish. (a) Its surficial watershed; (b) Its bedrock; (c) The depths of deposit above its bed‐ rock.

**Figure 9.** Discharge measurement structures for different runoff components from troughs with its location showing in Figure 8(a). 1- for rainfall; 2 and 3- for surface runoff; 4 and 5 – for interflow (30 cm below the soil surface); 6 – for interflow and groundwater flow (down to bedrock); 7 and 8 for total runoff (weirs are not shown). 1,3,5,6 and 8 are V-notch sharp-crested weir; 2, 4 and 7 are the full width rectangular sharp-crested weir; 9 is the tracking water head gauge. The picture design is the original concept, through renovations are now underway for gauge modernization.

## *The artificial catchment, Hydrohill*

has an altitude difference of 12.9 m with a surface slope of 6.7% to 17.1% at different directions. The brook gradient is 6.7%. The coverage during the watershed's construction in 1979 was natural grasses with small shrubs and a few Masson pines aged 5 to 6 years. Since then,

The main infrastructure includes: (1) measurement for energy flux, water flux, geochemical flux, gas flux; (2) changes of watershed storages of energy, water, gas and that of geochemical ions; (3) CZ tree experiment; and (4) variations of precipitation isotope fingerprints in flux and in storage. The surface and subsurface runoff are monitored directly via a trench, which extends upslope to capture subsurface and surface flow. This CZEB block has deep soils near the divide but only 1-meter depths near the outlet, making the block easy to close via a concrete wall installed to the bedrock at the outlet point. In this way, all surface and subsurface flow

**Figure 8.** The mini CZEB Nandadish. (a) Its surficial watershed; (b) Its bedrock; (c) The depths of deposit above its bed‐

**Figure 9.** Discharge measurement structures for different runoff components from troughs with its location showing in Figure 8(a). 1- for rainfall; 2 and 3- for surface runoff; 4 and 5 – for interflow (30 cm below the soil surface); 6 – for interflow and groundwater flow (down to bedrock); 7 and 8 for total runoff (weirs are not shown). 1,3,5,6 and 8 are V-notch sharp-crested weir; 2, 4 and 7 are the full width rectangular sharp-crested weir; 9 is the tracking water head gauge. The picture design is the original concept, through renovations are now underway for gauge modernization.

coverage has shifted to a dense forest with canopy height ca 8 m.

314 Current Perspectives in Contaminant Hydrology and Water Resources Sustainability

drain into discharge measuring structures (Figure 9).

rock.

The Hydrohill catchment with a drainage area of 490 m2 (horizontal projection by plane surveying) - 512 m2 for its inclined surface [49, 50] - is sited on a small andesitic hill. The entire hillslope was first excavated to bedrock to create a bare catchment of ca 4700 m2 . A concrete aquiclude was created above it and consisted of two intersecting slopes dipping towards each other at 10°with an overall downslope gradients of 16.9°(Figure 10a). Impermeable concrete walls enclose the catchment on all sides to prevent any flow of water across the catchment boundaries. Silt-loam soil was removed, layer by layer, from an agricultural field near the Hydrohill site and installed on the artificial aquiclude, layer-by-layer again. Also the bulk density of the soil was adjusted during filling to approximate the natural soil profile of the original agricultural field. Hence, the final 1-m soil 'profile' was identical, at least in compo‐ sition, to the natural soil profile. Grass was then planted over the surface and soil allowed to settle for 3 years (Figure 10b), after which time a central drainage trench was then excavated at the intersection of the two slopes and the water-sampling instrumentation was installed.

Five fiberglass troughs, each 40 cm wide and 40 m in length, were installed longitudinally in the trench. These troughs were stacked on top of each other to create a set of long zero-tension lysimeters (Figure 11a). Each trough has a 20 cm aluminium lip that extends horizontally into the soil layer to prevent leakage between layers. Water collected in troughs is routed through V-notch based logarithm weirs located in a gauging room (Figure 11b) under the hill where discharge is continuously monitored. Water samples are collected manually above the ponding at the weirs.

As illustrated in Figure 11a, the uppermost trough collects rain; the next lower trough collects surface runoff. The next three troughs collect subsurface flow from soil layers spanning the depths of 0-30, 30-60, and 60-100 cm. These troughs will be referred to as the 30 cm, 60 cm, and 100 cm troughs. The source of the water in these troughs (i.e., whether the water is derived from interflow or saturated flow) varies locally and during storms. The lowermost trough collects either saturated flow or interflow, depending on the height of the water table. When the water table is high, saturated flow may be collected in both of the lower two troughs.

A network of 21 aluminium alloy access tubes for neutron moisture gauges [51], a tensiometer scanner, and 22 wells for water-table measurements and groundwater sampling were installed (Figure 12). The wells were drilled to the aquiclude and are slotted along the lowermost 20 cm. After installation, the spaces around the slotted lengths were packed with sands to allow movement of groundwater to the well, space above the slotted lengths were packed with clay to prevent vertical drainage along the pipes (Figure 11c). The neutron probe access tubes (Figure 11c) and soil water potential tensiometers were positioned adjacent to the wells for soil water monitoring. Catchment evaporation was monitored by methods of water balance, energy budget and, soil water variations above zero flux potential. The plan view of the surface topography, the locations of wells and access tubes, and, the central stacked lysimeter troughs and that of energy budget set are showing in Figure 12.

**5. Selected findings**

**5.1. Runoff composition**

der renovation now.

than 1 km2

between components during an event.

The catchment response to precipitation consists of surface and subsurface components. These components were directly measured at Hydrohill with both isotopic and hydrochemical evidences (Figure 13) [53, 54]. Surface runoff was not always the dominant component during a rainfall-runoff event. The hydrograph composition of surface and subsurface flows was discriminated into four broad categories (Figure 14): 1) S type which is dominated by surface runoff; 2) SS type which is dominated by subsurface flows; 3) Intermediate type M, with largely equal contributions of surface and subsurface flows; and 4) Variation type V, with switching

Current Challenges in Experimental Watershed Hydrology

http://dx.doi.org/10.5772/55087

317

**Figure 12.** Hydrometric monitoring of Hydrohill. (a) Plan view of the surface topograph showing the locations of groundwater wells, access tubes for neutron moisture probe with potential scanner, the central stacked lysimeter troughs, the energy budget set. (b) Whole view of its surface watershed in operation during early years (1980's), un‐

The physical based distributed-parameter modeling appears to be a promising direction for watershed hydrology. However the scale of model elements, and their parameters, is a lingering problem. This has led to investigations into the establishment of a 'representative elementary area' (REP), 1 km2 [55]. In fact, because of the spatiotemporal heterogeneity of hydrological parameters, little is known regarding hydrologic parameterization in scales less

. Previously unimagined heterogeneities were observed at Hydrohill:

*Intrastorm isotopic heterogeneity of shallow groundwater.* During a storm in July 1989, the newly developed groundwater showed considerable spatial and temporal variability in theδ18O (Figure 15a), with values ranging from -12‰ to -6‰ while the groundwater derived from preevent soil water ranged from 0 to 100% at different times and locations in this catchment [56].

**5.2. Spatial and temporal variability in hydrological parameters**

**Figure 10.** The artificial catchment Hydrohill: (a) Its concrete aquiclude; (b) Three years waiting for the settling of filled soil to try to match close to the natural soil profile.

**Figure 11.** The artificial Hydrohill: (a) Schematic cross-section of rain, surface, and subsurface flow collectors at the catchment; (b) The V-notch based logarithm weirs for discharge measurement for these components; (c) The con‐ structions of wells for groundwater monitoring and access tubes for neutron moisture gauge. (b) and (c) are the origi‐ nal design, now under renovation;

The renovation phase of the main infrastructures for this physically modeled CZEB includes: (1) measurement for energy flux, water flux, geochemical flux, gas flux; (2) changes of watershed storage of energy, water, gas and that of geochemical ions; (3) a small separate laboratory within the gauging room for high-frequency, real time measurement of isotope concentrations in all hydrologic components using of a laser spectrometer LWIA [52]; (4) a removable rainfall simulator on the tracks capable of covering an area of more than 500 m2 , capable of simulating a range of rainfall intensities; and (5) a system for carbon balance of this block.

## **5. Selected findings**

## **5.1. Runoff composition**

**Figure 10.** The artificial catchment Hydrohill: (a) Its concrete aquiclude; (b) Three years waiting for the settling of filled

**Figure 11.** The artificial Hydrohill: (a) Schematic cross-section of rain, surface, and subsurface flow collectors at the catchment; (b) The V-notch based logarithm weirs for discharge measurement for these components; (c) The con‐ structions of wells for groundwater monitoring and access tubes for neutron moisture gauge. (b) and (c) are the origi‐

The renovation phase of the main infrastructures for this physically modeled CZEB includes: (1) measurement for energy flux, water flux, geochemical flux, gas flux; (2) changes of watershed storage of energy, water, gas and that of geochemical ions; (3) a small separate laboratory within the gauging room for high-frequency, real time measurement of isotope concentrations in all hydrologic components using of a laser spectrometer LWIA [52]; (4) a removable rainfall simulator on the tracks capable of covering an area of more than 500 m2

capable of simulating a range of rainfall intensities; and (5) a system for carbon balance of this

,

soil to try to match close to the natural soil profile.

316 Current Perspectives in Contaminant Hydrology and Water Resources Sustainability

nal design, now under renovation;

block.

The catchment response to precipitation consists of surface and subsurface components. These components were directly measured at Hydrohill with both isotopic and hydrochemical evidences (Figure 13) [53, 54]. Surface runoff was not always the dominant component during a rainfall-runoff event. The hydrograph composition of surface and subsurface flows was discriminated into four broad categories (Figure 14): 1) S type which is dominated by surface runoff; 2) SS type which is dominated by subsurface flows; 3) Intermediate type M, with largely equal contributions of surface and subsurface flows; and 4) Variation type V, with switching between components during an event.

**Figure 12.** Hydrometric monitoring of Hydrohill. (a) Plan view of the surface topograph showing the locations of groundwater wells, access tubes for neutron moisture probe with potential scanner, the central stacked lysimeter troughs, the energy budget set. (b) Whole view of its surface watershed in operation during early years (1980's), un‐ der renovation now.

## **5.2. Spatial and temporal variability in hydrological parameters**

The physical based distributed-parameter modeling appears to be a promising direction for watershed hydrology. However the scale of model elements, and their parameters, is a lingering problem. This has led to investigations into the establishment of a 'representative elementary area' (REP), 1 km2 [55]. In fact, because of the spatiotemporal heterogeneity of hydrological parameters, little is known regarding hydrologic parameterization in scales less than 1 km2 . Previously unimagined heterogeneities were observed at Hydrohill:

*Intrastorm isotopic heterogeneity of shallow groundwater.* During a storm in July 1989, the newly developed groundwater showed considerable spatial and temporal variability in theδ18O (Figure 15a), with values ranging from -12‰ to -6‰ while the groundwater derived from preevent soil water ranged from 0 to 100% at different times and locations in this catchment [56]. The groundwater samples for δ18O were collected at three times during the storm from the wells distributed in whole area as shown in Figure 15a. At the end of the storm, the ground‐ waterδ18O values showed a 4‰ range in composition. This range of compositions during the storm appears to be caused by the combination effects of intrastorm variation in rainδ18O (Figure 13) and spatial and temporal variability in subsurface flowpaths [56] – including the downward displacement of pre-storm water versus delivery of new water to the aquiclude via macropores [54].

**Figure 13.** Isotopic and hydrochemical evidences of various runoff components. Samples for isotopes were analysed in Menlo Park laboratory of USGS by C. Kendall (same hereinafter)

**Figure 15.** (a) Contour diagrams of δ18O in groundwater at three sampling times and (b) Typical spatial distribution of

Current Challenges in Experimental Watershed Hydrology

http://dx.doi.org/10.5772/55087

319

About ten years ago, Kirchner described a double paradox in catchment hydrology and geochemistry [31]. His objective was to raise questions rather than provide answers and a worthy challenge was then put forward (i.e., "need to take up the search for a unified theory", so that the paradoxes "would no longer seem paradoxes") [31]. For this challenge, perhaps catchments with directly observable surface and subsurface components could be useful in two ways: First, to identify how this double paradox will emerge in the subsurface and surface components individually, instead of in the integrated form at the stream outlet; and secondly,

*The delayed correspondency phenomena.* Three isolated rainfall events are presented from May 1987 at Hydrohill, and are ideal for establishing rainfall-runoff relationships as well as the unit hydrograph (Figure 16a). In short, the unit hydrography conceptualization is that rainfall event A produces the hydrograph peak A, rainfall events B and C produce the peaks B and C etc. It can be named as the concept of *'one-to-one correspondency'*. However, using isotope tracers, rainfall-runoff events at Hydrohill showed that rainfall event A corresponded to the hydrog‐ raphy peaks A1 and A2, though A2 only emerged later during event B. This pattern continued throughout subsequent storms with event B corresponding to hydrograph peaks B1 and B2, though B2 only emerged later during event C. This was termed as the phenomena of *delayed correspondency* [57]. In the case of Hydrohill, this delay appears very distinct. It is presumed that such a *delayed correspondency* will tend to diminish as drainage area increases. The diminishment of the *delayed correspondency* phenomena has been verified in a natural water‐

to trace surface and subsurface generation mechanisms using isotopic fingerprints.

the daily land evaporation rate.

**5.3. The double paradox**

**Figure 14.** Types of runoff composition summarized from observed data of Hydrohill during 1981- 1992.

*Annual heterogeneity of land evaporation rate* Land evaporation was measured via the variations of unsaturated soil water at different positions in the soil profile during its matrix potential distribution in that profile showing a zero flux plane. Intensive measurements of volumetric soil moisture by neutron probe at 4 points within each of the 21 access tubes (84 in total) were made during the wet and dry seasons. Representative examples of heterogeneity in the daily evaporation rate are shown in Figure 15b [51].

**Figure 15.** (a) Contour diagrams of δ18O in groundwater at three sampling times and (b) Typical spatial distribution of the daily land evaporation rate.

## **5.3. The double paradox**

The groundwater samples for δ18O were collected at three times during the storm from the wells distributed in whole area as shown in Figure 15a. At the end of the storm, the ground‐ waterδ18O values showed a 4‰ range in composition. This range of compositions during the storm appears to be caused by the combination effects of intrastorm variation in rainδ18O (Figure 13) and spatial and temporal variability in subsurface flowpaths [56] – including the downward displacement of pre-storm water versus delivery of new water to the aquiclude via

318 Current Perspectives in Contaminant Hydrology and Water Resources Sustainability

**Figure 13.** Isotopic and hydrochemical evidences of various runoff components. Samples for isotopes were analysed

**Figure 14.** Types of runoff composition summarized from observed data of Hydrohill during 1981- 1992.

*Annual heterogeneity of land evaporation rate* Land evaporation was measured via the variations of unsaturated soil water at different positions in the soil profile during its matrix potential distribution in that profile showing a zero flux plane. Intensive measurements of volumetric soil moisture by neutron probe at 4 points within each of the 21 access tubes (84 in total) were made during the wet and dry seasons. Representative examples of heterogeneity in the daily

in Menlo Park laboratory of USGS by C. Kendall (same hereinafter)

evaporation rate are shown in Figure 15b [51].

macropores [54].

About ten years ago, Kirchner described a double paradox in catchment hydrology and geochemistry [31]. His objective was to raise questions rather than provide answers and a worthy challenge was then put forward (i.e., "need to take up the search for a unified theory", so that the paradoxes "would no longer seem paradoxes") [31]. For this challenge, perhaps catchments with directly observable surface and subsurface components could be useful in two ways: First, to identify how this double paradox will emerge in the subsurface and surface components individually, instead of in the integrated form at the stream outlet; and secondly, to trace surface and subsurface generation mechanisms using isotopic fingerprints.

*The delayed correspondency phenomena.* Three isolated rainfall events are presented from May 1987 at Hydrohill, and are ideal for establishing rainfall-runoff relationships as well as the unit hydrograph (Figure 16a). In short, the unit hydrography conceptualization is that rainfall event A produces the hydrograph peak A, rainfall events B and C produce the peaks B and C etc. It can be named as the concept of *'one-to-one correspondency'*. However, using isotope tracers, rainfall-runoff events at Hydrohill showed that rainfall event A corresponded to the hydrog‐ raphy peaks A1 and A2, though A2 only emerged later during event B. This pattern continued throughout subsequent storms with event B corresponding to hydrograph peaks B1 and B2, though B2 only emerged later during event C. This was termed as the phenomena of *delayed correspondency* [57]. In the case of Hydrohill, this delay appears very distinct. It is presumed that such a *delayed correspondency* will tend to diminish as drainage area increases. The diminishment of the *delayed correspondency* phenomena has been verified in a natural water‐ shed with drainage area of 82.1 km2 (Figure 16b). Rainfall event A produces a runoff hydro‐ graph peak A, some of the constituents of rainfall event A were delayed to emerge during the rainfall event B. It is worth noting that the hydrograph produced by rainfall event A does not have the bell shape referred to in the unit hydrograph concept. In fact this *delayed correspond‐ ency* shows the formation of pre-event water ("old" water) during the runoff process. It follows that the current *'one-to-one correspondency',* used to conceptualize rainfall-runoff relationships in applied hydrology, will be associated with large uncertainty.

**Figure 17.** Processes of various runoff components, and that of their pre-event water of Nandadish. (a) for S type; (b)

Current Challenges in Experimental Watershed Hydrology

http://dx.doi.org/10.5772/55087

321

**Figure 18.** The Mg2+ and Cl- processes in event rainfall and runoff components of Hydrohill. Samples were analysed in

**5.4. Does diel signal of hydrochemical constituents emerge linkage among multi–processes?**

For the better understanding of the links between contaminants and multi-processes, explo‐ ration of diel signals in natural waters may yield "insight into the intricate linkages among hydrological, biological, and geochemical processes [23]". Diel variations of various hydro‐ logical constituents in surface and subsurface runoff responses during rainfall events were monitored in artificial catchments and monoliths with examples as following (Figure 19).

*Variations of pH.* The diel variations of pH in SR, IF, and GF of Hydrohill show their own individuations (solid lines in Figure 19a). Diel variation curve of SR to a large extent appears similar to that of rainfall. However, pH variation curve of IF is contrary to that of rainfall after 10 a.m., while the GF curve appears as the flattened IF curve. The pH variation curve of the interflow of monolith L1 and that of the flow of monolith LS are reasonably similar to that of

Atlanta laboratory of USGS by N.E. Peters (same hereinafter)

for SS type

**Figure 16.** The delayed correspondency phenomena.

*"Old" water paradox identified from both surface and subsurface flows.* Data from the 'mini CZEB' watershed, Nandadish, suggest that pre-event water ("old water") appeared in all runoff components including surface runoff, interflow and groundwater or saturated flow (Figure 17). Figure 17a and Figure 17b refer to the surface runoff dominated type (type S, from above), and subsurface flow dominated type (type SS) respectively. The surface runoff and subsurface runoff processes are shown separately with their corresponding proportions of pre-event water. For the type S, the pre-event water accounts for 9% and 24% of the total amount of surface runoff and of subsurface flow respectively while for the type SS, it becomes 11%and 89% respectively [58]. This reveals that even in a catchment with an average soil depth of 2.46 m, large volumes of pre-event water ("old" water) are stored and released promptly by event input.

*Hydrochemistry distribution paradox* This was termed by Kirchner as the "'variable chemistry of old water' paradox: although baseflow and stormflow are both composed mostly of 'old' water, they often have very different chemical signatures" [31]. However, the artificial catchment at Hydrohill shows paradoxically the distribution of hydrochemical compositions in different runoff components, including event rainfall. This '*hydrochemistry distribution paradox*' results largely from: (1) inorganic ions in event rainfall input emerge in all runoff components; (2) a strong similarity between rainfall and surface runoff but less similarity in subsurface components; and (3) the fact that the total amount of ions of event rainfall is sometimes much smaller than the sum of all runoff components [59]. Figure 18 shows the processes of Mg2+ and Cl in event rainfall, surface runoff, interflow and groundwater flow (saturated flow).

shed with drainage area of 82.1 km2

**Figure 16.** The delayed correspondency phenomena.

input.

processes of Mg2+ and Cl-

(saturated flow).

in applied hydrology, will be associated with large uncertainty.

320 Current Perspectives in Contaminant Hydrology and Water Resources Sustainability

(Figure 16b). Rainfall event A produces a runoff hydro‐

graph peak A, some of the constituents of rainfall event A were delayed to emerge during the rainfall event B. It is worth noting that the hydrograph produced by rainfall event A does not have the bell shape referred to in the unit hydrograph concept. In fact this *delayed correspond‐ ency* shows the formation of pre-event water ("old" water) during the runoff process. It follows that the current *'one-to-one correspondency',* used to conceptualize rainfall-runoff relationships

*"Old" water paradox identified from both surface and subsurface flows.* Data from the 'mini CZEB' watershed, Nandadish, suggest that pre-event water ("old water") appeared in all runoff components including surface runoff, interflow and groundwater or saturated flow (Figure 17). Figure 17a and Figure 17b refer to the surface runoff dominated type (type S, from above), and subsurface flow dominated type (type SS) respectively. The surface runoff and subsurface runoff processes are shown separately with their corresponding proportions of pre-event water. For the type S, the pre-event water accounts for 9% and 24% of the total amount of surface runoff and of subsurface flow respectively while for the type SS, it becomes 11%and 89% respectively [58]. This reveals that even in a catchment with an average soil depth of 2.46 m, large volumes of pre-event water ("old" water) are stored and released promptly by event

*Hydrochemistry distribution paradox* This was termed by Kirchner as the "'variable chemistry of old water' paradox: although baseflow and stormflow are both composed mostly of 'old' water, they often have very different chemical signatures" [31]. However, the artificial catchment at Hydrohill shows paradoxically the distribution of hydrochemical compositions in different runoff components, including event rainfall. This '*hydrochemistry distribution paradox*' results largely from: (1) inorganic ions in event rainfall input emerge in all runoff components; (2) a strong similarity between rainfall and surface runoff but less similarity in subsurface components; and (3) the fact that the total amount of ions of event rainfall is sometimes much smaller than the sum of all runoff components [59]. Figure 18 shows the

in event rainfall, surface runoff, interflow and groundwater flow

**Figure 17.** Processes of various runoff components, and that of their pre-event water of Nandadish. (a) for S type; (b) for SS type

**Figure 18.** The Mg2+ and Cl- processes in event rainfall and runoff components of Hydrohill. Samples were analysed in Atlanta laboratory of USGS by N.E. Peters (same hereinafter)

#### **5.4. Does diel signal of hydrochemical constituents emerge linkage among multi–processes?**

For the better understanding of the links between contaminants and multi-processes, explo‐ ration of diel signals in natural waters may yield "insight into the intricate linkages among hydrological, biological, and geochemical processes [23]". Diel variations of various hydro‐ logical constituents in surface and subsurface runoff responses during rainfall events were monitored in artificial catchments and monoliths with examples as following (Figure 19).

*Variations of pH.* The diel variations of pH in SR, IF, and GF of Hydrohill show their own individuations (solid lines in Figure 19a). Diel variation curve of SR to a large extent appears similar to that of rainfall. However, pH variation curve of IF is contrary to that of rainfall after 10 a.m., while the GF curve appears as the flattened IF curve. The pH variation curve of the interflow of monolith L1 and that of the flow of monolith LS are reasonably similar to that of IF and GF of Hydrohill respectively (dot dash lines in Figure 19a). This reveals that the variation of pH in interflow and saturated flow is not driven only by rainfall input but also by "the biological processes of photosynthesis and respiration [22, 60]." This is ascribed to the coupling of the three "reactors" mentioned above via their MC and IMC for transformation and exchange (Figure 4).

tions of the runoff responses from Hydrohill and monoliths (Figure 19) can be reasoned as the results of multi-coupling among hydrological, biological, and geochemical processes. It shows

Current Challenges in Experimental Watershed Hydrology

http://dx.doi.org/10.5772/55087

323

The measurements of various runoff components, within artificial systems containing controlled boundaries, provide the possibility to look inside the formation mechanisms of individual components. Isotopic and geochemical tracing can help to investigate these mechanisms but only if significant differences in isotopic compositions of components occur. The general mechanisms for these runoff components (i.e., the surface runoff, interflow and groundwater flow (saturated flow) ) are discussed in the following paragraphs [61,62].

*Surface runoff (SR).* Precipitation input is, of course, the essential condition for the generation of surface runoff. However, in order to actually generate runoff, there must be enough precipitation to form a thin, saturated soil layer (*Lsat*) at the ground surface. Saturation is key because unsaturated water movement is thought to be too slow to generate runoff. The thickness of *Lsat* at catchments, monoliths and plots was found to vary between 5 to 50 mm throughout events, increasing downward during rainfall events and receding once rainfall stopped, although these findings were complicated by the irregularities of the soil surface. SR was not generated until a *Lsat* was established at the surface. Once a *Lsat* was developed, regardless of its thickness, SR was generated immediately. *Overland flow* was only observed on impermeable surface (*DO*) at artificial plot and, on saturated surface (*SO*) from a special lysimeter designed for *Lsat* simulation. Intrastorm variation of isotopic composition of *DO* and *SO* can match that of event rainfall. In most cases however, SR is generated within the *Lsat*, with turbulent mixing of event and pre-event water stored in *Lsat*, i.e., the *saturated mixing surface flow* (*MS*); Alternatively, small amounts of event water act on the surface of *Lsat*, to force out pre-event water in the *Lsat*, termed here as the *saturated expelled surface flow* (*ES*). The isotopic composition of *MS* shows a mixing of event rainfall and pre-event water in *Lsat.* The isotopic

*Interflow (IF) in the unsaturated zone.* Three generation mechanisms are observed. (1) In cases there are soil layers with distinct bulk density and/or hydraulic conductivities, *IF* can be generated at the interface of soil layers. This was only observed in the 'mini CZEB', Nandadish, where IF occurred at the interface of the layers A (topsoil) and B (subsoil). This is termed as *layered interflow (LI)*. (2) During percolation, soil water moves downward and laterally towards the drainage interface and accumulated until saturation. The saturation will expand if the soil matrix potential ψ<0. Once ψ reaches zero, the accumulated saturated soil water will be released immediately at drainage interface. A small amount of infiltration water can trigger this process, to expel the accumulated soil water out as interflow. This process was observed

*(EI)*. (3) *Macropore interflow* (*MI*) was also observed in a special designed catchment without

*Groundwater flow (GF) from saturated zone.* In most cases, groundwater flow was due to a rising groundwater table resulted from event water recharge via the capillary fringe zone, and was

.

monolith. This is termed as *expelled interflow*

that the reactor II (Figure 4) provides a key operator in contaminant hydrology.

**5.5. Runoff generation**

composition of *ES* is similar to that of *Lsat*.

for most *IF* events in Hydrohill and the 32 m2

saturated zone with area of 4573 m2

**Figure 19.** Diel variations of pH, and hydrochemical constituents SO4 2-, Ca2+ in event rainfall and runoff responses of Hydrohill and monoliths. SR, IF, GF- runoff responses from Hydrohill; L1- monolith of 32 m2 serves as an element slope of Hydrohill; LS - saturated monolith of 1m2 serves as solum in reduction environment; EB3- Morning Glory Catchment of 4573 m2 serves as a catchment without saturated zone formed by debris only.

*Variation of ionic species.* (1) For anion SO4 2- : the diel variation curves of runoff responses of both Hydrohill and monoliths, with few exceptions, are contrary to that of rainfall. Mostly their peak concentrations happened during evening and midnight (Figure 19b). Different from the case of pH, SO4 2- of SR has a strong inversion with the event rainfall curve at night time. (2) For cation Ca2+: the diel curve of rainfall (Fig.19c) shows a small variation after sunrise. However, it triggers variations in runoff responses with their peak concentrations at both a.m. and before midnight. Highest peak happens to LI at afternoon. These diel variations in runoff perhaps are metabolism related. (3) Nimick et al [60] discussed diel cycles of dissolved trace metal concentration in a Rocky Mountain stream. They found that the anionic species "have their highest dissolved concentrations in the late afternoon" while the cationic species "have their highest dissolved concentrations shortly before sunrise"[23,60].

*Role of soil.* Variations of pH, SO4 2-, Ca2+ in total runoff of Morning Glory Catchment, a special designed catchment without soil but debris (EB3 in Figure 19 with broken lines), are very similar to that of rainfall curves after 20:00 (EB3 data are not enough before this time). The role of soil in the formation of ionic species in runoff is apparent. Even the variations of pH, SO4 2-, Ca2+ in total runoff of EB3 are triggered by the event rainfall, but the resultant variations appear much simpler than various curves of Hydrohill and monoliths. This implies that only a simple process (i.e., mainly the hydrological process is involved in EB3). So, the complex diel varia‐ tions of the runoff responses from Hydrohill and monoliths (Figure 19) can be reasoned as the results of multi-coupling among hydrological, biological, and geochemical processes. It shows that the reactor II (Figure 4) provides a key operator in contaminant hydrology.

## **5.5. Runoff generation**

IF and GF of Hydrohill respectively (dot dash lines in Figure 19a). This reveals that the variation of pH in interflow and saturated flow is not driven only by rainfall input but also by "the biological processes of photosynthesis and respiration [22, 60]." This is ascribed to the coupling of the three "reactors" mentioned above via their MC and IMC for transformation and

Hydrohill and monoliths. SR, IF, GF- runoff responses from Hydrohill; L1- monolith of 32 m2 serves as an element slope of Hydrohill; LS - saturated monolith of 1m2 serves as solum in reduction environment; EB3- Morning Glory Catchment

both Hydrohill and monoliths, with few exceptions, are contrary to that of rainfall. Mostly their peak concentrations happened during evening and midnight (Figure 19b). Different from

(2) For cation Ca2+: the diel curve of rainfall (Fig.19c) shows a small variation after sunrise. However, it triggers variations in runoff responses with their peak concentrations at both a.m. and before midnight. Highest peak happens to LI at afternoon. These diel variations in runoff perhaps are metabolism related. (3) Nimick et al [60] discussed diel cycles of dissolved trace metal concentration in a Rocky Mountain stream. They found that the anionic species "have their highest dissolved concentrations in the late afternoon" while the cationic species "have

designed catchment without soil but debris (EB3 in Figure 19 with broken lines), are very similar to that of rainfall curves after 20:00 (EB3 data are not enough before this time). The role of soil in the formation of ionic species in runoff is apparent. Even the variations of pH, SO4

Ca2+ in total runoff of EB3 are triggered by the event rainfall, but the resultant variations appear much simpler than various curves of Hydrohill and monoliths. This implies that only a simple process (i.e., mainly the hydrological process is involved in EB3). So, the complex diel varia‐

2- of SR has a strong inversion with the event rainfall curve at night time.

2-, Ca2+ in event rainfall and runoff responses of

2-,

2- : the diel variation curves of runoff responses of

2-, Ca2+ in total runoff of Morning Glory Catchment, a special

exchange (Figure 4).

**Figure 19.** Diel variations of pH, and hydrochemical constituents SO4

322 Current Perspectives in Contaminant Hydrology and Water Resources Sustainability

*Variation of ionic species.* (1) For anion SO4

*Role of soil.* Variations of pH, SO4

the case of pH, SO4

of 4573 m2 serves as a catchment without saturated zone formed by debris only.

their highest dissolved concentrations shortly before sunrise"[23,60].

The measurements of various runoff components, within artificial systems containing controlled boundaries, provide the possibility to look inside the formation mechanisms of individual components. Isotopic and geochemical tracing can help to investigate these mechanisms but only if significant differences in isotopic compositions of components occur. The general mechanisms for these runoff components (i.e., the surface runoff, interflow and groundwater flow (saturated flow) ) are discussed in the following paragraphs [61,62].

*Surface runoff (SR).* Precipitation input is, of course, the essential condition for the generation of surface runoff. However, in order to actually generate runoff, there must be enough precipitation to form a thin, saturated soil layer (*Lsat*) at the ground surface. Saturation is key because unsaturated water movement is thought to be too slow to generate runoff. The thickness of *Lsat* at catchments, monoliths and plots was found to vary between 5 to 50 mm throughout events, increasing downward during rainfall events and receding once rainfall stopped, although these findings were complicated by the irregularities of the soil surface. SR was not generated until a *Lsat* was established at the surface. Once a *Lsat* was developed, regardless of its thickness, SR was generated immediately. *Overland flow* was only observed on impermeable surface (*DO*) at artificial plot and, on saturated surface (*SO*) from a special lysimeter designed for *Lsat* simulation. Intrastorm variation of isotopic composition of *DO* and *SO* can match that of event rainfall. In most cases however, SR is generated within the *Lsat*, with turbulent mixing of event and pre-event water stored in *Lsat*, i.e., the *saturated mixing surface flow* (*MS*); Alternatively, small amounts of event water act on the surface of *Lsat*, to force out pre-event water in the *Lsat*, termed here as the *saturated expelled surface flow* (*ES*). The isotopic composition of *MS* shows a mixing of event rainfall and pre-event water in *Lsat.* The isotopic composition of *ES* is similar to that of *Lsat*.

*Interflow (IF) in the unsaturated zone.* Three generation mechanisms are observed. (1) In cases there are soil layers with distinct bulk density and/or hydraulic conductivities, *IF* can be generated at the interface of soil layers. This was only observed in the 'mini CZEB', Nandadish, where IF occurred at the interface of the layers A (topsoil) and B (subsoil). This is termed as *layered interflow (LI)*. (2) During percolation, soil water moves downward and laterally towards the drainage interface and accumulated until saturation. The saturation will expand if the soil matrix potential ψ<0. Once ψ reaches zero, the accumulated saturated soil water will be released immediately at drainage interface. A small amount of infiltration water can trigger this process, to expel the accumulated soil water out as interflow. This process was observed for most *IF* events in Hydrohill and the 32 m2 monolith. This is termed as *expelled interflow (EI)*. (3) *Macropore interflow* (*MI*) was also observed in a special designed catchment without saturated zone with area of 4573 m2 .

*Groundwater flow (GF) from saturated zone.* In most cases, groundwater flow was due to a rising groundwater table resulted from event water recharge via the capillary fringe zone, and was termed *recharged groundwater flow* (*RG*). RG was observed in most cases of both artificial and natural catchments. In addition, the *macropore-induced groundwater flow* (*MG*) is the only component not directly observed in this work but inferred from the isotope data in catchment *GF* during events.

Based on such problematic assumptions, one of the basic conclusions that "most stormflow is old water"[66], which resulted from two-component hydrograph separations in multiple

Current Challenges in Experimental Watershed Hydrology

http://dx.doi.org/10.5772/55087

325

*The physically ambiguous 'old' and 'new' water.* This two-component separation model in fact is based on the classic Horton infiltration theory with a very simple runoff generation concept that the soil surface acts as a sieve capable of separating rainfall into two basic components [67]. This 'old' and 'new' water separation model has led to many usages with a variety of defined components in addition to 'old' and 'new' water component of flow, e.g., 'pre-event and event component', 'pre-storm and storm component', 'pre-storm water and rain in storm run‐ off', 'pre-event and rainfall water', 'groundwater component and event water', 'groundwater and rainwater', 'groundwater discharge and surface runoff'[68], 'surface and subsurface runoff'[69]. The classic two-component method has also been extrapolated to incorporate three components (e.g. that of channel precipitation, soil water and groundwater [70]). Thus, it seems that success can be assured because there are no calibration constraints. As Kirchner [15] noted some models are "often good mathematical marionettes, they can dance to the tune of the

The 'old'(O) and 'new'(N) water actually correspond to different generation patterns of runoff components. In fact, multiple runoff mechanisms can result in the same proportions of old versus new water, a problem known as equi-finality. As seen in Hydrohill, surface and subsurface runoff patterns emerge from multiple runoff mechanisms (e.g., macropore flow can result in a large proportion of event water or a large proportion of pre-event water depending on antecedent wetness in the catchment). The runoff components corresponding to old water and new water can be very ambiguous. Applied to surface and subsurface runoff [68, 69], surface runoff is labeled new water [62]. This is a misconception as shown in Figure 20.

*The mass balance equations for this model are untenable.* The two-component tracer mixing model is stated by two mass balance equations for the composition of stormflow at any time which is used for hydrograph separation: Qs=Qn+Qo and Qsδs=Qnδn+Qoδ<sup>o</sup> where Q is streamflow, δ(Delta) is the D or 18O content, the subscripts s, n and o represent the stream, new and old respectively[56]. At any time t, the first equation is correct (i.e., discharge Qs(t) at the outlet of a watershed is equal to the sum of the unknown Qn(t) and Qo(t) right at the same outlet at the time t). During time t, the Qs(t) at the outlet and its isotopic concentrationδ<sup>s</sup> are measurable and are the value of known without problem, however both the Qn(t) and Qo(t) at the outlet of time t, are the results of confluence from somewhere within the watershed before time t, i.e., the results of a convolution integration. The Qn(t) and Qo(t) at the outlet of time t are concentrated from the separated isochrone area somewhere within the watershed with different time of concentration. So, the second equation with algebraic sum is physically unrealistic, hinging on the assumptions outlined above. For this linear algebraic sum to be physically realistic, the isotopic concentrations in the reservoirs of new and old water must be constant, without mixing or fractionating during the time interval of event hydrograph. To be physically realistic and representative of the mixing and non-conservative nature of water in natural catchments, it would likely take on the form of an integral or a finite-difference isochrone [63]. No affluxion

catchments with different drainage areas, appears suspect.

calibration data".

#### **5.6. Unreasonableness of current two-component isotope hydrograph separations**

Hydrograph separation using conservative, two-component mixing models has been done in various natural basins. However, data from these two catchments, including both the natural and artificial systems indicate that this technique is unacceptable for natural basins and artificial catchments, yielding misleading results due to the unreasonableness of most of the assumptions involved and of the unrealistic operation procedures [63].

*Violation of assumptions.* In studying natural processes it is inevitable to make simplifying assumptions. However, as Kennedy et al [64] indicated, when simplifying assumptions are made, one runs the risk of drawing misleading conclusions. That may have been the case for the five assumptions [65] summarized for the use of isotopic techniques in precipitation-runoff studies.

The first assumption, that "Groundwater and baseflow are characterized by a single isotopic content" [65], is problematic. Baseflow in natural basins can be recharged by the active zone, passive zone, each with differing flowpaths for recharge resulting in differences in the isotopic composition of baseflow itself, which changes with time and discharge. It is also unreasonable to expect a single isotopic content for groundwater which is subject to both temporal and spatial variations (Figure 16).

The second assumption, that "Rain or snowmelt can be characterized by a single constant isotopic composition or, the variations are documented" [65], is problematic, because, as shown in Figure 14a the intra-storm variability of isotopic composition of rainfall is very significant. High isotopic variability has also been observed in experimental watersheds at Maimai, Panola, etc [56]. Additionally, the spatial variability of isotopic composition of rainfall in a watershed exists and shouldn't be ignored. Isotopic data from a 82.1 km2 watershed showed that the largest difference in rainfallδ18O reached 8.2‰.

The third assumption, that "The isotopic composition of rain water is significantly different from that of groundwater/baseflow" [65], can be true. However it must be demonstrated by appropriate sampling [64].

The fourth assumption, that "Contributions from soil water are negligible, or the isotopic composition is identical to that of groundwater" [65] is problematic, because it is a miscon‐ ception. In a natural watershed, this can only occur if there is no unsaturated zone.

The fifth assumption, that "Contributions from surface water-bodies (such as ponds) are negligible" [65], is reasonable for most upland watersheds [64]. However, in catchments of hilly area with land use including series of ponds and paddy fields, the contributions are large enough to significantly influence results [62].

Based on such problematic assumptions, one of the basic conclusions that "most stormflow is old water"[66], which resulted from two-component hydrograph separations in multiple catchments with different drainage areas, appears suspect.

termed *recharged groundwater flow* (*RG*). RG was observed in most cases of both artificial and natural catchments. In addition, the *macropore-induced groundwater flow* (*MG*) is the only component not directly observed in this work but inferred from the isotope data in catchment

Hydrograph separation using conservative, two-component mixing models has been done in various natural basins. However, data from these two catchments, including both the natural and artificial systems indicate that this technique is unacceptable for natural basins and artificial catchments, yielding misleading results due to the unreasonableness of most of the

*Violation of assumptions.* In studying natural processes it is inevitable to make simplifying assumptions. However, as Kennedy et al [64] indicated, when simplifying assumptions are made, one runs the risk of drawing misleading conclusions. That may have been the case for the five assumptions [65] summarized for the use of isotopic techniques in precipitation-runoff

The first assumption, that "Groundwater and baseflow are characterized by a single isotopic content" [65], is problematic. Baseflow in natural basins can be recharged by the active zone, passive zone, each with differing flowpaths for recharge resulting in differences in the isotopic composition of baseflow itself, which changes with time and discharge. It is also unreasonable to expect a single isotopic content for groundwater which is subject to both temporal and

The second assumption, that "Rain or snowmelt can be characterized by a single constant isotopic composition or, the variations are documented" [65], is problematic, because, as shown in Figure 14a the intra-storm variability of isotopic composition of rainfall is very significant. High isotopic variability has also been observed in experimental watersheds at Maimai, Panola, etc [56]. Additionally, the spatial variability of isotopic composition of rainfall in a watershed exists and shouldn't be ignored. Isotopic data from a 82.1 km2 watershed

The third assumption, that "The isotopic composition of rain water is significantly different from that of groundwater/baseflow" [65], can be true. However it must be demonstrated by

The fourth assumption, that "Contributions from soil water are negligible, or the isotopic composition is identical to that of groundwater" [65] is problematic, because it is a miscon‐

The fifth assumption, that "Contributions from surface water-bodies (such as ponds) are negligible" [65], is reasonable for most upland watersheds [64]. However, in catchments of hilly area with land use including series of ponds and paddy fields, the contributions are large

ception. In a natural watershed, this can only occur if there is no unsaturated zone.

**5.6. Unreasonableness of current two-component isotope hydrograph separations**

assumptions involved and of the unrealistic operation procedures [63].

324 Current Perspectives in Contaminant Hydrology and Water Resources Sustainability

showed that the largest difference in rainfallδ18O reached 8.2‰.

*GF* during events.

studies.

spatial variations (Figure 16).

appropriate sampling [64].

enough to significantly influence results [62].

*The physically ambiguous 'old' and 'new' water.* This two-component separation model in fact is based on the classic Horton infiltration theory with a very simple runoff generation concept that the soil surface acts as a sieve capable of separating rainfall into two basic components [67]. This 'old' and 'new' water separation model has led to many usages with a variety of defined components in addition to 'old' and 'new' water component of flow, e.g., 'pre-event and event component', 'pre-storm and storm component', 'pre-storm water and rain in storm run‐ off', 'pre-event and rainfall water', 'groundwater component and event water', 'groundwater and rainwater', 'groundwater discharge and surface runoff'[68], 'surface and subsurface runoff'[69]. The classic two-component method has also been extrapolated to incorporate three components (e.g. that of channel precipitation, soil water and groundwater [70]). Thus, it seems that success can be assured because there are no calibration constraints. As Kirchner [15] noted some models are "often good mathematical marionettes, they can dance to the tune of the calibration data".

The 'old'(O) and 'new'(N) water actually correspond to different generation patterns of runoff components. In fact, multiple runoff mechanisms can result in the same proportions of old versus new water, a problem known as equi-finality. As seen in Hydrohill, surface and subsurface runoff patterns emerge from multiple runoff mechanisms (e.g., macropore flow can result in a large proportion of event water or a large proportion of pre-event water depending on antecedent wetness in the catchment). The runoff components corresponding to old water and new water can be very ambiguous. Applied to surface and subsurface runoff [68, 69], surface runoff is labeled new water [62]. This is a misconception as shown in Figure 20.

*The mass balance equations for this model are untenable.* The two-component tracer mixing model is stated by two mass balance equations for the composition of stormflow at any time which is used for hydrograph separation: Qs=Qn+Qo and Qsδs=Qnδn+Qoδ<sup>o</sup> where Q is streamflow, δ(Delta) is the D or 18O content, the subscripts s, n and o represent the stream, new and old respectively[56]. At any time t, the first equation is correct (i.e., discharge Qs(t) at the outlet of a watershed is equal to the sum of the unknown Qn(t) and Qo(t) right at the same outlet at the time t). During time t, the Qs(t) at the outlet and its isotopic concentrationδ<sup>s</sup> are measurable and are the value of known without problem, however both the Qn(t) and Qo(t) at the outlet of time t, are the results of confluence from somewhere within the watershed before time t, i.e., the results of a convolution integration. The Qn(t) and Qo(t) at the outlet of time t are concentrated from the separated isochrone area somewhere within the watershed with different time of concentration. So, the second equation with algebraic sum is physically unrealistic, hinging on the assumptions outlined above. For this linear algebraic sum to be physically realistic, the isotopic concentrations in the reservoirs of new and old water must be constant, without mixing or fractionating during the time interval of event hydrograph. To be physically realistic and representative of the mixing and non-conservative nature of water in natural catchments, it would likely take on the form of an integral or a finite-difference isochrone [63]. No affluxion happens to both the new and old water in the watershed. This can only happen in a small pond, operationally for a watershed, it's "the Procrustean bed [17]".

directly monitored hydrometrically; (4) isotopic and hydrochemical tracers are key to under‐ standing runoff generation mechanisms; and (5) account for interactions between the hydro‐

Current Challenges in Experimental Watershed Hydrology

http://dx.doi.org/10.5772/55087

327

Another kind of experimental basin is suggested going forward, namely the Critical Zone Experimental *Block* (CZEB), geologically a monolith-block within the Critical Zone. CZEB is a dynamic ecological and evolving system, coupled with various systems and united by hydrological process. The CZEB is a natural open system; both the energy and mass exchanges exist across its boundaries. It is a dissipative complex system with some degree of selforganization. The function of CZEB is threefold: mass, its material aspect; energy and force, its driven aspect and, organization and entropy, its thermodynamic/philosophical aspect.

To advance watershed hydrology and support development of a unified theory of hydrology, a two-way multi-scale experimental watershed system is suggested, including the natural system and the artificial system. Both have multi-scale subsystems from monolith, slope, subwatershed to watershed and follow the research idea of upwards and downwards routes. A trial for such a strategy is partly completed at the Chuzhou CZEB Experimental System.

To advance contaminant hydrology, the suggested two-way multi-scale experimental water‐ shed system may provide a key to unravel the complex mechanisms coupling hydrological, biological, and geochemical processes. Contaminant transformation and fate and their effects on regional degradation of groundwater basin involve highly complex mechanisms. This twoway multi-scale system calls for new models using both natural and artificial experimental basin results and using upwards and downwards approaches to multidisciplinary techniques,

All hydrological knowledge ultimately comes from observations, experiments, and measure‐ ments [15]. Progress in hydrology results mainly from challenges to prevailing approaches and concepts [71]. Hydrological experimentation, including CZEB experimental watershed studies, is the building block for development of a unified theory of hydrology including contaminant hydrology. However, it is important to remember Werner Heisenberg's warning that, "what we observe is not nature herself, but nature exposed to our method of questioning".

The renovation of the Chuzhou CZEB Experimental System is supported by the Hydrology Bureau of the Chinese Ministry of Water Resources and, the Nanjing Hydraulic Research Institutes. The adventurous plan for a two-way multi-scale hydrological experimental

Gu is deeply grateful to Jeffrey McDonnell and the group from USGS, Vance Kennedy, Carol Kendall, Norman (Jake) Peters for their kind help and support - the ever green cedars they planted by Hydrohill are now rooted and full of luxuriant foliage – as well as the wonderful water tracing methodology they taught, extending associations to new generation of collea‐ gues and students. The authors are also deeply grateful to those who have visited Hydrohill

watershed system and its realization are led by Acadimician Jian-Yun Zhang.

sphere and multiple watershed processes.

including isotope tracing.

**Acknowledgements**

**Figure 20.** Intrastorm isotopic variations of an event rainfall, and isotopic variations of observed surface runoff in dif‐ ferent watersheds. P–rainfall, SRI–surface runoff of the 'mini CZEB', Nandadish, with drainage area of 7897 m2, SRII– surface runoff of Hydrohill of 490 m2, SRIII- runoff of Morning Glory Catchment of 4573 m2.

## **6. Conclusion**

The river basin, watershed or catchment is a central concept in hydrology. Basin studies to assess watershed hydrology are approaching a period of transition and innovation. In fact, experimental watershed studies, the core of watershed hydrology, have tremendous and complicated challenges ahead.

The current challenges in basin experimental hydrology are mainly twofold. Advancing hydrologic science creates a fundamental challenge. Because the watershed system is a dynamic ecological system composed of a variety of biotic and abiotic processes driven by water and climatic processes, experimental watersheds should multi-couple these processes, organizing innovative measurements and approaches while continuing to support and test hydrological models. Anthropogenic hydrologic replumbing and natural climate oscillations are equally challenging. Field studies are the key to understanding and modeling the effect of hydrologic replumbing and climate change. There is a need for long term monitoring, systematic experimental facilities as well as data mining to reclaim historic data useful for determining baseline watershed metrics.

The main problems of the historical experimental basin approach are threefold: watershed surface laterality, hydrologic process laterality and downward components laterality. Lessons from fifty years of Chinese experimental basin studies are: (1) research facilities require natural conditions and also artificial controlled boundaries; (2) address not only the surface watershed, but also downward to the bedrock; (3) surface and subsurface runoff components should be directly monitored hydrometrically; (4) isotopic and hydrochemical tracers are key to under‐ standing runoff generation mechanisms; and (5) account for interactions between the hydro‐ sphere and multiple watershed processes.

Another kind of experimental basin is suggested going forward, namely the Critical Zone Experimental *Block* (CZEB), geologically a monolith-block within the Critical Zone. CZEB is a dynamic ecological and evolving system, coupled with various systems and united by hydrological process. The CZEB is a natural open system; both the energy and mass exchanges exist across its boundaries. It is a dissipative complex system with some degree of selforganization. The function of CZEB is threefold: mass, its material aspect; energy and force, its driven aspect and, organization and entropy, its thermodynamic/philosophical aspect.

To advance watershed hydrology and support development of a unified theory of hydrology, a two-way multi-scale experimental watershed system is suggested, including the natural system and the artificial system. Both have multi-scale subsystems from monolith, slope, subwatershed to watershed and follow the research idea of upwards and downwards routes. A trial for such a strategy is partly completed at the Chuzhou CZEB Experimental System.

To advance contaminant hydrology, the suggested two-way multi-scale experimental water‐ shed system may provide a key to unravel the complex mechanisms coupling hydrological, biological, and geochemical processes. Contaminant transformation and fate and their effects on regional degradation of groundwater basin involve highly complex mechanisms. This twoway multi-scale system calls for new models using both natural and artificial experimental basin results and using upwards and downwards approaches to multidisciplinary techniques, including isotope tracing.

All hydrological knowledge ultimately comes from observations, experiments, and measure‐ ments [15]. Progress in hydrology results mainly from challenges to prevailing approaches and concepts [71]. Hydrological experimentation, including CZEB experimental watershed studies, is the building block for development of a unified theory of hydrology including contaminant hydrology. However, it is important to remember Werner Heisenberg's warning that, "what we observe is not nature herself, but nature exposed to our method of questioning".

## **Acknowledgements**

happens to both the new and old water in the watershed. This can only happen in a small pond,

**Figure 20.** Intrastorm isotopic variations of an event rainfall, and isotopic variations of observed surface runoff in dif‐ ferent watersheds. P–rainfall, SRI–surface runoff of the 'mini CZEB', Nandadish, with drainage area of 7897 m2, SRII–

The river basin, watershed or catchment is a central concept in hydrology. Basin studies to assess watershed hydrology are approaching a period of transition and innovation. In fact, experimental watershed studies, the core of watershed hydrology, have tremendous and

The current challenges in basin experimental hydrology are mainly twofold. Advancing hydrologic science creates a fundamental challenge. Because the watershed system is a dynamic ecological system composed of a variety of biotic and abiotic processes driven by water and climatic processes, experimental watersheds should multi-couple these processes, organizing innovative measurements and approaches while continuing to support and test hydrological models. Anthropogenic hydrologic replumbing and natural climate oscillations are equally challenging. Field studies are the key to understanding and modeling the effect of hydrologic replumbing and climate change. There is a need for long term monitoring, systematic experimental facilities as well as data mining to reclaim historic data useful for

The main problems of the historical experimental basin approach are threefold: watershed surface laterality, hydrologic process laterality and downward components laterality. Lessons from fifty years of Chinese experimental basin studies are: (1) research facilities require natural conditions and also artificial controlled boundaries; (2) address not only the surface watershed, but also downward to the bedrock; (3) surface and subsurface runoff components should be

surface runoff of Hydrohill of 490 m2, SRIII- runoff of Morning Glory Catchment of 4573 m2.

**6. Conclusion**

complicated challenges ahead.

determining baseline watershed metrics.

operationally for a watershed, it's "the Procrustean bed [17]".

326 Current Perspectives in Contaminant Hydrology and Water Resources Sustainability

The renovation of the Chuzhou CZEB Experimental System is supported by the Hydrology Bureau of the Chinese Ministry of Water Resources and, the Nanjing Hydraulic Research Institutes. The adventurous plan for a two-way multi-scale hydrological experimental watershed system and its realization are led by Acadimician Jian-Yun Zhang.

Gu is deeply grateful to Jeffrey McDonnell and the group from USGS, Vance Kennedy, Carol Kendall, Norman (Jake) Peters for their kind help and support - the ever green cedars they planted by Hydrohill are now rooted and full of luxuriant foliage – as well as the wonderful water tracing methodology they taught, extending associations to new generation of collea‐ gues and students. The authors are also deeply grateful to those who have visited Hydrohill for their teaching, help and support, they are: Sklash M. from Canada; Geyh M.A., Plate E. and Seiler K-P from Germany; Gat J. from Israel; Shiklomanov I.A. from Russia; Verhagen B. from South Africa, Littlewood I. from UK, Kinzelbach W. from Switzerland. Thanks to post-doctoral researchers Ma Tao, Xu J-T for their figures. Many thanks to the Editor of this book, Paul Bradley, for his kind help and, encouragements.

ry. In: Xi R-Z, Gu W-Z, Seiler K-P (eds.) Research basins and hydrological planning: proceedings of the International Conference on Research Basins and Hydrological Planning, March 2004, Hefei, China. London: A.A. Balkema Publishers; (2004). ,

Current Challenges in Experimental Watershed Hydrology

http://dx.doi.org/10.5772/55087

329

[10] Toebes, C, & Ouryvaev, V. editors. Representative and Experimental Basins. Paris:

[11] Kelly, L. L, & Glymph, L. M. Experimental watersheds and hydrological research. In: Representative and Experimental Areas: proceedings of the International Symposi‐

[12] Kirby, C, Newson, M. D, & Gilman, K. editors. Plynlimon Research: The First Two

[13] Mcdonnell, J. J, Sivapalan, M, Vache, K, Dunn, S, Grant, G, Haggerty, R, Hinz, C, Hooper, R, Kirchner, J, Roderick, M. L, Selker, J, & Weiler, M. Moving beyond Heter‐ ogeneity and Process Complexity: A new Vision for Watershed Hydrology. Water

[14] Sivapalan, M. Pattern, Process and Function: Elements of a Unified Theory of Hy‐ drology at the Catchment Scale. In: Anderson MG (ed.) Encyclopedia of Hydrologi‐

[15] Kirchner, J. W. Getting the Right Answers for the Right Reasons: Linking Measure‐ ments, Analyses, and Models to Advance the Science of Hydrology. Water Resources

[16] NRCBasic Research Opportunities in Earth Science. Washington D.C.: The National

[17] Klemes, V. Dilettantism in Hydrology: Transition or Destiny? Water Resources Re‐

[19] Dooge, J. C. Looking for Hydrologic Laws. Water Resources Research (1986). S- 58S

[20] Eagleson, P. S. Climate, Soil, and Vegetation, 1. Introduction to Water Balance Dy‐

[22] Grayson, R. B, Moore, I. D, & Mcmahon, T. A. Physically Based Hydrologic Model‐ ing. 2. Is the Concept Realistic? Water Resources Research (1992). , 26(10), 2659-2666.

[23] National Research CouncilChallenges and Opportunities in the Hydrologic Sciences. Prepublication copy. Washington D.C.: The National Academies Press; (2012).

[18] Bras, R, & Eagleson, P. S. Hydrology, the Forgotten Earth Science. Eos (1987).

namics. Water Resources Research (1978). , 14, 705-712.

[21] Nash, J. E. Foreward. Journal of Hydrology (1988). v-viii.

[24] Vince, G. An Epoch Debate. Science (2011). , 334, 32-37.

um of Budapest, 28 September- 1 October. Paris: IASH; (1965).

Decades. Wallingford: Institute of Hydrology; (1991).

cal Science. John Wiley & Sons; (2005). , 193-219.

Resources Research (2007). W07301

Research (2006). W03S04

Academies Press; (2001).

search (1986). S-188S

22-31.

Unesco; (1970).

## **Author details**

Wei-Zu Gu1\*, Jiu-Fu Liu1 , Jia-Ju Lu1 and Jay Frentress2

\*Address all correspondence to: gweizu@163.com or gweizu@gmail.com

1 Institute for Hydrology and Water Resources, Nanjing Hydraulic Research Institutes, Nanjing, China

2 Oregon State University, Corvallis, USA

## **References**


ry. In: Xi R-Z, Gu W-Z, Seiler K-P (eds.) Research basins and hydrological planning: proceedings of the International Conference on Research Basins and Hydrological Planning, March 2004, Hefei, China. London: A.A. Balkema Publishers; (2004). , 22-31.

[10] Toebes, C, & Ouryvaev, V. editors. Representative and Experimental Basins. Paris: Unesco; (1970).

for their teaching, help and support, they are: Sklash M. from Canada; Geyh M.A., Plate E. and Seiler K-P from Germany; Gat J. from Israel; Shiklomanov I.A. from Russia; Verhagen B. from South Africa, Littlewood I. from UK, Kinzelbach W. from Switzerland. Thanks to post-doctoral researchers Ma Tao, Xu J-T for their figures. Many thanks to the Editor of this book, Paul

and Jay Frentress2

1 Institute for Hydrology and Water Resources, Nanjing Hydraulic Research Institutes,

[1] Rodda, J. C. Basin Studies. In: Rodda JC (ed.) Facets of Hydrology. London: John Wi‐

[2] Biswas, A. K. History of Hydrology. North-Holland Publishing Company; (1970).

[3] Unesco/WMO/IAHSThree Centuries of Scientific Hydrology. Paris: Unesco; (1974).

[4] Bates, C. G, & Henry, A. J. Forest and Stream-flow Experiment at Wagon Wheel Gap.

[5] Урываев ВА. Зкспериментальные Гидрологические Исследования на Валдае.

[6] Douglass, J. E, & Hoover, M. D. Forest Hydrology and Ecology at Coweeta. New

[7] Delfs, J, Friedrich, W, Kiesekamp, H, & Wagenhoff, A. Der Einfluss des Waldes und des Kahlschlages auf den Abflussvorgane den Wasserhaushalt und den Bodenab‐

[8] Jacquet, J. Les Etudes d'Hydrologie Analytique sur Bassins Versant Experimentaux.

[9] Gu, W-Z, Liu, C-M, Song, X-F, Yu, J-J, & Xia, J. Hydrological experimental system and environmental isotope tracing: a review on the occasion of the 50th Anniversary of Chinese basin studies and the 20th Anniversary of Chuzhou Hydrology Laborato‐

Bull. Du Centre de Recherche et d'Essais de Chatou. (1962). , 2, 3-25.

Bradley, for his kind help and, encouragements.

2 Oregon State University, Corvallis, USA

ley & Sons; (1976). , 257-297.

York: Springer-Verlag. (1988).

trag. Aus dem Walde. (1958). , 3, 223-325.

, Jia-Ju Lu1

328 Current Perspectives in Contaminant Hydrology and Water Resources Sustainability

\*Address all correspondence to: gweizu@163.com or gweizu@gmail.com

Monthly Weather Review. Supplement.(1928). , 30, 1-79.

Ленинград: Гидрометеорологическое Иэдательство; 1953

**Author details**

Nanjing, China

**References**

Wei-Zu Gu1\*, Jiu-Fu Liu1


[25] Alley, R. B, Clark, P. U, Keigwin, L. D, & Webb, R. S. Making Sense of Millennial-Scale Climate Change. In: Clark PU, Webb RS, Keigwin LD (eds.) Geophys. Monogr. Ser., Mechanisms of Global Climate Change at Millennial Time Scale. Washington D.C.: AGU; (1999). , 112

ration of the Critical Zone: Report of a workshop sponsored by the National Science

Current Challenges in Experimental Watershed Hydrology

http://dx.doi.org/10.5772/55087

331

[40] Report prepared by the CZO CommunityFuture directions for Critical Zone observa‐

[42] Klemes, V. Conceptualization and Scale in Hydrology. Journal of Hydrology (1983). [43] Seiler, K-P, & Lindner, W. Near Surface and Deep Groundwater. Journal of Hydrolo‐

[44] Arnold, R. W, Szabolcs, I, & Targulian, V. O. Global Soil Change. Laxenburg: Interna‐

[46] Weinberg, S. Dreams of a Final Theory: The Scientist's Search for the Ultimate Laws

[47] Anderson, S. P, Von Blanckenburg, F, & White, A. F. Physical and Chemical Controls

[48] Toffler, A. Forward: Science and Change. In: Prigogine I, Stengers I (eds.). Order out

[49] Gu, W-Z. Field research on surface water and subsurface water relationships in an artificial experimental catchment. In: Dahlblom P, Lindh G. (eds.) Interaction be‐ tween groundwater and surface water: proceedings of the international symposium,

[50] Gu, W-Z, & Freer, J. Patterns of surface and subsurface runoff generation. In: Lei‐ bundgut C (ed.) Tracer technologies for hydrological systems: proceedings of sympo‐ sium H4 at the XXI General Assembly of the International Union of Geodesy and

Geophysics, July 1995, Boulder, USA. IAHS Publ. (1995). , 1995(229), 265-273.

[51] Gu, W-Z. Measurements of spatial evapotranspiration characteristics of an experi‐ mental basin using a neutron probe. In: Isotope techniques in water resources devel‐ opment 1987: proceedings of the International Symposium on the Use of Isotope Techniques in Water Resources Development, 30 March- 3 April, Vienna. Vienna:

[53] Kendall, C, Mcdonnell, J. J, & Gu, W-Z. A Look Inside 'Black Box' Hydrograph Sepa‐ ration Models: A Study at the Hydrohill Catchment. Hydrological Processes (2001). ,

[54] Kendall, C, & Gu, W-Z. Development of isotopically heterogeneous infiltration wa‐ ters in an artificial catchment in Chuzhou, China. In: Isotope techniques in water re‐

30 May- 3 June, Ystad, Sweden. Lund: Lund University; (1988).

[45] Bateson, G. Mind and Nature: A Necessary Unity. New York: Dutton; (1979).

tional Institute for Applied Systems Analysis; (1990).

of Nature. London: Vintage Books; (1994).

of Chaos. Toronto: Bantam Books; (1984).

IAEA; (1987). , 1987, 789-793.

15, 1877-1902.

[52] Mcdonnell, J. Personal communications.

on the Critical Zone. Elements (2007). , 3, 315-319.

[41] National Critical Zone Observatory Program: http://criticalzoneorg/index.html

Foundation. (2005).

tory (CZO) science. (2010).

gy (1995). , 165, 33-44.


ration of the Critical Zone: Report of a workshop sponsored by the National Science Foundation. (2005).


[25] Alley, R. B, Clark, P. U, Keigwin, L. D, & Webb, R. S. Making Sense of Millennial-Scale Climate Change. In: Clark PU, Webb RS, Keigwin LD (eds.) Geophys. Monogr. Ser., Mechanisms of Global Climate Change at Millennial Time Scale. Washington

[26] Alley, R. B, Mayewski, P. A, Sowers, T, Stuiver, M, Taylor, K. C, & Clark, P. U. Holo‐ cene Climatic Instability: A Prominent, Widespread Event 8200 yr ago. Geology

[27] Adams, W. M. Green Development, Environment and Sustainability in the Third

[28] Scudder, T. A Sociological Framework for the Analysis of New Land Settlements. In: Cernea M (ed.) Putting People First: Sociological Variables in Rural Development.

[29] Gu, W-Z, & Lu, J-J. The Disposition of Water Resources of Arid Area with Special Reference to the Alxa Plateau, Inner Mongolia. In: Maldini D, Maher DM, Troppoli D, Studer M, Goebel J (eds.) Translating Scientific Results into Conservation Actions: New Roles, Challenges and Solutions for 21st Century Scientists. Boston: Earthwatch

[31] Kirchner, J. W. A Double Paradox in Catchment Hydrology and Geochemistry. Hy‐

[33] Steffen, W, Crutzen, P. J, & Mcneill, J. R. The Anthropocene: Are Humans Now Over‐

[35] Li, F-X, & Yao, J-H. Synthetical Research on the Economic Development and the En‐ vironmental Rectification for the He-Xi Corridor area (in Chinese). Beijing: Environ‐

[36] Lin, H. Earth's Critical Zone and Hydropedology: Concepts, Characteristics, and Ad‐

[37] Falkenmark, M, & Rockstrom, J. Balancing Water for Humans and Nature: The New

[38] Gu, W-Z. Experimental and representative basin stuides in China. In: Hooghart JC, Posthumus CWS, Warmerdam PMM. (eds.) Hydrological research basins and the en‐ vironment: proceedings of the international conference, September 1990, Wagenin‐

[39] Brantley, S. L, White, T. S, White, A. F, Sparks, D, Richter, D, Pregitzer, K, Derry, L, Chorover, J, Chadwick, O, April, R, Anderson, S, & Amundson, R. Frontiers in explo‐

[34] Smakhtin, V. U. Low Flow Hydrology: A Review. Journal of Hydrology (2001).

[30] Hall, P. Great Planning Disasters. London: Weidenfeld and Nicolson; (1980).

[32] Gu, W-Z. Isotope Hydrology (in Chinese). Beijing: Science Press; (2011).

whelming thr Great Forces of Nature? Ambio (2007). , 36, 614-642.

vances. Hydrology and Earth System Sciences (2010). , 14, 25-45.

Approach in Ecohydrology. London: Earthscan; (2004).

gen, The Netherland. Hague: CHO; (1990). , 24-28.

D.C.: AGU; (1999). , 112

(1997). , 25, 483-486.

Institute; (2007).

World. London: Routledge; (2001).

Oxford: Oxford University Press; (1991). , 148-167.

330 Current Perspectives in Contaminant Hydrology and Water Resources Sustainability

drological Processes (2003). , 17, 871-874.

mental Science Press; (1998).


sources development 1991: proceedings of the International Symposium on Isotope Techniques in Water Resources Development, 11-15 March, Vienna. Vienna: IAEA; (1992). , 1992, 61-73.

[65] Sklash, M. G, & Farvolden, R. N. The Use of Environmental Isotopes in the Study of High-runoff Episodes in Streams. In: Perry EC, Montgomery CW. (eds.) Isotope Studies of Hydrologic Processes. Dekalb: N. Illinois University Press; (1982). , 65-73.

Current Challenges in Experimental Watershed Hydrology

http://dx.doi.org/10.5772/55087

333

[66] Bishop, K. H. Episodic increases in stream acidity, catchment flow pathways and hy‐ drograph separation. PhD thesis. University of Cambridge, Dep. Geol., Jesus College,

[67] Chorley, R. J. The Hillslope Hydrological Cycle. In: Kirkby MJ. (ed.) Hillslope Hy‐

[68] Crouzet, E, & Hubert, P. Oliver Ph, Siwertz E. Le Tritium dans les Mesures d'Hydro‐ logie de Surface. Determination Experimentale du Coefficient de Ruissellement. Jour‐

[69] Wels, C, Cornett, R. J, & Lazerte, B. D. Hydrograph Separation: A Comparison of Ge‐ ochemical and Isotope Tracers. Journal of Hydrology (1991). , 122, 253-274.

[70] Dewalle, D. R, Swistock, B. R, & Sharpe, W. E. Three Component Tracer Model for Stormflow on a Small Appalachian Forested Catchment. Journal of Hydrology

[71] Yevjevich, V. Misconceptions in Hydrology and Their Consequences. Water Resour‐

drology. Chicheter: John Wiley & Sons; (1979). , 1-42.

nal of Hydrology (1970). , 11, 217-229.

ces Research (1968). , 4(2), 225-232.

Cambridge; (1991).

(1988). , 104, 301-310.


[65] Sklash, M. G, & Farvolden, R. N. The Use of Environmental Isotopes in the Study of High-runoff Episodes in Streams. In: Perry EC, Montgomery CW. (eds.) Isotope Studies of Hydrologic Processes. Dekalb: N. Illinois University Press; (1982). , 65-73.

sources development 1991: proceedings of the International Symposium on Isotope Techniques in Water Resources Development, 11-15 March, Vienna. Vienna: IAEA;

[55] Wood, E. F, Sivapalan, M, Beven, K. J, & Band, L. Effects of Spatial Variability and Scale with Implications to Hydrologic Modeling. Journal of Hydrology (1988). , 102,

[56] Kendall, C, & Mcdonnell, J. J. Effect of intrastorm isotopic heterogeneities of rainfall, soil water, and groundwater on runoff modeling. In: Tracers in hydrology: proceed‐ ings of the International Symposium on Tracers in Hydrology, 11-23 July 1993, Yoko‐

[57] Gu, W-Z. Experimental Research on Catchment Runoff Responses Traced by Envi‐ ronmental Isotopes. Advances in Water Science (1992). in Chinese with English ab‐

[58] Gu, W-Z, Shang, M-T, Zhai, S-Y, Lu, J-J, Jason, F, Mcdonnell, J. J, & Kendall, C. Rain‐ fall-runoff Paradox from A Natural Experimental Catchment. Advances in Water Sci‐

[59] Gu, W-Z, Lu, J-J, Zhao, X, & Peters, N. E. Responses of Hydrochemical Inorganic Ions in the Rainfall-runoff Processes of the Experimental Catchments and Its Significance for Tracing. Advances in Water Science (2007). in Chinese with English abstracts); ,

[60] Nimick, D. A, Cleasby, T. E, & Mccleskey, R. B. Seasonality of Diel Cycles of Dis‐ solved Trace Metal Concentrations in a Rocky Mountain Stream. Environmental Ge‐

[61] Gu, W-Z. Challenge on some rainfall-runoff conceptions traced by environmental isotopes in experimental catchments. In: Hotzl H, Werner A. (eds.) Tracer hydrology: proceedings of the 6th International Symposium on Water Tracing, 21-26 September

[62] Gu, W-Z. Various Patterns of Basin Runoff Generation Identified by Hydrological Experimentation and Water Tracing using Environmental Isotopes. Journal of Hy‐

[63] Gu, W-Z. Unreasonableness of current two-component isotopic hydrograph separa‐ tion for natural basins. In: Isotope techniques in water resources development 1995: proceedings of the International Symposium on Isotopes in Water Resources Devel‐

[64] Kennedy, V. C, Kendall, C, Zellweger, G. W, Wyerman, T. A, & Avanzino, R. J. De‐ termination of the Components of Stormflow using Water Chemistry and Environ‐ mental Isotopes, Mattole River Basin, California. Journal of Hydrology (1986). , 84,

1992, Karlsruhe, Germany. Rotterdam: A.A. Balkema; (1992). , 1992, 397-403.

draulic Engineering (1995). in Chinese with English abstracts); , 5, 9-17.

opment, 20-24 March, 1995, Vienna. Vienna: IAEA; (1996). , 1996, 261-264.

ence (2010). in Chinese with English abstracts); , 21(4), 471-478.

hama, Japan. IAHS Publ. (1993). , 1993(215), 41-48.

332 Current Perspectives in Contaminant Hydrology and Water Resources Sustainability

(1992). , 1992, 61-73.

stracts); , 3(4), 246-254.

ology (2005). , 47(5), 603-614.

18(1), 1-7.

107-140.

29-47.


## *Edited by Paul M. Bradley*

Human society depends on liquid freshwater resources to meet drinking, sanitation and hygiene, agriculture, and industry needs. Improved resource monitoring and better understanding of the anthropogenic threats to freshwater environments are critical to efficient management of freshwater resources and ultimately to the survival and quality of life of the global human population. This book helps address the need for improved freshwater resource monitoring and threat assessment by presenting current reviews and case studies focused on the fate and transport of contaminants in the environment and on the sustainability of groundwater and surface-water resources around the world. It is intended for students and professionals working in hydrology and water resources management.

Photo by JillianSuzanne / iStock

Current Perspectives in Contaminant Hydrology and Water Resources Sustainability

Current Perspectives in

Contaminant Hydrology

and Water Resources

Sustainability

*Edited by Paul M. Bradley*