Application of hydrological models in a snowmelt region of Aksu River Basin

This study simulated and predicted the runoff of the Aksu River Basin, a typical river basin supplied by snowmelt in an arid mountain region, with a limited data set and few hydrological and meteorological stations. Two hydrological models, the snowmelt-runoff model （SRM） and the Danish NedbФr-AfstrФmnings rainfall-runoff model （NAM）, were used to simulate daily discharge processes in the Aksu River Basin. This study used the snow-covered area from MODIS remote sensing data as the SRM input. With the help of ArcGIS software, this study successfully derived the digital drainage network and elevation zones of the basin from digital elevation data. The simulation results showed that the SRM based on MODIS data was more accurate than NAM. This demonstrates that the application of remote sensing data to hydrological snowmelt models is a feasible and effective approach to runoff simulation and prediction in arid unguaged basins where snowmelt is a major runoff factor.

A review of integrated surface-subsurface numerical hydrological models

Hydrological modeling,leveraging mathematical formulations to represent the hydrological cycle,is a pivotal tool in representing the spatiotemporal dynamics and distribution patterns inherent in hydrology.These models serve a dual purpose:they validate theoretical robustness and applicability via observational data and project future trends,thereby bridging the understanding and prediction of natural processes.In rapid advancements in computational methodologies and the continuous evolution of observational and experimental techniques,the development of numerical hydrological models based on physicallybased surface-subsurface process coupling have accelerated.Anchored in micro-scale conservation principles and physical equations,these models employ numerical techniques to integrate surface and subsurface hydrodynamics,thus replicating the macro-scale hydrological responses of watersheds.Numerical hydrological models have emerged as a leading and predominant trend in hydrological modeling due to their explicit representation of physical processes,heightened by their spatiotemporal resolution and reliance on interdisciplinary integration.This article focuses on the theoretical foundation of surface-subsurface numerical hydrological models.It includes a comparative and analytical discussion of leading numerical hydrological models,encompassing model architecture,numerical solution strategies,spatial representation,and coupling algorithms.Additionally,this paper contrasts these models with traditional hydrological models,thereby delineating the relative merits,drawbacks,and future directions of numerical hydrological modeling.

Integrated Hydrological Modeling of the Godavari River Basin in Maharashtra Using the SWAT Model: Streamflow Simulation and Analysis

Hydrological modeling plays a crucial role in efficiently managing water resources and understanding the hydrologic behavior of watersheds. This study aims to simulate daily streamflow in the Godavari River Basin in Maharashtra using the Soil and Water Assessment Tool (SWAT). SWAT is a process-based hydrological model used to predict water balance components, sediment levels, and nutrient contamination. In this research, we used integrated remote sensing and GIS data, including Digital Elevation Models (DEM), land use and land cover (LULC) maps, soil maps, and observed precipitation and temperature data, as input for developing the SWAT model to assess surface runoff in this large river basin. The Godavari River Basin under study was divided into 25 sub-basins, comprising 151 hydrological response units categorized by unique land cover, soil, and slope characteristics using the SWAT model. The model was calibrated and validated against observed runoff data for two time periods: 2003-2006 and 2007-2010 respectively. Model performance was assessed using the Nash-Sutcliffe efficiency (NSE) and the coefficient of determination (R2). The results show the effectiveness of the SWAT2012 model, with R2 value of 0.84 during calibration and 0.86 during validation. NSE values also ranged from 0.84 during calibration to 0.85 during validation. These findings enhance our understanding of surface runoff dynamics in the Godavari River Basin under study and highlight the suit-ability of the SWAT model for this region.

Flood Forecasting and Warning System: A Survey of Models and Their Applications in West Africa

