Probabilistic Analysis of Drought Effects on Rainfed Wheat Yield Using Copula Functions (Case Study: Tabriz Plain)

Authors

1 PhD. Student of Irrigation and Drainage, Department of Water Science and Engineering, Agriculture faculty, University of Tabriz. Tabriz. Iran.

2 Water Eng. Dept, Tabriz University

3 Professor, Department of Water science and engineering, Agriculture faculty, University of Tabriz.

Abstract

Extended Abstract
Background and Objectives
A combined assessment of drought risk and associated impacts on crop production based on a probabilistic approach seems appropriate to understand the multivariate nature of drought risk in agriculture. To overcome the problems caused by drought impact detection, several approaches have been developed in recent decades. Among the multivariate analysis approaches, copula functions are very popular. Copulas use univariate marginal distributions to form a joint distribution. The joint distribution can be described by the corresponding marginal distributions and copula functions that describe the dependency structure. In this research, using the statistical precipitation data in the Tabriz plain in eastern part of Lake Urmia basin, and yield of rainfed wheat in this area, a model based on copula functions was developed to determine the diffrent probabilities of yield risk and different drought conditions. Also, the application of copula functions related to rainfed wheat yield in this region was performed for the first time, and the presented method will be applicable to other areas and other crop cultivation.
Methodology
In the case of meteorological drought, the basis for calculating the degree of drought is determined by comparing the amount of precipitation with the long-term average or its normal values. The SPI index is considered to be an appropriate and powerful index to use as a time scale droughts monitoring. Basically, SPI was created to detect the lack of precipitation on multiple time scales. Among the reasons that make the use of this index so popular, we can mention the standard nature of this index as it can be used in regional studies and establish a temporal relationship between drought events in different parts of the same area. The SPI index is a dimensionless index and its more negative values show the more severe the drought.
Analysis of variables individually is easy and can be analyzed by statistical distribution functions; but statistical analysis joint variables is very complicated and impossible in most cases. If the correlation criterion of these variables is known, their joint probability distribution can be obtained using copula functions. Using copula functions for modeling has a high degree of flexibility as it is possible to choose different marginal distributions to create a multivariate model. Copulas are functions that form a bivariate or multivariate distribution based on two or more univariate marginal functions. Several copula functions can be used to construct a two-dimensional joint distribution of hydrological and agricultural variables, among which Archimedean and elliptical copula families are the most commonly used. In the present study, six copula functions are used and the parameters of the paired functions are determined using the two-stage maximum likelihood method, which estimates the parameters of the marginal distribution and the copula function by forming two likelihood functions. In order to investigate joint probability of rainfed wheat yield and drought index, the time series of rainfed wheat yield in the Tabriz Plain region and SPI index during the last 30 years from 1990 to 2020 was applicated in this study.
Findings
In general, the improvement in agricultural methods, investments and technological advances during this period have led to a continuous increase in crop yield, however, a sharp drop in crop yield is clearly evident at times during the reporting period. In this research, to ensure that the observed trend does not affect the results, the Copula model was built using the detrended time series data by removing the trend in the values and the variance in the yield data. Based on the results obtained, crop production decreases dramatically during severe droughts, so such sensitivity to moisture deficits caused by low rainfall after several wet years can be attributed to farmers' expectations and management policies driven by high productivity during the previous ones wet years were determined.
To use series of standardized yield and growing season SPI in copula, the most appropriate distributions were selected and used, logistic and normal, respectively. Also, according to the calculated Kendall correlation (0.35), the best fit joint was Claytons function with AIC = -11.10, RMSE = 0.036 and used to construct the joint probability distribution of the standardized yield series of rainfed wheat and SPI of growth Period. Results showed that the cumulative probability of yield risk event in mild, moderate, severe, and very severe drought was 18.42, 7.82, 4.26, and 2.32 percent, respectively,and for an overall rainfed wheat yield risk and drought condition is about 33%. Meanwhile, the probability of rainfed wheat yield risk and non-drought events is only about 7% and joint probability of yield risk and SPI>1 is so close to zero.
Conclusion
In the current research, the probability of occurrence of the rainfed wheat yield risk and drought conditions was extracted by copula functions. Based on the results obtained, most changes in yield occur when the drought index is in the range greater than -1.5. In other words, in the conditions of severe and very severe drought, the yield of rainfed wheat does not show noticeable changes in probability. The opposite situation can also be observed for the situation of very wet years, that when there is high amounts of precipitation and the drought index show values higher than 1.5, the yield of the crop does not show much change. Thus, the SPI must be greater than one in order to achieve the desired yield and not be at risk of rainfed wheat yield, which in this study is considered to be a standardized yield values greater than zero. Because in this range of SPI, the joint probability of yield risk is estimated to be very close to zero. Therefore, this threshold can be introduced as the rainfed wheat yield safety threshold, but depending on the situation of the region and the drought index, the drought threshold (SPI<0) contains cumulative probability about 33% of the rainfed wheat yield risk.

