Sensitivity Analysis of Multi-Temporal Sentinel-1 SAR Parameters to Crop Height and Canopy Coverage

The Polarimetric Synthetic Aperture Radar technique has provided various opportunities and challenges in agricultural activities mainly on crop management. The aim of this study is to investigate the sensitivity of 10 parameters derived from multi-temporal Sentinel-1 Synthetic Aperture Radar (SAR) data, to crop height and canopy coverage (CC) of maize, sunflower, and wheat. The correlation coefficient values indicate a high correlation for maize during the early growing stage. The coefficient determinations (R2) of 0.82 and 0.81 indicate that there is a strong relationship between the maize height and SAR parameters including VV + VH and VV, respectively. The maize CC is well correlated with VV parameter (R2 = 0.73), but it is observed that at the later growing stage the correlation became weaker. This means that the sensitivity decreases with increasing vegetation cover growth. Compared to maize, the sensitivity of SAR parameters to wheat variables is often good at the early stage. However, the highest correlation with wheat height represented by Alpha (α) decomposition parameter (R2 = 0.67). The sunflower height has an insignificant correlation with the majority of SAR parameters and only VH polarization shows low sensitivity (R2 = 0.31). The sunflower CC shows relatively higher correlation with VV polarization (R2 = 0.46) at the early stage while no considerable correlation is observed at the later stage. It is found that Sentinel-1 has a high potential for estimation of crop height and CC of the maize as a broad-leaf crop. The same is not true for sunflower as another broad-leaf crop.


Introduction
Spatio-temporal monitoring and evaluation of crop types as food resources has substantial importance in socio-economic issues of societies and it is essential for the sustainable management of agricultural activities. Thus, decision makers can have accurate and updated information regarding crops management; agricultural production planning, food security stability promotion, crop trading market etc. Besides, farmers also benefit from obtaining timely information of crop growth as well as yield estimation. The superiority of radar remote sensing techniques and capabilities of Synthetic observation, crop identification and yield estimation [7,16]. CC is used to estimate the influence area of the plant and defined as the percent of a fixed area covered by the crown of an individual plant species to the bare soil in a unit area. The total CC can reach a hundred percent (fully covered) since plants can overlap. Liao et al. [17] investigated the sensitivity of RADARSAT-2 polarimetric SAR (C-band) data using sixteen parameters to crop height and fractional vegetation cover (FVC) variables of corn and wheat. They explored that the corn height and FVC are strongly correlated with SAR parameters at the early growing stage but the correlation is low at the later growing stages. Furthermore, they observed that the sensitivity of SAR parameters to wheat variables (height and FVC) is very low. The study concluded that RADARSAT-2 polarimetric SAR (C-band) data has high potential in crop variables estimation for broad-leaf crops.
The Sentinel-1 mission is designed as a two-satellite constellation (Sentinel-1A and -1B) that equipped with a dual-polarized C-band SAR 20 m spatial resolution, and six days revisit time. Besides, the data can be acquired in different operational modes which enable users to access freely available long-term data archive for applications requiring long-range time series analyses. Abdikan et al. [3] utilized Sentinel-1 dual polarimetric C-band SAR images to study the efficiency of the sensor backscatter patterns for crop growth of maize in Turkey and revealed that the Sentinel-1 dual-polarized data provides useful information regarding plant development for agriculture monitoring and valid crop mapping. The potential use of Sentinel-1 radar data has been studied by Bousbih et al. [18] to estimate the soil characteristics and vegetation variables in agricultural areas. They reported that the sensitivity of the Sentinel-1 measurements to crop parameters decreases with increasing crop cover growth and variation in vegetation parameters.
This study investigates the potential of Sentinel-1 dual polarimetric SAR in estimation and monitoring of crop parameters, namely crop height and CC in an agricultural area. The objective of this study is to evaluate the sensitivity of different Sentinel-1 dual polarimetric SAR parameters to crop height and CC of maize, sunflower and wheat, and to investigate the changes in SAR backscatter arising from crop height and CC during crop phenological stages.

Materials and Methods
Within this framework, we have investigated the sensitivity of 10 parameters including linear polarization backscattering coefficients, H-A-α decompositions, polarization intensity ratios (VH/VV, VV/VH and VV − VH/VV + VH), Radar Vegetation Index (RVI) and intensity arithmetic calculations (VH − VV and VV + VH) derived from multi-temporal C-band Sentinel-1 SAR data, to crop height and Canopy Coverage (CC) of maize, sunflower, and wheat. For this purpose, field measurements were carried out simultaneously with SAR data acquisitions. For backscattering analysis four Sentinel-1 SAR images in descending orbit direction were acquired throughout the same growth season of the study area.

