Sensitivity of X-Band ( σ 0 , γ ) and Optical ( NDVI ) Satellite Data to Corn Biophysical Parameters

The objective of this work was to evaluate the sensitivity of three different satellite signals (interferometric coherence (γ), backscattering coefficient (σ0) and NDVI) to corn biophysical parameters (leaf area index, height, biomass and water content) throughout its entire vegetation cycle. All of the satellite and in situ data were collected during the Multi-spectral Crop Monitoring (MCM’10) experiment conducted in 2010 by the CESBIO Laboratory over eight different agricultural sites located in southwestern France. The results demonstrated that the NDVI is well adapted for leaf area index monitoring, whereas γ27.3° is much more suited to the estimation of the three other Biophysical Parameters throughout the entire crop cycle, with a coefficient of determination ranging from 0.83 to 0.99, using non-linear relationships. Moreover, contrary to the use of the NDVI or backscattering coefficients, the use of coherence exhibited a low sensitivity to the changes in vegetation and soil moisture occurring during senescence, offering interesting perspectives in the domain of applied remote sensing.


Introduction
Corn, which is primarily used for human and animal consumption, is the second-most-produced crop in the world (FAO, http://faostat.fao.org/).As food requirements increase in the future, the monitoring and manage-ment of corn crops to meet demand of resources will become a social issue.Satellites offer the possibility to monitor corn worldwide due to their high coverage and high temporal sampling.Among the many choices of operational satellite sensors, this work focuses on three high-resolution optical and radar imagers (pixel size < 20 m) that are able to monitor plots of corn ranging in size from a few to several hundred hectares.
In the optical domain, many studies have highlighted the interest in using satellite data for crop monitoring.These studies have primarily used the Normalized Difference Vegetation Index (NDVI) to estimate crop functioning parameters such as the Leaf Area Index (LAI), yield, soil water demand, fraction of cover, and fraction of active photosynthetic radiation [1]- [4].The primary drawbacks of these optical approaches are the limited periods of clear skies and the lack of sensitivity to crop geometry/architecture.Fewer studies have addressed crop monitoring in the microwave domain than in the optical domain, and majority of them used C-band data.These studies have explored a wide range of radar signal processing features, such as the backscattering coefficient (σ 0 ), polarimetric indexes (α angle, entropy, anisotropy) and/or interferometric coherence (γ), for crop classification and crop parameter estimation [5]- [12].In our knowledge, no studies focused on estimating corn biophysical parameters using SAR images acquired by satellite, particularly at X-band wavelengths.
The combined use of the optical and microwave domains within the electromagnetic spectrum began to increase at the beginning of the 21 th century with the launch of synthetic aperture radar (SAR) sensors (such as Envisat, Terrasar-X, Radarsat-2, and Alos), which were in limited use before 2003.Recent studies have demonstrated an interest in combining different sources of satellite data for crop management or monitoring, although there have only been a few studies on this topic [13]- [18].The primary limitation of such an approach is the lack of a well-adapted dataset that contains both in situ and satellite data to better explore the synergistic use of optical and radar remote sensing satellite data, particularly for corn [19] [20].
In this context, the objective of this work is to assess and compare the sensitivity of optical (NDVI) and X-band SAR (backscattering coefficients and interferometric coherence) satellite data to changes in corn biophysical parameters.Four parameters have been defined as target parameters-LAI, height, biomass and water content-which are all considered to be key parameters for food production and water resource management.Whereas the three first parameters can be used to derive the grain yield or leaf nitrogen concentration of the crop (based on empirical or semi-empirical approaches; [2] [21]- [23], the fourth generally serves to quantify the water stress of the plant [24]. This study is based on ground and satellite data acquired during the Multi-spectral Crop Monitoring experiment conducted by the CESBIO Laboratory in 2010 (Section 2) [20].The methodology is presented in Section 3. The results and discussion are provided in Sections 4 and 5.The conclusion and perspectives are presented in Section 6.

Study Site
The study site is located in southwestern France in the Midi-Pyrénées region (coordinates for the center of the area are 43˚29'36"N, 01˚14'14"E) (Figure 1).It is strongly anthropized, governed by a temperate climate and primarily conditioned by agricultural activities [20].During 2010, the annual rainfall reached 600 mm, and the mean daily air temperature ranged from 3.5 to 22˚C.Approximately 10% of the surface is attributed to urban, free water and forest areas.The remaining landscape is composed of grasslands and annual crops (e.g., corn, sunflower, wheat).The eight contrasted studied plots of corn are distributed inside the swath of the satellite images, as illustrated in Figure 1.The region of interest is extremely flat, with mean slopes of less than 0.4%, and the soil texture over the 8 studied plots is mainly defined as silty loam (silt = 54.4%,sand = 24.4%,clay = 21.2%),following the USDA classification scheme (Table 1).In 2010, corn was sown between April 2 nd and May 23 rd to be harvested at the beginning of the autumn (from September 7 th to October 27 th ).Plot boundaries were decreased by 10 m to reduce mixel effects.

Data Collection
Ground and satellite data were obtained from the Multi-spectral Crop Monitoring experiment (MCM'10) conducted in 2010 by the CESBIO Laboratory.2).Optical images were provided by high-resolution sensors onboard Formosat-2 and Spot-4/5 satellites in the visible to short-wave infrared domains (Table 2).The ground spatial resolutions were equal to 8, 20 and 10 m.All of the images were ortho-rectified using CNES ortho-rectification tools.Cloud detection and atmospheric corrections of the Formosat-2 images were applied, and the Spot images were atmospherically corrected using the CNES Kalideos processing chain.
Pre-processing of SAR and optical satellite data are described in Figure 2.

In Situ Data
In situ data included measurements of surface soil roughness, inter-row distance, Crop Height (CH), aboveground Crop Dry and Fresh Biomass (CDB and CFB, respectively), Crop Water Content (CWC) and Surface Soil Moisture (SSM).The two first parameters were considered stable over time (i.e., during the crop growth cycle).Soil roughness was measured using a 2-m pin profilometer, and the mean statistical parameters were calculated immediately after sowing, when the soil was smooth (h rms = 1.2 cm, l c = 4.5 cm).The inter-row distance was constant between the plots and was equal to 80 cm.
The four remaining parameters, which vary substantially over space and time, were collected throughout the vegetation cycle, from emergence to harvest.The dielectric constant of the soil was measured using mobile theta probe sensors, and the volumetric soil moisture of the first five centimeters was estimated using an in situ calibration function (R 2 = 0.75, RMSE = 4.1%, n = 403, [20]).For the 8 plots, the mean SSM value ranged between 9 and 29% m 3 •m −3 throughout the crop cycle.
The mean vegetation height was derived from measurements performed over more than ten samples for each plot.Plant height reached its maximum value (approximately 3.5 m) 100 days (on average) after sowing (Figure 3).
Regular biomass measurements were performed during the corn agricultural season over two plots (plot ID #7 and #8).For each date, five samples were weighed to obtain the fresh biomass, dried in an oven, and reweighed to determine the dry biomass.The total crop water content (including leaves, stems, and fruit) was derived from the difference between the CFB and the CDB, and the values were related to the phenological cycle of the crop.From emergence to flowering (approximately on DOY #220), the water content reached a maximum of 80% and decreased during senescence, reaching 50% only a few days before being harvested (the WC of the fruit was less than 30%).
All of the in situ variables were collected quasi-synchronously with satellite acquisitions (a mean time lag of less than 1 day was observed).

Geometric corrections NDVI calculation
Backscattering coefficient calculation

InSAR treatments
Coherence calculation

Opticalimages SAR images
Radiometric and atmospheric corrections

Methodology
Among

BP a c
where BP represents the biophysical parameter under consideration (LAI, CH, CFB, CDB, CWC and SSM), and a BP , b BP and c BP are the associated empirical parameters.
The relationships between biomass (wet and dry), crop water content, crop height, LAI surface soil moisture content and satellite signals were evaluated.The NDVI was estimated from three optical satellites (Formosat-2, Spot-4/5).The similarity of the NDVI values was determined by comparing the values estimated by each satellite from images acquired at the same time ( The backscattering coefficients of the Terrasar-X images (SL and SM modes) were estimated using the radiometric calibration procedure given in [16] and performed using NEST (Next ESA SAR Toolbox) software.Backscattering coefficients acquired at HH polarization ( 0 HH σ ) were derived from the digital number (DN) of the pixel (i), the calibration factor (K) and the incidence angle (θ) (Equation ( 2)).

( ) ( ) ( ) ( ) ( )
0 dB 20 log 10 log 10 log sin Geocoded interferometric coherence images were calculated from InSAR treatment processing based on the approach described in [25] and using open-source software, namely, DORIS (Delft Object-oriented Radar Interferometric) and ADORE (Automated DORIS Environment).Master and slave images were defined for each SAR mode (SM and SL) to calculate the coherence coefficients (γ) over the slave images, following Equation (3): where M i and S i represent the value of the pixel i of the master and slave images, respectively, * i M and * i S denote the complex conjugate matrix, and N denotes the total number of analysed pixels.
As described in [26], decorrelation sources can be separated in four factors: γ = thermal geometrical volume temporal γ γ γ γ ⋅ ⋅ ⋅ .Thermal decorrelation can be neglected for high signal-to-noise ratio [27] [28].Geometric decorrelation was reduced by applying a spectral filter, and the spatial baselines of all of the Insar pairs were verified as being lower than their critical baselines (Table 3).Volumetric and temporal coherences cannot be separated as a simple sum of independent coefficients.As demonstrated and modelled by [29], the second component is more affected by the vertical change in vegetation.Therefore, we used the coherence coefficient to estimate the crop biophysical parameters over time.
Ten and seven coherence images were generated in the StripMap and SpotLight modes, respectively.All of the images were geo-referenced using ortho-photos (spatial resolution of 50 cm) provided by the French National Geographic Institute, and the final location accuracy was less than 2 pixels, on average, considering the different products.Examples of NDVI, γ and 0 HH σ images are provided in Figure 4 for a sub-area of the study site from May 2010.

Influence of Surface Soil Moisture
The first results confirmed that no accurate empirical relationships could be estimated between the surface soil moisture changes and the satellite signals over the entire vegetation phenological cycle (0 < R 2 < 0.13) (Table 4).These weak correlations are attributed to the fact that the satellite signals were mainly acquired in the presence of vegetation.In these conditions, the X-band and optical electromagnetic waves did not penetrate through the vegetation layer but mainly interacted with the upper layer of the crop (i.e., leaves, stems, flowers) and not with the SSM.Nevertheless, according to [30] [31], stronger correlations can be observed when the soil is bare, from sowing to the first phenological stages of corn.

Sensitivity of Satellite Signal to LAI
The values of the statistical parameters associated with the relationships estimated between the LAI and satellite signals were disparate.The results associated with the optical signal confirmed that NDVI is well suited for LAI  retrieval, as previously noted in several studies [2].The radar signals (either backscattering coefficients or interferometric coherence) exhibited contrasting performances, with R 2 values ranging from 0.03 to 0.51 (Table 5).
According to the relationships based on the backscattering coefficients, the performances depended on the considered sensor configuration, with higher R 2 values associated with a high incidence angle (53.3˚ in this case).This result appears consistent with previous studies [32], showing that the sensitivity of the backscattering coefficients to soil and vegetation parameters depend on the considered incidence angle.

Sensitivity of Satellite Signal to Biomass and Crop Water Content
The performances associated with the relationships between corn biomass and satellite signals showed similar trends regardless of the considered fraction of the crop (CFB, CDB and CWC) (Table 6).The poorest values of R 2 were obtained for the optical index (NDVI) or backscattering coefficients (regardless of the incidence angle), with values never exceeding 0.37, whereas the highest R 2 values were obtained using interferometric coherences (regardless of the incidence angle), with values of R 2 close to 0.99.Despite these strong correlations, the application of these relationships is limited to the first phenological stages of corn due to the early saturation of the interferometric coherence (Figure 5).The results demonstrated that the saturation of the coherence occurred at a later phenological stage with the decrease in incidence angle, regardless of the considered crop parameter (CFB, CDB and CWC).For example, the saturation level was close to 1080 g•m −2 using γ 53.3˚ (beginning of stage 1: leaf development, BBCH), whereas this level reached 1738 g•m −2 using γ 27.3˚ (end of stage 3: stem elongation, BBCH) for the CFB.These observations can be interpreted based on the configuration of the X-band data and the associated sensitivity.The sensitivity of backscattering coefficients are influenced more by vegetation in images acquired at high incidence angles (>30˚) than at low incidence angles due to volume and double bounce contributions.At low incidence angles, backscattering coefficients are influenced more by a single scattering mechanism due to the soil contribution (associated with soil moisture and soil roughness changes).Thus, the correlation with crop height is lower at a low incidence angle than at a high incidence angle.According to the coherence, the lowincidence images are more strongly correlated with crop height, which is contrary to the results obtained from the backscattering coefficients (where soil moisture values more strongly affect the signal).An explanation of this phenomenon is that as the crop grows, the soil contribution decreases (R 2 < 0.13 between coherence and SSM) with the increase in vegetation contribution, which decorrelates more rapidly than soil, especially close to the nadir (27.3˚),where the vertical dynamic of the vegetation layer strongly affects the coherence.The decrease in the coherence over time can be attributed to the vegetation layer and not to the soil, which is expected to be coherent for an interval of a few days.The overall pattern of coherence decreases with the increase in crop height until saturation, inducing an increase in incoherent scattering mechanisms (Figure 6).The non-linear shape of the relationship can be explained by the progressive decrease in the soil contribution (for which the correlation is high), to the benefit of the incoherent vegetation contribution.
Figure 6 presents the best results obtained with optical and X-band satellite data.The relationships indicate comparable behaviors with two different periods: (i) a high sensitivity with an exponential decrease in the case of interferometric coherence and an exponential increase for the NDVI and backscattering coefficients, and (ii) a signal saturation phase for height values depending on the satellite signal under consideration.This level of saturation (highlighted in Figure 6 by the point located at 90% of the overall relationship saturation) is reached for corn heights near 75, 150 and 160 cm for the relationships based on backscattering coefficients, interferometric coherence and NDVI, respectively.Figure 6 also illustrates that surface changes (SSM, CWC, phenological stage), which occurred after 90% signal saturation, have a smaller effect on the coherence relationship than on the relationships based on the backscattering coefficients and/or NDVI.In addition, the results (figure not shown here) confirmed that parallel and/or perpendicular baselines values do not affect the values of coherence because they are inferior to the critical baseline.

Discussion
As expected, none of the tested satellite signals can be used to estimate the surface soil moisture changes throughout the vegetation cycle, especially when the presence of the vegetation layer masks the soil layer.This well-known phenomenon can be explained by the low penetration depth of the wavelength (both in optical and X-band frequencies) into the crop layer.The estimation of SSM under the crop layer remains challenging and is presently investigated using the POLSAR technique (giving access to the double bounce mechanism, partially associated with the SSM) or L-band signals (offering higher penetration depth and thus a stronger correlation with SSM) [33] [34].The results obtained in this paper confirm that techniques based only on backscattering coefficients, coherence or optical signals are not well adapted for monitoring SSM under the crop layer, and they reinforce the interest in developing new approaches such as those proposed by [34].However, the weak correla- The analyses of the relationships estimated between the satellite signals and the vegetation parameters revealed several important points.First, the results showed that among all of the tested satellite signals, optical data appear more suited to monitor the LAI of corn with sufficient confidence.Considering X-band signals, only the backscattering coefficient acquired at high incidence angle (53.3˚) could justifiably be used to estimate the LAI (with lower confidence than with optical data).Therefore, the monitoring of corn should be performed using a multi-temporal approach mainly based on optical data and partially composed of radar data, the latter being used only to replace optical data in case of cloud cover.Under these conditions, the combine use of optical and radar data could be used to accurately estimate some variables strongly correlated with the LAI, such as leaf transpiration, the interaction with the radiation by the leaf, and latent and sensible heat fluxes [35] [36].
The second part of the results demonstrated that the dry or fresh biomass and crop water content are strongly correlated over time.Consequently, considering that crop water content is a key parameter in the management of agricultural crops, the monitoring of dry or fresh biomass could be investigated to further estimate the water content of the crop, particularly to determine the maturity of the crop before harvest [37].Concerning the estimates of CDB, CFB and CWC, the best correlations were achieved using the coherence estimated from images acquired at 27.3˚ (R 2 = 0.99).Although these relationships quickly saturated (for values of fresh biomass close to 1738 g•m −2 ), the results (Figure 5) demonstrated that coherence is strongly sensitive to the first phenological stages of the crop (from emergence to the end of the stem elongation, as defined in the BBCH classification), contrary to other satellite signals, which were too strongly affected by SSM changes during the bare soil period or by the color changes and decrease in CWC during senescence.This accurate and early detection of the first phenological stages could be then implemented to estimate crop yields by using an agro-meteorological model, as demonstrated by [2] [38]- [42].
The estimate of CH demonstrated the potential of using the coherence instead of the other satellite signals, especially from images acquired at low incidence angle (27.3˚).Consequently, the coherence values could be implemented in empirical and/or semi-empirical methods for predicting yields as demonstrated by [21] [22].Moreover, the estimate of CH could be used as a key parameter in land-use mapping to more precisely differentiate similar crops where architectural differences mainly accrue from height differences [43].For example, the method could be investigated to improve corn and sorghum classification given their heights of approximately 1 and 3 m, respectively, over the study region.

Conclusions
This study evaluated the sensitivity of three sources of satellite data (X-band backscattering coefficients, interferometric coherence, and NDVI) to crop parameters over an entire vegetation cycle.The first part of the results demonstrated that the SSM changes were not correlated with satellite signals due to the presence of vegetation and the low penetration depth of the electromagnetic wave into the crop layer (regardless of the spectral domain).Consequently, we demonstrated that SSM did not significantly affect the relationships estimated between the satellite signals and the vegetation components (CDB, CFB, LAI and height).The second part of the results demonstrated that the LAI of corn can be estimated well using the NDVI (as was known from previous studies) and partially substituted by the use of backscattering coefficients acquired at a high incidence angle (53.3˚).The results also demonstrated the potential of using the coherence acquired at a low incidence angle (despite the early saturation of the coherence) to estimate the dry or fresh biomass and the water content, especially for the first phenological stages of the corn (R 2 = 0.99).The best corn height estimates were obtained by considering the interferometric coherence extracted from the images acquired at a low incidence angle (27.3˚) due to its higher decorrelation effects.In this case, crop height can be accurately estimated before 150 cm (R 2 = 0.83, RMSE = 56 cm, rRMSE = 32%).The use of backscattering coefficients at a high incidence angle and the use of NDVI were less accurate, even though the coefficients of determination were higher than 0.67 (the relative errors, which ranged between 78 and 92%, are still important).Moreover, the use of coherence is of particular interest for monitoring vegetation parameters because it is less sensitive to SSM and CWC than are σ 0 or NDVI during senescence.
In the near future, these relationships could be used in empirical or semi-empirical approaches to estimate crop yield at the landscape scale as presented in the "Discussion" section.Moreover, these results must be extended to other optical and X-band satellite missions such as Sentinel-2, COSMO-SkyMed, TamDEM-X and PAZ.

Figure 1 .
Figure 1.Locations of the eight monitored plots of corn and the climatic station of Lamasquère.

Figure 2 .
Figure 2. Flow chart of the satellite data pre-processing.

Figure 3 .
Figure 3. Temporal change in crop height over the 8 plots of corn (a) and in the fresh and dry biomass (b) measured over plot ID #7 (cross marker) and #8 (circle marker).
the wide choice of Biophysical Parameters (BP), the proposed methodology focuses on the relationships estimated between crop and soil parameters (LAI, CH, CFB, CDB, CWC, and SSM), and three satellite indicators for the entire crop cycle: NDVI, γ 27.3˚, γ 53.3˚,

Figure 5 .
Figure 5.Estimated relationships between the corn dry biomass (a), fresh biomass (b), water content (c) and interferometric coherence acquired at 27.3˚.Ninety percent of signal sensitivity is marked by vertical and horizontal dotted black lines for each sub-plot.

Figure 6 .
Figure 6.Estimated relationships between corn height and backscattering coefficients (a), interferometric coherence (b) and NDVI (c); 90% sensitivity is marked by vertical and horizontal dotted black lines for each sub-plot.

Table 1 .
Textural and topographic characteristics of the studied plots of corn.

Table 3 .
Summary of interferometric parameters for the two considered modes (SM and SL) for each pair of master and slave images: critical baseline (B critical ), perpendicular (B perp ), parallel baseline (B par ) and temporal baseline (B temp ).

Table 4 .
Values of the empirical parameters (a SSM , b SSM and c SMM ) used to determine the SSM, number of data points used (n), coefficients of determination ( 2 SSM R ), and absolute and relative root mean square errors (RMSE SSM and rRMSE SSM , respectively).

Table 5 .
Values of the empirical parameters (aLAI, bLAI and cLAI) used to determine the LAI, number of data points used (n), coefficients of determination ( 2 LAI R ), and absolute and relative root mean square errors (RMSELAI and rRMSELAI, respectively).

Table 6 .
Values of the empirical parameters (a, b and c) used to determine the fresh biomass, dry biomass and biomass water content, number of data used (n), coefficients of determination (R 2 ), and absolute and relative root mean square errors (RMSE and rRMSE, respectively).Linear relationships were applied when the value of c was not provided.

Table 7 )
. The relationships between CH and γ 27.3˚ exhibited the best statistical index value, with R 2 values greater than 0.8 and a maximum rRMSE of the corn height estimates of 32%.The performance of the other indicators decreased gradually for relations based on 0 σ −  .The associated R 2 and rRMSE values ranged from 0.71 to 0.04 and from 78 to 102%, respectively.

Table 7 .
Values of the empirical parameters (aCH, bCH and cCH) used to determine the corn height, number of data used (n), coefficients of determination ( 2 CH R ), and absolute and relative root mean square errors (RMSECH and rRMSECH, respectively).