Flood events occurrences and frequencies in the world are of immense worry for the stability of the economy and life safety. Africa continent is the third continent the most negatively affected by the flood events after Asia and Europe. Eastern Africa is the most hit in Africa. However, Africa continent is at the early stage in term of flood forecasting models development and implementation. Very few hydrological models for flood forecasting are available and implemented in Africa for the flood mitigation. And for the majority of the cases, they need to be improved because of the time evolution. Flash flood in Bamako (Mali) has been putting both human life and the economy in jeopardy. Studying this phenomenon, as to propose applicable solutions for its alleviation in Bamako is a great concern. Therefore, it is of upmost importance to know the existing scientific works related to this situation in Mali and elsewhere. The main aim was to point out the various solutions implemented by various local and international institutions, in order to fight against the flood events. Two types of methods are used for the flood events adaptation: the structural and non-structural methods. The structural methods are essentially based on the implementation of the structures like the dams, dykes, levees, etc. The problem of these methods is that they may reduce the volume of water that will inundate the area but are not efficient for the prediction of the coming floods and cannot alert the population with any lead time in advance. The non-structural methods are the one allowing to perform the prediction with acceptable lead time. They used the hydrological rainfall-runoff models and are the widely methods used for the flood adaptation. This review is more accentuated on the various types non-structural methods and their application in African countries in general and West African countries in particular with their strengths and weaknesses. Hydrologiska Byråns Vattenbalansavdelning (HBV), Hydrologic Engineer Center Hydrologic Model System (HEC-HMS) and Soil and Water Assessment Tool (SWAT) are the hydrological models that are the most widely used in West Africa for the purpose of flood forecasting. The easily way of calibration and the weak number of input data make these models appropriate for the West Africa region where the data are scarce and often with bad quality. These models when implemented and applied, can predict the coming floods, allow the population to adapt and mitigate the flood events and reduce considerably the impacts of floods especially in terms of loss of life.

Evapotranspiration estimation methods in hydrological models

Actual evapotranspiration is a key process of hydrological cycle and a sole term that links land surface water balance and land surface energy balance. Evapotranspiration plays a key role in simulating hydrological effect of climate change, and a review of evapotranspiration estimation methods in hydrological models is of vital importance. This paper firstly summarizes the evapotranspiration estimation methods applied in hydrological models and then classifies them into the integrated converting methods and the classification gathering methods by their mechanism. Integrated converting methods are usually used in hydrological models and two differences exist among them： one is in the potential evaporation estimation methods, while the other in the function for defining relationship between potential evapora tion and actual evapotranspiration. Due to the higher information requirements of the Pen- man-Monteith method and the existing data uncertainty, simplified empirical methods for calculating potential and actual evapotranspiration are widely used in hydrological models. Different evapotranspiration calculation methods are used depending on the complexity of the hydrological model, and importance and difficulty in the selection of the most suitable evapotranspiration methods is discussed. Finally, this paper points out the prospective de velopment trends of the evapotranspiration estimating methods in hydrological modeling.

Modelling Hydrological Consequences of Climate Change—Progress and Challenges

The simulation of hydrological consequences of climate change has received increasing attention from the hydrology and land-surface modelling communities. There have been many studies of climate-change effects on hydrology and water resources which usually consist of three steps： （1） use of general circulation models （GCMs） to provide future global climate scenarios under the effect of increasing greenhouse gases, （2） use of downscaling techniques （both nested regional climate models, RCMs, and statistical methods） for ＂downscaling＂ the GCM output to the scales compatible with hydrological models, and （3） use of hydrologic models to simulate the effects of climate change on hydrological regimes at various scales. Great progress has been achieved in all three steps during the past few years, however, large uncertainties still exist in every stage of such study. This paper first reviews the present achievements in this field and then discusses the challenges for future studies of the hydrological impacts of climate change.

Multi-Model Approach for Assessing the Influence of Calibration Criteria on the Water Balance in Ouémé Basin

Hydrological models are very useful tools for evaluating water resources, and the hydroclimatic hazards associated with the water cycle. However, their calibration and validation require the use of performance criteria which choice is not straightforward. This paper aims to evaluate the influence of the performance criteria on water balance components and water extremes using two global rainfall-runoff models (HBV and GR4J) over the Ouémé watershed at the Bonou and Savè outlets. Three (3) Efficacy criteria (Nash, coefficient of determination, and KGE) were considered for calibration and validation. The results show that the Nash criterion provides a good assessment of the simulation of the different parts of the hydrograph. KGE is better for simulating peak flows and water balance elements than other efficiency criteria. This study could serve as a basis for the choice of performance criteria in hydrological modelling.

GIS Application in Urban Flood Risk Analysis: Midar as a Case Study