Site Description
The Konya basin (38 • 40 N, 32 • 26 E) in central Turkey is selected for field measurements and satellite images collection ( Figure 1). The terrain of the study area is partly flat with a gently sloping (2%-6%) and, the smallest field area of approximately 0.5 and the largest 18 ha. in size. According to the Ministry of Agriculture and Forestry of Turkey, the distribution of major soils in the study region is Reddish Brown and Brownish soils Figure 2. The soil texture in the study area consists of clayed loam and loamy (medium structure), slightly alkaline, salt-free, and low contents of organic matter (1.30-2.08%). The region has an arid to semi-arid predominant weather conditions. The land use type is mainly agricultural land in the study site ( Figure 3). Therefore, the demand for water consumption for irrigation is increasing due to the extent and dense of agricultural activities. Maize, sunflower, and wheat are three investigated crops patterns as they have different structures. Maize and sunflower based on field campaign, generally are planted at the beginning of May, and harvested at the end of August or the beginning of September in this study area. Winter wheat is seeded in previous November and harvested at the end of July.   The region has an arid to semi-arid predominant weather conditions. The land use type is mainly agricultural land in the study site ( Figure 3). Therefore, the demand for water consumption for irrigation is increasing due to the extent and dense of agricultural activities. Maize, sunflower, and wheat are three investigated crops patterns as they have different structures. Maize and sunflower based on field campaign, generally are planted at the beginning of May, and harvested at the end of August or the beginning of September in this study area. Winter wheat is seeded in previous November and harvested at the end of July.  The region has an arid to semi-arid predominant weather conditions. The land use type is mainly agricultural land in the study site ( Figure 3). Therefore, the demand for water consumption for irrigation is increasing due to the extent and dense of agricultural activities. Maize, sunflower, and wheat are three investigated crops patterns as they have different structures. Maize and sunflower based on field campaign, generally are planted at the beginning of May, and harvested at the end

Field Surveys
In situ measurements were conducted for maize, sunflower and wheat fields in the springsummer agricultural season of the year 2016. Crop variables and parameters which indicate the growth rate including crop height and CC were collected simultaneously with the SAR data acquisition, and recorded during the field works. Field data collection includes measuring the row and plant cover by still tape in unit area, taking photographs by the camera and recording field characteristics such as soil properties and irrigation status. For calculating CC, photographs were taken in downward position and perpendicular to the ground with 100 cm distance from camera and canopy outmost in the tillering stage when the crop height was less than 100 cm. Considering the study area, 36 test fields that consist of 19 maize, 6 sunflower fields, and 11 wheat fields were dedicated for this research. The variability of development stages for maize, sunflower and wheat has been defined by applying the "Biologische Bundesanstalt, Bundessortenamt, and CHemische" (BBCH) [20] indicator for each field that generally consist of vegetative and reproductive stages. Seasonal maize, sunflower, and wheat crop calendar of the study area is presented in Figure 4. According to the defined height and CC thresholds, derived from SAR backscattering and BBCHscale, we call two stages for crops growth season; the early stage and the later stage.

