Landscape-Scale Crop Lodging Assessment across Iowa and Illinois Using Synthetic Aperture Radar (SAR) Images

: Crop lodging, the tilting of stems from their natural upright position, usually occurs after a heavy storm event. Since lodging of a crop seriously a ﬀ ects its yield, rapid assessment of crop lodging is valuable for farmers, policymakers, agronomists, insurance companies, and relief workers. Synthetic Aperture Radar (SAR) sensors have been recognized as valuable data sources for mapping lodging extent because of their good penetrating power and high-resolution remote sensing ability. Compared to other sources, SAR’s weather and illumination independence and large area coverage at ﬁne spatial resolution (3 m to 20 m) support frequent and detailed observations. Because of these advantages, SAR has the potential in supporting near real-time monitoring of lodging in ﬁelds when combined with automated image processing. In this study, a method based on change detection using modiﬁed Hidden Markov Random Field (HMRF) and Sentinel-1A data were utilized to identify lodging and map its extent. Results obtained have shown that when lodging occurs, the VH polarization’s backscatter ( σ VH) increases between the pre-lodging event image and the post-lodging event image. The increase in σ VH is due to the increase in volume scattering and vegetation-soil double bounce scattering resulting from the structural changes in the crop canopy. Using Sentinel-1A images and applying our proposed approach across several ﬁelds in Iowa and Illinois, we mapped the extent of the 2020 Derecho (wind storm) lodging disaster. In addition, we separated lodged regions into severely and moderately lodged areas. We estimated that approximately 2.56 million acres of corn and 1.27 million acres of soybean were lodged. Further analysis also showed the separation between un-lodged (healthy) ﬁelds and lodged ﬁelds. The observations in this study can guide future use of SAR-based information for operational crop lodging assessment. the VH mean polarization backscatter and 5 dB increase in the VV mean polarization backscatter for all the fields between the pre-lodging event date and the post- lodging event date. When we aggregated all the un-lodged fields and the lodged fields across Dallas county, we saw an increase of approximately 3 dB in the VH mean polarization backscatter and 1 dB in the VV mean polarization backscatter between the pre-lodging event date and the post-lodging event date. Taken together, these results suggest that differences in VH polarization and VV polarization can serve as useful lodging indicators at parcel- and landscape-level and enable rapid mapping of widespread lodging events using SAR data. Taken together, the results and the insights from this study demonstrate the practicability of using high resolution SAR remote sensing data for large scale identification of crop lodging. The approach employed is simple, effective, and can help decision-makers obtain critical information on lodging identification to support precision management in order to provide effective emergency relief. Future studies will evaluate if there is any improvement in the lodging detection performance


Introduction
With the increase in global population and the increase in food demand, the monitoring of agricultural activities has been of utmost importance. The increasing frequency and intensity of extreme weather events have also made monitoring of agricultural fields very critical. Lodging, the tilting of plant stems from their natural upright position, is a major yield limiting factor to crops such as corn, wheat, and barley [1][2][3]. Corn is vulnerable to lodging during its growing stages, particularly between early to late vegetative period [4]. Lodging in corn could be as a result (1) Understand the changes in backscatter over large-scale lodged fields and how to use the backscatter to classify lodging into severe or moderate categories. (2) Generate large-scale spatial extent maps of lodging using a change detection approach modified from our previous study [17] and determine the lodging rate (lodged crop per unit area) using the USDA's Crop Data Layer (CDL) map. (3) Qualitatively explore the relationship between high wind speed and lodged fields.
(4) Explore the capability of Sentinel-1A over an optical dataset like the one from Landsat-8. (5) Explore if lodged and un-lodged (healthy) fields differ between the pre-lodging event image and post-lodging event image.
In this paper, the proposed change detection method for lodging utilizes the concept of ratio image generation. The generation of ratio images suppresses background information while enhancing change information [17]. The generated ratio images were filtered using a non-local means filter [18,19] and classified into different classes using the Hidden Markov Random Field (HMRF). The HMRF considers the contextual information of neighboring pixels (i.e., a neighboring pixel is expected to have similar intensities and similar class labels) during classification. Published methods differ in their approach to extract change detection map. In the work done by Kasetkasem and Varshney [20], the authors used a MRF to model noiseless images for an optimal change image using the maximum a posteriori probability computation and the simulated annealing (SA) algorithm. In Zhao, et al. [21], the authors combined Voronoi Tessellation (VT) and HMRF based Fuzzy C-Means (FCM) algorithm (VTHMRF-FCM) for texture image segmentation. Similarly, in Yang and Yi [22], a novel method based on applying HMRF and generative adversarial network (GAN) on high-resolution SAR images was used for ship detection. To our knowledge and according to Chauhan, et al. [6], the use of satellite-based change detection for crop lodging is sparse. To date, there is no method available for corn and soybean lodging over large spatial areas.