The significance of this study lies in its exploration of the advanced applications of Geographic Information Systems (GIS) in assessing urban flood risks, with a specific focus on Midar, Morocco. This research is pivotal as it showcases that GIS technology is not just a tool for mapping, but a critical component in urban planning and emergency management strategies. By meticulously identifying and mapping flood-prone areas in Midar, the study provides invaluable insights into the potential vulnerabilities of urban landscapes to flooding. Moreover, this research demonstrates the practical utility of GIS in mitigating material losses, a significant concern in flood-prone urban areas. The proactive approach proposed in this study, centered around the use of GIS, aims to safeguard Midar’s population and infrastructure from the devastating impacts of floods. This approach serves as a model for other urban areas facing similar challenges, highlighting the indispensable role of GIS in disaster preparedness and response. Overall, the study underscores the transformative potential of GIS in enhancing urban resilience, making it a crucial tool in the fight against natural disasters like floods.

Effect of calibration data series length on performance and optimal parameters of hydrological model

In order to assess the effects of calibration data series length on the performance and optimal parameter values of a hydrological model in ungauged or data-limited catchments （data are non-continuous and fragmental in some catchments）, we used non-continuous calibration periods for more independent streamflow data for SIMHYD （simple hydrology） model calibration. Nash-Sutcliffe efficiency and percentage water balance error were used as performance measures. The particle swarm optimization （PSO） method was used to calibrate the rainfall-runoff models. Different lengths of data series ranging from one year to ten years, randomly sampled, were used to study the impact of calibration data series length. Fifty-five relatively unimpaired catchments located all over Australia with daily precipitation, potential evapotranspiration, and streamflow data were tested to obtain more general conclusions. The results show that longer calibration data series do not necessarily result in better model performance. In general, eight years of data are sufficient to obtain steady estimates of model performance and parameters for the SIMHYD model. It is also shown that most humid catchments require fewer calibration data to obtain a good performance and stable parameter values. The model performs better in humid and semi-humid catchments than in arid catchments. Our results may have useful and interesting implications for the efficiency of using limited observation data for hydrological model calibration in different climates.

Hydrological behaviour and water balance analysis for Xitiaoxi catchment of Taihu Basin

With the rapid social and economic development of the Taihu region, Taihu Lake now faces an increasingly severe eutrophication problem. Pollution from surrounding catchments contributes greatly to the eutrophication of water bodies in the region. Investigation of surface flow and associated mass transport for the Xitiaoxi catchment is of a significant degree of importance as the Xitiaoxi catchment is one of the major catchments within the Taihu region. A SWAT-based distributed hydrological model was established for the Xitiaoxi catchment. The model was calibrated and verified using hydrometeorological data from 1988 to 2001. The results indicate that the modeled daily and annual stream flow match the observed data both in the calibration period and the verification period, with a linear regression coefficient R2 and a coefficient e for modeled daily stream flow greater than 0.8 at Hengtangcun and Fanjiacun gauge stations. The results show that the runoff process in the Xitiaoxi catchment is affected both by rainfall and human activities （e.g., reservoirs and polder areas）. Moreover, the human activities weaken flood peaks more noticeably during rainstorms. The Water balance analysis reveals the percentages of precipitation made up by surface flow, evapotranspiration, groundwater recharge＇ and the change of soil storage, all of which are considered useful to the further understanding of the hydrological processes in the Xitiaoxi catchment. This study provides a good base for further studies in mass transport modeling and comparison of modeling results from similar hydrological models.

Spatial distribution characteristics and mechanism of nonhydrological time-variable gravity in China's Mainland

The purpose of this study is to explore nonhydrological mass transfer in China's Mainland.For this purpose,gravity recovery and climate experiment(GRACE)data were obtained to study the spatial distribution of time variant gravity signals in China's Mainland.Then,from auxiliary hydrological data processed according to the current hydrological model,a new more comprehensive hydrological model of China's Mainland was constructed.Finally,the time variant signals of this new hydrological model were removed from the time variant gravity field computed from GRACE data,thus obtaining a description of the nonhydrological mass transfer of China's Mainland.The physical sources and mechanisms of the resulting mass transfer are then discussed.The improved,more realistic,hydrological model used here was created by selecting the hydrological components with the best correlations in existing hydrological models,by use of correlation calculation,analysis,and comparison.This improved model includes water in soils and deeper strata,in the vegetation canopy,in lakes,snow,and glaciers,and in other water components(mainly reservoir storage,swamps,and rivers).The spatial distribution of the transfer signals due to nonhydrological mass in China's Mainland was obtained by subtracting the combined hydrological model from the GRACE time-variable gravity field.The results show that the nonhydrological signals in China's Mainland collected in GRACE data were mainly positive signals,and were distributed in the Bohai Rim and the northern and eastern parts of the Tibetan Plateau.The above nonhydrological mass transfer signals have been studied further and are discussed.The results show that the nonhydrological mass migration signals in the Bohai Rim region originate primarily from sea level change and marine sediment accumulation.The mass accumulation from Indian plate collision in the Tibetan Plateau appears to be the main reason for the increase in the residual gravity field in that region.

