Climate Change and Food Security: The Impact of Some Key Variables on Wheat Yield in Kazakhstan

Abstract

In such drought-prone regions as Kazakhstan, research on regional drought characteristics and their formation conditions is of paramount importance for actions to mitigate drought risks caused by climate change. This paper presents the results of research on the spatio-temporal patterns of atmospheric droughts as one of the most important factors hindering the formation of crop yields. The influence of several potential predictors characterizing teleconnection in the coupled “atmosphere–ocean” system and cosmic-geophysical factors affecting their formation is analyzed. The spatial relationships between atmospheric aridity at the individual stations of the investigated area and the wheat yield in Kazakhstan as well as its relationships with potential predictors were determined using econometric methods. High correlation was shown between wheat yield fluctuations and Multivariate El-Niño–Southern Oscillation (ENSO), galactic cosmic radiation, solar activity, and atmospheric drought expressed through the soil moisture index, which in turn depends on precipitation levels and temperatures. The model could be modified further so that the individual components could be forecasted into the future using various time series in an ARIMA model. The resulting integration of these forecasts would allow the prediction of wheat yields in the future. The obtained results can be used in the process of creating effective mechanisms for adaptation to climate change and droughts based on their early diagnosis.

1. Introduction

There is a strong sense of urgency in dealing with climate change. The latest IPCC report outlines that “Human activities are estimated to have caused approximately 1.0 °C of global warming above preindustrial levels, with a likely range of 0.8 °C to 1.2 °C. Global warming is likely to reach 1.5 °C between 2030 and 2052 if it continues to increase at the current rate” [1]. The United Nations announced 17 Sustainable Development Goals, which include SDG1 ‘No poverty’, SDG2 ‘Zero hunger’, SDG3 ‘Good health and wellbeing’, SDG11 ‘Sustainable cities and communities’, and SDG13 ‘Climate action’, most closely related to climate change and its impacts on food security [2]. Climate change is starting to have a significant impact on ecosystems and food production. The latest IPCC report asserts that “Climate models project robust differences in regional climate characteristics between present-day and global warming of 1.5 °C, and between 1.5 °C and 2 °C. These differences include increases in: mean temperature in most land and ocean regions, hot extremes in most inhabited regions, heavy precipitation in several regions, and the probability of drought and precipitation deficits in some regions” [1]. Anomalies in rainfall, changing seasons, droughts, and heat waves are appearing with increasing frequency, which is also affecting agricultural yields and prices. All the major crops, wheat, rice, coffee, cocoa, maize, sugarcane, tea, and grapes, will likely experience a change in yields as a result of changes in climate. IPCC assesses potential reductions in wheat yields due to climate change in the range of 2% per decade [3]. It is therefore important to identify and address various climatic, economic, and institutional factors and be able to forecast the likely changes to avoid food shortages or price spikes and protect the livelihoods of vulnerable communities. According to IPCC, “Populations at disproportionately higher risk of adverse consequences of global warming of 1.5 °C and beyond include disadvantaged and vulnerable populations, some indigenous peoples, and local communities dependent on agricultural or coastal livelihoods” [1].

Kazakhstan uses 80% of its territory for growing agricultural crops, thereby creating employment for over 7 million people. The production of wheat is one of the key strategically important sectors of Kazakhstan’s economy (70% of total grain production), which largely determines the country’s food security. Agriculture is most vulnerable to climate change, frequent extreme weather events, and water stress. Each additional °C in global temperatures is likely to reduce wheat yields by 6% [1]. Drought is one of the major natural disasters in Kazakhstan and exerts a considerable pressure on the agricultural sector, which is the most vulnerable sector of the economy to drought. As for the irrigated agriculture, glaciers in the mountains are the major storage of water resources and they are extremely sensitive to climate [4,5,6]. Both natural drought and climate change-induced variation in glaciers pose a great threat to the agriculture of Kazakhstan.

At the same time, the anthropogenic impact on the climatic system enters into complex interactions with the natural factors, influenced by the rotational cycling, determining the difficulty in predicting climatic changes. Every factor therefore must be viewed in a framework of diverse interactions with positive and negative feedback loops. The most famous impacts are those determined by atmospheric teleconnection, which leads to significant changes in seasonal precipitation, which, in turn, influences droughts in some regions and floods in others, breaking the regular weather patterns on Earth. The El Niño is defined by FAO as a “naturally occurring phenomenon characterized by the abnormal warming of sea surface temperature in the central and eastern equatorial Pacific Ocean. During El Niño episodes, normal patterns of tropical precipitation and atmospheric circulation are disrupted, triggering extreme climate events around the globe”. The effect of El Niño occurring every two to seven years and lasting up to 18 months could be predicted up to 5 months in advance. If such a teleconnection exists, the extreme consequences of adverse events could be adapted to, building resilience in the agricultural sector.

Accepting the 2030 Sustainable Development Agenda and the Paris Agreement on Climate Change, the international community took responsibility for designing a sustainable future. It should be noted that in order to not only eliminate hunger and poverty by 2030, but also tackle climate change, major structural changes in food production and agricultural systems will be needed all around the world [7].

One of the most important consequences of climatic changes is the increasing frequency of extreme natural events (droughts, floods, hail, heat waves, extreme frosts, etc.), which make a considerable impact on the sustainability of agricultural systems. In order to build a risk management system in relation to these events based on predictions and advance warnings, it is necessary to study in depth the mechanisms of their development, identify the most informative predictors, and build the management system accordingly. Droughts are one of the most significant and dangerous natural events for agriculture in many parts of the world. Due to the complex mechanisms of their occurrence, they are difficult to parameterize and predict [8,9,10].

This paper develops a comprehensive model of climate change impacts on drought and crop yields affecting food security in Kazakhstan, based on consideration of the insufficiently studied teleconnections in the ocean–atmosphere system and several other cosmic and geophysical factors. There is a significant gap in the literature on the climate–water–food nexus in the context of Kazakhstan experiencing growing impacts of climate change, and our paper will provide a new modelling approach allowing the assessment and prediction of impacts of changing climate on wheat yields. The originality of our paper stems from the inclusion of the potential explanatory variables: average monthly temperatures and precipitation levels for the studied region, changes in Total Solar Irradiance, changes in the Wolf index, cosmic rays, soil moisture, North Atlantic Oscillation (NAO), and the Multivariate El Niño–Southern Oscillation Index (MEI). All variables have been included with optimal time lags. The model seeks to explore optimal parameterization of atmospheric drought, and identification of the connection between the standardized moisture index and crop yield fluctuation, as well as analyzing the large-scale atmospheric processes related to drought affecting crop yields in Kazakhstan. The paper is structured as follows: Section 2 presents a detailed literature review of the state-of-the art in climate change impacts on crop yields, Section 3 presents a case study of Kazakhstan in greater depth, Section 4 presents our econometric modeling results, including the simultaneous equations model and the ARIMA forecasting model, Section 5 offers a discussion, and Section 6 concludes.

