It is essential to comprehend the relationship between agricultural yields and climatic conditions, especially concerning food security and the possible threats to crop output. Wheat is a crucial agricultural crop that covers a significant amount of rainfed production regions in Iran. This study utilized parametric and nonparametric approaches to assess rainfed wheat yield. The study centered on the Tabriz area in northwest Iran, examining precipitation patterns concerning rainfed wheat cultivation. The study focused on analyzing seasonal and distinct rainfall patterns during the cultivation period, utilizing widely recognized drought metrics such as the standardized precipitation index and the standardized precipitation–evaporation index. The study findings indicate a notable upward trend in rainfed wheat output over the analyzed period. The Mann–Kendall test resulted in a p-value of 0.031, indicating statistical significance for the observed rising trend. We conducted trend removal and normalized rainfed wheat yield figures based on seasonal precipitation to study the data more thoroughly. The second phase of the growing season was particularly notable, spanning from the completion of germination to the beginning of blooming. Instead of analyzing precipitation for the full growing season, concentrating on fall precipitation or the time from germination to blooming might improve yield forecasts and determinations.

  • This study utilized parametric and nonparametric approaches to assess rainfed wheat yield against drought indices.

  • The study findings indicate a notable upward trend in rainfed wheat output over the analyzed period.

The sustainable development of agricultural production is essential because of factors like population growth and the depletion of resources available for agriculture in highly populated countries. Reducing food losses as a result of extreme weather and climate change is one way to reach this goal. Nonetheless, yearly fluctuations in the climate pose a special risk to agricultural productivity (Gourdji et al. 2015; Wang et al. 2019; Eini et al. 2023b). However, as a result of the numerous physical, biological, and environmental systems that have been impacted by global warming brought on by human activity, agriculture has always been regarded as the most vulnerable industry (De Costa et al. 2007; Halkos & Tsilika 2017; Zhang et al. 2018; Zarei & Khaledi-Alamdari 2023).

One of the challenges facing the agricultural community is figuring out how to use weather information to implement risk management strategies, which will increase preparedness and decrease vulnerability to climate change (Meinke et al. 2003). Agriculture is a very risky business because of the significant weather fluctuations in arid and semi-arid regions, particularly in developing nations. To manage the impact of risks, climate risk management in these areas entails forecasting potential risks associated with climate change and offering solutions that can significantly reduce these risks (Ramesh 2010; Tomczyk et al. 2022).

To determine the optimal timing for agricultural operations and thereby maximize output, it is crucial to evaluate the risks posed by climate variability and extreme weather events to agricultural production (Andre et al. 2007). Therefore, researchers are increasingly interested in studying how climate-related factors impact the production of agricultural products. Consequently, there is a need to raise public awareness of the evaluation and prediction of these factors. Many factors that affect the climate are relevant to agricultural management and food security, particularly when it comes to rainfed agriculture in arid regions.

Among the most significant rainfed crops, wheat is a strategic plant vital to human nutrition and a staple food in many societies. Due to its widespread use and cultivation, this plant makes up the majority of all cultivations worldwide. As a result, its cultivated area is estimated to be around 17% of all cultivated area, and it provides 40% of the world's food supply (Giraldo et al. 2019; Tolwani & Shukla 2022). Climate variability has a major global impact on crop productivity. Remarkably, a number of researchers have reported that this variability causes a 60% change in yield (Ray et al. 2015; Matiu et al. 2017; Eini et al. 2023a).

Furthermore, rainfed agriculture is extremely susceptible to extreme weather events and changes (Zampieri et al. 2017; Ribeiro et al. 2019a), including drought, global warming, and heavy rainfall. As a result, extreme weather conditions like drought are thought to be a major threat to agricultural systems, with the effects being especially apparent in rainfed crops (Páscoa et al. 2017). Significant economic ramifications can result from quantitative and qualitative changes in variables linked to the success of strategic investments, which may also give rise to pervasive social issues. Rainfed crop yields are influenced by weather conditions and CO2, which are the most unpredictable (Szyga-Pluta et al. 2023; Marcinkowski & Piniewski 2024). Accurately estimating the weather in target areas aids farmers and planners in planting on schedule and providing for the plant's needs throughout the growing season. The development of quantitative and qualitative agricultural products may be particularly affected by this.

Furthermore, changes in agricultural crop yield in the region can be accurately predicted by determining and quantifying the impact of each climate variable. Plant performance, development, and growth conditions are significantly impacted by water stress. This is a severe issue in arid and semi-arid regions, where rainfall varies annually, and plants must withstand extended periods of water scarcity (Boyer 1982). As a result, it is evident that understanding how weather affects variations in crop yield is crucial (Faghih et al. 2021).

Analyzing climate trends using temperature and precipitation integration indices to identify and categorize droughts is essential for understanding how weather affects changes in crop yield (Vangelis et al. 2013). The most popular technique for researching droughts uses indicators that allow one to assess the drought's spatial and temporal characteristics as well as its relative intensity (Heim Jr 2002; Zargar et al. 2011). Agriculture's risk of drought, defined as the potential for crop failure due to drought conditions, frequently causes considerable financial losses (Skakun et al. 2016; Xie et al. 2018; Khaledi-Alamdari et al. 2023). Drought in the context of agriculture is characterized as a severe and ongoing lack of rainfall that drastically lowers agricultural output relative to normal (Warren et al. 1996), as well as meteorological drought, which is the outcome of the region's natural decline in rainfall over a specific period of time. It can occur over several years in a row or over a short period and is frequently represented by the difference between the annual and average precipitation totals. Compared to other natural disasters, drought is one of the most costly natural disasters in the world, impacting a large number of people and societies (Keyantash & Dracup 2002). Governments, regional and international organizations, and non-governmental organizations have focused especially on drought preparedness and policy during the past few decades (Wilhite 2000). Consequently, large-scale and occasionally local plant performance has been predicted using meteorological variables like temperature and precipitation (Bannayan et al. 2003; Palosuo et al. 2011). Farmers can take advantage of favorable climatic conditions to double production or make better decisions in the face of weather-related risks by accurately modeling multivariable weather distribution.

