Development and Evaluation of an Agricultural Drought Index by Harnessing Soil Moisture and Weather Data

A new agricultural drought index was developed for monitoring drought impacts on agriculture in Oklahoma. This new index, called the Soil Moisture Evapotranspiration Index (SMEI), estimates the departure of aggregated root zone moisture from reference evapotranspiration. The SMEI was estimated at five locations across Oklahoma representing different climates. The results showed good agreement with existing soil moisture-based (SM) and meteorological drought indices. In addition, the SMEI had improved performance compared to other indices in capturing the effects of temporal and spatial variations in drought. The relationship with crop production is a key characteristic of any agricultural drought index. The correlations between winter wheat production and studied drought indices estimated during the growing period were investigated. The correlation coefficients were largest for SMEI (r > 0.9) during the critical crop growth stages when compared to other drought indices, and r decreased by moving from semi-arid to more humid regions across Oklahoma. Overall, the results suggest that the SMEI can be used effectively for monitoring the effects of drought on agriculture in Oklahoma.


Introduction
Drought events occur frequently in the Southern Great Plains of the United States, negatively impacting agricultural production and sustainability. This is mainly due to scarce surface water resources in this region. Groundwater resources, such as the Ogallala aquifer, help producers mitigate drought impacts [1]. However, these resources are being depleted at an unsustainable rate, especially during dry periods [2]. In Oklahoma, recorded drought history dates to the start of the 20th century, when major droughts were experienced in the decadal spans of 1910, 1930, 1950, 1960, and 1970. The beginning of the 21st century brought severe drought-related losses to Oklahoma's agriculture. For instance, the drought episode of 2001-2002 cost $210 million to the state's economy due to significant yield loss of winter wheat, alfalfa, and hay [2].
A more recent period of extreme and exceptional drought in Oklahoma was between 2011 and 2015, with devastating impacts on irrigated agriculture. As a result of this episode and consequent declines in reservoir water levels, the Luger-Altus Irrigation District in southwest Oklahoma could not deliver irrigation water during the 2011-2014 period [2]. In addition, groundwater resources experienced greater-than-usual depletions. The Ogallala aquifer underlying the Oklahoma Panhandle  Oklahoma (1981Oklahoma ( -2010. Source: PRISM raster map. The locations of study sites and major rivers are also demonstrated.

Input Data
Each Mesonet station was equipped with automated sensors to record precipitation, air temperature, relative humidity, wind speed, solar radiation, and soil moisture. Soil moisture was estimated using heat dissipation sensors (model 229-L, Campbell Scientific, Inc., Logan, UT, USA; accuracy < ± 0.5 °C) installed at three depths of 5, 25, and 60 cm below the soil surface [31]. Each sensor was calibrated before installation to obtain an accurate relationship between heat dissipation and soil matric potential, which was then converted to volumetric water content (VWC) based on the van Genuchten [32] method. Measurements from undisturbed soil cores at each site were used for determining water contents at field capacity (FC) and wilting point (WP) as well as the van Genuchten parameters [33]. Using accurate estimates of soil parameters based on in-situ  Oklahoma (1981Oklahoma ( -2010. Source: PRISM raster map. The locations of study sites and major rivers are also demonstrated. Table 1. Annual average daily air temperature, total precipitation, and total reference evapotranspiration (ET r ) measured at each study site over the study period (2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016).

Input Data
Each Mesonet station was equipped with automated sensors to record precipitation, air temperature, relative humidity, wind speed, solar radiation, and soil moisture. Soil moisture was estimated using heat dissipation sensors (model 229-L, Campbell Scientific, Inc., Logan, UT, USA; accuracy < ± 0.5 • C) installed at three depths of 5, 25, and 60 cm below the soil surface [31]. Each sensor was calibrated before installation to obtain an accurate relationship between heat dissipation and soil matric potential, which was then converted to volumetric water content (VWC) based on the van Genuchten [32] method. Measurements from undisturbed soil cores at each site were used for determining water contents at field capacity (FC) and wilting point (WP) as well as the van Genuchten parameters [33]. Using accurate estimates of soil parameters based on in-situ measurements is a major requirement in obtaining reliable SM based drought indices, a point that was highlighted in previous studies [11][12][13]19]. Table 2 presents the soil texture and the VWC at FC and WP for the three sensor installation depths at each site. Daily VWC data for the study period (2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016) collected at each of the three depths from the five sites were used in the analysis. The Mesonet soil moisture database had some missing values, but they accounted for only 3% of all collected data at each site on average. In addition, more than half of missing values occurred during periods with no measured precipitation. Thus, the gaps were filled using a linear interpolation technique. Meteorological data such as precipitation, air temperature, relative humidity, wind speed, and solar radiation were also measured by the Mesonet station at each site and were retrieved for the study period (2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016).
To investigate the impact of drought on crop yield, the relationship between drought indices and winter wheat production-a major crop in Oklahoma-was analyzed. County wise data on winter wheat production were downloaded from the website of the U.S. Department of Agriculture, National Agriculture Statistics Service (https://quickstats.nass.usda.gov/). The reported values represented both irrigated and non-irrigated winter wheat in Oklahoma. Segregated datasets of non-irrigated and irrigated winter wheat were not available for the complete study period. Furthermore, almost 96% of the winter wheat harvested area was non-irrigated in Oklahoma between 2001 and 2009. Hence, the obtained data can be safely assumed to represent dryland production, which is more susceptible to drought. Wister station (see Figure 1) was not included in the analysis due to insufficient data. Tian et al. [34] conducted a similar study in the same region under limitation of data availability and found the combined dataset of winter wheat yield (irrigated plus non-irrigated) suitable for such analysis based on the low percentage of irrigated area.

Soil Moisture Based Drought Indices
Three previously developed SM based drought indices and a new index were calculated for each site during the study period. The existing indices included Soil Water Deficit Index (SWDI), Water Deficit Index (d), and Normalized Soil Moisture (NSM), and the new index was called Soil Moisture Evapotranspiration Index (SMEI). The SWDI was developed by Martínez-Fernández et al. [11] and was estimated as: where θ is the aggregated VWC of soil profile, θ FC is the VWC at FC, and θ AWC is the available water content estimated as the difference between VWC at FC and WP (all in m 3 m −3 ). The weighted average of θ was estimated using the following approach as suggested by [11]: where θ 5 , θ 25 , and θ 60 are VWC at 5, 25, and 60 cm soil layers.

Water Deficit Index
This index was suggested by Cammalleri et al. [19]: where n is an empirical exponent (unitless) and θ 50 is estimated by averaging VWC between soil moisture thresholds as described by [19]; θ was aggregated for the soil profile based on Equation (2).

Normalized Soil Moisture
The NSM was proposed by Dutra et al. [16] as: where θ m,y is the VWC for the month m and the year y (m 3 m −3 ), θ m is the mean monthly VWC (m 3 m −3 ), and σ m is the standard deviation for all study years. This index is a mathematical imitation of the standardized precipitation anomaly used by Jones and Hulme [35]. Modeled soil moisture estimates of the top 289 cm of soil were used by Dutra et al. [16]. In the present study, however, in-situ measurements of soil moisture were aggregated based on Equation (2) and used in calculating NSM at two different time scales-one month (NSM-1) and three months (NSM-3). For NSM-3, the average θ m,y over three consecutive months was used in Equation (4) along with corresponding θ m and σ m .

Soil Moisture Evapotranspiration Index
The SMEI is proposed in this study based on soil moisture and reference evapotranspiration (ET r ) estimates. The SMEI estimation follows a three-step approach. At the first step, the aggregated daily root zone water storage (in units of water depth, such as mm) is calculated using a zone-weighted approach: where RZWS i is the root zone water storage on day i and L1, L2, and L3 are the layers in the crop root zone represented by each soil moisture sensor. In this study, values of 10, 30, and 40 cm were used for L1, L2, and L3, respectively, assuming a root zone depth of 80 cm [36]. This aggregation of RZWS was in agreement with soil moisture sensor depths of Oklahoma Mesonet as discussed by Ochsner et al. [37]. The values of the L parameters can be modified in other applications based on the number and the depth of soil moisture sensors as well as the root zone.
In the second step, the difference between monthly average RZWS and monthly cumulative ET r (both in units of water depth) is estimated for each month during the study period: where D m,y is the difference between monthly average of RZWS i and monthly sum of ET ri for month m and year y, n is the number of days in each month, and ET ri is the reference ET on day i. The standardized approach of the American Society of Civil Engineers (ASCE) outlined in [38] was adopted in estimating ET ri using the REF-ET program [39], and Mesonet weather data were used as input. The ASCE standardized approach is based on the simplified ASCE Penman-Monteith equation for a tall agricultural crop as the reference surface. Detailed information on the required parameters and their units can be found in [38]. The final step involves standardizing D m,y using the average and standard deviation to estimate the drought index: where D m is the monthly average, and σ m is the monthly standard deviation of D m,y for the same month among all study years. In essence, SMEI is similar to Palmer's Z-Index where ET is subtracted from precipitation at a monthly time-step to model soil water deficit (explained later). In SMEI, however, precipitation is replaced with soil moisture since the latter parameter better represents the source of water available to agricultural crops [40]. This is because smaller amounts of precipitation evaporate from soil and plant surfaces before reaching the root zone, and larger amounts can generate runoff or deep percolation below the root zone. In either case, the precipitated water will not be available for crop consumption. The inclusion of a water consumption/demand parameter (ET) along with a water source parameter (soil moisture) in the same index can also augment the drought signal [41].
In estimating SMEI, the common approach of presenting drought indices as anomalies [42] was followed, and z-scores were calculated by normalizing D m,y to a mean of zero and a standard deviation of one. Therefore, the SMEI values can be interpreted similar to the values of Standardized Precipitation Index (SPI) and Standardized Precipitation Evapotranspiration Index (SPEI, explained later), where negative values represent drought condition. In addition to monthly SMEI (SMEI-1), 3 month SMEI (SMEI-3) was calculated in order to understand the characteristics of drought events as drought recurrence decreases and its duration increases at longer time scales [43]. The 3 month SMEI was estimated by using the corresponding parameters over three consecutive months in Equations (7) and (8) (similar to NSM-3). Longer time scales (9 and 12 month) were not included, since the focus of this study was exclusively on agricultural drought, and the 3 month time scale was considered sufficient for this purpose [44].
One of the main merits of SMEI is the use of estimated soil moisture, which has been assigned with the highest ranking in terms of transparency of a drought index by Keyantash and Dracup [23]. In addition, the computational approach of SMEI is relatively simple. The tractability of this index may seem lower because soil moisture sensors have not been widely used by agricultural producers in the past. However, soil moisture sensing devices are becoming increasingly affordable [45], and their inclusion in public weather networks provides a valuable resource to all end-users.

Meteorological Drought Indices
Seven widely used meteorological drought indices were estimated in this study, including Atmospheric Water Deficit (AWD), Standardized Precipitation Index (SPI), Palmer's Z-Index, Palmer Drought Severity Index (PDSI), Self-Calibrated PDSI (scPDSI), Standardized Precipitation Evapotranspiration Index (SPEI), and U.S. Drought Monitor (USDM). The calculation of SPI, Z-Index, PDSI, scPDSI, and SPEI requires time series data for a minimum of 30 years for a thorough drought analysis [46]. As mentioned before, the present study focused on the 17 year period of 2000-2016. Therefore, PRISM (Parameter-elevation Regressions on Independent Slopes Model) [47] and COOP (The National Weather Service Cooperative Observer Network) 4 × 4 km grid monthly datasets for precipitation and temperature were used for constructing the timeline before the year 2000. The PRISM data are available at https://wrcc.dri.edu/wwdt/time/. Schneider and Ford [48] found that the PRISM database was useful in developing "climate informed decision support" in case of non-availability of long-term precipitation data. The root mean square error (RMSE) for precipitation and temperature data collected from Mesnoet and PRISM-COOP for 2000-2016 was estimated as 18.33 mm and 0.70 • C, respectively.

Atmospheric Water Deficit
The AWD constitutes a simple calculation where weekly sums of precipitation are subtracted from weekly sums of ET [18]. In this study, the AWD was calculated reversely (precipitation-ET) similar to Martínez-Fernández et al. [11], and monthly sums were used as opposed to weekly. The ASCE standardized reference ET (ET r ) was used in AWD estimations.

Standardized Precipitation Index
The SPI is based on long-term monthly precipitation data and is calculated by fitting the time series to a probability distribution [43]. This distribution is eventually converted into normal distribution to translate the anomalies of precipitation into a score range. In this study, SPI was calculated for one and three month periods using the SPI_SL_6 program (downloaded from http://drought.unl.edu/archive/ Programs/SPI/spi_sl_6.exe). The purpose of the calculations that were based on different time scales was to compare the drought indices in regard to short-and long-term drought events [49].

Palmer's Z-Index
Palmer's Z-Index (hereinafter Z-Index) requires monthly precipitation and temperature along with the latitude and the available water capacity of the soil as input data. A two-layer soil model approach is used to calculate the water balance [22]. Hydrologic accounting is performed to compute the parameters including ET, recharge, run-off, and loss of moisture from soil along with their potential values. The potential values are further transformed using climate dependent weighting factors. These transformed values are then summed as an equivalent of precipitation required for maintaining the normal soil moisture, which is further subtracted from the actual precipitation in order to estimate the monthly moisture departure [50]. Subsequently, the moisture departure is multiplied with empirically derived climatic characteristic, and the monthly moisture anomaly index known as Z-Index is estimated. Z-Index is sensitive to abrupt variation in soil moisture and thus is considered an appropriate index for monitoring agricultural droughts [51].

Palmer Drought Severity Index
The PDSI relies heavily on Z-Index. With the purpose of integrating the drought magnitude with time, the empirically obtained duration factor was combined with Z-Index, and PDSI values were calculated. These values represented the drought severity and were suitable for quantifying long-term droughts.

Self-Calibrated PDSI
This index makes PDSI more consistent for spatial comparisons of drought through replacing the empirical constants of PDSI with location-specific values [50]. Z-Index, PDSI, and scPDSI were all calculated using the sc-PDSI program downloaded from the website of the GreenLeaf project at University of Nebraska-Lincoln (http://greenleaf.unl.edu/).

Standardized Precipitation Evapotranspiration Index
The SPEI works on the same principle as SPI. However, it includes the long-term variations in air temperature and follows the water balance approach by subtracting potential ET (based on the Thornthwaite method) from precipitation [52]. One and three month SPEI (SPEI-1 and SPEI-3) were obtained in this study using the SPEI calculator [53].

U.S. Drought Monitor
The USDM data were obtained by downloading weekly GIS layers from the web portal of the U.S. Drought Monitor (http://droughtmonitor.unl.edu/Data/GISData.aspx). The USDM reports drought severity in one of the following categories: abnormally dry (D0), moderate drought (D1), severe drought (D2), extreme drought (D3), and exceptional drought (D4). The drought category for each study sites was extracted from USDM layers. The categories were converted to numerical values of 0, 1, 2, 3, and 4 for D0, D1, D2, D3, and D4 categories, respectively. Finally, the median numerical value in each month was used as the monthly value.

Performance Analysis
In order to evaluate the performance of the newly developed index, its relationship with the existing SM and the meteorological indices was investigated. Correlation analysis is a widely used method for measuring the degree of association among two variables [54]. Pearson correlation analysis was employed, and the following categories of correlation coefficients (r) suggested by Mavromatis [55] were adopted in the present study: very high (>0.9), high (0.7-0.9), moderate (0.5-0.7), weak (0.3-0.5), and very weak (<0.3). The performance of the new index was also assessed by studying its ability to track temporal and spatial variations in drought as compared with several existing indices.
The last step of the performance analysis of the new drought index included investigating its potential for assessing the impact of agricultural drought on winter wheat production between 2000 and 2016. Winter wheat is Oklahoma's largest crop with an average planted area of 2 to 2.4 million hectares every year [56]. This crop is usually planted in October and harvested in June of the next year, followed by a fallow period or a short-season summer cover crop. The spring period is considered critical due to the occurrence of jointing, heading, and flowering in winter wheat [57], and any water stress in these growth stages can cause significant fluctuations in crop yield [58]. Hence, the correlation between annual winter wheat production and the averaged drought indices was estimated both for the growing season (October to June) and for the critical period related to drought stress (March and April) similar to previous drought studies comparing the magnitude of drought and crop production [14,17,44,59].

Comparison with SM Indices
Correlation coefficients of the new and the existing SM based indices for each study site are shown in Figure 2. The new index SMEI-1 had high to very high correlations with NSM-1, represented by r values ranging from 0.70 to 0.92. Moderate to high correlations were noted with SWDI and d, with r ranges of 0.62-0.75 and 0.60-0.75, respectively. In comparison to SMEI-1, SMEI-3 had similar or weaker correlations with existing SM indices, except NSM-3, where better correlations were found (most likely due to similarity in time scales).
shown in Figure 2. The new index SMEI-1 had high to very high correlations with NSM-1, represented by r values ranging from 0.70 to 0.92. Moderate to high correlations were noted with SWDI and d, with r ranges of 0.62-0.75 and 0.60-0.75, respectively. In comparison to SMEI-1, SMEI-3 had similar or weaker correlations with existing SM indices, except NSM-3, where better correlations were found (most likely due to similarity in time scales).

Comparison with Meteorological Indices
The correlation coefficients among all studied SM and meteorological indices are graphed in Figure 3. Compared to existing SM indices, SMEI-1 had similar or stronger correlations with meteorological indices. Weak to moderate correlations for SMEI-1 were found with SPI-1 (0.41-0.59) and SPEI-1 (0.45-0.69) at all sites. Increasing the time scale of these meteorological indices from one to three months (SPI-3 and SPEI-3) improved correlations (0.61-0.71 and 0.64-0.77, respectively). A similar trend was noted for SMEI-3 and other SM indices, where longer time scale SPI and SPEI had stronger correlations. In case of PDSI, the correlations were improved for SMEI-3 (0.61-0.78) compared to SMEI-1 (0.56-0.72), potentially because the PDSI represents medium to long-term drought conditions. Weak to moderate correlations were found between SMEI-1 and AWD (0.37-0.62), but the values were comparable or stronger than those with existing SM indices. Moderate correlations between SMEI-1 and Z-Index (0.60-0.68) were found at all study sites. These results were expected, as Z-Index accounts for moisture departure from normal on a monthly basis, mainly in terms of ET and soil moisture losses [60]. As a result, Z-Index is more suitable for detecting variation in soil moisture over a short period of time [51,61], and the same applies to SMEI-1.
To help facilitate a better comparison between new and previously developed SM indices, box and whisker plots of their correlation coefficients with all meteorological indices at all sites are shown in Figure 4. The average correlation between SMEI and meteorological drought indices was 1.3 times higher in comparison to existing SM indices. The SMEI-1 had the largest mean r (0.60) among all SM indices followed by NSM-1 with 0.52. The SMEI-1 also had one of the smallest ranges (0.37-0.76), indicating relatively better performance under variable climatic conditions. In general, the r values had a larger range for longer time scales such as NSM-3 and SMEI-3. This was expected, since indices estimated on a scale of three months are likely to be more strongly correlated with medium to long-term meteorological drought indices but not with short-term indices.
to three months (SPI-3 and SPEI-3) improved correlations (0.61-0.71 and 0.64-0.77, respectively). A similar trend was noted for SMEI-3 and other SM indices, where longer time scale SPI and SPEI had stronger correlations. In case of PDSI, the correlations were improved for SMEI-3 (0.61-0.78) compared to SMEI-1 (0.56-0.72), potentially because the PDSI represents medium to long-term drought conditions. The correlations between meteorological and existing SM drought indices found in this study were similar to those reported in previous studies. For example, the range of SWDI-AWD correlation coefficients was 0.42-0.58 in the present study. Martínez-Fernández et al. [11] found weak to moderate correlations (0.40-0.57) between weekly SWDI and AWD in a semi-arid region with heterogeneous soils. Mishra et al. [62] reported a correlation coefficient of 0.56 for the same indices under humid subtropical climate. In addition, moderate to high correlation coefficients (0.60-0.80) were reported for NSM-PDSI by Dutra et al. [16] for a Mediterranean region, and this study found similar r values ranging between 0.47 and 0.74. indices followed by NSM-1 with 0.52. The SMEI-1 also had one of the smallest ranges (0.37-0.76), indicating relatively better performance under variable climatic conditions. In general, the r values had a larger range for longer time scales such as NSM-3 and SMEI-3. This was expected, since indices estimated on a scale of three months are likely to be more strongly correlated with medium to longterm meteorological drought indices but not with short-term indices. The correlations between meteorological and existing SM drought indices found in this study were similar to those reported in previous studies. For example, the range of SWDI-AWD correlation coefficients was 0.42-0.58 in the present study. Martínez-Fernández et al. [11] found weak to moderate correlations (0.40-0.57) between weekly SWDI and AWD in a semi-arid region with heterogeneous soils. Mishra et al. [62] reported a correlation coefficient of 0.56 for the same indices under humid subtropical climate. In addition, moderate to high correlation coefficients (0.60-0.80) were reported for NSM-PDSI by Dutra et al. [16] for a Mediterranean region, and this study found similar r values ranging between 0.47 and 0.74. Spatial variations in correlation strength between SMEI-1 and meteorological drought indices at different sites were in accordance with the available literature. The strongest correlations of SMEI-1 and PDSI were found at Goodwell and Wister (r > 0.72), whereas Hollis, El Reno, and Pawnee had r values of 0.66, 0.58, and 0.56, respectively. Greater correlations with PDSI at arid and humid sites were most likely due to the inherent calibration of PDSI in extreme climates [50]. The SMEI-1 showed Spatial variations in correlation strength between SMEI-1 and meteorological drought indices at different sites were in accordance with the available literature. The strongest correlations of SMEI-1 and PDSI were found at Goodwell and Wister (r > 0.72), whereas Hollis, El Reno, and Pawnee had r values of 0.66, 0.58, and 0.56, respectively. Greater correlations with PDSI at arid and humid sites were most likely due to the inherent calibration of PDSI in extreme climates [50]. The SMEI-1 showed the strongest correlation with SPEI-3 (r = 0.76) at Goodwell and comparatively weaker correlation at Wister (r = 0.63). This may have been due to overestimation of potential evapotranspiration (PET) by the Thornthwaite equation for SPEI in humid regions [63]. The SMEI correlations with SPEI would probably have been higher if the Penman-Monteith equation was used to estimate the reference ET in SPEI calculations, as recommended by Beguería et al. [64].

Temporal Tracking of Drought
Due to the occurrence of multiple droughts of varying magnitudes between 2010 and 2015, this period was selected to further evaluate the performance of newly developed SMEI in tracking temporal variations in drought. Changes in SMEI-1 during the selected period were compared against fluctuations in USDM, two SM drought indices, and one meteorological drought index. The two SM indices were NSM-1 and SWDI. The d index was excluded from the comparison because it had the smallest average r (Figure 4). The single meteorological index used in comparison was the Z-Index, since it incorporates the soil-water balance in its calculations. The USDM timeline was included for comparison with the overall drought conditions. It should be noted that USDM estimates are based on a combination of meteorological drought indices, field observations, and stakeholder feedback [8].
According to USDM, the drought condition gradually started progressing between October 2010 and June 2011 at Goodwell, starting from D0 and reaching at D4 (Figure 5a). The drought continued for almost five years and eventually ended by June 2015. All drought indices showed multiple episodes of severe droughts during this period. The SMEI-1 and the Z-Index compared well with each other, and both indices captured the months of June and July 2011 as the worst drought-stricken months. The Z-Index values had a more rapid fluctuation, mainly due to the sensitivity of this index to unusually dry or wet months [50]. The NSM-1 had a relatively slow response, indicating wet conditions during peak drought months and gradually transforming into drought conditions. This shows the advantage of amplifying drought response by including ET in drought indices [41]. The SWDI trended similar to NSM-1 and signaled decreases in drought magnitude; however, it never confirmed a non-drought period. Almost all indices clearly indicated a significant wet period in May 2015, which marked the end of the drought. each other, and both indices captured the months of June and July 2011 as the worst drought-stricken months. The Z-Index values had a more rapid fluctuation, mainly due to the sensitivity of this index to unusually dry or wet months [50]. The NSM-1 had a relatively slow response, indicating wet conditions during peak drought months and gradually transforming into drought conditions. This shows the advantage of amplifying drought response by including ET in drought indices [41]. The SWDI trended similar to NSM-1 and signaled decreases in drought magnitude; however, it never confirmed a non-drought period. Almost all indices clearly indicated a significant wet period in May 2015, which marked the end of the drought.  During all major drought episodes at Hollis, El Reno, and Pawnee sites, the SMEI-1 values declined in correspondence with other drought indices (Figure 5b-d). In general, both SMEI-1 and Z-Index had similar trends in tracking drought intensification and reliefs during the 2010-2015 period, whereas NSM-1 and SWDI showed damped response towards short-term variations in drought condition. Since SMEI-1 is based on measured daily soil moisture and ET, its capability to identify the drought is more rigorous in comparison to Z-Index. In addition, SMEI relies on the ASCE Standardized Penman-Monteith method for estimating ET r , while Z-Index uses the less sophisticated but also less accurate Thornthwaite method [65].

Spatial Tracking of Drought
Two drought episodes in February 2006 and July 2011 were selected to investigate spatial variability in drought indices. The corresponding USDM maps were used in the analysis. According to Oklahoma Water Resources Board [66], the period of September 2005 to March 2006 was the driest cool growing season in Oklahoma since 1921. The departure from the normal precipitation during this season ranged from −366 mm for the Southeast division (Wister) to −105 mm for the Panhandle region (Goodwell), and the severity of drought gradually reduced from east to west. In February 2006, Wister and Pawnee were under D4 and D3 categories, respectively. The remaining sites were facing D2 category. The spatial variability of SMEI-1 was in agreement with USDM ( Figure 6), ranging from −3.34 at Wister to −1.14 at Goodwell. NSM and Z-Index also showed a similar pattern. However, SWDI had a nearly opposite trend, showing an increase in drought severity from east to west. This was most probably because this index is solely based on soil moisture availability in the root zone and does not include other parameters. In addition, it is not normalized based on the past data at each site. Despite having smaller departures of soil moisture from average, the magnitude of soil moisture was smaller in western sites compared to eastern ones due to their natural aridity. Hence, SWDI signaled a more severe drought. Keyantash and Dracup [23] mentioned that indices that are normalized are more appropriate for comparing across locations. data at each site. Despite having smaller departures of soil moisture from average, the magnitude of soil moisture was smaller in western sites compared to eastern ones due to their natural aridity. Hence, SWDI signaled a more severe drought. Keyantash and Dracup [23] mentioned that indices that are normalized are more appropriate for comparing across locations. An opposite spatial trend occurred in July 2011, and drought severity gradually decreased from Panhandle and Northwest regions towards the southeast. Based on USDM, Goodwell and Hollis experienced D4, El Reno was near the transition zone between D3 and D4, Pawnee was under D3, and Wister was under D1. The spatial pattern of SMEI-1 was in general compliance with the USDM, with the largest values of −0.41 estimated at Wister (Figure 7). The Z-Index followed a similar trend. However, the Z-Index had a value of −3.0 at Wister, indicating extreme drought conditions [51]. SWDI reflected an intensification of drought moving from Wister to Hollis but a less severe drought for Goodwell. The NSM showed non-drought condition at Goodwell, mainly due to the fact that aggregated soil moisture for July 2011 merely deviated from the long-term mean. Therefore, it was not able to detect the drought. An opposite spatial trend occurred in July 2011, and drought severity gradually decreased from Panhandle and Northwest regions towards the southeast. Based on USDM, Goodwell and Hollis experienced D4, El Reno was near the transition zone between D3 and D4, Pawnee was under D3, and Wister was under D1. The spatial pattern of SMEI-1 was in general compliance with the USDM, with the largest values of −0.41 estimated at Wister (Figure 7). The Z-Index followed a similar trend. However, the Z-Index had a value of −3.0 at Wister, indicating extreme drought conditions [51]. SWDI reflected an intensification of drought moving from Wister to Hollis but a less severe drought for Goodwell. The NSM showed non-drought condition at Goodwell, mainly due to the fact that aggregated soil moisture for July 2011 merely deviated from the long-term mean. Therefore, it was not able to detect the drought.

Comparison with Crop Production
When considering the entire winter wheat growing season, SMEI had the strongest correlations with crop yield in comparison to existing drought indices at all sites ( Figure 8). The ranges of r were 0.41-0.79, 0.20-0.73, 0.37-0.69, and 0.35-0.64 for SMEI-1, Z-Index, SWDI, and NSM-1, respectively. These results were in agreement with previous studies. For example, Z-Index had r values as high as 0.65 when compared with the wheat yield in varying climates using modeled data embedded in the 0.5° grid [67]. In addition, Carrão et al. [14] compared the remote sensing-based SM drought index with non-irrigated wheat yield from a humid region in South America for a complete growing season and found an r value of 0.82. Also, a general spatial trend was noted in examining correlation coefficients among drought indices and winter wheat yields, indicating a weakening of the correlation moving from the semi-arid Panhandle to more humid regions (Figure 8). This spatial variation was most probably due to the fact that winter wheat production is more sensitive to drought conditions in water limited environments, as root zone water availability is the major limiting factor for crop growth in rain-fed agriculture [68]. The spatial trend was in accordance with the findings of Tian et al. [34] who compared the winter wheat yield from Texas and Oklahoma with multiple drought indices (including Z-Index) and found stronger correlations (r > 0.5) in the western parts of the study area and weaker correlations towards the east.

Comparison with Crop Production
When considering the entire winter wheat growing season, SMEI had the strongest correlations with crop yield in comparison to existing drought indices at all sites ( Figure 8). The ranges of r were 0.41-0.79, 0.20-0.73, 0.37-0.69, and 0.35-0.64 for SMEI-1, Z-Index, SWDI, and NSM-1, respectively. These results were in agreement with previous studies. For example, Z-Index had r values as high as 0.65 when compared with the wheat yield in varying climates using modeled data embedded in the 0.5 • grid [67]. In addition, Carrão et al. [14] compared the remote sensing-based SM drought index with non-irrigated wheat yield from a humid region in South America for a complete growing season and found an r value of 0.82. Also, a general spatial trend was noted in examining correlation coefficients among drought indices and winter wheat yields, indicating a weakening of the correlation moving from the semi-arid Panhandle to more humid regions (Figure 8). This spatial variation was most probably due to the fact that winter wheat production is more sensitive to drought conditions in water limited environments, as root zone water availability is the major limiting factor for crop growth in rain-fed agriculture [68]. The spatial trend was in accordance with the findings of Tian et al. [34] who compared the winter wheat yield from Texas and Oklahoma with multiple drought indices (including Z-Index) and found stronger correlations (r > 0.5) in the western parts of the study area and weaker correlations towards the east. coefficients among drought indices and winter wheat yields, indicating a weakening of the correlation moving from the semi-arid Panhandle to more humid regions (Figure 8). This spatial variation was most probably due to the fact that winter wheat production is more sensitive to drought conditions in water limited environments, as root zone water availability is the major limiting factor for crop growth in rain-fed agriculture [68]. The spatial trend was in accordance with the findings of Tian et al. [34] who compared the winter wheat yield from Texas and Oklahoma with multiple drought indices (including Z-Index) and found stronger correlations (r > 0.5) in the western parts of the study area and weaker correlations towards the east. Comparatively larger correlations of SMEI-1 showed that the newly developed index is a useful tool for decision makers and growers to monitor drought variability and its effect on the production Comparatively larger correlations of SMEI-1 showed that the newly developed index is a useful tool for decision makers and growers to monitor drought variability and its effect on the production of winter wheat. However, the winter wheat farmers in Oklahoma would also be interested to know if they possibly could predict the anticipated yield in case of occurrence of droughts in those months that are critical for the crop growth, as they have to make important decisions by the spring period regarding grain production or utilizing the crop for grazing [28]. The correlation analysis conducted between winter wheat production and the spring period's (March and April) drought indices values showed that the SMEI-1 had the largest correlations ( Figure 9). For example, very high correlation was found at Goodwell with an r value of 0.92, whereas high (0.72) to moderate (0.61) correlations were noted for Hollis and El Reno, respectively, and weaker correlations (0.35) were noted for Pawnee. Narasimhan and Srinivasan [17] also analyzed the impact of drought on the winter wheat yield for the same crop growth period in Texas by using SM and ET deficit based indices and found high correlations (r > 0.75). These findings suggest that SMEI-1 could be suitable for predicting winter wheat production in drought-prone areas of Oklahoma. of winter wheat. However, the winter wheat farmers in Oklahoma would also be interested to know if they possibly could predict the anticipated yield in case of occurrence of droughts in those months that are critical for the crop growth, as they have to make important decisions by the spring period regarding grain production or utilizing the crop for grazing [28]. The correlation analysis conducted between winter wheat production and the spring period's (March and April) drought indices values showed that the SMEI-1 had the largest correlations ( Figure 9). For example, very high correlation was found at Goodwell with an r value of 0.92, whereas high (0.72) to moderate (0.61) correlations were noted for Hollis and El Reno, respectively, and weaker correlations (0.35) were noted for Pawnee. Narasimhan and Srinivasan [17] also analyzed the impact of drought on the winter wheat yield for the same crop growth period in Texas by using SM and ET deficit based indices and found high correlations (r > 0.75). These findings suggest that SMEI-1 could be suitable for predicting winter wheat production in drought-prone areas of Oklahoma.

Conclusions
Considering the need for an improved drought monitoring tool for Oklahoma, a new agricultural drought index-the Soil Moisture Evapotranspiration Index (SMEI)-was introduced and evaluated in this study. It was estimated as the normalized difference of monthly root-zone soil moisture and reference evapotranspiration (ETr). The SMEI was estimated along with several existing soil moisture (SM) based and meteorological indices during the 2000-2016 period at five sites across Oklahoma, U.S., with variable climatic conditions. The SMEI had strong correlations with existing SM indices. The correlations between SMEI and meteorological drought indices were 33% stronger on average when compared to existing SM indices. Higher correlation coefficients of SMEI were observed at all study sites, suggesting that SMEI performs well under variable climatic conditions experienced across Oklahoma. The correlation coefficients between existing SM and meteorological

Conclusions
Considering the need for an improved drought monitoring tool for Oklahoma, a new agricultural drought index-the Soil Moisture Evapotranspiration Index (SMEI)-was introduced and evaluated in this study. It was estimated as the normalized difference of monthly root-zone soil moisture and reference evapotranspiration (ET r ). The SMEI was estimated along with several existing soil moisture (SM) based and meteorological indices during the 2000-2016 period at five sites across Oklahoma, U.S., with variable climatic conditions. The SMEI had strong correlations with existing SM indices. The correlations between SMEI and meteorological drought indices were 33% stronger on average when compared to existing SM indices. Higher correlation coefficients of SMEI were observed at all study sites, suggesting that SMEI performs well under variable climatic conditions experienced across Oklahoma. The correlation coefficients between existing SM and meteorological indices found in the present study were similar to those reported previously.
The SMEI could capture temporal variations in drought and provided better response to the variation in drought magnitude in Oklahoma when compared to other indices. This can be attributed to the inclusion of ET with soil moisture, which improved the drought sensitivity of the index. In addition, SMEI had a better performance capturing the spatial variations in drought when compared to other indices and the maps developed by the U.S. Drought Monitor.
In general, the correlation of drought indices with the yield of winter wheat gradually reduced when moving from the semi-arid Panhandle to more humid southeast. In comparison to other indices, however, the SMEI had the strongest correlations during the critical wheat growth stages in spring (r > 0.9). These results suggest that SMEI can be used more effectively to demonstrate the progress of agricultural drought under varying climates and its impacts on the crop yields. Farmers could also potentially optimize their decision making for the best use of their crops according to the predicted yield beforehand during the early crop growth stages.
The newly introduced index (SMEI) in this study was calculated using in-situ soil moisture measurements and weather data. Measured soil moisture data may not be available for all places and periods of interest. However, widely available remotely sensed and/or modeled soil moisture data can be used in calculating SMEI if measured data are not available. Future studies must analyze the effectiveness and the sensitivity of SMEI over longer periods and under different climatic conditions than those of this study.