An Assessment of Drought Stress in Tea Estates Using Optical and Thermal Remote Sensing

: Drought is one of the detrimental climatic factors that affects the productivity and quality of tea by limiting the growth and development of the plants. The aim of this research was to determine drought stress in tea estates using a remote sensing technique with the standardized precipitation index (SPI). Landsat 8 OLI/TIRS images were processed to measure the land surface temperature ( LST ) and soil moisture index ( SMI ). Maps for the normalized difference moisture index ( NDMI ), normalized difference vegetation index ( NDVI ), and leaf area index ( LAI ), as well as yield maps, were developed from Sentinel-2 satellite images. The drought frequency was calculated from the classiﬁcation of droughts utilizing the SPI. The results of this study show that the drought frequency for the Sylhet station was 38.46% for near-normal, 35.90% for normal, and 25.64% for moderately dry months. In contrast, the Sreemangal station demonstrated frequencies of 28.21%, 41.02%, and 30.77% for near-normal, normal, and moderately dry months, respectively. The correlation coefﬁcients between the SMI and NDMI were 0.84, 0.77, and 0.79 for the drought periods of 2018–2019, 2019–2020 and 2020–2021, respectively, indicating a strong relationship between soil and plant canopy moisture. The results of yield prediction with respect to drought stress in tea estates demonstrate that 61%, 60%, and 60% of estates in the study area had lower yields than the actual yield during the drought period, which accounted for 7.72%, 11.92%, and 12.52% yield losses in 2018, 2019, and 2020, respectively. This research suggests that satellite remote sensing with the SPI could be a valuable tool for land use planners, policy makers, and scientists to measure drought stress in tea estates.


Introduction
Drought is one of the most deleterious climatic factors with regard to the changing climate that influences the growth of plants and ultimately affects the sustainable production of tea [1,2]. Drought influences plant physiological processes more than other environmental factors, altering cellular mechanisms and finally affecting the crop growth, productivity, and quality of tea [3,4]. The duration and magnitude of drought severity affects the yield and quality of tea, and the losses of crop yield caused by drought stress are becoming frequent and unpredictable because of climate change and increasing water stress in tea estates [5]. Tea (Camellia sinensis (L) O. Kuntze) is an evergreen perennial shrub cultivated in a wide range of soil and climatic conditions across tropical, subtropical and Mediterranean regions [6,7]. Both the quality and yield of tea depend on several factors, such as cultivar, cultural practices, and growth stages of plants, along with climate and season. Tea is a rainfed plantation crop, and its productivity depends mainly on weather conditions for optimal growth. Therefore, changes in climate, especially a reduction in rainfall, cause drought that affects the production of tea [8]. Bangladesh is one of the major tea-producing countries, ranking 9 th and producing 83 million kg of tea [9]. The information on land surface moisture content, which is an indispensable tool to measure drought severity [23].
Multiple studies have been conducted for assessing and monitoring drought using remote sensing technology. A study on NDVI-based agricultural drought assessment was performed to measure drought severity in the Palar Basin in Tamil Nadu, India, using IRS (Indian Remote Sensing Satellites)-1C & 1D [24]. Another study was performed in the western tracts of Tamil Nadu using satellite datasets of the National Oceanic and Atmospheric Administration (NOAA) [25]. Drought assessment was performed for hickory plantations in western Zhejiang, China, using a change detection approach from multitemporal Landsat 8 imagery [26]. A study was performed to investigate vegetation droughts in natural forests and rubber and oil palm plantations in peninsular Malaysia using MODIS datasets [27]. Another study was carried out in South Africa to measure the influence of drought on forest plantations using MODIS time series and climate data, suggesting that the normalized difference moisture index (NDMI) is a robust indicator to understand the interactions among precipitation, soil moisture, and water content in the plant canopy [28]. A comprehensive study on the impacts of drought and other disturbances on forest landscapes in North Carolina, USA, was conducted using MODIS and Landsat datasets [29]. In addition to these, several studies have been conducted for small-scale drought monitoring in tea plantations using wireless sensor networks (WSNs) [30][31][32]. The WSN-based drought monitoring system is not suitable for large and remote areas in many developing countries, such as Bangladesh, due to the disruption in power supply and internet networks, which are the major obstructions in the sensor network system. These studies provide information regarding the research gap in measuring drought severity for tea plantations to a greater extent using satellite remote sensing techniques.
To the best of our knowledge, no studies regarding drought monitoring for tea plantations on a large scale have been performed to identify drought severity in tea estates using satellite remote sensing technology. Furthermore, no standard drought index method has been developed for assessing drought severity. Therefore, a comprehensive study utilizing weather station-based precipitation datasets and multisatellite remote sensing datasets to develop a drought classification system was undertaken. The aim of the present study was to characterize the spatiotemporal pattern of drought in tea estates and to develop a drought classification system for tea plantations to ensure early warning for risk mitigation and sustainable tea production in Bangladesh.

