# Estimation of Dry Season Irrigation Water Abstraction in Lunsemfwa, Mulungushi, Mwomboshi, and Mkushi Subbasins from 2013 to 2017 in Zambia

Authors: Lamong Duke Tshenyego; Kamuti Mulonda; Isaac Nyambe Simate
Journal: Advances in Agriculture
 (2019)

---

## Abstract

Water abstraction depends on many variables that include the purpose for the abstraction, the location, the policies in place, and the type of water resources available for abstraction. The overall objective of this study was to estimate irrigation water abstraction from Mkushi, Mulungushi, Mwomboshi, and Lunsemfwa subbasins in Zambia. Reference evapotranspiration was determined using FAO ETo calculator and the results ranged from 6.84 mm/day to 7.02 mm/day. For this study the soils were set as described in the soil map of Zambia and put into the soil characteristic calculator to estimate their physical properties. The results estimate that a total maximum abstraction of 119,680,200 m3 was in 2013, and a minimum estimate of 74,951,400 m3 was in 2014. Wheat abstraction volumes (which were used to represent crops with higher water demand) were compared between catchments and significant differences exist when comparing Lunsemfwa catchment to Mkushi, Mulungushi, and Mwomboshi; thus there were no chances of similarity at an alpha level of 0.05. This means that Lunsemfwa catchment abstracted most irrigation water from 2013 to 2017 than the other three catchments as a result of having the largest proportion of irrigated area in the subbasin.


---

## Body

## 1. Introduction

Irrigation impact in Zambia depends on factors such as the level of demand and importantly on the prioritisation given to irrigation over hydropower production [1]. This is mostly the case in dry areas [2]. Agricultural irrigation is one of the global abstractors of water for approximately 70 % of water in the world. It is projected that irrigated land in developing countries will increase by 27 % in the next 20 years [3]. Water for irrigation is the only option to improve food production in low rain regions but it is affected by the changing weather conditions which may increase evapotranspiration and thus raise crop water requirement [4]. This also has posed some uncertainty about quantity and reliability of the supply [3], as irrigation can result in the reduction of upstream and downstream flows.Irrigation requirement is closely related to crop water requirement. Crop water requirement together with effective rainfall also includes additional water for leaching of salts and compensation for other water losses [5].The quantity of water to compensate for evapotranspiration losses from plant grown in a field is known as crop water requirement. It is normally assumed that the values of crop water requirement and evapotranspiration (ET) are identical [6]. Crop water requirement means the amount of water needed to be supplied to the crop for growth while crop evapotranspiration refers to the amount of water that is lost through evaporation and transpiration. Due to cost and difficulty in obtaining accurate field measurement ET is normally calculated using weather data. Many researchers have analyzed the fit of the various calculation methods in different locations. As an outcome of an Expert Consultation held in May 1990, the FAO Penman-Monteith method is now suggested as the standard method of calculating reference evapotranspiration (ETo) [7, 8]. Reliable estimation of ET is important in that it helps determine the size of canals, pumps, and dams when designing water abstraction for irrigation [9].In Zambia, the Mulungushi dam with a capacity of 2.60×108 m3 and the Mitta hills dam with a capacity of 6.79×108 m3 of water to generate hydropower [10] have been facing the problem of failure to fill up. There are two power stations, one at the Mulungushi dam and another at the Mitta hills dam. The Mulungushi dam power plant has the capacity of producing 31 MW and the Mitta hills power plant has the capacity of producing 54 MW [10]. The company produces power and sells to Zambia Electricity Supply Corporation Limited (ZESCO) which then distributes across the country [11]. This study will estimate the volume of water used in irrigation from the year 2013 to 2017 in the basins of the upper reaches of Lunsemfwa, Mulungushi, Mwomboshi, and Mkushi subbasins. This will help in the planning and harmonization of future water use.


## 2. Methodology

### 2.1. Simulation of Net Irrigation Water by AquaCrop Software

