Kinematic Storage Model (KSM) for Groundwater Development in Highly Permeable Hill Slope-Laboratory Study

Study of the groundwater table development and runoff generation is one of the most important parts of hydrology to develop a clear concept, especially in hill slope. The study is more complex in the real field rather than in the artificial system. The result in artificial systems developed and experimental observations may give good results. So, therefore, this study is aimed at modeling in the laboratory as artificial hill slope flows which include saturation excess surface runoff flows. The physical processes along with runoff generation depend on the factors—soil type, characteristic slope geometry, and initial soil conditions at the commencement of rainfall. The mechanisms involved in runoff formation process have been simulated successfully to compute hydrograph for hilly terrain and groundwater table development in highly permeable soil tested by kinematic storage model theory. The model represents the hill slope as a rectangular storage element of length 2.02 m, depth 0.15 m and width 1m With an impermeable bed making an angle of 10 degrees with the horizontal. The storage element is composed of two moisture zones: an unsaturated zone and a saturated zone. The result obtained is seemed good adjustment to the theory of hill slope model given by Nm Shakya, 1995. Aslo, the moisture profile variation in mixed sand profile was found immediately after the rainfall event. The result obtained shows that the timing and distribution of moisture over the depth where the maximum moisture content is 0.4 in mid of the depth which is more than in surface having a moisture level of 0.37.


Introduction
Hill slope stability is an important state of hill to avoid the land slide due to intense rainfall. In infinite slope length with highly permeable soil has high infiltration rate having the maximum volumetric capacity to store the infiltrated water. In such cases, the saturated conductivity of soil is always greater than the rainfall rate so therefore the rainfall water infiltrates into soil completely. If the depth of soil is low then the saturation level increases toward the slope surface and a time comes to reach the surface to be saturated and overland flow occurs. Hill slope having infinite length, the infiltrated water has high active pressure so that it may cause the land slide from a weak crack or fault line or the fault may be developed due to high shear stress development [1]. The initial moisture condition has significant role in the rise in moisture level and time for saturation. High antecedent moisture level causes fast surface saturation resulting starts of surface runoff as soon. Before the start of surface runoff, the subsurface flow is prominent to contribute the discharge in the river at the end of hill slope, the rate of which can be given by Darcy equation as: where K is the hydraulic conductivity of soil in cm/se, I is the hydraulic gradient (slope of the water pore pressure ) and A is the section area of the slope.
The seepage flow through subsurface is laminar in nature but may cause piping action from the ground level and may lead to the land slide [2]. The groundwater table depends on soil slope and type of soil. The moisture profile and its time for saturation completely depend on the rainfall intensity and its duration. Green and Ampt model gives the total volume infiltration and it's time for saturation. Which depends on the initial moisture deficit (IMD), the average hydraulic potential of soil, rainfall intensity and hydraulic conductivity. After the soil saturation, further rainfall cause rise in infiltrated volume due to the hydraulic conductivity which is less than the saturation conductivity. The rapid rise in moisture level causes the surface runoff earlier in the hill [3]. The infiltration rate can be estimated in ideal natural soil by different models among which the Hortan's equation (1930s) is more accurate. The rise in moisture with depth depends on the particles size distributions and pore size i.e. micro pore and macro pore interconnections of soil particles and it affects the capillary sanction and potential in soil [1]. The flow is in steady-state condition that follows to the kinematic wave approximation theory with Darcy law and the cause of flow is pore water pressure and gravitational force and the flow can be modeled by dynamic wave approximation theory too [4]. The interconnection of hill slope and river is the main source of water in the river due to subsurface flow and it may be three-dimensional or may be one-dimensional flow that can be analyzed by three-dimensional dynamic equation or assuming one-dimensional kinematic wave equation [5].
The moisture variation in hill slope and estimation of surface runoff is a challenging issue in hilly areas in case of under developed country as the massive where water is sufficiently available [6]. The vegetation, type of crop and land use pattern, is also a important factor affecting the infiltration process especially in clayey type soil slope which reduces to the surface runoff and enhances groundwater recharge [7]. So the study of the subsurface flow and its direction is felt to be carried out in an easy way which gives reliable results. Also, the aim of the study is to give validation for the existing water development process in hill slope despite the length of the slope but on the basis of initial moisture condition, the soil type and its slope. The study in real natural hill slope is a difficult process and the estimation of pore water will not be accurate. To give a clear concept of underground water development the laboratory approach is seemed mostly used even in literature. So for the study, a rainfall simulator is used in which an artificial rainfall of very low to intense rainfall can be made of short as well as long duration over the artificial catchment prepared at a certain slope. As the simulator is portable the slope angle can be made of the desired value. As set up is modeled a hill slope at an angle of 10˚ with the horizontal, the physical process occurred within the catchment represents the rainfall runoff process occurred in real hill terrain. The concept of underground water table, its hydrographs and the moisture profiles may represent or reflects the processes occurs in the real terrain.