Materials and Methods
Precipitation data collected from the Bangladesh Meteorological Department (BMD) were utilized to calculate the SPI and drought frequency. Both Landsat 8 OLI/TIRS and Sentinel-2 MSI satellite imagery were used as remote sensing datasets. Ground reference geographic positioning system (GPS) information was gathered during the field survey to extract the points over tea plantation areas to quantify the SMI, NDMI, LAI, and predicted yields.

Study Area
The study area is considered to lie within the tea-growing zone between latitudes 23 • 58 40 and 25 • 11 25 north and longitudes 91 • 16 09 and 92 • 30 17 east, and is located in the northeastern part of Bangladesh and consists of 21 subdistricts in the Sylhet Division ( Figure 1). The subdistricts include Moulvibazar, Kamalganj, Kulaura, Rajnagar, Sreemangal, Barlekha and Juri in the Moulvibazar district; Habiganj, Bahubal, Chunarughat, Madhabpur and Nabiganj in the Habiganj district; and Sylhet, Balaganj, Beanibazar, Fenchuganj, Golapganj, Gowainghat, Jaintiapur, Kanaighat and Zakiganj in the Sylhet district. The area is bordered by the Indian provinces of Meghalaya, Assam, and Tripura to the north, east, and south, respectively, subdistricts of Companiganj, Chatak, Jagannathpur, Baniachung and Lakhai to the west, and the Division of Chattogram to the southwest. There are 135 tea estates in the study area with 98,413 hectares of land, contributing approximately 85% of the total tea production of the country (63.42 million kg, 2018) [33]. The subdistricts include Moulvibazar, Kamalganj, Kulaura, Rajnagar, Sreemangal, Barlekha and Juri in the Moulvibazar district; Habiganj, Bahubal, Chunarughat, Madhabpur and Nabiganj in the Habiganj district; and Sylhet, Balaganj, Beanibazar, Fenchuganj, Golapganj, Gowainghat, Jaintiapur, Kanaighat and Zakiganj in the Sylhet district. The area is bordered by the Indian provinces of Meghalaya, Assam, and Tripura to the north, east, and south, respectively, subdistricts of Companiganj, Chatak, Jagannathpur, Baniachung and Lakhai to the west, and the Division of Chattogram to the southwest. There are 135 tea estates in the study area with 98,413 hectares of land, contributing approximately 85% of the total tea production of the country (63.42 million kg, 2018) [33]. This area is characterized by a sufficient amount of rainfall and favorable temperatures during the monsoon season that are congenial for tea production, but experiences drought stress during the winter season ( Figure 2). This area is characterized by a sufficient amount of rainfall and favorable temperatures during the monsoon season that are congenial for tea production, but experiences drought stress during the winter season ( Figure 2).

