Hydrology and Earth System Sciences www.hydrol-earth-syst-sci.net 1027-5606 1607-7938 11 2 2007 10.5194/hess-11-819-2007 http://www.hydrol-earth-syst-sci.net/11/819/2007/ http://www.hydrol-earth-syst-sci.net/11/819/2007/hess-11-819-2007.html http://www.hydrol-earth-syst-sci.net/11/819/2007/hess-11-819-2007.pdf 819 849 2007-02-05 Predictions of rainfall-runoff response and soil moisture dynamics in a microscale catchment using the CREW model H. Lee haksu.lee@noaa.gov E. Zehe M. Sivapalan School of Environmental Systems Enineering, The University of Western Australia, 35 Stirling Highway, Crawley, WA 6009, Australia Institute of Geoecology, University of Potsdam, Germany Departments of Geography & Civil and Environmental Engineering, University of Illinois at Urbana-Champaign, 220~Davenport Hall, 607 S. Mathews Avenue, Urbana, IL 61801, USA Predictions of catchment hydrology have been performed generally using either physically based, distributed models or conceptual lumped or semi-distributed models. In recognition of the disadvantages of using either of these modeling approaches, namely, detailed data requirements in the case of distributed modeling, and lack of physical basis of conceptual/lumped model parameters, Reggiani et al. (1998, 1999) derived, from first principles and in a general manner, the balance equations for mass, momentum and energy at what they called the Representative Elementary Watershed (or REW) scale. However, the mass balance equations of the REW approach include mass exchange flux terms which must be defined externally before their application to real catchments. Developing physically reasonable "closure relations" for these mass exchange flux terms is a crucial pre-requisite for the success of the REW approach. As a guidance to the development of closure relations expressing mass exchange fluxes as functions of relevant state variables in a physically reasonable way, and in the process effectively parameterizing the effects of sub-grid or sub-REW heterogeneity of catchment physiographic properties on these mass exchange fluxes, this paper considers four different approaches, namely the field experimental approach, a theoretical/analytical approach, a numerical approach, and a hybrid approach combining one or more of the above. Based on the concept of the scaleway (Vogel and Roth, 2003) and the disaggregation-aggregation approach (Viney and Sivapalan, 2004), and using the data set from Weiherbach catchment in Germany, closure relations for infiltration, exfiltration and groundwater recharge were derived analytically, or on theoretical grounds, while numerical experiments with a detailed fine-scale, distributed model, CATFLOW, were used to obtain the closure relationship for seepage outflow. The detailed model, CATFLOW, was also used to derive REW scale pressure-saturation (i.e., water retention curve) and hydraulic conductivity-saturation relationships for the unsaturated zone. Closure relations for concentrated overland flow and saturated overland flow were derived using both theoretical arguments and simpler process models. In addition to these, to complete the specification of the REW scale balance equations, a relationship for the saturated area fraction as a function of saturated zone depth was derived for an assumed topography on the basis of TOPMODEL assumptions. These relationships were used to complete the specification of all of the REW-scale governing equations (mass and momentum balance equations, closure and geometric relations) for the Weiherbach catchment, which are then employed for constructing a numerical watershed model, named the Cooperative Community Catchment model based on the Representative Elementary Watershed approach (CREW). CREW is then used to carry out sensitivity analyses with respect to various combinations of climate, soil, vegetation and topographies, in order to test the reasonableness of the derived closure relations in the context of the complete catchment response, including interacting processes. These sensitivity analyses demonstrated that the adopted closure relations do indeed produce mostly reasonable results, and can therefore be a good basis for more careful and rigorous search for appropriate closure relations in the future. Three tests are designed to assess CREW as a large scale model for Weiherbach catchment. The first test compares CREW with distributed model CATFLOW by looking at predicted soil moisture dynamics for artificially designed initial and boundary conditions. The second test is designed to see the applicabilities of the parameter values extracted from the upscaling procedures in terms of their ability to reproduce observed hydrographs within the CREW modeling framework. The final test compares simulated soil moisture time series predicted by CREW with observed ones as a way of validating the predictions of CREW. The results of these three tests, together, demonstrate that CREW could indeed be an alternative modelling framework, producing results that are consistent with those of the distributed model CATFLOW, and capable of ultimately representing processes actually occurring at the larger scale in a physically sound manner. Attinger, S.: Generalized coarse graining procedures for flow in porous media, Comp. Geosci., 7, 253–273, 2003. Beven, K. J.: Towards an alternative blueprint for a physically based digitally simulated hydrologic response modelling system, Hydrol. Processes, 16, 189–206, 2002. Beven, K. J. and Kirkby, M. J.: A physically-based variable contributing area model of basin hydrology, Hydrol. Sci. Bull., 24(1), 43–69, 1979. Beven, K. J. and Freer, J. E.: Equifinality, data assimilation, and uncertainty estimation in mechanistic modeling of complex environmental systems using the GLUE methodology, J. Hydrol., 249, 11–29, 2001. Binley, A., Elgy, J., and Beven, K.: A physically based model for heterogeneous hillslopes. 1. Runoff production, Water Resour. Res., 25(6), 1219–1226, 1989a. Binley, A., Beven, K., and Elgy, J.: A physically based model for heterogeneous hillslopes. 2. Effective hydraulic conductivities, Water Resour. Res., 25(6), 1219–1226, 1989b. Blöschl, G. and Zehe, E.: On hydrological predictability, Hydrol. Processes, 19(9), 3923–3929, 2005. Bras, R. L.: Hydrology: An Introduction to Hydrologic Science, Reading, Mass., USA: Addison-Wesley-Longman, 1990. Bresler, E. and Dagan, G.: Unsaturated flow in spatially variable fields. 2. Application of water flow models to various fields, Water Resour. Res., 19(2), 421–428, 1983. Chow, V. T., Maidment, D. R., and Mays, L. W.: Applied Hydrology, McGraw-Hill, Inc, 1988. Dagan, G.: Flow and Transport in Porous Formations, Springer-Verlag, New York, 1989. Duffy, C. J.: A two-state integral-balance model for soil moisture and groundwater dynamics in complex terrain, Water Resour. Res., 32(8), 2421–2434, 1996. Eagleson, P. S.: Climate, soil, and vegetation, 1. Introduction to water balance dynamics, Water Resour. Res., 14(5), 705–712, 1978a. Eagleson, P. S.: Climate, soil, and vegetation, 3. A simplified model of soil moisture movement in the liquid phase, Water Resour. Res., 14(5), 722–730, 1978b. Eagleson, P. S.: Climate, soil, and vegetation, 4. The expected value of annual evapotranspiration, Water Resour. Res., 14(5), 731–739, 1978c. Ichikawa, Y. and Shiiba, M.: Lumping of kinematic wave equation considering field capacity, 3rd Int. Conf. on Water Resources and Environmental Research, 1, 61–65, 2002. Kees, C. E., Band, L. E., and Farthing, M. W.: Effects of Dynamic Forcing on Hillslope Water Balance Models, Technical Report CRSC-TR04-12, Center for Research in Scientific Computation, Department of Mathematics, North Carolina State University, 2004. Kees, C. E., Farthing, M. W., Band, L. E., and Miller, C. T.: Choices of scale and process complexity in hillslope models, Eos Trans. AGU, 83(47), Fall Meet. Suppl., Abstract H62B-0838, 2002. Lee, H., Sivapalan, M., and Zehe, E.: Representative Elementary Watershed (REW) approach, a new blueprint for distributed hydrologic modelling at the catchment scale, in: Predictions in ungauged basins: international perspectives on state-of-the-art and pathways forward, edited by: Franks, S. W., Sivapalan, M., Takeuchi, K., and Tachikawa, Y., IAHS Press, Wallingford, Oxon, UK, IAHS Publications 301, Paper 15, 2005a. Lee, H., Sivapalan, M., and Zehe, E.: Representative Elementary Watershed (REW) approach, a new blueprint for distributed hydrologic modelling at the catchment scale: the development of closure relations, in: Predicting ungauged streamflow in the mackenzie river basin: today's techniques & tomorrow's solutions, edited by: Spence, C., Pomeroy, J. W., and Pietroniro, A., Canadian Water Resources Association (CWRA), Ottawa, Canada, 165–218, 2005b. Lunati, I., Attinger, S., and Kinzelbach, W.: Macrodispersivity for transport in arbitrary nonuniform flow fields: Asymptotic and pre-asymptotic results, Water Resour. Res., 38(10), 1187, doi:10.1029/2001WR001203, 2002. Maurer, T.: Physikalisch begrundete, zeitkontinuierliche Modellierung des Wassertransports in kleinen landlichen Einzugsgebieten, Universitat Karlsruhe, MitteiLungen Inst. F. Hydrologie u. Wasserwirtschaft, H. 61, Universitat Karlsruhe, 1997. Mualem, Y.: A new model for predicting the hydraulic conductivity of unsaturated porous media, Water Resour. Res., 12, 513–522, 1976. Philip, J. R.: General method of exact solution of the concentration-dependent diffusion equation, Aust. J. Phys., 13(1), 1–12, 1960. Press, W. H., Teukolsky, S. A., Vetterling, W. T., and Flannery, B. P.: Numerical Recipes in FORTRAN, Cambridge University Press, 1992. Reggiani, P. and Schellekens, J.: Modelling of hydrological responses: the representative elementary watershed approach as an alternative blueprint for watershed modelling, Hydrol. Processes, 17, 3785–3789, 2003. Reggiani, P., Sivapalan, M., and Hassanizadeh, S. M.: A unifying framework for watershed thermodynamics: balance equations for mass, momentum, energy and entropy and the second law of thermodynamics, Adv. Water Resour., 22(4), 367–398, 1998. Reggiani, P., Hassanizadeh, S. M., Sivapalan, M., and Gray, W. G.: A unifying framework for watershed thermodynamics: constitutive relationships, Adv. Water Resour., 23(1), 15–39, 1999. Reggiani, P., Sivapalan, M., and Hassanizadeh, S. M.: Conservation equations governing hillslope responses: exploring the physical basis of water balance, Water Resour. Res., 36(7), 1845–1863, 2000. Reggiani, P., Sivapalan, M., Hassanizadeh, S. M., and Gray, W. G.: Coupled equations for mass and momentum balance in a stream network: theoretical derivation and computational experiments, Proc. R. Soc. Lond., 457, 157–189, 2001. Robinson, J. S. and Sivapalan, M.: Catchment-scale runoff generation model by aggregation and similarity analyses, Hydrol. Processes, 9, 555–574, 1995. Robock, A., Vinnikov, K. Y., Srinivasan, G., Entin, J. K., Hollinger, S. E., Speranskaya, N. A., Liu, S., and Namkhai, A.: The Global Soil Moisture Data Bank. Bull. Amer. Meteor. Soc., 81, 1281–1299, 2000. Rogers, A. D.: The development of a simple infiltration capacity equation for spatially variable soils, B.E (Honours) thesis, Department of Civil and Environmental Engineering, The University of Western Australia, 1992. Schäfer, D.: Bodenhydraulische Eigenschaften eines Kleineinzugsgebiets – Vergleich und Bewertung unterschiedlicher Verfahren, PhD dissertation, Institute of Hydromechanics, University of Karlsruhe, Germany, 1999. Schulz, K., Seppelt, R., Zehe, E., Vogel, H. J., and Attinger, S.: Importance of spatial structures in advancing hydrological sciences, Water Resour. Res., 42(3), W03S03, doi:10.1029/2005WR004301, 2006. Sherman, L. K.: Comparison of F-curves derived by the methods Sharp and Holtan and of Sherman and Mayer, Trans. Am. Geophys. Union, 24, 465–467, 1943. Sivapalan, M. and Wood, E. F.: Spatial heterogeneity and scale in the infiltration response of catchments, in: Scale problems in hydrology, edited by: Gupta, V. K., Rodríguez-Iturbe I., and Wood, E. F., Reidel Publ., Dordrecht, pp 81–106, 1986. Sivapalan, M.: Linking hydrologic parameterizations across a range of spatial scales: hillslope to catchment to region, IAHS Publ. 212, Proc. Yokohama Symposium, 115–123, 1993. Sivapalan, M., Takeuchi, K., Franks, S.W., Gupta, V.K., Karambiri, H., Lakshmi, V., Liang, X., McDonnell, J.J., Mendiondo, E.M., O'Connel, P.E., Oki, T., Pomeroy, J.W., Schertzer, D., Uhlenbrook, S., and Zehe, E.: IAHS decade on Predictions of Ungauged Basins (PUB): Shaping an exciting future for the hydrological sciences, Hydrol. Sci. J., 48(6), 857–879, 2003. Sloan, P. G. and Moore, I. D.: Modeling subsurface stormflow on steeply sloping forested watersheds, Water Resour. Res., 20(12), 1815–1822, 1984. van Genuchten, M. T.: A closed-form equation for predicting the hydraulic conductivity of unsaturated soils, Soil Sci. Soc. Am. Jour., 44, 892–898, 1980. Viney, N. R. and Sivapalan, M.: A framework for scaling of hydrologic conceptualisations based on a disaggregation-aggregation approach, Hydrol. Processes, 18, 1395–1408, 2004. Vogel, H. J. and Roth, K.: Moving through scales of flow and transport in soil, J. Hydrol., 272, 95–106, 2003. Wagener, T., McIntyre, N., Lees, M. J., Wheater, H. S., and Gupta, H. V.: Towards reduced uncertainty in conceptual rainfall-runoff modelling: Dynamic identifiability analysis, Hydrol. Processes, 17, 455–476, 2003. Western, A. W. and Grayson, R. B.: The Tarrawarra data set: soil moisture patterns, soil characteristics and hydrological flux measurements, Water Resour. Res. 34(10), 2765–2768, 1998. Zehe, E., Maurer, T., Ihringer, J., and Plate, E.: Modeling water flow and mass transport in a Loess catchment, Phys. Chem. Earth (B), 26, 7–8, 487–507, 2001. Zehe, E. and Blöschl, G.: Predictability of hydrologic response at the plot and catchment scales – the role of initial conditions, Water Resour. Res., 40(10), W10202, doi:10.1029/2003WR002869, 2004. Zehe, E., Lee, H., and Sivapalan, M.: Derivation of closure relations and commensurate state variables for mesoscale hydrological models using dynamical upscaling, in: Predictions in ungauged basins: international perspectives on state-of-the-art and pathways forward, edited by: Franks, S. W., Sivapalan, M., Takeuchi, K., and Tachikawa, Y., IAHS Press, Wallingford, Oxon, UK, IAHS Publications 301, Paper 14, 2005a. Zehe, E., Becker, R., Bardossy, A., and Plate, E.: Uncertainty of simulated catchment runoff response in the presence of threshold processes: Role of initial soil moisture and precipitation, J. Hydrol., 315, 183–202, 2005b. Zehe, E., Lee, H., and Sivapalan, M.: Dynamical process upscaling for deriving catchment scale state variables and constitutive relations for meso-scale process models, Hydrol. Earth Syst. Sci., 10, 981–996, 2006.