Methodologies and Uncertainties in the Use of the Terrestrial Chlorophyll Index for the Sentinel-3 Mission

A methodology is described for the validation of Medium Resolution Imaging Spectrometer (MERIS) Terrestrial Chlorophyll Index (MTCI) data over heterogeneous land surfaces in an agricultural region in Southern Italy. The approach involves the use inverse canopy reflectance modeling techniques to derive maps of canopy chlorophyll content (CCC) and leaf area index (LAI) at fine spatial resolution. Indirect field measurements are used for validation of the fine spatial resolution data. Subsequently, these maps are aggregated based on a regular grid at 1 km spatial resolution to validate MERIS Level 2 MTCI (300 m). RapidEye satellite sensor data with a pixel size of 6.5 m are used for this purpose. Based on a set of independent ground measurements, fine spatial resolution maps achieved an R2 = 0.78 and RMSE = 0.39 for CCC and R2 = 0.76 and RMSE = 0.64 for LAI. The relationship between MERIS L2 MTCI and CCC [g∙m−2] achieved a coefficient of determination of 0.74 and it resulted to be extremely statistically significant (p-value < 0.001). Additionally, a relative validation of two other satellite products at medium resolution spatial scale, namely MERIS leaf area index (LAI) and Moderate Resolution Imaging Spectrometer (MODIS) LAI was performed by comparison with the fine spatial resolution LAI map. Results indicated a better accuracy in LAI estimation of MERIS (RMSE = 0.33) compared to MODIS (RMSE = 0.81) data.


Introduction
Medium spatial resolution Earth observation (EO) data are used routinely to monitor land surface characteristics from regional to global scales.Observations by NASA's Moderate Resolution Imaging Spectrometer (MODIS), CNES's Vegetation and ESA's Medium Resolution Imaging Spectrometer (MERIS) have been used to obtain systematic estimates of terrestrial biophysical variables such as Fraction of Absorbed Photosynthetically Active Radiation (fAPAR); leaf area index (LAI) [1,2]; canopy chlorophyll content (CCC) [3,4] and vegetation phenology [5].Estimates of these variables play an important role in ecosystem modeling and broader environmental studies and contribute to our understanding of biogeochemical fluxes and global climate.For instance, fAPAR is a measure of the photosynthetic activity and so an indicator of the amount and productivity of vegetation.For example, it is used in terrestrial ecosystem models to estimate the assimilation of carbon dioxide by vegetation.LAI characterizes the leaf surface available for energy and mass exchange between surface and atmosphere.It indicates the amount of leaf material in an ecosystem, which imposes an important control on photosynthesis, respiration, rain interception and other processes that link vegetation to climate [6].The canopy chlorophyll content is a measure of potential productivity and is a strong indicator of vegetation physiological status and health.For example, it has been used successfully for the detection of vegetation stress, photosynthetic capacity and productivity [7,8].
In the framework of the Global Monitoring for Environment and Security (GMES) initiative, algorithms to estimate land biophysical variables will be implemented in order to provide a wide range of operational services in applications such as understanding carbon and water cycle, estimating agricultural and forest productivity and monitoring environmental change.The GMES initiative will also guarantee a long-term sustainable observation ability capable of supporting the development of a family of satellites for regional, European and global scale monitoring.Within this context the European Space Agency (ESA) initiated the ocean and global land monitoring Sentinel-3 mission.This aims to provide continuity to operational services that have been developed via the European Remote Sensing (ERS) satellite and Environmental Satellite (Envisat) missions [9].Among other sensors, the Sentinel-3 payload will include the Ocean and Land Colour Instrument (OLCI).The OLCI will be an imaging spectrometer with a spectral sampling similar to that used for sea and land monitoring in recent years by Envisat-MERIS and SPOT-Vegetation.With a spatial resolution of 300 m (for land and coastal monitoring with 1.2 km in open ocean) and an ability to sample 21 narrow spectral intervals from ultraviolet to near-infrared, OLCI is designed to provide moderate spatial resolution biophysical products operationally [9].
In the terrestrial segment of the mission, two biophysical variables will be calculated and produced as Level 2B main products: OLCI Global Vegetation Index (OGVI) and OLCI Terrestrial Chlorophyll Index (OTCI).The OGVI is related to fAPAR whereas the OTCI is related to canopy chlorophyll content (CCC).Both biophysical variables are based on the MERIS heritage of operational land products [1,3] and will provide data continuity from Envisat's launch in 2002.The OTCI will be based on the MERIS Terrestrial Chlorophyll Index (MTCI) [3], which is already routinely distributed to users at regional and global scales as standard MERIS Level 2 (L2) products or eight days composites.The MTCI exploits the shift occurring in the red edge position of the spectral reflectance of vegetation.This shift in red edge position is related to changes in CCC [10][11][12], and has the advantage of remaining responsive even at very high CCC where traditional vegetation indices saturate.This CCC index has been applied successfully to many applications [13][14][15][16].The algorithm is fast, can be implemented easily with large data sets and is therefore suitable for operational purposes.
The assessment of accuracy of the data products is a pre-requisite for the provision of reliable and time-consistent datasets to monitor annual and seasonal vegetation trends.This is particularly important in the context of the GMES initiative, since the products will be used to provide value added services or will be distributed directly to users such as environmental agencies, decision and policy makers.
Sentinel-3 is scheduled for launch in 2013.However, the current available MERIS database offers an opportunity to investigate the potential of the OLCI Sentinel-3 mission and evaluate the performance of existing and new methodologies and instruments for land surface characterization.
The aim of this paper is to assess the performance of medium spatial resolution MERIS L2 MTCI products to estimate CCC.For ease of discussion we will call this 'validation'.The study also evaluates approaches to scaling biophysical variables from fine to medium spatial resolution and outlines the uncertainty associated with validation.