Hydrological characteristics and changes in the Nu-Salween River basin revealed with model-based reconstructed data

The Nu-Salween River(NSR),the longest free-flow river in Southeast Asia,plays an irreplaceable role in social development and ecological protection.The lower NSR region is particularly valuable as it is inhabited by approximately 6.7 million people.The basin has limited hydraulic conservancy infrastructure and insufficient ability to cope with climate change risks.Studying the hydrological characteristics and changes in the basin provides the scientific basis for rational protection and development of the basin.However,owing to the limitation of observation data,previous studies have focused on the local area and neglected the study of the lower reaches,which is not enough to reflect the spatial characteristics of the entire basin.In this study,the ECMWF 5th generation reanalysis data(ERA5)and Multi-Source Weighted-Ensemble Precipitation(MSWEP)were applied to develop a geomorphology-based hydrological model(GBHM)for reconstructing hydrological datasets(i.e.GBHM-ERA5 and GBHM-MSWEP).The reconstructed datasets covering the complete basin were verified against the gauge observation and compared with other commonly used streamflow products,including Global Flood Awareness System v2.1,GloFAS-Reanalysis dataset v3.0,and linear optimal runoff aggregate(LORA).The comparison results revealed that GBHM-ERA5 is significantly better than the other four datasets and provides a good reproduction of the hydrological characteristics and trends of the NSR.Detailed analysis of GBHM-ERA5 revealed that:(1)A multi-year mean surface runoff represented 39%of precipitation over the basin during 1980–2018,which had low surface runoff in the upstream,while areas around the Three Parallel Rivers Area and the estuary had abundant surface runoff.(2)The surface runoff and discharge coefficient of variations in spring were larger than those in other seasons,and the inter-annual variation in the downstream was smaller than that in the upstream and midstream regions.(3)More than 70%of the basin areas showed a decreasing trend in the surface runoff,except for parts of Nagqu,south of Shan State in Myanmar,and Thailand,where surface runoff has an increasing trend.(4)The downstream discharge has dropped significantly at a rate of approximately 680 million cubic metresper year,and the decline rate is greater than that of upstream and midstream,especially in summer.This study provides a data basis for subsequent studies in the NSR basin and further elucidates the impact of climate change on the basin,which is beneficial to river planning and promotes international cooperation on the water-and eco-security of the basin.

Light-absorbing Particles in Snow and Ice: Measurement and Modeling of Climatic and Hydrological impact

Light absorbing particles（LAP, e.g., black carbon, brown carbon, and dust） influence water and energy budgets of the atmosphere and snowpack in multiple ways. In addition to their effects associated with atmospheric heating by absorption of solar radiation and interactions with clouds, LAP in snow on land and ice can reduce the surface reflectance（a.k.a., surface darkening）, which is likely to accelerate the snow aging process and further reduces snow albedo and increases the speed of snowpack melt. LAP in snow and ice（LAPSI） has been identified as one of major forcings affecting climate change, e.g.in the fourth and fifth assessment reports of IPCC. However, the uncertainty level in quantifying this effect remains very high. In this review paper, we document various technical methods of measuring LAPSI and review the progress made in measuring the LAPSI in Arctic, Tibetan Plateau and other mid-latitude regions. We also report the progress in modeling the mass concentrations, albedo reduction, radiative forcing, and climatic and hydrological impact of LAPSI at global and regional scales. Finally we identify some research needs for reducing the uncertainties in the impact of LAPSI on global and regional climate and the hydrological cycle.

Efficiency and Feasibility of an Integrated Algorithm for Distributed Hydrological M odel Calibration

