Quantifying rainfall-runoff relationships on the Melkassa Hypo Calcic Regosol ecotope in Ethiopia

Droughts, resulting in low crop yields, are common in the semi-arid areas of Ethiopia and adversely influence the wellbeing of many people. The introduction of any strategy that could increase yields would therefore be advantageous. The objective of this study was to attempt to assess the influence of in-field rainwater harvesting (IRWH), compared to conventional tillage, on increasing the amount of water available to a crop like maize on a semi-arid ecotope at Melkassa situated in the eastern part of the Rift Valley. To achieve the objective of the study rainfall-runoff measurements were made during 2003 and 2004 on 2 m x 2 m plots provided with a runoff measuring system and replicated 3 times for each treatment. There were 2 treatments: conventional tillage (CT) on which hand cultivation was practised in a way that simulated the normal local CT; and a flat surface simulating the no-till, undisturbed surface of the IRWH technique (NT). Rainfall-runoff measurements were made over 2 rainy seasons during which there were 25 storms with > 9 mm of rain. From the 25 storms, only the 2nd season storms (8 storms) had runoff measurements. These storms were used for calibration and validation of the Morin and Cluff (1980) runoff model (MC Model). Appropriate values for final infiltration rate (If), surface storage (s) and for the crusting parameter (γ) were found to be: 6 mm∙h-1; 1.0 mm for NT and 6.0 mm for CT; 0.6 mm-1, respectively. The measured runoff (R) for the 2004 rainy season expressed as a fraction of the rainfall during the measuring period (P), i.e. R/P, gave values of 0.59 and 0.40 for the NT and CT treatments, respectively. There was a statistical difference between the runoff on the 2 treatments. Selected results from 7 years of field experiments with IRWH at Glen in South Africa were used together with measured maize yields and climate data over 16 seasons on the nearby Melkassa Experiment Station to estimate the yield benefits of IRWH compared to CT on the ecotope studied. The results ranged between 35 and 1 437 kg with a mean of 711 kg∙ha-1 over the 16 years. At Melkassa this was an estimated yield increase ranging from 13% to 49%. The mean increase was 33%.


CT
= conventional tillage D i = accumulated rain depth during time segment i (mm) DoY = day of year E s = evaporation from the soil surface (mm) ET = evapotranspiration (mm) ETo = reference evapotranspiration (mm) ET inc = increment in evapotranspiration (mm) (see Table 5) F Δti = total infiltration (mm) during time segment Δti with rainfall intensity I f = final infiltration rate (mm•h -1 ) I i = initial infiltration rate of the soil (mm•h -1 ) I t = instantaneous infiltration rate (mm•h -1 ) IRWH = in-field rainwater harvesting MC Model = Morin and Cluff (1980)  decreased by crusting; the extent to which the soil surface can store water before runoff starts which is described by a parameter termed surface detention (SD); and a crusting parameter (γ) describing its rate and extent of development.Their studies resulted in the formulation of a runoff model that satisfactorily predicted runoff from crusted soils in Arizona (Morin and Cluff, 1980), and in Israel (Morin et al., 1983).The model has been successfully used by Zere et al. (2005) for predicting the runoff measured by Du Plessis and Mostert (1965) over 18 years on a Tukulu form soil (Soil Classification Working Group, 1991) at Glen, South Africa.
The basis for the Morin and Cluff (1980) runoff model is provided by the following infiltration equation for crusted soils developed by Morin and Benyamini (1977): (1) Morin and Cluff (1980) showed that by integrating Eq. ( 1) with regard to time, and introducing changes in Pi over time segments of a storm (Δti), the following expression was valid: where: F Δti = total infiltration during time segment Δti with rainfall intensity Pi (mm) The other parameters are as defined for Eq. ( 1).
If the soil surface were such that it did not store any water before runoff occurred, then for any time segment during which Pi > I f , the runoff for each time segment Δti of a storm, i.e.Ri, could be calculated as: This is, however, not the case in practice as a soil surface always has some degree of surface roughness which will cause rain-water to accumulate, to an extent dependent on the degree and configuration of the roughness before runoff commences.Morin and Cluff (1980) deal with this factor by combining the Di of Eq. ( 2) and 'a surface detention' parameter (SD) into a parameter, SDm, termed 'maximum storage and detention'.By introducing this term into Eq.( 3) they showed that it was possible to compute the runoff of any storm, segment by segment, using the following equation: Substitution of the right-hand side of Eq. ( 2) into the F Δti term of Eq. ( 4) provides the complete Morin and Cluff (1980) runoff equation, i.e.Eq. ( 5):

