Mapping Evapotranspiration Coefficients in a Temperate Maritime Climate Using the Metric Model and Landsat Tm

The applicability of a land surface temperature (LST)-evapotranspiration (ET) regression model to estimate ET fraction (ETrF) was tested in the temperate maritime climate of Central Ireland. In this study, the Mapping ET at high Resolution and with Internalized Calibration (METRIC) model was applied to calculate evapotranspiration from a mixed land cover area in Central Ireland. The ETrF values estimated on a pixel-by-pixel basis using two different surface roughness maps derived from two different estimates of Leaf Area Index (LAI) showed no significant differences. The METRIC based ETrF values from a Landsat Enhanced Thematic Mapper Plus (ETM+) image from 12 July 2013 were in similar range with those obtained from a regression model using land surface temperature (LST)-ET relationship derived from a Landsat 5 Thematic Mapper (TM) image of 17 July 2006. Results indicated that ETrF values could be estimated by using LST values from Landsat TM, even when stress conditions were prevailing in the region.


Introduction
Regional estimates of Actual Evapotranspiration (ET a ) are one of the main components of the hydrological balance and an important factor affecting energy, hydrological, carbon and nutrient cycles.ET a can be estimated from physical observations near the ground over homogeneous surfaces using energy balance techniques, such as eddy covariance (EC) and Bowen ratio-energy balance (BREB), or by a soil water balance technique such as surface renewal (SR) [1].These hydrological and micrometeorological approaches [2][3][4][5] prevailed until the mid 1990s [6] where thermal infrared image data from satellites became available for ET a mapping by remote sensing using Soil-Vegetation-Atmosphere Transfer (SVAT) models [7].ET a models using remote sensing data enable development of real-time management applications for activities such as irrigation management that depend on changes in water status and energy balance at the earth surface through time and by location.Medium-or high-resolution daily evapotranspiration maps are more suitable for regional management than "point" data calculated using balance techniques because they convey the pattern and distribution from local to regional scale.
Satellite data can be used to estimate ET a from available energy and sensible heat flux using surface radiation balance and the residual method [8,9].Mapping EvapoTranspiration at high Resolution and with Internalized Calibration (METRIC) is a variant of the Surface Energy Balance Algorithm for Land (SEBAL) model [10][11][12] that implements the residual method to calculate ET a maps.The METRIC model uses the ET fraction (ET r F), which serves the same mathematical purpose as crop coefficient (K c ), defined as the ratio of satellite image evapotranspiration to ET r representing the crop type and the development stage of the crop.The concept of reference evapotranspiration (ET r ) was introduced to study the evaporative demand of the atmosphere independently of crop type, crop development and management practices.ET r values measured or calculated at different locations or in different seasons are comparable as they refer to the evapotranspiration from the same reference surface [13].ET r is considered a climatic parameter and can be computed from meteorological data assuming a theoretical ET calculation from an ideal crop of alfalfa.The FAO Penman-Monteith method is recommended as the best method for determining ET r [13].K c (ET r F hence) is an adjustment to the ET r for non-alfalfa crops taking into account the difference in evapotranspiration between the crop and ET r .There may be several ET r F values for a single crop depending on the crop stage of development.Crop coefficients are generally preferable for the estimation of water requirements and water use efficiency and this is the main motivation to use ET r F instead of evapotranspiration [13].Growth stages, plant physiology and rooting depth are used to define change in crop coefficient through the growing season.
SEBAL and METRIC have been successfully used for many studies around the world [9][10][11]14].SEBAL was originally used in Egypt, Spain and Niger [15,16] (predominantly hot, semi-arid to desert climates) and later applied in Mediterranean countries including Turkey [15] and more recently Greece [16,17] and Cyprus [18] (predominantly variants of Mediterranean climate).It has been used in other dry regions of the world; in parts of South America, USA, Europe, Africa, and Asia [19][20][21][22][23][24][25][26][27]; and, more recently, in cold, high altitude areas [28].These climates typically have prolonged clear dry weather with a tendency to water stress at times of the year.The objective of this work was to evaluate the application of METRIC to a completely different environment in terms of climate (temperate maritime) and land uses (peatland, permanent grassland, coniferous forest) as found in central Ireland.