Study Area
This study was carried out in Iowa and Illinois, spanning across 52 counties (Figure 1), covering the region affected by the 2020 Derecho lodging disaster which occurred on 10 August 2020. Iowa and Illinois are located in the Western part of the Corn Belt region of the United States and are part of the top agriculture production states. Iowa has a total agricultural land of 30.6 million acres, accounting for 92% of the state's total land area while Illinois has a total agricultural land of 27 million acres, accounting for 75% of the state's total land area. The two-dominant crops in both states are corn and soybean. Other crops like wheat and winter wheat are also planted but they account for a very small portion of all the cropland. In this study, the estimation of lodging rate (lodged crop per unit area) was only focused on corn and soybean. In Iowa and Illinois, the climate is cold and temperate. In Iowa, the wettest month with the most precipitation is June with an average precipitation of 107 mm and an average temperature of 82 • F, while the driest month is January with an average precipitation of 23 mm of rainfall and an average temperature of 31 • F. In Illinois, the wettest month with the most Remote Sens. 2020, 12,3885 4 of 15 precipitation is August in the Northern portion of Illinois with an average precipitation of 89 mm and an average temperature of 82 • F, while the driest month is February with an average precipitation of 49 mm of rainfall and an average temperature of 35 • F.

Field Data
In order to generate lodged crop per unit area, we used the 2018 CDL map because the 2020 CDL map is not yet available. Another reason for using the 2018 CDL map is that most farmers in Iowa and Illinois perform crop rotation. Fields where farmers consistently performed crop rotations for the last eight years were kept while fields where crop rotations were not consistent for the last eight years were removed. By doing this, we believe that corn and soybeans areas in 2018 should be similar to those of corn and soybean areas in 2020. The CDL data have a spatial resolution of 30 m, includes 132 detailed class labels, and was created from Landsat dataset using a decision tree algorithm trained on field samples [23]. The overall accuracy of the CDL dataset is about 95% for the United States Corn Belt. Apart from the CDL, we collected some ground truth data by field survey. Lodged fields were selected as samples for field observation. Mobile phone GPS was used for positioning the sampling points.

Remote Sensing Data
During 10 August 2020 and 11 August 2020, a derecho swept through Midwest that caused severe and widespread windstorms with some areas experiencing low-class tornadoes and heavy rain. As part of this, a windstorm was observed in Iowa and Illinois on 10 August 2020. The wind speed in both states ranged between 60 and 100 mph ( Figure 2). Wind speed data for this study were acquired from the National Weather Service (NWS). NWS used several sources like Iowa DOT, personal weather station, Automated Surface Observing Systems (ASOS), and Automated Weather Observing Systems (AWOS) to estimate the wind speed data. Since corn and soybean are already in their reproductive stages at the time of the windstorm, high wind speed led to severe lodging.

Field Data
In order to generate lodged crop per unit area, we used the 2018 CDL map because the 2020 CDL map is not yet available. Another reason for using the 2018 CDL map is that most farmers in Iowa and Illinois perform crop rotation. Fields where farmers consistently performed crop rotations for the last eight years were kept while fields where crop rotations were not consistent for the last eight years were removed. By doing this, we believe that corn and soybeans areas in 2018 should be similar to those of corn and soybean areas in 2020. The CDL data have a spatial resolution of 30 m, includes 132 detailed class labels, and was created from Landsat dataset using a decision tree algorithm trained on field samples [23]. The overall accuracy of the CDL dataset is about 95% for the United States Corn Belt. Apart from the CDL, we collected some ground truth data by field survey. Lodged fields were selected as samples for field observation. Mobile phone GPS was used for positioning the sampling points.