Introduction
More than 80% of Ethiopia's population is involved in agriculture, the backbone of the country's economy.Crop production is mostly under rain-fed conditions, most of which is marginalised by water stress (Ministry of Agriculture (MoA), 2000).This, and the frequent droughts, is a serious threat to those engaged in agriculture.The optimum utilisation of rain-water is therefore of utmost importance, requiring diligent adherence to the principle of 'more crop per drop', as appropriately stated recently by the former UN Secretary General, Kofi Annan.In scientific terms this means improving RWP, recently defined by Botha (2007) as the total long-term grain yield divided by total long-term rainfall.
One way of improving RWP is through the use of water harvesting.Many types of water conservation techniques that show significant crop yield increases have been tested worldwide (Berry and Mallet, 1988;Mwakalila and Hatibu, 1993;Kronen, 1994;Gicheru et al., 1998 andOjasvi et al., 1999).A technique that has given good results in a semi-arid area of South Africa is IRWH as described by Hensley et al. (2000).This technique is also known as mini-catchment runoff farming (Oweis et al., 1999).The technique is illustrated in Fig. 1.It combines the advantages of water harvesting from the notill, flat, crusted runoff strip, and decreased evaporation from the deeply infiltrating runoff water which accumulates in the mulched basin area.The technique led to maize yield increases of between 25% and 50% compared to conventional tillage practices, and resulted in significant increases in RWP.It was shown that the technique is suited to semi-arid areas with crusting soils that have a high water storage capacity (Botha et al., 2003).
Rainfall in semi-arid areas with fine-textured soils is mainly lost through evaporation from the soil surface (Es) and runoff (R).Under these conditions Es can be 60% to 70% of the annual rainfall (Bennie and Hensley, 2001), and R can vary between 8% and 49% of the annual rainfall depending on the prevailing conditions (Haylett, 1960;Du Plessis and Mostert, 1965;Bennie et al., 1994;Hensley et al., (2000) and Botha et al., 2003).Studies by Morin and Benyamini (1977) and Morin and Cluff (1980) showed that the most important factors influencing runoff in semi-arid areas were: rainfall intensity (Pi); the final infiltration rate of the soil (I f ), which is greatly

Hypothesis
• The in-field water harvesting technique described in Fig. 1 will result in increased crop yields compared with conventional tillage on certain semi-arid ecotopes of Ethiopia.• The MC Model will satisfactorily predict runoff on the chosen ecotopes.• It will be possible to make reasonable estimates of yield increases on the selected ecotopes using IRWH by predicting the extent of runoff collected in the basins, and therefore prevented from leaving the field and becoming unavailable to the crop.

Objectives
• To quantify rainfall-runoff relationships on the semi-arid Melkassa ecotope in Ethiopia over 2 rainy seasons.• To calibrate the MC Model for the Melkassa ecotope.
• To estimate for the Melkassa ecotope the maize yield benefits using the IRWH technique described in Fig. 1, compared to conventional tillage.Data from the first objective will be used to do this.

Procedure Study site
The study was carried out at Melkassa in one of the semi-arid regions of Ethiopia for 2 main rainy seasons during 2003 and 2004.Melkassa is located in the central part of the rift valley at longitude 39.31 o E and latitude 8.43 o N. The altitude is 1 550 m a. m. s. l. and the chosen site represents a gently sloping plain with a slope ranging from 0% to 5% comprising a foot slope of the rift valley.The ecotope is described by the geographic site name followed by the name of the soil.The soil is classified as a Hypo Calcic Regosol (WRB classification).The ecotope name is therefore Melkassa Calcic Fluvic Regosol.This soil covers about 10% of Ethiopia and about 16% of the rift valley (FAO, 1984;FAO, 1998b;Itanna, 2005).