Increasing complexity of distributed hydrological model(DHM)has lowered the efficiency of convergence.In this study,global sensitivity analysis(SA)was introduced by combining multiobjective(MO)optimization for DHM calibration.Latin Hypercube-once at a time(LH-OAT)was adopted in global parameter SA to obtain relative sensitivity of model parameter,which can be categorized into different sensitivity levels.Two comparative study cases were conducted to present the efficiency and feasibility by combining SA with MO(SA-MO).WetSpa model with non-dominated sorting genetic algorithm-Ⅱ(NSGA-Ⅱ)algorithm and EasyDHM model with multi-objective sequential complex evolutionary metropolis-uncertainty analysis(MOSCEM-UA)algorithm were adopted to demonstrate the general feasibility of combining SA in optimization.Results showed that the LH-OAT was globally effective in selecting high sensitivity parameters.It proves that using parameter from high sensitivity groups results in higher convergence efficiency.Study case I showed a better Pareto front distribution and convergence compared with model calibration without SA.Study case II indicated a more efficient convergence of parameters in sequential evolution of MOSCEM-UA under the same iteration.It indicates that SA-MO is feasible and efficient for high dimensional DHM calibration.

Hydrological Processes and Climate Change in Hydrographic Regions of Brazil

The objective of this work is to assess the impacts of IPCC AR5 climate change scenarios on water resources and hydrological processes across the entire Brazilian territory. Hydrological simulations are carried out in total drainage area of about 11,535,645 km2 and average stream flow of about 272,460 m3/s. The study area consists of different climates and land covers such as the Amazon Forest, Northeast Semiarid, Brazilian Savannah, Pantanal wetlands and temperate climate in the South. The atmospheric forcing to drive the large-scale hydrological model MGB-IPH is derived from the downscaling of two global climate models, HadGEM2-ES and MIROC5, by the Eta Regional Climate Model, at 20 km resolution. The Eta model provided the downscaling of the baseline (1961-1990) and three time-slices (2011-2040, 2041-2070 and 2071-2099). These projections adopted two emission scenarios, the RCP 4.5 and RCP 8.5. The change in the average and extremes of precipitation, evapotranspiration, rates of river discharge and soil moisture were assessed. The simulations showed the response of the hydrographic regions due to change of precipitation and potential evapotranspiration in the scenarios. Water availability decreases in almost the entire study area (exception for the South) and the major basins for hydroelectric power generation are affected. The Northwest, Amazon and a small area along the Northeast Atlantic coast exhibited intensification of the extremes discharges, where the anomaly is positive for high-flow (Q10) and negative for low-flow (Q95). The results highlight the most climatic sensitive regions in Brazil in terms of hydrological variables and water resources.

Contribution of Agro-Hydrological Modeling in the Evaluation of Water Availability of an Ungauged Basin Reservoir in Côte d’Ivoire: Case of the Loka Reservoir in Bouaké

The city of Bouaké, the second biggest city of Côte d’Ivoire, experienced a water shortage in 2018 that lasted four months due to the drying up of the Loka reservoir, which supplies two-thirds of the city. The challenge of the Loka reservoir is that it is located in an ungauged basin where very few hydrological studies have been carried out, despite the recurrent problems of access to drinking water. In the purpose to better understand the phenomena that caused this temporary drying of the dam, the methodology implemented was based on agro-hydrological modeling with SWAT using a regionalization technique of a nearby watershed. The model performance was assessed using three statistical indices (the Nash-Sutcliffe coefficient (NS), the coefficient of determination (R2) and the percentage of bias (PBIAS)) and the visual appreciation of hydrographs for monthly series. The statistical indices appear satisfactory with a NS and R2 ≥ 0.6 both for calibration and validation, and a PBIAS of -11.2 and -3.8 respectively for calibration and validation. The hydrological modeling of Loka basin has shown the impact of climate change already reported by some authors as well as anthropization. Thus, while the reservoir records a decrease in its water volume estimated at 384,604 m3 each year, the water demand undergoes an increase of 122,033 m3 per year.

Parameter Estimation of a Distributed Hydrological Model Using a Genetic Algorithm

