Next Article in Journal
Random Access Performance of Distributed Sensors Attacked by Unknown Jammers
Next Article in Special Issue
E2mC: Improving Emergency Management Service Practice through Social Media and Crowdsourcing Analysis in Near Real Time
Previous Article in Journal
Survey on the Performance of Source Localization Algorithms
Previous Article in Special Issue
How Well Can We Extract the Permanent Displacement from Low-Cost MEMS Accelerometers?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessing Crop Coefficients for Natural Vegetated Areas Using Satellite Data and Eddy Covariance Stations

1
Department of Civil and Environmental Engineering, Politecnico di Milano, Piazza Leonardo da Vinci, 32, 20133 Milan, Italy
2
Environmental Protection Agency of Aosta Valley, Climate Change Unit, 11100 Aosta, Italy
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(11), 2664; https://doi.org/10.3390/s17112664
Submission received: 30 September 2017 / Revised: 7 November 2017 / Accepted: 9 November 2017 / Published: 18 November 2017

Abstract

:
The Food and Agricultural Organization (FAO) method for potential evapotranspiration assessment is based on the crop coefficient, which allows one to relate the reference evapotranspiration of well irrigated grass to the potential evapotranspiration of specific crops. The method was originally developed for cultivated species based on lysimeter measurements of potential evapotranspiration. Not many applications to natural vegetated areas exist due to the lack of available data for these species. In this paper we investigate the potential of using evapotranspiration measurements acquired by micrometeorological stations for the definition of crop coefficient functions of natural vegetated areas and extrapolation to ungauged sites through remotely sensed data. Pastures, deciduous and evergreen forests have been considered and lower crop coefficient values are found with respect to FAO data.

1. Introduction

Accurate assessment of evapotranspiration is fundamental for efficient water resource management, for the design and management of water supply reservoirs, for designing and scheduling irrigation systems, and for environmental assessment [1]. Moreover, trends of modifications of water availability around the world, have prompted scientific investigations of the effect of climate change on the hydrological cycle [2,3,4,5]. In order to study the effects of climate change for several years in advance, the use of parsimonious models for the assessment of evapotranspiration is needed in order to maintain the simulation time within acceptable limits and because of limited availability of meteorological data at a spatial and temporal resolution sufficient to accurately capture the dynamics of hydrological processes.
Latent heat flux (LE), i.e., evapotranspiration (ET), is the joint term that rules the relationship between water and energy fluxes through vegetated soil and shallow atmosphere layer, since it is the only term present both in the water as well as in the energy balance. Despite its important role, direct measurements of ET are neither frequent nor simple, due to the complexity of the process, its spatial scale of representativeness and plant phenology. These issues often affect the reliability of direct measures of water losses, so that the scientific community proposed in the late forties onwards several methods to estimate ET based on meteorological measures [6]. A simplified method is based on the crop coefficient (kc) which embodies all the physiologic characteristics of a specific plant allowing one to relate a reference evapotranspiration of well irrigated grass to the potential evapotranspiration of each crop [7]. The kc is affected by growth stage, crop types, canopy conductance, soil characteristics, crop height and leaf area index.
This approach, which has been adopted by the Food and Agricultural Organization (FAO), is based on a revised version of the Penman-Monteith equation. Alternatives were proposed for reference evapotranspiration assessment that use Hargreaves-Samani equation [8] and its derivatives [9,10,11].
A database of kc values for a large number of agricultural crops has been created by [7] for different climates and [7] presented procedures for estimating kc values for natural vegetation as a function of leaf area index. However, local ground measurements are still needed to estimate a local corrected kc to account for local crop varieties and for specific weather conditions [12].
Potential evapotranspiration, which was usually measured with a weighting lysimeter as described in FAO paper no. 56, suffers from important limitations, and can also be defined with more innovative techniques. In particular micro-meteorological stations, based on eddy covariance, are regarded to be accurate enough for ET estimates and are now widely diffused [13,14,15,16] in cultivated fields as well as in natural vegetated areas.
Many works have been done in the last decades in order to calculate the crop coefficients in specific areas (more numerous in semi-arid and arid regions) for a large variety of crops [17,18,19,20,21]. The kc in these papers are locally obtained from lysimeter or eddy covariance measurements and result, in general, values that are lower than the FAO paper no. 56 ones.
Although the kc method has been widely used for crops, it has not been widely examined for natural ecosystems, such as forests and other perennial vegetation [22,23]. A few examples are available from [24] who analyzed small shrubs founding a kc between 0.50 and 0.85 while for well-developed shrubland a kc between 0.85 to 0.95 was determined. Reference [23] analyzed different vegetation covers obtaining that, except for evergreen forests, kc values have large seasonal variation as a function of crop growth and mainly precipitation.
Moreover in the last few years remote sensing data are also always more frequently used in support of evapotranspiration estimates defining vegetation parameters and crop coefficients in a distributed way for agricultural crops with empirically- or physically-based relationships [25,26,27,28,29,30]. The greater advantage of the use of remote sensing data is the possibility to have spatially distributed information. In particular, methods that combine vegetation index and crop coefficients have been developed in agricultural crops from 30 years ago [31].
Two different objectives are defined in this paper:
(a)
The definition of crop coefficient curves for natural area derived from eddy covariance data to be used in hydrological modelling to compute effective evapotranspiration.
(b)
Assessing the reliability and potentiality of using satellite data and conventional meteorological measurements for crop coefficient estimates in natural vegetated areas where eddy covariance stations are not available.
In this paper the proposed procedure is applied for three types of natural coverage: pastures, deciduous and evergreen forests.

2. Materials and Methods

