The three LSMs (VIC, Noah, and CLM) used in this study were forced with a common meteorological dataset that comprises daily precipitation, maximum and minimum temperatures, and wind speed. We used a 0.25∘ daily gridded precipitation product available for the period 1901–2015 from the India Meteorological Department (hereafter IMD; Pai et al., 2015), which was developed by IMD using data from 6995 gauge stations across India and an inverse distance weighting scheme (Shepard, 1984). In the gridded precipitation data from IMD, orographic and topographic features of precipitation are well captured along with the spatial variability associated with the Indian summer monsoon. We used 1∘ daily gridded maximum and minimum air temperatures from IMD (Srivastava et al., 2009), which were developed using the data from 395 observation stations across India. We re-gridded air temperatures from 1 to 0.25∘ using the method described in Maurer et al. (2002), which is based on a temperature lapse rate of 6.5 ∘C km-1 rise in elevation and the SYMAP algorithm. In re-gridding of air temperature, we used a 0.25∘ digital elevation model (DEM) that was resampled from the original 30 m elevation data from the Shuttle Radar Topography Mission (SRTM). The gridded precipitation and air temperature products have been used in many previous studies on drought and heat waves (Shah et al., 2017; Shah and Mishra, 2014; R. D. Shah and Mishra, 2016, 2015; Mishra et al., 2016).

We used simulated soil moisture from the three LSMs, the VIC, the Noah, and the CLM, to assess uncertainty in root-zone (60 cm) depth soil moisture and retrospective drought assessment. These three LSMs have been widely used for producing land-surface fluxes at global and regional scales (Rodell et al., 2004; Shah and Mishra, 2015; Unnikrishnan et al., 2013; see also Table S1 in the Supplement for a brief description of major hydrological processes). All three LSMs were forced with the same meteorological forcing from IMD at 0.25∘ resolution, and additional (radiation-related) forcing variables for the Noah and CLM were derived from the MTCLIM algorithm integrated in the VIC model (see Bohn et al., 2013). This will keep basic forcing data consistent across models. In all LSMs routines are there which disaggregate daily precipitation uniformly to a sub-daily timescale, while temperature and radiation are temporally disaggregated following the diurnal cycle. Regarding the usage of land-surface datasets, we note here that each model used a slightly different set of input datasets. Soil textural properties and resulting hydrologic parameters (like field capacity, wilting points, and available water) are mostly derived based on the Harmonized World Soil Database (HWSD) and the Food and Agriculture Organization (FAO) based soil maps. However, since land-surface models use soil parameters from two different sources that can lead to differences in available water (Fig. S2), we evaluated the difference in the simulated soil moisture anomalies from the VIC model using the soil parameters based on the FAO and HWSD datasets (Fig. S3). We do not find a significant difference in the simulated soil moisture anomalies based on the soil parameters from the two sources. It is worth noting that the underlying (soil-textural) dataset for the HWSD is mostly derived based on the FAO dataset. The different products may use different pedo-transfer functions to derive soil-related parameters. Nevertheless the derived (static) soil parameters do not induce significant differences in the temporal dynamics of simulated soil moisture anomalies (Figs. S2 and S3).

Vegetation characteristics specified within LSMs were derived from the Advanced Very High-Resolution Radiometer (AVHRR) and Moderate Resolution Imaging Spectroradiometer (MODIS) datasets. We evaluated the sensitivity of vegetation parameters derived from AVHRR and MODIS to simulated soil moisture using the Noah model (Fig. S4). In this case too we do not find any substantial differences in the dynamics of the simulated soil moisture anomalies due to differences in the underlying vegetation parameters estimated from the AVHRR and MODIS datasets (Fig. S4). We note that the analysis conducted here is at a 0.25∘ spatial resolution – as a result the differences among the model simulations due to fine-scale soil and vegetation characteristics may not be observed at a coarse scale.