2. Literature Review

There is considerable literature devoted to the study of factors influencing the development of droughts, early warnings, and their impacts on the yields of various agricultural crops. It has been noticed that climate change has indeed influenced agricultural productivity [11]. Crop prices and vegetable oil prices have also been found to respond to climatic variations [12]. The studies found great geographical variation among the magnitude and the nature of climatic teleconnections.

The first group of studies adopted a global perspective and revealed negative impacts of El Niño for approximately 22–24% of harvested areas in the world, including maize in the USA, China, Indonesia, and Mexico, wheat in the USA, China, Mexico, Australia, and Europe, rice in China, Tanzania, and Myanmar, and soybean in China and India. At the same time, it was found that El Niño makes a favorable impact on 30–36% of harvested areas, including wheat in Argentina, Kazakhstan, and South Africa, rice in China, Brazil, and Indonesia, soybean in Brazil and the USA, and maize in Argentina and Brazil [13]. The global studies identified Sahel, western and southern Africa, Australia and the western USA, South America and central Eurasia as experiencing flood anomalies, especially in La Niña years in western USA and El Niño years in central Eurasia [14]. Researchers have identified Galactic Cosmic Rays (GCR) as having a marked impact on the Earth’s cloud cover and climate variability [15,16].

A significant research cluster focused on the climate impacts for agriculture in Asia, the focus of the present study. Climate change impacts on crop yields in Asia were estimated to fluctuate +/−20% with south and southeast Asian countries more severely affected than east Asian countries [17]. The researchers showed that El Niño–Southern Oscillation (ENSO) had an influence on drought variation over the whole region of central Asia; NAO had significant impact on all the subregions. NAO and the Standardized Evapotranspiration Index (SPEI) were found to demonstrate significant correlations with yields in north Kazakhstan in 1974, 1983, 1992, 1996, 2005, and 2011 [8]. The studies focusing on wheat concluded that wheat productivity in China in 1978–2007 was affected by both climate change and adaptive management. Temperatures were not found to have a significant impact on wheat yield; however, the importance of the negative impact of precipitation as well as positive impacts of irrigation in China have been confirmed [18]. Research focused on wheat and maize identified a strong connection between the ENSO phases and crop yields in China and pointed towards high wheat yield probability distributions for both El Niño and La Niña years [19]. Researchers have found the Standardized Precipitation Index (SPI) and the Reconnaissance Drought Index (RDI) to be more appropriate than the SPEI for the arid regions in China [20]. In the context of India, rice productivity was shown to be correlated with monsoon fluctuations and reduced 4.3–13.8% during moderate and strong El Niño years [21]. Crop yield anomalies in India were revealed to be inversely correlated with Sea Surface Temperature (SST) anomalies from June to August over the Tropical Pacific Ocean (TPO) NINO 3 sector, rice being affected more than wheat with significant economic effects of up to USD 2.2 billion reductions in revenue in warm ENSO periods and up to USD 1.3 billion increases in revenue in the cold ENSO periods [22]. In the context of Pakistan, crop yields have been shown to be statistically correlated with solar radiation and temperatures [23]. At the same time, it was revealed that an increase of 1 °C in mean temperature during the germination and tillering stages is capable of reducing yields by up to 7.4% while a similar increase during the vegetative growth stage improved productivity by up to 6.4% [24]. In a study focused on China, Japan, and Australia, up to 73% of the yield variation was explained by biomass production, which was strongly positively correlated with intercepted radiation [25].

Studies focusing on North America found that monthly values of Pacific North American Pattern (PNA) and ENSO indices could predict summer weather in Canada up to 2–7 months in advance [26]. Crop yields have been shown to be higher in La Niña phases and lower than the trend immediately following La Niña years in south-eastern USA. Corn and tobacco yields as well as bell pepper harvests seem to be significantly affected by ENSO events, which doesn’t apply to tomatoes [27]. ENSO dynamics was found to significantly influence palm oil, soybean oil, sunflower seed oil, and rapeseed oil prices in the USA, where ENSO has been shown to include an auto-regressive nonlinearity [12]. In the context of Brazil, ENSO brings higher probability of increased rainfall in April and May in the Sao Paolo region, which could affect planting dates and yields [28]. Land use and land cover changes (LULCC) have also been found to influence albedo, humidity, precipitation, and temperature in South America [29]. In Africa, crop yields have been shown to correlate with ENSO events with El Niño having a negative impact on millet and a positive impact on maize and sorghum production with an indirect effect on temperature and rainfall during the short planting period in the Sahel [30]. Warm ENSO events in east Africa were found to coincide with drier conditions affecting surface temperatures and the Normalized Difference Vegetation Index (NDVI), while Pacific SST anomalies tend to affect air temperature in east Africa, according to Plisnier et al. [31], ENSO has been shown to be responsible for variations in temperature, wind, and humidity in Zambia, with three rainy seasons corresponding to La Niña, normal, and El Niño conditions [32].

In Europe, orange and tangerine yields have been shown to be strongly correlated with ENSO fluctuations, while wheat, orange, and lemon yields exhibited correlations with NAO index. In Spain, wheat, sunflower, barley, rye, olive, grapes, and tangerines were revealed to have lower yields during La Niña than El Niño phases. Moreover, there is an increasing interest in applying short-term climate forecasts to predict the risk of droughts for vulnerable crops [33]. The mean values of the forecast are based on atmospheric conditions, surface conditions, surface, and the depths of the ocean, which could all influence weather conditions. The Development of a European Multi-model Ensemble System for Seasonal-to-Inter-annual Prediction (DEMETER) is a project which is funded by the European Union and considers seven state-of-the-art general circulation models [34]. In detail, the models show seasonal simulations that cover the same period of 21 years. The state-of-the-art seasonal forecasting systems are presented for both maximum temperature and precipitation for the four seasons in Spain using a multimodel ensemble from the DEMETER project [34].

Atmospheric circulation is one of the determining factors of drought appearance [35]. A number of studies showed that annual variations of meteorological characteristics of the atmosphere in the northern hemisphere and hydrological variables in the ocean are associated with the ENSO and NAO phenomena [36,37,38,39,40,41]. During the El Niño phase, warmer and drier climate conditions are a crucial explanation for the negative impacts of El Niño on the crops, such as wheat in Australia, soybean in India, rice in Indonesia, and maize in Zimbabwe and Brazil [13,21,27,28]. To better assess drought in a warming climate [42] investigated the drought characteristics across different aridity zones with and without consideration of potential evapotranspiration (PET).

Numerous attempts have been made to find the ENSO response in the Atlantic–Eurasian sector [8,39,43]. El Niño contributes to a change in the localization and intensity of north Atlantic and Mediterranean storm tracks by changing the global midlatitude equator–pole temperature gradient during the ENSO period (Mokhov et al., 1995) [44]. A significant influence of ENSO on the NAO was detected in [43,45].