The main model (AquaCrop) used to simulate water abstraction in the basins and its data requirements was used as a guide to collect data. AquaCrop is a menu-driven program with a well-developed user interface. Before running a simulation; the following were specified in the main menu: the sowing date, the simulation period, and the appropriate environmental, initial, and off-season conditions [2]. The ETo calculator was used to calculate ETo. Variables which were used to calculate ETo are minimum and maximum air temperature (Tmin and Tmax, in °C); air humidity as mean relative humidity (RH mean in %); mean wind speed in meter per seconds measured at two (2) meters above ground; sunshine and net radiation (Rn) in MJ/m2.day. ETo was calculated in mm/day from January to December for the period 2013 to 2017. The main equation used by the ETo calculator to calculate ETo is given in (1). Figure 8 shows weather data used to run the ETo calculator from the three weather stations.(1)ETo=0.408ΔRn-G+γ900/T+273u2es-eaΔ+γ1+0.34u2whereE T o denotes reference evapotranspiration [mm day−1], R n denotes net radiation at the crop surface [MJ m−2 day−1], G denotes soil heat flux density [MJ m−2 day−1], T denotes mean daily air temperature at 2 m height [°C], u2 denotes wind speed at 2 m height [m s−1], e s denotes saturation vapour pressure [kPa], e a denotes actual vapour pressure [kPa], e s - e a denotes saturation vapour pressure deficit [kPa], Δ denotes slope vapour pressure curve [kPa °C−1], γ denotes constant [kPa °C−1].Soil Water Characteristics Hydraulic Properties Calculator was used to estimate soil water characteristics by texture as provided on the soil map of Zambia [9]. ETo, soil, and crop parameters were then fed into AquaCrop model by creating corresponding data file (shown in Table 1) for each subbasin to be used in crop growth simulation.Table 1
Information used to run the AquaCrop model.
IMPUT Climate Temp.: Taken from weather stations each month from 2013 to 2017 ETo: From ETo calculator each month from 2013 to 2017 Rainfall: Taken from weather stations each month from 2013 to 2017 CO2: Was left as default of model atmospheric CO2 concentration from 1902 -2099  Crop; Soybean Growing Cycle; day 1 after sowing 2 May to Maturity: 25 July Wheat Growing Cycle; day 1 after sowing 2 May to Maturity: 29 August  Management Determination of net irrigation water requirement  Soil  thickness (m) PWP FC SAT TAW Ksat (mm/hr) Soil Type % % % mm/m Mwomboshi 1.5 29.6 41.3 46.7 121 0.4 Sandy clay Lunsemfwa 1 6.6 11.3 40.3 47 59.9 Loamy Sand Mkushi 1 6.6 11.3 40.3 47 59.9 Loamy Sand Mulungushi 1 6.6 11.3 40.3 47 59.9 Loamy SandThis gave an estimation of water used by the crop from germination to harvest and was used to calculate gross irrigation water requirements under sprinkler irrigation systems. In this study the climate data was collected from the following weather stations: Lusaka Airport used for Mwomboshi catchment; Mulungushi in Kabwe used for Mulungushi and Lunsemfwa catchment; and Serenje used for Mkushi catchment. All the weather stations have digital data recording system and the information is available online under Southern African Science Service Centre for Climate Change and Adaptive Land Management (SASSCAL) program.From the meteorological data, ETo was then calculated using ETo calculator version 3.2 of September 2012.


### 2.2. AquaCrop Crop Data

Crop selection was based on crops grown commercially in the region with wheat as the main irrigated crop in the dry season. Soybean was selected as an alternative crop since it is mainly used in rotation as part of soil fertility management. Additionally, the two crops were representative of major field crops grown in the region and would provide a close estimate of dry season irrigation water abstractions. Table2 gives a summary of crop data as used in AquaCrop.Table 2
Crop parameters used in AquaCrop.
Wheat Soybean Length(days) Date Length (days) Date From day 1 after sowing  2-May  2-May To emergence 15 17-May 7 9-May To maximum Canopy 50 21-May 40 11-Jun To maximum root depth 65 6-Jun 40 11-Jun To start of canopy senescence 80 21-Jul 70 11-Jul To maturity 120 29-Aug 85 25-Jul To flowering 50 21-Jun 50 21-Jul Length building up HI 8 End 29 Jun 32 end 23 July Duration of flowering 15 End 6 Jul 10 end 1 July


### 2.3. Soil Data Use in AquaCrop

In AquaCrop the soil profile can be composed of up to five different horizons of variable depths, each with their own physical characteristics. The considered hydraulic characteristics are the hydraulic conductivity at saturation (Ksat) and soil water content at saturation (θsat), field capacity (θFC), and permanent wilting point (θPWP) [12]. The hydraulic property calculator was used and requires soil texture information in percent by mass of clay and sand (Dr. Keith E. Saxton, 2007). From this information, the type of soil and values for wilting point, field capacity, saturation, saturated hydraulic conductivity, and matric bulk density are estimated [13]. The soil profile is then specified in AquaCrop. For this study the soils were set as described in the soil map of Zambia (Figure 6) which was used at catchment scale to obtain the soil texture information including organic matter content. For this study, crop management was not considered. Monthly rainfall was put in the model for simulation of net irrigation, but mostly the period during simulation was having little or no rainfall (dry season).