All three LSMs were first spun up using 65 years (1951–2015) of data to establish initial conditions for the modeled states and fluxes. All three LSMs were manually calibrated to match simulated monthly streamflow with observed streamflow data obtained from the India-WRIS (http://www.india-wris.gov.in) at the gauging stations (H. L. Shah and Mishra, 2016) that are least affected by human interventions related to water diversion and water withdrawal for irrigation (see the Supplement, Fig. S1 and Table S2, for the geographical location of these stations). We identified the calibration parameter for each LSM based on prior studies (Cai et al., 2014; Hogue et al., 2005; Nijssen et al., 2001) and by performing a simple (one parameter at a time) sensitivity analysis. We used soil thickness also as calibration parameters following the success of calibrating the VIC using soil-layer thickness (Nijssen et al., 2001; R. D. Shah and Mishra, 2016). The calibration parameters were manually adjusted so as to match observed streamflow (see Table S2). Further, we evaluated the model skill by comparing simulated soil moisture with station and satellite-based soil moisture, and also by comparing the total column soil moisture changes with terrestrial water change based on GRACE products (see Sect. 3.1 for more details).

We used the Standardized Precipitation Index (hereafter SPI; McKee et al., 1993) and the Standardized Soil moisture Index (SSI; Hao and AghaKouchak, 2013) to represent meteorological and soil moisture (agriculture) droughts, respectively. We used 60 cm soil depth as a representative of root-zone soil moisture (Shah and Mishra, 2015). Since depths of root-zone and soil layers are different in all three LSMs (Table S1), we estimated 60 cm soil moisture for each grid cell and for each LSM, separately. A parametric (Gamma) distribution was fitted to precipitation and root-zone soil moisture to estimate SPI and SSI, respectively. For both SPI and SSI, the cumulative distribution functions obtained by fitting the Gamma distribution were mapped onto the normal distribution functions to represent a dimensionless index and to derive drought indices (see Shah and Mishra, 2015, and Appendix A for more details). We note that there are other approaches for estimating the soil moisture drought index – for example a non-parametric percentile-based drought index (Samaniego et al., 2013), but in this study we used a parametric (Gamma) distribution for estimating SSI so as to be consistent with the precipitation-based drought index (SPI).

Intensity–areal extent–frequency (IAF) curves for drought events were constructed to understand the frequency and severity of droughts in India for the period 1951–2015. The IAF curves were estimated using the root-zone soil moisture from the three LSMs (i.e. VIC, Noah, and CLM). We estimated drought severity using the 4-month SSI at the end of the monsoon and Rabi seasons (so as to represent the entire season) for the whole of India and for the Indo-Gangetic Plain region (longitude 75–90∘ E and latitude 23–30∘ N; Fig. S1) considering the climatological period of 1951–2015. The method to construct IAF curves has been described in detail in Mishra et al. (2016) and Mishra and Cherkauer (2010). IAF curves were estimated using the following steps: (i) for each year, the mean 4-month SSI value was estimated for all the grids for areal extents of 2, 5, 10, 20, 30, 40, 50, 60, 70, 80, 90, and 100 %; (ii) for each threshold of areal extent, mean severity of root-zone soil moisture drought was estimated for each year during the 1951–2015 period; (iii) the Generalized Extreme Value (GEV) distribution was fitted to the mean severity for the selected areal extents and parameters (shape, location, and scale) were estimated using the maximum likelihood method; (iv) drought severity was estimated for the selected return periods of 2, 5, 10, 20, 25, 50, 100, 200, and 500 years for each areal-extent threshold to construct IAF. Using IAF curves, mean intensity of drought can be estimated or for a given areal extent and frequency of drought. We evaluated the goodness of fit of the GEV distribution using QQ plots and a Chi-square goodness of fit test (Figs. S15–S17 and Tables S11–S13).

The three land-surface model (VIC, Noah, and CLM) parameterizations were manually constrained (calibrated) against observed streamflow across a set of 18 major river basins covering approximately the entire landmass of India (see Fig. S1 and Table S2). The performance of the LSMs in capturing the temporal dynamics of monthly streamflow during calibration and validation periods is quite satisfactory for most of the river basins (Table S2). The median correlation r (and Nash–Sutcliffe efficiency, NS) values estimated across these basins during the calibration period are around 0.91 (0.78), 0.90 (0.70), and 0.90 (0.70) for the VIC, Noah, and CLM, respectively. A similar level of (median) skill is also observed during the validation period (Table S2). The skill of the multi-model averaged streamflow of the three LSMs is comparatively better than that of individual models – with an overall median r (and NS) value estimated across all basins of 0.91 (0.80) and 0.94 (0.77) during the calibration and validation periods, respectively. The ensemble mean of simulated streamflow from the three LSMs performed better against the observations. We also notice a relatively poor skill for all three LSMs and the ensemble mean in the coastal basins (e.g., Cauvery and eastern coast basins – Table S2), which could be attributed to a number of factors including errors in forcing data and model parameterizations. Nevertheless, considering the wide range of hydro-climatic gradients across India, the efficiency of the three LSMs for capturing the observed streamflow can be considered reasonable.

Next we evaluated the skill of each model in capturing the observed dynamics of near-surface and 60 cm soil moisture (Figs. 1 and S3–S4). We used three different sources of soil moisture observations for this comparison purpose. The first set consisted of the weekly soil moisture observations taken at 18 IMD-based stations during the monsoon (JJAS) season for the period 2009–2013 (Unnikrishnan et al., 2013). The model simulated 60 cm soil moisture dynamics were compared against observations, which generally revealed a good skill for all three models (Fig. 1). Model simulated soil moisture showed a relatively higher correlation with observations in the northern and western regions as compared to those located in the southern coastal belt. Among models, the Noah simulated soil moisture exhibited higher correlation as compared to other two LSMs. For this setup, we also compared the calibrated vs. uncalibrated model runs to understand what improvements (if any) could be achieved by the parameter calibration of simulation of soil moisture anomalies. We find limited benefits of the model calibration in this case – only the VIC model benefited from the model calibration, mainly in the northern region locations and a few of the southern locations.

The second set of evaluations considered the continuous soil moisture observation datasets at an IIT Kanpur site available from the International Soil Moisture Network (ISMN: Dorigo et al., 2011). Although all three models exhibited a general bias in capturing absolute values of the observed soil moisture, their daily variability observed over the course of the year is well captured by all three models (Fig. S5). Moreover, the Noah and CLM models show improvements in terms of reducing overall bias as a result of model calibration.

Finally, our third set of model evaluations considered an assessment of the model skill for capturing the remote-sensing-based soil moisture available from the ESA-CCI product (Dorigo et al., 2012). Here we used the modeled top-layer (10–30 cm) annual soil moisture over the period 1979–2012 for the comparison; however, the ESA-CCI soil moisture product is more reliable after 2001 (Dorigo et al., 2017) (Fig. S6). Despite the limitation that the ESA-CCI soil moisture inference is for the top few centimeters of the earth's surface, we find a positive correlation with modeled soil moisture for all three models across a large part of India. A relatively higher correlation (more than 0.6) can be noticed for regions in the northwestern and southern peninsular parts of India.

We also evaluated the skills of LSMs for the terrestrial water storage (TWS) anomalies from the Gravity Recovery and Climate Experiment (GRACE – release v5.0; Landerer and Swenson, 2012) derived products (1∘ × 1∘) for the period 2002–2015. We used the ensemble mean of three available GRACE-TWS products (GeoForschungsZentrum, GFZ, Potsdam, Germany, Centre for Space Research at the University of Texas at Austin, USA, and Jet Propulsion Laboratory, USA) to reduce the noise (and scatter) among different TWS products. We compared ensemble mean GRACE-TWS against the monthly anomalies of modeled total column soil moisture from each of the three LSMs and their ensemble mean. The modeled total column soil moisture was aggregated to 1∘ spatial resolution to match the (coarse) resolution of the GRACE-TWS product. Overall, all three LSMs are able to capture the temporal dynamics of GRACE-TWS anomalies reasonably well across a large part of India (Fig. S7). The median correlation estimated across the modeled grid cells is more than 0.6 for all three LSMs – and the ensemble mean of simulated total column soil moisture anomalies showed an overall best (median) skill. All three LSMs (and the ensemble mean) exhibited a systematically lower performance in the northwestern part of India (Fig. S7), which is most probably related to groundwater pumping effects that are not modeled in LSMs but are captured in GRACE datasets (Asoka et al., 2017). Each LSM shows a slightly different area (grid cells) with the best skill score that motivates the use of a multi-model ensemble mean to capture the (GRACE-based) water storage anomalies across a large part of India.

We estimated areal extent of severe to exceptional droughts (SPI < -1.3) based on the 4-month SPI at the end of the monsoon season (representing accumulated precipitation for the June to September period) for the period of 1951–2015 (Fig. 2a). The top five monsoon season (JJAS) drought events occurred in 1987 (areal extent of severe to exceptional droughts: 35 %), 2002 (33.5 %), 1979 (27.7 %), 1972 (26.3 %), and 2009 (24.6 %) at the all-India level. The monsoon season droughts of 2015 (with an areal extent of 17.4 %) and 2014 (14.4 %) ranked 8th and 10th during the period of 1951–2015. Mishra et al. (2016) reported that the 2014–2015 monsoon season drought in the Indo-Gangetic Plain was the most severe during the history of 116 years, with a return period of 542 years.

We find that uncertainty in simulated areal extent of soil moisture drought (4-month SSI at the end of the monsoon season) estimated based on three LSMs is moderate during the monsoon season and has a year-to-year variability (Fig. 2a). For example, the 1 standard deviation representing the uncertainty in the simulated areal drought extent is on average estimated to be around 1.6 %. However, during the 1972 and 1979 monsoon seasons, uncertainty in areal extent of drought is 8 and 14 %, respectively (Fig. 2a; see also Table S3). We estimated the ensemble mean (ENS-SSI) areal extent of 60 cm soil moisture drought from the three LSMs for the monsoon season and found that 1987 (36.7 %), 2002 (35.9 %), 2009 (31.5 %), 1972 (29.8 %), and 1965 (23.4 %) are the top 5 drought years during the period of 1951–2015 (Fig. 2a). We notice that 1979 (areal extent: 20 %) ranked 6th on record, while 2015 (16.5 %) and 2014 (7.13 %) ranked 8th and 15th, respectively, based on ENS-SSI areal extent of the simulated 60 cm soil moisture drought during the monsoon season.

Since the Rabi season (NDJF) is the key crop-growing season in India, we estimated areal extents of 60 cm soil moisture drought from the three LSMs for the period of 1951–2015 (Fig. 2b). While year-to-year variability in the uncertainty of areal extent of droughts was found, the uncertainty in 60 cm soil moisture based on the areal extent of droughts was substantially higher (∼ 5 %) during the last decade of 1951–2015, which might be associated with frequent drought events during this period. We found a high uncertainty (∼ 7 %) in the areal extent of soil moisture drought in 2011 (Fig. 2b). Based on the ENS-SSI areal extent of 60 cm soil moisture drought, the top 5 drought years were 2003 (areal extent 30 %, uncertainty 9.6 %), 2001 (27.6, 12.0 %), 1966 (22.6, 4.7 %), 1973 (20.7, 8.9 %), and 1988 (20.6, 6.1 %). Droughts in the Rabi season can be driven by both monsoon season precipitation deficit and positive temperature anomalies during this season. Uncertainty in the top five drought events in the Rabi season was substantial (5–12 %), which underscores the need for multi-model drought assessment in the growing season.

The VIC, Noah, and CLM show an all-India median autocorrelation of 0.23, 0.37, and 0.71, respectively, at a 4-month lag (Fig. 3), indicating that the CLM has the highest persistence in the 60 cm soil moisture. Spatial differences in soil moisture persistence were also observed (Fig. 3). Our results are in agreement with the findings of Wang et al. (2009), who reported higher persistence for CLM modeled soil moisture, which can be attributed to its higher water holding capacity and thicker soil column (Fig. S2). Figure S2 shows that though there is not much difference in soil texture provided as input to three LSMs, the available water in the total column is much higher for the CLM. We find that soil layer thickness is strongly related to soil moisture persistence regardless of model calibration (Fig. S8). The same thing can be noted even considering the first three soil layers which cover a 60 cm soil column for all three LSMs, the soil-layer thickness of the CLM dominates over another two LSMs (Fig. S9).

Furthermore we find that all-India averaged mean monthly soil moisture is highest during July in the VIC and Noah LSMs, which is consistent with the seasonal cycle of all-India averaged precipitation (Fig. 3a–c). However, all-India averaged mean monthly 60 cm soil moisture reaches the highest level in the month of August in the CLM model (Fig. 3a–c). The 1-month lag between peak precipitation and peak 60 cm soil moisture from the CLM can be due to a relatively deeper soil column (Fig. S8) and higher total column water holding capacity, a higher number of soil layers, and differences in processes related to soil hydrology as discussed in Wang et al. (2009) and Xia et al. (2012).

To understand the relationship between meteorological and agricultural droughts, lagged correlation analysis was performed between the 4-month SPI at the end of the monsoon season and the 4-month SSI (at the end of JJAS, JASO, ASON, and so on). We find that the 4-month SSI at the end of the monsoon season for the VIC model showed the highest correlation with the 4-month SPI (JJAS), while the 4-month SSI from the Noah model showed the lowest correlation (Fig. S10a). These results indicate that the 60 cm soil moisture from the VIC model responds faster to the monsoon season precipitation than the other two LSMs, which can be associated with soil moisture persistence and model parameterization (Van Loon et al., 2012; Wang et al., 2009; Xia et al., 2012). However, we notice that the correlation between the 4-month SPI at the end of the monsoon season and the 4-month SSI declines rapidly after October (ONDJ, NDJF, and so on) for the VIC and Noah models (Fig. S10a). On the other hand, the CLM shows substantially higher persistence even for the March–June 60 cm soil moisture, which can be attributed to a deeper soil column and differences in the other processes related to soil hydrology. These results also indicate that the anomalous precipitation during the monsoon season can last longer and have a substantial influence on agriculture drought estimated using the 60 cm soil moisture from the CLM as reflected by the strength of the relationship between the 12-month SPI and 12-month SSI in the CLM (Fig. S10b).

Areal extent and severity of agricultural droughts estimated using the 60 cm soil moisture show a considerable uncertainty mainly due to the differences in soil moisture persistence characteristics among three LSMs (Fig. 3). We estimated areal extent of agriculture drought from the three LSMs considering the period that showed maximum correlation against the monsoon season precipitation (Fig. S10c). For instance, 4-month SSI at the end of October (JASO) showed the highest correlation with 4-month SPI at the end of September (JJAS) for the VIC and Noah models (Fig. S10a). On the other hand, the 4-month SSI at the end of November (ASON) showed the highest correlation with the 4-month SPI at the end of September (JJAS) for the CLM. Therefore, we considered root-zone soil moisture for JASO, JASO, and ASON periods from the VIC, Noah, and CLM, respectively (Fig. S10c) to understand the response of the monsoon season deficit in precipitation on agricultural drought. We find that the uncertainty in the areal extent of agricultural drought is substantially reduced considering the lagged response of the monsoon season precipitation and soil moisture (Fig. S10c and Table S3), indicating that the major source of uncertainty in areal extent of agricultural droughts is soil moisture persistence in the LSMs.

We reconstructed major monsoon season drought events over India using 60 cm soil moisture from the three LSMs for the period of 1951–2015 (Fig. 4). The meteorological and agricultural droughts were represented using the 4-month SPI and 4-month SSI, respectively, at the end of the monsoon season. We estimated the ensemble mean 4-month SSI from the three LSMs to understand whether the individual LSMs show a larger difference from the ensemble mean. Based on the 4-month SPI at the end of the monsoon season, we selected the top two most widespread drought events that occurred in 1987 and 2002 (Figs. 2a and 4). Moreover, we also selected a recent drought event (2015) that caused an enormous water crisis in the Indo-Gangetic Plain (Mishra et al., 2016). For all three major droughts (1987, 2002, and 2015) in the monsoon season, we compared areal extents of 60 cm soil moisture drought in the monsoon season estimated from the three LSMs.

Areal extents of the monsoon season droughts (meteorological and agricultural) in 1987, 2002, and 2015 show that droughts were mainly caused by the monsoon season precipitation deficits (Fig. 4). We notice positive air temperature anomalies in all 3 years (1987, 2002, and 2015); however, patterns of agricultural and meteorological droughts were largely similar (Fig. 4). Among the three drought events, the monsoon season air temperature anomaly (positive) was strongest in the 2015 monsoon season. Uncertainty in areal extent of agricultural droughts estimated using the 60 cm soil moisture is presented in supplemental Table S3. We notice large uncertainty in the areal extent of drought simulated from the three LSMs for the 2015 event (Table S3). The VIC model simulated areal extent of soil moisture drought was 14 % during the 2015 monsoon season, while the areal extent of drought simulated from Noah and the CLM was 21.2 and 18.1 %, respectively (Table S3).

Similar to the monsoon season droughts, we compared the spatial pattern of droughts simulated by the three LSMs for major droughts in the Rabi season, which occurred in 1966, 1973, 2001, and 2003 (Fig. S11). We notice that major droughts in the Rabi season were also largely driven by the precipitation deficit, and the role of positive air temperature anomalies was relatively minor (Fig. S11). Overall, the VIC model shows the lesser intensity of drought during the Rabi season as compared to the Noah and CLM for all the years, which can be attributed to differences in soil moisture persistence in the three LSMs (Fig. 3a). We find higher uncertainty in areal extent of drought during the 2001 Rabi season (Fig. S11) than other years of major drought events, which may be due to higher impacts of air temperature on drought during 2001. The overall uncertainty in the areal extent of agricultural droughts estimated using the 60 cm soil moisture estimated from the three models is presented in supplemental Table S4. In 2001, the areal extent of droughts simulated by the VIC, Noah, CLM, and their ensemble mean were 17.2, 40.7, 24.7, and 26.1 %, respectively (Table S4). Large differences in the areal extent of 60 cm soil moisture drought simulated from the three LSMs were also noted for the years 1966, 1973, 2001, and 2003 (Table S4).

Uncertainty in the multi-model drought estimates was evaluated using IAF curves (Fig. 5). Drought intensity associated with the selected areal extent thresholds was estimated for the return periods of 10, 20, 50, 100, 200, and 500 years (Fig. 5). We find that multi-model-based drought intensity for a selected areal extent has a much larger uncertainty when the 95 % confidence interval of the GEV parameters was considered (Table S6). However, uncertainty in drought intensity was lower when only mean values of the GEV parameters for drought intensity and areal extents were considered. Uncertainty in drought intensity appeared to grow with an increase in the return period (Fig. 5). Based on the IAF curves, we find that the 2002 monsoon season drought has a return period of 50 years. Moreover, uncertainty in drought intensity from the three LSMs was larger for smaller areal extents (Fig. 5). For instance, a drought of 50 % areal extent and 50-year return period has intensities of -1.66, -1.91, and -1.53 simulated from the VIC, Noah, and CLM, respectively (Table S5). On the other hand, a drought of 5 % areal extent and 50-year return period has intensities of -2.89, -3.79, and -3.09 simulated from the VIC, Noah, and CLM, respectively (Table S5). Most of the return periods and areal-extent drought intensities were higher for the Noah model and lower for the VIC model (Table S5), which can be associated with the differences in 60 cm soil moisture persistence in the three LSMs.

A considerably higher uncertainty in IAF curves during the Rabi season was noticed (Fig. S12; Tables S7 and S8). For instance, a drought of 50 % areal extent and a 50-year return period can have intensities of -1.23, -1.62, and -1.50 for 60 cm soil moisture obtained from the VIC, Noah, and CLM, respectively (Fig. S12 and Table S7). A drought of 5 % areal extent and the 50-year return period have intensities of -2.17, -3.71, and -3.25 simulated from the VIC, Noah, and CLM, respectively (Table S7). Higher uncertainty in drought intensities during the Rabi season can be attributed to the response of soil moisture in the three LSMs to meteorological forcing (precipitation and air temperature).

As the Indo-Gangetic Plain is one of the most intensive crop-growing regions in the world, we evaluated the uncertainty in the IAF curves constructed using the 60 cm soil moisture from the three LSMs (Fig. S13). We used the 12-month SSI at the end of December based on the areal averaged mean annual 60 cm soil moisture over the Indo-Gangetic Plain. Similar to IAF curves for the all-India averaged 60 cm soil moisture, a large uncertainty was found due to differences in the GEV parameters (Fig. S13, Tables S9 and S10). Moreover, uncertainty in IAF curves of the Indo-Gangetic Plain increases with the return period and declines with an increase in areal extent of droughts. For instance, for an aerial extent of 5 % and a return period of 50 years, drought intensities were -2.59, -3.06, and -2.48 for the VIC, Noah, and CLM (Table S9). However, when areal extent increases to 50 %, drought intensities are -1.42, -1.49, and -1.45 for the VIC, Noah, and CLM, respectively (Table S9, Fig. S11). Overall, uncertainty in drought intensity is higher for localized droughts that have higher return periods. These results further indicate that 60 cm soil moisture drought characteristics can have large uncertainty arising from different LSMs, which can be reduced by considering the multi-model ensemble agricultural drought assessments.

We evaluated the differences in the coupling of 60 cm soil moisture (SSI) with monsoon season precipitation and air temperature (Figs. 6 and S14). We find that a 4-month SSI at the end of the monsoon season from the VIC model is strongly coupled (correlation coefficient = 0.90) with the monsoon season precipitation over India. On the other hand, the 60 cm SSI showed correlation coefficients of 0.79 and 0.74 against the monsoon season precipitation for the Noah and CLM, respectively (Fig. 6c and e). These results show differences in the response of 60 cm SSI against changes in the monsoon season precipitation. However, the all-India averaged 60 cm SSI showed stronger coupling with the monsoon season air temperature for the Noah and CLM models (correlation = -0.65 and -0.67) than that of the VIC model (correlation = -0.53) (Fig. 6b, d and f). Interestingly, the coupling between the 60 cm SSI for the Rabi and monsoon season precipitation is stronger for the CLM (correlation = 0.76) and Noah (correlation = 0.66) than that of the VIC model (correlation = 0.55). These results indicate that the monsoon season precipitation deficit can have a larger influence on the Rabi season drought in the CLM and Noah models compared to that of the VIC model. Similarly, the Rabi season air temperature showed a stronger relationship with the 60 cm SSI for the CLM (correlation = -0.51) and Noah (-0.37) than that of the VIC model (-0.31). Similar differences in the Rabi and monsoon season 60 cm SSI for the Indo-Gangetic Plain were observed for the three LSMs (Fig. S14), indicating that the drought indices based on the 60 cm soil moisture may show different sensitivity to atmospheric forcing, which can show a substantial variation across regions and seasons.