1. Introduction
The Curve Number (CN) is one of the most frequently used methods for predicting runoff, and determining rainfall–runoff relationships [
1,
2]. Curve Numbers represent average antecedent runoff conditions for urban, cultivated agricultural, and arid and semi-arid range land uses. Curve Numbers can be determined from soil type, land use, soil moisture conditions, and hydrologic conditions. In hydrology, CN values are used to determine how much rainfall is lost (initial abstractions and infiltration) and how much becomes runoff. A high CN means high runoff and low losses as with urban areas, whereas a low CN means low runoff and high infiltration. The Soil Conservation Service (SCS) CN method is one of the most popular methods for computing the runoff volume from a rainstorm. It is popular because it is simple, easy to understand and apply, stable, and accounts for most of the runoff producing watershed characteristics. The SCS CN is a function of land use and hydrologic soil groups and is calculated using a table relating the two [
3,
4]. For each land use type, there is a corresponding CN for each of the SCS hydrologic soil types A, B, C, and D, with values ranging theoretically from 0–100 [
2]. Selection of the CN is generally accomplished using lookup tables and site-specific information based on soil groups, land use, vegetative cover, and antecedent moisture. Previous studies have focused on CN identification and utilization, and have suggested adjusting CNs for watershed conservation [
5,
6,
7,
8].
Several models exist that can calculate watershed runoff [
9]. Each component of a model is an aspect of the rainfall–runoff process, such as precipitation, losses, runoff transformation, and river storage and routing. Representing one of these components requires a set of parameters, which specify the particular characteristics of the watershed in terms of their mathematical and empirical relation to physical processes. Two such models that can predict watershed runoff that are used in this study are the HEC-1 Model, and the Rainfall Runoff Model (RRM).
Surface water resources in Jordan are mainly used for agriculture because they are mixed with treated water and cannot be used for domestic purposes. Groundwater is also used for domestic purposes but it is not enough to satisfy the growing demand for domestic and industrial sectors [
10]. Surface water resources are distributed among 15 major basins in Jordan, with the Zarqa River Basin (ZRB) considered as one of the most significant basins due to its economic, social and agricultural importance as the second main tributary to the Jordan River. The Zarqa River Basin has been taken as a case study because it suffers from several water problems ranging from poor water management, water quality, and from water conflicts over allocation amongst various sectors [
10,
11]. It is considered to be the area most affected by water scarcity because of its limited water resources and high demand, which is due to the fact that the basin houses about 52% of Jordan’s population [
12,
13]. The heavy utilization of the Zarqa Basin has resulted in the reduction of the base flow of the Zarqa River from 5 m
3/s to less than 1 m
3/s and the discharge of springs from an average 317 million cubic meter (MCM) prior to 1985 to less than 130 MCM after 2000.
The Zarqa Basin drains 4120 square kilometers and extends from Syria to Amman, and then to the Jordan River. The highest point in the basin is located near the Syrian town of Salkhad in Jebel Al-Arab at an elevation of 1460 m. The basin discharges its water at the confluence of the Zarqa River with the Jordan River at an elevation of −350 m below sea level. The steep bed slopes of the river and its tributaries are in the range of 0.7 to 1.5 percent on average. The stream flow of the Zarqa River is impounded by the King Talal Dam (KTD), south-west of the town of Jaresh at an altitude of about 120 m above sea level, which has a storage capacity of 75 MCM, the largest dam in Jordan. The flood (runoff) gauge station of New Jerash Road Bridge is the only gauge in the Zarqa Basin and has a database that is long enough to be suitable for evaluation. The station near Ain Ghazal, important for monitoring the runoff from Amman, stopped operating in 1979. The Zarqa River Basin spans five governorates: Amman, Balqa’a, Zarqa, Jerash, and Mafraq [
14].
Figure 1 shows the location of the Zarqa River Basin and the main cities it contains.
The soils of the ZRB differ widely according to various slopes. In the most humid west, reddish to brownish clay and clay loams prevail. Toward the east, the soil becomes silty loam to loamy in texture, yellowish brown to strong brown with very high carbonate content. Soil erosion on the steeper slope is the main cause of alluvial processes where sheet erosion is more pronounced in the eastern part [
15]. These changes have taken place in a matter of decades starting in the early 1960s.
Land use in the ZRB has undergone considerable changes [
15]. The recent expansion of Amman and surrounding towns has been extensive; the area where large expanses of grazing land and more fertile agricultural lands could once be found between Amman and other towns has now developed into a sprawling urban conglomerate. The Basin has regions of natural forests occurring in the mountainous area that are composed of oak, pine, juniper, wild olive and cypress trees. Agriculture is scattered within the Basin comprising rain-fed orchards, olive, and field crops, with irrigated agriculture on the river basin banks. The major industrial activities in the Basin include the Hussein thermal plant, petrol refining, the textile industry, paper processing, leather production, and intermediate petrochemical processing. The ZRB which was once full of scattered forest, rangelands and a haven for wild and domestic animals, has suffered from desertification caused by human activities and climate change.
Therefore, the critical situation of the Zarqa Basin requires proper management to enhance watershed restoration efforts. For this purpose, the identification of CN values based on existing land uses and hydrologic soil groups (HSGs) classification is an important first step for developing a model dealing with watershed management problems. This paper presents possible solutions for watershed restoration by using the HEC-1 and the RRM models, with the aim of estimating the baseline CN values and proposing new CN values that would represent restoration scenarios. The prediction of CNs for the average rainy season of 2001/2002 will be in parallel with the prediction of runoff using the HEC-1 and the RRM models. The HEC-1 Model was developed by the US Army Corps of Engineers [
16], while the RRM was implemented as a web-based distributed client-server system and was provided by the Environmental Software Systems-ESS, Austria [
17,
18].
2. Methodology
The data used for running these models includes elevation maps; precipitation amounts; temperature; flow; runoff; satellite imagery; soil maps, and basin characteristics. This data is used for the 2001/2002 water year from 1 October 2001 to 30 September 2002. The water year of 2001/2002 is taken because it was an average rainy season with a good distribution of rainfall amounts. For that year, all data related to daily rainfall and runoff records was available. The gauging station of New Jerash Road Bridge is the most essential gauge with a complete record and is used for the HEC-1 and the RRM model calibrations. This allows for the simulated flow from the HEC-1 Model and the RRM to be compared with the actual data to calibrate parameters and verify the accuracy of the model.
2.1. Watershed Modeling System (WMS): HEC-1 Model
The Watershed Modeling System (WMS) is a comprehensive graphical modeling environment for all phases of watershed hydrology and hydraulics. The WMS includes tools to automate model parameterization such as automated basin delineation, geometric parameter calculations, GIS overlay computations, and cross-section extraction from terrain data. The WMS software supports many hydrologic models, including the HEC-1 [
19], and can be used for the estimation of the similar parameters used by the RRM.
The HEC-1 Model is one of the most widely-used lumped parameter models available. It is designed to simulate the surface runoff response of a river basin to precipitation by representing the basin as an interconnected system of hydrologic and hydraulic components. It simulates a single storm event, and includes several options for modeling and computing total runoff volume for a rain storm event. Precipitation loss is calculated based on the CN values and on the initial surface moisture storage capacity. The CN and surface moisture storage capacity are related to the total runoff depth for a storm. The certainty of the HEC-1 application was deduced from the compatibility of the modeled flow to the actual flow data. For instance, the availability of daily runoff data allows for the calibration of the modeled flow with a known data series. The limitation of the HEC-1 Model is its simulation to a single storm due to the fact that provisions are not made for soil moisture recovery during periods of no precipitation. Therefore, the total seasonal runoff volume was considered in simulating the watershed.
2.2. WaterWare System: Rainfall–Runoff Model (RRM)
The WaterWare system is a model-based information and decision support system for water resources management. The system is designed to support the implementation of the Water Framework Directive or any similar national legislation. The system is a modular client-server system that smoothly integrates several databases, GIS, and a range of dynamic simulation and optimization models for a wide range of Integrated Water Resources Management (IWRM) tasks. The system is implemented in a fully web-enabled environment, supporting the seamless integration of databases, simulation and optimization models, and analytical tools into a common, easy-to-use framework. This includes a multi-media user interface with Internet access, a GIS with hierarchical map layers, object databases, time series analysis, reporting functions, and an embedded expert system for estimation, classification and impact assessment tasks. Monitoring data acquisition, automatic forecast runs and reporting, is driven by a real-time rule-based expert system. The components of the WaterWare System can be summarized as follows: a- Database management for all components of a basin including monitoring time series. b- Embedded GIS functionality. c- A set of dynamic simulation models in the base system. The WaterWare system provides the following: a- a scenario analysis; b- continuous monitoring, now-casting and scheduled forecasting; c- event-based operation triggered by user-defined events.
The WaterWare on-line system and databases include the following: the Water Resources Model (WRM), the Rainfall–Runoff Model (RRM), the Irrigation Water Demand Model, the Water Quality Model, the WRM Optimization Model, the Water Technologies Database, the Discrete Multi Criteria Decision Support System (DSS), the River Basin Object Database, the Water Institutions Database; and the Water Issues Questionnaire [
18].
The WaterWare system was developed through a series of applications, where different institutions from several Mediterranean countries (Austria, Italy, Greece, Cyprus, Malta, Turkey, Lebanon, Jordan, Palestine, Tunisia, and Morocco) were involved in the OPTIMA Project (Optimization for Sustainable Water Management), sponsored by the EU under the Sixth Framework Programme (FP6) with contract No. INCO-CT-2004-50909 [
20,
21]. The Rainfall–Runoff Model is one of the components of the WaterWare system [
18]. It is a dynamic, spatially lumped, basin scale water budget model that operates at a daily time step. It is designed for un-gauged basins, estimating their outflow as a function of basin characteristics including elevation distribution, land use distribution, rainfall in mm and average air temperature in degree Centigrade as dynamic inputs.
Editing a scenario for the Rainfall–Runoff Model involves major components:
Basic basin characteristics data such as basin size, and channel characteristics;
Defining basin morphometry, and geometry
Providing specific model parameters including the initial conditions
Land use data such as forests, pastures, agricultural practices, built up areas
Temperature and precipitation data as a daily time series
Daily runoff model results data shown together with the observed runoff data that can also be used for model calibration
Daily main output data from the Rainfall–Runoff Model as well as a stand-alone model for analysis, or generation of the RRM data in a file format that is compatible with the Water Resources Model (WRM) in order to be used as input for the WRM
2.3. Soil and Land Use Classification for Predicting Curve Numbers
Soil and land use classifications are used for predicting the baseline CN values of the Zarqa River Basin. This is done by overlaying the classification of Hydrologic Soil Groups with the land use maps of the Zarqa River Basin. The identification of the existing CN for the average rainy year will be in parallel with the prediction of runoff using the HEC-1 Model, and the Rainfall–Runoff Model (RRM).
These classifications are based on USDA Soil Conservation Service [
2], which considers factors such as soil texture, soil depth, infiltration rate, and layers that impede the downward movement of water. Four hydrologic soil groups A, B, C and D are considered as the basic classification of soils of the watershed as defined in
Table 1. The soil textures of the Zarqa River Basin differ according to slope and soil depth in terms of percentages of silty clay loam, clay, silty clay, clay loam, sandy clay loam, and others. The properties of the representative soil profiles of each soil map unit and the description of each soil map unit are classified into their hydrologic soil groups.
The land use classifications from the September 2001 Landsat Enhanced Thematic Mapper (ETM) imagery of the Zarqa Basin include 6 different classes established by using Environment for Visualizing Images (ENVI) software for satellite data analysis. The classes are represented in percentages.
The CN values for each unique area are found by identifying the soil type and assigned land use in the lookup tables of the Soil Conservation Service (SCS) [
2]. An example of the values which are considered as a reference for determining the Zarqa Basin CNs is shown in
Table 2. The weighted value of CN for the Zarqa Basin is determined according to the Equation (1).
CN = Weighted curve number.
CNi = Curve number for each unique soil and land use.
Ai = Area with curve number CNi, and A = Total area.
The CN for the average antecedent moisture condition AMC II can be converted to the dry condition AMC I and the wet condition AMC III using Equations (2) and (3) [
3]:
3. Results and Discussion
3.1. HEC-1 Model, and Rainfall Runoff Model
Prior to running the HEC-1 model, and the RRM, all parameters are checked for consistency. The profile of each soil map unit is classified and results in two Hydrologic Soil Groups (HSGs), HSG B and HSG C, whereas, the land use (existing condition of the land) of the September 2001 Landsat Enhanced Thematic Mapper (ETM) imagery of the Zarqa Basin has been classified into six different classes. The different land use classes are water body, forest and orchard, mixture of grass, weeds and desert shrubs, crops, urban areas, and the bare soil.
Table 3 shows the land use classification of the Zarqa River Basin in percentages. The Zarqa River Basin hydrologic soil group map and the land use classification map are overlaid to determine the CNs of the entire watershed using the HEC-1 Model.
Figure 2 Shows the Zarqa River Basin’s soil and land use classification maps.
The standard SCS Table [
2] is used to estimate the initial CN values, then, the HEC-1 Model runs the simulations and adjusts the CN values until the measured runoff is computed within a close amount.
Table 4 shows the weighted CNs for the Zarqa Basin where the average weighted CNs for the Zarqa River Basin are 59 (dry), 78 (normal) and 89 (wet). During the growing season, the average curve numbers are 52, 72 and 86. These weighted baseline CNs are the keys to modeling/simulating the runoff at the Zarqa Basin’s outlet; Runoff Gauge at Jerash Bridge, and the expected CN adjustments for restoration scenarios.
The HEC-1 Model is run using the Watershed Modeling System (WMS) [
16,
17,
18,
19]. The results of running the HEC-1 Model show that the annual calibrated (modeled) runoff volume is 36.4 MCM from the total rainfall storms, which totaled 911 MCM. The measured (actual or observed) total runoff volume of the Zarqa Basin outlet was 36.6 MCM. Thus, the simulated runoff with a value of 36.4 MCM almost matches the figures of the gauged values from the Ministry of Water and Irrigation. It should be noted that this excludes the base flow. The base flow consists of the effluent of the wastewater treatment plant of khirbet es-Samra which is about 2.0 m
3/s. The calibrated runoff value of 36.4 MCM is achieved after adjusting the CN values (model calibration), and shows that the land uses are 0.3% water body, 9.3% forest and orchard, 71% mixture of grass, weeds, and desert shrubs, 7.0% crops, 4.0% urban areas, and 8.4% bare soil. Thus, the amount of runoff in the HEC-1 Model is a function of the parameters, primarily, the composite CN of the Zarqa River Basin for a specific year. Using the land use percentages of the HEC-1 Model with the RRM will also simulate the flow for the same year 2001/2002. The RRM incorporates the existing conditions by setting up the basin’s parameters as rainfall, temperature, geometry, morphometry, and land use distributions.
The Rainfall–Runoff Model (RRM) is run online using the huge capability of the Environmental Software and Services (ESS) server [
17,
18].
Figure 3 shows the Zarqa River Basin’s flow simulation. It also shows the time series of precipitation (mm/day), the average daily temperature (°C), and the measured and modeled flow for the studied year (1 October 2001 to 30 September 2002). The results of running the RRM show that the average modeled flow (runoff and base-flow) is about 2.97 m
3/s or about 93.7 MCM/year, which matches the figures of the Ministry of Water and Irrigation 2002, with an average value of 2.96 m
3/s or about 93.4 MCM/year. The simulated flow is highly compatible to the current (measured) flow due to the simulation capability of RRM, which has been set up based on the existing conditions of the basin, primarily the land use for that year.
Comparing the results of the HEC-1 Model in predicting the Zarqa River Basin runoff, with the RRM’s predicted water flow of the Zarqa River Basin proved their simulation efficiency. Nevertheless, the HEC-1 Model uses a single storm event and is based on values of curve number, while the RRM is a dynamic, spatially lumped model that operates at a daily time step, and estimates the flow as a function of the basin characteristics. The models’ predictions for the real runoff and flow definitely establish the modeling certainty and help water resource developers to incorporate different basin features for watershed simulation and restoration.
3.2. CN Applications for Watershed Restoration
Using the baseline result from
Section 3.1, the proposed changes in land use can be simulated via the HEC-1 Model and provide a mechanism to evaluate conservation mechanisms for Zarqa River Basin restoration planning. Two scenarios are developed here as tools for making decisions on the land class [
22] of mixtures of grass, weeds, and desert shrubs.
Table 5 shows the weighted CNs for the Zarqa Basin restoration, simulated using the HEC-1 Model.
First scenario: This scenario proposes the use of contours to lands which are mixtures of grass, weeds, and desert shrubs. The average weighted CNs for the Zarqa River Basin are 58, 76 and 88 under dry, normal and wet conditions respectively. During the growing season, the average CNs are 51, 71 and 85. The results of running the HEC-1 Model show that the annual average predicted runoff volume is 33.1 MCM from the total rainfall.
Second scenario: This scenario proposes the use of contours, terraces and crop residues to lands which are mixtures of grass, weeds, and desert shrubs. The average weighted CNs for the Zarqa River Basin are 57, 75 and 87. During the growing season, the average CNs are 50, 70 and 84. The results of running the HEC-1 Model show that the annual predicted runoff volume is 28.4 MCM from the total rainfall based on the average rainy season of 2001/2002.
Herein, the results of running the HEC-1 Model and the RRM as baseline scenarios for runoff prediction are compared. The HEC-1 Model shows that the average calibrated runoff is about 36.4 MCM (excluding the base-flow), whereas the RRM shows that the average modeled flow (runoff + base-flow) is about 2.97 m3/s; about 93.7 MCM/year. Land use change scenarios are affecting CNs as the predicted runoff decreases from 36.4 to 33.1 by applying contours, and then to 28.4 by applying contours, terraces and crop residues. The land use changes will help in enhancing the infiltration of water into the soil to increase soil moisture storage and recharge groundwater, while maintaining sufficient surface water downstream of the basin. This will lead to the gradual restoration of the watershed with time. Land use changes have been demonstrated in protected areas in the adjacent watershed. Thus, the models can present, simulate, predict the flow, and provide watershed restoration options for decision makers.
3.3. Sustainability of the Zarqa River Basin
The area of the basin before 1900 was thriving with scattered forests, bushes and pastures and was grazing land for domestic and wild animals such as camels, sheep, deer and Zebras. It was also a haven for migrating birds. Since the beginning of the 20th century, the area has suffered from gradual deforestation, overgrazing, land degradation and soil erosion resulting in bare lands which are subject to wind and water erosion. In addition, since the middle of last century, the basin has witnessed the expansion of urban areas associated with an increase in the construction of roads and other infrastructure. This further increases land degradation and loss of land cover. These conditions also enhance runoff and erosion and reduce water storage in the soil profile.
The suggested interventions of land restoration such as contouring, terracing and crop residues has changed land use and thus reduced CN values (
Table 5). Runoff has decreased from 36.4 MCM to 28.4 MCM using the above measures, indicating that an additional 8.0 MCM of water will be stored in the soil profile. Over the area of the basin of 4120 km
2, an additional 50 mm of water could help in supporting plant cover. This could be considered as an element of sustainability [
23] by maintaining the basin’s integrity, thus improving food production.