Objectives of the Study
The study aims the following objectives  To carry out the low and intense artificial rainfall drop over the artificial catchment to observe its rainfall and runoff process in highly permeable sand slope.  To observe the rise in pore water pressure along the length of hill slope so that the hydraulic grade line can be marked.  The development of the rise in moisture level which represents the water table development process in highly permeable soil and the time for saturation at which the surface runoff starts. The moisture profile immediate after a rainfall event can be observed at different time interval. And to validate, the groundwater development process in sandy soil by kinematic storage model (KSM).

Infiltration in the Ideal & Natural Soil
Infiltration is one of the most important factors in hydrology which influences to the groundwater and its estimation is complex [8]. The ideal soil is considered to be one which is homogeneous throughout the profile, and in which all of the pores are interconnected by capillaries. In addition, it is assumed that the applied rainfall falls uniformly over the soil surface. Because the movement of water into the soil is really uniform, the infiltration process can be considered to be one dimensional. For this ideal case, perhaps the most important factors, which affect the infiltration capacity, are soil type and moisture content. The soil type determines the size and number of the capillaries through which the water must flow. The moisture content determines the capillary potential in the soil and the relative conductivity as shown for a typical soil in Figure 1. For low moisture content, the suctions are high and the conductivity low. At high moisture contents, near saturation, the soil suction is low, and the relative conductivity is high. One can now examine a typical infiltration event-rainfall falling on a dry soil. Because the initial moisture content is low the relative conductivity is small. This implies that the moisture level must be higher before water will move further into the soil mass. Because of this, a distinct wetting front will form with the moisture content ahead of the front still low, and the soil behind the front virtually saturated. A considerable capillary suction exists at the wetting front. In case of highly permeable soil having sufficient depth if the rainfall rate is lower than the infiltration rate or field capacity of soil the infiltrated volume will be greater than at the saturation condition in shallow depth. And for this soil slope, the rise in water table due to a long duration rainfall can be estimated by the well-known equation kinematic storage Model in hill slope [1]. The rainfall intensity its duration is the most important for the moisture storage in hill slope [9]. The infiltration curve goes exponentially in dry soil always [10] as shown in Figure 1. Also, Figure 2 is the hypothetical moisture profile during rainfall event in soil. Figure 3 is the corresponding infiltration curve of Figure 2 during the rainfall event.

Kinematic Storage Model
The moisture storage depends on the soil type, its slope and hill slope geometry [11]. It also depends on the rainfall duration and the time interval. It varies for long time interval record and immediate after the rainfall occurs at shallow and deep soil [12]. There are various analytical models for soil moisture study in both saturated and unsaturated conditions. kinematic wave equation and Hill slope storage Boussinecsq based on darcy law are mostly used model [13]. The hill slope model as shown in Figure 4 represents the hill slope as a rectangular storage element of length L 2.0 m and depth D 0.15 m and width 1 m, making an angle of α (10 degree) with the horizontal. The continuity equation for the saturated zone may be expressed as: where u θ is the average volumetric water content in the unsaturated zone, r is the precipitation input rate, i is the vertical input rate to the saturated zone subscript 1 and 2 represents refer to the beginning and end of the time step, respectively. Sloan and Moore (1984) assumed that the vertical input rate is a function of the volumetric water content, i.e.