The influence of ENSO on air temperature over the Eurasian continent was explored by Kryjov and Park (2007) [46]. The influence works through the so-called annular mode of the northern hemisphere, an alternation of periods of an intense circumpolar vortex in the lower stratosphere and corresponding western transport in the upper troposphere to the north of 45° N. and the periods of weakening of the circumpolar vortex and the western transport. Positive air temperature anomalies are noted over all temperate and circumpolar latitudes of Eurasia in the first case, with negative temperatures observed otherwise. Significant connections appear in the second half of the 20th and beginning of the 21st centuries, whereas the analysis of longer time series (100 years and more) did not show such interactions between processes in the equatorial Pacific Ocean and the North Atlantic [25].

In the spring–summer period of El Niño (La Niña) year, the atmospheric circulation in the north Atlantic is weakening (intensifying). This contributes to the formation of positive (negative) anomalies in the SST. In spring, the biggest differences in SST in the El Niño and La Niña years are observed off the coast of Africa (5–15° N, 15–30° W), and in summer—east of Newfoundland (45–55° N, 35–45° W).

Nesterov [47] shows that the El Niño–Southern Oscillation phenomenon has the greatest impact on atmospheric circulation in the Atlantic–European region in winter, coinciding with the phase of maximum development of the phenomenon, still having an effect next winter. In the first winter of the El Niño (La Niña) period, a positive (negative) phase of the east Atlantic is energized. Next winter—a positive (negative) phase of the north Atlantic oscillation takes place. The negative phase of both oscillations is more clearly expressed on composite maps of geopotential and surface pressure. The La Niña year is characterized by more variability of the indices of the main fluctuations in the atmospheric circulation in the Atlantic–European region.

The North Atlantic Oscillation (NAO) index peaked in the early 1990s, and has been dropping steadily to the present. The NAO index exhibits a 60–70-year cycle and demonstrates a significant positive correlation with TS in the northern hemisphere [41]. Periods with positive NAO indices are characterized by intense western transport of air masses and notable warming of the most part of extra-tropical latitudes of the northern hemisphere, especially strong in the winter–spring season. NAO has been considered to be one of the possible factors of climatic change; numerous studies are devoted to the effect of NAO on weather and climate in Europe, and the links between the NAO and the state of the stratospheric polar vortex. Currently, there is widespread understanding of how El Niño–Southern Oscillation events play a crucial role in determining agricultural production. Furthermore, problems arising from variability of water availability and soil degradation are currently major challenges to agriculture globally. These problems are likely to be exacerbated in the future if global climate change reaches 1.5 °C and 2.0 °C [1].

The key publications showing the results of relevant empirical research on climate change impacts on agriculture in different global regions and specific countries are presented and analyzed in

Table 1. Comparative analysis of literature on the climatic effects on crop production.

