1. Introduction
Impact modeling, such as agricultural modeling and monitoring activities, is heavily based on the availability of daily meteorological variables, such as precipitation and mean temperature. In regions with sparse meteorological stations and/or with restricted data accessibility and availability, reanalysis products offer a convenient way to work. It is therefore essential to evaluate their behavior as well as their intrinsic uncertainties. Here, we focus on precipitation in the core monsoon region of Asia (Noska and Misra 2016; Ding and Chan 2005), a key agricultural region of the world. This region is characterized by high population and intensive agricultural land use. Precipitation variability within the monsoon season can cause adverse impacts on agricultural yields (Singh et al. 2014, and references therein). For example, prolonged dry spells during the Indian summer monsoon season can substantially reduce yields of Kharif (monsoon) crops when these events occur during soil preparation, transplanting, or at critical crop growth phases. The Indian summer monsoon is responsible for 85% of the annual precipitation (Singh et al. 2014) and it is vital for Indian agriculture. Sensitivity of agricultural production to rainfall dynamics is generally higher over rain-fed areas (e.g., Ceglar et al. 2016). Although irrigation is intensively practiced in monsoonal Asia (Hatcho et al. 2010), irrigation water availability still depends on precipitation and could become an issue in a changing climate context (e.g., Elliott et al. 2014). Monitoring droughts during the monsoon season is therefore necessary to manage natural resources and it has high importance for policy makers (Shah and Mishra 2014).
Several evaluation studies have been performed on monsoon precipitation by using different observational and reanalysis datasets. Recently, Shah and Mishra (2014) reported positive biases in the Indian monsoon precipitation affecting the ECMWF interim reanalysis (ERA-Interim, ERA-I; Dee et al. 2011), the NASA reanalysis MERRA (Rienecker et al. 2011), and the NCEP reanalysis CFSR (Saha et al. 2010). Huang et al. (2016) reported overall good performance of reanalyses in reproducing the climatology and the variability of the East Asia summer monsoon precipitation. Rana et al. (2015) analyzed reanalyses (CFSR and ERA-I) and rain gauge– and satellite-based datasets; reanalyses reproduced the seasonality in the Indian subcontinent well, especially in low-elevation areas, while more uncertainties remained in high-elevation areas and areas with poor weather station coverage. It is also worthwhile highlighting that, as precipitation is a diagnostic variable of reanalyses, model physics and simulated circulation play an important role (Bosilovich et al. 2008).
As for agricultural monitoring and modeling, precipitation characteristics at different time scales (from daily to seasonal) are relevant and must be evaluated. As discussed by Huang et al. (2016), there are large uncertainties in nonrainfall and heavy rainfall classes among different reanalyses in eastern Asia. Therefore, this study aims at evaluating four latest-generation reanalysis datasets (some of them offering bias-corrected precipitation data) in terms of various precipitation characteristics at different time scales. Seasonality, distributional characteristics of daily precipitation, precipitation extremes, and drought during the monsoon season are analyzed by using recently proposed methods and approaches.
2. Data and methods
Four reanalysis datasets and two observational gridded datasets (one rain gauge based and another one satellite and rain gauge based) have been used for this study (Table 1). These datasets have been selected due to their spatial resolution and their spatiotemporal coverage. The analyzed area extends over 0°–40°N, 65°–125°E, covering the core Asian monsoon areas, except the northeastern part of the region affected by the East Asian summer monsoon (Ding and Chan 2005). To perform a quantitative comparison, all datasets have been regridded on the same 0.75° reference grid (i.e., the one of ERA-I) by applying a conservative remapping procedure (Chen and Knutson 2008, and references therein). Among the chosen reanalyses, ERA-Interim/Land (the land version of ERA-I, hereafter ERA-I-Land; Balsamo et al. 2015) and AgMERRA (the agricultural version of MERRA; Ruane et al. 2015) provide bias-corrected precipitation. ERA-I-Land is a single 32-yr simulation with the latest ECMWF land surface model driven by meteorological forcing from ERA-I (Dee et al. 2011). It is characterized by a monthly bias adjustment of precipitation based on the Global Precipitation Climatology Project (GPCP) monthly precipitation dataset (Adler et al. 2003). AgMERRA (Ruane et al. 2015) is a NASA reanalysis based on MERRA and MERRA-Land reanalysis (Rienecker et al. 2011; Reichle et al. 2011). It includes a monthly bias adjustment of precipitation based on the Climate Research Unit (CRU) wet days dataset (Harris et al. 2014), the Global Precipitation Climatology Centre (GPCC) dataset (Schneider et al. 2014), the University of Delaware precipitation dataset (WM; Willmott and Matsuura 1995), and satellite-based products [TRMM (Huffman et al. 2007), CMORPH (Joyce et al. 2004) and PERSIANN (Hsu and Gao 1997)]. The other two reanalyses, ERA-I (Dee et al. 2011) and JRA-55 (Kobayashi et al. 2015), are not bias-corrected. Concerning the gridded observations, APHRODITE (Yatagai et al. 2012) and Climate Hazards Group Infrared Precipitation with Stations (CHIRPS; Funk et al. 2014) have been retrieved. APHRODITE is a daily precipitation dataset based on rain gauge data developed by a consortium of the Research Institute for Humanity and Nature and the Meteorological Research Institute of the Japan Meteorological Agency. CHIRPS is a quasi-global satellite/rain gauge dataset developed by the Climate Hazards Group.
Precipitation datasets.
As the temporal availability of the datasets is not the same, the common period from 1981 to 2007 has been selected for the analysis. Moreover, the period from June to August (JJA) has been chosen as the monsoon period. The onset of the monsoon depends on both the year and the region, but JJA can be thought as a conservative choice to avoid the complexity and the uncertainties of dealing with the spatiotemporal nonstationarity.
a. Seasonality