Experimental design
The experiment was carried out at Melkassa Agricultural Research Center (MARK) research field with a slope of 1%.There were 2 treatments and 3 replications in a random ised complete block design.The plot size was 2 m x 2 m.The treatments were: • Conventional tillage (CT), i.e. the normal/traditional farmer's tillage practice • No tillage on a flat surface (NT), i.e. simulating the runoff strip of IRWH.
Both treatments were uncropped and weeds were controlled by hand weeding.The lower side of each plot was equipped with a runoff collecting device.Each plot was surrounded by a galvanised iron sheet protruding 20 cm to 30 cm above the surface of the soil, and inserted about 20 cm deep into the soil.This 'wall' served to isolate each plot hydraulically.Runoff was collected in a gutter at the lower side of the plot.The gutter channelled the runoff water into a 200 ℓ barrel buried at the side of each plot.Runoff data were collected for each rainfall event.The MC Model describes a rainstorm as a group of rain segments for which the breaks in the rain are less than 24 h.Huff (1967) defines a storm as a rain period separated from a preceding and succeeding rainfall event by 6 h or more.The latter definition was used.Runoff was simply measured by recording the height of the runoff inside the barrel.
Rainfall amount and intensity was measured by an automatic tipping bucket rain gauge (Hobo Event (C) Onset Computer Corp, Model No. 7, Version No. 4) installed at the experimental site to store detailed data for every storm.Each bucket tip measures 0.2 mm in a time interval determined by the intensity of the rainfall.The rain gauge is capable of measuring 0.2 mm in 0.01 s.The rain gauge was equipped with a data-logger with memory capacity of 32 768 bytes.The data were downloaded to a laptop computer, and then prepared to 1 min intensity.The record included the starting date and time, as well as the terminating date and time of each storm.The data collected were analysed to characterise each rainstorm at the Melkassa ecotope during the measuring period.

Climate
The Melkassa Hypo Calcic Regosol ecotope is located about 15 km south-east of Nazret City.The main rainy season is during the months June to September, during which 68% of the annual rainfall occurs (Table 1).The measured Class A pan evaporation data (Eo) and the potential evapotranspiration (ETo), calculated using the Penman-Monteith equation, correlate well with R 2 of 0.92.The highest evaporative demand occurs during the months of March, April and May.During these months, the mean maximum temperature (Txm) is around 30 o C while the mean relative humidity (RHm) drops to 51%.During the main crop growing season of June to September conditions are more favourable with Txm and RHm approximately 27 o C and 64% respectively.(FAO, 1998a) According to the recent agroecological zones classification of Ethiopia (MoA, 2000), the Melkassa Hypo Calcic Regosol ecotope falls in the zone termed hot to warm semi-arid lowlands (SA1).This belt exhibits 2 growing seasons of 50 d and 100 d in length for the 1 st and 2 nd seasons, respectively, and has an annual rainfall (P) and potential evapotranspiration (PET) of about 772 mm and 1 994 mm, respectively.The aridity index (AI) (P/PET) of 0.39 (Table 1) identifies this as a semi-arid area.

Soil
A profile pit was dug to a depth of 3 000 mm.The soil profile was described and classified as follows: Hypo Calcic Regosol according to the World Resource Base System (FAO, 1998b); Etosha Vetkuil (2111) according to the South African System (Soil Classification Working Group, 1991); Regosol according to the FAO system (FAO, 1984).A soil map (FAO, 1998a) of the Rift Valley in this vicinity shows the dominance of the Regosols (Fig. 2).Important characteristics of the Melkassa soil are a favourable clay loam texture of the fine earth throughout the profile, with high silt content.The topsoil is strongly crusting.The water holding capacity of the potential root zone for maize is considered to be high.
Determinations of the following physical properties were made: drainage curve; soil water retention curves; bulk density; initial (I i ) and final (I f ) infiltration rates.Detailed results are presented in Welderufael (2006).

Calibration and validation of the MC Model
The measured rainfall and runoff data during the main season of 2004 were used to calibrate and validate the MC Model.Half of the data were used for calibration and the other half for validation.The data were used together with the determined values of I i and I f to run the model.The remaining parameters in the model, i.e. maximum surface detention (SD m ), and γ were fixed using a sensitivity analysis to obtain 'best fit' values.Model calibration was carried out by changing the values of γ between 0.1 and 0.9 and SD m between 0 and 10 mm, while keeping the measured and first approximation I i and I f values fixed.Once the optimum values for γ and SD m were obtained, the sensitivity analysis was conducted manually and expertly (Madsen et al., 2002) to improve the I f value until the performance evaluation functions had reached their optimum level, and the observed and simulated runoff values matched reasonably well.Once the models were calibrated and the parameters fixed, validation was carried out on the remaining data using the procedure of Willmott (1981).

