Multivariate probabilistic estimates of heat stress for rice across China

Heat stress, a major threat to rice (Oryza sativa) production across China, would tend to increase in frequency and intensity under warming climate. Unlike probabilistic analysis via a univariate character, heat stress events, characterized by three variables (i.e., duration, peak and accumulated detrimental intensity), were identified in the past years. Nine distribution functions (i.e., Beta, Cauchy, Logistic, Normal, Exponential, Gamma, Lognormal, Weibull and Generalized Extreme Value) were firstly introduced and compared to select the best-fit marginal distribution of univariable by using Kolmogorov–Smirnov test, and seven copula functions (i.e., Normal and t, Gumbel–Hougaard, Clayton, Frank, Joe, Ali-Mikhail-Haq) were applied in the distributions of multivariables by Akaike Information Criterion statistics. It was obvious that higher magnitude was in the eastern parts in the context of heat stress frequency and characteristic variables. Critical values of heat stress variables corresponding to the certain return periods (i.e. 5, 10, 20 and 50 years) successively expanded in intensity and spatial scope. Inter-correlations of heat stress variables were significant, enlightening the importance of copula in connecting heat stress variables. The combined and co-occurrence bivariate and trivariate return period at certain univariate value corresponding to the given return periods, were consistent at the spatial scale. Accordingly, it was highlighted that eastern parts, especially Zhejiang, central-northern Fujian and eastern Jiangxi, were prone to heat stress, as a consequence of not only univariate but also multivariate probabilistic analysis. These results can be helpful in quantitatively assessing the vulnerability of rice to heat stress and provide us desired information of prevention strategies for heat stress.