### 2.4. Catchments Delineation

The catchment area for the upper Lunsemfwa basin was defined using shape file fromhttps://earthexplorer.usgs.gov/ under United States Geological Survey (USGS). For this study, the 3 sec DEM (30 m x 30 m) data of USGS HydroSHEDS were collected and used for the derivation of catchment and drainage flow directions.The map that covers the area under the study was obtained from Short Topography Mission which is downloadable as a Geotiff file [14]. The downloaded footprint Geotiff files were collected from the Quantum Geography Information System (QGIS) software 2.18.19 which is free to download online [15]. The tiles that form the area under the study were corrected using the Zambia authority ID of European Petroleum Search Group EPSG (32535) latitude and of World Geographical System (WGS) 72BE/UTM zone 35S coordinate referencing system. This was achieved by using spatial reference online to find the right EPSG code, which covers the study area. The delineated area is represented in Figure 1.Figure 1
Catchment orientation and river networks.


### 2.5. Irrigation Area Estimation

For this study, the USGS earth explorer website was used to obtain satellite imagery. The procedures for downloading the free images included registering for USGS accounts and searching the desired area of study and then downloading the Geotiff image from Landsat 8 OLI/TIRS archives [14]. Path 171-171 and 172-172 row 70-70 and 69-69, 10 percent clouds cover day for months of May, July, and September 2013 to 2017 was set.The landsat images were set to the date of early May to August of each year from 2013 to 2017 to cover the period when irrigation is carried out. Landsat imagery (L8OLI/TIRS) was then clipped to the shape files of the delineated catchments in order to classify irrigated farms. The then clipped catchments were classified as seen from Landsat imagery (L8OLI/TIRS) distinguished by different colours and spatial geographical features. Areas under irrigation were calculated with the aid of field calculator in the attribute table of every catchment file in QGIS.For each year, four Landsat imageries (L8OLI/TIRS) for each month of May, July, and August were used, then composited to one image, and clipped to the catchment shape files. This process was repeated for the years 2013 to 2017 in each catchment. Irrigated areas (Figure7) were results of the average of all irrigation classified area for the months of May, July, and August of each year under the study period. Table 3 shows the average irrigated area classified from the QGIS software and the percentage irrigated with respect to the catchment area over five (5) periods.Table 3
Proportion of the various sub-basins under irrigation (wheat/soybean) from 2013 to 2017.
Catchments Study Area (ha) Av. irrigated (ha) % irrigated Mkushi 348000.0 2076.7 0.6 Lunsemfwa 779400.5 13840.0 1.78 Mulungushi 444500.0 2720.0 0.61 Mwomboshi 304300.0 2303.3 0.76The areas were used to calculate the volume of water used after simulation of the net irrigation requirements for each crop in the AquaCrop.


## 2.1. Simulation of Net Irrigation Water by AquaCrop Software

The main model (AquaCrop) used to simulate water abstraction in the basins and its data requirements was used as a guide to collect data. AquaCrop is a menu-driven program with a well-developed user interface. Before running a simulation; the following were specified in the main menu: the sowing date, the simulation period, and the appropriate environmental, initial, and off-season conditions [2]. The ETo calculator was used to calculate ETo. Variables which were used to calculate ETo are minimum and maximum air temperature (Tmin and Tmax, in °C); air humidity as mean relative humidity (RH mean in %); mean wind speed in meter per seconds measured at two (2) meters above ground; sunshine and net radiation (Rn) in MJ/m2.day. ETo was calculated in mm/day from January to December for the period 2013 to 2017. The main equation used by the ETo calculator to calculate ETo is given in (1). Figure 8 shows weather data used to run the ETo calculator from the three weather stations.(1)ETo=0.408ΔRn-G+γ900/T+273u2es-eaΔ+γ1+0.34u2whereE T o denotes reference evapotranspiration [mm day−1], R n denotes net radiation at the crop surface [MJ m−2 day−1], G denotes soil heat flux density [MJ m−2 day−1], T denotes mean daily air temperature at 2 m height [°C], u2 denotes wind speed at 2 m height [m s−1], e s denotes saturation vapour pressure [kPa], e a denotes actual vapour pressure [kPa], e s - e a denotes saturation vapour pressure deficit [kPa], Δ denotes slope vapour pressure curve [kPa °C−1], γ denotes constant [kPa °C−1].Soil Water Characteristics Hydraulic Properties Calculator was used to estimate soil water characteristics by texture as provided on the soil map of Zambia [9]. ETo, soil, and crop parameters were then fed into AquaCrop model by creating corresponding data file (shown in Table 1) for each subbasin to be used in crop growth simulation.Table 1
Information used to run the AquaCrop model.
IMPUT Climate Temp.: Taken from weather stations each month from 2013 to 2017 ETo: From ETo calculator each month from 2013 to 2017 Rainfall: Taken from weather stations each month from 2013 to 2017 CO2: Was left as default of model atmospheric CO2 concentration from 1902 -2099  Crop; Soybean Growing Cycle; day 1 after sowing 2 May to Maturity: 25 July Wheat Growing Cycle; day 1 after sowing 2 May to Maturity: 29 August  Management Determination of net irrigation water requirement  Soil  thickness (m) PWP FC SAT TAW Ksat (mm/hr) Soil Type % % % mm/m Mwomboshi 1.5 29.6 41.3 46.7 121 0.4 Sandy clay Lunsemfwa 1 6.6 11.3 40.3 47 59.9 Loamy Sand Mkushi 1 6.6 11.3 40.3 47 59.9 Loamy Sand Mulungushi 1 6.6 11.3 40.3 47 59.9 Loamy SandThis gave an estimation of water used by the crop from germination to harvest and was used to calculate gross irrigation water requirements under sprinkler irrigation systems. In this study the climate data was collected from the following weather stations: Lusaka Airport used for Mwomboshi catchment; Mulungushi in Kabwe used for Mulungushi and Lunsemfwa catchment; and Serenje used for Mkushi catchment. All the weather stations have digital data recording system and the information is available online under Southern African Science Service Centre for Climate Change and Adaptive Land Management (SASSCAL) program.From the meteorological data, ETo was then calculated using ETo calculator version 3.2 of September 2012.