Measurements
Rainfall amounts (P) and intensities (Pi) were measured during the main rainy seasons of 2003 and 2004.Runoff (R) measurements were taken only for 2004.Melkassa storms generally exhibited intense rainfall during the 1 st and 2 nd quartiles of the events.Huff (1967) in his study at Illinois in the USA also found a similar pattern.In the 2003 rainy season the total amount of the rainfall for events ≥ 9 mm was 297 mm.It was uniformly distributed throughout the season.There were six storms in July, seven in August and four in September, with 103 mm, 114 mm and 80 mm rainfall amounts, respectively (Table 2 A).After calibration the MC Model predicted R/P as 0.3 and 0.16 on NT and CT plots respectively.Most of R (57%) came from the 3 big storms on DoY 199, 236 and 249.For 2004, rainfall events ≥ 9 mm totalled 210 mm.Rainfall for the entire season amounted to 251 mm producing R/P values of 0.6 and 0.4 on the NT and CT plots respectively (  2003, the rainfall distribution in 2004 was non-uniform with 2 large storms in July followed by 1 large storm in August, and 1large storm each in September and October.This pattern would have caused a shortage of water during the flowering and maturity stage of the cropping period.During 2004, where measured runoff data were available, there was a significant difference at the 0.05 probability confidence level between the runoff on the 2 cultivation practices, with an overall mean per storm of 15.5 mm and 10.5 mm on the NT and CT plots respectively.The significant difference is attributed to the larger SD m values of the CT plots, presumably due to the 2 cultivation practices carried out on them during the season, and also due to the relatively few heavy storms capable of producing a similar crusted surface to that on the NT plots.The CT plot was cultivated at the beginning of the study and after the storm on DoY 212.This left a rough surface with considerable depressions that persisted for a longer period throughout the rainy season.
The MC Model was calibrated and validated using the rainfall-runoff measurements for 2004.The validated model was then used to predict R for each storm of the 2003 season.
A similar procedure for model calibration and validation to that used at Dera was followed (Welderufael, 2006).Results are presented in Tables 3 A and B and Table 4. Appropriate values for I f and γ were found to be 6 mm•h -1 and 0.6 mm -1 , respectively.These are the same as those selected for Dera (Welderufael, 2006).
The calibration procedure revealed that the s values (= SD m ) which gave the best results with the MC Model were 1 mm and 6 mm for the NT and CT plots, respectively.The following criteria were used for making the decision: D-index, R 2 , and RMSEu/RSME as close as possible to 1.0.A high value of the latter parameter is of particular importance since it indicates that the error is mainly not of a systematic nature.
These values, and their assessment parameters, are printed in bold (Table 3).Results of the validation test using these values, and the I i , I f and γ values ( Once the model had been calibrated and validated, it was used to simulate the runoff of each storm during both years (Table 2).The good agreement between the measured and predicted values is reflected by the R 2 values of 0.86 and 0.94 for the NT and CT plots, respectively, for all the storms during the 2 seasons (P≥9 mm).

Well-simulated storms
Included are storms that start with intense rainfall (P i > I f ), and those storms that acquire high intensities (P i > I f ) later than in the 1 st quartile, and continued with P i > I f for sufficient time to fulfil the sorptivity and SD m demand of the soil.A study of the rainfall vs. time graphs for a number of storms of this type indicates that about 4 mm of rain is needed to satisfy the requirements of sorptivity and SD m (1 mm) on the NT plots.Therefore, subtracting 4 mm from the cumulative rainfall value at the point where P i becomes less than I f , or at the point where the steepest part of the cumulative rainfall line terminates, will directly give an estimated amount of runoff on the NT plots.
Figures 3 and 4 show storms that begin with high intensity, during 2003 and 2004 respectively.For the storm on DoY 199 of year 2003, the point where P i < I f is indicated by an arrowed line.The P i > I f part of the storm lasted for about 48 min.The arrow gives a value of 18 mm on the y-axis of the cumulative rainfall.Therefore, subtracting 4 mm (the value of sorptivity + SD m ) from 18 mm will give 14 mm.This is a similar result to the amount of runoff simulated by the Morin and Cluff (1980) runoff model for NT plot which equals 13.9 mm (Table 2).Similarly, for CT plots if we subtract 3 mm plus the value of SD m for CT plots (3 + 6 = 9), we will obtain 9 mm of runoff.Again the result is very close to the one estimated by the model as 10.2 mm (   The storm on DoY 252 of 2004 (Fig. 4) had Pi > I f throughout its 22 min duration giving 22 mm of cumulative rainfall.Subtracting 4 mm from 22 mm gives an expected runoff of 18 mm on NT plots.The measured runoff was 16.4 mm, while the model simulated 18.2 mm.Using the same calculation as for storm on DoY 199 the expected runoff on CT plots is 13 mm.The measured and simulated values were 11.4 mm and 13.8 mm respectively.Storms with P i < I f were also well simulated, in all cases giving zero runoff.
The 2 nd group of storms that were well simulated were characterised by P i > I f for a certain period during the middle of the storm's duration (2 nd , 3 rd or 4 th quartiles).It can be assumed that the sorptivity and SD m demand for these storms was satisfied by the rain that fell before the intense part started, or else by bursts of intense rains (P i > I f ) that occurred before or after the major intense period (m.i.p.).Huff (1967) defined 'burst' as a cessation in rainfall or an abrupt, persistent change in rainfall rate.But here 'burst of intense rain' was taken as part of the storm that showed P i > I f for a short time interval compared to the m.i.p. of P i > I f .At Melkassa this kind of storm was rare. Figure 5 shows one of these storms (DoY 277) during 2004.In this storm the starting point and end-point of the m.i.p. are indicated by the arrowed lines giving 14 mm and 6 mm of rainfall on the cumulative rainfall y-axis.Thus, by subtraction, 14 mm minus 6 mm, gives 8 mm of expected runoff, the same as the measured value.The model simulated 8.2 mm.Similarly, for storm on DoY 249 of 2003 (Fig. 6) the arrowed lines indicate 26 mm and 7 mm of the cumulative rainfall as boundary values of the m.i.p.This gives 19 mm of expected runoff on NT plots, while the model simulated 18.2 mm.
It is clear that a long dry period between storms will increase the sorptivity of the soil.In addition, high SD m values were encountered when storms occurred immediately after cultivation on CT plots..Both these factors will influence the accuracy of simulations.Unlike the Dera Calcic Fluvic Regosol ecotope (Welderufael, 2006), the CT plots on the Melkassa Hypo Calcic Regosol ecotope retained an almost similar SD m value (6 mm) throughout the 2004 rainy season.This may be due to the smaller number of intense rain events after the 2 nd cultivation practice carried out on DoY 215.

Examples of storms not well simulated
Storm on DoY 223 of 2004 gave an exceptionally high measured runoff value on the NT plot of 25.9 mm whereas the simulated value was only 12.2 mm (Fig. 7).The high measured R was probably due to the fact that it occurred 48 h after 2 continuous storms on DoY 220 and 221.Although these storms produced little R (1.8 mm and 2.3 mm from NT) they probably contributed enough water to leave the soil surface wet after 48 h.As a result the demand for sorptivity was minimised (an indication of how the model could be improved).The other relevant factor was the occurrence of continuous small bursts of P i > I f that lasted for about 133 min, between 49 min and 182 min (Fig. 7).They covered approximately 3 quarters of the storm's duration.These bursts may not have been considered by the model as significant enough time segments to produce runoff.Similarly, storms on DoY 232 and 251 of 2004 and 2003, respectively, were under-simulated by the model.

Estimating yield increases using IRWH
Empirical procedures were followed to estimate the benefit of IRWH to maize production on the Melkassa Hypo Calcic Regosol ecotope.Maize yields and climate data from the Melkassa Agricultural Research Center (MARC) for 16 growing seasons (1988 to 2003) were used.The average maize yield for this period using conventional tillage was found to be 2 115 kg•ha -1 (Table 5).The climate data included rainfall (P), temperature (T), relative humidity (RH), sunshine hours (SH) and wind speed (WS).
The CROPWAT programme developed by FAO was used for making Es + T estimates.(Es + T) is termed ET in the programme.The detailed description of this programme and the calculations used are given in Welderufael, 2006.The programme first makes use of the climatic data needed (i.e.P, T, RH, SH, and WS) to calculate ETo for each day of each growing season using the Penman-Monteith equation.In the 2 nd step the programme combines a crop coefficient (Kc) with ETo to estimate the potential ET of maize (Kc*ETo) for each growth stage, i.e. the amount of water it would require for ET to attain maximum yield.For each of the 16 growing seasons the ratio of the final grain yield to the sum of the ET values for the season yielded the water productivity (WP ET ) for that particular season (Table 5).The mean WP ET value over all the seasons was found to be 6.5kg•ha -1 •mm -1 (Table 5).
To proceed further it is necessary to have an estimate of the fraction of the extra water produced by runoff on the NT plots that will become used for increasing yield, i.e. in this case used specifically for ET.The results obtained by Hensley et al. (2000) and Botha (2007) for field experiments comparing the IRWH and CT production techniques with maize on the Glen/Bonheim ecotope, over 7 growing seasons, was employed as follows to estimate this fraction.The following information was extracted for each growing season: • Infield runoff (R inf ) from the IRWH treatment with a bare runoff area • The difference in water used for ET on IRWH compared to CT (ET IRWH -ET CT = ΔET) • The ratio of ΔET/ R inf The average value of ΔET/ R inf over the 7 seasons was 0.62.This indicates that on average, on the Glen Bonheim ecotope with maize, ET IRWH can be expected to be increased to the extent of 0.62 x R if above the ET of maize with conventional tillage, i.e.ΔET ≈ 0.62 * R if .A comparison of the runoff characteristics of the Melkassa ecotope and the Glen/Bonheim ecotope shows that they have similar characteristics.Their I f values are also the same (6 mm•h -1 ).It is therefore a reasonable 1 st approximation to employ the calculated ET:R if relationship for the Glen/ Bonheim ecotope on the Melkassa Hypo Calcic Regosol ecotope.
The following is the description of the procedure used to estimate the expected maize yield increment with IRWH.Results are presented in  * 1 R if = estimate of infield runoff based on measurements made on this ecotope during 2003and 2004(Welderufael, 2006)).* 2 ET inc = estimate of the increase in ET due to R if * 3 WP ET = water productivity calculated based on evapotranspiration (ET) NT plot runoff to rainfall events > 9 mm was developed.The equation is: R = 0.714P -6.8959 (R 2 = 0.87) where: R is the estimated runoff in mm P is the amount of rain (> 9 mm) of the rainfall event (Welderufael, 2006) Applying this equation to each rainfall event during each growing season from 1988 to 2003 (Table 5) provides an estimate of what the runoff (R if ) would have been from the runoff strip (Fig. 1) had IRWH been employed.Multiplication of this value by 0.62 gives an estimate of the expected ET increment (ET inc ).The multiplication of ET inc *WP ET provides a logical estimate of the increased yield with IRWH.
Results are presented in Table 5.Values vary between 13% and 49%.The mean increase is shown to be 33%, which represents an estimated average annual yield increase of 711 kg•ha -1 .

Conclusions
The 3 objectives of the study were achieved.Firstly, the Morin and Cluff runoff model was successfully calibrated and validated.Appropriate values for the 3 parameters needed by the model for use on the Melkassa Hypo Calcic Regosol were determined, i.e.I f = 6 mm•h -1 ; s for NT and CT were 1 mm and 6 mm respectively; and γ = 0.6 mm -1 .Secondly, rainfallrunoff relationships on the Melkassa ecotope during 2004 were quantified giving values of 0.59 and 0.40 for the NT and CT treatments, respectively.The significant difference between the runoff on the 2 treatments during 2004 was caused mainly by the 2 cultivation operations on the CT plots that caused large SDm and had a major influence on runoff.Thirdly, maize yield benefits using the IRWH technique instead of conventional tillage on this ecotope were estimated to be between 35 and 1 437 kg•ha -1 .
The study shows how crop yields in semi-arid regions of Sub-Saharan Africa could be increased significantly by employing in-field rainwater harvesting rather than conventional tillage.Since it is expected that the technique will only be successful on ecotopes with specific properties, prior detailed characterisation of these relevant properties is recommended.

Figure 2
Figure 2Soil map of the Rift Valley showing the experimental site(FAO, 1998a) Figure 4 Storm on DoY 252 of 2004 on the Melkassa Hypo Calcic Regosol ecotope Figure 5 Storm on DoY 277 of 2004 on the Melkassa Hypo Calcic Regosol ecotope

TABLE 2 Measured and simulated runoff during the 2003 (A) and 2004 (B) rainy seasons at Melkassa (A) A
Table 3) are presented in Table 4. Significant aspects for both NT and CT are the relatively low RMSE values, very high D-index and R 2 values, and very high RMSEu/RMSE values.These are also compatible to the good overall correlation finally obtained between measured and predicted runoff values during 2004 for both the NT and CT treatments; R 2 values were 0.86 and 0.94, respectively.