This study examined the variables that better justify the changes in rainfed wheat yield using meteorological and climatic parameters, including precipitation and temperature. The variables used were analyzed parametrically and nonparametrically to clarify the effects. The role of droughts during the studied period was assessed using the widely recognized standardized precipitation index (SPI) and the standardized precipitation–evaporation index (SPEI), which were compared to rainfed wheat yield. For the first time, this type of analysis provides very important and significant results for agricultural decision-making and management, as well as the discussion about food security.

Study area

As a part of the Lake Urmia Basin (Delavar et al. 2024), the Tabriz Plain is situated between 37°45′ N and 38°28′ N latitude and between 45°48′ E and 46°51′ E longitude. About 5,481 km2 of this study area is divided into 1,703 km2 of flat land and 3,778 km2 of mountainous terrain. Figure 1 shows the geographic location of the study area. According to the Emberger climate index, the Tabriz Plain has a semi-arid (cold and dry) climate with an average temperature of 13 °C and an average rainfall of about 250 mm. The Iran Meteorological Organization data center (https://data.irimo.ir) provided the daily meteorological data, and the Agricultural Organization of Tabriz Province, Iran, provided the wheat yield data. Figure 2 shows a diagram of the workflow.
Figure 1

Tabriz Plain location in Iran and Lake Urmia Basin.

Figure 1

Tabriz Plain location in Iran and Lake Urmia Basin.

Close modal
Figure 2

Flowchart of the current study.

Figure 2

Flowchart of the current study.

Close modal

The growth cycle of rainfed wheat and the calculation of growing degree days

The method described by Zadoks et al. (1974) identified 10 fundamental phases in the growth process of wheat: germination, seedling growth, tillering, stem elongation, booting, ear emergence, flowering, milk development, dough development, ripening, and ultimately, the plant reaches its full maturity stage. Each plant needs specific heat units to complete its distinct phenological stages. The length of the rainfed wheat growing season was measured in this study using the growing degree days (GDD) method, which is calculated using Equation (1) (Hundal et al. 1997).
(1)
where GDD is the growing degree days (cumulative heat), Tmin and Tmax are the minimum and maximum daily temperature in degrees Celsius, Tb is the temperature of the plant base in degrees Celsius, and finally, a and b are the starting time and end of the desired phenological stage, respectively, which is taken into account for calculating the entire wheat growing period according to the date of sowing and harvesting. It should be noted that the base temperature is the lowest temperature at which the plant is considered not to grow at lower temperatures. For wheat, this value is Tb = 0, so if the average daily temperature has a value equal to or less than Tb, the value of GDD = 0 is considered (Sharma et al. 2004).

To evaluate the role of each phase of wheat growth in relation to yield, the values of heat units during six phases are used, which consist of (1) germination, (2) completion of germination to onset of flowering, (3) flowering stage, (4) completion of flowering to the start of seed filling, (5) seed filling stage, and (6) total growing season.

Drought indices

Meteorological droughts occur before other types of droughts (Eini et al. 2023c), and their occurrence is considered a necessary condition for other droughts (Mirabbasi et al. 2012; Eini et al. 2023d). In the case of meteorological drought, the basis for calculating the degree of dryness is determined by comparing the amount of precipitation with the long-term average or its normal value. To assess and monitor climatic droughts, researchers worldwide use different methods, which are not accessible everywhere and are sometimes not efficient enough for other regions.

In the current study, the SPI was used to assess meteorological droughts. The SPI proposed by McKee et al. (1993) was calculated for each region based on the long-term precipitation occurring in each region (Edwards 1997). Basically, the SPI was developed to detect the lack of precipitation at multiple time scales. These time scales show the specific impacts of drought on the accessibility of different water sources. Instead of using just precipitation, as the SPI does, there seems to be a recent trend toward using the SPEI, which makes use of both temperature and precipitation (Dikshit et al. 2021). The SPEI is one of the indicators of drought, which can be determined at different scales by the difference between precipitation (P) and potential evapotranspiration or reference transpiration (ET0). This includes the influence of temperature on soil moisture, which affects plant growth.

Yield and trend analysis

To check the presence of trends in the variables under study such as temperature, precipitation, and yield, the Mann–Kendall test (Mann 1945; Kendall 1948) is considered the most widely used nonparametric testing method (Nikhil Raj & Azeez 2012; Zarch et al. 2015; Leng & Hall 2019; Szyga-Pluta et al. 2023).

Kendall's τ is a coefficient that indicates ordinal correlation and is used to determine the dependence criterion between random variables. This coefficient is generally used to determine copula functions.

Kendall's τ coefficient is defined as the difference between the probability of alignment or concordance and the probability of misalignment or discordance:
(2)
where X1X2 and (Y1Y2) are random vectors with uniform and independent distribution that have the same common cumulative distribution function (CDF) HXY(x, y). P[(X1X2)(Y1Y2) > 0] is the probability of concordance and P[(X1X2) s(Y1Y2) < 0] is the probability of discordance (Nelsen 2007). The pairs (xi, xj) and (yi, yj) are concordance if (xixj)(yiyj) > 0 and discordance if (xixj)(yiyj) < 0. Furthermore, when (xixj)(yiyj) = 0, these two pairs are neither concordance nor discordance. The Kendall coefficient τ is limited to the interval [−1, 1], where values 1 indicate complete concordance, −1 indicates complete discordance, and 0 indicates zero concordance. If (x1y1) and (x2y2) are two data pairs of sample space values of size n, Kendall's τ correlation can be calculated as follows:
(3)
where nc is the number of concordance pairs, nd is the number of discordance pairs, and n is the number of data. In case of duplicate data, the following equation is used:
(4)
where ti is the number of repeated observations in the ith group, nx is the number of repeated observations for the variable x, ui is the number of repeated observations in the ith group, and ny is the number of repeated observations for the variables x y. I there are no repetitions in the x and y data, then τ = τb, and the equation becomes the previous equation.

Precipitation status of the region

According to the trend recorded from the precipitation data of the Tabriz synoptic station, each type of precipitation impacts the region's climate and weather conditions. Therefore, all recorded rainfall events were included in the present study. The trend of changes in annual precipitation and long-term average monthly precipitation in the studied area are shown in Figures 3 and 4, respectively.
Figure 3

Long-term annual precipitation changes in the Tabriz Plain during the studied period.

Figure 3

Long-term annual precipitation changes in the Tabriz Plain during the studied period.

Close modal
Figure 4

Average monthly precipitation of the Tabriz Plain during the study period.

Figure 4

Average monthly precipitation of the Tabriz Plain during the study period.

Close modal

Rainfall at different growing seasons

The only source of water supply for rainfed wheat in the study area is the occurrence of precipitation at different times when most of the precipitation in this region falls during the growing season of rainfed wheat. In the area studied, according to available statistics, on average, more than 90% of precipitation occurs in each crop year and between October and May. Therefore, by extracting the precipitation that occurs only during the growing season, its influential role in the yield of rainfed wheat can be studied. Figure 5 shows annual growing season precipitation over the studied period from 1989 through the 2020 growing season.
Figure 5

Annual growing season precipitation (P) over the studied period.

Figure 5

Annual growing season precipitation (P) over the studied period.

Close modal
According to Figure 5, during the rainfed wheat growing period, significant changes can be seen in the state of annual precipitation, which has a very slight increasing trend; however, the trend test of Mann–Kendall (Mann 1945) confirms no significant trend at the 5% level (α = 0.05) with respect to growing season precipitation (p-value = 0.284), so it can be recognized that the observed trend for rainfall during the wheat growing season in the Tabriz Plain is not significant. The range of precipitation recorded in the studied area ranges from 102 mm in the lowest occurrence period (the growing season in the 2007–2008 crop year) to 348.5 mm in the highest occurrence (the growing season in the 1992–1993 crop year). The average rainfall during the growing season is estimated to be around 240 mm. Rainfed wheat goes through different stages during its growing period, including germination (S1), completion of germination until the beginning of flowering (S2), the flowering phase (S3), the end of flowering to the beginning of seed filling (S4), the seed filling phase (S5), and the entire growing season (GP). The precipitation occurring in each stage has a different contribution. In the studied area, the main contribution ranges from germination to seed filling, which is related to the second stage, i.e., from the completion of germination to the beginning of the flowering period. Figure 6 shows the percentage of precipitation in different phases of the wheat growing season as a long-term average.
Figure 6

The contribution of precipitation related to different stages of rainfed wheat growth in the Tabriz Plain.

Figure 6

The contribution of precipitation related to different stages of rainfed wheat growth in the Tabriz Plain.

Close modal

According to Figure 6, the proportion assigned to each stage, in descending order, refers to the following: the stages from the completion of germination to the beginning of flowering, from the completion of flowering to the beginning of seed filling, the flowering phase, the seed filling phase, and germination.

Rainfed wheat yield

Considering the precipitation status during the growing season and changes in rainfed wheat yield, these two variables have a good correlation (Figure 7). Therefore, to study the role of each seasonal rainfall in forming the final wheat yield, the relationship between performance and seasonal rainfall was established and presented after separating the rainfall of each season.
Figure 7

Changes in yield and seasonal precipitation during the studied period (yellow: fall; blue: winter; green: spring).

Figure 7

Changes in yield and seasonal precipitation during the studied period (yellow: fall; blue: winter; green: spring).

Close modal

Elimination of crop yield trends and anomalies

To study the yield of rainfed wheat with the drought index, in this research, the time series of the yield of rainfed wheat in the Tabriz Plain region during the last 30 years from the crop year 1991–1990 to 2019–2020 was studied, and its values are shown in Figure 7. In general, it can be seen that the improvements in agricultural methods, investments, and technological advancements have led to a continuous increase in crop yield during this period, but a sharp decline in yield is visible at times during the studied period. The Mann–Kendall trend test (Mann 1945) confirms a statistically significant positive trend at 5% probability (α = 0.05) in rainfed wheat yield (p-value = 0.031) to ensure that the observed trend does not influence the results. The data analysis was performed using the time-series data without trend and by detrending values and variance (Figure 8).
Figure 8

Yield of rainfed wheat in the Tabriz Plain region and detrended values.

Figure 8

Yield of rainfed wheat in the Tabriz Plain region and detrended values.

Close modal
Crop yield data, by its nature, has positive values, and standardizing these values can identify and reveal anomalies and changes; therefore, the yield data are extracted in a standardized manner, and the results are shown in Figure 9. As can be seen from Figure 9, the standardized yield of rainfed wheat in the Tabriz Plain reached the negative range in some periods of the studied period, and the negative values under these conditions are expressed as the risk of reduced crop yield (Leng & Hall 2019).
Figure 9

Standardized yield during the study period in the Tabriz Plain.

Figure 9

Standardized yield during the study period in the Tabriz Plain.

Close modal

Growing season SPI drought

Using the rainfed wheat growing season precipitation data, the SPI drought values were estimated for each growing season, and the results are shown in Figure 10. According to the results, the study area experienced mild, moderate, severe, and even very severe droughts at times. The most severe drought is associated with the 2007–2008 crop year. In the period studied (from the 1990–1991 crop year to the 2020–2019 crop year), drought periods occurred in approximately half of the growing seasons, and the values for the entire period ranged between −3.14 and 1.54.
Figure 10

SPI values in the rainfed wheat growth season in the Tabriz Plain.

Figure 10

SPI values in the rainfed wheat growth season in the Tabriz Plain.

Close modal

Standardized precipitation–evaporation index

The SPEI determines the drought status of the region by taking into account simultaneous precipitation and evapotranspiration. This index was calculated for the studied area over a period of 30 years (growing seasons 1990–1991 to 2019–2020), and its results are shown in Figure 11. The annual scale SPEI or the length of the growing period is calculated as a 12- or 9-month index, and its value is calculated for the wheat growing period, taking into account the results from the last month of that period (Wang et al. 2014; An et al. 2020). The advantage of this index is the use of the evapotranspiration parameter, and according to the results, the highest value of the SPEI drought is associated with the 2007–2008 growing year. During the studied period (from 1991–1990 to 2019–2020), 17 growing seasons experienced drought periods, and the values for the entire period ranged from −2.41 to 2.04.
Figure 11

SPEI values in the rainfed wheat growth season in the Tabriz Plain.

Figure 11

SPEI values in the rainfed wheat growth season in the Tabriz Plain.

Close modal

Study of precipitation and yields of rainfed wheat

The changes in annual precipitation, especially during the growing season of rainfed wheat, are evident when comparing the yield of rainfed wheat in each growing season during periods when precipitation decreased significantly. The yield of rainfed wheat also reacted to this and fell significantly. Figure 7 shows the yield compared to growing season precipitation during the study period. In the growing season 2007–2008, when the lowest precipitation was recorded during the growing season, the yield of rainfed wheat (detrended) in this region was also lower at the lowest level, 103 kg/ha. This is despite the fact that in the case of the highest precipitation of the growing season in the region, which was measured at 348.5 mm, the yield of rainfed wheat was not the highest, which could indicate a very strong response of the yield of rainfed wheat to small amounts of atmospheric precipitation.

Study of SPI and rainfed wheat yield

Figure 12 shows the time series of rainfed wheat yield without trend and compares it with the SPI of the Tabriz Plain region during the 30-year period from 1990 to 2020. As can be seen, crop production decreases dramatically during severe drought periods, so such sensitivity to moisture deficits due to poor precipitation can lead to significant economic losses after several wet years due to general expectations as well as management policies resulting from adequate productivity in rainy years (Mpelasoka et al. 2008). In the 1995–1996 and 2016–2017 growing seasons, it is clearly visible that after several consecutive years, the occurrence of a mild drought led to a significant decline in yield. Awareness of the risk posed by the impact of drought on rainfed wheat yield can largely prevent this damage. Madadgar et al. (2017) also pointed out that correlation coefficient analysis shows that precipitation (SPI) strongly correlates with annual crop yield.
Figure 12

Time series of rainfed wheat yield without trend compared to the SPI.

Figure 12

Time series of rainfed wheat yield without trend compared to the SPI.

Close modal

Study of the SPEI and rainfed wheat yield

The changes in rainfed wheat yield compared to the changes in the SPEI at the growing season scale of rainfed wheat for the study area during the 30-year period (1990–1991 to 2019–2020) are shown in Figure 13. The greatest decline in yield occurred in 2007–2008, which also saw the most severe drought situation according to the SPEI. In many other cases, the occurrence of drought caused yield decline, but, in some cases, the yield status of rainfed wheat was not proportional to these events despite the occurrence of drought or drought. Many factors can mitigate the impact of drought on farmers' fields. For example, agricultural management such as irrigation can significantly mitigate crop yield response to lack of rainfall (Tafoughalti et al. 2018; Eini et al. 2023e). Nevertheless, assessing drought risk to crop yield provides valuable information for adaptation and targeted risk reduction (Leng & Hall 2019).
Figure 13

Time series of rainfed wheat yield without trend compared to the SPEI.

Figure 13

Time series of rainfed wheat yield without trend compared to the SPEI.

Close modal

A recent study focused on the response of wheat and barley plants to different drought indicators in Spain. It showed the differing efficiencies of several drought indicators and highlighted the importance of the multiscale feature of droughts, especially the SPEI (Peña-Gallardo et al. 2019). The study by Ribeiro et al. (2019b) demonstrates the efficiency and appropriateness of using the SPEI in jointly analyzing the risk of crop yields and drought events.

Correlation and statistical analysis of parameters

Since many parameters play a role in determining the yield of rainfed wheat, it is essential to specify the parameters that play a significant role in the variations of rainfed wheat yield in the region. In Table 1, according to the parameters mentioned in the previous sections, their parametric and nonparametric correlations with the yield of rainfed wheat in three different states of normal yield, detrended yield, and standardized yield are given, to which the values refer to annual rainfall. Precipitation during six growing seasons of rainfed wheat, including germination (S1), completion of germination to start of flowering (S2), flowering stage (S3), completion of flowering to start of seed filling (S4), seed filling stage (S5), and the entire growing season (GP), is given in the table above. Based on the results obtained, the highest parametric and nonparametric (ordered) correlation of precipitation is related to the total growing season precipitation at normal yield and detrended yield. The Kendall correlation values for growing season precipitation relative to normal yield and detrended yield are estimated at 0.396 and 0.366, respectively, both significant at the 1% level.

Table 1

Parametric and nonparametric correlation of precipitations, SPI, and SPEI time series with rainfed wheat yield in the Tabriz region

Statistics
Yield
Standardized yield
Detrended yield
MeanMinimumMaximumKendallPearsonKendallPearsonKendallPearson
Precipitation Annual 256.91 146.3 362.23 0.387* 0.586* 0.356* 0.554* 0.375* 0.562* 
Growing season (GP) 240.43 102.08 348.51 0.396* 0.602* 0.347* 0.583* 0.366* 0.591* 
Fall 67.41 7.18 144.97 0.345* 0.541* 0.370* 0.515* 0.370* 0.514* 
Winter 72.59 18.41 168.31 0.327** 0.430** 0.209 0.344 0.228 0.376** 
Spring 100.9 27.92 181.95 0.041 0.137 0.076 0.198 0.076 0.186 
Summer 16.02 0.72 40.22 −0.12 −0.334 −0.145 −0.379** −0.163 −0.383** 
S1 4.24 24.01 0.245 0.301 0.168 0.261 0.196 0.271 
S2 181.22 87.28 297.74 0.396* 0.577* 0.320** 0.503* 0.338* 0.512* 
S3 18.5 0.41 47.67 0.041 0.066 0.025 0.096 0.025 0.098 
S4 30.27 1.03 64.77 0.138 0.126 0.259 0.117 0.254 
S5 6.21 35.05 0.229 0.237 0.24 0.261 0.249 0.265 
Drought index SPI −3.14 1.55 0.396* 0.665* 0.347* 0.652* 0.366* 0.658* 
SPEI −0.01 −2.41 2.04 0.253** 0.510* 0.361* 0.605* 0.361* 0.588* 
Statistics
Yield
Standardized yield
Detrended yield
MeanMinimumMaximumKendallPearsonKendallPearsonKendallPearson
Precipitation Annual 256.91 146.3 362.23 0.387* 0.586* 0.356* 0.554* 0.375* 0.562* 
Growing season (GP) 240.43 102.08 348.51 0.396* 0.602* 0.347* 0.583* 0.366* 0.591* 
Fall 67.41 7.18 144.97 0.345* 0.541* 0.370* 0.515* 0.370* 0.514* 
Winter 72.59 18.41 168.31 0.327** 0.430** 0.209 0.344 0.228 0.376** 
Spring 100.9 27.92 181.95 0.041 0.137 0.076 0.198 0.076 0.186 
Summer 16.02 0.72 40.22 −0.12 −0.334 −0.145 −0.379** −0.163 −0.383** 
S1 4.24 24.01 0.245 0.301 0.168 0.261 0.196 0.271 
S2 181.22 87.28 297.74 0.396* 0.577* 0.320** 0.503* 0.338* 0.512* 
S3 18.5 0.41 47.67 0.041 0.066 0.025 0.096 0.025 0.098 
S4 30.27 1.03 64.77 0.138 0.126 0.259 0.117 0.254 
S5 6.21 35.05 0.229 0.237 0.24 0.261 0.249 0.265 
Drought index SPI −3.14 1.55 0.396* 0.665* 0.347* 0.652* 0.366* 0.658* 
SPEI −0.01 −2.41 2.04 0.253** 0.510* 0.361* 0.605* 0.361* 0.588* 

*Significant in 0.01 level.

**Significant in 0.05 level.

Furthermore, the studied indices have an appropriate and acceptable correlation with respect to rainfed wheat yield, and Kendall's correlation for the SPI and SPEI with respect to normal yield are 0.396 and 0.253, respectively. These values suggest a detrending yield of 0.366 and 0.361, respectively. Based on the results presented in Table 1, four time series of total growing steps, seasonal precipitation, and SPI and SPEI were obtained.

This study provides a novel approach to rainfed wheat yield estimation by analyzing the specific impacts of precipitation during different growth stages, which is crucial for regions that rely heavily on rainfed agriculture. The key finding is the identification of the second growth phase – from the completion of germination to the beginning of flowering – as critically influential, with fall precipitation emerging as a particularly strong predictor of final yield. This early-stage precipitation monitoring offers a new perspective for agricultural planning and food security management.

The study's use of SPI and SPEI to correlate drought events with yield changes further contributes to the understanding of how seasonal precipitation patterns influence rainfed wheat production. This research is innovative in its detailed dissection of precipitation's role across different growth stages and its implications for more accurate yield forecasting in the face of climate variability.

The novelty of this manuscript lies in its focused examination of precipitation effects across distinct growth stages of rainfed wheat, particularly emphasizing the critical role of early-stage precipitation. This study introduces a novel methodology by separating precipitation occurrences based on the GDD and using both parametric and nonparametric correlation analyses to evaluate their impact. The implications of these findings are significant, offering actionable insights for early-stage agricultural interventions and more effective food security strategies in regions where rainfed agriculture is predominant.

Numerous factors influence rainfed wheat yield, and understanding these factors is crucial for improving crop yield estimation and forecasting. By analyzing precipitation during different growth stages using commonly used drought indicators, researchers can determine the impact of precipitation on rainfed wheat yield. This study highlights the significant role of precipitation during specific growth stages, particularly from the completion of germination to the beginning of flowering. However, it is important to acknowledge that all growth stages are critical; insufficient precipitation in the early stages cannot be compensated by more precipitation later on, underscoring the importance of the early stages.

The study also examined the seasonality of precipitation – specifically fall, winter, and spring precipitation – during the growing season of rainfed wheat in the Tabriz region. The results indicate that autumn precipitation plays a more significant role than other seasons, showing strong correlations with the final rainfed wheat yield. Early monitoring of fall rainfall can, therefore, be a valuable tool for predicting wheat yield, which is essential for planning and managing food security.

To account for drought events during the growing season, the study utilized the SPI and SPEI, both of which effectively explained variations in rainfed wheat yield in the region. Given that rainfed crops rely heavily on precipitation, drought events directly impact yield, and changes in these indices correspond to changes in crop yield.

Future research could benefit from comparing the findings of this study with similar studies in other regions or climates to assess the broader applicability of the results. Additionally, exploring new or improved methods for yield prediction, such as integrating remote sensing data or advanced machine learning techniques, could provide more accurate forecasts. Studies could also investigate the potential of other environmental factors, technological advancements, and socio-economic variables in influencing rainfed wheat yield, offering a more comprehensive understanding of the factors at play.

The authors thank the University of Tabriz for supporting this research. M.R.E. was supported by National Science Centre (Narodowe Centrum Nauki – NCN, PRELUDIUM BIS-1 project, UMO-2019/35/O/ST10/04392, Poland).

M.K.-A. and A.F.-F. conceptualized, analyzed, and prepared the figures. M.K.-A. and A.M.-H. prepared and provided a literature review and methodology investigation. A.F.-F. and M.R.E. supervised and validated the results. All authors reviewed and participated in the writing of the manuscript.

Data cannot be made publicly available; readers should contact the corresponding author for details.

The authors declare there is no conflict.

An
Q.
,
He
H.
,
Nie
Q.
,
Cui
Y.
,
Gao
J.
,
Wei
C.
,
Xie
X.
&
You
J.
(
2020
)
Spatial and temporal variations of drought in Inner Mongolia, China
,
Water
,
12
,
1715
.
André
R. G. B.
,
Garcia
A.
,
De Melo á Abreu
J. P.
,
Nieto
R.
,
Prasad
P. V. V.
,
White
D.
&
Garcia
A.
(
2010
)
Climate and weather risk assessment for agricultural planning. In: World Meteorological Organisation
.
Guide to agricultural meteorological practices
,
Geneva: World Meteorological Organisation, pp. 7.1–7.33
.
Bannayan
M.
,
Crout
N.
&
Hoogenboom
G.
(
2003
)
Application of the CERES-wheat model for within-season prediction of winter wheat yield in the United Kingdom
,
Agronomy Journal
,
95
,
114
125
.
Boyer
J. S.
(
1982
)
Plant productivity and environment
,
Science
,
218
,
443
448
.
De Costa
W.
,
Weerakoon
W.
,
Chinthaka
K.
,
Herath
H.
&
Abeywardena
R.
(
2007
)
Genotypic variation in the response of rice (Oryza sativa L.) to increased atmospheric carbon dioxide and its physiological basis
,
Journal of Agronomy and Crop Science
,
193
,
117
130
.
Delavar
M.
,
Raeisi
L.
,
Eini
M. R.
,
Morid
S.
,
Mohammadi
H.
&
Abbasi
H.
(
2024
)
Assessing the effectiveness of water-saving plans at the farm and basin level using agrohydrological modeling and water-accounting approaches
,
Journal of Irrigation and Drainage Engineering
,
150
,
04024009
.
Dikshit
A.
,
Pradhan
B.
&
Huete
A.
(
2021
)
An improved SPEI drought forecasting approach using the long short-term memory neural network
,
Journal of Environmental Management
,
283
,
111979
.
Edwards
D. C.
(
1997
)
Characteristics of 20th century drought in the United States at multiple time scales. Monterey, CA: Storming Media
.
Eini
M. R.
,
Motehayeri
S. M. S.
,
Rahmati
A.
&
Piniewski
M.
(
2023b
)
Evaluation of the accuracy of satellite-based rainfed wheat yield dataset over an area with complex geography
,
Journal of Arid Environments
,
212
,
104963
.
Eini
M. R.
,
Najminejad
F.
&
Piniewski
M.
(
2023c
)
Direct and indirect simulating and projecting hydrological drought using a supervised machine learning method
,
Science of the Total Environment
,
898
,
165523
.
Eini
M. R.
,
Rahmati Ziveh
A.
,
Salmani
H.
,
Mujahid
S.
,
Ghezelayagh
P.
&
Piniewski
M.
(
2023d
)
Detecting drought events over a region in Central Europe using a regional and two satellite-based precipitation datasets
,
Agricultural and Forest Meteorology
,
342
,
109733
.
Eini
M. R.
,
Salmani
H.
&
Piniewski
M.
(
2023e
)
Comparison of process-based and statistical approaches for simulation and projections of rainfed crop yields
,
Agricultural Water Management
,
277
,
108107
.
Faghih
H.
,
Behmanesh
J.
,
Rezaie
H.
&
Khalili
K.
(
2021
)
Climate and rainfed wheat yield
,
Theoretical and Applied Climatology
,
144
,
13
24
.
Giraldo
P.
,
Benavente
E.
,
Manzano-Agugliaro
F.
&
Gimenez
E.
(
2019
)
Worldwide research trends on wheat and barley: a bibliometric comparative analysis
,
Agronomy
,
9
,
352
.
Gourdji
S.
,
Läderach
P.
,
Valle
A. M.
,
Martinez
C. Z.
&
Lobell
D. B.
(
2015
)
Historical climate trends, deforestation, and maize and bean yields in Nicaragua
,
Agricultural and Forest Meteorology
,
200
,
270
281
.
Halkos
G. E.
&
Tsilika
K. D.
(
2017
)
Climate change effects and their interactions: an analysis aiming at policy implications
,
Economic Analysis and Policy
,
53
,
140
146
.
Heim
R. R.
Jr
(
2002
)
A review of twentieth-century drought indices used in the United States
,
Bulletin of the American Meteorological Society
,
83
,
1149
1166
.
Hundal
S.
,
Singh
R.
&
Dhaliwal
L. K.
(
1997
)
Agro-climate indices for predicting phenology of wheat (Triticum aestivum) in Punjab
,
Indian Journal of Agriculture Science
,
7
(
6
),
265
268
.
Kendall
M. G.
(
1948
)
Rank correlation methods. London: Charles Griffin & Company Limited
.
Keyantash
J.
&
Dracup
J. A.
(
2002
)
The quantification of drought: an evaluation of drought indices
,
Bulletin of the American Meteorological Society
,
83
,
1167
1180
.
Khaledi-Alamdari
M.
,
Majnooni-Heris
A.
,
Fakheri-Fard
A.
&
Russo
A.
(
2023
)
Probabilistic climate risk assessment in rainfed wheat yield: copula approach using water requirement satisfaction index
,
Agricultural Water Management
,
289
,
108542
.
Madadgar
S.
,
AghaKouchak
A.
,
Farahmand
A.
&
Davis
S. J.
(
2017
)
Probabilistic estimates of drought impacts on agricultural production
,
Geophysical Research Letters
,
44
,
7799
7807
.
Mann
H. B.
(
1945
)
Nonparametric tests against trend
,
Econometrica
,
13
(
3
),
245
259
.
https://doi.org/10.2307/1907187
.
McKee
T. B.
,
Doesken
N. J.
&
Kleist
J.
(
1993
)
The relationship of drought frequency and duration to time scales
,
Proceedings of the 8th Conference on Applied Climatology
,
Boston
, pp.
179
183
.
Meinke
H.
,
Abawi
Y.
,
Stone
R. C.
,
Hammer
G. L.
,
Potgieter
A. B.
,
Nelson
R.
,
Howden
S. M.
,
Baethgen
W.
&
Selvaraju
R.
(
2003
)
Climate risk management and agriculture in Australia and beyond: linking research to practical outcomes
,
The International Conference on Total Disaster Risk Management
, p.
4
.
Mirabbasi
R.
,
Fakheri-Fard
A.
&
Dinpashoh
Y.
(
2012
)
Bivariate drought frequency analysis using the copula method
,
Theoretical and Applied Climatology
,
108
,
191
206
.
Mpelasoka
F.
,
Hennessy
K.
,
Jones
R.
&
Bates
B.
(
2008
)
Comparison of suitable drought indices for climate change impacts assessment over Australia towards resource management
,
International Journal of Climatology
,
28
,
1283
1292
.
Nikhil Raj
P.
&
Azeez
P.
(
2012
)
Trend analysis of rainfall in Bharathapuzha River basin, Kerala, India
,
International Journal of Climatology
,
32
,
533
539
.
Nelsen
R. B.
(
2007
)
An introduction to copulas.
Berlin: Springer Science & Business Media
.
Palosuo
T.
,
Kersebaum
K. C.
,
Angulo
C.
,
Hlavinka
P.
,
Moriondo
M.
,
Olesen
J. E.
,
Patil
R. H.
,
Ruget
F.
,
Rumbaur
C.
&
Takáč
J.
(
2011
)
Simulation of winter wheat yield and its variability in different climates of Europe: a comparison of eight crop growth models
,
European Journal of Agronomy
,
35
,
103
114
.
Páscoa
P.
,
Gouveia
C.
,
Russo
A.
&
Trigo
R.
(
2017
)
The role of drought on wheat yield interannual variability in the Iberian Peninsula from 1929 to 2012
,
International Journal of Biometeorology
,
61
,
439
451
.
Peña-Gallardo
M.
,
Vicente-Serrano
S. M.
,
Domínguez-Castro
F.
&
Beguería
S.
(
2019
)
The impact of drought on the productivity of two rainfed crops in Spain
,
Natural Hazards and Earth System Sciences
,
19
,
1215
1234
.
Ramesh
K.
(
2010
)
More efficient use of agricultural inputs as part of adoption of preparedness strategies: multiple cropping
. In: Stigter, K. (ed.)
Applied Agrometeorology
,
Springer
, pp.
435
440
.
Ray
D. K.
,
Gerber
J. S.
,
MacDonald
G. K.
&
West
P. C.
(
2015
)
Climate variation explains a third of global crop yield variability
,
Nature Communications
,
6
,
1
9
.
Ribeiro
A. F.
,
Russo
A.
,
Gouveia
C. M.
&
Páscoa
P.
(
2019a
)
Modelling drought-related yield losses in Iberia using remote sensing and multiscalar indices
,
Theoretical and Applied Climatology
,
136
,
203
220
.
Ribeiro
A. F.
,
Russo
A.
,
Gouveia
C. M.
,
Páscoa
P.
&
Pires
C. A.
(
2019b
)
Probabilistic modelling of the dependence between rainfed crops and drought hazard
,
Natural Hazards and Earth System Sciences
,
19
,
2795
2809
.
Sharma
A.
,
Sood
R.
&
Kalubarme
M.
(
2004
)
Agrometeorological wheat yield forecast in Himachal Pradesh
,
Journal of Agrometeorology
,
6
,
153
160
.
Skakun
S.
,
Kussul
N.
,
Shelestov
A.
&
Kussul
O.
(
2016
)
The use of satellite data for agriculture drought risk quantification in Ukraine. Geomatics
,
Natural Hazards and Risk
,
7
,
901
917
.
Szyga-Pluta
K.
,
Tomczyk
A. M.
,
Piniewski
M.
&
Eini
M. R.
(
2023
)
Past and future changes in the start, end, and duration of the growing season in Poland
,
Acta Geophysica
,
71
,
3041
3055
.
Tafoughalti
K.
,
El Faleh
E.
,
Moujahid
Y.
&
Ouargaga
F.
(
2018
)
Wheat yield vulnerability: relation to rainfall and suggestions for adaptation
,
Journal of Degraded and Mining Lands Management
,
5
,
1251
.
Tolwani
S.
&
Shukla
R.
(
2022
)
Role of genetic variability, heritability and genetic advance for various traits of bread wheat (Triticum aestivum L.) under variable sowing times
,
International Journal of Bio-Resource & Stress Management
,
13
,
122
130
.
Tomczyk
A. M.
,
Piniewski
M.
,
Eini
M. R.
&
Bednorz
E.
(
2022
)
Projections of changes in maximum air temperature and hot days in Poland
,
International Journal of Climatology
,
42
(
10
),
5242
5254
.
Vangelis
H.
,
Tigkas
D.
&
Tsakiris
G.
(
2013
)
The effect of PET method on reconnaissance drought index (RDI) calculation
,
Journal of Arid Environments
,
88
,
130
140
.
Wang
Q.
,
Wu
J.
,
Lei
T.
,
He
B.
,
Wu
Z.
,
Liu
M.
,
Mo
X.
,
Geng
G.
,
Li
X.
&
Zhou
H.
(
2014
)
Temporal-spatial characteristics of severe drought events and their impact on agriculture on a global scale
,
Quaternary International
,
349
,
10
21
.
Wang
B.
,
Feng
P.
,
Chen
C.
,
Li Liu
D.
,
Waters
C.
&
Yu
Q.
(
2019
)
Designing wheat ideotypes to cope with future changing climate in South-Eastern Australia
,
Agricultural Systems
,
170
,
9
18
.
Warren
A.
,
Sud
Y.
&
Rozanov
B.
(
1996
)
The future of deserts
,
Journal of Arid Environments
,
32
,
75
89
.
Wilhite
D. A.
(
2000
)
Drought as a natural hazard: concepts and definitions. In: Wilhite, D. A. (ed.) Drought: a global assessment, vol. I, London: Routledge, pp. 3–18
.
Xie
W.
,
Xiong
W.
,
Pan
J.
,
Ali
T.
,
Cui
Q.
,
Guan
D.
,
Meng
J.
,
Mueller
N. D.
,
Lin
E.
&
Davis
S. J.
(
2018
)
Decreases in global beer supply due to extreme drought and heat
,
Nature Plants
,
4
,
964
973
.
Zampieri
M.
,
Ceglar
A.
,
Dentener
F.
&
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
,
064008
.
Zadoks
J. C.
,
Chang
T. T.
&
Konzak
C. F.
(
1974
)
A decimal code for the growth stages of cereals
.
Weed Research
,
14
,
415
421
.
https://doi.org/10.1111/j.1365-3180.1974.tb01084.x
Zarei
G.
&
Khaledi-Alamdari
M.
(
2023
)
Zoning the heating and cooling requirements of several major greenhouse plants in Iran
,
Journal of Agricultural Meteorology
,
11
,
27
35
.
Zargar
A.
,
Sadiq
R.
,
Naser
B.
&
Khan
F. I.
(
2011
)
A review of drought indices
,
Environmental Reviews
,
19
,
333
349
.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY 4.0), which permits copying, adaptation and redistribution, provided the original work is properly cited (http://creativecommons.org/licenses/by/4.0/).