Precipitation Data and the SPI Calculation
The daily precipitation data for the local weather stations, Sylhet and Sreemangal, which represent the rainfall scenario of the study area, were collected from the Bangladesh Meteorological Department to obtain the monthly precipitation sums. Drought monitoring was performed employing the monthly precipitation sums to calculate the standardized precipitation index (SPI) for January 2018 to March 2021. Positive SPI values indicate wet conditions, and negative values denote drought periods with less precipitation than normal [34]. The precipitation data gathered for the local weather stations were utilized to classify wet and dry conditions [35]. Drought characterization was performed using the SPI to measure the severity, duration, and monthly coverage of drought of each category [36]. The accumulated SPI values were utilized to classify the drought category. The SPI value, ≤−1, refers to SPI-1, SPI-3, SPI-6, and SPI-12. Thus, the negative SPI values represent the corresponding drought severity for quantitative assessment of drought [37][38][39]. The SPI-1 computation was performed for 2018-2019, 2019-2020, and 2020-2021 using monthly precipitation data to investigate future drought events. The monthly precipitation sums could be used further to calculate the SPI utilizing precipitation records of the Bangladesh Meteorological Department. Drought classification was performed according to Razali et al. 2016 [27].
The SPI in a cumulative distribution is expressed using the gamma function. The gamma distribution is used to define the probability density or the function of frequency [39]. The gamma distribution can be represented as follows: Here, α > 0 and β > 0 refer to the shape parameter and scale parameter, respectively [40].
x denotes the amount of precipitation in millimeters in consecutive months i, and Г (α) is the gamma function. When xi = 0, the cumulative gamma distribution is unknown, and the cumulative probability is encountered. The SPI estimation also involves a

