Utilization of Landsat-8 data for the estimation of carrot and maize crop water footprint under the arid climate of Saudi Arabia

The crop Water Footprint (WF) can provide a comprehensive knowledge of the use of water through the demarcation of the amount of the water consumed by different crops. The WF has three components: green (WFg), blue (WFb) and grey (WFgr) water footprints. The WFg refers to the rainwater stored in the root zone soil layer and is mainly utilized for agricultural, horticultural and forestry production. The WFb, however, is the consumptive use of water from surface or groundwater resources and mainly deals with irrigated agriculture, industry, domestic water use, etc. While the WFgr is the amount of fresh water required to assimilate pollutants resulting from the use of fertilizers/agrochemicals. This study was conducted on six agricultural fields in the Eastern region of Saudi Arabia, during the period from December 2015 to December 2016, to investigate the spatiotemporal variation of the WF of silage maize and carrot crops. The WF of each crop was estimated in two ways, namely agro-meteorological (WFAgro) and remote sensing (WFRS) methods. The blue, green and grey components of WFAgro were computed with the use of weather station/Eddy covariance measurements and field recorded crop yield datasets. The WFRS estimated by applying surface energy balance principles on Landsat-8 imageries. However, due to non-availability of Landsat-8 data on the event of rainy days, this study was limited to blue component (WFRS-b). The WFAgro of silage maize was found to range from 3545 m3 t-1 to 4960 m3 t-1; on an average, the WFAgro-g, WFAgro-b, and WFAgro-gr are composed of < 1%, 77%, and 22%, respectively. In the case of carrot, the WFAgro ranged between 297 m3 t-1 and 502 m3 t-1. The WFAgro-g of carrot crop was estimated at <1%, while WFAgro-b and WFAgro-gr was 67% and 32%, respectively. The WFAgro-b is occupied as a major portion in WF of silage maize (77%) and carrot (68%) crops. This is due to the high crop water demand combined with a very erratic rainfall, the irrigation is totally provided using groundwater delivered by center pivot irrigation systems. On the other hand, the WFRS-b estimated using Landsat-8 data was varied from 276 (±73) m3 t-1 (carrot) and 2885 (±441) m3 t-1 (silage maize). The variation (RMSE) between WFRS-b and WFAgro-b was about 17% and 14% for silage maize and carrot crops, respectively.