№	Author, Year	Region and Country	Research Focus	Method	Variables	Conclusions
1.	Iizumi et al., 2014 [13]	Global: USA, China, Kazakhstan, Mexico, Tanzania, Australia, India, etc.	Maize, soybean, rice, wheat	Statistical analysis: polynomial regression	Sea surface temperature (SST) anomaly, Oceanic Niño Index, La Niña, surface air temperature, and soil moisture	El Niño results in negative impacts on yields in 22–24% of harvested areas worldwide: maize in USA, China, Mexico, and Indonesia; soybean in India and China; rice in China, Myanmar, and Tanzania; and wheat in China, USA, Australia, Mexico, and parts of Europe. In contrast, El Niño impacts positively on 30–36% of harvested areas worldwide: wheat in Argentina, Kazakhstan, and South Africa; rice in China, Indonesia, and Brazil; soybean in the USA and Brazil; and maize in Brazil and Argentina.
2.	Ward et al., 2014 [14]	Global except Antarctica and Greenland	Flood risk	Cascade of hydrological and hydraulic models; statistical analyses: Spearman’s rank	Daily gridded discharge and flood volume at a horizontal resolution of 0.5° ×0.5° forced by daily meteorological fields; El Niño and La Niña	Flood anomalies revealed during El Niño and La Niña years: 34% of Earth surface excluding Antarctica and Greenland in El Niño years and 38% in La Niña years, respectively. At the regional scale strong flood anomalies are observed in the Sahel, western and southern Africa, Australia, the western United States (esp. in La Niña years), South America, and central Eurasia (especially during El Niño years).
3.	Svensmark, 2000 [15]	Global	Cosmic rays and Earth’s climate	Analytical review	Cosmic ray flux, solar cycle, Earth’s climate; clouds; Earth’s temperature; volcanic dust; El Niño–Southern Oscillation	During the last solar cycle, changes in the Earth’s cloud cover are correlated with the Galactic Cosmic Ray variations. It was found that on the time horizon over 10 years the sun has a significant influence on climate variations based on correlation between isotopes and Earth’s temperature for the last 1000 years assessed by proxy data.
4.	Ormes, 2018 [48]	Global	Cosmic rays and climate	Analytical review	Solar irradiance and climate; The Little Ice Age; cosmic ray ion intensity; direct effects of cosmic rays on climate: the Svensmark Effect; cloud classification	Galactic Cosmic Ray ionization effect in the troposphere could amount to 5–10%, whereas at latitudes above 50, this could be as high as 15–20%, affecting cloud formation through stimulating cloud condensation nuclei.
5.	Sidorenkov, 2016 [49]	Global	Synchronization of atmospheric processes with frequencies of the Earth–moon–sun system	Analytical review; comparative analysis	The movement of the poles of the Earth; quasi-two-year cyclic (QTC) zonal wind; the Chandler movement of the Earth’s poles (CMP); Multanovsky’s Natural synoptic periods (NSP); decade climate change	The author demonstrates how the 35-year Brückner cycle is correlated with the beat of the annual solar (365 days) and annual lunar (355 days) changes in meteorological parameters. A hypothesis is proposed on how the lunisolar tides influence the air temperature via the radiation conditions in the atmosphere due to fluctuations in the cloud cover. The author concludes that climatic characteristics, temperature, monsoons, masses of ice sheets, etc., correlate with changes in the speed of rotation of the Earth.
6.	Guo et al., 2018 [8]	Central Asia	Spatial and temporal characteristics of droughts in central Asia during 1966–2015	Principle Components Analysis (PCA); varimax rotation method; the Sen’s slope and the Modified Mann–Kendall method (MMK); continuous wavelet analysis	Climatic indices: North Atlantic Oscillation, Siberian High Index, Tibetan Plateau Index; PCA regionalization; Standardized Precipitation Evapotranspiration Index; drought	The drought patterns differ substantially in central Asia with most significant droughts observed in 1973–1979 (NW, NK, SW, and SE); 1983–1988 (SE and HX), and 1997–2003 (NE, NK, SW, SE, and HX) and an overall wetting tendency between 1966 to 2015. Most regions except north Kazakhstan exhibit a drying trend in 2003–2015. ENSO has an influence on drought variation over the whole region of central Asia; NAO has significant impact everywhere except HX. NAO and SPEI demonstrate short-period strongly significant correlations in northern Kazakhstan for 1974, 1983, 1992, 1996, 2005, and 2011.
7.	Min et al., 2014 [18]	China	Wheat yield	Exploratory spatial data analysis (ESDA): spatial autocorrelation and spatial heterogeneous dynamic analysis of wheat yield; Durbin model	Wheat yield; the average monthly precipitation; the average monthly temperature during; input of fertilizer; input of seed; effective irrigation area of wheat; disaster area of wheat	Wheat productivity in China in 1978–2007 is both affected by climate change and adaptive management. Temperature does not have a significant impact on wheat yield and precipitation has a significant negative effect. Wheat yield is significantly influenced by irrigation and damaged area during its growth. The authors infer that the adaptive management has played a more important role than climate change for wheat yield.
8.	Liu et al., 2014	[19] China	Wheat and maize	Statistical analysis: linear regression, Student’s t-test, Kolmogorov–Smirnov (KS) test. Agricultural production systems simulator (APSIM)	Daily weather data, data on yields of winter wheat and summer maize, SST over the tropical Pacific (ENSO index)	Crop yields on the North China Plain were affected by ENSO phases. The wheat yield probability distributions were exceptionally high for El Niño and La Niña years, although in Nanyang region the differences between El Niño and La Niña events were significant only at the 5% level. For maize the lowest probability distributions were seen in El Niño years and reached the confidence level of 95% in La Niña years.
9.	Xu et al., 2015 [20]	China	Climate variability and trends at a national scale	Analytical review; regression model; China Meteorological Data Sharing System; fractal dimensions and trend indices of each climatic factor for 1960–2013 in the 579 meteorological stations; Mann–Kendall tests for trends; Hurst index of meteorological parameters; mapping by ArcGIS10 software	3 months Standardized Precipitation Index (SPI3), 3 months Reconnaissance Drought Index (RDI3), 3 month Standardized Precipitation Evapotranspiration Index (SPEI3)	High spatial resolution daily climate variability data from 1960 to 2013 were used in conjunction with a 3 months Standardized Precipitation Index (SPI3), 3 months Reconnaissance Drought Index (RDI3), and a 3 months Standardized Precipitation Evapotranspiration Index (SPEI3). The most severe droughts of 1962–1963 and 2010–2011 in China affected over half of the nonarid regions. These droughts were the strongest in regions from the North China Plain to the downstream of the Yangtze River. SPI and RDI were found to be more appropriate than SPEI in the arid regions.
10.	Subash and Gangwar, 2014 [21]	India	Kharif rice yield	Interannual variability, standard deviation, coefficient of variation, Mann–Kendall nonparametric test	Monthly rainfall data, production and productivity of kharif rice, Oceanic Niño Index (ONI)	Rice productivity was found to be related to monsoon fluctuations. Six out of eight El Niño years exhibited a reduction in rainfall during the monsoon season, ranging from −20.3% in 2002 to −5.5% in 1991. A higher deficit was observed in July and September. During 5 out of 8 moderate and strong El Niño years the kharif rice productivity reduced −4.3% in 1986 and −13.8% in 2002. There exists a wide spatial variability in rice productivity. July rainfall influenced 71% of the variations in rice productivity.
11.	Selvaraju, 2003 [22]	India	Rice, sorghum, maize, pigeon pea, blackgram wheat, chickpea	Hamming-type filter, three-point weighted moving average using Hanning weights, correlation coefficients, Kruskal–Wallis (KW) H-test	Time series of actual total food grain production, NINO3 SST index, averaged SST anomalies, summer monsoon rainfall (SMR)	Crop yield anomalies in India were found to be inversely and statistically significantly related to SST anomalies from June to August over TPO NINO3 sector. Rice tends to be more affected than wheat. Total yields decreased by 1–15% during warm ENSO phase. Economic impacts were found to be significant, reducing revenue by up to USD 2.2 billion during a warm ENSO period and increasing revenues by USD 1.3 billion in a cold ENSO period. NINO3 ENSO index was concluded to be a good policy predictor for agriculture in India.
12.	Iglesias et al., 1996 [17]	India, Bangladesh, Pakistan, Japan, south China	Soybeans, maize, wheat, roots, tubers, fruit, vegetables	General circulation models (GCMs); climate scenarios; IBSNAT crop models	Temperature and precipitation, doubled CO2	The estimated effects of climate change on crop yields in Asia measure up to +/- 20%. South and southeast Asian countries were found to be more affected; east Asia was considered to be less vulnerable. Fluctuations in ENSO event frequency and severity were shown to be having an impact on agricultural production in Asia.
13.	Ahmed and Hassan, 2011 [23]	Pakistan	Wheat yield	Statistical analysis; scatter plot regression model	Crop data and climate variables (T1, T2, SR1, SR2, PTQ1, and PTQ2)	Yields were found to be statistically related to solar radiation and temperature.
14.	Ahmad et al., 2014 [24]	Pakistan	Wheat yield	Regression model	Solar radiation, crop yields	Climatic change is having an impact on wheat yields in Pakistan. An increase of 1 °C in the mean temperature during the germination and tillering stages is capable of reducing yields by 7.4%. Conversely, a 1 °C increase in the mean temperature during the vegetative growth stage improves productivity by 6.4%.
15.	Yulihastin et al., 2008 [50]	The Indonesian Maritime Continent (IMC)	Rainfall	Dipole mode index	El Niño–Southern Oscillation (ENSO) and the Indian Ocean dipole (IOD) mode, SST	ENSO and IOD tend to influence rainfall in Indonesia. During a combined strong El Niño phase and a strong positive IOD, Indonesia experiences negative anomalies of precipitations during June–August, September–November, and December–February.
16.	Diamond et al., 2012 [51]	Southwest Pacific	Southwest Pacific tropical cyclone	20th-century reanalysis data and the coupled ENSO index (CEI)	Mean sea level pressure difference between Papeete, French Polynesia and Darwin, Australia; El Niño and La Niña atmospheric circulation	Positive correlations were found between Tropical Cyclones (TCs), sea surface temperature, and atmospheric circulation. The period of 1991–2010 exhibited higher TC frequency as opposed to the 1970–1990 period. TCs were found to be connected to sea surface temperature (SSTs) in the SW Pacific basin and atmospheric circulation in specific regions.
17.	Huang et al., 2016 [25]	Japan, Australia, and China	Rice yield	Analytical review; linear regression; box diagrams	Rice; grain yield; biomass production; harvest index; radiation use efficiency (RUE); solar radiation; intercepted radiation	Up to 73% and 6% of the yield variation was explained by biomass production and the harvest index, respectively. Biomass production was strongly positively correlated to intercepted radiation.
18.	Garnett et al., 1998 [26]	Canada	Summer weather forecasting	Multiple regression analysis; linear correlation coefficients	PNA and ENSO indices, summer mean temperature, and mean total precipitation	Monthly values of PNA and ENSO indices could predict summer weather in Canadian prairies 2–7 months in advance, which could be beneficial for agriculture.
19.	Zhao et al., 2019 [52]	USA	Cereals, roots, vegetables, and fruits	SIMPLE crop model: model calibration, sensitivity analysis; simulation	13 parameters; climate impact	The SIMPLE crop model includes 13 parameters including crop type, daily weather data, crop management, and soil water holding parameters, and was calibrated for 14 crops. The limitations of the model include lack of response to vernalization and photoperiod effect on phenology, and the model doesn’t include nutrient dynamics.
20.	Ubilava and Holt, 2013 [12]	USA	Corn production	Smooth Transition Vector Error Correction Models (STVECM); Smooth Transition Autoregressive Models (STAR)	Panel of crop yield, temperature, precipitation, and ENSO data (anomaly, Niño 3.4)	ENSO dynamics were found to include an autoregressive nonlinearity. Palm oil, soybean oil, sunflower seed oil, and rapeseed oil are significantly affected by climate anomalies driven by ENSO events.
21.	Hansen et al., 1999 [27]	Southeastern USA	Peanut, tomato, cotton, tobacco, corn, and soybean	Canonical correlation analysis; Fourier analysis; spectral analysis, analysis of variance (ANOVA)	Historical records of yields, area harvested, price, and total value of production, sea surface temperature anomalies	Crop yields were shown to be higher in La Niña phases and lower than the trend in years immediately following La Niña phases. Corn and tobacco yields were found to be statistically significantly affected by ENSO phases. Bell pepper prices exhibited a correlation with ENSO events, which didn’t apply to tomatoes.
22.	Soler et al., 2010 [28]	Brazil	Maize yield	CSM-CERES-Maize model, DSSAT model	El Niño, La Niño daily weather data; soil characterization data; cultivar coefficients; crop management information	A higher probability of increased rainfall in April and May in Sao Paulo during an El Niño phase was found. The study explored the significance of ENSO phases on planting dates and yields.
23.	Beltrán- Przekurat et al., 2012 [29]	Argentina	Land-use and land-cover change	Coupled atmospheric- biospheric regional climate GEMRAMS model	The near-surface energy balance, temperature, humidity, and precipitation	Land-use/land-cover changes (LULCC) were found to influence the albedo, near-surface energy balance, humidity, precipitation, and temperature over southern South America.
24.	Okonkwo and Demoz, 2014 [30]	Sahel (Africa)	Maize, millet, and sorghum	Wavelet analysis (CWT, XWT), correlation analysis	Niño 3.4, SST, precipitation, temperature, soil moisture	Statistically significant correlation between ENSO and crop yield in the Sahel region was confirmed. El Niño events appear to have a negative impact on millet yield and a positive impact on maize and sorghum production. The indirect impact of ENSO on crop yield is through its effect on temperature and available rainfall during the short planting period in Sahel.
25.	Plisnier et al., 2000 [31]	East Africa	East African ecosystems, Lake Tanganyika	Statistical analysis: moving average (MA) filter	SOI, rainfall, air humidity, NDVI, Ts	A positive and significant correlation between the air temperature in east Africa and Pacific SST anomalies was found. Warm ENSO events tend to coincide with drier conditions. Correlation coefficients between NDVI and Ts anomalies and the ENSO index were found to be highly significant.
26.	Makaudze, 2014 [53]	Southern Africa, Zimbabwe	Maize yield	The Decision Support System for Agricultural Technology program, DSSAT v4	Crop inputs, daily weather data: solar radiation, rainfall, temperature evapotranspiration, El Niño, La Niño	Climate forecasts have shown to lead to smallholder farmers in Zimbabwe recording higher yield gains (28%) compared to those without forecasts.
27.	Kanno et al., 2015 [32]	Zambia	Rainy seasons and rainfall simulations	Weather Research and Forecasting Model (WRF) ver. 3.2	Air temperature, relative humidity, precipitation, wind direction and speed	The three rainy seasons in Zambia corresponded in turn with La Niña, normal, and El Niño conditions. The intra-seasonal variations in precipitation appeared to be driven by ENSO. Variations in temperature, wind, and humidity over time appeared to be highly sensitive to large-scale atmospheric movements, such as ENSO.
28.	Gimeno et al., 2002 [33]	Spain	Wheat, rye, barley, oats, sunflower seeds, olives, grapes, citrus fruits	ANOVA analysis; cross-correlation functions (CCF)	ENSO and NAO phases; crop yields	Orange and tangerine yields were found to be highly correlated to fluctuations in ENSO. Wheat, orange, and lemon showed correlation with the NAO index. Wheat, sunflower, barley, rye, olive, grape, and tangerine tend to exhibit lower yields during La Niña phases than during El Niño phases. At the same time, oats, rye, wheat, and citrus yields are higher in positive NAO phases in Spain.
29.	Frias et al., 2010 [34]	Spain	Precipitation and temperature seasonal forecasts	Multimodel DEMETER system; statistical downscaling method	Seven global atmosphere–ocean coupled models, each running from an ensemble of nine initial conditions, ENSO events	ENSO is found to be affecting dry events in spring in the south and Mediterranean coast and hot events in the southern part of Spain during El Niño years. La Niña events manifest themselves in dry events in winter in the northwest of Spain and hot events in the summer in the south and Mediterranean coast.