## 2.2. AquaCrop Crop Data

Crop selection was based on crops grown commercially in the region with wheat as the main irrigated crop in the dry season. Soybean was selected as an alternative crop since it is mainly used in rotation as part of soil fertility management. Additionally, the two crops were representative of major field crops grown in the region and would provide a close estimate of dry season irrigation water abstractions. Table2 gives a summary of crop data as used in AquaCrop.Table 2
Crop parameters used in AquaCrop.
Wheat Soybean Length(days) Date Length (days) Date From day 1 after sowing  2-May  2-May To emergence 15 17-May 7 9-May To maximum Canopy 50 21-May 40 11-Jun To maximum root depth 65 6-Jun 40 11-Jun To start of canopy senescence 80 21-Jul 70 11-Jul To maturity 120 29-Aug 85 25-Jul To flowering 50 21-Jun 50 21-Jul Length building up HI 8 End 29 Jun 32 end 23 July Duration of flowering 15 End 6 Jul 10 end 1 July


## 2.3. Soil Data Use in AquaCrop

In AquaCrop the soil profile can be composed of up to five different horizons of variable depths, each with their own physical characteristics. The considered hydraulic characteristics are the hydraulic conductivity at saturation (Ksat) and soil water content at saturation (θsat), field capacity (θFC), and permanent wilting point (θPWP) [12]. The hydraulic property calculator was used and requires soil texture information in percent by mass of clay and sand (Dr. Keith E. Saxton, 2007). From this information, the type of soil and values for wilting point, field capacity, saturation, saturated hydraulic conductivity, and matric bulk density are estimated [13]. The soil profile is then specified in AquaCrop. For this study the soils were set as described in the soil map of Zambia (Figure 6) which was used at catchment scale to obtain the soil texture information including organic matter content. For this study, crop management was not considered. Monthly rainfall was put in the model for simulation of net irrigation, but mostly the period during simulation was having little or no rainfall (dry season).


## 2.4. Catchments Delineation

The catchment area for the upper Lunsemfwa basin was defined using shape file fromhttps://earthexplorer.usgs.gov/ under United States Geological Survey (USGS). For this study, the 3 sec DEM (30 m x 30 m) data of USGS HydroSHEDS were collected and used for the derivation of catchment and drainage flow directions.The map that covers the area under the study was obtained from Short Topography Mission which is downloadable as a Geotiff file [14]. The downloaded footprint Geotiff files were collected from the Quantum Geography Information System (QGIS) software 2.18.19 which is free to download online [15]. The tiles that form the area under the study were corrected using the Zambia authority ID of European Petroleum Search Group EPSG (32535) latitude and of World Geographical System (WGS) 72BE/UTM zone 35S coordinate referencing system. This was achieved by using spatial reference online to find the right EPSG code, which covers the study area. The delineated area is represented in Figure 1.Figure 1
Catchment orientation and river networks.


## 2.5. Irrigation Area Estimation