Validation of EO-Based Medium Spatial Resolution CCC Products
The validation of CCC products derived from remotely sensed data is particularly challenging since the signal observed by satellite sensors integrates the effects of leaf biochemical components and canopy structure (i.e., foliage amount, spatial arrangement and leaf orientation).Moreover, it is influenced by the intercepted radiation from other elements such as non-photosynthetic materials (branches, stems and shoots) and underlying soils [17,18].Determination of CCC in the field requires the estimation of two variables: Leaf area index and leaf chlorophyll concentration.The product of these two variables provides an approximation of chlorophyll content at canopy level that can be compared with that observed from satellite sensors at medium spatial resolutions [19].Several field methods have been proposed for the estimation of LAI [20] and leaf chlorophyll concentration [21,22].Whatever procedure is adopted, each field estimate of LAI or leaf chlorophyll concentration should refer to a common elementary sampling unit (ESU) that is representative, wherever possible of one vegetation type.The spatial sampling strategy (e.g., transects, grid sampling, random, etc.) and intensity (number of measurements per ESU) can vary depending on the equipment, canopy architecture and site homogeneity.Moreover, the extent of each ESU should correspond to the spatial resolution of remotely sensed products under validation.This is particularly complex to achieve at the scale of medium and coarse resolution satellite sensor data due to land cover heterogeneity and poor spatial correspondence between ESU and image pixel as a result of geolocation errors [23][24][25].Direct comparisons between estimates of a variable in the field and products derived from a satellite sensor data are feasible for large homogenous areas [26][27][28].However, where the land surface landscape is heterogeneous, alternative validation approaches are required.Several regional and global validation studies [26,29] have used fine spatial resolution remotely sensed data (10-30 m) to map biophysical variables as an intermediate step in scale between field biophysical estimates and medium or coarser spatial resolution satellite sensor imagery.For instance, validation of MERIS biophysical variables by the VALERI group (http://www.avignon.inra.fr/valeri/)involved first, establishing statistical 'transfer functions' between fine spatial resolution (e.g., 20 m) satellite sensor images and field biophysical estimates.Then, fine spatial resolution biophysical variable maps are produced for comparison with medium spatial resolution MERIS products.This intermediate step (between field biophysical estimates to fine spatial resolution satellite sensor data) provides a more appropriate representation of spatial variability in the land surface and reduces the mismatch between field biophysical estimates (often acquired at a spatial resolution of few meters) and satellite sensor data acquired at spatial resolutions ranging from 300 m to 1 km or more.This validation approach was presented in the framework of the Land Product Validation Subgroup of the Committee on Earth Observation Satellites (CEOS) [29] as the 'bottom-up' approach to validation of LAI products at a global scale from field biophysical measurements.
Indirect (and typically non-destructive) measurement techniques have been used to estimate LAI and leaf chlorophyll concentration in the field because they are simple and time efficient when compared to direct (and typically destructive) measurement techniques.Indirect measurement techniques when used adequately can provide estimates that are close to direct measurements for a range of vegetation types and conditions and can be used speedily over large areas [30].For instance, LAI estimates can be derived using the Li-Cor LAI-2000 sensor (Li-Cor, Inc.USA); digital hemispherical photography (DHP) [20]; quantum sensors, such as TRAC (3 rd Wave Engineering, ON, Canada) and Decagon AccuPAR (Decagon Devices, Inc.) or ceptometers (Delta-T devices Ltd, UK).However, hemispherical sensors (such as LAI-2000 and DHP) have been shown to be more robust than other sensors [31,32] in terms of their limited sensitivity to environmental conditions, operational procedures and assumptions, for example in relation to canopy architecture.The LAI-2000 instrument uses a fish-eye lens to measure the canopy light interception in multiple directions at five angles.LAI and the foliage mean tilt angle (MTA) are then estimated directly by inverting a simple radiative transfer model for vegetative canopies [33].The digital hemispherical photography technique is based on upward-looking fish-eye photographs taken from beneath the plant canopy that are processed using dedicated software to estimate LAI [34,35].A detailed review of field techniques for the estimation of LAI can be found in [20,[30][31][32].
Non-destructive estimates of leaf chlorophyll concentration can be obtained from portable leaf chlorophyll meters such as a Minolta SPAD-502 (Minolta Camera Co. Ltd., Japan) or an Opti-Sciences CCM-200 meter (Opti-Sciences, Inc., USA).These instruments use transmitted radiation in selected visible wavelengths to estimate chlorophyll concentration in leaves.Although such contact devices are quicker than destructive sampling, the measurements of transmission represent relative values of chlorophyll concentration and are very sensitive to leaf structural properties such as thickness, density, cuticle thickness and pubescence [36][37][38].However quantitative leaf chlorophyll concentration estimates can be obtained using a conversion equation derived from destructive extraction procedures [21] and contemporaneous contact non-destructive measurements.While a useful approach, [39][40][41] have reported curvilinear relationship with sensor sensitivity declining at higher levels of leaf chlorophyll concentration [42][43][44][45][46].An overview of the scatter around relationships between SPAD-502 readings and leaf chlorophyll concentration is presented in Figure 1.Based on this dataset acquired in our past experimental campaigns linear relationships are observed for the majority of leaf types.

Figure 1.
An overview of the scatter around relationships between SPAD-502 readings and leaf chlorophyll concentration extracted in the laboratory for nine leaf types.This dataset is part of our past experimental campaigns and destructive chlorophyll extraction was performed according to the procedures described in [47,48].

Overview
In this study, high spatial resolution maps of LAI and CCC were derived using inverse canopy reflectance modeling with RapidEye imagery.The resulting maps were aggregated at 1 km spatial resolution to validate MERIS L2 MTCI.Field measurement techniques were used during an experimental campaign to estimate LAI and leaf chlorophyll concentration and to validate fine spatial resolution LAI and CCC maps.Generalized and, where possible, canopy specific conversion equations were used to estimate leaf chlorophyll concentration from SPAD-502 measurements.A schematization of the validation processing is presented in Figure 2. Particular attention was dedicated to LAI; a relative validation of two other satellite products at medium resolution spatial scale, namely MERIS LAI and MODIS LAI, was performed by comparison with the fine spatial resolution LAI map.

Figure 2. Schematization of the validation process for the Medium Resolution Imaging
Spectrometer (MERIS) L2 Terrestrial Chlorophyll Index (MTCI) data.Similarly, this validation scheme was used for MERIS and Moderate Resolution Imaging Spectrometer (MODIS) leaf area index (LAI) data.

Field Campaign
This study was undertaken at the 560 km 2 'Piana del Sele' site in the Campania region of Southern Italy (Figure 3) (Lat.40.52°N,Long.15.00°E).The study site, one of the largest agricultural areas in the region, is characterized by irrigated agriculture (mainly forage crops, fruit trees and vegetables) with an average field size of about 2 ha [47].Within this large and relatively flat region, LAI and leaf chlorophyll concentration were estimated within 400 m 2 elementary sampling units (ESU) that were relatively homogeneous in terms of both vegetation type and growth stage.The center of each ESU was geo-located by means of a GPS, with an accuracy of ±3-5 m.
LAI and leaf chlorophyll concentration were estimated using a random sampling scheme within each of 36 ESUs.These ESUs were located in the study site according to field accessibility but were distributed within the area of interest to cover a variety of crop conditions.A LAI map derived from a SPOT fine spatial resolution image acquired on 12 August 2009 was used to locate possible fields for LAI and CCC sampling.The SPOT-LAI map was derived using a semi-empirical relationship based on the Weighted Difference Vegetation Index (WDVI).A logarithmic relationship was used to estimate LAI from WDVI, which was calibrated during several field campaigns in the study site in the  LAI was estimated using the Plant Canopy Analyser LAI-2000.The instrument has four well-documented limitations and two of these were dealt with explicitly in this study.Measurements were made during early morning and late afternoon under diffuse light conditions to minimize the effects of direct sunlight and thereby minimize LAI underestimation [49].A view cap of 180 degrees was used to physically limit the sensor field-of-view and thus to reduce interference due to the presence of an operator.However, it was accepted that the sensor did not distinguish photosynthetically active leaf tissue from other plant elements, such as stems, flowers or senescent leaves and so LAI would be overestimated, especially in woody species.By contrast, the 'clumping effect', i.e., non-random positioning of canopy elements, was also neglected causing LAI underestimation.Different definitions of LAI have been used in the literature depending on vegetation type and measurement device [31].These have included green LAI [50], effective LAI [51] and plant area index [52].The 'LAI' recorded using the LAI-2000 sensor represents a measure close to the effective Plant Area Index ('PAIe') [53,54] for reasons discussed above.The average LAI, resulting from two replications of one above-canopy and nine below-canopy measurements, was taken as representative measure for each ESU.The two replications (for a total of 18 measurements) were sampled randomly within the ESU.LAI values within the 36 ESUs ranged between 0.5 and 5.5 with a mean of 2.2 and a standard deviation of 1.3.
Leaf chlorophyll concentration was estimated using the SPAD-502 chlorophyll meter.At each ESU, 30 random measurements of leaves in different canopy layers were averaged to produce one value.Destructive sampling was not possible during the campaign.Therefore crop-specific conversion equations were taken from the literature (peach, apricot, plum trees: [55]; maize: [56]; and a generalized conversion equation for the other vegetation types: [57]) and used to convert SPAD-502 values to leaf chlorophyll concentration [g•m −2 ] (see Table 1).
SPAD-502 values within the 36 ESUs ranged between 30 and 68 with a mean of 48 and a standard deviation of 10.The total CCC for an indicative m 2 of ground was obtained by multiplying leaf chlorophyll content with the corresponding LAI of each ESU.The mean and standard deviation for LAI and CCC field estimates are provided in Table 2. Though care was taken to distribute the ESUs within the study site, their locations were restricted to fields with good access, which could have introduced a bias in the field data collection.

Fine Spatial Resolution Data from RapidEye Constellation
Multispectral data were acquired on 17 August 2009 (at 10:35 UTC) from Choros, a satellite in the RapidEye constellation.The sensor records radiance in five broad spectral bands corresponding to blue, green, red, red-edge and near-infrared part of the electromagnetic spectrum at a spatial resolution of 6.5 m.Four scenes, acquired within a few seconds, with a maximum across-track incident angle of 5° were adequate to cover the study site (about 560 km 2 ).Radiometrically calibrated Level 3A data were provided with a geometric accuracy of 13.95 m (root-mean-square-error, RMSE = 6.50 m).Further geometric correction was performed using Ground Control Points (GCPs) and the scenes were resampled using the nearest neighbor method, resulting in a final geo-location accuracy of about 3 m, which is consistent with the area of the ESU (20 × 20 m) and allows a small cluster of pixels to be averaged and compared with field estimates.One scene, containing identifiable ground targets (i.e., asphalt, sea water, concrete and sand), was atmospherically corrected by using ATCOR-2 [58].The known spectral reflectance of these targets was used for the parameterization of the atmospheric radiative transfer model.Subsequently, an empirical line method was applied to normalize the other three scenes to the first image.For this purpose, uniform areas in the overlapping regions between adjacent images were considered: Twenty zones of about 200 m 2 representing dark and bright surfaces were selected for each image and correction functions were derived for each spectral band.

Medium Spatial Resolution Data from MERIS
The MERIS L2 image used in this study was acquired on 17 August 2009 (at 10:33 UTC) at a spatial resolution of approximately 300 m (full resolution) and made available through the ESA EO-LISA archive.The MERIS L2 data are accompanied by different types of geophysical information (according to the level of geophysical processing for 'Land', 'Water' and 'Cloud'), Rayleigh-corrected reflectances and bottom-of-atmosphere vegetation indices, such as MTCI.The image geometric correction was performed using the map projection tool (nearest neighbor resampling method) implemented in the Basic ENVISAT Toolbox for (A)ATSR and MERIS (BEAM).
The MERIS top-of-canopy (TOC) vegetation algorithm (implemented as a plug-in 'TOC-VEG' in the ESA BEAM VISAT software) was used to estimate LAI from MERIS L2 reflectance data.As detailed in the algorithm description, TOC-VEG [27,59] is a model-inversion technique applied to a coupled leaf (PROSPECT, [60]) and canopy (SAIL, [61]) model.By means of a neural network four vegetation biophysical variables are derived: fAPAR, fCover, LAI and LAI×chlorophyll content.

Medium Spatial Resolution Data from MODIS
The MODIS LAI product (tile h19v04) was collection 5 (MOD15A2), 8-day composite (from 13 to 20 August 2009), atmospherically-corrected, and at a 1 km spatial resolution.Data were made available through the Earth Resources Observation System (EROS) Data Centre Distributed Active Archive Center.The MODIS LAI product was derived using a Look-Up-Table (LUT), called the 'main algorithm', which provides an inversion of a three-dimensional radiative transfer (RT) model [2] within six canopy biomes: Grasses and cereal crops, shrubs, broadleaf crops, savannah, broadleaf and needle-leaf forests [62].However, the main LAI algorithm is not applied when there are large uncertainties in surface reflectance data, misclassifications in the biome map or limitations in the LUTs [63].Should this happen, a backup algorithm based on empirical relations between the Normalized Difference Vegetation Index (NDVI) and LAI associated with a biome classification map is triggered in order to estimate LAI.A quality flag associated to the product indicates the algorithms used for producing LAI maps.The quality flag (SFC_QC) has five possible values; '0' and '1' refer to estimates by the main algorithm (LUT and RT) with '0' being the 'best' possible result.'2' and '3' refer to estimates by the backup algorithm (NDVI) and '4' for those that could not be estimated and so no data are given.The LAI values with the SFC_QC = 1 are referred to saturation LAI data and they usually do not represent the LAI on the surface even though they do provide a lower limit of the 'true' LAI.As would be anticipated, saturation usually occurred in areas with dense canopies [64].

Deriving LAI and CCC Maps from Fine Spatial Resolution Data from RapidEye Imagery
Two groups of techniques have been used to scale-up estimates of biophysical variables at the field scale in order to derive maps of those variables at fine spatial resolution [29].The first are based on empirical equations, such as regressions between vegetation indices and field data.The second are based on physical models (i.e., inverse canopy reflectance modeling, [26]).Empirical methods are site-and sensor-specific that require accurate field data [65] and most employ vegetation indices using data in only two wavebands, usually red and near infrared [23,66].Such simple ratios are rarely able to account for the complexity of the canopy reflected signal at fine spatial resolution as this depends also on canopy geometry (leaf inclination and spatial arrangement, plant row orientation and spacing), leaf and soil optical properties, sun position and viewing position [67,68].A physical modeling technique was used in this study.Top-of-Canopy (atmospherically corrected) reflectance data in four wavebands (green, red, red edge, near-infrared) were used in the inversion of a leaf/canopy reflectance model to produce maps of LAI and CCC.Field data were not required for the estimation of biophysical variables and could be used exclusively for product validation, therefore avoiding circularity in the next stage of the process when an empirical relationship is derived between MTCI and CCC using chlorophyll maps derived with the same combination of wavebands.
The widely used coupled PROSPECT+SAILH model ('PROSAIL') was applied here.PROSAIL is a combination of the leaf model PROSPECT-4 [69] and the canopy model SAILH [61,70,71].It calculates the bi-directional reflectance of homogeneous canopies as a function of several structural and biophysical variables, soil reflectance, illumination and viewing geometry.The PROSAIL model inversion was based on a LUT approach.The model was parameterized using the technique described in [72] and field measurements used to validate the accuracy of resulting LAI and CCC maps.

Comparing Fine Spatial Resolution CCC with Medium Spatial Resolution MERIS L2 MTCI Data
Fine spatial resolution CCC data derived in Section 2.6 was compared with MERIS L2 MTCI.To minimize overlap errors due to pixel misalignment between the fine and the medium spatial resolution data, the fine spatial resolution map was re-projected to the same coordinate system as the medium resolution image [29].The analysis was performed at approximately 1 km (3 × 3 full resolution MERIS pixels) spatial resolution.A regular grid was generated as an overlay to MERIS pixels.Fine spatial resolution CCC data and MERIS L2 MTCI were averaged within each grid cell.

Comparing Fine Spatial Resolution LAI with Medium Spatial Resolution MERIS and MODIS LAI Products
LAI plays an important role in the validation of CCC and so an additional validation phase was used.In this step, the fine spatial resolution RapidEye satellite sensor LAI map derived in Section 2.6 was compared with the MERIS LAI and the MODIS LAI products.In the case of MERIS, the spatial aggregation was conducted as for CCC based on the 1 km grid.In the case of MODIS, a regular grid based on the position of MODIS 1 km pixels was used as the base on which to average both fine spatial resolution LAI and MODIS LAI data.The intercomparison was performed only for pixels for which there was an LAI estimate derived using the main MODIS algorithm.The selection of those pixels was achieved using the quality flag (SFC_QC = 0) which accompanies the LAI product.

Estimation of LAI and CCC from Fine Spatial Resolution Data from RapidEye Imagery: Model Inversion Results
The accuracy with which LAI and CCC can be mapped was quantified using the coefficient of determination (R 2 ), the root-mean-square-error (RMSE) and the relative RMSE (RMSE/mean value of field estimates) between satellite sensor and ground based estimations.Estimations of LAI using the PROSAIL LUT model inversion approach resulted in an accuracy of R 2 = 0.76.For CCC, the estimation accuracy was R 2 = 0.78.Correlations between both LAI and CCC estimated on the ground and from satellite sensors are presented in Figure 4.The RMSE and the relative RMSE for satellite sensor estimates of LAI and CCC variables and corresponding variables estimated in the field are presented in Table 3.

Validation of Medium Spatial-Resolution MERIS L2 MTCI Data
To compare the fine and medium spatial resolution data sets, the CCC data derived from RapidEye satellite sensor data were aggregated based on a 1 km grid using the approach described in Section 2.7.Cells in which less than 80% of the area was classified as vegetation by the fine spatial resolution data were not considered further.Maps of CCC at fine spatial resolution and medium spatial resolution MERIS L2 MTCI are given in Figure 5.

Comparison with Medium-Resolution MERIS and MODIS LAI Data
MERIS LAI was derived using the MERIS TOC-VEG processor implemented in BEAM VISAT software and compared to high spatial resolution RapidEye satellite sensor LAI data.The aggregation at 1 km spatial resolution was performed as described in Section 2.8.The relationship between MERIS LAI and LAI estimated from fine spatial resolution RapidEye satellite sensor data is presented in Figure 7.The coefficient of determination was high (R 2 = 0.80) and significant (p-value < 0.001), the RMSE resulted in 0.33 (no. of pixels 504).The evaluation of the MODIS LAI was performed only for those pixels flagged as best quality and processed by the main algorithm.The relationships between RapidEye satellite sensor LAI data and MODIS LAI estimated according to the MODIS three-dimensional model inversion approach is shown in Figure 8.The comparison is presented following the MODIS biome classification, which is based on six canopy architectural types (grasses and cereal crops, shrubs, broadleaf crops, savannah, broadleaf forest, and needle-leaf forest).RMSE and number of pixels for each biome type and all biomes combined are provided in Table 4.
Table 4. RMSE and number of pixels for the relationships MODIS LAI vs. RapidEye LAI for the six biomes of the MODIS classification scheme.RapidEye LAI data were aggregated at 1 km based on the position of MODIS pixels.The RMSE of our validation ranged from 0.30 (broadleaf crops) to 1.57 (broadleaf forest) with a combined RMSE for the study area of 0.81.The accuracy reported for MODIS LAI (MODIS land team, http://landval.gsfc.nasa.gov/) is 0.66 (RMSE units).If the broadleaf forest class is excluded the resulting accuracy is 0.50.However, this accuracy is affected by (i) uncertainties in the input of the algorithm (surface reflectance and land cover) and (ii) model uncertainties, namely the consistency between RT simulations stored in the LUT and the corresponding MODIS surface reflectance [73].Our study highlighted misclassification in the land cover product used for the LAI algorithm.For instance, 236 pixels-out of a total of 494 pixels-were classified as savannah (biome 4), whereas the majority of the land cover is irrigated agriculture and needle leaf forest (large evergreen conifer forests are present in the coastal area of the study site).The broadleaf forest class (48 pixels), which is also rarely represented in the study area, and savannah are the two land covers with lower LAI estimation accuracy.This uncertainty is probably due to misclassification and the resultant selection of an erroneous biome type for LAI modeling.

Robustness of MERIS L2 MTCI-CCC Relationships
To evaluate the robustness of MERIS L2 MTCI and CCC linear equations, different sources of uncertainty were considered.Validation of the CCC map used field estimates of LAI and leaf chlorophyll concentration, derived using indirect sampling methods.Regarding LAI field estimates, errors can be up to 10% depending on the degree of crop heterogeneity [54].Earlier studies of LAI estimation have reported that indirect methods underestimated LAI between 25% and 30% for most woody plant canopies [74].
A major source of uncertainty in the estimation of leaf chlorophyll concentration is the definition of conversion equations used to convert field readings into leaf chlorophyll concentration.An example of the scatter around relationships between SPAD-502 readings and leaf chlorophyll concentration extracted in the laboratory is presented in Figure 1 for nine leaf types.The varying nature of the SPAD-leaf chlorophyll concentration relationships (linear and exponential) and the absence of instrument-and crop-specific conversion functions necessitate the destructive sampling of a small number of leaves in order to convert relative leaf chlorophyll concentration readings into absolute values per unit fresh weight.Using the conversion equations (not shown) derived from relationships in Figure 1 (using a minimum of 24 leaves) the best estimate of leaf chlorophyll concentration can be achieved with an error ranging from 7% to 28% (with a mean error of 12%).The robustness of the conversion equations for each of the nine crops is reported in Table 5.The sampling procedure is an important consideration as leaf chlorophyll concentration varies across the vertical profile of the vegetation and across different portions of the same leaf.It is recommended that the leaf surface is dry and clean (free from dust or honeydew secreted by aphids and other insects).Large variation in the measurements can be due to leaf veins and changes in thickness due to age and position of the leaf in the vertical profile.Previous studies reported that SPAD-502 readings and leaf chlorophyll concentration may also vary with plant growth stage [75], growing conditions [76], growing season [22] and genotype [77].Thus, leaf measurement position should remain invariant within crop types.To represent the chlorophyll content in the canopy, individual leaf readings should be taken across the vertical profile of multiple plants within the same LAI sampling unit.In addition, alternative leaf chlorophyll concentration sampling strategies have been reported to increase the accuracy of CCC estimation.For instance, Gitelson et al. [78] estimated CCC throughout a growing season by measuring leaf chlorophyll concentration of only the upper leaves along with the LAI of the canopy.The CCC was estimated with an accuracy (RMSE) of 0.22 g•m −2 (R 2 = 0.97) in maize and 0.23 g•m −2 (R 2 = 0.92) in soybean.
An additional source of uncertainty is in the up-scaling procedures from field measurements to fine spatial resolution biophysical maps.In this study, the high-resolution CCC map was generated with a RMSE of 0.39 g•m −2 relative to the ground measurements from 36 ESUs (Relative RMSE = 34%).This error might be due to limitations in the biophysical variable retrieval algorithm and to the ill-posed inverse problem.
Land surface heterogeneity and different plant species composition might also influence the validation results.For example, an earlier field experiment [28] provided a slightly different function compared to Equation (1) (the two equations differ with a p-value = 0.022): CCC In [28] field data were acquired for six vegetation species in large (>25 ha) homogenous fields.LAI and leaf chlorophyll concentration at 25 random sampling locations in each field were used to derive the chlorophyll content per field and then related to the mean MTCI value for full resolution MERIS pixels within each field.

Conclusions
In this study, the performance of medium spatial resolution MERIS Terrestrial Chlorophyll Index Level 2 product (MERIS L2 MTCI) to estimate canopy chlorophyll content (CCC) was addressed for heterogeneous land surfaces in an agricultural region in Southern Italy.
We applied a methodology to produce a fine spatial resolution map of CCC using the inverse modeling techniques with RapidEye satellite sensor data.The validation of this map was performed with independent ground measurements.Subsequently, the CCC map (6.5 m) was aggregated to medium spatial resolution based on a regular grid (1 km) to validate MERIS L2 MTCI.A relationship between MTCI and CCC was established for land surfaces with a vegetation cover >80% at a 1 km pixel scale.Different sources of uncertainty from limitations in field measurements to scaling issues were considered and discussed.Additionally, the study concentrated on the validation of medium spatial resolution MERIS and MODIS leaf area index (LAI) products considering the same scaling approach.
In summary, results from the validation and uncertainty lead to the following conclusions: • Medium spatial resolution MERIS L2 MTCI was validated based on fine spatial resolution CCC data derived from RapidEye satellite sensor data.The relationship between MERIS L2 MTCI and CCC [g•m −2 ] achieved a coefficient of determination of 0.74 and it resulted extremely statistically significant (p-value < 0.001).• Inverse modeling techniques based on a Look-Up-Table (LUT) approach were used to derive fine spatial resolution maps (6.5 m) of CCC with R 2 = 0.78 and RMSE = 0.39.Similarly, LAI achieved an accuracy of R 2 = 0.76 and RMSE = 0.64.The validation was based on independent LAI and leaf chlorophyll concentration estimated in 36 sampling units covering different crop species.Leaf chlorophyll concentration was estimated using the SPAD-502 chlorophyll meter.
• Based on indirect SPAD-502 measurements and crop and instrument specific conversion functions, leaf chlorophyll concentration can be estimated with an error ranging from 7% to 28% (with a mean percentage error of 12%).Wherever possible, the study suggested the use of instrument-and crop-specific conversion functions.
• The fine spatial resolution LAI map was compared to well-established medium resolution products namely MERIS LAI (based on the TOC-VEG processor implemented in the ESA BEAM VISAT software) and MODIS LAI (MOD15A2).Results indicated a better accuracy in LAI estimation of MERIS (R 2 = 0.80 and RMSE = 0.33) compared to MODIS (R 2 = 0.27 and RMSE = 0.81).It was not surprising that MERIS data had a similar performance to our fine spatial resolution LAI product since both maps are based on the inversion of the same radiative transfer model.Moreover, MERIS benefits of a finer spectral and higher spatial configuration compared to MODIS.On the contrary, misclassification in the MODIS land cover biomes used for LAI modeling might be the source of uncertainty for LAI estimation.
This work provides the framework to evaluate the robustness of MERIS L2 MTCI equations and its accuracy to predict CCC.The validation protocol has been described and it will be used in future field campaigns as a standard approach to assess the performance of MTCI.Validation methodologies and quantification of uncertainties are particularly relevant in the framework of the GMES (Global Monitoring for Environment and Security) initiative, since land surface products (such as MTCI) will be distributed directly to users or will be used to provide value added services.In particular, the Sentinel-3 payload will include the Ocean and Land Colour Instrument (OLCI) and will provide the OLCI Terrestrial Chlorophyll Index (OTCI) as Level 2B main product based on the MERIS L2 MTCI heritage.
Acquisition of fine spatial resolution (6.5 m) RapidEye satellite sensor data Geometric correction and estimation of Top-Of-Canopy reflectance Retrieval of LAI and CCC maps using canopy reflectance model inversion Validation of LAI and CCC maps with indirect ground measurements Aggregation of CCC and MERIS L2 MTCI maps at a common 1 km grid and validation Acquisition of medium spatial resolution (300 m) MERIS L2 data and extraction of MTCI data last years (R 2 = 0.64) [48].LAI and leaf chlorophyll concentration were sampled within ESUs over a period of three days (from 23 to 25 August 2009) for a range of crops including forage crops (1 alfalfa and 11 maize ESUs), fruit trees (1 plum, 2 apricot, 5 kiwi and 11 peach ESUs) and vegetables (1 aubergine, 1 pepper, 1 artichoke ESU).In addition two high-biomass poplar plots were also sampled.

Figure 3 .
Figure 3.The study was carried out in Campania region in Southern Italy (left).The test site was located in the 'Piana del Sele' (right), one of the largest agricultural areas in the region.The gridded box represents the footprint of the 4 RapidEye image tiles used in the study.

Figure 4 .
Figure 4. LAI (left) and CCC (right) estimated from field measurements (X axis) versus RapidEye satellite sensor data (Y axis) estimates.Symbols correspond to: '•' fruit trees, '□' maize and '◊' other vegetation types.Error bars in LAI indicate standard deviations of the LAI estimate.

Figure 5 . 1 )Figure 6 .
Figure 5. Fine spatial resolution CCC [g•m −2 ] map from RapidEye satellite sensor data (left) and medium spatial resolution MERIS L2 MTCI data (right).The two maps were aggregated and compared at 1 km spatial resolution based on the position of MERIS pixels.

Figure 7 .
Figure 7.Comparison between MERIS LAI (derived using the TOC-VEG processor, BEAM VISAT) and LAI estimated from RapidEye satellite sensor data aggregated at 1 km spatial resolution.

Table 1 .
Conversion equations used to convert relative SPAD-502 values to leaf chlorophyll concentration.

Table 2 .
Mean values (and standard deviation) of LAI and CCC field estimates.

Table 3 .
RMSE and relative RMSE for RapidEye satellite sensor estimates of LAI and CCC variables and corresponding variables estimated in the field.

Table 5 .
[47,48]ient of determination (R 2 ), standard error and number of observations of the equations (not shown) derived from destructive sampling to convert SPAD-502 values in leaf chlorophyll concentration values.The data set is part of past experimental campaigns and destructive chlorophyll extraction was performed according to[47,48].