Methods
The study area is County Laois, located in the Midlands Region of Ireland (centred on 53 • 0 N, 7 • 24 W; Figure 1).The centre of the county is Slieve Bloom Mountains with a maximum elevation of around 530 m above sea level (ASL) and a plateau to the south-east of the county at 180 to 245 m ASL influence on the climate of the area, particularly temperature and rainfall distribution [29][30][31].According to Köppen and Geiger the climate can be characterized as temperate maritime (Cfb).December and January are usually the wettest months, while March, April and June are usually the driest.However, even the driest month still has abundant rainfall compared to Mediterranean and arid climates.At the Durrow weather station (52 • 50 43 N, 7 • 23 48 W), there are around eleven raindays per month with 50-60 mm per month in the summer, while the average annual temperature is 9.5 • C with an average maximum of 15.1 • C in July and August and average minimum of 4.4 • C in January.The coldest months are January and February, while July and August are the warmest [30].The prevailing wind direction is from the southeast to the west.Agriculture is the predominant land-use in the county and that was one of the main reasons for its selection.
The predominant land use is managed grassland, accounting for about 70% of the land area with the remaining land being urban, tillage farming, plantation forestry, uplands including scrub and peatland and low level peatlands [31].
The METRIC model was used to produce ET a and ET r F maps for 17 July 2006, utilizing a clear, high quality Landsat 5 TM image with a 30 m spatial resolution (path/row 206/23).Image data were acquiring from the United States Geological Survey (USGS) GloVis tool.These data were chosen because the image was virtually cloud free (rare in Ireland) and represented a time when a water stress had built up due to dry weather in the previous weeks.It should be noted that water stress in the context of a temperate maritime climate is very different from that experienced in semi-arid, Mediterranean and desert climates.Meteorological data from adjacent synoptic weather stations [30] were also used.Radiometric and geometric correction were applied to the initial TM image deploying the related header file given by USGS [32][33][34][35].Spectral radiance and reflectance image files for each TM band were derived using Digital Numbers (DNs) and Gain and Bias levels for each channel as described at Landsat 7 Science Data Users Handbook [36].The geometrically corrected images were projected to the Irish map conversion system using the Irish 1965 datum, based on a Transverse Mercator projection and the Airy Modified 1849 Spheroid.The same system is used throughout this study.Each image was finally cropped to the area of interest using the appropriate county shape file.The results of these processes were corrected radiance values for each TM band for the study region.The METRIC model eliminates the need for refined atmospheric correction of surface temperature and albedo due to its internal calibration [37,38], but these steps do improve the overall result [15, [39][40][41].
METRIC is centred on the equation of surface energy balance (Equation (1)): where λET is the latent heat flux (W/m 2 ), Rn is net radiation (W/m 2 ), G is soil heat flux (W/m 2 ) and H is sensible heat flux (W/m 2 ).Net radiation was computed from the land surface radiation balance (Equation ( 2)): where α is surface albedo, RS ↓ (W/m 2 ) is the incoming solar radiation, RL ↓ (W/m 2 ) is incoming long wave radiation, RL ↑ (W/m 2 ) is the outgoing long wave radiation and εo is the broad band surface emissivity.Soil heat flux G (W/m 2 ) was empirically estimated using the function of Bastiaanssen [15], which is based on albedo, surface temperature and NDVI (Equation (3)): where LST is surface temperature (K) and NDVI is the widely used Normalized Difference Vegetation Index.The sensible heat flux, H (W/m 2 ) for neutral atmospheric conditions was calculated by (Equation (4)): where ρ is air density (kg/m 3 ) related to atmospheric pressure, cp is air specific heat capacity (1004 J•kg −1 •K −1 ), dT (K) is the temperature difference (T1 − T2) between two reference heights (z1 and z2), and rah is aerodynamic resistance to heat transport (s•m −1 ).Equation (4) has two unknowns, rah and dT, which is normally difficult to solve.The first unknown parameter (rah) was computed for neutral stability as (Equation ( 5)): Radiometric and geometric correction were applied to the initial TM image deploying the related header file given by USGS [32][33][34][35].Spectral radiance and reflectance image files for each TM band were derived using Digital Numbers (DNs) and Gain and Bias levels for each channel as described at Landsat 7 Science Data Users Handbook [36].The geometrically corrected images were projected to the Irish map conversion system using the Irish 1965 datum, based on a Transverse Mercator projection and the Airy Modified 1849 Spheroid.The same system is used throughout this study.Each image was finally cropped to the area of interest using the appropriate county shape file.The results of these processes were corrected radiance values for each TM band for the study region.The METRIC model eliminates the need for refined atmospheric correction of surface temperature and albedo due to its internal calibration [37,38], but these steps do improve the overall result [15, [39][40][41].
METRIC is centred on the equation of surface energy balance (Equation ( 1)): where λET is the latent heat flux (W/m 2 ), R n is net radiation (W/m 2 ), G is soil heat flux (W/m 2 ) and H is sensible heat flux (W/m 2 ).Net radiation was computed from the land surface radiation balance (Equation ( 2)): where α is surface albedo, R S↓ (W/m 2 ) is the incoming solar radiation, R L↓ (W/m 2 ) is incoming long wave radiation, R L↑ (W/m 2 ) is the outgoing long wave radiation and ε o is the broad band surface emissivity.Soil heat flux G (W/m 2 ) was empirically estimated using the function of Bastiaanssen [15], which is based on albedo, surface temperature and NDVI (Equation ( 3)): where LST is surface temperature ( • C) and NDVI is the widely used Normalized Difference Vegetation Index.The sensible heat flux, H (W/m 2 ) for neutral atmospheric conditions was calculated by (Equation ( 4)): where ρ is air density (kg/m 3 ) related to atmospheric pressure, c p is air specific heat capacity (1004 J•kg −1 •K −1 ), dT (K) is the temperature difference (T 1 − T 2 ) between two reference heights (z 1 and z 2 ), and r ah is aerodynamic resistance to heat transport (s•m −1 ).Equation ( 4) has two unknowns, r ah and dT, which is normally difficult to solve.The first unknown parameter (r ah ) was computed for neutral stability as (Equation ( 5)): where z 1 and z 2 are heights in meters above the zero plane displacement, u * is the friction velocity (m/s), and k is von Karman's constant (0.41).Friction velocity (u * ) quantifies the turbulent velocity fluctuations in the air using the logarithmic wind law (Equation ( 6)) for neutral atmospheric conditions: where k is von Karman's constant, u x is the wind speed (m/s) at height z x , and z om is the momentum roughness length (m), which is a measure of interaction between the surface and the adjacent layer of air above.Momentum transfer is mainly caused by turbulent drag on roughness elements [42,43].
Wind speed at 200 m above the weather station (u 200 ) was calculated then using Equation ( 6).The height of 200 m is referred in METRIC as the "blending height" [12] (Equation ( 7)): where u * is the friction velocity at the weather station.This calculation was done on a spreadsheet [40].
Wind speed (u x ) was acquired from the weather station of the study, while momentum roughness length (z om ) at the station was empirically estimated using the following formula (Equation ( 8)) [42]: where h is the vegetation height (m).
Friction velocity u * for each pixel was finally computed using Equation ( 9): where z om is the particular momentum roughness length for each pixel, while u 200 is assumed to be constant for all the pixels of the image since it is unaffected by surface features.The momentum roughness length (z om ) for each pixel was computed using a land-use map and Leaf Area Index (LAI) [44].For agricultural areas, z om can be calculated as a function of LAI (Equation ( 10)): Non-agricultural areas, have known z om values [40] and can be assigned using a land-use map (Figure 1b) which has been already created for this study utilizing Corine Land-Cover raster data [45].
METRIC suggests linear change in dT with LST [10] meaning that Equation (4) has only two unknown parameters (a and b).This can be expressed as (Equation ( 11)): where a and b are correlation coefficients, and dT is defined as the difference between the air temperature very near the surface (at 0.1 m above the zero plane displacement height) and the air temperature at 2.0 m above the zero plane displacement height [40].To estimate a and b, the so-called "cold" and "hot" pixels are selected from the TM image.H values can be reliably estimated for these pixels which serve as boundary conditions for the total methodology.A "cold" pixel in METRIC is generally defined as one representing a wet, well-irrigated crop surface with full vegetation cover [10].A pixel with high value of LAI and low LST is a very good candidate to be a "cold" pixel.
A "hot" pixel is representing a dry, exposed soil, typical of a harvested agricultural field, where all the energy is used to heat the surface (ET = 0).A pixel with low values of LAI and high LST is a very good candidate to be a "hot" pixel.The selection of "hot" pixels is more difficult than the selection of "cold" pixel due to the wide range of temperatures for hot pixel candidates.For the assessment of "cold" and "hot" pixels, a pseudo-colour image was created for each of LAI, and LST maps in order to acquire a better visualization of the prevailing differences.Finally, cold and hot pixel candidates are recorded and then carefully selected for each case according to the above criteria.H was finally computed fitting a line according to Equation ( 11) and solving Equation ( 4) assuming neutral atmospheric conditions, meaning that no atmospheric stability or instability was existed.An iterative process was finally used to correct for possible atmospheric instability caused by buoyancy effects of surface heating by applying Monin-Obukhov theory [10,46].New dT values are then computed for the "cold" and "hot" pixels according to the estimated values of a and b, and a new corrected value for H is repeatedly computed while the stability correction is applied until H stabilizes.More details about atmospheric stability correction equations can be found in Bastiaanssen [10,11].Finally, reference ET fraction values (ET r F) can be calculated using the standardized ASCE Penman-Monteith equation applied using local meteorological observations assuming ET r F is constant over the day (Equation ( 12)): where ET ainst is the instantaneous ET a (mm/day) and ET r (mm/day) is the evapotranspiration of a known reference crop, usually alfalfa for METRIC.A new image file where each pixel denotes a unique ET r F value, with 30 m × 30 m spatial resolution, is the final product of the model.The whole METRIC model is summarized at Figure 2.
Water 2017, 9, 23 5 of 14 pixel candidates are recorded and then carefully selected for each case according to the above criteria.H was finally computed fitting a line according to Equation (11) and solving Equation ( 4) assuming neutral atmospheric conditions, meaning that no atmospheric stability or instability was existed.An iterative process was finally used to correct for possible atmospheric instability caused by buoyancy effects of surface heating by applying Monin-Obukhov theory [10,46].New dT values are then computed for the "cold" and "hot" pixels according to the estimated values of a and b, and a new corrected value for H is repeatedly computed while the stability correction is applied until H stabilizes.More details about atmospheric stability correction equations can be found in Bastiaanssen [10,11].Finally, reference ET fraction values (ETrF) can be calculated using the standardized ASCE Penman-Monteith equation applied using local meteorological observations assuming ETrF is constant over the day (Equation ( 12)): where ETainst is the instantaneous ETa (mm/day) and ETr (mm/day) is the evapotranspiration of a known reference crop, usually alfalfa for METRIC.A new image file where each pixel denotes a unique ETrF value, with 30 m × 30 m spatial resolution, is the final product of the model.The whole METRIC model is summarized at Figure 2. As it can be seen from Figure 2, METRIC model has an obvious complexity, consisting of many sub-models and routines, but, as a residual method, is based on Equation (1) which has a very simple meaning: evapotranspiration can be computed based on Rn, G and H. Sensible Heat Flux (H) is a very important parameter and is generally difficult to compute, but this is where METRIC serves as a solution, utilizing the theory of "boundary" pixels.However, one of the most important steps for the As it can be seen from Figure 2, METRIC model has an obvious complexity, consisting of many sub-models and routines, but, as a residual method, is based on Equation (1) which has a very simple meaning: evapotranspiration can be computed based on R n , G and H. Sensible Heat Flux (H) is a very important parameter and is generally difficult to compute, but this is where METRIC serves as a solution, utilizing the theory of "boundary" pixels.However, one of the most important steps for the computation of H is the computation of roughness length, playing a very critical role in the computation of friction velocity resistance (Equations ( 4)-( 6)), and this is where this study mainly contributes [47,48].The original METRIC model suggests computing roughness length using leaf area index (LAI, m 2 •m −2 ) multiplied by 0.018, where LAI is computed as (Equation ( 13)): SAVI is computed as (Equation ( 14)): where Ch x refers to the reflectance expressed as digital number of the Landsat TM channel number (x) and L is an empirical parameter that is difficult to compute for Ireland due to the lack of available images [40].These relationships (Equations ( 13) and ( 14)) were developed for ID, USA, which is very different from the study region so an alternative method for computing LAI was used derived from Enhanced Vegetation Index 2 (EVI2) [49][50][51][52][53], using the formula [50] (Equation ( 15)): and LAI was finally derived from EVI2 using the method of O'Connell et al. [52,53] (Equation ( 16)): The described methodology has the ability to compute LAI using EVI2 avoiding the use of SAVI at the study region which is difficult to compute [40,52,54].Eventually, LAI affects the computation of roughness length, which in turn is affecting the computation of friction velocity resistance (Equations ( 4)-( 6)), and finally the computation of H (Equation ( 4)).The new methodology of LAI computation was incorporated into the total METRIC model, and the new results were compared with the METRIC ET r F.
Furthermore, LST values generated using Landsat's Thematic Mapper (TM) thermal band [40,55] were correlated with ET r F values, considering that ET r F values are well correlated with thermal units [56][57][58], while Simplified Surface Energy Balance (SSEB) assumes that actual ET varies linearly with LST [59].LST is a governing parameter in water and energy exchange between land and atmosphere and is also an essential input to most of the irrigation water management activities [60].Additionally, Sun et al. [61] found linear relationships between ET and LST values, while there is a strong indication that maybe there is a significant correlation even if stress conditions are prevailing over regions having maritime climatic conditions [62,63].A linear equation was finally derived between ET r F and LST values for 17 July 2006 at this study.
For the validation of the derived equation a new image from the same region under the same climatic conditions had to be tested, and that was during the 12 July 2013.The same methodology was followed as that with the initial 2006 image, but with the exception that Landsat 5 was in no operation any more, thus, Landsat ETM+ imagery has been utilized instead.ETM+ has almost the same spectral characteristics as TM [36].A new METRIC ET r F map as well as a new LST map based on low gain band 6 from ETM+ were then generated for the 12 July 2013 using the described model.It is crucial to note that a new pair of cold and hot pixel has been used within the model of METRIC.The previous derived linear equation was then applied to the latest case and a new modelled ET r F map was generated.METRIC ET r F values and modelled ET r F values were compared, applying statistics between the two corresponding images.

Results and Discussion
ET r F image file product was firstly derived according to METRIC (Figure 3a).The computation of METRIC was implemented again using the new EVI2 based LAI values, and a new ET r F map was generated (Figure 3b).Examining Figure 3         Equation ( 17) was then tested by applying to the initial LST image file, producing a new modelled ET r F image as shown in Figure 6.The values of the new map have an average value of 0.64 (std = 0.18) while the initial METRIC map has an average of 0.72 (std = 0.16).The total correlation between the initial and final image map applied to a number of 1,909,074 pixels is high (R = 0.63).differencing defined as the subtraction of the pixel digital values of the first image recorded at one date from the corresponding pixel values of the second image was applied then [64][65][66][67].The results showed that the difference image file between METRIC ET r F map and the modelled one has an average of 0.02 (std = 0.06).
Water 2017, 9, 23 8 of 14 Equation ( 17) was then tested by applying to the initial LST image file, producing a new modelled ETrF image as shown in Figure 6.The values of the new map have an average value of 0.64 (std = 0.18) while the initial METRIC map has an average of 0.72 (std = 0.16).The total correlation between the initial and final image map applied to a number of 1,909,074 pixels is high (R = 0.63).Image differencing defined as the subtraction of the pixel digital values of the first image recorded at one date from the corresponding pixel values of the second image was applied then [64][65][66][67].The results showed that the difference image file between METRIC ETrF map and the modelled one has an average of 0.02 (std = 0.06).A validation of Equation ( 17) was finally implemented using the new case of 12 July 2013.Two ETrF maps were generated applying METRIC (average value 0.61) and the model (Equation ( 17)) respectively (average value 0.60), and the results are shown in Figure 7.The total correlation between METRIC and model based ETrF values is high (R = 0.88).A validation of Equation ( 17) was finally implemented using the new case of 12 July 2013.Two ET r F maps were generated applying METRIC (average value 0.61) and the model (Equation ( 17)) respectively (average value 0.60), and the results are shown in Figure 7.The total correlation between METRIC and model based ET r F values is high (R = 0.88).A validation of Equation ( 17) was finally implemented using the new case of 12 July 2013.Two ETrF maps were generated applying METRIC (average value 0.61) and the model (Equation ( 17  The mean value for the difference image obtaining after subtracting the two previous images (Figure 8a) is 0.005 (std = 0.07).It is clear that larger differences (in red colour) can be seen over regions defined as inland water at Figure 1.For that reason, separate regression models were developed for three land cover types.More specifically, grass, forest and agricultural were merged into one class (GFA), while artificial (AR) and inland wetlands (IW) used as they were.The result of these regressions is summarized by the following equations (R 2 values > 0.90): ET r F =19.161 − 0.0610 × LST (GFA) ET r F = 19.944− 0.0637 × LST (AR) ET r F = 20.288− 0.0642 × LST (IW) (20) Comparisons among model outputs using Equations ( 17)- (20) were then conducted on a pixel-by-pixel basis for the three classes.The analysis showed that the resulted ET r F values were not significant different between GFA and AR classes, but a remarkable variation appeared when using Equation (20).It was therefore decided to utilize Equation ( 17) throughout the LST image except the pixels belonging to IW class.A mask was created for those pixels, and a separate equation (Equation ( 20)) was used instead.The whole ET r F image map was a result of merging the two previous images resulting from Equations ( 17) and (20) into one.After applying the previous methodology to the new modelled ET r F image file, a new difference image file between METRIC and modelled ET r F was generated (Figure 8b).The results finally showed a clear improvement (difference mean value = 0.004, std = 0.06) as might have been expected.
Interestingly, the evaluation of the energy balance components, such as the net radiation and the sensible and latent heat fluxes, is characterized by significant uncertainties related to surfaces heterogeneity [68].Therefore, a new pixel oriented parameterization would be required for IW surfaces [69].This issue must be further considered as new sensors in higher spatial and temporal resolutions become available, contributing to the development of new land use classifications identical for energy balance applications.pixels belonging to IW class.A mask was created for those pixels, and a separate equation (Equation ( 20)) was used instead.The whole ETrF image map was a result of merging the two previous images resulting from Equations ( 17) and ( 20) into one.After applying the previous methodology to the new modelled ETrF image file, a new difference image file between METRIC and modelled ETrF was generated (Figure 8b).The results finally showed a clear improvement (difference mean value = 0.004, std = 0.06) as might have been expected.Interestingly, the evaluation of the energy balance components, such as the net radiation and the sensible and latent heat fluxes, is characterized by significant uncertainties related to surfaces heterogeneity [68].Therefore, a new pixel oriented parameterization would be required for IW surfaces [69].This issue must be further considered as new sensors in higher spatial and temporal resolutions become available, contributing to the development of new land use classifications identical for energy balance applications.
Summarizing, METRIC model showed reasonable results estimating the ETrF values on a pixel by pixel basis, but the main weakness of this study remains the absence of real measured evapotranspiration data for the appropriate validation.However, METRIC is very well established between other energy balance models for many years now [9].It also noteworthy that Earth Engine Evapotranspiration Flux (EEFlux) which is a web-based tool operating on the Google Earth Engine and computational cloud uses METRIC model as foundation [70].Summarizing, METRIC model showed reasonable results estimating the ET r F values on a pixel by pixel basis, but the main weakness of this study remains the absence of real measured evapotranspiration data for the appropriate validation.However, METRIC is very well established between other energy balance models for many years now [9].It also noteworthy that Earth Engine Evapotranspiration Flux (EEFlux) which is a web-based tool operating on the Google Earth Engine and computational cloud uses METRIC model as foundation [70].
Another limitation of the study is that the model has only been applied to a small scale.However, after establishing the appropriate regression equations of the model, the same equations might be applied to other sensors [71].A well-known procedure for taking advantage of a better temporal resolution is called downscaling [72,73].It has been suggested [73][74][75] that is feasible to apply various downscaling methods to combine MODIS and Landsat imagery in order to obtain both high temporal and high spatial ET resolutions.Using a daily temporal resolution like MODIS [75,76] it is easy to exploit new images and apply the model to a new data set, strengthening the applicability of the procedure.Furthermore, introducing hyperspectral data is possible to improve the accuracy of LST computation, while the use of microwave data during the days affected by weather could make an LST estimation feasible.An integrated use of multi-source, multi-scale remote sensing data to expand temporal scale of ET r F based on the characteristics of study area, is a breakthrough for future research [77][78][79].On such basis, this study provides the framework for future analysis elsewhere.

Conclusions
In this study, METRIC model was applied to compute ET r F values in a temperate maritime climate.In order to examine the proposed model, Laois Co. central Ireland has been selected mainly because it is a typical regional agricultural region.Landsat TM, ETM+ and meteorological data have been used together in order to compute energy balance parameters necessary for the final computation of ET r F values.The results show reasonable spatial variations in correspondence with the already known land-use classes of the region.Then, METRIC model was implemented using locally adapted LAI values, given as a function of EVI2.LAI is much easier to compute using the proposed formula, since it is not dependent on time-series of remotely sensed physical parameters which is difficult to obtain.There was a significant positive correlation.
Furthermore, an approach based on the correlation between LST and ET r F values has been introduced, which shows high performance in estimating ET r F values after applying the related regression model (R 2 = 0.89).The approach was based on a Landsat TM image taken on 17 July 2006.
A linear equation explained that relation was validated using a new Landsat ETM+ image taken on 12 July 2013.Subtracting the new modelled ET r F image from METRIC ET r F image, the findings were reasonable.However the performance was not ideal especially on inland water pixels, leading to a different approach, where those pixels were masked and excluded from the initial regression.Another regression has been performed only for those pixels belonging to IW class, and then the two results were merged.The final difference image map showed a clear improvement in comparison with the first approach.
The most important limitation of the study is the fact that there is no available in-situ validation data.Another limitation is that the model remains local in nature as it uses a locally derived LAI equation and a regression model derived from one image.However, when using downscaling techniques and different sensors in higher spatial and temporal resolutions, it is feasible to take advantage of the derived regression equations.Future work will concentrate on those downscaling techniques, combining with the use of new thermal in-situ sensors.Further studies, which take land-use classification into account, will also need to be undertaken, because this study showed a remarkable sensitivity in that aspect.
In conclusion, a regression between METRIC ET r F and LST values is suggested using Landsat imagery, applied in a temperate maritime climate.Although METRIC is a promising tool for mapping ET r F values accurately worldwide, additional evaluation is needed under a variety of crop, land-use and weather conditions to fully assess its capability to accurately estimate spatially distributed ET r F values.

Figure 8 .
Figure 8. Difference between METRIC and modelled ET r F (12 July 2013): (a) initial case; and (b) after masking IW class.