For this study, the USGS earth explorer website was used to obtain satellite imagery. The procedures for downloading the free images included registering for USGS accounts and searching the desired area of study and then downloading the Geotiff image from Landsat 8 OLI/TIRS archives [14]. Path 171-171 and 172-172 row 70-70 and 69-69, 10 percent clouds cover day for months of May, July, and September 2013 to 2017 was set.The landsat images were set to the date of early May to August of each year from 2013 to 2017 to cover the period when irrigation is carried out. Landsat imagery (L8OLI/TIRS) was then clipped to the shape files of the delineated catchments in order to classify irrigated farms. The then clipped catchments were classified as seen from Landsat imagery (L8OLI/TIRS) distinguished by different colours and spatial geographical features. Areas under irrigation were calculated with the aid of field calculator in the attribute table of every catchment file in QGIS.For each year, four Landsat imageries (L8OLI/TIRS) for each month of May, July, and August were used, then composited to one image, and clipped to the catchment shape files. This process was repeated for the years 2013 to 2017 in each catchment. Irrigated areas (Figure7) were results of the average of all irrigation classified area for the months of May, July, and August of each year under the study period. Table 3 shows the average irrigated area classified from the QGIS software and the percentage irrigated with respect to the catchment area over five (5) periods.Table 3
Proportion of the various sub-basins under irrigation (wheat/soybean) from 2013 to 2017.
Catchments Study Area (ha) Av. irrigated (ha) % irrigated Mkushi 348000.0 2076.7 0.6 Lunsemfwa 779400.5 13840.0 1.78 Mulungushi 444500.0 2720.0 0.61 Mwomboshi 304300.0 2303.3 0.76The areas were used to calculate the volume of water used after simulation of the net irrigation requirements for each crop in the AquaCrop.


## 3. Results

### 3.1. Climatic Conditions for Mkushi, Lunsemfwa, Mulungushi, and Mwomboshi Catchments

Reference evapotranspiration ETo was determined using FAO ETo calculator [16] for each catchment. Monthly average ETo for five years was plotted as shown in Figure 2. The results show that ETo is highest at an average of 7.02 mm/day in Serenje and Mkushi while the least was in Mulungushi at an average of 6.84 mm/day.Figure 2
Five-year average monthly ETo Mkushi, Mwomboshi, Lunsemfwa, and Mulungushi catchment.


### 3.2. Soil Types Used for Irrigation Water Estimation

Catchment soil was analyzed using the soil water characteristic software [17] and results are shown in Table 1. Table 4 shows the parameter of loamy sands soil found in the Lunsemfwa, Mulungushi, and Mkushi catchment area.Table 4
Soil physical properties.
Soil Parameter Catchements Mkushi Lunsemfwa Mulungushi Mwomboshi Wilting Point ( % ) 6.6 6.6 6.6 29.6 Field Capacity ( % ) 11.3 11.3 11.3 41.3 Saturation ( % ) 40.3 40.3 40.3 46.7 Organic Matter ( % ) 0.5 0.5 0.5 1.2 Saturated Hydraulic Conductivity (mm/hr) 59.9 59.9 59.9 0.4 Soil Type Loamy sand Loamy sand Loamy sand Sandy Clay


### 3.3. AquaCrop Simulation and Net Irrigation Estimation

The results from the catchments and soil characteristics were used to set the AquaCrop model simulation of net irrigating requirements [1] for wheat and soybean. Wheat gave higher net irrigation water equipment than soybean. This helps in estimating the minimum and maximum water that could be abstracted due to irrigation. Net irrigation water requirement is represented in Figure 3. ETo of the crop varies as simulated by the ETo calculator giving alternating values of net irrigation from 2013 to 2017. Mulungushi and Lunsemfwa show similar values as a result of same weather and soil type [18]. From Figure 3 cumulative ETo for wheat and soybean are different because of differences in length of growing season. Wheat recorded higher values for cumulative ETo as it was simulated for longer period (120 days) while soybean was simulated for shorter periods (85 days). ETo influenced the net irrigation water, due to alternating solar radiation, humidity, and wind speeds recorded from weather station resulting in alternating ETo. It was also due to variability in air temperature affecting both ETo and the length of crop growing season (Allen 2008).Figure 3
Yearly analysis of net irrigation per hectare and evapotranspiration from soybean and wheat in Mkushi, Lunsemfwa, Mulungushi, and Mwomboshi catchments.


### 3.4. Volumes of Water Abstraction