Introduction
Although a significant amount of water is being consumed in the industrial and domestic sectors, the agricultural sector is considered as the largest consumer of water with 80% of water consumption worldwide [1]. Therefore, unless the current water management practices become dramatically wiser; many parts of the world will face rigorous competition for water among agriculture, energy, industry and civil activities [2].
Water scarcity, which refers to the lack of satisfactory available water resources to meet the water needs within a particular region, can be classified into two types, physical and economic [3]. The physical water scarcity, as per Srinivasan et al. [3], results from natural water resources that are inadequate to meet the demand of a certain area; however, the economic water scarcity results from poor handling of accessible water resources. Water scarcity is considered as one of the most critical problems facing many societies worldwide, and is directly interlinked to the food sector as over 80% of the world water withdrawal is to meet the requirements of the increasing population and the continuing development [4]. Due to water scarcity, many countries are forced to import food from abroad to meet the demand of its own people.
Water Footprint (WF), originated from the conception of Virtual Water (VW) [5], can be defined as being the total volume of freshwater used to produce the goods and services consumed by the individual or community or produced by the business, and is measured as volume of water consumed (evaporated) and/or polluted per unit of time [6,7]. Therefore, WF is an effective means that is used to quantify stress on water resources to address global, regional, national and local water scarcities. In this text, WF of a crop is defined as being the volume of freshwater used to produce a unit of a specific agricultural product (m 3 t -1 ) throughout the period of its production [7,8]. Crop WF provides an effective tool to investigate the linkage between food and water resources as a function of climate, soil and agricultural practices [9]. In general, WF analysis connects a wide range of sectors and issues and provides a multi-disciplinary scope for efficient management of water resources [10].
The three components of the WF that provide a thorough insight of the use of water through the demarcation of the consumed water are the green, the blue, and the grey WF. The blue WF (WF b ) refers to the consumption of blue water resources (surface and groundwater) by a product throughout its water demand stages. Water consumption is defined as being the loss of water from the available ground and/or surface water in a catchment area. Losses occur when water evaporates, returns to another catchment area (i.e. the sea) or is incorporated into a product [7]. The green WF (WF g ) is the consumption of green water resources (rainwater as far as it is not considered as run-off). The grey WF (WF gr ), however, refers to the pollution and is defined as being the volume of freshwater required to absorb the load of pollutant concentrations in a given natural environment.
The blue water resources are generally scarcer and costlier compared to the green water resources, which are valid reasons to place most of the emphasis and focus on only the blue water footprint. However, the green water resources are also limited and, therefore, are considered scarce, which provides an argument to account for the green WF as well. In addition, the green water can be substituted by the blue water and, in agriculture, the reverse can also take a place indicating that a complete knowledge can be obtained by accounting for both types of water.
After the International Expert Meeting on virtual water (VW) Trade held in December 2002 and the special session on VW Trade and Geopolitics during the Third World Water Forum held in March 2003, many studies have highlighted WF as a tool for global, regional and national water savings [8,[11][12][13][14][15][16], others have calculated the WF on a global scale [6,17]. Many studies have focused on WF and VW trade among countries with primary crops at subnational/local levels in conjunction with variations in climate, soil and other factors [18]. Studies have been accomplished, at the national level, on the use of water for crop production, municipal water consumption and the flow of water in countries, such as China [19,20], India [18], the Netherlands [21], the UK [22], Indonesia [23] and Nepal [24]. However, such studies are very limited in Saudi Arabia [25][26][27], and are not available even at the level of fields or farms, which are experiencing dynamic climates, scarce water reserves, and limited rainfall.
On the other hand, existing methods calculate the WF using data from national statistics, reports and climatic databases [7,19,24,25]. Remote sensing techniques provide the possibility of mapping the WF of a specific crop through the use of models, such as crop growth, productivity, and crop water use models. However, the use of remote sensing data in the quantification of WF is limited [28][29][30]. Therefore, this study was designed to bridge the gap in knowledge existed in the area of WF in the Kingdom of Saudi Arabia by quantifying and analyzing the spatial variation of WF of maize and carrot crops, cultivated during the period from December 2015 to December 2016, in the Eastern region of the Kingdom. The WF of silage maize and carrot crops was estimated in two ways, namely: agro-meteorological (WF Agro ) and remote sensing (WF RS ) methods. The green (WF Agro-g ), blue (WF Agro-b ) and grey (WF Agro-gr ) components of WF Agro were computed with the use of weather station/Eddy covariance measurements and field recorded crop yield datasets. The WF RS estimated by applying surface energy balance principles on Landsat-8 imageries.
Landsat images are widely used in ET mapping and in the estimation of the WF g and the WF b . However, the sparse and erratic rainfall conditions over the study site, resulted in a lack of Landsat-8 images coinciding with the rainy days. Hence, this study was limited to the blue component (WF RS-b ) of the WF RS approach. Moreover, the WF Agro estimates were used as a reference for the accuracy assessment of Landsat-8 derived WF RS-b in conjunction with the climatic conditions and cropping patterns.

Study area
The study was carried out in six 50 ha agricultural fields that were part of the 47 fields of the Todhia Arable Farm (TAF) located about 250 km Southeast of Riyadh, the capital city of Saudi Arabia, at coordinates of 24˚11 0 00@ E and 48˚56 0 14.6@ N (Fig 1). TAF was in an arid climate with hot summers (40 ± 2˚C) and cold to moderate winters (15 ± 3˚C), with a mean air temperature of 35˚C. The annual rainfall was about 90 mm, most of which occurred in the period from November to February. Due to the high crop water demand combined with highly erratic rainfall, irrigation is totally provided using groundwater delivered by center pivot irrigation systems. The major crops cultivated in the study area were forages (alfalfa, Rhodes grass and corn) and vegetable crops (carrot and lettuce). Depending on the demand, vegetable crops were cultivated throughout the year.

Field data
For the determination of WF of carrot and maize crops, a field survey was conducted to understand the cropping pattern of the experimental farm and to develop the sampling strategy (permission to conduct the study was issued by the Farm Manager, Mr. Alan King). Carrot crop was cultivated throughout the year; while, maize crop was cultivated twice a year (March-June, and July-November). Out of the 47 fields, maize crop was cultivated in 23 fields and carrot was cultivated in only seven fields during the period from December 2015 to December 2016. For this study, six center pivot irrigated fields were considered as sample fields, four fields (TE-2, TE-9, TE-11, PAL) were designated for silage maize and the remaining two for carrot (3)(4)(5). Datasets pertaining agricultural practices, such as sowing and harvesting dates, crop growth stages, amount of applied irrigation water, agro-climatic data and crop yields were obtained for the whole crop growth period. Except for the agro-climatic data which was extracted from an Eddy Covariance system installed in the farm, the data sets were taken from the records of the experimental farm.

Water footprint of agricultural crops
The methodological flow of determining the WF of agricultural crops is provided in Fig 2. The agro-meteorology (WF Agro ) and remote sensing (WF RS ) methods were applied to estimate the WF of silage maize and carrot crops. The WF Agro is computed with the use of empirical equations. Weather station/Eddy covariance system measured parameters (temperature, precipitation, evapotranspiration, etc.,) and the field measured crop yield datasets were used as input parameters. The WF Argo-g and WF Argo-b were considered as the amount of water consumed by the crop, while the WF Argo-gr is the amount of water needed for leaching requirement and to assimilate possible applied agro-chemicals [16].

Agro-meteorological (empirical) approach
Crop water requirement. Crop water requirement (CWR) was calculated by multiplying the reference crop evapotranspiration (ET o ) by the crop coefficient (K c ), as in Eq (1) described by Savva and Frenken [31]. In this study, the CWR was assumed to be equal to the actual crop evapotranspiration (ET c ) and there were no water limitations to crop growth so that the CWR was fully met. However, the crop type, variety and the developmental stage of the crop significantly affect the ET c . The K c values listed in FAO [32] were used in this study (Table 1).
where, ET c is the crop evapotranspiration (mm d -1 ), K c is the crop coefficient and ET o is the reference crop evapotranspiration (mm d -1 ).
Climatic variations, such as temperature, sunshine, wind speed, and humidity are the controlling factors for K c and ET o . The ET o (mm d -1 ) was estimated using the formula described in ASCE-EWRI [33] as shown in Eq (2).
where, R n is the measured net irradiance at the crop canopy (MJ m -2 d -1 ), G is the soil heat flux density (MJ m -2 ), T is the mean daily air temperature (˚C), U 2 is the mean daily wind speed at 2 m height (m s -1 ), e s is the saturated vapor pressure (kPa), e a is the mean actual vapor pressure (kPa), Δ is the slope of the saturation vapor pressure-temperature-pressure curve (kPa˚C -1 ) and γ is the psychometric constant (kPa˚C -1 ). Effective rainfall (P eff ) refers to the percentage of rainfall, which becomes available to crops. A number of factors can influence the effective rainfall including the amount and density of the rainfall, soil texture and the bulk density and the topography and the slope of the study area. The P eff was calculated using Eq (3) following Kuo et al. [34].
where, P eff is the effective rainfall (mm) and P tot is the total rainfall (mm). Soil salinity is a common aspect of irrigated agriculture over hyper-arid regions salinity due to evaporation of irrigation water. Therefore, to improve/protect soil quality, leaching out the accumulated salts from the soil profile is essential by applying an excess amount of water at the beginning of the growing season. The amount of water required for leaching, the leaching requirement (LR), is calculated from the Eq (5) as described in FAO [32].
where, EC w is salinity of the applied irrigation water (dS m -1 ), EC e is soil salinity measured from the soil saturation extract. In addition to LR, based on the amounts of CWR and P eff , the agricultural water requirement for each field was estimated using Eq (5) as described by Chowdhury et al. [35].
where, Q is the monthly agricultural water requirement of a field or irrigation scheme (m 3 d -1 ), i is the crop index, A i is the field area (ha), ET c is the crop evapotranspiration (mm d -1 ) and P eff is the effective rainfall (mm d -1 ). Fully automated sprinkler irrigation system was used to achieve high irrigation application efficiency. The efficiency of the sprinkler system was fixed at 70% as described in Al-Zeid et al. [36] for the application of irrigation water to the crops in the experimental fields.
Crop water use (CWU). The green and blue components of the CWU (m 3 ha -1 ) were calculated from the accumulation of the daily evapotranspiration (ET, mm d -1 ) over the entire growing period. They were calculated using Eqs (6) and (7) according to Chapagain and Hoekstra [37] as follows: where, lgp is the length of the growth period of the studied crop; ET green is the green water evapotranspiration and ET blue is the blue water evapotranspiration; CWU green is the total rainwater evaporated from the field during the growing period and CWU blue is the total irrigation water evaporated from the field.
ET green ¼ min ðET c ; P eff Þ ½length=time ð8Þ where, the ET green , i.e. evapotranspiration of rainfall water, is equal to the minimum difference between the total ET c and P eff . The ET blue , however, is equal to the maximum difference between the total ET c and P eff ). The grey assimilation water use (AWU grey , m 3 ha -1 ), the amount of water is used for the leaching (wash out) of salts from the rooting zone. The AWU grey was calculated using Eq (10), by dividing the pollutant load (L leached , in mass/time) by the difference between the ambient water quality standard for that pollutant (the maximum acceptable concentration C max , mg l -1 ) and its natural concentration in the receiving water body (C nat , mg l -1 ) [26].
The accumulated salts/pollutants generally consist of fertilizers, pesticides, and insecticides. In this study, the Nitrate (N) was the representative pollutant, and α was valued at 10% flat rate based on the study of Hoekstra and Chapagain [8]. Subsequently, C nat was considered as being equal to zero and C max as being equal to 11.5 mg l -1 as per the Saudi water quality standards [27]. The magnitude of N leaching, however, depends on soil conditions (irrigation frequency, rainfall pattern, soil texture, percolation rate, etc.) and methods of fertilizer application, including rate, time and agronomical practices [37]. Based on the chemical application rate per hectare (AR, kg ha -1 ) and the leaching fraction (α), the L leached can be estimated using Eq (11) according to Chapagain and Hoekstra [37].
Assessment of water footprint (WF Agro ). Water Footprint (WF) of both maize and carrot crops was calculated, based on the framework explained in Chapagain and Hoekstra [37], as the ratio of the total water used (m 3 ha -1 ) and the crop yield (t ha -1 ). The three components, WF Agro-b originated from irrigated water, WF Agro-g originated from rainwater and, the WF Agro-gr related to leaching of salts and wash out agro-chemicals [8,19] estimated and subsequently, the total water footprint (WF Agro ) of silage maize and carrot crops was calculated as in Eq (12). The obtained WF of silage maize and carrot crops was expressed in m 3 t -1 .
The WF g , which is the WF components produced from the green portion (i.e. rainfall) of crop water use, was calculated by dividing the component in crop water use (CWU green , m 3 ha -1 ) by the crop yield (Y, t ha -1 ) as in Eq (13). The WF b and the WF g were calculated in a similar way as in Eqs (14) and (15).
where, WF b and WF g are the WF components produced from the portion of the irrigated water and the chemicals assimilated water, respectively.