Remote Sensing Data
During 10 August 2020 and 11 August 2020, a derecho swept through Midwest that caused severe and widespread windstorms with some areas experiencing low-class tornadoes and heavy rain. As part of this, a windstorm was observed in Iowa and Illinois on 10 August 2020. The wind speed in both states ranged between 60 and 100 mph ( Figure 2). Wind speed data for this study were acquired from the National Weather Service (NWS). NWS used several sources like Iowa DOT, personal weather station, Automated Surface Observing Systems (ASOS), and Automated Weather Observing Systems (AWOS) to estimate the wind speed data. Since corn and soybean are already in their reproductive stages at the time of the windstorm, high wind speed led to severe lodging.
In order to cover the agricultural areas damaged by the windstorm, we acquired six Sentinel-1A images in Interferometric Wideswath (IW) instrument mode. Our Sentinel-1A images were acquired between 29 July 2020, and 22 August 2020. Each Sentinel-1A image collection had a resolution of 10 m, a dual polarization (VH and VV; V = Vertical, H = Horizontal), and consisted of Level-1 Ground Range Detected (GRD) scenes. As shown in Table 1 and Figure 3, three images were acquired before the lodging event (pre-lodging event) and three images were acquired after the lodging event (post-lodging event).
Remote Sens. 2020, 12, 3885 5 of 15 rain. As part of this, a windstorm was observed in Iowa and Illinois on 10 August 2020. The wind speed in both states ranged between 60 and 100 mph ( Figure 2). Wind speed data for this study were acquired from the National Weather Service (NWS). NWS used several sources like Iowa DOT, personal weather station, Automated Surface Observing Systems (ASOS), and Automated Weather Observing Systems (AWOS) to estimate the wind speed data. Since corn and soybean are already in their reproductive stages at the time of the windstorm, high wind speed led to severe lodging.   In order to cover the agricultural areas damaged by the windstorm, we acquired six Sentinel-1A images in Interferometric Wideswath (IW) instrument mode. Our Sentinel-1A images were acquired between 29 July 2020, and 22 August 2020. Each Sentinel-1A image collection had a resolution of 10 m, a dual polarization (VH and VV; V = Vertical, H = Horizontal), and consisted of Level-1 Ground Range Detected (GRD) scenes. As shown in Table 1 and Figure 3, three images were acquired before the lodging event (pre-lodging event) and three images were acquired after the lodging event (postlodging event).

Methods
A graphical overview of the framework used for this study is illustrated in Figure 4. Major steps of the framework include preprocessing, ratio image formation, and change detection classification approach. All the steps in this paper can be reproduced quickly and all analyses for the study were

Methods
A graphical overview of the framework used for this study is illustrated in Figure 4. Major steps of the framework include preprocessing, ratio image formation, and change detection classification approach. All the steps in this paper can be reproduced quickly and all analyses for the study were carried out using python programming language.
Remote Sens. 2020, 12, x FOR PEER REVIEW 6 of 15 In this paper, our proposed approach requires some parameters to be set beforehand. The parameters that need to be set include (i) the neighborhood size of the non-local means filtering step; (ii) the kernel size of the majority filter; (iii) the structuring element of the morphological filter; and, finally, (iv) the maximum number of allowed change classes. Please note that while we identified optimal settings for these parameters, we found that the performance of our algorithm does not critically depend on the exact choice for these variables. This is true for the following reasons: (i) as non-local means filtering is performed very early in the workflow, the impact of changes in the neighborhood size is mitigated by subsequent processing steps such as the application of mathematical morphology. Hence, we found that varying the neighborhood size from its optimal value changed system performance only slowly; (ii) the increase or decrease in the kernel size of the majority filter slowly decreases our change detection performance, yet this reduction of performance does not become significant unless the kernel size is increased tremendously; (iii) from an analysis of a broad range of data from different change detection projects we found that (1) a 5 × 5 pixel-sized structuring element of the morphological filter led to the most consistent results; and that (2) change detection performance changed slowly with deviation from the 5 × 5 pixel setting. Hence, while 5 × 5 pixel was found to be optimal, the exact choice of the window size is not critical for change detection success; finally, (iv) the maximum number of allowable change classes is a very uncritical variable as it merely sets an upper bound for a subsequent algorithm that automatically determines the number of distinguishable classes in a data set. By presetting this variable to 3 classes we ensure that changes as a result of lodging are captured.

