Advances in Agriculture
Volume 2019 (2019), Article ID 8979837, 9 pages
https://doi.org/10.1155/2019/8979837

Research Article

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

Lamong Duke Tshenyego, Kamuti Mulonda, and Isaac Nyambe Simate

Department of Agricultural Engineering, University of Zambia, Lusaka, Zambia

Correspondence should be addressed to Lamong Duke Tshenyego; dukelamong@gmail.com

Received 6 February 2019; Revised 15 April 2019; Accepted 17 April 2019; Published 14 July 2019

Academic Editor: Gábor Kocsy

Copyright © 2019 Lamong Duke Tshenyego et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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 m³ was in 2013, and a minimum estimate of 74,951,400 m³ 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.

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×10⁸ m³ and the Mitta hills dam with a capacity of 6.79×10⁸ m³ 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/m².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.

where

  denotes reference evapotranspiration [mm day⁻¹],  denotes net radiation at the crop surface [MJ m⁻² day⁻¹], G denotes soil heat flux density [MJ m⁻² day⁻¹], T denotes mean daily air temperature at 2 m height [°C], u₂ denotes wind speed at 2 m height [m s⁻¹],  denotes saturation vapour pressure [kPa],  denotes actual vapour pressure [kPa],  denotes saturation vapour pressure deficit [kPa],  denotes slope vapour pressure curve [kPa °C⁻¹],  denotes constant [kPa °C⁻¹].

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

CO₂:
Was left as default of model atmospheric CO₂ 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 Sand

This 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. Table 2 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 (), 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 from https://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 (Figure 7) 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.76

The 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 (Figure 7) 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 m³. The maximum volume was 119,680,200 m³ overall in 2013, while the minimum was estimated to be in 2014 at 74,951,400 m³. In simple terms the maximum irrigation water abstraction decreased from 2013 (119,680,200 m³) to 2015 (94,686,566.66 m³), increased to 118,193,500 m³ in 2016, and reduced in 2017 to 100,715,566 m³. The minimum irrigated water abstraction drop from 98,272,900 m³ in 2013 to 77,944,500 m³ in 2014, then increased to 93,979,933.33 m³ in 2015, and lastly dropped in 2016 and 2017 to 74,951,400 m³. The summary is included in Table 5 and Figure 5.

Table 5: Parameter of study area.

Parameter
Catchments

Mkushi
Lunsemfwa
Mulungushi
Mwomboshi

Catchment area (km²)
3480.0
7794.48
4445.0
3043.0

Irrigated area (km²)
20.8
138.4
27.2
23.0

Average Yearly Rainfall (mm)
1032.68
1165.54
1165.54
882.64

Max Volume abstracted (m³)
66413317
3.40E+08
6.60E+07
6.81E+07

Min Volume abstracted (m³)
46737683
4.67E+07
5.48E+07
4.10E+07

Soil Type
Loamy sand
Loamy sand
Loamy sand
Sandy Clay Loam

Figure 4: Average wheat and soybean water abstraction from Lunsemfwa.

(a)

(b)

Figure 5: (a) Total annual abstraction volumes and (b) average annual for each crop.

Figure 6: Soil map of Zambia which was used to obtain the soil characteristics of the four catchments under study https://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 (Karam et 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 Table 5 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×10⁷ m³ in the five-year period followed by an averages of 1.36×10⁷ m³ in Mwomboshi catchment and Mkushi at 1.32 ×10⁷ m³ and lowest average in Mulungushi at 1.31×10⁷ m³. 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.

Appendix

A. Input Weather Data Used to Calcuate ETo in ETo Calculator

See Figures 2, 3, 4, 5, 6, and 7 and Tables 4 and 5.

Weather Data. See Figure 8.

B. SPSS Least Significant Different Comparisons Results

Based on observed means, the error term is mean square (error) = 47340650502465.200.”

The mean difference is significant at the 0.05 level.

See Table 6.

Data Availability

The data used to support the findings of this study are included as the supplementary information file.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

I am highly obliged to my sponsors, Mobility to Enhance Training of Engineering Graduates in Africa (METEGA), first for awarding me the scholarship to pursue a MEng-Agricultural Engineering at the University of Zambia and secondly covering the tuition fees and the other necessary support such as medical insurance and a monthly stipend. Without these, this paper would not have been possible. Honest appreciation goes to my supervisors Mr. K. Mulonda and Dr. Isaac N. Simate for their insightful contribution that helped me stay on track and finish this paper within the time frame allocated to the study.

Supplementary Materials