, explaining the research focus, key methods, variables, and research conclusions.

3. Case Study

Kazakhstan, a vast country of 2.7 million km², is almost equidistantly positioned between the Atlantic, Pacific, and Indian Oceans at the center of the Eurasian continent, which determines its continental climate. The territory of Kazakhstan is represented by four climatic zones: deserts (44%), semideserts (14%), steppe (26%), and forest-steppe (6%), and is vulnerable to climatic change.

The geography of central Asia is complex and consists of high mountains including the Tianshan and Karakoram ranges, low endorheic basins including the Tarim Basin, vast deserts including the Gobi, Kara Kum, Kyzyl Kum, and Taklimakan, and treeless grasslands in the northern part of Kazakhstan. The elevation ranges from less than −130 m in the Turpan Depression and the Karagiye Depression to over 7500 m in the Pamir and the Karakoram mountains. Due to the inland position and remoteness from oceans, central Asia is characterized by a semiarid and arid climate with low relative humidity and persistent soil moisture deficit making this region one of the most arid in the world.

As a result of climate change and human activities, central Asia is experiencing an increasingly severe deficit of water resources with increasing temperature and variable precipitation pattern [8]. Lands of Kazakhstan are featuring a variety of soils: a large part of the forest-steppe zone is occupied by chernozem (black earth), dark brown, light brown, and brown soils can be found to the south of this area. Desert and semidesert soils are represented in the form of sierozem (gray-brownish soil).