Remote sensing approach
A total of 32 cloud-free Landsat-8 satellite images for the period from December 2015 to December 2016 were downloaded from the portal of the USGS Earth Explorer (http:// earthexplorer.usgs.gov). The acquired images (Path 164 and 165; Row 45) covered the entire growth period of carrot and silage maize crops. In order to generate a reliable WF of carrot and silage maize crops, Landsat-8 data were analyzed for soil-adjusted vegetation index (SAVI), land surface temperature, a fraction of ET (ET f ) and actual ET (ET a ). Initially, Landsat-8 data were geo-referenced to the Universal Transverse Mercator (UTM) map projection with the World Geodetic System 84 (WGS84) datum. Subsequently, Landsat-8 digital numbers were converted to spectral radiance and transformed to Top-Of-Atmosphere (TOA) reflectance [38,39]. The spectral reflectance and land surface temperature products were generated by executing the "ATCOR" module of Geomatica software program (ver. 2015). Landsat-8 derived evapotranspiration (ET a ). Simplified surface energy balance (SSEB) model was used for estimating ET a (Eq 16) of carrot and silage maize in the study area. Predefined hot and cold conditions of each pixel were used, in Eq 15, to compute the ET f as described in Senay et al. [40].
where ET o is the reference ET (mm d -1 ), ET f is the ET fraction (unitless), and k is a scaling coefficient (unitless). The ET o was obtained from the weather datasets of the Eddy covariance (EC) system. In this study, k value was set to 1 as described in Senay et al. [40] and ET fraction (ET f ) was computed as per Eq 17.
where, T h is the reference hot pixel temperature (K), T s is the land surface temperature (K) obtained from Landsat-8 images, and T c is the reference cold pixel temperature (K). Satellite image analysis was carried out with the help of Geomatica software program (ver. 2015). Each individual Landsat-8 scene was processed separately for computing ET a on the day of satellite overpass. All the available daily ET a images were used for upscaling to the entire season level in order to compute ET a with respect to the length of the growth period (lgp) of each crop (Eq 18).
where, ET lgp is the crop specific lgp total ET (mm d -1 ), ET oi is the reference ET (mm) for period i (days) and ET fi is the representative ET f (unitless) for the period i. The mean monthly ET, derived from Landsat-8 products, was used in the estimation of the WF of carrot and silage maize crops. Landsat-8 estimated water use (ET a ) and crop yield were compared with the actual values (WF Agro ), which were obtained through the agrometeorological method.

Landsat-8 based water footprint (WF RS ).
After computing the ET lgp of each crop, based on the precipitation days and irrigation, the ET lgp was attempted to segregate into blue and green components of ET a . Due to the fact that the study area was low in its annual precipitation (~90 mm y -1 ), the green component of ET a was at its minimum values. The availability of satellite images on rainy days is limited. Hence, the Landsat-8 based WF estimation was limited to only the blue component of ET a (i.e. WF RS-b ). Subsequently, the WF RS-b was computed as the ratio between the crop water use (i.e. Landsat-8 estimated ET a of blue portion = CWU blue ) and the predicted yield (Y P ) as in Eq 14.

Accuracy assessment
The accuracy of Landsat-8 derived WF RS-b was assessed against the agrometeorological computed WF Agro-b . The performance indicators that were used for the accuracy assessment included Pearson correlation coefficient (R 2 ), root means square error (RMSE), mean bias error (MBE) and Nash-Sutcliff Efficiency (NSE).

Meteorological conditions
The average monthly meteorological data for the period from December 2015 to December 2016 are presented in Table 2. On the average, the monthly minimum, maximum and mean annual air temperatures were 19.3˚C, 34.0˚C and 26.5˚C, respectively. The total amount of rainfall during the study period was recorded at 13.6 mm. In addition, the average monthly wind speed was 5.0 m S -1 ; however, the average monthly ET o values during the study period ranged between 155 mm (November 2016) and 530 mm (July 2016), with an average daily ET o of 11.0 mm.

Electrical conductivity (EC) for soil and irrigation water
Electrical conductivity (EC) can provide accurate estimates of the number of salts presented in soil and water. The EC of both the soil and irrigation water correlates significantly with other agricultural field properties that affect crop productivity such as soil texture, cation exchange capacity (CEC), drainage conditions, organic matter level. To maintain the quality of soil, it is necessary to leach out salts from the root zone by means of additional irrigation water. For the determination of the Leaching factor (α) as in Eq (4), the soil EC (EC e ) and irrigation water EC (EC w ) were used as inputs ( Table 3). The values of EC e for carrot fields (3-5-S and 5-5-S) were determined at 1.43 (±0.33) dS m -1 and 2.76 (±1.11) dS m -1 , respectively. On the other hand, the EC e of silage maize fields varied from 2.27±0.48 dS m -1 (PAL field) to 5.21±1.26 dS m -1 (TE-2 field). The mean EC w of irrigation water in the experimental fields ranged between 1.48 dS m -1 (3-5-S field) and 2.09 dS m -1 (TE11 field).

Crop water requirements
Since groundwater is the main source of irrigation, crops are cultivated throughout the year depending on the demand and price. As part of agrometeorological estimation of WF Agro , the computation of crop water requirement (CWR) is essential for accurate scheduling of irrigation water. For the application of irrigation water through the sprinkler system, the irrigation interval and the amount of water applied to the experimental crops (carrot and silage maize) were calculated using CROPWAT software program (Ver. 8.0). As illustrated in Table 4, the CWR varies across the crops and lengths of growth period (lgp). The temporal dynamics of ET o and salinity of soil and water plays a key role in the variability of CWR. The CWR together with Leaching requirement (LR) for the summer grown silage maize crop, was estimated as high (1622 mm) as compared to spring grown the crop (1359 mm). Similarly, the CWR+LR of carrot crop grown in summer months of 3243 mm was higher than that of the winter-grown carrot of 620 mm. The agro-climatic variables, such as wind speed and temperature, significantly influenced the ET o and the lgp of a crop and the respective K c values, which were the main cause of the increase in CWR of crops in summer.

Crop yield
Silage maize and carrot yields predictive regression models were developed based on the values of the Soil Adjusted Vegetation Index (SAVI). The SAVI was computed for each crop throughout its growth period from seedling to seven days prior to the harvest. Season and crop wise datasets, about 6 to 7 Landsat-8 images of multiple dates were analyzed for this purpose (S1 File). Subsequently, the obtained SAVI is correlated with the field-recorded yields (Y A ) for the prediction of yield (Y P ). A linear relationship between the Landsat-8 derived SAVI and Y A of tested crops were developed, and the best-fit model was used in the prediction of Y P (Fig 3). A total of 100 samples, i.e. 40 (from two carrot fields) and 60 (from two silage maize fields) sampling locations were identified for each season and the Y A measurements was recorded. Of which, about 60 points (60%) used for the development of the Y P models, while the remaining 40% (40 samples) were used for the cross-validation of the models. The SPSS statistical software (Ver. 18.1) was used for the development and cross-validation of the Y P models.
As shown in Table 5, the obtained models were validated for their accuracy against the insitu yields using performance indicators, such as Pearson correlation coefficient (R 2 ), root  mean square error (RMSE), mean bias error (MBE) and Nash-Sutcliff Efficiency (NSE). The best relationship between crop yield and SAVI was obtained when the crops were in mid-stage of development. In the case of silage maize, the best response was observed on the Julian days of 162 (SAVI = 0.456) and 272 (SAVI = 0.462) for summer and spring seasons, respectively, when the crops were at their peak growth stage. The best response for carrot fields was observed on the Julian days of 68 (winter) and 146 (summer), when the values of SAVI were estimated at 0.416 and 0.398, respectively. The average predicted yield (Y P , DM t ha -1 ), for silage maize and carrot, was 31.98 and 37.65, respectively. The variation between Y A and Y P was found to be 15% (i.e. 4.5 t ha -1 ) for carrots and 17% (5.4 t ha -1 ) for silage corn.

Crop water use (CWU)
The CWU of carrot and silage maize crops was obviously observed to vary across the seasons. The pattern of the CWU (green, blue and grey) and the leaching fraction (pollutant load for the dilution of the applied nutrient salts) across the experimental fields were presented, along with the actually applied irrigation, in Table 6. Although the study farm was keen to achieve the best yields through effective management of irrigation water, there was a reduction ranging between 5.12% and 39.8% in actually applied irrigation water compared to the estimated CWR for the two experimental crops ( Table 4). The highest reduction was observed during the summer season, where 39.8% and 19% reduction occurred for silage maize and carrot, respectively. The green, blue and grey components of CWU (mean) for carrot fields were 0.4%, 68.1% and 31.8%, respectively. For silage maize crop, these values were 0.81%, 77.25% and 22.54%. In order to understand the spatial variation in WF, Landsat-8-based CWU RS-b (ET) of the two crops was predicted (Figs 4 to 7). The predicted CWU RS-b for silage maize was ranged from 712 (winter) up to 1108 (summer), while for carrot crop it was about 543 and 1382 mm for both summer and winter seasons, respectively. The RMSE (mm) between CWU Agro-b and CWU RS-b was ranged from 75.49 (carrot) to 234.37 (maize).

Water footprint (WF)
The patterns of WF (green, blue and grey) are presented in Table 7. Among the four silage maize investigated fields, filed number TE-11 was the lowest in WF Agro (3545 m 3 t -1 ); while the highest WF Agro (4960 m 3 t -1 ) was recorded for field number PAL. On the average, the contribution of WF Agro-g , WF Agro-b and WF Agro-gr for silage maize was estimated at 0.75%, 77.33% and 22.48%, respectively. The WF Agro of carrot crop ranged was ranged between 297 m 3 t -1 (field number 3-5N) and 502 m 3 t -1 (field number 5-5N), with an average WF Agro value of 396 m 3 t -1 over the entire study period. The WF Agro-g , WF Agro-b and WF Agro-gr of carrot were contributed with 0.6%, 66.9% and 32.9%, respectively. The remote sensing (RS) based WF RS-b varied across the crops from 276 (±73) m 3 kg -1 (carrot) to 2884 (±441) m 3 t -1 (silage maize). The RS based yield (Y P ), WF Agro and WF RS-b are provided in Table 7. The variation between WF RS-b and WF Agro-b was about 17% and 14% for silage maize and carrot, respectively. The WF Agro composed with the highest fraction of the WF Agro-b from 67% for carrot and 77% for maize crop. The grey water, however, corresponded to local salinity and crop salt tolerance and ranged from 22% (maize) to 33% (carrot). The current results for carrot crop indicated that the obtained WF was 30% lower than that reported by Multch et al [26]. Results of silage maize, however, indicated that the obtained WF was 2.5 times lower than that reported by Chowdhury et al. [36]. The total WF Agro values of silage maize (4014 m 3 t -1 ) and carrots (396 m 3 t -1 ) obtained in this study were lower than that of earlier studies (Table 8). However, the WF Agro-gr determined in this study was observed to be higher than the global averages stated in Mekonnen and Hoekstra [16] for both silage maize and carrot crops. Results of this study indicated that WF Agro was relatively higher than the global WF statistics [16]. This may be due to the fact that global statistics utilized both irrigated and arid crops for the compilation of global WF.
The empirical (agro-meteorological) approach based WF Agro-b , calculated from actual insitu data, was used as a reference against which the accuracy of the Landsat-8 determined WF RS-b was evaluated. A comparative analysis of the blue component of WF (WF b ) estimated by both WF Agro-b and WF RS-b approach is illustrated in Fig 8. Results revealed a highly significant linear relationship between the empirical (WF Argo-b ) and the Landsat-8 derived WF RS-b . In the case of silage maize, the recorded R 2 was 0.82 (P>F = 0.010) with the RMSE value of 501 m 3 t -1 (17%) and an MBE value of 218 m 3 t -1 (8%). For the carrot crop, the R 2 found to be  [40] found that the SSEB model was able to capture the seasonal ET a well with a strong correlation, where the R 2 value ranged between 0.76 (carrot) and 0.97 (silage maize) for the relationship between model estimated and actual ET. However, in terms of magnitude, the specific agreement in overestimation or underestimation of ET values may depend on the seasonal dynamics of climate and vegetation cover. Particularly, the extreme ET is usually reliable since the subjective selection of hot and cold reference pixels is eliminated and the only variable from season to season is the land surface temperature.
On the other hand, irrigation efficiency at the border areas of the center pivots depends mainly on the system efficiency along with the soil type, slope, and surface roughness. In this study, the boarder-effect is minimal as the fields are having slopes of less than 5%. However, the deposition of fine sand as a result of sand dunes and the prevailing winds from the surroundings of the "desert area" resulted in the variation of the WF at the border areas.

Conclusions
A field study was conducted to investigate the Water Footprint (WF) for carrot and silage maize crops cultivated in Saudi Arabia during the period from December 2015 to December 2016. The specific conclusions drawn from the study are as follows: • This study demonstrated the methods of estimating water footprint using Eddy covariance (WF Agro ) and remote sensing-based approach for WF mapping from Landsat-8 multispectral imagery.
• Due to the limited availability of satellite images on rainy days, the Landsat-8 based estimates were restricted to only the blue component of WF (i.e. WF RS-b ).  Utilization of Landsat-8 data for the estimation of crop water footprint • The utility of Landsat-8 data in mapping CWU RS-b showed reliable seasonal estimates of 1199 mm (summer) and 761 mm (spring/winter), which were in accordance with the Eddy covariance measured ET at 939 mm and 750 mm for summer and spring/winter growth periods, respectively.
• Among the six experimental fields, Landsat-8 determined WF RS-b of silage maize varied from 2526 to 3822 m 3 t -1 . For carrot, WF RS-b values were estimated at 242 to 302 m 3 t -1 .
• The Landsat-8 derived WF RS-b showed a highly significant linear relationship with the empirical WF Agro-b approach (R 2 = 0.77, P>F = 0.001).
• Feasible water footprint assessment system of agricultural crops for the efficient use of available water resources was developed by the integration of remote sensing technology (Landsat-8 satellite images) and weather data from the agro-meteorological station (Eddy Covariance system).