( )
where k is hydraulic conductivity.
The continuity equation for the saturated storage element may be expressed as: Or in finite difference form where S is the drainable water stored in the saturated zone per unit with and q is the discharge from the profile per unit width.
where h is the depth of the water table at the outlet and θ drain is the drainable porosity of the soil. With the assumption that the hydraulic gradient is equal to the bed slope, the discharge per unit width at the outlet is given by Darcy's law as: where v is the velocity of flow at the out let and is given by where K s is the saturated hydraulic conductivity and α is the hill slope angle. The water table outlet depth at the end of any time interval ∆t may be determined explicitly by:

Method of Study and Material
The rainfall simulator (Rig) product of Arm field company UK ( Figure 5) has been used for the study having the capacity to generate rainfall event of small to intensive intensity of constant vertical drop size distributed in time and space over the catchment prepared. The setup is useful for the experimental study on sub-surface and surface flow due to a rainfall event. Also, the simulator has a good capacity to study the groundwater table development and effects on hydraulic gradient due to infiltration of water from the soil to the groundwater.
The setup is so made that the change in water pressure can be observed on the     The material used for the simulator should be of 0.2 mm to 2 mm size for this study. The fine and mixed sand can be used for the test material forming the catchment area. The chart in Figure 8 shows the particle size for simulator S12 [15].

Catchment Preparations
As per the material required was collected by sieving the sand in the mixed form of 0.2 to 2 mm and was collected in the tank to construct the catchment. It was filled as per the specification of the instruction manual and marked line within the tank of total area 0.202 m 2 and volume of 0.24 m 3 . The slope was made at 10 degree with the horizontal plane. After the reconnaissance, detailed leveling was done by trial and error method so that accuracy was maintained. As all facilities for controlling, collecting, and runoff reading facilities are already associated became easier to observe, collect and measure runoff. Time observa-Open Journal of Civil Engineering tion was done by stopwatch provided. After the completion of the profile preparation, the test was carried out and data recording were done simultaneously.
As our main aim was to observe the instantaneous runoff discharge, time of concentration and peak discharge due to the different rainfall intensities from 3 l/min to 13 l/min was recorded. The manometer readings were observed at the one-minute intervals with the help of stopwatch which showed the rise of sub-surface water table in soil. Figure 9 below is the constructed soil slope profile within the rig.

Results and Discussion (Model Performance)
The catchment prepared represented the hill slope as a rectangular storage element of Length 2.02 m and Depth of 0. 15 m (Figure 4 Micro topography) with impermeable bed making an angle α (10 degree) with the horizontal plane. The soil is comprised of unsaturated with initial moisture content of 0.23 and saturated moisture content of 0.45. The residual moisture content for sand lays 5% -10%. So 10% was taken for the prediction at the uniform rainfall intensity of 0.25 cm/minute. Permeability of the soil is assumed as 0.3 cm/min During the rainfall, kinematic wave approximation was used for the prediction of the possible velocity. The depth estimation and observation were done at the time of one-minute interval. From the observation after 15-minute rainfall at the rate of 5 ltr/min, the manometer reading was stable at 127 mm. Which indicated the subsurface flow was up to 15 minutes and afterwards overland flow started. Up to that time period, the discharge was observed 0.1125 ltr/se while from the kinematic storage model it was predicted 0.15 ltr/se. Immediate after that rainfall was stopped and the recession value was observed at the one minute time intervals and made comparison with the predicted from kinematic storage model theory. As Figure 10 shows the rainfall runoff simulation and rise in pore water pressure in manometers attached with the simulator As the initial moisture content of the sand was 23%, the change in the moisture content during the rainfall was primary analysis as it contributed to the infiltration rate of the water into the bottom layer of the soil (Figure 11). So the moisture profile from the kinematic storage model was obtained as in Table 1 and Figure 12 for light rainfall period and after rainfall stopped. From the result, it seemed clear raise in the moisture content on the sand maximum of 39% which was near about the saturation moisture content tested and was an acceptable value. After rain, the moisture profile was reverse of rainfall period.
This validated the general principle of the moisture content phenomena in sandy soil due to rainfall. Figure 10. Photos of the rainfall events made, outlet discharges and soil water pressure heads.

