Cultivation Features Using Meteorological and Satellite Data from 2001 to 2010 in Dakhla Oasis, Egypt ()
1. Introduction
Agriculture in Egypt is strongly correlated to distance from the Nile River [1] . Figure 1 illustrates the distribution
Figure 1. Distribution of NDVI, which indicates the vegetation amount, in Egypt on June 2012 from satellite data acquired by MODIS. Red circles indicate the locations of the oases.
of the Normalized Difference Vegetation Index (NDVI), which indicates vegetation amount, in Egypt in June 2012 using the MODe rate resolution Imaging Spectroradiometer (MODIS) satellite image data set (details will be discussed later). NDVI is high only along the Nile River and in the delta area, and not much vegetation can be seen far from the Nile River because of the hyperarid climate. However, some vegetation is found in the desert area far from the Nile River, although the NDVI there is little compared to that of the delta area. These oasis areas have been developed from time immemorial.
Cultivation in these oasis areas depends on fossil water (artesian ground water) originating from the Nubian aquifer [2] . Exploitable ground water resources are very large. However, because there is essentially no rainfall in this region, water is not rechargeable; therefore, the water wells are destined to dry up with time. The ground water levels of aquifers differ according to the underground geological conditions; thus, new water wells continue to be exploited, and old water wells and the land are abandoned. Cultivation in oasis areas is therefore vulnerable from the viewpoint of sustainable water use.
Traditional agriculture has focused on date palm in Dakhla oasis. Systematic agriculture has been conducted since the 1960s using a government well and rotating cultivation, as instructed by the government. In this study, we focus on cultivation features throughout the government and local well districts using meteorological and satellite data from 2001 to 2010, and the interview in Rashda Village, Dakhla Oasis.
2. Study Site and Analysis Methods
2.1. Site Description and Analysis Method
Dakhla, a major oasis in the Western Desert of Egypt (Figure 1), is located on low-level ground, at 0 to 200 m above sea level. The aridity index used by the United Nations Environment Program [3] indicates that this district lies in a hyperarid region, with annual rainfall of nearly 0 mm. Consequently, the surrounding area is bare of vegetation, except for very isolated clumps of Acacia or Tamarix that survive on the shallow moisture held in blown sand, and moisture from dew, isolated showers, and frontal storms [4] . The climate of Dakhla is controlled mainly by the Harmattan, which is a hot wind resulting from high pressure in the Sahara. Although wind directions are mostly north, a slight westerly skewing of the maximum probably reflects topographic in draft of northern winds into the Dakhla depression [4] .
Meteorological data between 2001 and 2010 sourced from National Climatic Data Center (NCDC) Climate Data Online of the Dakhla Meteorological Station (station number 624,320, 25.48˚N, 29.0˚E, 117 m) was used to analyze the meteorological features and reference evapotranspiration (ET0) defined by [5] . Daily ET0 (m×md−1) as expressed by [5] is
(1)
where Δ is the slope vapor pressure curve (kPa×˚C−1), Rn is the net radiation (MJ×m−2×day−1), G is the soil heat flux (MJ×m−2×day−1), γ is the psychrometric constant (kPa×˚C−1), T is the air temperature (˚C), U is the wind speed (m×s−1), es is the saturation vapor pressure (kPa), and ea is the actual vapor pressure (kPa). The procedure for calculating ET0 is described in detail by [5] . The effect of G is ignored for actual daily calculations. For Rn calculation, the data of sunshine duration is needed. Since it was not observed at the Dakhla Meteorological Station, we referred to the measured relative sunshine duration (ratio of actual sunshine duration to possible sunshine duration) for 1990 to 2005 at the Kharga Meteorological Station (station number 624,350, 25.45˚N, 30.53˚E, 73 m), which is 120 km from Dakhla [6] . The annual mean percentage of relative sunshine duration was high at 87% (minimum at 82% in January, maximum at 90% in June). The observational height of U in Equation (1) is defined as 2 m. Because the actual observational height is 10 m, we converted the observed wind speed using following formula [5] :
(2)
where U2 is the wind speed at 2 m (m×s−1) and Uz is the observed wind speed at z m (m×s−1).
Daily crop evapotranspiration (ET) (m×md−1) was calculated by multiplying ET0 by crop coefficient Kc (dimensionless):
(3)
Table 1 lists the Kc value [5] over each growth stage of the main cultivated crops in Rashda Village. Since these Kc values are determined under typical irrigation management and soil wetting, ET in Equation (3) is assumed as the evapotranspiration for non-stressed conditions. The seasonal variation of Kc was determined from the growing schedule using the interview (Table 1), the actual growth pattern using the NDVI, and the Kc value in Table 1.
LANDSAT-5/TM and LANDSAT-7/ETM+ images were used to detect the amount of vegetation and water in the government and local well districts. We obtained 108 images during the analysis period between 2001 and 2010. NDVI provides a measure of the amount of vegetation on the land surface [7] . In general, higher NDVI values indicate greater amounts of vegetation. NDVI produces values of ?1 to 1. The negative values indicate the presence of water bodies (e.g., oceans, rivers, lakes, and snow). Small positive values indicate bare ground with little or no green vegetation, while values approaching 1 are a sign of lush green vegetation. [8] indicated that 0.05 ≤ NDVI < 0.2 is bare land or bare land surface with less vegetation, based on the [9] . The NDVI can be estimated from the reflectances of bands 3 and 4 using the following formula:
(4)
(5)
where ρ4 and ρ3 are the reflectances for bands 4 and 3, λ is the band number, Lλ is the sensor observed radiance
Table 1. Cultivation schedule, crop coefficient value (Kc) in the initial, mid-season and end-season periods, and evapotranspiration of main cultivated crops in Rashda.
for band λ (W×m−2×ster−1×μm−1), ESUNλ is the mean solar exoatmospheric irradiance for band λ (W×m−2×ster−1×μm−1), θ is the solar incident angle normal to the surface, and dr is the inverse squared relative distance between the Earth and the Sun.
The Normalized Difference Water Index (NDWI) is a satellite-derived index acquired from the NIR and SWIR reflectances, and provides a measure of the amount of water in vegetation canopies [10] . NDWI gives a positive value for green vegetation and a negative one for dry vegetation. It is calculated using the following formula, which is the same as that used to determine NDVI:
(6)
where ρ5 is the reflectance for band 5.
The government well district (area 70 ha, center location 25˚35'35''N, 28˚55'17''E, 97 m) and the local well district (area 22 ha, center location 25˚35'05''N, 28˚56'04''E, 99 m) were classified polygonally using the image in Figure 2, and the averaged NDVI and NDWI values were calculated. Each of the 34 farmers in the government well district owned land in the North, South, and West sub-districts. Averaged NDVI and NDWI values were calculated for each sub-district in the government well district and the local well district. Bands 3, 4, and 5 had a resolution of 30 m.
(a)(b)
Figure 2. The North (green), South (red), and West (blue) sub-districts in the government well district (a) and the local well district (b). Red circles indicate the locations of the wells. These maps were made using Google map.
We also used the MODIS-Terra monthly 1 km NDVI product (MODIS/Aqua Vegetation Indices Monthly L3 Global 1 km SIN Grid V005*, http://modis.gsfc.nasa.gov/) to clarify the spatial distribution of vegetation amount in Egypt (Figure 1).
2.2. Irrigation and Cultivation System in the Government and Local Well Districts in Rashda Village
The water source in Dakhla Oasis and throughout the Western Desert is the Nubian sandstone aquifer, which is the world’s largest fossil water aquifer system. Therefore, the irrigation and cultivation system in Dakhla Oasis differs from that in the Nile delta and valley, which depends completely on the Nile River as its water source. The ground water of the Nubian aquifer was recharged 14,000 to 30,000 years ago, and the last rainy season related to the recharge occurred 8000 years ago, according to isotopic analysis [2] . Currently, rainfall provides little recharge water to the recharge area (northwest part of Sudan and Uwein at Upland located at the boundary of Egypt, Libya, and Sudan), so the irrigation water in Dakhla Oasis is mainly fossil water having ancient roots [11] .
In Dakhla Oasis, the underground water is extracted using a well. The estimated number of wells and springs is 70 in Rashda Village. [1] provided detailed information on these wells. Government wells are utilized by the government (Ministry of Irrigation) and are owned and managed by the Ministry of Irrigation. They generally have a depth of 1200 m and are supposed to last 50 years. There are 12 wells of this type in Rashda. Local wells are artesian wells that are dug in the traditional way, and are owned and managed collectively by local cultivators. Their depth is usually 85 m or less. The well leader maintains a list of cultivators who have access to the water source. The water source is supposed to last 20 years. There are 29 wells of this type in Rashda. Figure 2(a) indicates the locations of five government wells, and Figure 2(b) indicates that of one local well.
Irrigation water flows from wells through main and branch canals into parcels of cultivated land through the main and branch distributors. Irrigation is then managed according to the basin irrigation system in both government and local well districts. The interview indicates that in the government well district, 22 m3 of water is irrigated to each feddan (1 feddan = 0.42 ha) once every 12 days from October to April, and 40 m3 of water is irrigated to each feddan once every 12 days from May to September. In the local well district, 30 m3 of water is irrigated to each feddan once every 10 days through the year.
The prevailing crops in the government well district are wheat or clover hay (berseem) in winter, and rice in summer under the following rotation:
Dates are generally produced on land associated with the local well. The land associated with governmental wells does not have palm trees. Water management, crop types grown, and their rotation are decided by the Department of Agriculture in the government well district. However, this rule is not strictly applied in areas irrigated by local wells.
3. Results and Discussion
3.1. Meteorological Features in Dakhla Oasis
Figure 3 indicates the average seasonal variation of daily average temperature, wind speed, vapor pressure, and daily ET0 from 2001 to 2010. The annual daily average temperature was 24.4˚C (highest 34.9˚C in July and lowest 13.0˚C in January). The annual average daily maximum temperature was 32.8˚C (highest 42.5˚C in July and lowest 20.9˚C in January). The annual average daily minimum temperature was 16.0˚C (highest 26.4˚C in July and lowest 4.6˚C in January). The annual average daily vapor pressure was 1.1 kPa (highest 1.5 kPa in August and lowest 0.7 kPa in February). Dakhla is comparatively humid because irrigation farming is widespread. The annual daily average wind speed was 1.9 m∙s−1 (strongest 2.9 m∙s−1 in April and weakest 1.0 m∙s−1 in November).
3.2. ET0 and Crop Evapotranspiration Using Meteorological Data
The average seasonal variation of daily ET0 from 2001 to 2010 is indicated in Figure 3. The annual average ET0
Figure 3. Seasonal variation of daily averaged temperature, daily maximum temperature, daily minimum temperature, daily averaged wind speed, daily averaged vapor pressure, and daily averaged ET0 between 2001 and 2010.
was 5.0 m∙md−1 (greatest 7.9 m∙md−1 in July and least 2.0 m∙md−1 in December). The annual ET0 was 1816 m∙my−1. If actual sunshine duration was assumed to be the same as possible sunshine duration, annual ET0 was 1884 m∙my−1. The difference was 68 m∙my−1 (0.19 m∙md−1), which is within observational error [12] . Thus, possible sunshine duration can be used to calculate ET0 in Dakhla Oasis. Since ET0 in winter is less than in summer, cultivation in winter makes sense from the viewpoint of water saving, as discussed in the next section.
Evapotranspiration for each crop during the growing season is indicated in Table 1. It should be noted that this evapotranspiration was calculated assuming sufficient wet conditions; therefore, actual evapotranspiration will be less than these values. Under the meteorological conditions in Dakhla, water requirements are highest for date palm, and second highest for rice paddy. According to the interview, farmers preferred to grow wheat or clover hay rather than date palm in government well districts because the cultivation of date palm requires an abundance of water. However, date palm is still cultivated in local well districts. Water demand for winter wheat is comparatively low. However, because most nutrients in the soil are consumed in its cultivation, a fallow period and clover hay cultivation are necessary after wheat production to recover land conditions. It is noteworthy that water requirements for winter wheat are almost the same as for clover hay.
3.3. Seasonal Variation of NDVI for Government and Local Well Districts
Figure 4 indicates the seasonal NDVI variation (date palm) for each sub-district in the government and local well districts between 2001 and 2010. In Figure 5, the solid line represents a sixth-order polynomial equation, and four polynomial lines meet. In the government well district, the NDVI value reached maximum (0.6) in mid-February and decreased drastically until the end of May. Although the NDVI value decreased in the order of South, North, and West sub-districts during the winter season (November to May), their differences were small. A significant change in NDVI can be seen among the sub-districts during the summer season (June to October). After harvesting of winter plants (winter wheat or clover hay), some crops (e.g., rice) were planted and cultivated from the beginning of November in the South sub-district, although the scale of cultivation was
Figure 4. Seasonal variation in the NDVI for each sub-district of the government well district and the local well district between 2001 and 2010. The solid line indicates a sixth-order polynomial equation.
small compared to that of the winter season; however, the land remained fallow in the North and West sub-dis- tricts most years (Figure 4). Here, the fallow period was defined as when the NDVI was below 0.2 and NDWI was below 0, as discussed in Section 2.1. The high irrigation water volume required for cultivation in the summer climate means that fallow or small-scale cultivation is more suitable for Rashda Village during the summer from the viewpoint of water saving. Additionally, the government regulates rice paddy cultivation, based on the amount of water available for irrigation in the summer, under current irrigation laws. Such a rotation system in the government well districts could be confirmed even in the interview.
The NDVI value of date palm in the local well district was almost constant. However, a slight increase began in May and reached its maximum (0.5) in December. According to the interview, the date palm harvesting period was from September to November, and the average yield was 120 kg per tree (Figure 5).
3.4. Seasonal Variation of NDWI for Government and Local Well Districts
Figure 6 indicates the seasonal NDWI variation for each sub-district in the government well district and the local well district (date palm) between 2001 and 2010. In Figure 7, the solid line represents a sixth-order polynomial equation, and four polynomial lines. The seasonal variation trend of NDWI was similar to that of NDVI. The NDWI value of date palm was low in summer and began to increase toward the harvesting period. Since the NDWI provides a measure of the water amount in vegetation canopies, a low NDWI value might indicate water stress conditions for date palm in summer. The irrigation amount was constant through the year in the local well district (30 m3 of water irrigated to each feddan once every 10 days). The irrigation amount might be low for date palm in summer, when evapotranspiration requirements are large.
If the accumulated day-to-day NDWI value is assumed to be the total water requirement for each crop, differences in water requirements among crops can be estimated from the satellite data only. Figure 7 indicates that the accumulated NDWI was 109 for date palm in the local well district, and 42 for winter plants (winter wheat or clover hay) in the government well district, based on the averages of the North (42), South (45), and West (38) sub-districts during the winter season (November to May). Accumulated NDWI of date palm was 2.6 times that of winter plants. The given irrigation amounts were 257 mm for date palm and 102 mm for winter plants according to the interview; thus, irrigation amounts for date palm were 2.5 times those for winter plants. Table 1 also indicates that evapotranspiration of date palm was 2.5 times that of winter plants (average of winter
Figure 5. Seasonal variation in the NDVI using a sixth-order polynomial equation (indicated in Figure 4).
Figure 6. Seasonal variation in the NDWI for each sub-district of the government well district and the local well district between 2001 and 2010. The solid line indicates a sixth-order polynomial equation.
wheat and clover hay). It should be noted that this evapotranspiration was calculated assuming sufficient wet conditions; thus, actual evapotranspiration will be less than these values. Actually, irrigation amounts were much smaller than the calculated evapotranspiration in Table 1. These results may indicate the possible use of NDWI to detect differences in water requirements among crops. It will be important to compare the accumulated NDWI with actual evapotranspiration for future development of this concept.
4. Conclusions
Dakhla Oasis is located in the heart of the western desert in Egypt, 190 km west of Kharga Oasis. It contains highly fertile land and is rich in water, and it supports a higher population than Kharga Oasis. The Dakhla
Figure 7. Seasonal variation in the NDWI using a sixth-order polynomial equation (indicated in Figure 6).
depression extends 155 km in the east-west direction. The area of the oasis that is suitable for agriculture is 155 km long and 60 km wide. In 2012, the number of wells increased 30% compared to that in 2008. In recent years, irrigation farming has spread significantly in Dakhla Oasis. To extend the use of the limited groundwater resources, efficient irrigation farming is crucial.
This study focuses on cultivation features throughout the government and local well districts using the meteorological and satellite data from 2001 to 2010, and the interview in Rashda Village, Dakhla Oasis. It is one of a limited number of interdisciplinary studies involving the natural and social sciences, and will contribute to fill the lack of research on the cultivation features in the oasis region.
1) Since solar radiation and temperature are suitable for crop production, cultivation intensity depends on available groundwater in Dakhla. Since ET0 in winter is less than in summer, cultivation in winter makes sense from the viewpoint of water saving.
2) Evapotranspiration is highest in the order of date palm, rice paddy, winter wheat, and clover hay under sufficient wet conditions in Rashda Village.
3) Approaches to detect cultivation features and water requirements were examined using NDVI and NDWI for the government and local well districts in Rashda Village. There was a distinct difference in cultivation systems between the summer and winter seasons in the government well district. Fallow or small-scale cultivation was applied in summer because of the high evapotranspiration demand and irrigation water saving. Water requirements for date palm were 2.6 times those for winter plants (winter wheat or clover hay) using the accumulated day-to-day NDWI value.
Our analysis has the following limitations.
1) Data on solar radiation were lacking in Dakhla Oasis. The western desert of Egypt was dominated by high pressure on the synoptic scale (2000 to 10,000 km) throughout the year. Thus, we assumed that the solar radiation condition in Dakhla Oasis was similar to that in Kharga Oasis 120 km from Dakhla Oasis.
2) Actual evapotranspiration data on each crop were not obtained. The proposed method to detect the water requirements using NDWI should be verified using the actual observed data of evapotranspiration for future development.
However, the observation system will be expensive and difficult to establish for a variety of reasons. Data obtained by satellite and the interview used in this study will also become a useful tool for analyzing cultivation features in arid regions, including the oasis area in Egypt.