Field Surveys
In situ measurements were conducted for maize, sunflower and wheat fields in the spring-summer agricultural season of the year 2016. Crop variables and parameters which indicate the growth rate including crop height and CC were collected simultaneously with the SAR data acquisition, and recorded during the field works. Field data collection includes measuring the row and plant cover by still tape in unit area, taking photographs by the camera and recording field characteristics such as soil properties and irrigation status. For calculating CC, photographs were taken in downward position and perpendicular to the ground with 100 cm distance from camera and canopy outmost in the tillering stage when the crop height was less than 100 cm. Considering the study area, 36 test fields that consist of 19 maize, 6 sunflower fields, and 11 wheat fields were dedicated for this research. The variability of development stages for maize, sunflower and wheat has been defined by applying the "Biologische Bundesanstalt, Bundessortenamt, and CHemische" (BBCH) [20] indicator for each field that generally consist of vegetative and reproductive stages. Seasonal maize, sunflower, and wheat crop calendar of the study area is presented in Figure 4. According to the defined height and CC thresholds, derived from SAR backscattering and BBCH-scale, we call two stages for crops growth season; the early stage and the later stage. Appl. Sci. 2019, 9, x FOR PEER REVIEW 6 of 19 Four field surveys were conducted to obtain accurate ground measurements in late May-mid June, early July, late July, and late August 2016. General growth stages categories (leaf development, stem elongation, heading and flowering), are defined according to the BBCH-scale (Table 1). For calculating the crop height of each test site, five height measurements were obtained and their mean value was calculated to represent the crop height of the relevant test field. From the test sites, plant cover and row were measured and also photographs were captured to determine and evaluate the CC percentage. A synopsis of the maize, sunflower and wheat of different growing stages is given in Figure 5.    Four field surveys were conducted to obtain accurate ground measurements in late May-mid June, early July, late July, and late August 2016. General growth stages categories (leaf development, stem elongation, heading and flowering), are defined according to the BBCH-scale (Table 1). For calculating the crop height of each test site, five height measurements were obtained and their mean value was calculated to represent the crop height of the relevant test field. From the test sites, plant cover and row were measured and also photographs were captured to determine and evaluate the CC percentage. A synopsis of the maize, sunflower and wheat of different growing stages is given in Figure 5.   Four field surveys were conducted to obtain accurate ground measurements in late May-mid June, early July, late July, and late August 2016. General growth stages categories (leaf development, stem elongation, heading and flowering), are defined according to the BBCH-scale (Table 1). For calculating the crop height of each test site, five height measurements were obtained and their mean value was calculated to represent the crop height of the relevant test field. From the test sites, plant cover and row were measured and also photographs were captured to determine and evaluate the CC percentage. A synopsis of the maize, sunflower and wheat of different growing stages is given in Figure 5.    In this study site, the BBCH-scale is considered as 53 when the maize height was in range of 120-150 cm. This growth stage of the maize is at the inflorescence emergence and heading stage. When maize height was greater than 220 cm, the BBCH-scale was 69 and represented the end of flowering. Once the sunflower height was greater than 92 cm, the BBCH-scale is found as 79, indicating the end of flowering and the inflorescence reaches full size. The BBCH-scale is considered as 59 after that wheat height reached to 53 cm and inflorescence fully emerged. We observed that different wheat height could have the similar wheat BBCH due to variation of wheat growth conditions which cause to distinction even though they are at the same phenology, in agreement with the study of Liao et al. [17]. Figure 6 shows the relationship between crop height and the BBCH-scale corresponding to the each crop principal growth stages.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 7 of 19 In this study site, the BBCH-scale is considered as 53 when the maize height was in range of 120-150 cm. This growth stage of the maize is at the inflorescence emergence and heading stage. When maize height was greater than 220 cm, the BBCH-scale was 69 and represented the end of flowering. Once the sunflower height was greater than 92 cm, the BBCH-scale is found as 79, indicating the end of flowering and the inflorescence reaches full size. The BBCH-scale is considered as 59 after that wheat height reached to 53 cm and inflorescence fully emerged. We observed that different wheat height could have the similar wheat BBCH due to variation of wheat growth conditions which cause to distinction even though they are at the same phenology, in agreement with the study of Liao et al. [17]. Figure 6 shows the relationship between crop height and the BBCH-scale corresponding to the each crop principal growth stages. Scatterplots (d-f) in Figure 6 show SAR response to crop heights. Note that the sensitivity of the SAR backscatter to CC of three different crops has been changed in different height and CC in each crop. For maize, the variation of correlation was determined when the maize CC threshold is 75%. This threshold is 85% and 60% for sunflower and wheat respectively.

Sentinel-1 Dual Polarimetric SAR Data Statement and Processing
We used Single Look Complex (SLC) Sentinel-1A interferometric wide swath (IW) data in descending pass direction ( Table 2). Sentinel-1 satellites are equipped with C-band sensor with an incidence angle range between 29.1° and 46.0° at 5405 MHz radar frequencies. The satellite obtains data with dual polarization (VV and VH) backscatter where electromagnetic waves have polarized vertically (V) for transmission and H/V for reception [21]. Data can be acquired in both ascending and descending orbit pass directions, which means that any object on the earth surface can be evaluated as two independent sights. The difference is that due to steeper incident angle of the ascending orbit therefore the mean backscatter values are rather higher in compare with descending orbit [22]. For time series analysis of SAR backscatter, four C-band polarimetric SAR data during the growth stages of the selected crops from mid-June to late August were acquired for investigation of the sensitivity of Sentinel-1 data to the crops.  Scatterplots (d-f) in Figure 6 show SAR response to crop heights. Note that the sensitivity of the SAR backscatter to CC of three different crops has been changed in different height and CC in each crop. For maize, the variation of correlation was determined when the maize CC threshold is 75%. This threshold is 85% and 60% for sunflower and wheat respectively.