Runoff Volume
From the comparative analysis, it seemed the total drainage volume of the water from the catchment of mixed sand was 72 liter whereas the kinematic storage model showed its value of 78.26 liters. But the rainfall rate was of 5l pm, rainfall period of 15 minutes. So the total rainfall volume was 75 liters. However, the rained water was more than 75 liter because even after switching off the sump motor it had drawn water already so there was rain after the stop of the motor run and rainfall stop took time a bit more. Figure 13 was the result obtained.
Hence, from the analysis, the variation of the prediction and observation is not so unrealistic and the experimentation conducted seemed under the KSM theory and obeyed it in terms of volume.

Runoff Hydrograph
After saturation of the soil, the overland flow started at a rainfall rate of 5 lpm and runoff at lps was recorded at the time interval of one minute with the help of stop watch. At the time of concentration (Tc) of 18 minutes, the peak discharge 0.13 lps were recorded. Then after the rainfall was stopped and recession was recorded. Figure 14 below shows the hydrograph obtained on the sand plot soil slope. At the peak discharge, the runoff coefficient was obtained within the range   of 0.4 to 0.52 by rational method for the catchment which is within the standard value as 0.05 to 0.9. Figure 15 shows the variation of runoff coefficients for the slope profile at different rainfall rate according to the rational formula for small catchment producing surface runoff after saturation.

Moisture Profiles
The comparison of moisture profile for sand of sizes 0.2 mm, and 2 mm (mix and fine), over the depth was found as shown in Figure 16. It was observed under the intense rainfall rate of 15 lpm for 5 minute at which the overland flow occurred. Since the infiltration capacity of the fine and mixed sand was high due to the larger size macrospores, significant moisture transaction took place within the pores resulting in the complete adjustment of the moisture. The percolating rate in the larger size of the pore was high so the high moisture content was measured at the mid of the mixed sand. The moisture contents were measured just after the rain in soil slopes. Owing to the low depth, the moisture profiles obtained are not exactly alike as found in the hill slope scale runoff process as explained by Shakya (1995) and S, Scheres (1996). However; the result obtained shows that the timing and distribution of moisture over the depth are a good agreement with that produced by Shakya (1995) and Scherer (1996) for immediate after rainfall occurs in hill slope. Open Journal of Civil Engineering Figure 15. Variation of runoff coefficient with the rainfall rate in sand soil slope. Figure 16. Moisture profile with depth due to intense rainfall in sandy soil immediate after rainfall event.

Conclusions
 Validation by the Hill Slope runoff Process with Kinematic Storage Model to the observation for the sand slope has been ascertained that the experiments obeyed the theory of hill slope model and applicability of the KSM theory for the hill slope runoff analysis with permeable soil. It is observed a realistic result between the observed and predicted values which also shows a good adjustment between KSM theory and physical hill slope runoff process. The peak values for the estimated was found 15 cm whereas observed value was 12.7 cm. The total infiltrated volume was 78.26 liters by the model estimation whereas the rainfall gave was of 75 liter which is considerable values as there is a small variation between estimated and observed.  Owing to the low depth, the moisture profiles obtained are not exactly alike as found in the hill slope scale runoff process as explained by Shakya (1995) and S. Scheres (1996). However; the result obtained shows that the timing and distribution of moisture over the depth are good agreement with that produced by Shakya (1995) and Scherer 1996.