According to the [48], most of the territory of Kazakhstan is exposed to droughts of varying intensity. Scarcity and uneven distribution of water resources is amplifying the problem of droughts and leading to the emergence of sands over 30 million ha and saline land over 127 million ha.

The disturbance in the seasonal characteristics of soil formation under the effect of drought creates favorable conditions for land degradation. The weak formation of the land cover and vegetation also leads to desertification. These natural features of Kazakhstan are a cause for poor resilience of the environment to human pressures (it is estimated that about 75% of the country is under an increased risk of ecological destabilization). The analysis of the impact of adverse meteorological phenomena on agricultural crop yields showed that the share of atmospheric and soil drought is about 80%, heavy rain and hail—14%, frosts—2%, and severe frosts and strong winds—1% each [54,55].

The occurrence of drought in Kazakhstan is due to the specifics of the general circulation of the atmosphere. Drought sets almost across the entire territory of Kazakhstan, when the anticyclones of the Azores origin move from east to west, creating a band of high pressure covering the whole territory of the republic [56]. Geographic distribution of atmospheric droughts is demonstrated during invasion of arctic air from the north or from the northwest (from the water area of the Barents and Kara Seas) and the formation of a powerful anticyclone. If the Arctic air arrives from the Kara Sea to western Siberia, a stationary anticyclone is formed over central and eastern Kazakhstan. Consequently, an atmospheric drought is observed in the east of Kazakhstan. Western parts of the country at this time are exposed to cyclones. If the Arctic air invades from the Barents Sea in the western part of Russia, the center of a stationary anticyclone is located over the Urals. Accordingly, drought is observed in the west of the country. Particularly dangerous are blocking processes in the atmosphere. Due to the large area of Kazakhstan, the distribution of droughts on its territory and their intensity has considerable variation.

All three types of droughts occur in Kazakhstan: atmospheric, soil and atmospheric-soil, or general droughts. The latter represents the greatest threat to agricultural crops. The main feature of atmospheric drought is considered to be a stable anticyclone weather with long dry periods, high temperatures, and considerable air aridity. Soil drought occurs as a consequence of atmospheric drought, when deposits of soil moisture decrease rapidly during increased evaporation and become insufficient for the normal growth and development of plants.

Atmospheric-soil drought is a combination of conditions that characterize atmospheric and soil drought. They are divided into strong, medium, and weak in terms of intensity and area coverage. In terms of time of occurrence, droughts are separated into spring, summer, autumn, and winter [57]. Spring drought is characterized by low temperatures and low air humidity, low deposits of yielding moisture in the soil, and dry winds. Summer drought is usually distinguished by high fever and hot dry winds causing increased evaporation of water from soil and intensive transpiration of plants. Autumn drought is characterized by low temperatures and low deposits of yielding moisture in the root soil horizons. Winter drought occurs amid snow cover with a lack of moisture in the root soil horizons and at temperatures above 0 °C, when transpiration of some plants, especially winter crops, is resumed aggravated by sun and windy weather.

Figure 1 shows the calculated standardized moisture index (Kσ), which indicates the distribution of moisture zones in Kazakhstan.

About 94% of the territory lies in the following zones: hyperarid zone extends to the southern half of Kazakhstan, covering Atyrau, Mangystau, southern half of Aktobe, Zhambyl, Karaganda, and Almaty regions; arid zone extends along the northern part of Kazakhstan, as well as a narrow band in the southeast; semiarid zone is situated in the flat area of the north Kazakhstan and Akmola region, on the territory of the Tarbagatay mountain range, in low and middle mountain areas of the East Kazakhstan, South Kazakhstan, Almaty, and Zhambyl regions. Subhumid and humid zones have the essentially high elevations of the Tianshan Mountains. Combating desertification, therefore, presents an important challenge for Kazakhstan [58].

Especially strong droughts in Kazakhstan tend to arise in the negative phase of the MEI. Catastrophic droughts were observed in the Kostanay region in 1975, 1983, 1991, 1995, 1998, and 2004, coinciding with strong MEI and El Niño phases. A significant number of drought cases in the area are associated with positive NAO index values.

4. Methodology

The main purpose of the present study was to identify reliable predictors for wheat yields in Kazakhstan. To achieve this, we used an econometric model developed by the lead author at Environment Europe in Oxford using OxMetrics software. In modelling the dynamics of wheat yield in the Kostanay region of Kazakhstan we employed the cointegration and error-correction methodology, where explanation of the long-run changes is combined with short-term dynamics. To find adequate predictors that have a significant impact on drought formation and crop yield in Kazakhstan we analyzed the role of changes in soil moisture resulting from precipitation and temperature dynamics, cosmic radiation, and solar activity. First, an attempt was made to find a significant signal of the influence of the ENSO phenomenon in the formation of the weather regime in the northern regions of Kazakhstan.

The data for this research comes from the National Oceanic and Atmospheric Administration, the Ministry of Education and Research of the Republic of Kazakhstan, the Ministry of Agriculture of the Republic of Kazakhstan, Kazhydromet, Novosibirsk Neutron Monitor, NASA, and other sources.

The research shows that interannual climate variability in Kazakhstan is affected by several factors: El Niño and North Atlantic Oscillation (ENSO/NAO), Solar Irradiance (SI), and Galactic Cosmic Rays (GCR) [8,16,59].

In this paper we have chosen the Multivariate ENSO Index (MEI) based on the six main observed variables: sea level pressure (SLP), sea surface temperature (SST), zonal and the outgoing long wave radiation over the Tropical Pacific basin sky [60,61].

To study drought in Kazakhstan in this article we used the soil moisture index (K), which shows the ratio of annual precipitation to annual evaporation:

$$K=\raisebox{1ex}{$R$}\!\left/ \!\raisebox{-1ex}{${\overline{E}}_0$}\right.$$

(1)

R—sum of precipitation for the year, mm; ${\overline{E}}_0$—evaporation for the year, mm.

If the value of K > 1, then there is an excess of moisture, K < 1 indicates a drought.

To assess the severity of droughts we used the standardized moisture index (K_σ), calculated as the ratio of the anomaly to the standard deviation of K:

$$K_{\sigma }=\raisebox{1ex}{$\left(K_i-\overline{K}\right)$}\!\left/ \!\raisebox{-1ex}{$\sigma $}\right.$$

(2)

K_i—soil moisture index K in i-year; $\overline{K}$—average value of moisture index for the period from 1971 to 2011; σ—standard deviation of the soil moisture index.

The results obtained correspond to the conclusions of [20] who showed that both drought and aridity indicate an imbalance in water availability. While drought is a natural temporal hazard, aridity is a constant climatic feature.

