1. Introduction
Hailstorms are a frequent natural disaster, but a challenge to forecast, in the Canadian Prairies that can cause catastrophic damage to field crops. The specific atmospheric and geographic features of Western Canada are conducive to hailstorm activity, resulting in ‘hail belt’ regions [
1,
2]. In the warmest months of June and July, the cool dry air flowing from the Rocky Mountains encounters the hot moist surface air on the Prairies, resulting in an updraft of warm air. When the strong updrafts carrying moisture hit the high-altitude cold air, water droplets condense and freeze, creating hailstones. The Prairies represent 80% of Canadian farmland and suffer significant hail-to-crop damage and economic losses annually. The average annual agriculture loss by hail is about CAD 200 million [
3]. In 2020, the Canadian Crop Hail Association (CCHA) reported the total hail insured payment on 12,100 claims was almost CAD 250 million [
4]. Julian et al. [
1] investigated anthropogenic climate change over North America in the present and future. The study anticipated a more frequent occurrence of large hailstorms over the Canadian Rockies and Northern Plains in the future. With such numerous hailstorms, assessment of hail damage is challenging to the conventional on-farm inspector system, which requires assessing damage on individual plants within multiple locations of each damaged field.
In North America, the hail damage to field crops is generally categorized: (1) direct plant stand loss, (2) leaf defoliation, (3) stem cut-off, and (4) grain loss [
5,
6,
7]. The first damage assessment is normally conducted 7–10 days after hail occurs and involves counting plants and identifying the damage on individual plants at multiple locations within each field. The hail assessment requires careful, accurate, and reliable analysis to help farmers make decisions for replanting, harvest, and arriving at fair economic compensation for insurers. Conventional hail evaluation is typically used to meet the requirement, but it is considered time-consuming and labor-intensive. For these reasons, a less labor-intensive approach to estimate crop hail losses should benefit policymakers, insurance companies, and farmers.
Previous studies have found potential for using vegetation indices for estimation of hail damage in crops [
8,
9]. Hail damage causes physical defoliation and reduction in leaf chlorophyll, carotenoid, and polyphenol content [
10]. Spectral indices, such as Normalized Difference Vegetation Index (NDVI), Normalized Difference Water Index (NDWI), and Plant Senescence Radiation Index (PSRI), were found to be sensitive to crop stress [
11,
12,
13]. Sosa et al. [
14] used five Sentinel-1 microwave indices and five Sentinel-2 multispectral indices to detect the changes in vegetation caused by hail.
Satellite remote sensing provides a planet perspective to detect the dynamic changes in the field to regional scales, which could provide an inexpensive and more reliable approach to estimate hail damage in crops. In recent years, Sentinel-2 data have gained the attention of the remote sensing community for cropland mapping due to: (1) high spatial (10 m), temporal (5 days), and spectral resolution (13 bands); (2) free and open access; and (3) available for cloud computing (Google Earth Engine, GEE) [
15].
The primary remote sensing techniques commonly used in hail damage evaluation include change detection and time-series data analysis [
16,
17,
18]. Change detection analyzes the variability in two or multiple images for a specific area over a distinct time period. Spectral indices before and after an event are extracted and compared. The index difference is computed; a classification method, such as thresholding, is used to classify and generate damage maps. For example,
from Landsat 8 imagery was used to detect crop hail damage [
17]. Time-series data analysis has been used intensively in recent years for change detection. Time-series data analysis involves temporally dense monitoring of land surface dynamics over a defined period [
19]. In crop damage assessment, the difference of temporal patterns of damage and non-damage (temporal integration) is computed to estimate hail damage and its severity. Sosa et al. [
14] calculated the rate of index change (standard deviation) by differentiating each of the time-series obtained in each plot. The k-mean unsupervised classification on vegetation indices was used to detect the homogeneous hail damage zone on 91 plots. There was 87% damaged area strictly classified to similar degrees of damage. The area under curve (AUC) is frequently used to combine time-series disease progress into a single value and estimate disease progression’s effect on crop yield [
20]. Our research group [
21] has used AUC to measure the flowering duration and estimate crop yield in canola. Hail damage is a spatially discrete event and causes further different crop damages. So, it should be possible to assess hail damage severity by computing the area under the curve for a time-series of vegetative indices in hail-damaged fields, whereby measuring the realistic crop damage, crop growth restriction, and crop recovery.
Despite the increasing interest in using remote sensing for hail damage mapping, the major focus was on hail-strike locations, patterns, and damaged areas. Studies on modeling hail damage severity compared with the actual crop yield loss at a field scale are limited. Previous studies have assessed hail damage severity by using the estimated biomass and predicted yield models to quantify the yield loss by hail in the simulating hailstorm field experiments [
22,
23].
The main objectives of the study were (1) to explore NDVI and associated temporal profiles of the dominant crop types in the hail-damaged area, (2) to estimate crop damage using series data analysis and change detection methods, (3) to identify a suitable time frame on hail event for crop damage assessment, and (4) to compare the accuracy of AUC and delta temporal analysis for hail-to-crop damage severity assessment.
4. Discussion
The current study used Sentinel-2 to estimate hail damage in 54 fields (3402 hectares) of canola, wheat, and lentil crops near Assiniboia (49.6328° N 105.9921° W) on 4 July 2020. The crop damage due to hail was variable among the three crops. The effect of hail on vegetation indices was longer in canola and wheat compared to lentil. Wheat was found to recover at 33 days after the hail; whereas canola and lentil showed a delayed recovery.
Differential crop tolerance could be the reason for varying vegetation index correlation between crops. The correlation of the vegetative indices and crop damage showed no optimal vegetative index that best estimated all hail-crop damage with NDVI, NDWI, and PSRI all having the greatest r and the lowest RMSE in canola, wheat, and lentil, respectively. All vegetation indices use NIR-reflectance, which has been reported to be highly sensitive to the changes in the canopy structure and surface wetness due to hailstorm [
41].
It is interesting, though perhaps not unexpected, that the mesophyll in the leaf structure is sensitive to hail. Both the NDWI and PSRI reflect the changes of the mesophyll in the plant leaf. PSRI tracks the carotenoid to chlorophyll ratio in plants [
13]. In the plant leaves, carotenoid absorbs the blue-green to green wavelengths; although chlorophyll absorbs blue and red wavelengths. Since the leaf defoliation or scarred leaves by hail reduced the chlorophyll absorption of blue and red wavelengths, the ratio of carotenoid to chlorophyll in the chloroplast in the mesophyll increased; the PRSI would be raised after hail [
42]. Hail also affects the water content of the spongy mesophyll, which is the site of gaseous exchange. Tartachnyk et al. [
43] found that the stomata in the leaf closed during hail and remained closed for several hours, which further decreased water evapotranspiration rate. Araki et al. [
44] showed that rice had sharply reduced CO
2 exchange rate and mesophyll conductance in a flooded environment, and inhibited gas exchange resulted in high water content in the spongy mesophyll. The amount of water available in the internal leaf structure essentially controls the spectral reflectance in the SWIR interval of the electromagnetic spectrum, and SWIR reflectance is negatively related to leaf water content [
45]. Therefore, the NDWI index increased as hail appeared, and the leaf stomata remained closed in plants that were highly damaged.
NDVI has been used in previous hail damage assessments using remote sensing [
17,
46,
47]. The index is sensitive to plant health and is used to identify vegetation stress. Ref. [
8] detected hail areas by: (1) classifying the anomaly pixels in the NDVI layer compared to the background and (2) passing them through kernel filters of the histogram to determine damaged or undamaged areas. However, NDVI is general to plant health; thus, confounding factors, such as flooding stress, leaf damage by plant disease, nutrient deficiency, and others could influence this index [
48]. Moreover, NDVI may not be suitable for hail occurring near crop maturity.
Several studies have used methods such as quantifying the change of indices between pre- and post-hail or undamaged vs. damaged area to map the hailed area and separated different damage severity classes [
41,
49,
50]. Zhao et al. [
46] calculated the delta indices of pre- and post-hailstorm NDVI from HJ-1 multispectral images (30 m) to map hailed areas and, subsequently, categorized the damage into broad classifications of serious, moderate, and light, with a Kappa coefficient (κ) of 0.78. Choudhary [
41] calculated the percentages changes in NDVI and Normalized Difference Tillage Index (NDTI) between pre- and post- hail and estimated the affected area to five damage severity classes with an accuracy of 86.7%. Our study calculated the delta NDVI from a similar calculation but computed a finer level of damage severity with a coefficient correlation of 0.77~0.89. Additionally, the 12-day
NDVI indices also showed more accurate estimation than the 12-day AUC in canola and lentils. Corresponding to the fact that the canola and lentil did not recover as rapidly as wheat, the delta index was useful to assess the direct hail damage.
The comparison between the different hail damage approaches provided some interesting points. The present study calculated the AUC of the changed vegetative indices to estimate the hail severity. The crop damage estimation showed that the 32-day AUC had the highest r and lowest RMSE compared with the 12-day AUC, 12-day delta index, and 7-day delta index. The protracted AUC analysis captures both the direct damage and plant recovery in the temporal profile; in contrast, the delta approach only measured the direct damage. Another relative improvement of the AUC approach is fitting various vegetation indices. The difference of RMSE amongst the vegetation indices in the 12-day AUC assessment is minor compared to the 12-day delta assessment. Hence, the AUC approach’s estimation accuracy varies less with the vegetation indices.
Instead of comparing the plant reflectance variation by hail, another hail damage assessment compares the estimated plant biomass and yield with the potential yield from the previous years. Gobbo et al. [
51] estimated the biomass and yields in 90 hectares of hail-damaged corn with a precise −4.9~3.4% difference compared to the insurance company’s estimates. The author found the spatial variability within the field confounded the model and underestimated the actual biomass by 15.1%. The predicted yield model could only evaluate severe damage (no grain production), whereas the yield loss for moderate and low damage could not be estimated in the large-scale fields (8574 hectares) [
46]. The results of the study showed that changes in vegetation indices had a reasonably good correlation with crop inspector assessments. We believe that the ground survey from an insurance inspector is still essential and, when combined with index-based methodologies, more accurate estimations are possible.
We think that the ground survey from an insurance inspector is still essential for monitoring crop growth stages, evaluating potential yield, and providing aid to farming distress. The crop-damage assessment by satellite-based imagery can assist inspectors having quick overview on damage extent, as well as severity. Suitable vegetation indices explored from this study could support decision making in locating survey sites. Combining ground surveys with satellite imagery techniques makes more accurate and timely compensation possible.
Although crop damage by hail can be assessed, the study still has some limitations. First, the sample size is small, especially on lentil (six data points). As field data were collected from a single hail event, hail event time associated with crop growing stage has not been taken into consideration. Secondly, satellite data were collected from one source (Sentinel-2), so spatial and temporal resolution imagery can still be increased; environmental factors that contributed to the development of crops over time have not been included. Therefore, further research should focus on field data collection so that more sophisticated models, such as machine learning, can be used. This is possible as field data collection on crop damage by hail is conducted annually by the Municipal Hail Insurance. More satellite data sources, such as Sentinel-1, PlanetScope, Landsat, will be included to enhance data quality, as well as prediction power. Some satellite-based data on precipitation, temperature, soil moisture, evapotranspiration will be included so that damage prediction on a larger spatial scale is possible.