Much literature discusses “potential” and “actual” ET, the former being the maximum evapotranspiration which would occur if no factors other than energy supply and atmospheric demand limiting the ET rate [32]. Conversely, “actual” ET is the amount of water evapotranspirated in the case that water availability becomes a limiting factor [33].
Potential evapotranspiration can be derived from ground measurements, lysimeter or eddy covariance stations, or computed using a wide range of equations that are very well assessed in the research community, such as the Prietsley-Taylor equation [34], Hargreaves [8] or Penman Monteith [6]. In this study the Penman-Monteith equation has then been applied.

2.1. Penman-Monteith Equation

Penman-Monteith’s equation, that combines the energy balance with the mass transfer method, allows computing a potential evapotranspiration (ETPPenMon (mm·day−1) as:
E T P P e n M o n = [ Δ ( R n G ) + c p ρ a ( e 0 e a ) / r a ] Δ + γ ( 1 + r c / r a )
where Rn is the net radiation (MJ·m−2·d−1), G is the soil round heat flux (MJ·m−2·d−1), (e0ea) is the vapour pressure deficit of the air (kPa), ρa is the mean air density (kg·m−3), cp is the specific heat of the air (MJ·kg−1·°C−1), Δ is the slope of the saturation vapour pressure (kPa·°C−1), γ is the psychrometric constant (kPa·°C−1) and rc is the crop resistance (s·m−1) and ra is the aerodynamic resistance (s·m−1).
In particular net radiation is computed as sum of longwave and shortwave radiation:
R n = ( 1 α ) × R s R n l
where α is albedo, Rs is the incoming shortwave radiation (MJ·m−2·d−1) and Rnl is the net longwave radiation (MJ·m−2·d−1).
The aerodynamic resistance is a function of wind velocity (u) (m·s−1) and canopy height (hc) (m) [35]:
r a = ln ( z u 0.667 h c 0.123 h c ) ln ( z t 0.667 h c 0.0123 h c ) 0.168 × u
where zu and zrh are respectively the measurement heights of wind velocity and relative air humidity (m).
The crop resistance is computed as [36]:
r c = r s min 0.5 × L A I
where rsmin is the minimum stomatal resistance (s·m−1) and LAI is leaf area index. G is approximated as a fraction (0.1 during daytime [7]) of net radiation.

2.2. Sensitivity Analysis

A sensitivity analysis on kc is performed to understand how much it deviates from its mean value. A sensitivity coefficient is computed in respect to kc for each parameter included in crop coefficient estimates (e.g., LAI, albedo and vegetation height) and then is adimensionalized so that the results are comparable. The coefficient of sensitivity is computed as:
S K c = lim Δ k c 0 ( Δ V i / V i Δ k c / k c ) = V i k c k c V i
where SKc is the sensitivity coefficient and Vi is the ith parameter. A positive value of SKc means that the parameter will increase as kc increases and the opposite behaviour for negative values of SKc.

2.3. Eddy Covariance Technique

The eddy covariance method determines the surface fluxes as the sum of turbulent fluxes, measured above the surface, and the flux divergence between the surface and the eddy covariance measurement level [37]. The basic equations to estimate latent and sensible heat fluxes are:
L E = λ ρ w H 2 O ¯
H = ρ C p w T ¯
where λ is the vaporization latent heat, ρ the air density and w H 2 O ¯ the covariance between vertical wind velocity and concentration of vapor in the air. Cp is the specific heat at constant pressure and w T ¯ is the covariance between vertical wind velocity and temperature. The symbol (ʹ) indicates instantaneous fluctuation from the time-averaged values of a specific variable in according with Reynolds decomposition of a meteorological stochastic signal [38].
The reliability of flux measurements depends on different theoretical assumptions of the eddy covariance technique [17,39], among which the most important are the horizontal homogeneity, the stationarity and the mean vertical wind speed equal to zero during the averaging period. A number of corrections needs to be applied to obtain high quality fluxes flowing procedures which are now well assessed in literature [40,41].

2.4. The FAO Crop Coefficient

According to FAO Paper no. 56 [7], the crop coefficient is given by:
k c = E T P E T 0
where ET0 is the reference potential evapotranspiration (mm·day−1) and ETP is a generic potential evapotranspiration (mm·day−1). According to FAO Paper no. 56, three values of kc are defined over time in accordance to the crop development stage: (i) initial stage (kc-ini), (ii) mid-season stage (kc-mid) and the harvest stage (kc-end).
The reference evapotranspiration, ET0, is a potential evapotranspiration which derives from Penman-Monteith’s equation and is defined for a reference surface defined as an “hypothetical crop with an assumed height of 0.12 m having a surface resistance of 70 s·m−1 and an albedo of 0.23, closely resembling the evaporation of an extension surface of green grass of uniform height, actively growing and adequately watered” [7].
ET0 (mm·day−1) is computed as:
E T 0 = 0.408 Δ ( R n G ) + γ 900 T + 273 u 2 ( e 0 e a ) Δ + γ ( 1 + 0.34 u 2 )
where u 2 is wind speed at 2 m above ground surface (m·s−1) expressed as:
u 2 = 4.87 ln ( 67.8 z u 5.42 ) u
For natural vegetation [7] proposed to assess crop coefficient as a function of leaf area index (LAI)
k c F A O = k c , min + ( k c , max k c , min ) ( 1 e 0.7 L A I )
where:
k c , max = k c , h + [ 0.04 ( u ¯ 2 2 ) 0.004 ( R H min 45 ) ] ( h 3 ) 0.3
k c , h = 1 + 0.1 h
where u ¯ 2 is the mean value for wind speed at 2 m height during the mid-season (m·s−1), R H min is the mean value for minimum daily relative humidity during mid-season (%), and h is the mean maximum plant height (m).

2.5. Crop Coefficients

In this paper, crop coefficient is assessed in three ways that are different as regards ETP estimates:
k c e d d y = E T P e d d y E T 0
k c P e n M o n , s a t = E T P P e n M o n , s a t E T 0
k c P e n M o n , g r o u n d = E T P P e n M o n , g r o u n d E T 0
where ETPeddy is the potential evapotranspiration derived by eddy covariance stations, ETPPenMon,sat is the potential evapotranspiration estimated with Penman-Monteith’s equation using leaf area index and albedo retrieved from remotely sensed images and ground measured meteorological forcings, and ETPPenMon,ground is the potential evapotranspiration estimated with Penman-Monteith’s equation using data coming from ground measurements.
Since the crop coefficient approach proposed by FAO is intended for potential evapotranspiration estimate, kceddy is significant only when vegetation can evapotranspirate with no limitation in unlimited water conditions. As it can be expected that, when available moisture of the soil is enough, ET from eddy covariance station should be similar to the one calculated with the Penman-Monteith’s equation (Equation (1)), in the present work kceddy is computed only when the condition E T P e d d y E T P P e n M o n , g r o u n d 1 (±0.05) is satisfied.
In addition to above mentioned methods for crop coefficient assessment, Equation (12) is applied (kcFAO). A schematic of the methodology followed in this research is illustrated in Figure 1.

2.6. Sites and Data

Monitoring took place in four sites representative of three natural ecosystems: the first one is located in a grass field in the mountain area of Torgnon in Italy [45°50′40′′ N, 7°34′41′′ E, 2160 m above sea level], the second, Chestnut Ridge in the U.S. [35.9′ N, 84.3′ W, 286 m a.s.l.], and the third, Duke Forest in the U.S. [35.9′ N, 79.1′ W, 168 m a.s.l.], are located in the deciduous broadleaf forest region, while the fourth is in the Black Hills in the U.S. [44.1′ N, 103.6′ W, 1718 m a.s.l.] characterized by an evergreen forest.

2.6.1. Torgnon

The study site of Torgnon [42,43] is an abandoned pasture located in northwestern Italian Alps (Aosta Valley, Italy) at an elevation of 2160 m a.s.l. The site is part of the FLUXNET network. This site is characterized by an intra-Alpine semi-continental climate, with an annual mean temperature of 3.1 °C and a mean annual precipitation of about 880 mm. From the end of October to late May, the site is covered by a thick snowcover (90–120 cm) which limits the growing period, and hence plants evapotranspiration, to an average of five months. Dominant vegetation is composed by matgrass (Nardus stricta L., Festuca nigrescens All., Arnica montana L., Carex sempervirens Vill., Geum montanum L., Anthoxanthum alpinum L.L., Potentilla aurea L., Trifolium alpinum L..) and the soil is classified as Cambisol (WRB classification).
The eddy covariance method is used to measure the fluxes of H2O between the ecosystem and the atmosphere. Measurement of wind speed in the three components is performed by a CSAT3 three-dimensional sonic anemometer (Campbell Scientific, Logan, UT, USA), while H2O vapor air densities were measured by a LI-7500 open-path infrared gas analyzer (LICOR, Lincoln, NE, USA). Instruments were placed at 2.5 m above the ground. Data from both the anemometer and the IRGA were measured at 10 Hz. A weather station is installed close to the eddy covariance tower to continuously measure additional meteorological variables. Air and soil temperatures are measured respectively by a HMP45 (Vaisala, Helsinki, Finland) and with temperature probes type therm107 (Campbell Scientific) at different depths. Soil water content is assessed with soil water reflectometers model CS616 (Campbell Scientific). The variables of interest for this study, i.e., soil heat flux and net radiation are measured respectively by HFP01 plates (Hukseflux, Delft, The Netherlands) and by a CNR4 (Kipp & Zonen, Delft, The Netherlands) net radiometer. Snow height is measured with a SR50A sonic snow depth sensor (Campbell Scientific).
Since the grassland is unmanaged, canopy structure varies during the season following the phenological development. Measurements of canopy height and Leaf Area Index (LAI) at site are available for the study years, in detail: canopy height is measured and phytomass is sampled every two weeks at 12 selected plots of 40 cm × 40 cm in the study area. Total LAI and canopy height for the 12 plots are averaged to obtain a site-level value.

2.6.2. Chestnut Ridge

Chestnut Ridge is a site included in the FLUXNET network (ORNL DAAC, 2011) located in Tennessee (USA) in a temperate climate and managed by Oak Ridge National Laboratory in Tennessee. The vegetation type below the station is deciduous broadleaf forest and in particular oaks and hickories are present. The tower height is 60 m. The station is active from 2006 and the data analysed in this paper are from the year 2007.
The micro-meteorological station is supplied with the following sensors: a Young-81000 3D sonic anemometer (Young, Traverse City, MI, USA), a LI-7500 open path gas analyzer (LICOR), a CNR1net radiometer (Kipp & Zonen), a Gcap Temperature/Humidity Probe (ATDD/NOAA, College Park, MD, USA), soil moisture sensor delta-T probes at different depths, Soil Temperature Probe and a HFP01heat flux plate (Hukseflux) and a TB3 Tipping Bucket Rain Gauge (Hydrological Services, Lake Worth, FL, USA). Data are available every 30 min. A more detailed description of this station can be found in [44,45].

2.6.3. Duke Forest

The Duke Forest C-H2O Research Site is located at the Blackwood Division of the Duke Forest near Durham (North Carolina, USA) and is part of the FLUXNET network. The forest is characterized by oak-hickory trees of about 33 m. The tower, which is 40 m tall, is equipped with a CSAT3 3D sonic anemometer (Campbell Scientific), a LI-6262 open path gas analyzer (LICOR), a NR-LITE net radiometer (Kipp & Zonen), a HMP35C Temperature/Humidity Probe (Vaisala) and soil moisture sensor delta-T probes at different depths. Data are available every 30 min and data from 2001 and 2002 are used in this paper. The papers that describe the setup of the station are [46,47].

2.6.4. Black Hills

The Black Hills eddy covariance station is located in South Dakota (USA) in a temperate climate and is managed by Oak Ridge National Laboratory in Tennessee (USA). It is part of the FLUXNET network. The tower is 24 m height and the vegetation coverage is composed by evergreen needleleaf forest, in particular ponderosa pines. The station is equipped with a R3sonic anemometer (Gill), a LI-7500 open path gas analyzer (LICOR), a CNR1 net radiometer (Kipp & Zonen), a HMP50Y Temperature/Humidity Probe (Vaisala), a HFP01 soil heat flux plate (Hukseflux), a ThetaProbe ML2 soil moisture sensors (Delta-T, Cambridge, UK) and soil temperature probes (ATDD/NOAA). A more detailed description of this station can be found in [47]. Data are available every 30 min and data from 2006 and 2007 are used in this paper.

2.6.5. Satellite Data

Leaf area index and albedo are retrieved from satellite data measured from MODIS radiometers on board TERRA and AQUA NASA satellites (https://earthdata.nasa.gov/). These satellites fly over the same area twice a day, providing a night and a day image. The values of LAI [48] and albedo [49] are available products which are averaged over a range of eight days to avoid problems of cloudy conditions. These parameters can be considered constant during this period of time.

3. Results and Discussion

kc values are computed at daily scale from measured eddy covariance data at 30 min or 1 h. A day of data is considered only if at least the 80% of daytime data is available (when Rn > 0), otherwise the entire day is deleted. Other gaps are also present in the observed series due to the identification of potential evapotranspiration conditions. In fact days characterized by soil water stress conditions are eliminated from the data series.
Finally in order to compute the difference between the computed crop coefficient values, RMSE is evaluated as:
R M S E = [ i = 1 n ( X s i m i X o b s i ) 2 n ] 0.5
where Xsimith is the ith simulated variable, Xobsith is the ith measured variable, n the sample size, and X o b s ¯ the average observed variable.

3.1. Results of the Sensitivity Analysis

In Figure 2 the sensitivity coefficients for kc are reported for LAI, albedo and vegetation height considering a variation from the mean value of ±50%, using the meteorological variables to compute energy fluxes measured by an eddy covariance station located in a maize field in Northern Italy [50].
Albedo and vegetation height have a low influence on crop coefficient estimates; in fact SKc remains almost constant and equal to 0.1 for variations of albedo and vegetation height from −50% to +50%. Instead the crop coefficient is highly influenced by the variability of leaf area index and in particular decreases in exponential way from −50% variation of LAI with SKc equal to 0.44 till values of 0.13 with a variation of +50% of LAI, so that the crop coefficient will be more sensitive during the growing period of the vegetation. These results are also confirmed by [51] who found that the most sensitive parameter in ET0 estimation is rc.

3.2. Pasture: Torgnon

In Figure 3 the computed crop coefficient values are reported for 2009 and 2010. The selected data belong to snow-free periods from day 183 to 305 for 2009 and from 142 to 302 for 2010. Results for kc are summarized in Table 1, while in Table 2 the lengths of the growth stages are reported.
During the year 2009, two crop coefficient values from eddy covariance measurements can be identified: the first one is relative to a completely developed grass with a mean values of 0.8, while the second one is relative to the senescence phase when grass starts to become yellow with a mean value of 0.47. A similar distribution is found also during 2010, with the first kceddy equal to 0.88 and the second to 0.43.
When the crop coefficient is computed with Penman-Monteith’s equation, a kcPenMon,ground of 0.76 is obtained during 2010 for the first period and 0.39 for the senescence period. These values are quite similar to those obtained from eddy covariance data, showing that for pasture correct kc values can be defined using meteorological data from basic network stations which are more available than from eddy covariance towers.
The analyses are then performed using satellite information of LAI and kcPenMon,sat are equal to 0.76 for the first period and 0.49 for the second period during 2009, while 0.76 and 0.36 for 2010. kcPenMon,sat are similar to kcPenMon,ground and kceddy, showing that this technique can be applied in a distribute way in river basins coupling basic network stations for meteorological data and satellite data for vegetation information.
Finally, RMSE values have been computed between kcPenMon,sat and kceddy and RMSE values equal to 0.15 and 0.19 are found for 2009 and 2010 respectively. When kcPenMon,ground and kceddy are compared for 2010, a RMSE of 0.18 is obtained.

3.3. Deciduous Forest

3.3.1. Chestnut Ridge

In Figure 4 crop coefficient values are reported for the deciduous forest of Chestnut Ridge. The data belong to the whole year, so that the complete seasonal cycle can be analyzed differencing the behaviors during periods with or without leaves (Table 1). In Table 2 the lengths of the growth stages are reported.
During the winter period, kceddy has low values equal to 0.19. Similar results are obtained during the autumn season when the leaves are not present as well. kc from eddy covariance tower is then equal to 0.2. When these periods are analyzed using satellite information, much lower values are obtained with kcPenMon,sat for the initial period equal to 0.04 and kc-end to 0.05 for both years. This behaviour is linked to values near zero or zero of LAI from satellite during these winter–autumn periods.
The growth and decrease branches of the two kc curves are quite in accordance, according to Figure 4, as well as for the summer period when kc-mid from eddy covariance station is equal to 0.47 and from satellite data to 0.48. Overall, analyzing the estimates accuracy for the whole year, RMSE obtained from the comparison between kc computed from eddy covariance data and from satellite data is found to be equal to 0.18.

3.3.2. Duke Forest

Crop coefficient results for the years 2001 and 2002 are reported in Figure 5. kceddy during 2001 is low during winter with a mean value of 0.11 while kcPenMon,sat is 0.07. Similar values are obtained for 2002 with kcPenMon,sat being 0.04 and kceddy 0.09. These results are obtained also during autumn (Table 1). During the summer season, when trees are completely lush, as expected, crop coefficient values are higher and equal to 0.51 from eddy covariance data and to 0.44 from satellite data for 2001; while for 2002 kc values of 0.43 and 0.44 are reached for kceddy and kcPenMon,sat respectively. Results for the two analyzed year are quite in accordance with each other. The lengths of the growth stages are similar between the two analyzed years (Table 2).
When the yearly RMSE is evaluated between kc from eddy covariance data and from Penman-Monteith’s equation with satellite data, values of 0.05 for 2001 and of 0.04 for 2002 are obtained.

3.3.3. Intercomparison

A comparison between Chestnut Ridge and Duke Forests results is performed in order to define a general crop coefficient curve which is representative for deciduous forests. According to Table 1, results of crop coefficient from eddy covariance data are quite similar between the two stations for the initial stage value for all the analysed year ranging between 0.09 and 0.19, as well as for the final stage with values between 0.12 and 0.2. kc-mid is also quite defined with values between 0.43 and 0.48.

3.4. Evergreen Forest: Black Hills

In Figure 6 crop coefficient values are reported for 2006 and 2007 for the Black Hills station. Even though the station is located in an evergreen forest, evapotranspiration substantially decreases during the winter period and so the crop coefficient does.
In Table 1 kceddy has low values equal to 0.05 during the winter period for both years. Similar valus are obtained for kcPenMon,sat. The final stage has similar low values equal to 0.03 and 0.07 during 2006 from eddy covariance data and from satellite data respectively, while for 2007 kceddy is equal to 0.03 while kcPenMon,sat to 0.04.
The growth and decrease branches of the two kc curves are quite in accordance, according to Figure 4 and Table 2, as well as for the summer period when kc-mid from eddy covariance station is equal to 0.17 and from satellite data to 0.2 for 2006 and to 0.18 and 0.2 respectively for 2007. Overall, analyzing the estimates accuracy for the whole year, RMSE between kc from eddy covariance data and from satellite data is equal to 0.07 in 2006 and to 0.05 in 2007.

4. Conclusions

In this paper, the potentiality of using evapotranspiration measurements from micro-meteorological stations for definition of crop coefficient functions of natural vegetated area has been evaluated. Sensitivity analysis shows that albedo and vegetation height have a low influence on crop coefficient estimates, while these are highly influenced by the variability of leaf area index.
For pasture, crop coefficient obtained from eddy covariance data are similar to the values obtained with ground measured data and those obtained using leaf area index and albedo retrieved from remotely sensed images. This leads to two important results: (1) measurements from widespread available stations provide robust estimate of crop coefficient, and (2) satellite data are powerful to compute crop coefficient in a distribute way in river basins wherever basic ground measured meteorological data are available.
The potentials of remotely sensed data have been confirmed using data available from the FLUXNET network both on deciduous and evergreen forests. The crop coefficient values for the considered coverage of pastures, deciduous and evergreen forests are found to be lower than the FAO values during the mid-season stage, showing the necessity to check FAO data before applying them to areas with climatic and vegetation characteristics different form the ones considered by FAO.

Acknowledgments

This work was supported by SIM project “Smart Irrigation from Soil Moisture Forecast Using Satellite and Hydro-Meteorological Modelling ” [http://www.sim.polimi.it], and INNOMED project “Innovative Options for Integrated Water Resources Management in the Mediterranean” funded by the Ministry for Education, University and Research of Italy within the EU Water Joint Programming Initiative. Authors thank Stefano Morella, Paride Gianoli for their help in elaborating data. Thanks also to G. Katul, R. Oren, T. Meyers for providing Fluxnet data.

Author Contributions

C.C., G.R., M.M. and E.C. conceived and designed the experiments; C.C., G.R. and M.G. performed the experiments; C.C. and G.R. analyzed the data; C.C., G.R. and M.G. wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

References

  1. Xu, C.Y.; Singh, V.P. Evaluation of three complementary relationship evapotranspiration models by water balance approach to estimate actual regional evapotranspiration in different climatic regions. J. Hydrol. 2005, 308, 105–121. [Google Scholar] [CrossRef]
  2. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration—Guidelines for Computing Crop Water Requirements—FAO Irrigation and Drainage Paper 56; Food and Agriculture Organization of United Nations: Rome, Italy, 1998. [Google Scholar]
  3. Beniston, M. Climatic Change and Its Impacts. An Overview Focusing on Switzerland; Kluwer Academic Publishers: Dordrecht, The Netherlands, 2004. [Google Scholar]
  4. Bocchiola, D.; Diolaiuti, G.; Soncini, A.; Mihalcea, C.; D’Agata, C.; Mayer, C.; Lambrecht, A.; Rosso, R.; Smiraglia, C. Prediction of future hydrological regimes in poorly gauged high altitude basins: The case study of the upper Indus, Pakistan. Hydrol. Earth Syst. Sci. 2011, 15, 2059–2075. [Google Scholar] [CrossRef] [Green Version]
  5. Diolaiuti, G.A.; Maragno, D.; D’Agata, C.; Smiraglia, C.; Bocchiola, D. Glacier retreat and climate change: Documenting the last 50 years of Alpine glacier history from area and geometry changes of Dosdè Piazzi glaciers (Lombardy Alps, Italy). Prog. Phys. Geogr. 2011, 35, 161–182. [Google Scholar] [CrossRef]
  6. Nijssen, B.; O’Donnell, G.M.; Hamlet, A.F.; Lettenmaier, D.P. Hydrologic Sensitivity of Global Rivers to Climate Change. Clim. Chang. 2001, 50, 143–175. [Google Scholar] [CrossRef]
  7. Monteith, J.L. Evaporation and environment. Symp. Soc. Exp. Biol. 1965, 19, 205–234. [Google Scholar] [PubMed]
  8. Hargreaves, G.H.; Samani, Z.A. Reference crop evapotranspiration from temperature. Appl. Eng. Agric. 1985, 1, 96–99. [Google Scholar] [CrossRef]
  9. Ravazzani, G.; Corbari, C.; Morella, S.; Gianoli, P.; Mancini, M. Modified Hargreaves-Samani equation for the assessment of reference evapotranspiration in Alpine river basins. J. Irrig. Drain. Eng. 2012, 138, 592–599. [Google Scholar] [CrossRef]
  10. Trajkovic, S. Temperature-based approaches for estimating reference evapotranspiration. J. Irrig. Drain. Eng. 2005, 131, 316–323. [Google Scholar] [CrossRef]
  11. Vanderlinden, K.; Giraldez, J.V.; Van Meirvenne, M. Assessing Reference Evapotranspiration by the Hargreaves Method in Southern Spain. J. Irrig. Drain. Eng. 2004, 130, 184–191. [Google Scholar] [CrossRef]
  12. Pereira, L.S.; Allen, R.A.; Smith, M.; Raes, D. Crop evapotranspiration estimation with FAO56: Past and future. Agric. Water Manag. 2015, 147, 4–20. [Google Scholar] [CrossRef]
  13. Baldocchi, D.; Falge, E.; Gu, L.; Olsen, R.; Hollinger, D.; Running, S.; Anthony, P.; Bernhofer, C.; Davis, K.; Evans, R.; et al. FLUXNET: A new tool to study the temporal and spatial variability of ecosystem scale carbon dioxide, water vapor, and energy flux densities. Am. Meteorol. Soc. 2001, 82, 2415–2434. [Google Scholar] [CrossRef]
  14. Foken, T. Micrometeorology; Springer: Berlin, Germany, 2008; p. 306. [Google Scholar]
  15. Wilson, K.; Goldstein, A.; Falge, E.; Aubinet, M.; Baldocchi, D.; Berbigier, P.; Bernhofer, C.; Ceulemans, R.; Dolman, H.; Field, C.L.; et al. Energy balance closure at FLUXNET sites. Agric. For. Meteorol. 2002, 113, 223–243. [Google Scholar] [CrossRef]
  16. Papale, D.; Reichstein, M.; Aubinet, M.; Canfora, E.; Bernhofer, C.; Kutsch, W.; Longdoz, B.; Rambal, S.; Valentini, R.; Vesala, T.; et al. Towards a standardized processing of net ecosystem exchange measured with eddy covariance technique: Algorithms and uncertainty estimation. Biogeosciences 2006, 3, 571–583. [Google Scholar] [CrossRef]
  17. Bezerra, B.G.; Bergson, G.; da Silva, B.B.; Bezerra, J.R.C.; Sofiatti, V.; dos Santos, C.A.C. Evapotranspiration and crop coefficient for sprinkler-irrigated cotton crop in Apodi Plateau semiarid lands of Brazil. Agric. Water Manag. 2012, 107, 86–93. [Google Scholar] [CrossRef]
  18. Casa, R.; Russell, G.; Lo Cascio, B. Estimation of evapotranspiration from a field of linseed in central Italy. Agric. For. Meteorol. 2000, 104, 289–301. [Google Scholar] [CrossRef]
  19. Čereković, N.; Todorović, M.; Snyder, R.L.; Boari, F.; Pace, B.; Cantore, V. Evaluation of the crop coefficients for tomato crop grown in a Mediterranean climate. Opt. Méditerr. A 2010, 95, 91–94. [Google Scholar]
  20. Facchi, A.; Gharsallah, O.; Corbari, C.; Masseroni, C.; Mancini, M.; Gandolfi, C. Determination of maize crop coefficients in humid climate regime using the eddy covariance technique. Agric. Water Manag. 2013, 130, 131–141. [Google Scholar] [CrossRef]
  21. Piccinni, G.; Ko, J.; Marek, T.; Howell, T. Determination of growth-stage-specific crop coefficients (kC) of maize and sorghum. Agric. Water Manag. 2009, 96, 1698–1704. [Google Scholar] [CrossRef]
  22. Glenn, E.P.; Neale, C.M.U.; Hunsaker, D.J.; Nagler, P.L. Vegetation index-based crop coefficients to estimate evapotranspiration by remote sensing in agricultural and natural ecosystems. Hydrol. Process. 2011, 25, 4050–4062. [Google Scholar] [CrossRef]
  23. Liu, C.; Sun, G.; McNulty, S.G.; Noormets, A.; Fang, Y. Environmental controls on seasonal ecosystem evapotranspiration/potential evapotranspiration ratio as determined by the global eddy flux measurements. Hydrol. Earth Syst. Sci. 2017, 21, 311–322. [Google Scholar] [CrossRef]
  24. Irmak, S.; Kabenge, I.; Rudnick, D.; Knezevic, S.; Woodward, D.; Moravek, M. Evapotranspiration crop coefficients for mixed riparian plant community and transpiration crop coefficients for Common reed, Cottonwood and Peach-leaf willow in the Platte River Basin, Nebraska-USA. J. Hydrol. 2013, 481, 177–190. [Google Scholar] [CrossRef]
  25. D’Urso, G.; Menenti, M. Mapping crop coefficients in irrigated areas from Landsat TM images. Proc. SPIE 1995, 2585, 41–47. [Google Scholar]
  26. Er-Raki, S.; Chehbouni, A.; Guemouria, N.; Duchemin, B.; Ezzahar, J.; Hadria, R. Combining FAO-56 model and ground-based remote sensing to estimate water consumptions of wheat crops in a semi-arid region. Agric. Water Manag. 2007, 87, 41–54. [Google Scholar] [CrossRef]
  27. Calera, A.; Gonzalez-Piqueras, J.; Meli, J. Monitoring barley and corn growth from remote sensing data at field scale. Int. J. Remote Sens. 2004, 25, 97–109. [Google Scholar] [CrossRef]
  28. Gonzalez-Dugo, M.P.; Neal, C.M.U.; Mateos, L.; Kustas, W.P.; Prueger, J.H.; Anderson, M.C.; Li, F. A comparison of operational remote sensing-based models for estimating crop evapotranspiration. Agric. For. Meteorol. 2009, 149, 1843–1853. [Google Scholar] [CrossRef]
  29. Heilman, J.L.; Heilman, W.E.; Moore, D.G. Evaluating the crop coefficient using spectral reactance. Agron. J. 1982, 74, 967–971. [Google Scholar] [CrossRef]
  30. Hunsaker, D.J.; Pinter, P.J., Jr.; Barnes, E.M.; Kimball, B.A. Estimating cotton evapotranspiration crop coefficients with a multispectral vegetation index. Irrig. Sci. 2003, 22, 95–104. [Google Scholar] [CrossRef]
  31. Neale, C.M.U.; Bausch, W.; Heermann, D. Development of reflectance-based crop coefficients for corn. Trans. ASAE 1989, 32, 1891–1899. [Google Scholar] [CrossRef]
  32. Zhou, M.C.; Ishidaira, H.; Takeuchi, K. Comparative study of potential evapotranspiration and interception evaporation by land cover over Mekong basin. Hydrol. Process. 2008, 22, 1290–1309. [Google Scholar] [CrossRef]
  33. Rodriguez-Iturbe, I. Ecohydrology: A hydrologic perspective of climate soil-vegetation dynamics. Water Resour. Res. 2000, 36, 3–9. [Google Scholar] [CrossRef]
  34. Priestley, C.H.B.; Taylor, R.J. On the Assessment of Surface Heat Flux and Evaporation Using Large-Scale Parameters. Mon. Weather Rev. 1972, 100, 81–92. [Google Scholar] [CrossRef]
  35. Thom, A.S. Momentum, Mass and Heat Exchange of Plant Communities. In Vegetation and Atmosphere; Monteith, J.L., Ed.; Academic Press: London, UK, 1975; pp. 57–110. [Google Scholar]
  36. Jarvis, P.G. The interpretation of the variations in leaf water potential and stomatal conductance found in canopies in the field. Philos. Trans. R. Soc. B 1976, 273, 593–610. [Google Scholar] [CrossRef]
  37. Barr, A.; Morgenstern, K.; Black, T.; McCaughey, J.; Nesic, Z. Surface energy balance closure by the eddy-covariance method above three boreal forest stands and implications for the measurement of CO2 flux. Agric. For. Meteorol. 2006, 140, 322–337. [Google Scholar] [CrossRef]
  38. Garratt, J. The Atmospheric Boundary Layer; Cambridge University Press: Cambridge, UK, 1993; p. 316. [Google Scholar]
  39. Kaimal, J.C.; Finnigan, J.J. Atmospheric Boundary Layer Flows-Their Structure and Measurement; Oxford University Press: New York, NY, USA, 1994; p. 289. [Google Scholar]
  40. Aubinet, M.; Grelle, A.; Ibrom, A.; Rannik, U.; Moncrieff, J.; Foken, T.; Kowalski, A.S.; Martin, P.H.; Berbigier, P.; Bernhofer, C.; et al. Estimates of the annual net carbon and water exchange of forests: The euroflux methodology. Adv. Ecol. Res. 2000, 30, 113–175. [Google Scholar]
  41. Mauder, M.; Foken, T. Documentation and Instruction Manual of the Eddy Covariance Software Package TK2; University Abt Mikrometeorol: Arbeitsergebn, Bayreuth, Germany, 2004; pp. 26–42. [Google Scholar]
  42. Migliavacca, M.; Galvagno, M.; Cremonese, E.; Rossini, M.; Meroni, M.; Sonnentag, O.; Cogliati, S.; Manca, G.; Diotri, F.; Busetto, L.; et al. Using digital repeat photography and eddy covariance data to model grassland phenology and photosynthetic CO2 uptake. Agric. For. Meteorol. 2011, 151, 1325–1337. [Google Scholar] [CrossRef]
  43. Galvagno, M.; Wohlfahrt, G.; Migliavacca, M.; Cremonese, E.; Rossini, M.; Filippa, G.; Manca, G.; Siniscalco, C.; Morra di Cella, U.; Colombo, R.; et al. Phenology and carbon dioxide source/sink strength of a subalpine grassland in response to an exceptionally short snow season. Environ. Res. Lett. 2013, 8, 025008. [Google Scholar] [CrossRef]
  44. Baldocchi, D.; Black, T.A.; Curtis, P.S.; Falge, E.; Fuentes, J.D.; Granier, A.; Gu, L.; Knohl, A.; Pilegaard, K.; Schmid, H.P.; et al. Predicting the onset of net carbon uptake by deciduous forests with soil temperature and climate data: A synthesis of FLUXNET data. Int. J. Biometeorol. 2005, 49, 377–387. [Google Scholar] [CrossRef] [PubMed]
  45. Hollinger, D.Y.; Ollinger, S.V.; Richardson, A.D.; Meyers, T.P.; Dail, D.B.; Martin, M.E.; Scott, N.A.; Arkebauer, T.J.; Baldocchi, D.D.; Clark, K.L. Albedo estimates for land surface models and support for a new paradigm based on foliage nitrogen concentration. Glob. Chang. Biol. 2010, 16, 696–710. [Google Scholar] [CrossRef]
  46. Stoy, P.; Katul, G.G.; Siqueira, M.B.S.; Juang, J.Y.; Novick, K.A.; McCarthy, H.R.; Oishi, A.C.; Uebelherr, J.M.; Kim, H.S.; Oren, R. Separating the effects of climate and vegetation on evapotranspiration along a successional chronosequence in the southeastern U.S. Glob. Chang. Biol. 2006, 12, 1–21. [Google Scholar] [CrossRef]
  47. Stoy, P.; Katul, G.G.; Siqueira, M.B.S.; Juang, J.Y.; Novick, K.A.; Uebelherr, J.M.; Oren, R. An evaluation of models for partitioning eddy covariance-measured net ecosystem exchange into photosynthesis and respiration. Agric. For. Meteorol. 2006, 141, 2–18. [Google Scholar] [CrossRef]
  48. Myneni, R.B.; Hoffman, S.; Knyazikhin, Y.; Privette, J.L.; Glassy, J.; Tian, Y.; Wang, Y.; Song, X.; Zhang, Y.; Smith, G.R.; et al. Global products of vegetation leaf area and fraction absorbed PAR from year one of MODIS data. Remote Sens. Environ. 2002, 83, 214–231. [Google Scholar] [CrossRef]
  49. Schaaf, C.; Gao, F.; Strahler, A.; Lucht, W.; Li, X.; Tsung, T.; Strugnell, N.; Zhang, X.; Jin, Y.; Muller, J.-P.; et al. First operational BRDF, albedo and nadir reflectance products from MODIS. Remote Sens. Environ. 2002, 83, 135–148. [Google Scholar] [CrossRef]
  50. Masseroni, D.; Corbari, C.; Mancini, M. Validation of theoretical footprint models using experimental measurements of turbulent fluxes over maize fields in Po Valley. Environ. Earth Sci. 2014, 72, 1213–1225. [Google Scholar] [CrossRef]
  51. Zhao, P.; Lüers, J. Parameterization of Evapotranspiration Estimation for Two Typical East Asian Crops. Atmosphere 2017, 8, 111. [Google Scholar] [CrossRef]
Figure 1. Methodology to compute crop coefficients with the four different techniques.
Figure 1. Methodology to compute crop coefficients with the four different techniques.
Sensors 17 02664 g001
Figure 2. Sensitivity coefficient of kc at LAI, albedo and vegetation height variations.
Figure 2. Sensitivity coefficient of kc at LAI, albedo and vegetation height variations.
Sensors 17 02664 g002
Figure 3. Crop coefficient values for Torgnon eddy covariance station over pasture for: (a) 2009; and (b) 2010.
Figure 3. Crop coefficient values for Torgnon eddy covariance station over pasture for: (a) 2009; and (b) 2010.
Sensors 17 02664 g003
Figure 4. Crop coefficient values for Chestnut Ridge eddy covariance station for 2007.
Figure 4. Crop coefficient values for Chestnut Ridge eddy covariance station for 2007.
Sensors 17 02664 g004
Figure 5. Crop coefficient values for Duke forest eddy covariance station for: (a) 2001; and (b) 2002.
Figure 5. Crop coefficient values for Duke forest eddy covariance station for: (a) 2001; and (b) 2002.
Sensors 17 02664 g005
Figure 6. Crop coefficient values for Black Hills forest eddy covariance station for: (a) 2006; and (b) 2007.
Figure 6. Crop coefficient values for Black Hills forest eddy covariance station for: (a) 2006; and (b) 2007.
Sensors 17 02664 g006
Table 1. Crop coefficient values for the four sites: kceddy, kcPenMon,sat, kcPenMon,ground, kcFAO.
Table 1. Crop coefficient values for the four sites: kceddy, kcPenMon,sat, kcPenMon,ground, kcFAO.
kc-inikc-midkc-end
Pasture
Torgnon kcFAO0.750.750.75
2009kceddy-0.80.47
kcPenMon,sat-0.760.49
2010kceddy-0.880.43
kcPenMon,sat-0.760.36
kcPenMon,ground-0.760.39
Deciduous Forest
Chestnut Ridge2007kcFAO0.150.910.15
kceddy0.190.470.2
kcPenMon,sat0.040.480.05
Duke Forest2001kcFAO0.150.80.15
kceddy0.110.510.12
kcPenMon,sat0.070.440.06
2002kcFAO0.150.90.15
kceddy0.090.43-
kcPenMon,sat0.040.440.07
Evergreen Forest
Black Hills2006kcFAO0.150.790.15
kceddy0.05 0.17 0.04
kcPenMon,sat0.04 0.20 0.07
2007kcFAO0.150.780.15
kceddy0.05 0.18 0.03
kcPenMon,sat0.05 0.20 0.04
Table 2. Lengths of crop development stages (Julian day of beginning and end) respect to kc-ini, kc-mid and kc-end.
Table 2. Lengths of crop development stages (Julian day of beginning and end) respect to kc-ini, kc-mid and kc-end.
Days kc-iniDays kc-midDays kc-end
Pasture
Torgnon2009 183–305
2010 143–303
Deciduous Forest
Chestnut Ridge20071–105139–269328–365
Duke Forest20011–71149–297320–365
20021–80121–281306–365
Evergreen Forest
Black Hills20061–98142–231281–365
20071–98189–247301–365

Share and Cite

MDPI and ACS Style

Corbari, C.; Ravazzani, G.; Galvagno, M.; Cremonese, E.; Mancini, M. Assessing Crop Coefficients for Natural Vegetated Areas Using Satellite Data and Eddy Covariance Stations. Sensors 2017, 17, 2664. https://doi.org/10.3390/s17112664

AMA Style

Corbari C, Ravazzani G, Galvagno M, Cremonese E, Mancini M. Assessing Crop Coefficients for Natural Vegetated Areas Using Satellite Data and Eddy Covariance Stations. Sensors. 2017; 17(11):2664. https://doi.org/10.3390/s17112664

Chicago/Turabian Style

Corbari, Chiara, Giovanni Ravazzani, Marta Galvagno, Edoardo Cremonese, and Marco Mancini. 2017. "Assessing Crop Coefficients for Natural Vegetated Areas Using Satellite Data and Eddy Covariance Stations" Sensors 17, no. 11: 2664. https://doi.org/10.3390/s17112664

APA Style

Corbari, C., Ravazzani, G., Galvagno, M., Cremonese, E., & Mancini, M. (2017). Assessing Crop Coefficients for Natural Vegetated Areas Using Satellite Data and Eddy Covariance Stations. Sensors, 17(11), 2664. https://doi.org/10.3390/s17112664

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop