Evaluation of Different Radiative Transfer Models for Microwave Backscatter Estimation of Wheat Fields

: This study aimed to analyze existing microwave surface (Oh, Dubois, Water Cloud Model “WCM”, Integral Equation Model “IEM”) and canopy (Water Cloud Model “WCM”, Single Scattering Radiative Transfer “SSRT”) Radiative Transfer (RT) models and assess advantages and disadvantages of different model combinations in terms of VV polarized radar backscatter simulation of wheat ﬁelds. The models are driven with ﬁeld measurements acquired in 2017 at a test site near Munich, As vegetation descriptor for the canopy models Leaf Area Index (LAI) was used. The effect of empirical model parameters is evaluated in two different ways: (a) empirical model parameters are set as static throughout the whole time series of one growing season and (b) empirical model parameters describing the backscatter attenuation by the canopy are treated as non-static in time. The model results are compared to a dense Sentinel-1 C-band time series with observations every 1.5 days. The utilized Sentinel-1 time series comprises images acquired with different satellite acquisition geometries (different incidence and azimuth angles), which allows us to evaluate the model performance for different acquisition geometries. Results show that total LAI as vegetation descriptor in combination with static empirical parameters ﬁt Sentinel-1 radar backscatter of wheat ﬁelds only sufﬁcient within the ﬁrst half of the vegetation period. With the saturation of LAI and/or canopy height of the wheat ﬁelds, the observed increase in Sentinel-1 radar backscatter cannot be modeled. Probable cause are effects of changes within the grains (both structure and water content per leaf area) and their inﬂuence on the backscatter. However, model results with LAI and non-static empirical parameters ﬁt the Sentinel-1 data well for the entire vegetation period. Limitations regarding different satellite acquisition geometries become apparent for the second half of the vegetation period. The observed overall increase in backscatter can be modeled, but a trend mismatch between modeled and observed backscatter values of adjacent time points with different acquisition geometries is observed. coupled performance of and models how of backscatter attenuation through the canopy inﬂuence the total backscatter calculation for different vegetation stages. two novelties of this are the evaluation of different combinations of widely used surface and canopy RT models on one test site and the analysis over time of empirical model parameters coe f (SSRT) and B (WCM) describing the backscatter attenuation through the canopy T By using a backscatter static calibration performance results in the form non-static The validation accuracy for the non-static parameter approach was evaluated with ubRMSE and ranges for all model combinations between 1.82 and 2.22dB. Validation results with SSRT as canopy model show better in combination with all surface models when compared to respective combinations using WCM canopy has shown backscatter

• investigate different RT model combinations and asses their advantages and disadvantages; • evaluate if different radar acquisition geometries are modeled adequately with the used RT models; and • serve as preliminary work for future synergistic retrieval approaches of SAR and optical sensors with a focus on high spatial and temporal resolutions.
Section 2 presents the used dataset. Section 3 summarizes the used RT models. In Section 4, calibration and validation results are shown and discussed. Finally, the main conclusions are drawn in Section 5.

Study Area
The study area, Munich North Isar (MNI), is located in southern Germany (Bavaria), near Munich (48 • 13'N-48 • 20'N, 11 • 39'E-11 • 45'E, Figure 1). Since 2014, almost every year from spring until autumn, different field campaigns targeting agricultural purposes were carried out [65][66][67][68]. From March to September 2017, an intensive field campaign focusing on maize and wheat fields for validation of soil and vegetation parameter retrievals from Sentinel-1, Sentinel-2 and the future EnMAP satellite was conducted. MNI is characterized by intensive agriculture with wheat, maize, and grassland as main crop types. In close vicinity (<10 km) to the test site two meteorological stations, Freising (470 m a.s.l.) and Eichenried (475 m a.s.l.) managed by the Bavarian State Research Institute (LFL) and one meteorological station Munich-airport (446 m a.s.l.) managed by the German Meteorological Service (DWD) are situated. The annual measured mean temperature for 2017 ranges between 9 • C (Freising) and 9.3 • C (Eichenried). The average annual precipitation for 2017 reached 753 mm (Munich-airport) to 853 mm (Eichenried). The used data for this study includes field campaign data of wheat fields from 2017 (Section 2.2) and Sentinel-1 satellite data (Section 2.3).

Field Data
During the MNI field campaign of 2017, weekly field measurements of different biophysical parameters (Table 1) were conducted. The total LAI was measured with a LI-COR Biosciences LAI-2200C device (LI-COR Biosciences Inc., Lincoln, NE, USA) as an average of 14 measurements from the same area. The measurements were taken within each test field at three different locations ( Figure 1). The accuracy of LAI in terms of mean standard deviation of repeated measurements ranges within the fields between 0.45 and 0.52. The monitoring period started end of March and ended shortly before the fields were almost simultaneously harvested in mid of July. Additionally, Decagon TM5 soil moisture sensors using the capacity method were installed permanently within the first five centimeters of the soil surface. Soil moisture changes were monitored with a time interval of 10 min. Information about the soil was provided by earlier campaigns when soil samples were taken from the fields, and the soil properties were analyzed in the laboratory ( Table 2). The soils bulk density with 1.45 ± 0.13 g/cm 3 , and the clay content with 7.38 ± 1.8%, show no high variability between the different fields. The sand content, on the other hand, shows higher variability, with 24.08 ± 10.46%.

Satellite Data
For this study, C-band Sentinel-1 SLC data of Sentinel-1A/B was used. The Sentinel data was pre-processed with ESA's SNAP Toolbox Version 7.0.3. An overview of all applied pre-processing steps is given in Figure 2. For the geometric correction, SRTM data with 1 arc-second resolution was chosen as digital elevation model input. Afterward, the radiometric correction method of Kellndorfer et al. [69] was applied. In a second pre-processing step, a multi-temporal Lee-sigma filter was used for speckle reduction. The temporal filter was applied on each image with information of 6 other images (three before the target and three after) with a spatial window size of 5 × 5 pixels, a sigma of 0.9, and a target window size of 3 × 3 pixels. For the period of the field campaign in 2017 (March to July), in total, 78 Sentinel images covering the study area are available. Considering images with different orbit directions (ascending and descending) and different incidence angles (ranges from 35 • to 45 • ) a revisit time of 1.5 days was archived. The spatial resolution of the processed data was 10 × 10 m. A more detailed overview of the used Sentinel-1 dataset and the image properties is given in Table 3. The primary acquisition mode of Sentinel-1 provides data with polarization VV and VH. For our study, the focus was set on polarization VV due to findings that, for retrieving soil moisture, the usage of VH alone or in addition to VV is not suitable for well-developed vegetation [45,70].  Schematic overview of SNAP pre-processing steps to retrieve geometric and radiometric corrected images from Sentinel-1 SLC data (left). Pre-processing steps for speckle reduction using a multi-temporal speckle filter (right).

Empirical Water Cloud Model (WCM Surface Part)
The WCM, often referred to as the tau-omega model, was developed by Attema and Ulaby in 1978 [22]. For a given polarization pq (pq = HH, VV, or HV), the surface contribution σ 0 s pq of WCM to the backscattered radar signal in dB scale is defined as with empirical fitted soil parameters C pq and D pq and soil moisture content mv. C pq is an empirical calibration constant, whereas D pq , as a calibration factor, indicates the sensitivity of soil moisture on the received radar signal. The WCM surface part is a purely empirical model thus no additional information about surface roughness or incidence angle is needed. Therefore, the empirical parameters have to be calibrated for each test site separately.

Semi-Empirical Oh Model 1992 (Oh92)
In 1992, Oh et al. [26] developed an approach for the retrieval of soil moisture and soil surface roughness by empirical determined functions. Based on scatterometer measurements and various ground measurements of soil moisture and surface roughness, two functions for the co-(p = and with θ as local incidence angle, k as radar wave number (k = 2π/λ), where λ is the wavelength, and s as rms height. R 0 is the Fresnel reflectivity coefficient at nadir given by where r is the relative dielectric constant. The VV polarized backscatter coefficient σ 0 s VV is further defined as with the Fresnel coefficients for horizontal R h and vertical R v polarization where µ r is the relative permittivity. Furthermore, the backscatter coefficients σ 0 s HH and σ 0 s HV are given with respect to σ 0 s VV , p, and q by The model in its original form can be applied for the retrieval of soil moisture or soil surface roughness for bare soil conditions at several frequencies (X-to L-Band) and a broad range of incidence angles (10-70 • ). Because it is a semi-empirical model, the validity range of the model in terms of soil moisture and soil surface roughness is defined with 0.1 < ks < 6, and 9 Vol.% < mv < 31 Vol.%.

Semi-Empirical Oh Model 2004 (Oh04)
In 2004, Oh [27] revised and simplified his original approach to use only soil moisture (mv) as an independent variable rather than R 0 and r (Section 3.1.2). Thus, if using mv as input variable, no additional information about the soil properties (bulk density, sand and clay content) is needed. The model is defined by Oh04 is optimized for bare soils with 0.13 < ks < 6.98, 4 Vol.% < mv < 29.1 Vol.% and 10 • < θ < 70 • .
The Dubois model was optimized for bare soil conditions and has a validity range for soil moisture of mv ≤ 35 Vol.% and soil surface roughness of ks ≤ 2.5.

Physical Integral Equation Model (IEM)
The IEM is a theoretical backscattering model and was developed by Fung et al. [29] in 1992. Since then, Fung and colleagues extended the IEM to bistatic scattering [71]. The general co-polarized backscatter coefficient σ 0 pp for pp = VV or HH is defined as where I n pp is I n pp = (2 ks cosθ) n f pp e −k 2 s 2 cos 2 θ + (ks cosθ) n F pp , with W n as the Fourier transform of the nth power of the surface correlation function p(x, y). Furthermore, the backscatter at cross polarization σ 0 s HV is calculated as The Kirchhoff coefficients f HH , f VV and complementary field coefficients F HH , F VV , F HV are given as with Fresnel coefficient at horizontal R h (Equation (6)) and vertical R v (Equation (7)) polarization thus R is described by The Fourier transform of the nth power of the surface correlation coefficient W n (a, b) is calculated by The distribution of surface correlation function p(x, y) can be described for low surface roughness as exponential and high surface roughness values as Gaussian by p(x, y) = e −( |x|+|y| L ) (exponential), (25) p(x, y) = e −( with L as correlation length.

Empirical Water Cloud Model (WCM)
The WCM [22] with respect to surface σ 0 s pq and canopy σ 0 c pq contribution, as well as two-way attenuation T 2 pq , is defined as whereas the canopy part σ 0 c pq in linear scale and the two-way attenuation T 2 pq are written as where θ is the local incidence angle, V1 and V2 are empirical vegetation descriptors, and A pq , as well as B pq , are fitted parameters of the model which depend on the vegetation properties and the radar configuration. For σ 0 s pq in Equation (27), each of the described surface models in Section 3.1 can be used.

Semi-Empirical Single Scattering Radiative Transfer (SSRT) Model
The SSRT model used by De Roo [35] and Ulaby [17] is a semi-empirical first-order scattering model. The model is defined as where with T p and T q as attenuation of the canopy for different polarizations and σ 0 s pq describing the pure surface scattering mechanism. Similar to the definition of WCM in Section 3.2.1, all surface models described in Section 3.1 can be used for calculating the surface contribution σ 0 s pq within SSRT. Furthermore, the p polarized one way transmittivity of the canopy T p is defined as with τ p as the p polarized attenuation of the canopy given by whereas H represents the canopy height. The extinction coefficient k p e which accounts for the absorption and scattering losses of the electromagnetic wave through the canopy is defined as In general, a canopy consists of leaves, stalks, and branches with different shapes and orientations, which do not show a linear distribution in the vertical. However, in the applied SSRT, it is assumed that k p e , k p a , and k p s follow a uniform distribution in the vertical as a function of z within the canopy layer. In addition to the extinction coefficient, the scattering part k p s of k p e can be derived by where ω represents the single scattering albedo. For the direct backscattering contribution of the canopy σ 0 c pq , Attema and Ulaby's [22] water cloud approach of identical scatterers, which are uniformly distributed within the volume, is used. Thus, multiple scattering effects are ignored. As a consequence, the volume backscattering coefficient σ back V pq of the vegetation medium is defined as with N v as the number of scattering particles per unit volume and σ back pq as the pq polarized backscattering cross section of a single particle. Finally, the pq polarized canopy backscattering coefficient σ 0 c pq within Equation (30) can be obtained from Furthermore, ground/canopy (σ 0 gc pq ) and canopy/ground (σ 0 cg pq ) scattering contributions are defined as where H is the canopy height, σ bist v pq is the bi-static scattering cross section of a single leaf or stalk, and R p describes the p polarized Fresnel reflectivity (Equations (6) and (7)). Thus, the total canopy ground contribution σ 0 cgt pq within Equation (30) as the sum of σ 0 gc pq and σ 0 cg pq can be written as Furthermore, ground canopy ground contribution (σ 0 gcg pq ) within Equation (30) is defined as

Practical Considerations
Each described model requires a different set of input parameters. A summary for the different RT models is given in Table 4. For some parameters, field measurements (Section 2.2) or literature values (s and ω, Table 5) are used, and other site dependent parameters have to be fitted. The analyzed wheat fields were sown in autumn of 2016. By the starting point of the observation period (end of March 2017), the soil surface was already smoothed out by rain and the ground was covered by wheat plants (height > 10 cm). Marzahn et al. [72] showed that, for wheat fields with the above mentioned state, only minor changes in terms of roughness throughout the vegetation period are observable. Previous studies regarding periodic features and roughness changes [73,74] found out that changes in surface roughness due to soil rows as periodic feature within wheat fields are essential if the viewing angle is nearly perpendicular to the row orientation, but, for other viewing angles, the changes within surface roughness are negligible. The viewing angle and the row orientation of the wheat fields in our study are always <75 • . Therefore, changes due to periodic soil rows are assumed to be negligible. Typical roughness measurements of various winter wheat fields suggest rms values between 1.0 and 1.3 [72,73,[75][76][77][78]. With the assumption of only minor roughness changes throughout the vegetation period, a literature value for surface height s of 1.2 cm was chosen. The value of 1.2 cm was used due to former field campaign measurements of wheat fields in Germany [72,73]. For single scattering albedo ω, a common literature value of 0.03 [79] was set. The conversion from soil moisture field measurements to the required model input of dielectric constant for models Oh92, Dubois95, and IEM was performed by using a dielectric mixing model for soils after Dobson et al. [80]. The required soil information about sand and clay content, bulk density (Table 2), and soil moisture (Table 1) were provided by field measurements and laboratory results. Additionally, for some parameters, in their original form, adjustments were made within this study. In particular, DeRoo et al. [35] used for parameterizing of the extinction coefficient k p e a combination of an empirical parameter, vegetation water mass and vegetation height. In our study, the vegetation water mass, and the vegetation height were replaced by total LAI. Therefore, k p e is defined by with coef as an empirical parameter. To reduce the required parameters for model IEM, a well-established approach of Baghdadi et al. [31,45] was used. The correlation length L was replaced by a fitted parameter Lopt which is dependent on s, θ, and the polarization. Lopt for C-band VV polarization data and the Gaussian correlation function after Baghdadi et al. [45] is defined as The refined version of Baghdadi is hereinafter referred to as IEM_B. A schematic illustration of the RT model calibration approach is shown in Figure 3. All RT model combinations (surface + canopy) are driven by field measurements and the required empirical parameters ( Table 4). The fitting of the empirical parameters was carried out by minimizing the sum of the squared error between modeled and measured (Sentinel-1) radar backscatter values. For the measured Sentinel-1 backscatter value of each measurement point shown in Figure 1, the mean backscatter of 5 × 5 pixel (50 × 50 m) around the measurement location was chosen. In a first fitting approach, all parameters shown in Table 4 were defined to be static for the entire vegetation period. In a second fitting approach, WCM parameters C, D, and A were set to be static; therefore, the mean (Table 5) of all model results of the static approach was used. The used values for the different parameters for validation are shown in Table 5. Additionally, the attenuation of the backscatter through the canopy was defined to be variable throughout the time series. More specifically, parameters coe f within SSRT and B within WCM were fitted for each time step individually by taking three observations before and after into account. By making only coe f or B variable, changes within the results can then clearly related to changes of the attenuation of the radar backscatter signal by the canopy. Field measurements used as model input parameters show multidimensionally unstructured inter-and intra-field correlations. Therefore, measurement point independency is assumed, and a validation of the parametrized RT models is performed using a leave-one-out cross-validation approach. Hereby, the parameter mean of the calibration results of eight measurement points is validated with the remaining measurement point.

Differences Between Applied Models
WCM is a purely empirical model and depends, therefore, only on the calibrated empirical coefficients. It is applicable under almost all surface/vegetation conditions, but it has to be calibrated for different test sites separately. The transferability to other test sites and or other surface/vegetation conditions is not possible. Surface models, like Oh92, Oh04, Dubois95, and IEM, were originally only developed for bare soil and or sparse vegetation conditions. Nevertheless, different studies replaced the surface component within WCM and SSRT with Oh [17,35,43,47], Dubois95 [51,54,55], or IEM/IEM_B [17,45,53] models. One key advantage of the semi-empirical-based (Oh92, Oh04, Dubois95) or theoretical-based (IEM/IEM_B) surface models in comparison to WCM is their better transferability to other test sites and surface/vegetation conditions. Models, like Oh92, Oh04, or Dubois95, are based on a hybrid construction with experimental data guided by trends predicted by theoretical models [17]. Theoretical models, like IEM/IEM_B, on the other hand, have a theoretical foundation, whereas, for used mathematical approximations within the model, various assumptions to retrieve a analytical solution are made [17]. The Oh model of 1992 was developed based on a single experiment with information about only four different soil surfaces [26]. For the model version of 2004, Oh used information of approximately 40 bare soil fields conducted over seven experiments [27]. Furthermore, the usage of and R in model Oh92 was refined to the usage of mv in model Oh04. Model Dubois95 is the only model which was developed only for co-polarized backscatter data (HH or VV), whereas the other models can calculate co-and cross-polarized backscatter values (HH, VH, VV). Differences between vegetation models WCM and SSRT exists in the form of their simplicity. WCM calculates only the volume backscattering component, whereas, within SSRT, additional backscatter components, such as plant-ground and ground-plant interaction, as well as ground-plant-ground scattering contributions, are considered. In our implementation, the vegetation descriptor of WCM is LAI, whereas the vegetation descriptor of SSRT consists of LAI and canopy height. In general, the computational time and the required additional input parameters are increasing from empirical to semi-empirical and theoretical models. A summary of the different models with information about type, validity range, site dependency, required parameters, and used polarization are given in Table 4.

Static Empirical Parameters
In a first calibration approach, the empirical parameters for the different models were treated as static throughout the entire time series. The modeled backscatter was then compared to the measured backscatter from Sentinel-1. Table 6 shows the mean of RMSE and R 2 of all analyzed sample points ( Figure 1) for different surface and canopy model combinations. The retrieved RMSE of the calibration results ranges from 1.92 to 2.25 dB with R 2 of 0.08 to 0.34, respectively. A more detailed picture regarding differences between modeled and Sentinel-1 backscatter during the time series is shown for field 508 in Figure 4. While, for the first half of the vegetation period, all model combinations show a relatively good fit to Sentinel-1 backscatter data, in the second half, deviations are obvious. Furthermore, from the beginning of June where LAI reaches its saturation point and the maximum plant height is almost reached, no significant change over time within the modeled backscatter can be observed. The analyzed soil models (different colors) show small differences, whereas a clear separation between analyzed canopy models (solid vs dashed lines), especially for later vegetation stages, is noticeable. As described in Section 2.3, Sentinel-1 data of four different overpasses, and therefore with different incidence and azimuth angles, were used for this analysis. In Figure 4, every fourth point (same icon) of the Sentinel-1 backscatter time series represents the same satellite acquisition geometry (same incidence and azimuth angle). The incidence angle is implemented within the used RT models, whereas the models do not account for difference azimuth angles. The black line (Sentinel-1 backscatter) in Figure 4 shows that the observed backscatter values differ with varying incidence angles and changes in soil moisture. The model predictions (different colors) illustrates that the models can account for varying incidence angles and changes in soil moisture only until the end of May. Furthermore, the good correlation of modeled and Sentinel-1 backscatter values suggests that the effect of different azimuth angles on backscatter values seems to be negligible until the end of May. The main change within the wheat fields in June and July in terms of phenology are the flowering, the development of the fruit, and, later on, the ripening [75,81]. For these phenology stages, the increase in backscatter is caused by higher sensitivity of the radar signal to the ground contribution due to water loss within the vegetation [81]. Mattia et al. [76] identified the heading period as turning point where the sensitivity of the radar backscatter to above-ground biomass decreases, whereas the sensitivity to soil surface increases. The temporal evolution of modeled ground contribution to the total backscatter for different model combinations (different colors) with the observed Sentinel-1 total backscatter (black line) is shown in Figure 5 (top part). Until the increase in canopy height at the beginning of May (Table 4), the modeled ground scattering part seems to be the main contributor to the total backscatter, whereas the canopy part is negligible. With increase, especially in canopy height, the ground contribution drops significantly. The expected decrease in backscatter of the ground contribution due to a bigger canopy layer (increase of canopy height and LAI) can be more clearly observed for SSRT than WCM. The differences between SSRT and WCM might be found by using different canopy descriptors (SSRT: LAI and canopy height; WCM: LAI). Differences in ground contributions between Dubois95 and the other surface models are related to differences in the modeled attenuation through the canopy T (Figure 4, bottom part). T is regulating the contribution intensity of the ground and canopy part for the total backscatter calculation. The temporal evolution of T with a value range from 0 (dominant canopy contribution) to 1 (dominant ground contribution) is shown in Figure 5 (bottom part). All model combinations show a similar temporal shape with slightly higher values for model Dubois95. T decreases from April (dominant ground contribution) to mid-May (dominant canopy contribution) and stays at its minimum after mid-May. The expected increase [76,81] of the modeled ground contribution due to higher ground sensitivity in June and July (phenology: flowering, development of the fruit, and ripening) is not observed within the modeled data. In our case, T is mainly driven by static empirical parameter B (WCM) and coe f (SSRT), as well as non-static vegetation descriptors LAI and/or canopy height. Therefore, with almost no changes in LAI and canopy height in June and July (Figure 4, middle part), the two-way attenuation T stays near zero, which indicates a dominant canopy model contribution to the total backscatter calculation. By definition, LAI is defined as the one-sided leaf area per measured ground unit [82]. By the time wheat plants reach their maximum height, the leaves are fully developed. Changes within the wheat plants, especially during vegetation stages of flowering, fruit development, or ripening (Figure 4), are based mainly on increasing biomass within grains and stems, as well as changes of the vegetation water content. However, biomass changes in grains and stems, as well as vegetation water content loss, especially during the ripening stage, are not reflected within the LAI. Therefore, almost no information about the increased biomass and the water loss due to ripening of the plants is given within this model configuration. Plant moisture reduction affects the attenuation of the radar signal by the canopy in a way that the canopy is more transparent for the radar wave [83]. Therefore, the sensitivity of the radar signal to the canopy should decrease, whereas the sensitivity to the surface increases. The almost non-existent deviation between Sentinel-1 and modeled backscatter in early vegetation stages suggests that the interaction between surface and canopy model, and therefore the attenuation of the backscatter signal by the canopy, described by static empirical parameters and LAI, can be modeled sufficiently only in early vegetation stages. For good backscatter model results during later vegetation stages, the backscatter changes due to water loss within the plants have to be considered. The effects of utilizing non-static empirical parameters to account for these shortcomings are discussed in the next sections.

Non-Static Empirical Parameters
Due to the results in Section 4.1.1, which showed that LAI and height cannot account for the observed backscatter changes in June and July (increase of Sentinel-1 backscatter), a second calibration approach with non-static empirical parameter where chosen. As already mentioned, the observed increase in backscatter is caused by a higher sensitivity of the radar signal to the ground contribution. The attenuation of the backscatter T, more precisely B, for WCM and coe f for SSRT were identified as the main drivers for the increase or decrease of the ground contribution; therefore, a non-static approach for these parameters was tried. The other empirical parameters of WCM (surface: C, D; canopy: A) were set as static values (mean of the retrieved static values during the calibration approach shown in Section 4.1.1) to clearly relate the observed changes to the attenuation of the backscatter. The modeled backscatter was compared to Sentinel-1 observations and the statistics for different model combinations are shown in form of mean RMSE and R 2 of all analyzed field measurement points in Table 7. The retrieved RMSE ranges from 1.13 to 1.60dB with R 2 of 0.45 to 0.82, respectively. By comparing different model combinations, almost no differences between the two analyzed canopy models can be observed. A different picture is shown by comparing surface models only. WCM seems to outperform all others, whereas differences between WCM and Oh92, Oh04, and IEM_B are smaller than differences of these models to Dubois95. Similar to the approach with static empirical parameters, the evolution over time of modeled and Sentinel-1 backscatter is shown in Figure 6 for one measurement point of field 508. Like the results in Section 4.1.1, modeled backscatter results of the first half of the vegetation period fit well to the observed Sentinel-1 backscatter data. Contrary to the static parameter approach, the second half of the vegetation period shows high correlations between Sentinel-1 and modeled backscatter. Unlike the results of the static approach ( Figure 5, top part), the expected increase of the ground contribution at the end of the vegetation period can be observed in Figure 7 (top part). The increase of the ground contribution to the total backscatter is also reflected within the two-way attenuation by the canopy T shown in Figure 7 (bottom part). Compared to the static approach ( Figure 5, bottom part), the values of T from April to the beginning of June are very similar for all model combinations, except Dubois95. The temporal changes of non-static parameters B and coe f are shown for the validation results in Figures 8 and 9 and will be further discussed in the validation Section 4.2. Besides the modeled backscatter increase for the second half of the vegetation period, different trends between Sentinel-1 and modeled backscatter for individual time steps are observed. This mismatch might occur for several reasons. As already mentioned in Section 4.1.1, Sentinel-1 data with different incidence and azimuth angles were used. A closer look at every fourth modeled point (same incidence and azimuth angle) in Figure 6 (same icon) shows a steady increase of modeled backscatter values at the end of the vegetation period. Therefore, an overall trend of an increasing backscatter can be seen for Sentinel-1 and modeled data. Due to non contradictory trends during the first half of the vegetation period, it seems that the influence on backscatter due to incidence and azimuth angles increases for the second half of the vegetation period. Thus, higher canopy heights and LAI values may also increase the impact of different incidence and azimuth angles on the backscatter behavior. In addition, the acquisition time (Table 3) of the Sentinel-1 images might play another role for different attenuation effects by the canopy [81]. The acquisition time of the satellite differs for different overpasses. For the MNI test site, Sentinel-1 data was acquired during early morning or late afternoon (Table 3). This might lead to differences in the observations due to dew [85,86] or different plant alignments towards the sun which, in theory, leads to different backscatter attenuations through the canopy [87][88][89], which are not accounted for within the models. Despite the different trends during the temporal evolution of modeled and Sentinel-1 backscatter for the second half of the vegetation period, the observed overall increase in backscatter at the end of the vegetation period can be modeled well.

Model Validation Results
The validation for the non-static approach was performed by a leave-one-out cross-validation method. Means of RMSE, ubRMSE, and R 2 are shown in Table 7. The results for the non-static approach show ubRMSE values between 1.92 and 2.22dB. All validation (ubRMSE) results are in a range of 0.65 to 0.95dB poorer than the calibration (RMSE) results. Furthermore, for the validation results R 2 ranges between 0.49 and 0.64, whereas the poorest results are archived with surface model Dubois95. A comparison of ubRMSE for canopy model SSRT and WCM yields for all surface models slightly better results for canopy model SSRT. A more detailed overview of the validation results separated by models and measurement points is given in Figure 10 for canopy model SSRT and in Figure 11 for canopy model WCM. Comparing different fields, as well as field measurement points, a diverse picture is drawn. In general, field 508, especially 508-1, shows the best results, with R 2 values mainly higher than 0.6 and ubRMSE lower than 2.2. Field 542, and especially 542-1, show the worst results, with R 2 around 0.4 and ubRMSE higher than 2.2. Comparing different surface model combinations, IEM_B and WCM show, for both canopy models, slightly better results than the others. It is noticeable that the results of Dubois95 are, in general, poorer than the results of the other surface models, with the exception of point 542-1 (measurement point with the poorest overall results) and partly 301-2. Differences between Oh92 and Oh04 for the different measurement points are present but, compared to the other models, very low. A closer look into each model combination, and especially at the used non-static empirical parameters coe f and B, and thus at the change of the attenuation of the canopy in time, are given with show, for all surface models (except Dubois95), for parameter coe f , a almost steady decrease from the beginning to the end of the vegetation period and a similar evolution for parameter B, except for a short period of increasing values to a relative maximum from mid-May to the end of May. The differences in shape between coe f and B might be explained due to different model input data. SSRT uses the canopy height and LAI for the description of the canopy, whereas WCM only uses LAI. The increase of parameter B correlates very well with the measured increase of the canopy height. This suggests that parameter B compensates for possible shortcomings of WCM due to missing information about the canopy height. Furthermore, not only the absolute values but also the spread of coe f and B for different models are higher at the beginning of the vegetation period and strongly decrease over time. This means that coe f and B show higher differences between fields and field measurement points at the beginning of the vegetation period. This is in line with observations made during the field campaign, where higher differences within one field and between fields could be observed more easily at early vegetation stages. At the end of the vegetation period, such differences were not detectable anymore. Another indicator of low field differences regarding the ripening stage of the wheat plants was the almost simultaneous harvesting date. Another factor to be considered if looking at higher variability of coe f and B at the beginning of the vegetation season is the model implementation itself. B and coe f are influencing the attenuation of the backscatter by the canopy and therefore how strong each model compartment (surface or canopy) contributes to the total backscatter calculation. At the beginning of the growing season, LAI and canopy height are low; therefore, notable differences within coe f and B might have not such a high impact on total backscatter predictions. Differences in the range of coe f (0 to 2.5) and B (from 0 to 0.35) are most likely based on differences between the model definitions and the required input parameters.

Conclusions
Modeled backscatter results of wheat field time series data using different RT model combinations (surface: Oh92, Oh04, Dubois95, IEM_B, WCM; canopy: WCM, SSRT) were compared to observed C-band data from Sentinel-1. Differences between the models were analyzed. The used dataset was acquired by an intense field campaign throughout one vegetation period in 2017. The analysis focused on coupled performance of surface and canopy models and especially on how changes of backscatter attenuation through the canopy influence the total backscatter calculation for different vegetation stages. The two novelties of this study are the evaluation of different combinations of widely used surface and canopy RT models on one test site and the analysis over time of empirical model parameters coe f (SSRT) and B (WCM) describing the backscatter attenuation through the canopy T.
Results show that, for total LAI, as vegetation descriptor, a static parameter influencing the backscatter attenuation through the canopy is suitable for the first half but not for the second half of the vegetation period. By using a non-static parameter approach, the backscatter increase at the end of the vegetation period can be modeled. The static calibration performance results in the form of RMSE improved from 1.92-2.25dB to 1.13-1.60dB for the non-static approach. The validation accuracy for the non-static parameter approach was evaluated with ubRMSE and ranges for all model combinations between 1.82 and 2.22dB. Validation results with SSRT as canopy model show better results in combination with all surface models when compared to respective combinations using WCM for the canopy part. Furthermore, it has been shown that the modeled backscatter results highly depend on the non-static empirical parameter. The evolution of the empirical parameter is similar for all surface models except Dubois95. At the beginning of the vegetation period, high values are observed, which decrease during the vegetation season to a minimum shortly before harvesting. Differences between the canopy model SSRT and WCM are noticeable in the form of higher variability (more outliers) of the empirical parameter from canopy model WCM. Furthermore, the empirical parameter for canopy model WCM has a relative maximum at the end of June. This increase of WCM's empirical parameter B can most likely be explained by not including information about the canopy height. Overall, the results of this study indicate that more complex models, like IEM_B as surface and SSRT as canopy model, provide the best results in our setup regarding ubRMSE and R 2 . It should be mentioned that the disadvantage of more complex models, by requiring more input parameters, were set to a minimum because of model parameter reduction by model adjustments and the use of literature values for rms height and single scattering albedo. Therefore, based on this study, we suggest using surface model IEM_B in combination with canopy model SSRT.
To accomplish a very dense time series of satellite acquisitions with a revisit time of 1.5 days, Sentinel-1 images with different incidence and azimuth angles were used. The models can account for backscatter changes due to different acquisition geometries only during the first half of the vegetation period. During the second half, a trend mismatch between Sentinel-1 and modeled backscatter is apparent in all model results. Therefore, it has to be stated that the used models in this study are only partially able to handle differences due to changes in radar acquisition geometries.
To take full advantage of dense time series provided, e.g., by incorporating different sensors on top of varying radar acquisition geometries, extended research on the suitability of certain RT models is required. Respectively, arising radiometric differences added to the disparity in radar acquisition geometries must be investigated, for that matter. To advance synergistic retrieval methods of SAR and optical data for soil moisture estimation, further research on incorporating non-static empirical parameter sets retrieved from optical sensors, to account for vegetation phenological states, is needed.