The volume of water was estimated from irrigated areas (Figure7) and gross irrigation water requirements for wheat and soybean. The assumption was that the crops were irrigated by sprinkler with an efficiency of 0.75. Figure 4 shows that all irrigation water abstracted from the study area would increase depending on the crop irrigated. It was estimated that the average abstraction in all catchments was 107,727, 813 m3. The maximum volume was 119,680,200 m3 overall in 2013, while the minimum was estimated to be in 2014 at 74,951,400 m3. In simple terms the maximum irrigation water abstraction decreased from 2013 (119,680,200 m3) to 2015 (94,686,566.66 m3), increased to 118,193,500 m3 in 2016, and reduced in 2017 to 100,715,566 m3. The minimum irrigated water abstraction drop from 98,272,900 m3 in 2013 to 77,944,500 m3 in 2014, then increased to 93,979,933.33 m3 in 2015, and lastly dropped in 2016 and 2017 to 74,951,400 m3. The summary is included in Table 5 and Figure 5.Table 5
Parameter of study area.
Parameter Catchments Mkushi Lunsemfwa Mulungushi Mwomboshi Catchment area (km2) 3480.0 7794.48 4445.0 3043.0 Irrigated area (km2) 20.8 138.4 27.2 23.0 Average Yearly Rainfall (mm) 1032.68 1165.54 1165.54 882.64 Max Volume abstracted (m3) 66413317 3.40E+08 6.60E+07 6.81E+07 Min Volume abstracted (m3) 46737683 4.67E+07 5.48E+07 4.10E+07 Soil Type Loamy sand Loamy sand Loamy sand Sandy Clay LoamFigure 4
Average wheat and soybean water abstraction from Lunsemfwa.Figure 5
(a) Total annual abstraction volumes and (b) average annual for each crop.
(a) (b)Figure 6
Soil map of Zambia which was used to obtain the soil characteristics of the four catchments under studyhttps://www.google.co.zm/search?q=soil+map+of+zambia&rlz=1C1CHBF.Figure 7
Classified landsat images from study area showing classified irrigated area.Figure 8
Weather data from Serenje, Kabwe, and Lusaka airport Zambia.


## 3.1. Climatic Conditions for Mkushi, Lunsemfwa, Mulungushi, and Mwomboshi Catchments

Reference evapotranspiration ETo was determined using FAO ETo calculator [16] for each catchment. Monthly average ETo for five years was plotted as shown in Figure 2. The results show that ETo is highest at an average of 7.02 mm/day in Serenje and Mkushi while the least was in Mulungushi at an average of 6.84 mm/day.Figure 2
Five-year average monthly ETo Mkushi, Mwomboshi, Lunsemfwa, and Mulungushi catchment.


## 3.2. Soil Types Used for Irrigation Water Estimation

Catchment soil was analyzed using the soil water characteristic software [17] and results are shown in Table 1. Table 4 shows the parameter of loamy sands soil found in the Lunsemfwa, Mulungushi, and Mkushi catchment area.Table 4
Soil physical properties.
Soil Parameter Catchements Mkushi Lunsemfwa Mulungushi Mwomboshi Wilting Point ( % ) 6.6 6.6 6.6 29.6 Field Capacity ( % ) 11.3 11.3 11.3 41.3 Saturation ( % ) 40.3 40.3 40.3 46.7 Organic Matter ( % ) 0.5 0.5 0.5 1.2 Saturated Hydraulic Conductivity (mm/hr) 59.9 59.9 59.9 0.4 Soil Type Loamy sand Loamy sand Loamy sand Sandy Clay


## 3.3. AquaCrop Simulation and Net Irrigation Estimation

The results from the catchments and soil characteristics were used to set the AquaCrop model simulation of net irrigating requirements [1] for wheat and soybean. Wheat gave higher net irrigation water equipment than soybean. This helps in estimating the minimum and maximum water that could be abstracted due to irrigation. Net irrigation water requirement is represented in Figure 3. ETo of the crop varies as simulated by the ETo calculator giving alternating values of net irrigation from 2013 to 2017. Mulungushi and Lunsemfwa show similar values as a result of same weather and soil type [18]. From Figure 3 cumulative ETo for wheat and soybean are different because of differences in length of growing season. Wheat recorded higher values for cumulative ETo as it was simulated for longer period (120 days) while soybean was simulated for shorter periods (85 days). ETo influenced the net irrigation water, due to alternating solar radiation, humidity, and wind speeds recorded from weather station resulting in alternating ETo. It was also due to variability in air temperature affecting both ETo and the length of crop growing season (Allen 2008).Figure 3
Yearly analysis of net irrigation per hectare and evapotranspiration from soybean and wheat in Mkushi, Lunsemfwa, Mulungushi, and Mwomboshi catchments.


## 3.4. Volumes of Water Abstraction

The volume of water was estimated from irrigated areas (Figure7) and gross irrigation water requirements for wheat and soybean. The assumption was that the crops were irrigated by sprinkler with an efficiency of 0.75. Figure 4 shows that all irrigation water abstracted from the study area would increase depending on the crop irrigated. It was estimated that the average abstraction in all catchments was 107,727, 813 m3. The maximum volume was 119,680,200 m3 overall in 2013, while the minimum was estimated to be in 2014 at 74,951,400 m3. In simple terms the maximum irrigation water abstraction decreased from 2013 (119,680,200 m3) to 2015 (94,686,566.66 m3), increased to 118,193,500 m3 in 2016, and reduced in 2017 to 100,715,566 m3. The minimum irrigated water abstraction drop from 98,272,900 m3 in 2013 to 77,944,500 m3 in 2014, then increased to 93,979,933.33 m3 in 2015, and lastly dropped in 2016 and 2017 to 74,951,400 m3. The summary is included in Table 5 and Figure 5.Table 5
Parameter of study area.
Parameter Catchments Mkushi Lunsemfwa Mulungushi Mwomboshi Catchment area (km2) 3480.0 7794.48 4445.0 3043.0 Irrigated area (km2) 20.8 138.4 27.2 23.0 Average Yearly Rainfall (mm) 1032.68 1165.54 1165.54 882.64 Max Volume abstracted (m3) 66413317 3.40E+08 6.60E+07 6.81E+07 Min Volume abstracted (m3) 46737683 4.67E+07 5.48E+07 4.10E+07 Soil Type Loamy sand Loamy sand Loamy sand Sandy Clay LoamFigure 4
Average wheat and soybean water abstraction from Lunsemfwa.Figure 5
(a) Total annual abstraction volumes and (b) average annual for each crop.
(a) (b)Figure 6
Soil map of Zambia which was used to obtain the soil characteristics of the four catchments under studyhttps://www.google.co.zm/search?q=soil+map+of+zambia&rlz=1C1CHBF.Figure 7
Classified landsat images from study area showing classified irrigated area.Figure 8
Weather data from Serenje, Kabwe, and Lusaka airport Zambia.


## 4. Discussion

### 4.1. Catchment Evapotranspiration

The findings for the period of study suggests that Serenje weather station, as used for weather to calculate ETo for Mkushi catchment [11], has highest ETo while the lowest was in Mwomboshi catchment. This is in agreement with geographical locations of the catchment as Serenje (located in the northern end of Mkushi catchment) is further north and closer to the equator while Kabwe (used for Mwomboshi catchment) is in the southern end and slightly far from the equator. Figure 3 reviews a close relationship between seasonal net irrigation and ETo for both wheat and soybean. Thus net irrigation is influenced by cumulative ETo during the growing season [19] as evidenced by relatively high water demands per hectare in Mkushi catchment. For the same crops and soil type simulation in Mkushi, Lunsemfwa, and Mulungushi, the differences in weather influenced seasonal ETo and consequently net irrigation requirements (Karamet al., 2007). Despite using catchment averaged soil and weather characteristics, the resulting irrigation water requirements are in an acceptable range of 410 mm–510 mm for soybeans and 450 mm–710 mm for wheat. The FAO guideline gives a range from 450 mm to 700 mm and 450 mm to 650 mm for soybeans and wheat, respectively.


### 4.2. Irrigation Areas and Volumes in the Basins

According to the findings, Lunsemfwa catchment is the largest irrigated area as affirmed by [20] with Mkushi catchment showing low irrigated areas as demonstrated in Figure 7. Although other catchments have availability of water somehow the Lunsemfwa area has more access to roads and electricity as well as good soil for growing crops leading to high irrigated area. Table 6 shows the variation in mean irrigated area for all catchments. Although the classification process was good in identifying most irrigated areas under pivot irrigation systems, poor resolution of free online satellite images has limited accuracy.Table 6
Multiple Comparisons of wheat Net irrigation water [mm/ha] 2013 to 2017 LSD.
(I) Catchment (J) Catchment Mean Difference (I-J) Std. Error Sig. 95% Confidence Interval Lower Bound Upper Bound Lunsemfwa Mkushi -146.66 44.045 0.01 -242.625 -50.694 Mwomboshi -96.42 44.045 0.05 -192.385 -0.454 Mkushi Lunsemfwa 146.66 44.044 0.01 50.694 242.625 Mwomboshi 50.24 44.044 0.28 -45.725 146.205 Mwomboshi Lunsemfwa 96.42 44.044 0.05 0.454 192.385 Mkushi -50.24 44.044 0.28 -146.205 45.725 “Based on observed means: The error term is Mean Square (Error) = 4849.887.” ∗. The mean difference is significant at the 0.05 level.Further study or research should be carried out using remote sensed images with higher resolutions.Volume of irrigation abstraction is greatly influenced by area under irrigation. It can be seen from Table5 that Lunsemfwa basin has the highest volume while Mkushi has the lowest and this same trend is observed for irrigation areas, respectively. Although amount of areas irrigated has impact on the amount of volume abstracted it was not the case with Mwomboshi which has second highest abstraction and third highest irrigated area (Table 5). This was due to the different soil type (Sandy clay) from the other three catchments as well as the thickness of the soil horizon (Table 1).