Image Preprocessing
Sentinel-1A preprocessing was carried out using SeNtinel Application Platform (SNAP) software version 6.0 (https://step.esa.int/main/toolboxes/snap/), an open source common architecture provided by European Space Agency (ESA). The preprocessing step includes orbit file correction, GRD border noise removal, thermal noise removal, calibration, filtering using refined Lee filter, radiometric terrain correction (RTC), and geometric terrain correction (GTC). For more details on the preprocessing step and the importance of RTC, the reader is referred to Ajadi, Meyer and Webley [17].

Logarithmic Scaling and Ratio Image Formation
While both VH polarization and VV polarization were sensitive to crop lodging assessment, we only used the VH polarization of Sentinel-1A dataset in this study because it depicts crop phenology very well and, moreover, it is very sensitive to crop canopy structure. In order to increase the detectability of lodging and to suppress background information from SAR data, ratio images ( = , , ) were formed ( Figures 5 and 6) using the pre-lodging event image and of the post-lodging event image of similar geometry, respectively (Table 1 and Figure 3). Note that, due to the performed radiometric correction steps, images are not required to come from identical geometries. Also, because the pre-lodging and post-lodging event images are logarithmically scaled, the creation of ratio images is performed as a subtraction operation. Afterwards, a fast-non-local means filtering procedure was applied to all ratio images in order to filter out the speckle noise while In this paper, our proposed approach requires some parameters to be set beforehand. The parameters that need to be set include (i) the neighborhood size of the non-local means filtering step; (ii) the kernel size of the majority filter; (iii) the structuring element of the morphological filter; and, finally, (iv) the maximum number of allowed change classes. Please note that while we identified optimal settings for these parameters, we found that the performance of our algorithm does not critically depend on the exact choice for these variables. This is true for the following reasons: (i) as non-local means filtering is performed very early in the workflow, the impact of changes in the neighborhood size is mitigated by subsequent processing steps such as the application of mathematical morphology. Hence, we found that varying the neighborhood size from its optimal value changed system performance only slowly; (ii) the increase or decrease in the kernel size of the majority filter slowly decreases our change detection performance, yet this reduction of performance does not become significant unless the kernel size is increased tremendously; (iii) from an analysis of a broad range of data from different change detection projects we found that (1) a 5 × 5 pixel-sized structuring element of the morphological filter led to the most consistent results; and that (2) change detection performance changed slowly with deviation from the 5 × 5 pixel setting. Hence, while 5 × 5 pixel was found to be optimal, the exact choice of the window size is not critical for change detection success; finally, (iv) the maximum number of allowable change classes is a very uncritical variable as it merely sets an upper bound for a subsequent algorithm that automatically determines the number of distinguishable classes in a data set. By presetting this variable to 3 classes we ensure that changes as a result of lodging are captured.

Image Preprocessing
Sentinel-1A preprocessing was carried out using SeNtinel Application Platform (SNAP) software version 6.0 (https://step.esa.int/main/toolboxes/snap/), an open source common architecture provided by European Space Agency (ESA). The preprocessing step includes orbit file correction, GRD border noise removal, thermal noise removal, calibration, filtering using refined Lee filter, radiometric terrain correction (RTC), and geometric terrain correction (GTC). For more details on the preprocessing step and the importance of RTC, the reader is referred to Ajadi, Meyer and Webley [17].

Logarithmic Scaling and Ratio Image Formation
While both VH polarization and VV polarization were sensitive to crop lodging assessment, we only used the VH polarization of Sentinel-1A dataset in this study because it depicts crop phenology very well and, moreover, it is very sensitive to crop canopy structure. In order to increase the detectability of lodging and to suppress background information from SAR data, ratio images (X Ris = [X Ri1 , X Ri2 , X Ri3 ]) were formed ( Figures 5 and 6) using the pre-lodging event image and of the Remote Sens. 2020, 12, 3885 7 of 15 post-lodging event image of similar geometry, respectively (Table 1 and Figure 3). Note that, due to the performed radiometric correction steps, images are not required to come from identical geometries. Also, because the pre-lodging and post-lodging event images are logarithmically scaled, the creation of ratio images is performed as a subtraction operation. Afterwards, a fast-non-local means filtering procedure was applied to all ratio images in order to filter out the speckle noise while preserving the details. The fast-non-local means uses redundant information to reduce noise and restore the original noise-free image by performing a weighted average of pixel values, considering the spatial and intensity similarities between pixels [17].

Change Detection Classification Approach
In this study, we employed the Hidden Markov Random Field (HMRF) approach to perform our lodging-based change detection classification. The HMRF approach fully utilizes and enhances our previous change detection classification approach described in [17] by improving its robustness and sensitivity to false alarms. The method in [17] employed the Finite Gaussian Mixture (FGM) model for image classification. For example, if a ratio image contains dimensional vector of pixels with = 1, 2, … … . . . , being the set of pixel indices, then for each pixel in a class label is inferred using the conditional probability as shown in [17]. Each pixel in the FGM is independent from their neighboring pixels, meaning they do not consider the relationship of pixels within its neighboring system. To improve the FGM, we modified the HMRF approach proposed by Zhang, et al. [24] and employed it. By assuming a Gaussian distribution, the HMRF model is given by

Change Detection Classification Approach
In this study, we employed the Hidden Markov Random Field (HMRF) approach to perform our lodging-based change detection classification. The HMRF approach fully utilizes and enhances our previous change detection classification approach described in [17] by improving its robustness and sensitivity to false alarms. The method in [17] employed the Finite Gaussian Mixture (FGM) model for image classification. For example, if a ratio image contains dimensional vector of pixels with = 1, 2, … … . . . , being the set of pixel indices, then for each pixel in a class label is inferred using the conditional probability as shown in [17]. Each pixel in the FGM is independent from their neighboring pixels, meaning they do not consider the relationship of pixels within its neighboring system. To improve the FGM, we modified the HMRF approach proposed by Zhang, et al. [24] and employed it. By assuming a Gaussian distribution, the HMRF model is given by

Change Detection Classification Approach
In this study, we employed the Hidden Markov Random Field (HMRF) approach to perform our lodging-based change detection classification. The HMRF approach fully utilizes and enhances our previous change detection classification approach described in [17] by improving its robustness and sensitivity to false alarms. The method in [17] employed the Finite Gaussian Mixture (FGM) model for image classification. For example, if a ratio image X Ri contains N dimensional vector of pixels with I = {1, 2, . . . . . . . . . , N} being the set of pixel indices, then for each pixel i in X Ri a class label x i is inferred using the conditional probability as shown in [17]. Each pixel in the FGM is independent from their neighboring pixels, meaning they do not consider the relationship of pixels within its neighboring system. To improve the FGM, we modified the HMRF approach proposed by Zhang, et al. [24] and employed it. By assuming a Gaussian distribution, the HMRF model is given by where q(l| x Ni ) is a conditional probability mass function (pmf) for a class label given that x Ni are N neighbors for a pixel x i . The segmentation process of HMRF requires an initial estimate for the class labels x 0 and initial parameters (θ 0 = mean(µ 0 ), variance σ 2 0 ). In this research, we used K-means clustering approach to provide these initial labels and initial parameters. The update for parameter θ 0 was estimated iteratively using our Expectation-maximization (EM) algorithm [17]. The initial labels were updated iteratively using the maximum a posteriori (MAP) algorithm. The updated label is now used to solve forx that minimizes the total posterior energŷ where χ is the set of all possible configurations of labels,x is the estimated class label, U(X LR |x) is the likelihood energy, and U(x) is the prior energy function which is defined by clique potentials. A clique is a subset of nodes in which every node is connected to every other node. The clique potential was defined on pairs of neighboring pixels. We assumed that one pixel has at most eight neighboring pixels. The MAP algorithm iteration stops when Equation (2) converges or when the maximum iteration we set is reached. In this study, we classified each ratio image in X Ris into three classes namely no change, moderate change, and severe change. Finally, the classified ratio image was filtered using a majority filter with mathematical morphology to remove small isolated misclassified pixels. Majority filter with mathematical morphology works by first replacing the pixels in a neighborhood using majority of their adjacent neighboring pixels and then applying opening by reconstruction followed with closing by reconstruction. The order of first doing opening by reconstruction followed by closing by reconstruction was designed to reduce noise while preserving the geometric details in the image.
In this study, we set the kernel size of our majority filter to 3 and the size of our morphological kernel to 5. In the final classification map, we removed the no change class, and we multiplied each changed class (moderate change and severe change) by each crop (corn and soybean) mask derived from the CDL to estimate lodged crops per unit area.

Sentinel-1A Backscatter Analysis for Lodging Detection
The first instances of lodging were observed in the fields on 11 August 2020 due to the windstorm on 10 August 2020. Because of the temporal frequency of Sentinel-1, lodging assessment was done on acquisitions after the lodging event. As seen in the ratio images ( Figures 5 and 6), the VH polarization's backscatter (σVH) increased between the pre-lodging image and the post-lodging image. We predict the increase in value of σVH is based on two conditions: (1) The first condition is due to volume scattering from the crops. In a healthy crop, the orientation of the crop canopy is erect. When crop lodging occurs, the orientation of the canopy elements changes, the ears and the stems bend downwards and incline against each other. Since σVH is sensitive to structural changes, when lodging occurs, the volume scattering increases.
(2) The second condition is due to the high sensitivity between the interaction of vegetation and soil (double bounce effect). When a field is lodged, the vegetation-soil double scattering increases and the σVH increases as well.
The approach for the mapping of lodging as detailed in Section 3 captures the increase in backscatter as a result of the lodging event. When the lodging is severe, condition 1 and condition Remote Sens. 2020, 12, 3885 9 of 15 2 hold and changes in σVH are very high but when lodging is moderate, condition 1 holds and changes in σVH are moderate. The reader should note that both severe and moderate lodging are very detrimental and they will lead to yield loss. In this study, ratio images in X Ris were classified into three classes. The first class, which depicts no change, comprises mostly of urban areas, healthy crops, and standing water. The second class, which depicts moderate change, contains moderately lodged fields while the third class, which depicts severe change, contains severely lodged fields (Figures 7  and 8). After classification of the lodging areas, we then quantified the amount of corn and soybean that were moderately and severely lodged (

Qualitative Relationship between High Wind Speed and Lodged Fields
To qualitatively check if there is a relationship between high wind speed and lodged fields, we

Qualitative Relationship between High Wind Speed and Lodged Fields
To qualitatively check if there is a relationship between high wind speed and lodged fields, we

Qualitative Relationship between High Wind Speed and Lodged Fields
To qualitatively check if there is a relationship between high wind speed and lodged fields, we overlaid Figures 5b and 7b with the wind speed map (Figure 2). These overlay is shown in Figure 9a,b respectively. The result showed that our estimated lodged fields followed areas with increased wind speed (Figure 9a,b). In particular, areas with wind speed of 80+ mph captured most of our estimated severely and moderately lodged fields while areas with wind speed of 60+ mph captured few lodged fields.

Reliability of the Approach Employed in Mapping Lodging
To qualitatively access the performance of the Derecho lodged area mapped, we acquired eight ground truth data from field observation along the interstate I-80W. Due to limited amount of ground truth data collected, we were unable to perform extensive quantitative analysis of our method. While overlaying the eight ground truth points acquired with our lodged map, we observed that each point overlaid perfectly on an area that was lodged ( Figure 10). Furthermore, we compared an example of the lodged change detection map generated from

Reliability of the Approach Employed in Mapping Lodging
To qualitatively access the performance of the Derecho lodged area mapped, we acquired eight ground truth data from field observation along the interstate I-80W. Due to limited amount of ground truth data collected, we were unable to perform extensive quantitative analysis of our method. While overlaying the eight ground truth points acquired with our lodged map, we observed that each point overlaid perfectly on an area that was lodged ( Figure 10). Furthermore, we compared an example of the lodged change detection map generated from Sentinel-1A ratio image and optical Landsat-8 ratio image by using our classification approach. The Landsat-8 ratio image was generated from the green band of pre-lodging event image acquired on 10 July 2020, and the post-lodging event image acquired on 11 August 2020 ( Figure 11).
We chose the green band over every other band of Landsat-8 because the green band has a larger dynamic range which allowed for better separation of lodged and un-lodged areas. As can be seen from Figure 12, the lodging map generated from Landsat-8 ratio image and the Sentinel-1A ratio image were able to clearly detect lodging. While our approach was able to classify the Sentinel-1A ratio image into severely lodged and moderately lodged fields, we were unable to do so in the Landsat-8 ratio image. This observation showed that Sentinel-1A radar sensor has the advantage in reflecting the structural changes of lodging fields, while the Landsat-8 optical sensor only has the advantage in detecting the biochemical changes of lodging fields. All in all, for this study, we believe that the obvious feature for lodging is the structural change and the distribution pattern of lodged fields (severe and moderate) from Sentinel-1A image were similar to the distribution pattern of the lodged fields in the Landsat-8 image.
speed map in Figure 2.

Reliability of the Approach Employed in Mapping Lodging
To qualitatively access the performance of the Derecho lodged area mapped, we acquired eight ground truth data from field observation along the interstate I-80W. Due to limited amount of ground truth data collected, we were unable to perform extensive quantitative analysis of our method. While overlaying the eight ground truth points acquired with our lodged map, we observed that each point overlaid perfectly on an area that was lodged ( Figure 10). Furthermore, we compared an example of the lodged change detection map generated from Sentinel-1A ratio image and optical Landsat-8 ratio image by using our classification approach. The Landsat-8 ratio image was generated from the green band of pre-lodging event image acquired on 10 July 2020, and the post-lodging event image acquired on 11 August 2020 ( Figure 11).  ground truth data from field observation along the interstate I-80W. Due to limited amount of ground truth data collected, we were unable to perform extensive quantitative analysis of our method. While overlaying the eight ground truth points acquired with our lodged map, we observed that each point overlaid perfectly on an area that was lodged ( Figure 10). Furthermore, we compared an example of the lodged change detection map generated from Sentinel-1A ratio image and optical Landsat-8 ratio image by using our classification approach. The Landsat-8 ratio image was generated from the green band of pre-lodging event image acquired on 10 July 2020, and the post-lodging event image acquired on 11 August 2020 ( Figure 11).  We chose the green band over every other band of Landsat-8 because the green band has a larger dynamic range which allowed for better separation of lodged and un-lodged areas. As can be seen from Figure 12, the lodging map generated from Landsat-8 ratio image and the Sentinel-1A ratio image were able to clearly detect lodging. While our approach was able to classify the Sentinel-1A ratio image into severely lodged and moderately lodged fields, we were unable to do so in the Landsat-8 ratio image. This observation showed that Sentinel-1A radar sensor has the advantage in reflecting the structural changes of lodging fields, while the Landsat-8 optical sensor only has the advantage in detecting the biochemical changes of lodging fields. All in all, for this study, we believe that the obvious feature for lodging is the structural change and the distribution pattern of lodged fields (severe and moderate) from Sentinel-1A image were similar to the distribution pattern of the lodged fields in the Landsat-8 image.

Spatial-Extent of Lodging and Lodging Rate
Based on the reports by USDA Risk Management Agency (RMA), 57 counties in Iowa and Illinois were in the path of the storm which led to the Derecho lodging disaster. In the report by RMA, there are roughly 14 million acres of insured crops within those 57 counties. Also, based on the Storm Prediction Center preliminary storm reports and assessment from MODIS satellite imagery, the Iowa Department of Agriculture and Land Stewardship estimated that about 3.57 million acres of corn and 2.5 million acres of soybeans were likely to have been impacted by the severe wind on 10 August 2020. Based on our findings using Sentinel-1A data across 52 counties, we estimated that a total of approximately 2.56 million acres of corn and approximately 1.27 million acres of soybean were impacted during the Derecho lodging disaster (Figure 8 and Table 2). Furthermore, we observed that

Spatial-Extent of Lodging and Lodging Rate
Based on the reports by USDA Risk Management Agency (RMA), 57 counties in Iowa and Illinois were in the path of the storm which led to the Derecho lodging disaster. In the report by RMA, there are roughly 14 million acres of insured crops within those 57 counties. Also, based on the Storm Prediction Center preliminary storm reports and assessment from MODIS satellite imagery, the Iowa Department of Agriculture and Land Stewardship estimated that about 3.57 million acres of corn and 2.5 million acres of soybeans were likely to have been impacted by the severe wind on 10 August 2020. Based on our findings using Sentinel-1A data across 52 counties, we estimated that a total of approximately 2.56 million acres of corn and approximately 1.27 million acres of soybean were impacted during the Derecho lodging disaster (Figure 8 and Table 2). Furthermore, we observed that out of the 52 counties, only 32 counties have more than 10 thousand acres per county impacted by the storm. Within these 32 impacted counties, we leveraged our proprietary yield prediction capability to quantify the yield potential prior to the storm, and we were able to quantify the total bushels lost following the windstorm. We observed the following in our findings: (a) Corn acres impacted:~2.36 million acres impacted out of~4.90 million acres planted (48% acre impact). (b) Corn bushels impacted:~442.37 million bushels impacted out of~918.36 million bushels expected prior to the storm (48% reduction). (c) Soybean acres impacted:~1.27 million acres impacted out of~4.69 million acres planted (28% acre impact). (d) Soybean bushels impacted:~80.90 million bushels impacted out of~279.00 million bushels expected prior to the storm (29% yield reduction) As shown above and in Table 2, fields with moderate lodging were more frequent than fields with severe lodging and corn fields lodged more than soybean fields. We believe the reason for this difference is due to relatively more resistance of soybean to wind than corn.

Temporal Behavior of Un-Lodged (Healthy) and Lodged Fields throughout the Observation Period
The main difference in a field before and after lodging is the reduction of plant height. Therefore, the backscattering coefficients for lodged and un-lodged field will be different because for a lodged field, either one of the two conditions stated in Section 4.1 will hold. In this study, we compared nine randomly selected un-lodged fields to nine randomly selected lodged fields ( Figure 13) using both the VH mean polarization backscatter and VV mean polarization backscatter. While observing Figure 13, on the one hand, un-lodged fields showed no change at all in both the VH and VV mean polarization backscatter. On the other hand, we saw an increase of approximately 6 dB in the VH mean polarization backscatter and 5 dB in the VV mean polarization backscatter for all the fields between the pre-lodging event date and the post-lodging event date respectively. Even though there is a slightly higher difference in the VH mean polarization backscatter, the sensitivity of VH polarization to lodging is similar to the sensitivity of VV polarization to lodging. In this study, the increase in both the VH and VV mean polarization backscatters follows the reported observations in [8]. In Chauhan, et al. [8], the authors observed a clear linear trend of increasing VH polarization backscatter and VV polarization backscatter with the increase in the lodging severity. It is worth mentioning that corn and soybean in all the fields are in the grain-filling period or plateau stage. At this stage, the height and greenness of the crop should be fairly constant, and this is why we see no change in the un-lodged fields. To further show how lodged fields and un-lodged fields differ on a large scale, we generated the VH and VV mean polarization backscatter over Dallas county in Iowa ( Figure 14). By qualitatively observing the mean backscatter between the pre-lodging event date and the post-lodging event date in Figure 14, we saw an increase of approximately 3 dB in the VH mean polarization backscatter and 1 dB in the VV mean polarization backscatter. Based on the difference observed between the lodged fields and un-lodged fields in Dallas county, it is quite evident that our approach is sensitive enough to detect lodging. large scale, we generated the VH and VV mean polarization backscatter over Dallas county in Iowa ( Figure 14). By qualitatively observing the mean backscatter between the pre-lodging event date and the post-lodging event date in Figure 14, we saw an increase of approximately 3 dB in the VH mean polarization backscatter and 1 dB in the VV mean polarization backscatter. Based on the difference observed between the lodged fields and un-lodged fields in Dallas county, it is quite evident that our approach is sensitive enough to detect lodging.

Conclusions
The potential and feasibility of using SAR data for monitoring the Derecho lodging disaster was demonstrated in this study. Though crop lodging assessment using SAR data have been shown in a few studies earlier, studies that have assessed crop lodging using satellite-based data at large spatial scale are still sparse and knowledge relating to the changes in SAR signatures was lacking in literature. With the advent of dense SAR time series from Sentinel-1 and proposed missions like NASA ISRO Synthetic Aperture Radar (NISAR), there is an increased interest in exploring SAR for lodging identification. The main conclusions of this study are summarized below.
(1) The modified change detection approach used was shown to be capable of providing near realtime monitoring of the Derecho lodging disaster by generating detailed parameters, such as backscatter changes, lodging extent, and lodging rate in corn and soybean. The generated lodging extent maps from SAR showed both severely and moderately damaged fields. The use

Conclusions
The potential and feasibility of using SAR data for monitoring the Derecho lodging disaster was demonstrated in this study. Though crop lodging assessment using SAR data have been shown in a few studies earlier, studies that have assessed crop lodging using satellite-based data at large spatial scale are still sparse and knowledge relating to the changes in SAR signatures was lacking in literature. With the advent of dense SAR time series from Sentinel-1 and proposed missions like NASA ISRO Synthetic Aperture Radar (NISAR), there is an increased interest in exploring SAR for lodging identification. The main conclusions of this study are summarized below.
(1) The modified change detection approach used was shown to be capable of providing near real-time monitoring of the Derecho lodging disaster by generating detailed parameters, such as backscatter changes, lodging extent, and lodging rate in corn and soybean. The generated lodging extent maps from SAR showed both severely and moderately damaged fields. The use of CDL also allowed the estimation of lodged crop per unit area, showing relatively more lodging in corn fields than soybean fields. We believe the sensitivity of corn to lodging was caused by its unique structural characteristics (long vertical orientation of its stalk). We estimated that a total of approximately 2.56 million acres of corn and approximately 1.27 million acres of soybean were impacted during the Derecho lodging disaster. (2) The modified change detection approach used was reliable and the reliability can be seen by the similar distribution patterns of lodged fields in the Sentinel-1A imagery and Landsat-8 imagery. Furthermore, the generated lodged field maps show correlation with areas of extreme wind speed. (3) The backscatter difference between the timeseries of lodged and un-lodged (healthy) fields differ.
Our analyses from nine fields showed almost no change between the pre-lodging event image and post-lodging event image of an un-lodged field while we noted an approximately 6 dB increase in the VH mean polarization backscatter and 5 dB increase in the VV mean polarization backscatter for all the fields between the pre-lodging event date and the post-lodging event date. When we aggregated all the un-lodged fields and the lodged fields across Dallas county, we saw an increase of approximately 3 dB in the VH mean polarization backscatter and 1 dB in the VV mean polarization backscatter between the pre-lodging event date and the post-lodging event date. Taken together, these results suggest that differences in VH polarization and VV polarization can serve as useful lodging indicators at parcel-and landscape-level and enable rapid mapping of widespread lodging events using SAR data.
Taken together, the results and the insights from this study demonstrate the practicability of using high resolution SAR remote sensing data for large scale identification of crop lodging. The approach employed is simple, effective, and can help decision-makers obtain critical information on lodging identification to support precision management in order to provide effective emergency relief. Future studies will evaluate if there is any improvement in the lodging detection performance when SAR and multispectral datasets are fused. We will also undertake quantitative measures to further validate our results. Using the validated crop lodging area in this study as a reference, we will investigate the use of deep learning models to generate a standard near-real time framework for crop lodging assessment. We will also explore the use of Gamma distribution for fitting the EM algorithm in the HMRF rather than Gaussian distribution. With adequate ground truth information and further validation, large-scale crop lodging assessment can improve crop yield forecasting and crop insurance, ultimately benefitting food production and food security initiatives.