Water is a vital resource, and can also sometimes be a destructive force. As such, it is important to manage this resource. The prediction of stream flows is an important component of this management. Hydrological models are very useful in accomplishing this task. The objective of this study is to develop and apply an optimization method useful for calibrating a deterministic model of the daily flows of the Miramichi River watershed (New Brunswick). The model used is the CEQUEAU model. The model is calibrated by applying a genetic algorithm. The Nash-Sutcliffe efficiency criterion, modified to penalize physically unrealistic results, was used as the objective function. The model was calibrated using flow data (1975-2000) from a gauging station on the Southwest Miramichi River (catchment area of 5050 km2), obtaining a Nash-Sutcliffe criterion of 0.83. Model validation was performed using flow data (2001-2009) from the same station (Nash-Sutcliffe criterion value of 0.80). This suggests that the model calibration is sufficiently robust to be used for future predictions. A second model validation was performed using data from three other measuring stations on the same watershed. The model performed well in all three additional locations (Nash-Sutcliffe criterion values of 0.77, 0.76 and 0.74), but was performing less well when applied to smaller sub-basins. Nonetheless, the relatively strong performance of the model suggests that it could be used to predict flows anywhere in the watershed, but caution is suggested for applications in small sub-basins. The performance of the CEQUEAU model was also compared to a simple benchmark model (average of each calendar day). A sensitivity analysis was also performed.

Comparison of two different methods for determining flow direction in catchment hydrological modeling

Digital elevation models （DEMs） are widely used to define the flow direction in distributed hydrological models for simulation of streamflow. In recent decades, numerous methods for flow direction determination have been applied successfully to mountainous regions. Nevertheless, some problems still exist when those methods are used for flat or gently sloped areas The present study reviews the conventional methods of determining flow direction for such landscapes and analyzes the problems of these methods. Two different methods of determining flow direction are discussed and were applied to the Xitiaoxi Catchment, located in the Taihu Basin in southern China, which has both mountainous and flat terrain. Both the agree method and the shortest path method use drainage networks derived from a remote sensing image to determine the correct location of the stream. The results indicate that the agree method provides a better fit with the DEM for the hilly region than the shortest path method. For the flat region where the flow has been diverted and rerouted by land managers, both methods require observation of the drainage network to determine the flow direction. In order to clarify the applicability of the two methods, both are employed in catchment hydrological models conceptually based on the Xinanjiang model and implemented with PCRaster. The simulation results show that both methods can be successfully applied in hydrological modeling. There are no evident differences in the modeled discharge when using the two methods at different spatial scales.

Hydrological assessment of TRMM rainfall data over Yangtze River Basin

High-quality rainfall information is critical for accurate simulation of runoff and water cycle processes on the land surface. In situ monitoring of rainfall has a very limited utility at the regional and global scale because of the high temporal and spatial variability of rainfall. As a step toward overcoming this problem, microwave remote sensing observations can be used to retrieve the temporal and spatial rainfall coverage because of their global availability and frequency of measurement. This paper addresses the question of whether remote sensing rainfall estimates over a catchment can be used for water balance computations in the distributed hydrological model. The TRMM 3B42V6 rainfall product was introduced into the hydrological cycle simulation of the Yangtze River Basin in South China. A tool was developed to interpolate the rain gauge observations at the same temporal and spatial resolution as the TRMM data and then evaluate the precision of TRMM 3B42V6 data from 1998 to 2006. It shows that the TRMM 3B42V6 rainfall product was reliable and had good precision in application to the Yangtze River Basin. The TRMM 3B42V6 data slightly overestimated rainfall during the wet season and underestimated rainfall during the dry season in the Yangtze River Basin. Results suggest that the TRMM 3B42V6 rainfall product can be used as an alternative data source for large-scale distributed hydrological models.

Simulating streamflow and water table depth with a coupled hydrological model

A coupled model integrating MODFLOW and TOPNET with the models interacting through the exchange of recharge and baseflow and river-aquifer interactions was developed and applied to the Big Darby Watershed in Ohio, USA. Calibration and validation results show that there is generally good agreement between measured streamflow and simulated results from the coupled model. At two gauging stations, average goodness of fit （ R2 ）, percent bias （ PB ）, and Nash Sutcliffe efficiency （ENS） values of 0.83, 11.15%, and 0.83, respectively, were obtained for simulation of streamflow during calibration, and values of 0.84, 8.75%, and 0.85, respectively, were obtained for validation. The simulated water table depths yielded average R2 valuesof0.77 and 0.76 for calibration and validation, respectively. The good match between measured and simulated streamflows and water table depths demonstrates that the model is capable of adequately simulating streamflows and water table depths in the watershed and also capturing the influence of spatial and temporal variation in recharge.