## 4.1. Catchment Evapotranspiration

The findings for the period of study suggests that Serenje weather station, as used for weather to calculate ETo for Mkushi catchment [11], has highest ETo while the lowest was in Mwomboshi catchment. This is in agreement with geographical locations of the catchment as Serenje (located in the northern end of Mkushi catchment) is further north and closer to the equator while Kabwe (used for Mwomboshi catchment) is in the southern end and slightly far from the equator. Figure 3 reviews a close relationship between seasonal net irrigation and ETo for both wheat and soybean. Thus net irrigation is influenced by cumulative ETo during the growing season [19] as evidenced by relatively high water demands per hectare in Mkushi catchment. For the same crops and soil type simulation in Mkushi, Lunsemfwa, and Mulungushi, the differences in weather influenced seasonal ETo and consequently net irrigation requirements (Karamet al., 2007). Despite using catchment averaged soil and weather characteristics, the resulting irrigation water requirements are in an acceptable range of 410 mm–510 mm for soybeans and 450 mm–710 mm for wheat. The FAO guideline gives a range from 450 mm to 700 mm and 450 mm to 650 mm for soybeans and wheat, respectively.


## 4.2. Irrigation Areas and Volumes in the Basins

According to the findings, Lunsemfwa catchment is the largest irrigated area as affirmed by [20] with Mkushi catchment showing low irrigated areas as demonstrated in Figure 7. Although other catchments have availability of water somehow the Lunsemfwa area has more access to roads and electricity as well as good soil for growing crops leading to high irrigated area. Table 6 shows the variation in mean irrigated area for all catchments. Although the classification process was good in identifying most irrigated areas under pivot irrigation systems, poor resolution of free online satellite images has limited accuracy.Table 6
Multiple Comparisons of wheat Net irrigation water [mm/ha] 2013 to 2017 LSD.
(I) Catchment (J) Catchment Mean Difference (I-J) Std. Error Sig. 95% Confidence Interval Lower Bound Upper Bound Lunsemfwa Mkushi -146.66 44.045 0.01 -242.625 -50.694 Mwomboshi -96.42 44.045 0.05 -192.385 -0.454 Mkushi Lunsemfwa 146.66 44.044 0.01 50.694 242.625 Mwomboshi 50.24 44.044 0.28 -45.725 146.205 Mwomboshi Lunsemfwa 96.42 44.044 0.05 0.454 192.385 Mkushi -50.24 44.044 0.28 -146.205 45.725 “Based on observed means: The error term is Mean Square (Error) = 4849.887.” ∗. The mean difference is significant at the 0.05 level.Further study or research should be carried out using remote sensed images with higher resolutions.Volume of irrigation abstraction is greatly influenced by area under irrigation. It can be seen from Table5 that Lunsemfwa basin has the highest volume while Mkushi has the lowest and this same trend is observed for irrigation areas, respectively. Although amount of areas irrigated has impact on the amount of volume abstracted it was not the case with Mwomboshi which has second highest abstraction and third highest irrigated area (Table 5). This was due to the different soil type (Sandy clay) from the other three catchments as well as the thickness of the soil horizon (Table 1).


## 5. Conclusions and Recommendations

The volumes of irrigation water abstracted from the catchments were estimated to be highest in the Lunsemfwa catchment at average of 6.78×107 m3 in the five-year period followed by an averages of 1.36×107 m3 in Mwomboshi catchment and Mkushi at 1.32 ×107 m3 and lowest average in Mulungushi at 1.31×107 m3. The large abstraction is due to bigger irrigation area in the Lunsemfwa catchment. The following are the recommendations:(1)
Doing bathymetric studies on the dams upstream to ascertain the amount of water farmers are capable of storing. This is to make abstraction studies more accurate and to improve information availability for future planning.(2)
To compliment remote sensing applications, field surveys and experiments need to be carried out to validate the software results.



---

*Source: 8979837-2019-07-14.xml*