Keywords


Angelidis P, Maris F, Kotsovinos N and Hrissanthou V, 2012. Computation of drought index SPI with alternative distribution functions. Water Resource Management 26(9): 2453–2473.
Arbenz P, 2013. Bayesian copulae distributions, with application to operational risk management—some comments. Methodology and Computing in Applied Probability 15:105–108.
Bahremand A, Alvandi E, Bahrami M, Dashti Marvili M, Heravi H.,Khosravi GR, Kornejady A,Samadi Arghini H,Tajiki M and Teimouri M, 2016. Copula functions and their application in stochastic hydrology. Journal of Conservation and Utilization of Natural Resources 4(2):1-20. ( in Persian with English abstract)
Ben-Ari T, Adrian J, Klein T, Calanca P, Van der Velde M and Makowski D, 2016. Identifying indicators for extreme wheat and maize yield losses. Agricultural and Forest Meteorology 220: 130-140.
Bokusheva R, Kogan F, Vitkovskaya I, Conradt S and Batyrbayeva M, 2016. Satellite-based vegetation health indices as a criteria for insuring against drought-related yield losses. Agricultural and Forest Meteorology 220:200-6.
Challinor AJ, Ewert F, Arnold S, Simelton E and Fraser E, 2009. Crops and climate change: progress, trends, and challenges in simulating impacts and informing adaptation. Journal of Experimental Botany 60: 2775-2789.
Chiew FH, Piechota TC, Dracup JA and McMahon TA, 1998. El Nino/Southern Oscillation and Australian rainfall, streamflow and drought: Links and potential for forecasting. Journal of Hydrology 204(1-4):138-49.
Edwards DC and McKee TB, 1997. Characteristics of 20th century drought in the United States at multiple time scales. Atmospheric Science Paper 634:1-30.
Elliott J, Müller C, Deryng D, Chryssanthacopoulos J, Boote KJ, Büchner M, Foster I, Glotter M, Heinke J, Iizumi T and Izaurralde RC, 2015. The global gridded crop model intercomparison: Data and Modeling Protocols for Phase 1 (V1. 0). Geoscientific Model Development 8(2):261-77.
Fang HB, Fang KT and Kotz S, 2002. The meta-elliptical distributions with given marginals. Journal of Multivariate Analysis 82:1-16.
Folberth C, Elliott J, Müller C, Balkovic J, Chryssanthacopoulos J, Izaurralde RC, Jones CD, Khabarov N, Liu W, Reddy A and Schmid E, 2016. Uncertainties in global crop model frameworks: effects of cultivar distribution, crop management and soil handling on crop yield estimates. Biogeosciences Discussions [preprint] 1-30.
Frank MJ, 1979. On the simultaneous associativity of F (x, y) and x + y − F (x, y). Aequationes Mathematicae 19:194–226.
Genest C, Favre AC, Béliveau J and Jacques C, 2007. Metaelliptical copulas and their use in frequency analysis of multivariate hydrological data. Water Resources Research 43(9):1-12.
Genest C and Rivest LP, 1993. Statistical inference procedures for bivariate Archimedean copulas. Journal of the American Statistical Association 88(423):1034-43.
Godfray HCJ, Beddington JR, Crute IR, Haddad L, Lawrence D, Muir JF, Pretty J, Robinson S, Thomas SM and Toulmin C, 2010. Food security: the challenge of feeding 9 billion people. Science 327(5967):812-818.
Gumbel EJ, 1960. Distributions of extreme values in several dimensions. Paris Institute of Statistics 9:171-173.
Hernandez-Barrera S, Rodriguez-Puebla C and Challinor AJ, 2017. Effects of diurnal temperature range and drought on wheat yield in Spain. Theoretical and Applied Climatology 129(1):503-19.
Hlavinka P, Trnka M, Semeradova D, Dubrovský M, Žalud Z and Možný M, 2009. Effect of drought on yield variability of key crops in Czech Republic. Agricultural and Forest Meteorology 149 (3–4):431–442.
Huang J, Zhuo W, Li Y, Huang R, Sedano F, Su W, Dong J, Tian L, Huang Y and Zhu D, 2020. Comparison of three remotely sensed drought indices for assessing the impact of drought on winter wheat yield. International Journal of Digital Earth 13: 504-526.
Huang X and Wang Z, 2018. Probabilistic spatial prediction of categorical data using elliptical copulas. Stochastic Environmental Research and Risk Assessment 32:1631-1644.
Lee T, Modarres R and Ouarda TB, 2013. Data‐based analysis of bivariate copula tail dependence for drought duration and severity. Hydrological Processes 27:1454-1463.
Leng G and Hall J, 2019. Crop yield sensitivity of global major agricultural countries to droughts and the projected changes in the future. Science of the Total Environment, 654:811-21.
Lesk C, Rowhani P and Ramankutty N, 2016. Influence of extreme weather disasters on global crop production. Nature 529 (7584): 84–87.
Li Y, Gu W, Cui W, Chang Z and Xu Y, 2015. Exploration of copula function use in crop meteorological drought risk analysis: a case study of winter wheat in Beijing, China. Natural Hazards 77:1289-1303.
Li C, Singh VP and Mishra AK, 2013. A bivariate mixed distribution with a heavy-tailed component and its application to single-site daily rainfall simulation. Water Resources Research, 49:767–789.
Lobell DB, Roberts MJ, Schlenker W, Braun N, Little BB, Rejesus RM and Hammer GL, 2014. Greater sensitivity to drought accompanies maize yield increase in the US Midwest. Science 344(6183):516-9.
Madadgar S, AghaKouchak A, Farahmand A and Davis SJ, 2017. Probabilistic estimates of drought impacts on agricultural production. Geophysical Research Letters 44 (15): 7799–7807.
Mann HB, 1945. Non-parametric tests against trend. Econometrica 13:163–171.
Matiu M, Ankerst DP and Menzel A, 2017. Interactions between temperature and drought in global and regional crop yield variability during 1961–2014. PLoS One 12 (5): e0178339.
McKee TB, Doeskin NJ and Kleist J, 1993. The relationship of drought frequency and duration to time scales. Pp.179-184. Proceedings of the 8th Conference on Applied Climatology. 17-22 January, Anaheim, California.
Mirabbasi R, Fakheri-Fard A and Dinpashoh Y, 2012. Bivariate drought frequency analysis using the copula method. Theoretical and Applied Climatology 108:191-206.
Mpelasoka F, Hennessy K, Jones R and Bates B, 2008. Comparison of suitable drought indices for climate change impacts assessment over Australia towards resource management. International Journal of Climatology: A Journal of the Royal Meteorological Society 28(10):1283-92.
Nelsen RB, 2006. An Introduction to Copulas. Springer Science & Business Media. Springer New York, NY.
Páscoa P, Gouveia CM, Russo A and Trigo RM, 2017. The role of drought on wheat yield interannual variability in the Iberian Peninsula from 1929 to 2012. International Journal of Biometeorology 61(3):439-51.
Potopova V, Boroneant C, Boincean B and Soukup J, 2016. Impact of agricultural drought on main crop yields in the Republic of Moldova. International Journal of Climatology 36(4): 2063–2082.
Quiring SM, 2009. Developing objective operational definitions for monitoring drought. Journal of Applied Meteorology and Climatology 48: 1217-1229.
Ribeiro AF, Russo A, Gouveia CM and Páscoa P, 2019a. Modelling drought-related yield losses in Iberia using remote sensing and multiscalar indices. Theoretical and Applied Climatology 136(1):203-20.
Ribeiro AF, Russo A, Gouveia CM, Páscoa P and Pires CA, 2019b. Probabilistic modelling of the dependence between rainfed crops and drought hazard. Natural Hazards and Earth System Sciences 19(12): 2795-2809.
Rosenzweig C, Elliott J, Deryng D, Ruane AC, Müller C, Arneth A, Boote KJ, Folberth C, Glotter M, Khabarov N and Neumann K, 2014. Assessing agricultural risks of climate change in the 21st century in a global gridded crop model intercomparison. Proceedings of the National Academy of Sciences 111(9):3268-73.
Salvadori G and De Michele C, 2004. Frequency analysis via copulas: Theoretical aspects and applications to hydrological events. Water Resources Research 40(12): 1-17.
Sanainejad SH, Ansari H, Davari K and Morid S, 2003. Monitoring and assessment of drought severity in mashhad at different time scales using standardized precipitationindex (SPI). Iranian Journal of Soil Research 17(2): 201-209. ( in Persian with English abstract) 
Shi W and Tao F, 2014. Vulnerability of African maize yield to climate change and variability during 1961–2010. Food Security 6 (4): 471–481.
Sklar M, 1959. N-dimensional distribution functions and their margins. Publications of the Statistical Institute of the University of Paris 8:229-231.
Srinivas S, Menon D and Meher Prasad A, 2006. Multivariate simulation and multimodal dependence modeling of vehicle axle weights with copulas. Journal of Transportation Engineering 132(12):945-55.
Tao F, Zhang Z, Liu J and Yokozawa M, 2009. Modelling the impacts of weather and climate variability on crop productivity over a large area: a new super-ensemblebased probabilistic projection. Agricultural and Forest Meteorology 149 (8):1266–1278.
Tebaldi C and Lobell D, 2008. Towards probabilistic projections of climate change impacts on global crop yields. Geophysical Research Letters 35 (8):1-6.
Tilman D, Balzer C, Hill J and Befort BL, 2011. Global food demand and the sustainable intensification of agriculture. Biological Sciences 108 (50): 20260–20264.
Troy T, Kipgen C and Pal I, 2015. The impact of climate extremes and irrigation on US crop yields. Environmental Research Letters 10(5): 1-10.
Van Dijk AI, Beck HE, Crosbie RS, De Jeu RA, Liu YY, Podger GM, Timbal B and Viney NR, 2013. The millennium drought in southeast Australia (2001–2009): Natural and human causes and implications for water resources, ecosystems, economy, and society. Water Resources Research 49(2):1040-1057.
Yu C, Li C, Xin Q, Chen H, Zhang J, Zhang F, Li X, Clinton N, Huang X, Yue Y and Gong P, 2014. Dynamic assessment of the impact of drought on agricultural yield and scale-dependent return periods over large geographic regions. Environmental Modelling & Software 62:454-464.
Zampieri M, Ceglar A, Dentener F and Toreti A, 2017. Wheat yield loss attributable to heat waves, drought and water excess at the global, national and subnational scales. Environmental Research Letters 12 (6): 1-11
Zipper SC, Qiu J and Kucharik CJ, 2016. Drought effects on US maize and soybean production: spatiotemporal patterns and historical changes. Environmental Research Letters 11 (9): 1-11