Figure 2 shows ENSO MEI series from 1950 to 2017. During our study period between 1971 to 2011, MEI is characterized by considerable variability in time and predominance of negative values is observed in the last few decades, which means the warm phase of ENSO. Strong El Niño occurred in 1982–1983, 1991–1993, 1997–1998, and 2010. During these years, the value K_σ was below 0 and the Kostanay region of Kazakhstan observed droughts of varying intensity.

The amount of cloud cover and the intensity of Galactic Cosmic Rays (GCR) registered on the surface of the Earth were shown to be statistically correlated as demonstrated by satellite data [15,16]. Marsch and Svensmark [62] revealed that the interannual variability of ISCCP-D2 low-cloud properties observed in 1983–1994 was highly statistically correlated to both Galactic Cosmic Rays and El Niño–Southern Oscillation. It was shown as early as 1962 that cosmic radiation is likely to have an impact on plant growth [63]. This has been proven with recent data in a study focusing on tree growth in Britain [64]. The physics of Galactic Cosmic Rays—cloud cover interdependence can be explained as follows: at the height of 10–15 km from the Earth’s surface, the altitude of the tropopause separating the troposphere and the stratosphere, the Galactic Cosmic Rays reach the top of the clouds and lose a significant portion of their energy. Ref. [62] found that during two 11-year solar activity cycles, the correlation coefficient between the variations of the Galactic Cosmic Rays’ fluxes and the cloud cover on Earth was extremely high.

Based on the research outlined in the literature review, combining variables that are rarely found together, we have built a complex of econometric models, covering time series data between 1971 and 2009. Our approach comprised a simultaneous equations model of wheat yield and soil moisture in Kazakhstan as functions of a range of climatic and astronomical variables as well as the forecasting model capable of predicting likely changes in wheat yield into the future. The next section will deal with our econometric modeling experiments in greater depth.

5. Modelling Results

One of our hypotheses in this study was that soil moisture is affecting agricultural productivity under climatic change. Analysis showed that sufficiently good coherence was observed between temporal dependence of standardized moisture index and crop yield fluctuation, i.e., high values of Kσ corresponded to high yield values, and negative values to low-yielding years. These results are in agreement with the Second National Report of Climate Changes of the Republic of Kazakhstan [4]. We were able to model soil moisture based on the precipitation and temperature variables.

The second hypothesis that we tested was that a combination of El-Niño-related variables, characteristics of the solar activity, and cosmic radiation coupled with soil moisture will be enough to explain the variations in wheat yields. First, we estimated the short-run difference model of wheat yields. We then modelled the wheat yields dynamics in levels and combined the short-term and long-term elements in the yield model to arrive at a final equation. The fitted values and the data for wheat yields in Kazakhstan are presented in Figure 3.

The model coefficients can be seen in

Table 2. Econometric model of wheat yield dynamics in the Kostanay region, Kazakhstan.

Variables	Coefficient	Std. Error	t-Value	t-Prob	Part.R2
Change in YIELD	0.484347	0.07903	6.13	0.0000	0.5643
Change in Temperature [June]	−0.255337	0.1347	−1.90	0.0680	0.1103
Change in Total Solar Irradiancet−1	−6.48107	1.893	−3.42	0.0019	0.2880
Change in Wolft−1	0.0465606	0.01537	3.03	0.0051	0.2403
Cosmic Rays [December] t−1	0.00117928	2.894 × 10−5	40.7	0.0000	0.9828
Soil Moisture	1.18548	0.2816	4.21	0.0002	0.3792
MEI [June] t−1	−1.20924	0.3088	−3.92	0.0005	0.3459
ECM t−1	0.775867	0.09002	8.62	0.0000	0.7192

Source: author calculations.

.

At the same time, we modeled soil moisture separately in an additional equation based on the connections with precipitation and temperatures. Figure 4 presents the performance of our model for soil moisture. It is interesting to note the years that came up as significant outliers in the model equation: 1987, 1991, 1993, 1994, and 2005. During the years of 1987 and 1991 moderate El Niño events were observed; 1994 and 2005 were the years of weak El Niño.

It is shown in

Table 3. Econometric model of soil moisture dynamics in the Kostanay region, Kazakhstan.

Variables	Coefficient	Std. Error	t-Value	t-Prob	Part.R2
Precipitation [June]	0.00348383	0.0006696	5.20	0.0000	0.4915
Precipitation [July]	0.00358584	0.0006185	5.80	0.0000	0.5456
Precipitation [August]	0.00416834	0.0008602	4.85	0.0000	0.4561
Precipitation [October]	0.00316804	0.001026	3.09	0.0045	0.2539
Temperature [July]	−0.118263	0.01035	−11.4	0.0000	0.8233
I1987	1.45623	0.3843	3.79	0.0007	0.3389
I1991	−0.998864	0.3810	−2.62	0.0140	0.1971
I1993	1.18435	0.3848	3.08	0.0046	0.2528
I1994	1.22896	0.3822	3.22	0.0033	0.2697
I2005	−1.03599	0.3944	−2.63	0.0138	0.1977

Source: author calculations.

that soil moisture depends on precipitation in June, July, August, and October and temperatures in July.

Finally, we obtained a resulting two-equation model (

Table 4. Simultaneous equation for soil moisture.

Variables	Coefficient	Std. Error	t-Value	t-Prob
Precipitation [June]	0.00334953	0.0006376	5.25	0.0000
Precipitation [July]	0.00402316	0.0005775	6.97	0.0000
Precipitation [August]	0.00378489	0.0007966	4.75	0.0001
Precipitation [October]	0.00278728	0.0009193	3.03	0.0069
Temperature [July]	−0.117935	0.01046	−11.3	0.0000
I1987	1.37304	0.3439	3.99	0.0008
I1991	−0.750801	0.3380	−2.22	0.0387
I1993	1.23502	0.3581	3.45	0.0027
I1994	1.37440	0.3549	3.87	0.0010
I2005	−0.851416	0.3532	−2.41	0.0262

Source: author calculations.

and

Table 5. Simultaneous equation for wheat yield.

Variables	Coefficient	Std. Error	t-Value	t-Prob
Change in YIELD	0.512376	0.07265	7.05	0.0000
Change in Temperature [June]	−0.307918	0.1263	−2.44	0.0247
Change in Total Solar Irradiancet−1	−6.50333	1.767	−3.68	0.0016
Change in Wolft−1	0.0508515	0.01450	3.51	0.0024
Cosmic Rays [December] t−1	0.00118554	3.080 × 10−5	38.5	0.0000
MEI [June] t−1	−1.29832	0.2851	−4.55	0.0002
Soil Moisture	1.32344	0.2979	4.44	0.0003
ECM t−1	0.834856	0.08633	9.67	0.0000

Source: author calculations.

) capable of explaining the dynamic of wheat yields and the soil moisture in response to changes in climatic variables.