b. Distributional analysis



c. Precipitation events




d. Drought
To characterize droughts, the Standardized Precipitation Index (SPI; Mckee et al. 1993) is here applied. SPI is a dimensionless index based on the probability distribution of precipitation. It can be generated for different time scales, such as monthly and seasonal. Drought events are categorized as follows: extreme (SPI less than −2), severe (SPI between −2 and −1.5), and moderate (SPI between −1.5 and −1.0). For the assessment of monsoon season droughts, the summer season SPI is used here (i.e., SPI on three months’ time scale for August).
3. Results
a. Seasonality
As observational gridded datasets also suffer from uncertainties and can fail in reproducing precipitation field over areas with poor data coverage, it is worthwhile to compare APHRODITE and CHIRPS. In terms of seasonality, APHRODITE shows a more uniformly distributed precipitation through the year in western China w.r.t. CHIRPS (Fig. 1). In northwestern India/Pakistan, CHIRPS shows a broader area with precipitation being concentrated in a few months (Fig. 1). Concerning reanalyses, all of them show more uniformly distributed monthly precipitations when compared to CHIRPS (Fig. 2). However, differences can be observed in the dimensionless seasonality index S (Fig. 3) as, for instance, a broader and more pronounced monsoonal behavior in Myanmar in ERA-I (Fig. 3). When compared to APHRODITE, JRA-55 and ERA-I overall show a similar spatial pattern in terms of seasonality with precipitation being more uniformly distributed during the year, except in some areas such as northwestern India/Pakistan and western China (Fig. 2). However, ERA-I is characterized by a broader and more pronounced monsoonal behavior in the northeastern side of the Bay of Bengal and along the southern side of the Himalaya, while JRA-55 overall shows a less pronounced monsoonal behavior (Fig. 3). Both AgMERRA and ERA-I-Land show less pronounced differences (Fig. 2), although they both have higher monsoonal activity over Myanmar.
Relative entropy (DX) of (a) APHRODITE and (b) CHIRPS and (c) their difference, and the dimensionless seasonality index (DSI) of (d) APHRODITE and (e) CHIRPS and (f) their difference.
Citation: Journal of Climate 30, 2; 10.1175/JCLI-D-16-0227.1
Difference in relative entropy between the four different reanalyses and (top) APHRODITE and (bottom) CHIRPS observational datasets.
Citation: Journal of Climate 30, 2; 10.1175/JCLI-D-16-0227.1
Difference in dimensionless seasonality index between the four different reanalyses and (top) APHRODITE and (bottom) CHIRPS observational datasets.
Citation: Journal of Climate 30, 2; 10.1175/JCLI-D-16-0227.1
b. JJA differences
During JJA, CHIRPS shows higher accumulated precipitation (especially in the area under the influence of both the Indian summer monsoon and the west Pacific summer monsoon) except over the northern part of the domain and along the western coast of India (Fig. 4). This behavior is almost stationary during the investigated period and differences are always more pronounced in the zonal band comprised between 10° and 15°N (Fig. 4). As for reanalyses, AgMERRA is generally characterized by the smallest differences in summer mean precipitation cumulates, except over Myanmar and along the southeastern side of the Tibetan Plateau where it significantly overestimates JJA precipitation (w.r.t. APHRODITE; Fig. 5). Large significant positive biases affect ERA-I especially in the eastern part of the domain (Fig. 5) whereas these biases are considerably reduced in ERA-I-Land, especially over southeastern China. JRA-55 has a completely different behavior w.r.t. the other reanalyses, being largely affected by significant underestimation of JJA accumulated precipitation (Fig. 5). These differences contribute to explain and better understand the identified behavior in the dimensionless seasonality index discussed in the previous subsection. In terms of spatial correlation of JJA precipitation, AgMERRA achieves the best performance, being characterized by correlation values higher than 0.9 with both APHRODITE and CHIRPS (Fig. 6). JRA-55 and ERA-I-Land have similar performance (with correlation values varying around 0.8), but ERA-I-Land is affected by a sudden drop in 1990. ERA-I is the worst in terms of spatial correlation with values always lower than 0.8 (Fig. 6). Concerning the temporal evolution of the identified JJA-biases, it is worthwhile to highlight the sudden decrease in the positive bias affecting ERA-I (w.r.t. APHRODITE) in the mid-1990s over the zonal band comprised between 10° and 15°N (Fig. 7) as well as the changes affecting ERA-I-Land. The other two reanalyses do not show remarkable temporal changes (Fig. 7).
Mean summer precipitation (mm) of (a) APHRODITE and (b) CHIRPS. (c) Differences in the mean summer accumulated precipitation of APHRODITE and CHIRPS datasets. Regions where differences are statistically significant (p < 0.05) are denoted with white crosses. (d) Mean zonal differences in the summer accumulated precipitation (mm) of APHRODITE and CHIRPS from 1981 to 2007.
Citation: Journal of Climate 30, 2; 10.1175/JCLI-D-16-0227.1
Differences in the mean summer accumulated precipitation of four different reanalyses and (top) APHRODITE and (bottom) CHIRPS observational datasets. Regions where differences are statistically significant (p < 0.05) are denoted with white crosses.
Citation: Journal of Climate 30, 2; 10.1175/JCLI-D-16-0227.1
Spatial correlations between summer accumulated precipitation of the different reanalyses and the two observational datasets, calculated for each year separately: (left) APHRODITE and (right) CHIRPS.
Citation: Journal of Climate 30, 2; 10.1175/JCLI-D-16-0227.1
Mean zonal differences in the summer accumulated precipitation of the four reanalyses w.r.t. (top) APHRODITE and (bottom) CHIRPS. Values (in mm) are presented for period between 1981 and 2007.
Citation: Journal of Climate 30, 2; 10.1175/JCLI-D-16-0227.1
Besides the identified difference in the JJA precipitation of the two observational datasets, the (rescaled) distributions of APHRODITE and CHIRPS are still significantly different (although the differences are not too large), with the former being characterized by fatter tails (not shown). Concerning the reanalyses, they all show more complex differences w.r.t. observations than a simple bias in the mean. Indeed, they all tend to underestimate the (rescaled) distributions and ERA-I shows the largest differences (except over the Tibetan Plateau; Fig. 8).
Distributional comparison of daily rescaled summer precipitation of the four different reanalyses w.r.t. (top) APHRODITE and (bottom) CHIRPS. Shown is the statistic A [Eq. (3)]. Higher values are associated with large differences in the distributions. White areas denote nonsignificant differences.
Citation: Journal of Climate 30, 2; 10.1175/JCLI-D-16-0227.1
c. Spatial and temporal variation in the skill of daily precipitation from reanalyses
Figure 9 shows the FSS as a function of the spatial scale for the two selected thresholds. The skill increases with increasing size of spatial scale, since rainfall events are more accurately captured over larger areas and models often tend to misplace localized intense rainfall events (Roberts and Lean 2008). Moreover, higher skill can be observed for lower precipitation thresholds (here, 75th percentile), including less intense and spatially broader events that are more accurately captured by reanalyses. Higher skill at the model grid scale can be observed for AgMERRA; however, for the 90th percentile the skill is generally below 0.5. With respect to APHRODITE, the FSS becomes comparable between ERA-I, JRA-55, and AgMERRA at spatial scales higher than 2.25°. When using CHIRPS as the reference dataset, AgMERRA and JRA-55 show similar skill throughout all spatial scales, whereas slightly lower performance is identified for ERA-I. The skill of ERA-I-Land is below the skill of all the other reanalyses, most probably due to sudden change after 1990, as described in the previous subsection. It should be emphasized that before that year the skill of ERA-I-Land is higher than the ones estimated for both ERA-I and JRA-55, especially when considering the 75th percentile threshold.
The distribution of FSS as a function of spatial scale for (left) 75th and (right) 90th percentiles chosen as precipitation thresholds. The box plot at each of the spatial scales represents the temporal distribution of FSS (i.e., the skill is assessed for each day of all summer seasons).
Citation: Journal of Climate 30, 2; 10.1175/JCLI-D-16-0227.1
Results suggest that skillful spatial scale is temporally varying. A scale of 2.25° appears to be a skillful scale for AgMERRA using the 75th percentile threshold. For ERA-I and JRA-55, the skillful spatial scale varies between 2.25° between 1990 and 2000 and 3.75° during the rest of the analyzed period. Considering the 90th percentile, the skillful spatial scale of 3.75° prevails in ERA-I, JRA-55, and AgMERRA; ERA-I-Land has a comparable skillful spatial scale before 1989, but it drops substantially after that year.
d. Drought
In terms of interannual drought variability, APHRODITE and CHIRPS tend to agree except over Myanmar, Cambodia, and the northwestern part of the domain (not shown). Consistency between drought-affected areas can be found for events such as the one in 1987 (northwestern India, Pakistan, and the Tibetan Plateau), 1997 (central India, central-western part of China, and parts of southeastern Asia), and 2002 (northwestern India, Pakistan, central China, and the North China Plain) (e.g., Zou et al. 2005; Wang 2006). The drought events shown by APHRODITE in 1999 (North China Plain, northwestern India, and Indonesia) and 2006 (central and southern India, Myanmar, and southern China) are not well reproduced by CHIRPS. Figure 10 shows the time series of the area affected by extreme, severe, and moderate drought (as defined in the previous section). Overall, drought-affected areas are well represented by all the reanalyses. However, higher differences from observations after 1990 can be observed in the time series of ERA-I-Land (especially for moderate events).
Difference in the area (expressed as the proportion of grid cells over land) affected by extreme (SPI less than −2), severe (SPI between −2 and −1.5), and moderate (SPI between −1.5 and −1) drought events, as represented by the reanalyses, with respect to the two observational datasets: (left) APHRODITE and (right) CHIRPS.
Citation: Journal of Climate 30, 2; 10.1175/JCLI-D-16-0227.1
Taylor diagrams shown in Fig. 11 represent the estimated differences in the spatial patterns of drought events. AgMERRA clearly outperforms the other reanalyses. Its spatial correlation with the observational datasets generally exceeds 0.8. Lower performance can be seen in Fig. 11 associated with ERA-I and JRA-55, which are characterized by similar behavior, while worse spatial representation is observed for ERA-I-Land. In this case, two well-separated groups of years can be identified. Spatial patterns tend to be better represented before 1990 (when ERA-I-Land outperforms ERA-I). Afterward, the spatial representativeness of drought patterns tends to be poor, with nonsignificant correlations and RMS, when compared to both observational datasets. The Taylor skill score S exhibits high interannual variability (not shown) but generally highly resembles the spatial correlation patterns of Fig. 11. A substantial drop of the skill score characterizes ERA-I-Land after 1990. ERA-I and JRA reach lower skill than AgMERRA, varying respectively between 0.2 and 0.4 and between 0.4 and 0.7. The skill score of ERA-I-Land drops from 0.4 to 0.1 after 1990. Generally, the skill drops in AgMERRA, ERA-I, and JRA after 1999.
Taylor diagrams for summer SPI over monsoon Asia, using (left) APHRODITE and (right) CHIRPS as reference datasets. Each point represents a specific year and colors are associated with the different reanalyses.
Citation: Journal of Climate 30, 2; 10.1175/JCLI-D-16-0227.1
To better understand and characterize the different drought representation, a specific analysis of severe drought events over selected regions of monsoon Asia is performed. Figure 12 shows areal averages of SPI over central India (region 1), Southeast Asia (region 2), and the North China Plain (region 3). Severe drought events occurred during the analyzed period in central India in 1987 and 2002 according to the observational datasets. The 1987 drought event is well represented by all the reanalysis datasets. On the contrary, the 2002 event is underestimated by ERA-I and ERA-I-Land. Reanalyses and CHIRPS also point to a drought event in 1992, whereas no drought occurrence is shown by JRA-55 and APHRODITE. Southeast Asia is characterized by a highly variable drought signal among different datasets. The most notable is the underestimation of SPI values in ERA-I and simultaneous overestimation of SPI in ERA-I-Land after 2000. Better agreement among different datasets can be observed over the North China Plain. An exception is ERA-I-Land, where the drought signal is not well reproduced after 2000. The most intensive drought events after 1990 (1997 and 2002) are well reproduced by AgMERRA and JRA-55, whereas ERA-I and ERA-I-Land fail to capture those events.
Time series of summer SPI over central India (region 1), Southeast Asia (region 2), and the North China Plain (region 3). Colors are associated with the different datasets.
Citation: Journal of Climate 30, 2; 10.1175/JCLI-D-16-0227.1
4. Discussion and conclusions
The comparison between the two observational precipitation datasets has shown remarkable differences (e.g., in the represented seasonality of precipitation as well as in the distribution of JJA daily precipitation). This highlights the uncertainties affecting all gridded observational datasets associated both with the methodological interpolation procedures and data spatiotemporal homogeneity and availability. Significant differences have been also identified in the comparison of all the reanalyses with the two observational dataset. An interesting sudden change in the 1990s seems to affect ERA-I in a specific zonal band (between 10° and 15°N) and ERA-I-Land (affected by another change at the beginning of the twenty-first century). In terms of droughts, ERA-I-Land fails to reproduce drought occurrence patterns after 1990, although its performance is higher than JRA-55 and ERA-I before that year. In terms of spatial correlation of JJA-accumulated precipitation, a drop is clearly affecting ERA-I-Land in 1990 (Fig. 6). However, the spatial correlation of seasonal precipitation of ERA-I-Land with the observational datasets remains still higher than the one observed between ERA-I and observations. This sudden change can be attributed to a change in the ERA-I assimilation rain-affected radiances from SSM/I (Geer et al. 2008), but also the observed interdecadal shift that occurred over the region at the beginning of the 1990s (Huang et al. 2016, and references therein) could have played a role. Given that the spatial correlation between observations and ERA-I remains stable, the shifts in JJA precipitation might have been amplified by modified land modeling component in ERA-I-Land. This shift in precipitation translates into an amplified drought shift. Indeed, Balsamo et al. (2015) have shown that the changes in various simulated components of the hydrological cycle over Asia are mainly a result of model changes, whereas the bias correction in precipitation has been shown to be rather neutral over Eurasia.
Overall, AgMERRA outperforms the other three reanalyses. This performance of AgMERRA is somewhat expected due to the high level of statistical postprocessing (involving several different observational datasets) of the MERRA reanalysis.
This study has shown that the use of reanalyses for impact modeling and monitoring over monsoonal Asia should be based on an accurate evaluation process and choices should be tailored to the specific needs. Indeed, AgMERRA achieves the best performance, but it is not a near–real time (or almost near real time) updated dataset. ERA-I-Land shows good spatial performance as well, but when the interest is on the temporal evolution (such as drought analysis) JRA-55 should be preferred, as it exhibits the most stable temporal behavior.
Acknowledgments
We acknowledge three anonymous reviewers whose comments and suggestions greatly improved this study. We thank B. Casati for useful discussions on forecast verification and the FSS.
REFERENCES
Adler, R. F., and Coauthors, 2003: The version-2 Global Precipitation Climatology Project (GPCP) Monthly Precipitation Analysis (1979–present). J. Hydrometeor., 4, 1147–1167, doi:10.1175/1525-7541(2003)004<1147:TVGPCP>2.0.CO;2.
Balsamo, G., and Coauthors, 2015: ERA-Interim/Land: A global land surface reanalysis dataset. Hydrol. Earth Syst. Sci., 19, 389–407, doi:10.5194/hess-19-389-2015.
Bosilovich, M. G., J. Chen, F. R. Robertson, and R. F. Adler, 2008: Evaluation of global precipitation in reanalyses. J. Appl. Meteor. Climatol., 47, 2279–2299, doi:10.1175/2008JAMC1921.1.
Ceglar, A., A. Toreti, R. Lecerf, M. Van der Velde, and F. Dentener, 2016: Impact of meteorological drivers on regional inter-annual crop yield variability in France. Agric. For. Meteor., 216, 58–67, doi:10.1016/j.agrformet.2015.10.004.
Chen, C.-T., and T. Knutson, 2008: On the verification and comparison of extreme rainfall indices from climate models. J. Climate, 21, 1605–1621, doi:10.1175/2007JCLI1494.1.
Dee, D. P., and Coauthors, 2011: The ERA-Interim reanalysis: Configuration and performance of the data assimilation system. Quart. J. Roy. Meteor. Soc., 137, 553–597, doi:10.1002/qj.828.
Ding, Y., and J. C. Chan, 2005: The East Asian summer monsoon: An overview. Meteor. Atmos. Phys., 89, 117–142, doi:10.1007/s00703-005-0125-z.
Elliott, J., and Coauthors, 2014: Constraints and potentials of future irrigation water availability on agricultural production under climate change. Proc. Natl. Acad. Sci. USA, 111, 3239–3244, doi:10.1073/pnas.1222474110.
Feng, X., A. Porporato, and I. Rodriguez-Iturbe, 2013: Changes in rainfall seasonality in the tropics. Nat. Climate Change, 3, 811–815, doi:10.1038/nclimate1907.
Funk, C., and Coauthors, 2014: A quasi-global precipitation time series for drought monitoring. U.S. Geological Survey Data Series 832, 4 pp., doi:10.3133/ds832.
Funk, C., and Coauthors, 2015: The climate hazards infrared precipitation with stations—A new environmental record for monitoring extremes. Sci. Data, 2, 150066, doi:10.1038/sdata.2015.66.
Geer, A. J., P. Bauer, and P. Lopez, 2008: Lessons learnt from the operational 1D + 4D-Var assimilation of rain- and cloud-affected SSM/I observations at ECMWF. Quart. J. Roy. Meteor. Soc., 134, 1513–1525, doi:10.1002/qj.304.
Genovese, C., and L. Wasserman, 2004: A stochastic process approach to false discovery control. Ann. Stat., 32, 1035–1061, doi:10.1214/009053604000000283.
Harris, I., P. Jones, T. Osborn, and D. Lister, 2014: Updated high-resolution grids of monthly climatic observations—The CRU TS3.10 dataset. Int. J. Climatol., 34, 623–642, doi:10.1002/joc.3711.
Hatcho, N., S. Ochi, and Y. Matsuno, 2010: The evolution of irrigation developments in monsoon Asia and historical lessons. Irrig. Drain., 59, 4–16, doi:10.1002/ird.542.
Hirota, N., Y. N. Takayabu, M. Watanabe, and M. Kimoto, 2011: Precipitation reproducibility over tropical oceans and its relationship to the double ITCZ problem in CMIP3 and MIROC5 climate models. J. Climate, 24, 4859–4873, doi:10.1175/2011JCLI4156.1.
Hsu, K., and X. Gao, 1997: Precipitation estimation from remotely sensed information using artificial neural networks. J. Appl. Meteor., 36, 1176–1190, doi:10.1175/1520-0450(1997)036<1176:PEFRSI>2.0.CO;2.
Huang, D.-Q., J. Zhu, Y.-C. Zhang, Y. Huang, and X.-Y. Kuang, 2016: Assessment of summer monsoon precipitation derived from five reanalysis datasets over East Asia. Quart. J. Roy. Meteor. Soc., 142, 108–119, doi:10.1002/qj.2634.
Huffman, G. J., and Coauthors, 2007: The TRMM Multisatellite Precipitation Analysis (TMPA): Quasi-global, multiyear, combined-sensor precipitation estimates at fine scales. J. Hydrometeor., 8, 38–55, doi:10.1175/JHM560.1.
Joyce, R. J., J. E. Janowiak, P. A. Arkin, and P. Xie, 2004: CMORPH: A method that produces global precipitation estimates from passive microwave and infrared data at high spatial and temporal resolution. J. Hydrometeor., 5, 487–503, doi:10.1175/1525-7541(2004)005<0487:CAMTPG>2.0.CO;2.
Kobayashi, S., and Coauthors, 2015: The JRA-55 Reanalysis: General specifications and basic characteristics. J. Meteor. Soc. Japan, 93, 5–48, doi:10.2151/jmsj.2015-001.
McKee, T. B., N. J. Doesken, and J. Kleist, 1993: The relationship of drought frequency and duration time scales. Eighth Conf. on Applied Climatology, Anaheim, CA, Amer. Meteor. Soc., JP1.7. [Available online at http://ccc.atmos.colostate.edu/relationshipofdroughtfrequency.pdf.]
Noska, R., and V. Misra, 2016: Characterizing the onset and demise of the Indian summer monsoon. Geophys. Res. Lett., 43, 4547–4554, doi:10.1002/2016GL068409.
Pascale, S., V. Lucarini, X. Feng, A. Porporato, and S. ul Hasson, 2015: Analysis of rainfall seasonality from observations and climate models. Climate Dyn., 44, 3281–3301, doi:10.1007/s00382-014-2278-2.
Rana, S., J. McGregor, and J. Renwick, 2015: Precipitation seasonality over the Indian subcontinent: An evaluation of gauge, reanalyses, and satellite retrievals. J. Hydrometeor., 16, 631–651, doi:10.1175/JHM-D-14-0106.1.
Reichle, R. H., R. D. Koster, G. J. M. De Lannoy, B. A. Forman, Q. Liu, S. P. P. Mahanama, and A. Touré, 2011: Assessment and enhancement of MERRA land surface hydrology estimates. J. Climate, 24, 6322–6338, doi:10.1175/JCLI-D-10-05033.1.
Rienecker, M. M., and Coauthors, 2011: MERRA: NASA’s Modern-Era Retrospective Analysis for Research and Applications. J. Climate, 24, 3624–3648, doi:10.1175/JCLI-D-11-00015.1.
Roberts, N., and H. Lean, 2008: Scale-selective verification of rainfall accumulations from high-resolution forecasts of convective events. Mon. Wea. Rev., 136, 78–96, doi:10.1175/2007MWR2123.1.
Ruane, A. C., R. Goldberg, and J. Chryssanthacopoulos, 2015: Climate forcing datasets for agricultural modeling: Merged products for gap-filling and historical climate series estimation. Agric. For. Meteor., 200, 233–248, doi:10.1016/j.agrformet.2014.09.016.
Saha, S., and Coauthors, 2010: The NCEP Climate Forecast System Reanalysis. Bull. Amer. Meteor. Soc., 91, 1015–1057, doi:10.1175/2010BAMS3001.1.
Schneider, U., A. Becker, P. Finger, A. Meyer-Christoffer, M. Ziese, and B. Rudolf, 2014: GPCC’s new land surface precipitation climatology based on quality-controlled in situ data and its role in quantifying the global water cycle. Theor. Appl. Climatol., 115, 15–40, doi:10.1007/s00704-013-0860-x.
Shah, R., and V. Mishra, 2014: Evaluation of the reanalysis products for the monsoon season droughts in India. J. Hydrometeor., 15, 1575–1591, doi:10.1175/JHM-D-13-0103.1.
Singh, D., M. Tsiang, B. Rajaratnam, and N. S. Di, 2014: Observed changes in extreme wet and dry spells during the South Asian summer monsoon season. Nat. Climate Change, 4, 456–461, doi:10.1038/nclimate2208.
Taylor, K. E., 2001: Summarizing multiple aspects of model performance in a single diagram. J. Geophys. Res., 106, 7183–7192, doi:10.1029/2000JD900719.
Toreti, A., and P. Naveau, 2015: On the evaluation of climate model simulated precipitation extremes. Environ. Res. Lett., 10, 014012, doi:10.1088/1748-9326/10/1/014012.
Wang, B., 2006: The Asian Monsoon. Springer Science, 787 pp.
Wang, B., and Q. Ding, 2008: Global monsoon: Dominant mode of annual variation in the tropics. Dyn. Atmos. Oceans, 44, 165–183, doi:10.1016/j.dynatmoce.2007.05.002.
Willmott, C. J., and K. Matsuura, 1995: Smart interpolation of annually averaged air temperature in the United States. J. Appl. Meteor., 34, 2577–2586, doi:10.1175/1520-0450(1995)034<2577:SIOAAA>2.0.CO;2.
Yatagai, A., K. Kamiguchi, O. Arakawa, A. Hamada, N. Yasutomi, and A. Kitoh, 2012: APHRODITE: Constructing a long-term daily gridded precipitation dataset for Asia based on a dense network of rain gauges. Bull. Amer. Meteor. Soc., 93, 1401–1415, doi:10.1175/BAMS-D-11-00122.1.
Zhang, L., and Y. Qian, 2003: Annual distribution features of precipitation in China and their inter-annual variations. J. Meteor. Res., 17, 146–163.
Zou, X., P. Zhai, and Q. Zhang, 2005: Variations in droughts over China: 1951–2003. Geophys. Res. Lett., 32, L04707, doi:10.1029/2004GL021853.