Precipitation Data and the SPI Calculation
The daily precipitation data for the local weather stations, Sylhet and Sreemangal, which represent the rainfall scenario of the study area, were collected from the Bangladesh Meteorological Department to obtain the monthly precipitation sums. Drought monitoring was performed employing the monthly precipitation sums to calculate the standardized precipitation index (SPI) for January 2018 to March 2021. Positive SPI values indicate wet conditions, and negative values denote drought periods with less precipitation than normal [34]. The precipitation data gathered for the local weather stations were utilized to classify wet and dry conditions [35]. Drought characterization was performed using the SPI to measure the severity, duration, and monthly coverage of drought of each category [36]. The accumulated SPI values were utilized to classify the drought category. The SPI value, ≤−1, refers to SPI-1, SPI-3, SPI-6, and SPI-12. Thus, the negative SPI values represent the corresponding drought severity for quantitative assessment of drought [37][38][39]. The SPI-1 computation was performed for 2018-2019, 2019-2020, and 2020-2021 using monthly precipitation data to investigate future drought events. The monthly precipitation sums could be used further to calculate the SPI utilizing precipitation records of the Bangladesh Meteorological Department. Drought classification was performed according to Razali et al. 2016 [27].
The SPI in a cumulative distribution is expressed using the gamma function. The gamma distribution is used to define the probability density or the function of frequency [39]. The gamma distribution can be represented as follows: Here, α > 0 and β > 0 refer to the shape parameter and scale parameter, respectively [40].
x denotes the amount of precipitation in millimeters in consecutive months i, and Γ (α) is the gamma function. When x i = 0, the cumulative gamma distribution is unknown, and the cumulative probability is encountered. The SPI estimation also involves a matching density function of the gamma likelihood to the frequency distribution of precipitation.
For an undefined gamma function, when x i = 0, the value of G (x i ) is expressed as follows: H( where q is the probability of zero. Then, the cumulative probability H(x i ) is transformed to the standard normal distribution for the SPI value [39]. To calculate the SPI, the cumulative probability distribution is transformed into a normal distribution that can be expressed as follows: When 0 < H(x) < 0.5, the following expression is true.
The drought frequency was calculated according to the following formula: Drought frequency (%) = Number of months affected by drought Total number of months × 100 (8)

Field Survey and Ground Reference Information
Ground reference information was collected during the 2019 field survey. GPS locations of tea estates were recorded around the study area using a handheld GPS locator (eTrex, Garmin, Olathe, AR, USA). Furthermore, this GPS information was used to extract point values from the remote sensing data employing the Spatial Analyst Tools in ArcGIS ( Figure 3). According to the field survey, among the 135 tea estates in the study area, 91 were located in Moulvibazar, 25 in Habiganj and 19 in Sylhet District ( Figure 4).

Remote Sensing Data
Multisatellite image data from the Landsat 8 OLI/TIRS and Sentinel-2 MSI were collected for drought estimation from December to February for the years 2018-2019, 2019-2020 and 2020-2021 ( Figure 3). Landsat 8 is equipped with the OLI (Operational Land Imager) and TIRS (Thermal Infrared Sensor) instruments. The OLI sensor provides image data in the coastal aerosol, visible (blue, green and red), NIR (near-infrared), and SWIR (short-wave infrared) spectra with a spatial resolution of 30 m and in a panchromatic band with a spatial resolution of 15 m. The TIRS sensor provides image data in two thermal infrared bands with a spatial resolution of 100 m. The TIRS dataset, registered to the OLI data, consists of geometrically, radiometrically, and terrain-corrected 12-bit products. Landsat 8 C1 level-1 images with path 136 and row 43 were processed to develop LST and SMI maps. The collected Landsat 8 image data were processed at a resolution of 30 m using the resampling technique in ArcGIS followed by subsequent mosaicking and masking. All selected bands from Landsat 8 were resampled to a 30 m resolution for uniformity in cell size and data. Furthermore, an algebraic operation, radiometric calibration, and atmospheric and geometric corrections were performed on the raster images. The average reflectance for the raster images was calculated using the raster calculator tool to minimize the spatial variability. All selected bands from Sentinel-2 image data resampled to a 20 m resolution were used to develop NDMI, LAI, and yield maps. To calculate the LST, the NDVI was retrieved measuring the reflectance of near-infrared, and red bands from the data obtained from the optical sensors of Landsat 8 OLI, followed by the retrieval of the proportion vegetation (PV), as well as land surface emissivity (LSE). The spectral radiance was derived from the long-wave infrared (LWIR) bands obtained from the TIRS of Landsat 8, followed by measuring the brightness tem-

Remote Sensing Data
Multisatellite image data from the Landsat 8 OLI/TIRS and Sentinel-2 MSI were collected for drought estimation from December to February for the years 2018-2019, 2019-2020 and 2020-2021 ( Figure 3). Landsat 8 is equipped with the OLI (Operational Land Imager) and TIRS (Thermal Infrared Sensor) instruments. The OLI sensor provides image data in the coastal aerosol, visible (blue, green and red), NIR (near-infrared), and SWIR (short-wave infrared) spectra with a spatial resolution of 30 m and in a panchromatic band with a spatial resolution of 15 m. The TIRS sensor provides image data in two thermal infrared bands with a spatial resolution of 100 m. The TIRS dataset, registered to the OLI data, consists of geometrically, radiometrically, and terrain-corrected 12-bit products. Landsat 8 C1 level-1 images with path 136 and row 43 were processed to develop LST and SMI maps. The collected Landsat 8 image data were processed at a resolution of 30 m using the resampling technique in ArcGIS followed by subsequent mosaicking and masking. All selected bands from Landsat 8 were resampled to a 30 m resolution for uniformity in cell size and data. Furthermore, an algebraic operation, radiometric calibration, and atmospheric and geometric corrections were performed on the raster images. The average reflectance for the raster images was calculated using the raster calculator tool to minimize the spatial variability. All selected bands from Sentinel-2 image data resampled to a 20 m resolution were used to develop NDMI, LAI, and yield maps. To calculate the LST, the NDVI was retrieved measuring the reflectance of near-infrared, and red bands from the data obtained from the optical sensors of Landsat 8 OLI, followed by the retrieval of the proportion vegetation (PV), as well as land surface emissivity (LSE). The spectral radiance was derived from the long-wave infrared (LWIR) bands obtained from the TIRS of Landsat 8, followed by measuring the brightness temperature (BT); and finally, the LST map was developed incorporating the values of PV and LSE with BT. On the other hand, the optical sensors from Sentinel-2 contributed to determining the NDMI and LAI, as well as predicted yield utilizing the visible (red), near-infrared, and short-wave infrared (SWIR) bands. All remote sensing indices were calculated from satellite data with less than 2% cloud cover for three consecutive years and analyzed in the ArcGIS ® environment.

Estimation of Land Surface Temperature (LST)
Land surface temperature is a controlling parameter of canopy temperature and water content in vegetation that can be used effectively to monitor soil moisture in vegetated land [41]. LST was calculated for the tea estates using temporal information from the Landsat 8 OLI/TIRS dataset with less cloud coverage [42,43]. There were two steps involved in the estimation of the LST. First, the NDVI data were retrieved using two wave bands, near-infrared (B5) and red (B4), according to the following formula: PV was calculated using the NDVI, which can be expressed as follows: Here, NDVI min and NDVI max were obtained from the respective NDVI images. The LSE was retrieved from the PV and calculated according to the following expression [44]: In the second step, the thermal bands were transformed to digital numbers to calculate the radiance. The spectral radiance was calculated from the thermal bands (B10 and B11), which can be expressed as follows: Lλ = ML + Q CAL + AL Lλ = 0.0003342 * Band10 + 0.1 and Lλ = 0.0003342 * Band11 + 0.1 (12) where Lλ denotes the spectral radiance of the TOA (top of atmosphere) at the sensor aperture, ML is the multiplicative rescaling factor of the corresponding bands that were derived from the metadata, Q CAL refers to the pixel values (DN) of the quantized and calibrated standard products, and AL is the additive rescaling factor of the specific band.
The BT was calculated using the following expression: where BT denotes the satellite brightness temperature in Celsius. K1 and K2 refer to the calibration constants in Kelvin derived from the metadata of the thermal conversion constants of the respective bands. Finally, the LST was calculated using the following formula:

Estimation of Soil Moisture Index (SMI)
Soil moisture is a potential determinant to understand drought severity and is primarily driven by precipitation through the process of infiltration as well as evapotranspiration [45,46]. Soil moisture is variable both spatially and temporally due to differences in soil properties and drainage [47]. Among various indices, the LST-based SMI is a widely used index to monitor soil moisture in vegetated land. The SMI is the proportion of the difference between the current soil moisture and the permanent wilting point to the field capacity as well as the residual soil moisture [48]. The value of the SMI ranges between 0 and 1, where 0 refers to very dry conditions and 1 denotes extremely wet conditions. The SMI is associated with the derived LST of a particular area that can be calculated using the following expression: where LST max and LST min refer to the maximum and minimum surface temperatures of the study area, respectively.

Estimation of Normalized Difference Moisture Index (NDMI)
The normalized difference moisture index determines the liquid water molecules in plant canopies that interact with solar radiation, and thus, water deficit results in a decrease in photosynthesis affecting the growth and development of plants [49]. The NDMI is a measure of water content in the spongy mesophyll tissues of vegetation canopies, which acts as a functioning indicator of climate and soil properties regulating water availability. The NDMI responds to variations in both water content (absorption of B8a) and spongy mesophyll tissues in plant canopies (reflectance of B11) [50,51]. The NDMI was calculated from Sentinel 2 satellite images using the following formula [52]: where NIR n is the near-infrared (narrow) reflectance of B8a with wavelengths ranging from 855 to 875 nm and SWIR 1 is the shortwave infrared reflectance of B11 with wavelengths of 1565-1655 nm.

Statistical Analysis
Mean values per pixel of individual estates were used for drought stress assessment in tea estates. Pearson correlation analyses were performed to assess the relationship between the soil moisture and vegetation canopy moisture. Simple Linear Regression analyses were also performed to show the association of the LAI and time-series yield information for yield prediction in the drought season.

Validation of Yield in the Drought Period
Yield data for 2018, 2019, and 2020 were collected from the Bangladesh Tea Board. The normalized difference vegetation index (NDVI) for the drought period was calculated from the reflectance of B4 (red) and B8 (near-infrared) of Sentinel-2 satellite imagery using the following expression.
The leaf area index (LAI) is an important determinant of tea yield and was also calculated from the NDVI using geospatial techniques. The LAI was calculated according to the equation used in previous studies based on the correlation between the LAI and NDVI utilizing the least square method [53,54].
Yield data were compared with the corresponding extracted LAI values in the scatter plot using regression analysis. The predicted yield maps in the drought period were developed using the linear regression equation in the ArcGIS environment. Finally, the extracted values from the yield maps were compared with the actual yield in the drought period to measure the yield loss in tea.

Standardized Precipitation Index (SPI) and Drought Classification
The results of drought frequency observed in the 39-month period for the Sylhet station demonstrated that 15 months belonged to the near-normal, 14 months normal, and 10 months moderately dry categories, which accounted for 38.46%, 35.90%, and 25.64%, respectively. In contrast, the Sreemangal station showed that the frequencies for near-normal, normal, and moderately dry conditions were 28.21%, 41.02%, and 30.77%, respectively ( Figure 5 and Table 1).

Standardized Precipitation Index (SPI) and Drought Classification
The results of drought frequency observed in the 39-month period for the Sylhet station demonstrated that 15 months belonged to the near-normal, 14 months normal, and 10 months moderately dry categories, which accounted for 38.46%, 35.90%, and 25.64%, respectively. In contrast, the Sreemangal station showed that the frequencies for near-normal, normal, and moderately dry conditions were 28.21%, 41.02%, and 30.77%, respectively ( Figure 5 and Table 1).

Land Surface Temperature (LST)
The average LST in the drought months (December to February) for the three consecutive years derived from the Landsat 8 remote sensing data showed that the surface temperature ranged between 14.25 and 27.68, 9.95 and 26.51, and 16.31 and 29.53, respectively ( Figure 6).

Land Surface Temperature (LST)
The average LST in the drought months (December to February) for the three consecutive years derived from the Landsat 8 remote sensing data showed that the surface temperature ranged between 14.25 and 27.68, 9.95 and 26.51, and 16.31 and 29.53, respectively ( Figure 6).

Soil Moisture Index (SMI)
The

Normalized Difference Moisture Index (NDMI)
The average NDMI in plant canopies observed in the drought period ranged between

Soil Moisture Index (SMI)
The

Normalized Difference Moisture Index (NDMI)
The average NDMI in plant canopies observed in the drought period ranged between

Spatial Correlation between the SMI and NDMI
The correlation coefficients between the SMI and NDMI were 0.84, 0.77, and 0.79 in the three consecutive years, respectively. This spatial correlation provides evidence that there is a strong interrelationship between moisture in soil and plant biomass during the drought period (Figure 9).

Validation of Tea Yield in the Drought Period
The validation of tea yield was performed using the NDVI-and LAI-based yield prediction technique (Figures 10-14).

Spatial Correlation between the SMI and NDMI
The correlation coefficients between the SMI and NDMI were 0.84, 0.77, and 0.79 in the three consecutive years, respectively. This spatial correlation provides evidence that there is a strong interrelationship between moisture in soil and plant biomass during the drought period (Figure 9).

Spatial Correlation between the SMI and NDMI
The correlation coefficients between the SMI and NDMI were 0.84, 0.77, and 0.79 in the three consecutive years, respectively. This spatial correlation provides evidence that there is a strong interrelationship between moisture in soil and plant biomass during the drought period (Figure 9).

Validation of Tea Yield in the Drought Period
The validation of tea yield was performed using the NDVI-and LAI-based yield prediction technique (Figures 10-14).

Validation of Tea Yield in the Drought Period
The validation of tea yield was performed using the NDVIand LAI-based yield prediction technique (Figures 10-14). The trendline obtained from the regression analysis between the observed LAI and time-series yield information during the drought period in tea estates showed a strong relationship with R 2 values, 0.66, 0.71, and 0.68, for the three consecutive years, respectively ( Figure 12). The regression analysis between LAI and the corresponding yield in the drought period of the three consecutive years shows a strong relationship between LAI and tea yield. According to the trajectories of predicted yield during the drought period, a lower yield was observed than the actual yield for 83 estates in 2018, whereas the number of responsible estates was 81 each in 2019 and 2020. According to the predicted yield, it was also evident that the tea estates in the study area demonstrated yield losses of 7.72%, 11.92%, and 12.52% in the three consecutive years, respectively ( Figure 14).  The trendline obtained from the regression analysis between the observed LAI and time-series yield information during the drought period in tea estates showed a strong relationship with R 2 values, 0.66, 0.71, and 0.68, for the three consecutive years, respectively ( Figure 12). The regression analysis between LAI and the corresponding yield in the drought period of the three consecutive years shows a strong relationship between LAI and tea yield. According to the trajectories of predicted yield during the drought period, a lower yield was observed than the actual yield for 83 estates in 2018, whereas the number of responsible estates was 81 each in 2019 and 2020. According to the predicted yield, it was also evident that the tea estates in the study area demonstrated yield losses of 7.72%, 11.92%, and 12.52% in the three consecutive years, respectively ( Figure 14). The trendline obtained from the regression analysis between the observed LAI and time-series yield information during the drought period in tea estates showed a strong relationship with R 2 values, 0.66, 0.71, and 0.68, for the three consecutive years, respectively ( Figure 12). The regression analysis between LAI and the corresponding yield in the drought period of the three consecutive years shows a strong relationship between LAI and tea yield. According to the trajectories of predicted yield during the drought period, a lower yield was observed than the actual yield for 83 estates in 2018, whereas the number of responsible estates was 81 each in 2019 and 2020. According to the predicted yield, it was also evident that the tea estates in the study area demonstrated yield losses of 7.72%, 11.92%, and 12.52% in the three consecutive years, respectively ( Figure 14).

Discussion
The use of the standardized precipitation index (SPI) with remote sensing datasets for drought monitoring was important for this study [55]. The combined utilization of the optical and thermal infrared remote sensing technology from Landsat 8 OLI/TIRS and Sentinel-2 was useful to fulfill the aim of this study. The thermal infrared sensors along with the optical sensors from Landsat 8 were used to measure the land surface temperature, and soil moisture content. On the other hand, the optical remote sensing technology utilized visible (red), near-infrared, and short-wave infrared sensors to measure the NDMI, NDVI, and LAI from Sentinel-2. This study also focused on assessing the relationship between soil and plant canopy moisture during the drought season, where both the optical and thermal infrared sensors contributed. In addition to this, the optical sensors from Sentinel-2, which are considered robust tools for providing high resolution remote sensing data, were used to perform yield loss assessment for validation of drought.
The optical and thermal infrared bands of remote sensing satellite data are predominantly used to monitor soil moisture, despite their susceptibility to the influence of clouds [56]. In contrast, microwave satellite remote sensing data with a higher temporal

Discussion
The use of the standardized precipitation index (SPI) with remote sensing datasets for drought monitoring was important for this study [55]. The combined utilization of the optical and thermal infrared remote sensing technology from Landsat 8 OLI/TIRS and Sentinel-2 was useful to fulfill the aim of this study. The thermal infrared sensors along with the optical sensors from Landsat 8 were used to measure the land surface temperature, and soil moisture content. On the other hand, the optical remote sensing technology utilized visible (red), near-infrared, and short-wave infrared sensors to measure the NDMI, NDVI, and LAI from Sentinel-2. This study also focused on assessing the relationship between soil and plant canopy moisture during the drought season, where both the optical and thermal infrared sensors contributed. In addition to this, the optical sensors from Sentinel-2, which are considered robust tools for providing high resolution remote sensing data, were used to perform yield loss assessment for validation of drought.
The optical and thermal infrared bands of remote sensing satellite data are predominantly used to monitor soil moisture, despite their susceptibility to the influence of clouds [56]. In contrast, microwave satellite remote sensing data with a higher temporal resolution can be used to measure soil moisture, recognizing the changes in the dielectric properties of surface soil water; however, this is not reliable due to the lower spatial resolution [57]. Therefore, the optical and thermal remote sensing bands of Landsat 8 OLI/TIRS were used to monitor drought for the tea plantations during the cloud-free winter season. Land surface temperature and vegetation indices such as NDMI are important indicators to measure soil moisture [58]. The utilization of Sentinel-2 datasets for NDMI measurements was useful because of its higher spatial resolution to monitor plant water content by measuring the reflectance that served the aim of this study.
It is evident that the study area receives adequate rainfall in the active growing period during the monsoon season but faces drought after the monsoon season; thus, this hampers productivity in tea during the drought months in the winter season. According to the drought frequency classification, this area is characterized by normal and near-normal conditions for 8 to 9 months and moderate drought conditions for 3 to 4 months in a year. As tea is sensitive to soil moisture stress, drought for 3 to 4 months causes yield loss, affecting both tea production and quality [1,6]. To validate the drought for tea estates, the LAI-based yield prediction model was adopted using remote sensing technique. The growth and development of tea leaves are reduced in the drought season due to the reduction in soil moisture. Another important reason for the reduction in leaf area of tea is that the tea bushes are heavily pruned at this time. These might be the reasons for the large differences between the actual yield and the predicted yield of tea.
The SMI-NDMI model based on optical and thermal infrared remote sensing is efficient and easy to calculate; however, it is susceptible to adverse weather conditions and cannot provide valid data in the presence of clouds [59]. Therefore, it was not possible to obtain year-round cloud-free data, which was an obstacle to preparing the year-round trajectories. To overcome this situation, integrated optical and thermal infrared remote sensing, as well as microwave remote sensing techniques, should be adopted to obtain year-round land surface soil moisture and plant water content, which would be an indispensable tool for drought monitoring in tea estates. This study utilized the precipitation-based SPI, as well as the remote sensing-based SMI and NDMI, to understand the water relationships under drought conditions rather than considering evaporation and transpiration losses of water from vegetated land, as the surface soil moisture distribution depends on the soil-vegetation-atmosphere of drought conditions. These issues should be addressed in future studies for precise monitoring of drought in tea estates.

Conclusions
This study suggests a drought stress assessment and classification system for tea plantations in Bangladesh using a standardized precipitation index and optical and thermal infrared remote sensing techniques. According to the results of drought frequency, the ratios of near-normal, normal, and moderately dry months for the Sylhet station were 38.46%, 35.90%, and 25.64%, respectively. In contrast, the Sreemangal station showed frequencies of 28.21%, 41.02%, and 30.77% for near-normal, normal, and moderately dry conditions, respectively. The spatial correlations between the SMI and NDMI were 0.84, 0.77, and 0.79 in the three consecutive years, respectively, indicating that there was a strong interrelationship between soil and vegetation canopy moisture during the drought period. This study also predicts that 83, 81, and 81 tea estates among the 135 estates showed lower yields than the actual yield during the drought period; thus, they incurred yield losses, which accounted for 7.72%, 11.92%, and 12.52% in 2018, 2019, and 2020, respectively. The results of the drought stress assessment, as well as drought classification for tea plantations in Bangladesh, could be significant for policy makers, land use planners, and scientists in the decision-making process in ensuring the early warning of drought conditions and thus would help tea growers to understand future drought events and to take the necessary initiatives for increasing production.