The final two-equation model (Figure 5) is composed of the two interrelated equations for soil moisture and wheat yield changing over time. The first equation explains variation in soil moisture through precipitation in June, July, August, and October (moisture tends to increase with higher precipitation levels), July average temperatures (moisture tends to decrease with higher temperatures), strong positive impacts in 1987 (moderate El Niño), 1993 (weak El Niño), and 1994 (moderate El Niño), and negative impacts in 1991 (strong El Niño year) and 2005 (weak El Niño year). All variables in this equation are statistically significant at levels below 1% and only two variables, I1991 and I2005, are significant at 5%.

The second equation explains variation in wheat yields through a range of variables: there is an acceleration effect of a change in yield (with a positive coefficient), change in average June temperatures (yields tend to be lower with higher June temperatures), and change in total solar irradiance in the previous year (yields tend to be lower the greater the change); change in the Wolf number of the previous year (higher number of sunspots tends to increase yields); December amount of cosmic rays hitting the Earth (the higher the amount of cosmic rays hitting the Earth, the higher the yields); MEI in June of the previous year (yields tend to be lower if El Niño was observed in June of the previous year and higher if La Niña was observed in the same period), soil moisture (higher moisture in the soil tends to produce higher yields), and an error-correction term. It should be added that most variables are statistically significant at a level lower than 1% with one exception of the change in June temperatures, where it is significant at 5%.

The visual observation of the moisture and yield curves in Figure 5 points to the fact that the model covering the 36 observations between 1973 and 2008 fits the data well. The explanatory model has been transformed into a forecasting model via ARIMA time series modelling and final integration. The resulting model is capable of forecasting fluctuations in yields up to the year 2030. As Figure 6 illustrates, we should have expected a drop in yields in 2013 to a level of just above 7 t/ha and should expect gradual increases to a maximum yield level of approximately 11 t/ha in 2019, after which the yields continue fluctuating at around 10 t/ha.

In our final model, wheat yield depends on soil moisture; there is a significant El Niño response in the form of MEI6_t−1, and it also depends on the CR12_t−1 (cosmic rays in December of the previous year). There is a significant response to the changes in TSI_t−1 (total solar irradiance) and the changes in solar activity (Wolf_t−1) as well as changes in the temperatures in June and changes in wheat yield itself.

The estimation sample is: 1973–2008.

We can see that most model specification tests have been satisfied apart from the slight issues with autoregression in individual equations. All vector tests for normality of residuals, autocorrelation of residuals, and residual homoscedasticity have been met.

6. Discussion

Our results reflect the most recent discoveries in the climatic drivers for crop yields. We have identified several effects to be significant in the model of wheat yield in Kazakhstan; most importantly, the connection between El Niño and wheat yield, the role of changes in cosmic radiation, solar activity, and soil moisture, resulting from precipitation and temperature dynamics. Our paper has confirmed a link between cosmic rays, solar activity, and wheat yield, which are rarely studied together with traditional climatic variables such as precipitation and temperatures. The model could be modified further so that the individual components could be forecasted into the future using various forms of a time series ARIMA model. The resulting integration of these forecasts allows the prediction of wheat yields into the future (Figure 6).

This model will have a significant practical value for predicting crop yields in Kazakhstan under changing climatic conditions and will be tested further as new data becomes available. Currently it can forecast wheat yields into the future up to 2030 under the “ceteris paribus” assumption (everything else being equal). Both our hypotheses have therefore been confirmed.

7. Conclusions

This paper presents recent research highlighting the impacts of climatic factors on droughts and yields of wheat in northern Kazakhstan. Drought is one of the major manifestations of climate change in Kazakhstan and exerts a considerable pressure on the agricultural sector, which is the most vulnerable sector to drought in the economy. To mitigate its impact, a warning system to deliver timely information about the possible occurrence of droughts is needed. This is possible only if the influence of large-scale atmospheric processes, formed due to the impact of cosmic-geophysical conditions and interdependent factors in the ocean–glacier–land-atmosphere systems with direct and reverse links, is taken into account.

Climatic change is clearly a process determined by the complex interaction of natural and anthropogenic causes. Without explaining the exact mechanism of anthropogenic impacts on the climatic system, we explore the impact of a range of indicators, some of them experiencing clear anthropogenic pressures, to explain the variability in wheat yields in Kazakhstan.

In modelling the dynamics of wheat yield in the Kostanay region of Kazakhstan we employed the cointegration and error-correction methodology, where explanation of the long-run changes has been combined with short-term dynamics. The results can be summarized as follows.

We have conducted the zoning of the territory of Kazakhstan based on the standardized moisture index. It was demonstrated that the region under study, the Kostanay region, is characterized by almost all zones from arid to moist. Most parts of Kazakhstan, about 94% of the territory, lie in two zones: arid and semiarid. Catastrophic droughts were observed in the Kostanay region in Kazakhstan in 1975, 1983, 1991, 1995, 1998, and 2004, coinciding with strong MEI and El Niño phases.

Our analysis showed that sufficiently good coherence was observed in the relationship between temporal variation of standardized moisture index and crop yield fluctuation, i.e., high values of moisture index correspond to high yield values, and negative values to low-yielding years.

The article found an important connection between the Multivariate ENSO Index (MEI) and wheat yields in Kazakhstan. Our research found significant correlation between Galactic Cosmic Rays (GCR) received on Earth and wheat yield in the Kostanay region of Kazakhstan. If the impact of the cosmic rays on droughts is real, the mechanism is attributed to the ionization caused by GCR influencing droughts through the microphysics of the cloud formation. Our research has confirmed the effect of total solar irradiance (TSI) on wheat yield.

The econometric model proposed in this paper explains changes in wheat yield as a function of changes in soil moisture, changes in June air temperatures, changes in total solar irradiance, changes in TSI_t−1 (total solar irradiance), and changes in solar activity (Wolf_t−1). In our model, there is a significant El Niño response in the form of MEI6_t−1, and wheat yield also depends on the CR12_t−1 (cosmic rays in December of the previous year) and changes in wheat yield itself.

Our results reflect the most recent discoveries in the climatic drivers for crop yields. We have identified several effects to be significant in the model of wheat yield in Kazakhstan, most importantly, the connection between El Niño, the role of changes in cosmic radiation, solar activity, and soil moisture, resulting from precipitation and temperature dynamics. The model proposed in our paper could be modified further so that the individual components could be forecasted into the future using various forms of a time series ARIMA model. The resulting integration of these forecasts allows the prediction of wheat yields into the future, thereby building resilience and strengthening food security in the face of a changing climate.

One of the potential limitations of this study stems from the fact that Kazakhstan has vastly different climatic conditions in various regions and new modeling experiments would need to be conducted to confirm the found dependencies for other regions.

Our future research plans include more in-depth analysis of the climate–water–food nexus in relation to other sensitive agricultural crops experiencing significant impacts of climate change, such as rice, coffee, etc.