Excel file used to support the findings of this study. (Supplementary Materials)

References

R. Stricevic, M. Cosic, N. Djurovic, B. Pejic, and L. Maksimovic, “Assessment of the FAO AquaCrop model in the simulation of rainfed and supplementally irrigated maize, sugar beet and sunflower,” Agricultural Water Management, vol. 98, no. 10, pp. 1615–1621, 2011.

C. F. Yamba, H. Walimwipi, H. Kling et al., “Water supply and demand scenarios for the Zambezi river,” Basin Climate Change and Upstream Development Impacts on New Hydropower Projects in the Zambezi Project Report for Climate & Development Knowledge Network Lead Author: Randall Spalding-Fecher, 2014.

J. K. Kiprop, J. K. Lagat, P. Mshenga, and A. M. Macharia, “Determining the economic value of irrigation water in Kerio Valley Basin (Kenya) by residual value method,” Journal of Economics and Sustainable Development, vol. 6, no. 7, pp. 102–108, 2015.

T. Uhlendahl, P. Salian, C. Casarotto, and J. Doetsch, “Good water governance and IWRM in Zambia: Challenges and chances,” Water Policy, vol. 13, no. 6, pp. 845–862, 2011.

F. N. Gichuki, “Hydrological impacts of flood storage and management on irrigation water abstraction in upper Ewaso Ng’iro River Basin, Kenya,” Water Resources Management, pp. 1859–1879, 2008.

J. Osorio, J. Jeong, K. Bieger, and J. Arnold, “Influence of Potential Evapotranspiration on the Water Balance of Sugarcane Fields in Maui, Hawaii,” Journal of Water Resource and Protection, pp. 852–868, 2014.

R. G. Allan, L. S. Pereira, D. Raes, and M. Smith, “Crop evapotranspiration-Guidelines for computing crop water requirements-FAO Irrigation and drainage paper 56,” 1998.

R. G. Allen, L. S. Pereira, D. Raes, M. Smith, and W. Ab, “Crop evapotranspiration - Guidelines for computing crop water requirements,” FAO Irrigation and Drainage Paper 56, 1998.

K. T. Zeleke and L. J. Wade, “Evapotranspiration estimation using soil water balance, weather and crop data,” in Evapotranspiration - Remote Sensing and Modeling, 2009.

B. Lawson, “Mulungushi hydropower project – Environmental and Social Impact Assessment (ESIA),” Mulungushi Hydropower Project Draft Scoping Report, Environmental Resources Management, 2013.

M. Everett, “Muchinga Hydropower Project – Environmental and Social Impact Assessment (ESHIA) Muchinga Hydropower Project Draft Scoping,” Report Version 2, 2013.

D. Raes, P. Steduto, T. C. Hsiao, and E. Fereres, “Calculation procedures AquaCrop,” 2012.

D. Raes, P. Steduto, T. C. Hsiao, and E. Fereres, “Users guide AquaCrop,” 2012.

F. Akram, M. G. Rasul, M. M. K. Khan, and M. S. I. I. Amir, “Automatic Delineation of Drainage Networks and Catchments using DEM data and GIS Capabilities: a case study,” 2012.

T. Athan, P. Ersts, and O. D. Werner Macho, “QGIS User Guide,” 2018. https://docs.qgis.org/2.18/pdf/en/QGIS-2.18-UserGuide-en.pdf.

D. Raes, “The ETo Calculator Table of Contents,” 2012.

K. E. Saxton and W. J. Rawls, “Soil water characteristic estimates by texture and organic matter for hydrologic solutions,” Science Society of America Journal, vol. 70, no. 5, pp. 1569–1578, 2006.

A. Chapman, H. Walimwipi, and H. Kling, “The vulnerability of hydropower production in the Zambezi River Basin to the impacts of climate change and irrigation,” Mitigation and Adaptation Strategies for Global Change, vol. 21, no. 5, pp. 721–742, 2016.

S. Er-Raki, A. Chehbouni, N. Guemouria, B. Duchemin, J. Ezzahar, and R. Hadria, “Combining FAO-56 model and ground-based remote sensing to estimate water consumptions of wheat crops in a semi-arid region,” Agricultural Water Management, vol. 87, no. 1, pp. 41–54, 2007.

I. Malasha, “Food security and poverty alleviation through improved valuation and Governance of river fisheries in africa the governance of small scale fisheries in Zambia,” Submitted to the Research Project on Food Security and Poverty Alleviation through Improve 1, 2016, http://pubs.iclarm.net/wfcms/file/bmz/Zambia%20governance.pdf.