Sentinel-1 Dual Polarimetric SAR Data Statement and Processing
We used Single Look Complex (SLC) Sentinel-1A interferometric wide swath (IW) data in descending pass direction ( Table 2). Sentinel-1 satellites are equipped with C-band sensor with an incidence angle range between 29.1 • and 46.0 • at 5405 MHz radar frequencies. The satellite obtains data with dual polarization (VV and VH) backscatter where electromagnetic waves have polarized vertically (V) for transmission and H/V for reception [21]. Data can be acquired in both ascending and descending orbit pass directions, which means that any object on the earth surface can be evaluated as two independent sights. The difference is that due to steeper incident angle of the ascending orbit therefore the mean backscatter values are rather higher in compare with descending orbit [22]. For time series analysis of SAR backscatter, four C-band polarimetric SAR data during the growth stages of the selected crops from mid-June to late August were acquired for investigation of the sensitivity of Sentinel-1 data to the crops.

SAR Backscatter
Several studies on the processing of radar data indicate that the phenological stages of the plant have an effect on the backscatter of the signal, and there is a significant correlation between the biophysical parameters of the plants, including height, leaf area index, vegetation mass, plant water content, and radar signal backscatter [15,[23][24][25]. Radar backscattering from vegetation is a function of both wavelengths polarization and frequency. Different frequencies and polarizations enable one to infer various and supplementary information from the single object. In agricultural radar applications, combination of polarizations (co-polarization and cross-polarization) allows analyst to extract extra information about crop characteristics.
The polarization of backscattering microwaves indicates the target structural properties and visualizes scattering characteristics of observed features [1,5]. A majority of the space borne radar systems often transmit only one polarization and receives both polarization giving rise to dual polarimetric SAR data (e.g., Sentinel-1 with VH and VV polarizations), while some collecting full polarimetric so-called quad polarization (HH, VV, HV, and VH) imagery (e.g., PALSAR, TerraSAR-X, and RADARSAT-2). Fully polarimetric SAR data is acquired using the H and V polarizations which extracted from SLC data and can be represented by a 2 × 2 scattering matrix S (Equation (1)) including polarimetric information for each individual resolution cell [26][27][28].
Scattering matrix which depends on incident and the scattered field, has four components, each representing the received and transmitted polarizations [29,30]. The scattering matrix consists of information on the nature and characteristic of the observed media and features. Full polarimetric SAR data set which is described as scattering matrix is foundation for several coherent polarimetric decomposition and analysis. For polarimetric analysis an alternative procedure is derived from a covariance matrix (C 3 ) that represents the average polarimetric information extracted from a set of neighboring pixels to produce the mean polarimetric response. The covariance matrix C 3 , (Equation (2)), is determined from the outer element of the vector form of the scattering matrix with its Hermitian conjugate, KC [17,27]. The averaged target vector (covariance matrix) for fully polarimetric data is given by Equation (3).
Appl. Sci. 2019, 9, 655 where ensemble averaging is shown by the | | represents the modulus, the * indicates complex conjugation and the complex conjugate transpose shown by superscript H. For natural targets, in case of monocratic radar, S HV ≈ S VH when Srt indicates the complex scattering amplitude for received and transmitted polarization (r, t ∈ {h, v}) for horizontal and vertical polarization and the scattering matrix is defined by three-element complex target vector, K C = S HH √ 2S HV S VV T , where superscript T indicates the matrix transpose [29,30]. In the covariance matrix, diagonal elements (C 11 = σ 0 HH , C 22 = σ 0 HV , and C 33 = σ 0 VV ) define backscattering coefficients and the upper or lower triangular components represent complex numbers. The backscattering coefficients have correlation with the structural characteristics of the features [17,31].
In comparison to the quad polarization, dual polarimetric SAR sensors collect a fraction of total (precisely half of the scattering matrix components) polarimetric information involved in fully polarimetric imagery [27]. It means that each resolution cell at each time point is defined by a 2 × 2 covariance matrix (C 2 ) that is obtained from C 3 . The resulting covariance matrix which is for dual polarization (e.g., Sentinel-1) is represented by Equation (4).
Since dual polarization has only diagonal elements, the matrix with off-diagonal components are set to zero and do not follow a complex Wishart distribution; however, the two diagonal blocks (1 by 1) do [30,32].
Polarimetric Synthetic Aperture Radar (PolSAR) technique has resulted many different investigations and improvements in crop growth monitoring, yield estimation, crop disaster prediction and prevention and in more general terms providing accurate information for precision farming. PolSAR products, such as Entropy (H), Alpha (α) and Anisotropy (A) decompositions are calculated from the covariance matrix. The H-α-A decompositions are used to extract average parameters from experimental data suggested by Cloude and Pottier [33]. This approach is based on second-order statistics using a smoothing algorithm [34]. Natural measure of the inherent reversibility of the backscattering data is defined by entropy (H), and indicates the randomness of the scatter, while the underlying average scattering mechanisms, scattering type (surface, double-bounce and volume scattering) can be identified using Alpha parameters. The relative power the second and third eigenvectors is described by Anisotropy (A), which represents being of different properties in different directions when measured along different axes [33,35]. The Entropy (H) decomposition parameter has more sensitivity to the crop parameters and the density and randomness of some vegetation canopy than Alpha and Anisotropy [28,36].
In agricultural radar monitoring, Radar Vegetation Index (RVI) is a method for observation of the level of the vegetation growth in time series data analysis as an alternative to NDVI (Normalized Difference Vegetation Index) method used in optical image processing studies [37]. Ranging between 0 and 1, RVI is used for measuring the randomness of scattering in microwave signal [38]. It is close to 0 for a smooth bare surface and as vegetation grows the value increases till the crop reaches to the end of growth cycle and it is affected by vegetation water content and sensitive to the biomass [39]. RVI calculation needs quad-polarized data, thus for full polarization, RVI is retrieved by the Equation (5).
where σ 0 HH and σ 0 VV are co-polarized backscattering coefficients and σ 0 HV is cross-polarized backscattering coefficient in power units. According to the Charbonneau et al. [40] the assumption that supposes σ 0 HH ≈ σ 0 VV then Equation (5) can be reduced to the form as Equation (6).
Melanie et al. [41] studied the RVI and concluded that RVI HH is useful when just two polarizations are available and can be an appropriate approximation of the surface scattering if the interaction between the surface plane and vegetation is insignificant.
Since Sentinel-1 is dual polarization and has VH and VV polarizations, following Charbonneau et al. [40] assumption of possibility to modification of RVI in case of availability of two polarizations we assume an alternative to RVI for dual polarization as shown in Equation (7).
The index is about the contribution of volume scattering which is indicated by cross-polarized response. Pre-processing steps of satellite images were carried out using open source tools of Sentinel Application Platform (SNAP) software [42]. Mean backscatter values and temporal variation of backscatter for each field for three different crop types relying on backscatter statistic results are extracted using rasterstats zonal statistics and interpolated point queries function with a Python module [43]. The module is used to summarize geospatial raster datasets to extract information based on vector geometries [44]. Quantum GIS [45], an open source GIS software is applied to draw the region of interest (ROI) polygons as vector geometries used in rasterstats zonal statistics function. To decrease the effects of mix of the classes, the polygons are set at proper interval from the edges of the field boundary and homogeneous pixels are selected for the evaluation. Figure 7 shows the flowchart of Sentinel-1 dual polarimetric SAR data processing. Appl. Sci. 2019, 9, x FOR PEER REVIEW 11 of 19 Figure 7. Flowchart of Sentinel-1 dual polarization SAR data processing.

Results
Different Sentinel-1 dual polarimetric SAR parameters including linear polarization backscattering coefficients, H-A-α decompositions, polarization intensity ratios (VH/VV, VV/VH and VV − VH/VV + VH), Radar Vegetation Index (RVI), intensity arithmetic calculations and crop variables (crop height and CC) are investigated. Mean backscattering coefficient values for each sample site for different crops are defined and the correlation between the in situ measurements are analyzed. The results show that responses of polarimetric SAR parameters to the crop variables (crop height and CC) vary for different crop types at different phenological stages of the crops. For maize, SAR backscatter coefficient is more sensitive to crop height at the early stage of crop growth compared to the later stage, while the same situation is not true for sunflower and wheat.

SAR Parameters and Their Correlation with Crop Height
Almost in all visited fields, the mean backscatter values indicated similar tendency in four stages of crop growth. In Figure 8, mean backscatter coefficient (σVV and σVH and VV + VH) values of each sample site are determined and the correlation between the field measurements are presented. It is observed that the maize had high sensitivity at the beginning but starts to decrease when the maize height is higher than 150 cm at inflorescence emergence and heading stages. Sunflower is not sensitive to the crop height only when the sunflower height is greater than 90 cm. At the end of flowering and the inflorescence, it shows low sensitivity compared to its early and later stages. In comparison with maize and sunflower, wheat had relatively good sensitivity at the end of heading. Regarding wheat full inflorescence which is varying in different heights, the correlation between SAR

Results
Different Sentinel-1 dual polarimetric SAR parameters including linear polarization backscattering coefficients, H-A-α decompositions, polarization intensity ratios (VH/VV, VV/VH and VV − VH/VV + VH), Radar Vegetation Index (RVI), intensity arithmetic calculations and crop variables (crop height and CC) are investigated. Mean backscattering coefficient values for each sample site for different crops are defined and the correlation between the in situ measurements are analyzed. The results show that responses of polarimetric SAR parameters to the crop variables (crop height and CC) vary for different crop types at different phenological stages of the crops. For maize, SAR backscatter coefficient is more sensitive to crop height at the early stage of crop growth compared to the later stage, while the same situation is not true for sunflower and wheat.

SAR Parameters and Their Correlation with Crop Height
Almost in all visited fields, the mean backscatter values indicated similar tendency in four stages of crop growth. In Figure 8, mean backscatter coefficient (σ VV and σ VH and VV + VH) values of each sample site are determined and the correlation between the field measurements are presented. It is observed that the maize had high sensitivity at the beginning but starts to decrease when the maize height is higher than 150 cm at inflorescence emergence and heading stages. Sunflower is not sensitive to the crop height only when the sunflower height is greater than 90 cm. At the end of flowering and the inflorescence, it shows low sensitivity compared to its early and later stages. In comparison with maize and sunflower, wheat had relatively good sensitivity at the end of heading. Regarding wheat full inflorescence which is varying in different heights, the correlation between SAR backscattering (VV VH and σ 0 VV ) and wheat height is considerably less as in inferred from Figure 8c,f in growing stage. backscattering (VV VH and σ⁰VV) and wheat height is considerably less as in inferred from Figure 8 (c and f) in growing stage.  Table 3 gives the coefficient of determination (R 2 ) between Sentinel-1 dual polarimetric SAR parameters and measured crop height for maize, sunflower and wheat. Both σ⁰VV and σ⁰VH showed high correlation (R 2 = 0.81 and 0.80 respectively) and VV + VH had highest correlation (R 2 = 0.82) with maize height at the early stage. Whereas the correlation between SAR parameters to sunflower height is very low (R 2 = 0.31) at the same stage. Almost all SAR parameters show relatively good correlation at the early stage of wheat. For instance, a good negative correlation with VH/VV intensity ratio and σ⁰VH are obtained at the early stage of wheat (R 2 = 0.66, and 0.65 respectively). However, all of the crops have very weak correlation or even they are not correlated and sensitive to plant height in the later stage.   Table 3 gives the coefficient of determination (R 2 ) between Sentinel-1 dual polarimetric SAR parameters and measured crop height for maize, sunflower and wheat. Both σ 0 VV and σ 0 VH showed high correlation (R 2 = 0.81 and 0.80 respectively) and VV + VH had highest correlation (R 2 = 0.82) with maize height at the early stage. Whereas the correlation between SAR parameters to sunflower height is very low (R 2 = 0.31) at the same stage. Almost all SAR parameters show relatively good correlation at the early stage of wheat. For instance, a good negative correlation with VH/VV intensity ratio and σ 0 VH are obtained at the early stage of wheat (R 2 = 0.66, and 0.65 respectively). However, all of the crops have very weak correlation or even they are not correlated and sensitive to plant height in the later stage.
Amongst the H-Alpha decomposition parameters, Alpha (α) decomposition parameter represents the highest correlation (R 2 = 0.67) with wheat height during the early growth stage. The maize height is relatively correlated with H-Alpha decomposition parameters at the early stage, although no considerable correlation is observed during the early stage of sunflower and at the later stage of wheat. Figure 9 shows the correlation between the H-Alpha decomposition parameters and crops height of maize, sunflower and wheat. Amongst the H-Alpha decomposition parameters, Alpha (α) decomposition parameter represents the highest correlation (R 2 = 0.67) with wheat height during the early growth stage. The maize height is relatively correlated with H-Alpha decomposition parameters at the early stage, although no considerable correlation is observed during the early stage of sunflower and at the later stage of wheat. Figure 9 shows the correlation between the H-Alpha decomposition parameters and crops height of maize, sunflower and wheat. (g) (h) (i) Figure 9. Correlation between Entropy (a-c), Anisotropy (d-f) and Alpha (g-i) decompositions of maize, sunflower and wheat with their height during growing stages.

SAR Parameters and Their Correlation with Crop Coverage
The CC measured in field includes measuring the row and plant cover by steel tape in unit area. For validation of the CC calculated from measured data, photographs were taken using the camera in downward position and perpendicular to the ground with 100 cm distance from the camera lens. The CC extraction process is done by application of Python glob image processing package [46]. The photos are first converted to HSV (Hue, Saturation, and Value) model and used as a detector for the

SAR Parameters and Their Correlation with Crop Coverage
The CC measured in field includes measuring the row and plant cover by steel tape in unit area. For validation of the CC calculated from measured data, photographs were taken using the camera in downward position and perpendicular to the ground with 100 cm distance from the camera lens. The CC extraction process is done by application of Python glob image processing package [46]. The photos are first converted to HSV (Hue, Saturation, and Value) model and used as a detector for the type and shape and to do color constancy processing by grouping or classifying the image. Following the HSV conversion green mask is applied to slice the green areas as white ratio and black to the bare soil. The sample classified green area for wheat is given in Figure 10. The unit area from the photograph is calculated using the ground sampling distance (GSD) formula and setting calibration derived from images. The contributing parameters for determination of ground resolution are the camera's height above the ground, the camera's pixel size, and the lens' focal length.
(g) (h) (i) Figure 9. Correlation between Entropy (a-c), Anisotropy (d-f) and Alpha (g-i) decompositions of maize, sunflower and wheat with their height during growing stages.

SAR Parameters and Their Correlation with Crop Coverage
The CC measured in field includes measuring the row and plant cover by steel tape in unit area. For validation of the CC calculated from measured data, photographs were taken using the camera in downward position and perpendicular to the ground with 100 cm distance from the camera lens. The CC extraction process is done by application of Python glob image processing package [46]. The photos are first converted to HSV (Hue, Saturation, and Value) model and used as a detector for the type and shape and to do color constancy processing by grouping or classifying the image. Following the HSV conversion green mask is applied to slice the green areas as white ratio and black to the bare soil. The sample classified green area for wheat is given in Figure 10. The unit area from the photograph is calculated using the ground sampling distance (GSD) formula and setting calibration derived from images. The contributing parameters for determination of ground resolution are the camera's height above the ground, the camera's pixel size, and the lens' focal length. The coefficient of determination (R 2 ) between Sentinel-1 polarimetric SAR parameters and the measured CC for maize, sunflower and wheat for early and later growing stages are given in Table  4. According to the principle growth stage [20] when the maize reaches to the end of its stem elongation, leaves completely unfold and have full size. When the BBCH-scale is greater than 39 for maize, CC is measured 75%. The CC for sunflower and wheat is measured as 85% and 60%, respectively when the BBCH-scale is greater than 51 and 49.  The coefficient of determination (R 2 ) between Sentinel-1 polarimetric SAR parameters and the measured CC for maize, sunflower and wheat for early and later growing stages are given in Table 4. According to the principle growth stage [20] when the maize reaches to the end of its stem elongation, leaves completely unfold and have full size. When the BBCH-scale is greater than 39 for maize, CC is measured 75%. The CC for sunflower and wheat is measured as 85% and 60%, respectively when the BBCH-scale is greater than 51 and 49.
Considering Sentinel-1 dual polarimetric SAR parameters and measured CC, the highest correlation is obtained for the VV polarization (R 2 = 0.73) and VV+VH (R 2 = 0.73) when the CC is lower than 75% for maize while with increasing CC the correlation decreased (R 2 = 0.49) for both crops. This is contrary to the findings of Liao et al. [17] where the sensitivity of RADARSAT-2 polarimetric SAR and its correlation with fractional vegetation cover (FVC) in HV polarization is high for maize, and they suggest that VV polarization is not a useful parameter for monitoring broad-leaf crops. Their findings may have conflict with our results due to the sensors properties. Regarding the sunflower, in its early stage (i.e., CC < 85%), again the higher coefficient of determinations are obtained for VV polarization (R 2 = 0.46) and VV + VH (R 2 = 0.47). However, comparing the maize and the sunflower as two different broad-leaf crops in our study, R 2 of the sunflower is lower than the maize. The discrepancy in correlations might be due to the difference in leaves geometry of maize and sunflower. In other respects, similar to the findings of Liao et al. [17], no correlation is observed for wheat at both stages in our study. This may be due to the wheat structure and leaves geometry as narrow-leaf crop. Wheat can reach its full development at early stage and penetrating from biomass occurs due to narrow leaves, stem affectation and contributing of underlying soil resulting attenuation of SAR backscatter [47]. Figure 11 depicts the parameters which have highest correlation with CC during the growth stage for maize and sunflower.

Conclusions
In this study, three different crop types that have "broad" and "narrow" leaves were selected to investigate the sensitivity of the Sentinel-1 dual polarimetric SAR parameters to plant height and CC. Maize and sunflower is selected as broad leaf and wheat as narrow leaf crop. The reason for studying two broad-leave crops was to validate the high sensitivity of SAR polarimetric data to early stages of broad-leave crops and then comparing the findings for wheat as narrow-leave crop. We demonstrate that maize presents higher correlation during the early stages of the crop growth when the crop height is less than 150 cm. It is strongly correlated with the SAR parameters including VV + VH (R 2 = 0.82), VV (R 2 = 0.81), and VH (R 2 = 0.80). Besides, CC of maize was well correlated with VV polarization (R 2 = 0.73) at the early stage before the heading stage although at the later growing stage the correlation becomes weaker after the heading stage. From the backscatter analysis, the same result is not observed in the sunflower. The sunflower height has very low correlation with the most of SAR parameters. Only VH polarization shows slightly better sensitivity when its height is below 92 cm before the ending of flowering and the inflorescence reaches to full size. The sunflower CC is relatively correlated with VV polarization at the early stage (during the flowering stage) while any considerable correlation between SAR parameters and sunflower height and CC is observed at the later stage. The sensitivity of SAR parameters to wheat variables is often low compared to maize and sunflower. The high, but negative correlations are related to the VH/VV intensity ratio and VH polarization. However, Alpha (α) decomposition parameter shows highest correlation at the beginning stage and represents absence of SAR parameters sensitivity with wheat height at the later growing stage and CC at both stages. The results we have obtained reveals that Sentinel-1 dual polarimetric SAR (C-band) has a high potential for identifying growth stages and estimation of crop height, canopy coverage of maize as a broad-leaf crops. However, this is not proven for sunflower,

Conclusions
In this study, three different crop types that have "broad" and "narrow" leaves were selected to investigate the sensitivity of the Sentinel-1 dual polarimetric SAR parameters to plant height and CC. Maize and sunflower is selected as broad leaf and wheat as narrow leaf crop. The reason for studying two broad-leave crops was to validate the high sensitivity of SAR polarimetric data to early stages of broad-leave crops and then comparing the findings for wheat as narrow-leave crop. We demonstrate that maize presents higher correlation during the early stages of the crop growth when the crop height is less than 150 cm. It is strongly correlated with the SAR parameters including VV + VH (R 2 = 0.82), VV (R 2 = 0.81), and VH (R 2 = 0.80). Besides, CC of maize was well correlated with VV polarization (R 2 = 0.73) at the early stage before the heading stage although at the later growing stage the correlation becomes weaker after the heading stage. From the backscatter analysis, the same result is not observed in the sunflower. The sunflower height has very low correlation with the most of SAR parameters. Only VH polarization shows slightly better sensitivity when its height is below 92 cm before the ending of flowering and the inflorescence reaches to full size. The sunflower CC is relatively correlated with VV polarization at the early stage (during the flowering stage) while any considerable correlation between SAR parameters and sunflower height and CC is observed at the later stage. The sensitivity of SAR parameters to wheat variables is often low compared to maize and sunflower. The high, but negative correlations are related to the VH/VV intensity ratio and VH polarization. However, Alpha (α) decomposition parameter shows highest correlation at the beginning stage and represents absence of SAR parameters sensitivity with wheat height at the later growing stage and CC at both stages. The results we have obtained reveals that Sentinel-1 dual polarimetric SAR (C-band) has a high potential for identifying growth stages and estimation of crop height, canopy coverage of maize as a broad-leaf crops. However, this is not proven for sunflower, that may be due to the crop structure and leaves geometry since they may change the SAR backscatter value in any stages of crops. This study also demonstrates that Sentinel-1 dual polarimetric SAR data can be a good alternative to other commercial data which enables users to access freely available of a constant long-term data archive for applications requiring long-range time series.
Although the use of single orbit direction (descending) provides possibility of monitoring and investigating agricultural growth stages, future studies should use both the ascending and descending orbits to have independent results from different angles of view. In addition, the relationship between backscatter values of multi-temporal Sentinel-1 data and vegetation variables will be incorporated to improve the crops mapping and classification accuracy as a future work.