Abstract
Climate change constitutes a significant threat to agricultural productivity, particularly in India, where nearly half of the population relies on agriculture for their livelihoods. The sector is increasingly exposed to rising temperatures, irregular precipitation patterns and the intensification of extreme climatic events. While previous studies have identified climate-induced yield reductions in staple crops, much of the existing evidence is based on short-term analyses and lacks a comprehensive evaluation across diverse crop groups. This study addresses this limitation by employing a panel dataset covering 20 crops across 20 Indian states from 1966 to 2016 to assess both historical impacts and projected outcomes under shared Socioeconomic Pathways (SSP) climate scenarios. The analysis employs crop-specific, non-linear panel regression models that incorporate key climatic variables, including temperature, precipitation, wind speed and evapotranspiration, alongside critical non-climatic factors such as mechanisation, irrigation and fertiliser application. The findings indicate that crops like rice and wheat are most adversely affected by rising temperature, whereas crops like pearl millet and cotton exhibit greater resilience, though cotton remains highly sensitive to precipitation variability. Pigeon pea and sesamum benefit from increased evapotranspiration up to a threshold, beyond which yields decline. By the end of the twenty-first century, the combined effects of rising temperature and precipitation are projected to substantially reduce the yields of pigeon pea and maize. The adoption of High-Yielding Varieties (HYVs) demonstrates notable productivity gains, especially for pearl millet and rice. These findings underscore the urgent need for targeted adaptation strategies, including climate-resilient crop varieties, water-saving technologies and diversified cropping systems to enhance agricultural resilience.
JEL
1. Introduction
Climate change poses one of the most formidable threats to global agriculture, with projected productivity declines of up to 15.9 per cent by the 2080s and even sharper drops of 19.7 per cent in developing countries if global warming persists (IPCC, 2022; Lal, 2011; Ortiz-Bobea et al., 2021). In India, while the Green Revolution achieved food grain self-sufficiency, it triggered environmental costs like soil degradation and groundwater overuse (Harwood, 2021; Lele & Goswami, 2017; Eliazer Nelson et al., 2019), now worsened by climate change. Rising temperatures, heat waves, droughts, heavy rains and cyclones threaten Indian agriculture, which employs 42.86 per cent of the country’s workforce (Arora, 2019; Birthal & Hazrana, 2019; Singh & Singh, 2017; World Bank, 2024). Extreme weather caused crop losses worth 0.25 per cent of India’s GDP in 2018–2019 (Singh et al., 2020). With the global population projected to reach 11 billion by 2100 (Ruane & Rosenzweig, 2018), it is imperative to sustain and increase yields. Therefore, it is crucial to investigate the specific ways in which climate change affects agricultural productivity to inform strategies that will ensure food security. Hence, this study addresses the research question of how climate change has affected agricultural productivity in India over the past five decades, focusing on non-linear, diverse crop-specific impacts and their potential change under future climate scenarios.
Projections indicate that climate change could severely affect global crop yields, with end-of-century losses estimated for wheat (9–12 per cent), corn (5.6–6.3 per cent), rice (18.1–19.4 per cent) and soyabeans (15.1–16.1 per cent) (Cui, 2020; Hussain et al., 2020). In India, alterations in temperature and precipitation regimes are anticipated to have profound implications for agricultural productivity. Over the past five decades, temperatures have exhibited a persistent upward trajectory, while long-term precipitation trends have remained relatively stable, though marked by notable interannual variability (see Figure 1). These trends could reduce rice and wheat yields by 8.1 per cent and 6.51 per cent by 2080, respectively (Wang et al., 2022). Rising temperatures may disrupt crop growth, intensify water stress, and heighten pest and disease pressures. Although precipitation remains stable, it may not offset these effects, especially as rainfall during the rice-growing season has declined by 7 per cent in the last 30 years (Daloz et al., 2021). Without adaptation, rain-fed rice yields could fall by 20 per cent by 2050 and by 47 per cent by 2080, while irrigated rice may see declines of 3.5 per cent and 5 per cent (Wang et al., 2018). Understanding these impacts necessitates examining both climatic and non-climatic factors, as temperature and rainfall directly influence crop productivity (Anderson et al., 2020; Aryal et al., 2020; Kumar & Sharma, 2022).

Changes in temperature and precipitation, along with shifts in evapotranspiration and wind speed, have altered agro-ecological zones, affecting cereals, pulses, oilseeds and cash crops (Meena et al., 2015). In some cases, crop responses vary widely; for instance, rice yields rise by 5.21 per cent with increased rainfall but drop with higher wind and evapotranspiration; bajra benefits from temperature rise but declines with evapotranspiration; chickpea, groundnut and sugarcane also show similar sensitivities (Mohapatra et al., 2022). While most research centres on staples like rice and wheat, according greater attention to pulses, oilseeds and cash crops is vital for India’s food security, economic stability and rural livelihoods. A broader focus is essential to address the future implications of climate change across this wide array of crops through inclusive and targeted adaptation strategies.
Agricultural reforms in the 1990s, supported by better access to advanced seeds, fertilisers and financial resources, significantly improved crop yields across India. These policy reforms led to major increases in cereal production in northern India and oilseed output in the south (see Figure 2). 1 Key inputs like improved seeds, irrigation and mechanisation played a critical role in boosting productivity and reducing labour dependence (Gulati & Juneja, 2020). These input-driven gains contributed to nationwide agricultural growth (Nadkarni, 2022). However, the long-term sustainability of these gains remains uncertain, highlighting the need to understand how agricultural inputs interact with climate variability and external stresses.

Many studies examining the intersection of climate change and agriculture rely on temporal scopes of less than 30 years. While these shorter time frames effectively capture near-term climate variability and extreme weather events, they often miss the cumulative effects of gradual climate change on agronomic systems. For instance, Chen and Gong’s (2021) 35-year panel study, Chandio et al.’s (2020) 32-year analysis and Da et al. (2024) offer valuable insights into adaptation patterns, yet remain limited in capturing long-duration climatic shifts. Extending the temporal scale of analysis is vital to uncover long-term resilience thresholds and adaptive mechanisms as agricultural systems confront escalating climate risks.
India’s pivotal role in global food security underscores the importance of understanding how climate change affects its agricultural productivity. With the global population projected to reach 11 billion by 2100, sustaining food supply is crucial (Patel et al., 2020). Under high-emissions scenarios like SSP5-8.5, monsoon rainfall could increase by 24.3 per cent; around 0.33 mm/day or 5.3 per cent per °C of warming by the end of the century (Katzenberger et al., 2021; Yadav et al., 2024). These shifts may severely disrupt rice and wheat yields (Cui & Xie, 2022). To ensure resilience, it is equally vital to assess climate impacts on pulses and oilseeds, which are critical for dietary diversity and nutritional security.
Projected climatic changes in India include temperature increases of 2.3°C–3.9°C and rainfall gains between 61.5 mm and 98.2 mm by the 2080s. Eastern regions may see even sharper changes, with temperature rises of 2°C–4°C and rainfall increasing by up to 15 per cent, while coastal areas could face wind speed reductions of about 1 m/s (Karmacharya et al., 2017; Mishra et al., 2022). This regional variability necessitates adaptive strategies tailored to India’s diverse agro-ecological zones. Our study builds on these insights by integrating Intergovernmental Panel on Climate Change (IPCC) projections across multiple Shared Socioeconomic Pathways (SSP1-2.6 to SSP5-8.5), offering forward-looking analysis of crop productivity in India. While earlier studies, such as Choudhary et al. (2020), have analysed climate impacts on staple crop yields using extensive historical data over six decades, they have largely overlooked forward-looking projections that consider evolving climate scenarios. By bridging this gap, our study offers a critical framework to aid policymakers and stakeholders in anticipating and mitigating risks to agricultural productivity under future climate conditions, addressing both current impacts and long-term adaptation needs.
This research advances existing literature by addressing several key gaps. First, it undertakes a longitudinal analysis of crop–climate interactions using a panel dataset that spans 50 years (1966–2016), a period marked by profound transitions in Indian agriculture, including the Green Revolution, successive policy shifts and technological innovations. This extended time frame enables a robust assessment of how crop yields have responded to both climatic variability and non-climatic pressures over time. Second, the study expands the scope of prior research by including a diverse set of crops across four major categories, such as cereals, pulses, oilseeds and cash crops, many of which are typically overlooked in empirical assessments. This broad coverage allows for a detailed examination of heterogeneous crop responses to climate change, offering comparative insights that are essential for designing targeted adaptation strategies. Third, the analysis integrates forward-looking yield projections based on SSP climate scenarios, thereby providing a future-oriented perspective on the potential impacts of climate change on Indian agriculture. Fourth, the study employs advanced econometric techniques, including non-linear panel regression models, to address methodological challenges such as endogeneity, omitted variable bias and heterogeneity in crop response functions. Finally, by identifying empirically grounded yield thresholds, the research offers crop-specific and policy-relevant adaptation strategies that can inform agricultural planning and resilience-building efforts in a changing climate.
The structure of the paper is as follows. After the introductory Section 1, Section 2 outlines the theoretical background of the study. Section 3 details the data sources and the empirical methods used for the analysis. Section 4 presents the results and discussion, while Section 5 concludes with policy implications.
2. Theoretical Framework
To analyse the impact of climate change on Indian agricultural crops, this study conceptualises agricultural production through a weather-augmented framework where climatic variables act as environmental efficiency shifters, where the crop yield Yit for state i at time t is determined by a vector of non-climatic inputs Xit and a set of climatic variables Wmit. Unlike conventional inputs that typically exhibit monotonic returns, environmental variables such as temperature and rainfall follow a non-linear pattern, contributing to productivity up to an optimal threshold before becoming detrimental. To capture this multiplicative relationship where climate scales the effectiveness of all other inputs, we model the production function using an exponential specification. In this framework, the production function is defined as follows:
Here, the term exp(a) captures the Total Factor Productivity (TFP) characterised by constant returns to scale, while exp(XitØ) represents the contribution of non-weather control variables Xit. The coefficients Ø represent the semi-elasticities, indicating that a one-unit change in non-climatic inputs results in a proportional percentage change in output. The term X it represents the net environmental efficiency derived from climatic conditions. To account for the non-linear impacts of weather, we specify X it as a quadratic function of the m weather variables.
The coefficients γm and θm capture the linear and quadratic effects of the mth weather variable, respectively. This form allows the marginal effect of the weather variable Wmit to vary with its intensity, capturing both the beneficial and adverse deviations from ideal conditions.
By taking the natural logarithm of this production function, we derive the semi-logarithmic specification employed in our empirical analysis. The resulting equation is
Equation (3) corresponds directly to Equation (6) in Section 3. Furthermore, the inclusion of the idiosyncratic error εit and unobserved random variables uit [as detailed in Equation (7)] accounts for stochastic shocks to production, ensuring that the theoretical model is fully consistent with the econometric strategy of minimising bias from unobserved heterogeneity.
To further elaborate the relationship between climate and yield, we examine the marginal impact of climatic variables. The marginal effect of a specific weather variable Wmit on the log of crop yield is given by the partial derivative with respect to Wmit:
Equation (4) underscores that the impact of climate is not constant but depends on the current level of the weather variable. We anticipate a concave relationship where θm < 0, indicating that the marginal benefits of weather variables diminish as they increase. This allows us to identify the optimal climatic condition, or ‘inflection point’, where yield is maximised. By setting the marginal effect to zero, we can solve for the optimal level of the weather variable W*
mit
:
Values of Wmit below this threshold indicate that increases in the weather variable are beneficial, while values above it suggest that further increases lead to yield losses.
3. Materials and Methods
3.1 Data Sources
This study utilises a longitudinal dataset covering a 50-year period from 1966 to 2016, drawing on ICRISAT (2024) data from 20 major agricultural states in India that collectively represent the principal contributors to national food grain production. The analysis focuses on 20 key crops, which are categorised into cereals, pulses, oilseeds and cash crops (see Appendix B, Table B1). In order to ensure data consistency and suitability for analysis, comprehensive data cleaning techniques were employed, resulting in a balanced panel dataset. This preparation is essential for conducting robust panel regression analysis and deriving reliable empirical insights.
Climatic variables, including temperature, precipitation, wind speed and evapotranspiration, were sourced from the Centre for Monitoring Indian Economy (CMIE, 2024), which provides consistent annual state-level data. Non-climatic variables, such as irrigation coverage, mechanisation and fertiliser use, were obtained from crop-specific data available through ICRISAT (2024). Since crop-wise fertiliser consumption data are not directly available, estimates were derived by prorating total fertiliser use across crops in accordance with their respective shares of cultivated area, following the method employed by Verma et al. (2020). In the absence of crop-disaggregated labour data, total agricultural labour in each state was proportionally allocated across crops based on their cultivated area. Data on electric pump sets and tractors were sourced from ICRISAT VDSA’s (2024) agricultural census. Due to the lack of crop-specific information for these inputs, the same proportional distribution method was applied, in line with the strategies adopted by Gupta et al. (2014), Padakandla (2016) and Mohapatra et al. (2022). This comprehensive methodological approach ensures consistency in estimating crop-level input intensities across states and over time.
Scatter plot analysis confirms a non-linear relationship between crop yields and various climate variables, 2 revealing a possible quadratic association, as shown in Figure 3. While certain climate variables show strong correlations with yields, the relationship remains non-linear. For this study, two primary climate variables, temperature and precipitation, are emphasised, as these will serve as the basis for projecting crop yields under future climate scenarios.

3.2 Empirical Methodology
This study employs a semi-log production function due to its econometric tractability and interpretability in capturing percentage changes in output. The functional form for state θ and time t is represented as: yit = f(Xit, Wit). After accounting for all variables, the production function becomes
Here, α captures the intercept in a production function characterised by constant returns to scale, implying linearity and homogeneity of degree one. Yit represents the total crop yield for state i at time t. Xjit are non-weather control variables, coefficient Øj are semi-elasticities; a one-unit change in Xjit implies an approximate 100 Øj per cent change in Yit. Wmit are weather variables (e.g., temperature and rainfall), γm and θm capture the linear and quadratic effects of the mth weather variable, allowing non-linear responses. dt are year fixed effects, xit are state-specific linear time trends, and fit is the idiosyncratic error. The analysis applies crop-specific regressions to capture the nuanced responses to climatic and non-climatic factors.
The adoption of the semi-log functional form, as indicated by Sargan’s criteria (see Appendix B, Table B3), rejected the null hypothesis, supporting this structure’s necessity (see Appendix A, Note A1 for more details). We analyse the relationship between crop yield and both weather and non-weather factors using panel regression models, with several pre-estimation tests to ensure model robustness. The Fisher unit-root test (see Table B2) confirmed variable stationarity, while the Variance Inflation Factor (VIF) test (see Table B5) ruled out multicollinearity concerns.
To address unobserved random variables like farm inputs, denoted as γ, the regression equation is further expressed as (see Note A2):
However, applying the Ordinary Least Squares (OLS) method leads to bias and inconsistency in the βˆ estimator, as f(x, c) ≠ 0, indicating a correlation between unobserved variables and independent variables. The Hausman test suggests that γ follows a random effects pattern across states, ruling out the pooled OLS method. Additionally, the Pesaran test (see Table B4) confirmed significant cross-sectional dependence, while the Wooldridge test identified first-order autocorrelation, indicating both cross-sectional and serial autocorrelation across models (see Note A3).
Although the Feasible Generalised Least Squares (FGLS) method is widely used for addressing serial correlation and cross-sectional dependence, it is often limited in accurately capturing heteroscedasticity in crop-specific regressions (see Table B9). In contrast, the Beck–Katz estimator, which employs Panel Corrected Standard Errors (PCSE), is more effective in adjusting for these issues without transforming the original data structure. PCSE provides consistent and reliable standard error estimates, particularly in datasets with a large number of cross-sectional units and relatively short time dimensions. As highlighted by Famanta et al. (2024), PCSE offers greater efficiency and accuracy in estimating model parameters in such cases, making it more suitable than FGLS for robust inference. Accordingly, the final model adopts the PCSE methods and is structured as follows:
In this context, i indicates the total number of panels, while t corresponds to the time frame of the study within each panel. The symbol ε refers to the idiosyncratic error term, which may exhibit autocorrelation over time (t) or contemporaneous correlation across panels (i).
The error covariance matrix Ω
NT
is structured as Ω
NT
= ΣX π:
Here, Σ is the error covariance matrix and π is a temporal correlation matrix.
The PCSE estimator is computed as follows:
Here, x~ and y~ represent the Prais-transformed observations of the explanatory and dependent variables, respectively. The PCSE method accounts for both cross-sectional dependence and serial autocorrelation, making it well-suited for our analysis of crop yields. The PCSE estimator for the covariance matrix is
In this study, we applied PCSE regression to account for both cross-sectional dependence and serial autocorrelation, resulting in more reliable forecasts of future crop yields under changing climate conditions. To assess the marginal effects of climatic variables such as temperature and precipitation, we calculate partial derivatives based on the estimated panel regression coefficients. The net impact of these weather variables on crop yield is derived by summing their respective marginal effects. The equation for this marginal effect is expressed as follows:
(See Appendix A, Note A4 for more details.)
To assess crop tolerance under different climate scenarios, we determine the inflection points, thresholds beyond which further changes in weather may reduce yields. If θm ≠ 0; the turning point (value of Wmit that maximises/minimises lnYit) is:
With the methodology established, Section 4 delves into the historical impacts of climate variables on crop yields in India, highlighting the effects across various crop categories, including cereals, pulses, oilseeds and cash crops.
4. Results and Discussion
The following three sub-sections present the empirical results. Section 4.1 examines historical climate impacts on crop yields, comparing climatic and non-climatic effects across crop categories. Section 4.2 explores future yield projections under alternative IPCC climate scenarios, analysing crop-specific responses to combined increases in temperature and precipitation over the twenty-first century. For robustness, Section 4.3 compares the estimation results from the FGLS and PCSE methods, respectively.
4.1 Historical Climate Effects on Crop Yields
4.1.1 Cereals
The elasticity analysis (see Appendix B, Table B8) for cereals highlights the complex interplay between yield, climate variables and non-climate factors. Temperature emerges as a major determinant, with rice and wheat exhibiting strong negative elasticities of −24.58 and −17.50, respectively, indicating that yield declines as temperatures rise. Wheat appears particularly vulnerable to heat stress, with a critical threshold at 17.63°C beyond which yields decrease sharply. This trend aligns with the study by Jyothsna et al. (2024), who report that high night temperatures (>25°C) can reduce rice yield by 10–15 per cent, while day temperatures (>35°C) can cause a 10–20 per cent decline due to increased respiratory losses. In contrast, pearl millet and finger millet show positive elasticities of 6.78 and 20.77, respectively, supporting their characterisation as climate-resilient, next-generation crops (Neupane et al., 2022). However, their squared terms reveal that the positive impact tapers off beyond certain thresholds.
Wind speed also significantly influences cereal yields, though effects vary across crops. Rice benefits from moderate wind, showing a positive elasticity of 6.50, possibly aiding pollination. In stark contrast, wheat is highly sensitive, displaying a sharp negative elasticity of −77.54, likely reflecting wind-induced physical damage. Sorghum appears more resilient (see Figure 4), though its squared term indicates a turning point at 2.05 m/s, beyond which wind becomes detrimental. This resilience strengthens the view of DeFries et al. (2023) that sorghum could be a climate-resilient substitute for wheat under future climate scenarios.

The role of precipitation’s role is equally critical, especially when considering its interaction with irrigation. Rice shows a negative elasticity of −0.42, suggesting susceptibility to excess rainfall and related risks like waterlogging. However, Ashkra et al. (2023) find a positive link between rainfall and rice yield, indicating a 1.58 per cent increase in yield per additional millimetre of rainfall, highlighting the urgent need for temperature-tolerant, water-efficient rice varieties. Maize and sorghum, with positive elasticities of 0.0023 and 0.3197, respectively, thrive under increased rainfall, particularly in rain-fed regions. For finger millet, precipitation benefits decline beyond 1,000 mm.
Among the non-climate variables, fertiliser application positively affects yields across cereals, except in the case of wheat, highlighting the importance of nutrient inputs. Mechanisation, particularly tractor usage, strongly enhances rice and sorghum yields, with rice showing a regression coefficient of 5.99 (see Table B6). Irrigation infrastructure, reflected by electric pump set usage, also proves critical. Rural literacy shows mixed effects, likely boosting yields through better adoption of technology. However, labour elasticity varies, influenced by differing mechanisation levels. The spread of High-Yielding Varieties (HYVs) significantly improves yields, especially in the case of rice and pearl millet, underlining the essential role of improved seed technologies in strengthening climate resilience.
4.1.2 Pulses
Chickpea shows a modest positive temperature coefficient of 0.2305, suggesting a slight benefit from warmer conditions when only climate factors are considered (see Table B7). This positive effect grows significantly with the inclusion of non-climate variables, such as irrigation, where elasticity rises to 10.32. This finding aligns with the study by Neelima and Kumar (2023), who emphasise that technological advancements can enhance crop resilience to climate variability. However, Bera (2021) highlights that erratic rainfall and extreme weather introduce challenges for pulse production, such as delayed sowing and heat-stress-induced yield losses, potentially offsetting temperature-related gains.
For groundnut, higher temperatures initially result in a negative coefficient of −0.4200 under the climate-only model, indicating heat vulnerability. Yet, when non-climate factors are considered, elasticity becomes positive (23.06), suggesting that management practices can effectively mitigate temperature stresses. This supports the view of Gull et al. (2020), who found that pulses are highly sensitive to heat stress, particularly during flowering, though targeted interventions can help buffer some impacts.
Wind speed generally negatively affects pulse crops, with chickpea and pigeon pea showing strong negative elasticities of −51.51 and −58.24, respectively, likely due to physical damage or increased evapotranspiration. Groundnut is less affected (−26.92) but remains vulnerable. Precipitation initially benefits groundnut yields (0.00077), but when non-climate factors are introduced, the precipitation deviation coefficient shifts to −2.90E−04 (see Figure 5), suggesting that excess rainfall can become detrimental under compounding stressors.

Other factors, such as rural literacy and agricultural labour availability, positively influence pulse yields, though minor pulses exhibit slight declines (−0.89), reflecting varying labour needs across pulse varieties.
4.1.3 Oilseeds and Other Seeds
Temperature affects oilseed yields in diverse ways, underscoring the varied sensitivities of different crops to climatic changes. Castor shows a strong positive elasticity of 56.48, indicating its ability to thrive under warmer conditions, likely due to its adaptability and heat tolerance (see Appendix B, Table B8). In contrast, rapeseed and mustard, and linseed display negative elasticities of −6.68 and −7.61, respectively, suggesting that higher temperatures inhibit their productivity. The case of rapeseed and mustard is particularly noteworthy, as the significant negative square term of −0.0100 (see Figure 6) implies that temperatures exceeding an optimal 21°C lead to a sharp reduction in yields. These results align with Pillai and Walia (2024), who emphasise the heightened vulnerability of mustard to climate change, with both temperature and rainfall negatively impacting its resilience.

Wind speed presents distinct challenges across oilseed varieties. Sesamum records a negative elasticity of −14.21, reflecting its sensitivity to high wind speeds that can elevate evapotranspiration and cause physical damage. Conversely, castor demonstrates a positive elasticity of 51.08, suggesting greater resilience, possibly due to its robust physiology or cultivation practices in more sheltered environments. Moreover, Maity et al. (2023) highlight that self-pollinated crops, such as several oilseeds, suffer greater yield reductions under temperature stress, indicating that wind may amplify yield variability.
Evapotranspiration further influences yields in complex ways. Sunflower exhibits a turning point at 828 mm, with an elasticity of −4.06 beyond this threshold, while safflower, despite a positive elasticity of 231.60, faces declining yields beyond 1,071 mm of evapotranspiration.
Precipitation effects are equally nuanced. Castor benefits from rainfall, showing a positive elasticity of 0.95, whereas sesamum and linseed respond negatively, with elasticities of −0.07 and −0.25, respectively. Notably, sesamum yields sharply decline beyond 1,759.25 mm of precipitation.
Socio-agronomic variables also play critical roles. Irrigation significantly boosts sesamum yields (elasticity 23.65), though excessive irrigation reduces rapeseed and mustard productivity (−19.48). Fertiliser use enhances safflower yields (elasticity 0.0110), highlighting the importance of nutrient management.
4.1.4 Cash Crops
Sugarcane and cotton display distinct sensitivities to temperature. For sugarcane, a 1 per cent increase in average temperature boosts yield by 7.01 per cent, indicating a strong positive response. Cotton similarly exhibits a positive elasticity of 7.65 per cent, though the effect tapers at higher temperatures, suggesting vulnerability to extreme heat. This pattern aligns with Jans et al. (2021), who project a 50 per cent increase in cotton output under RCP8.5 by 2,099, despite lower water availability due to CO2-driven effects. Farooq et al. (2023) also warn that rising temperature and precipitation variability threaten the resilience of cotton.
Evapotranspiration similarly affects these crops in nuanced ways. For sugarcane, the positive influence of evapotranspiration declines beyond 782.05 mm, as indicated by negative squared terms (see Figure 7). Cotton’s critical threshold is higher, at 1,291.67 mm, after which further evapotranspiration reduces yields.

Precipitation presents a dual impact. Sugarcane benefits from rainfall up to 3,750 mm, beyond which yields decline, likely due to waterlogging or nutrient leaching. Cotton, conversely, shows a negative elasticity of −0.24 per cent, highlighting its sensitivity to excess moisture.
Non-climatic factors are equally crucial. Irrigation supports both sugarcane and cotton crops, though sugarcane experiences diminishing returns when combined with heavy rainfall. Cotton benefits from irrigation but requires careful moisture regulation. Mechanisation also shows contrasting effects: sugarcane yields fall with tractor use (elasticity of −121.74), while cotton gains significantly from mechanisation, with a positive elasticity of 1,009.17, reflecting its greater dependence on efficient mechanised practices.
4.2 Future Climate Scenarios
To analyse the future projections of crop yields under climate change scenarios, we use findings from SSP1-2.6 for precipitation and SSP5-8.5 for temperature across three periods: 2021–2040, 2041–2060 and 2081–2100 (see Figure 8). 3 These projections highlight varying levels of crop resilience and sensitivity to climatic changes, revealing both adverse and favourable trends.

Among cereals, wheat is the most negatively impacted, with yield losses increasing from 18.74 per cent to 27.35 per cent by the end of the century. This sharp decline is consistent with Daloz et al. (2021), who reported wheat yield reductions ranging from 4 per cent to 36 per cent under rising temperatures and reduced irrigation. In contrast, pearl millet displays resilience, with projected yield gains ranging from 17.99 per cent to 24.09 per cent, suggesting its adaptability under warming and erratic rainfall conditions. Rice experiences moderate yield reductions between −3.51 per cent and −8.71 per cent, particularly in later periods, supporting Singh et al. (2024), who found similar declines in rice yields in Uttar Pradesh under comparable SSP scenarios. Maize exhibits significant vulnerability, with yields falling from −12.78 per cent to −28.44 per cent, highlighting its sensitivity to temperature extremes. Conversely, finger millet performs better, with projected yield gains increasing from 12.57 per cent to 20.78 per cent. These contrasts suggest that future climatic conditions may favour pearl millet and finger millet while posing challenges for wheat, rice and maize cultivation.
The pulses group shows mixed outcomes. Pigeon pea is highly vulnerable, with yields declining steeply from −15.09 per cent to −33.51 per cent, reflecting its pronounced sensitivity to temperature rises as observed by Bera (2021). Chickpea also records negative trends, with yield reductions ranging from −4.22 per cent to −11.18 per cent, indicating its susceptibility to heat stress and reduced water availability. However, groundnut presents a strong positive response, with projected yield improvements ranging from 29.17 per cent to 38.04 per cent. This finding aligns with the study by Kadiyala et al. (2021), who reported similar yield gains under a conductive climate effect, highlighting the greater adaptability of groundnut relative to other pulses.
Among oilseeds, rapeseed and mustard display consistent yield reductions, worsening from −11.38 per cent to −18.22 per cent, due to temperature-induced damage to flowering and seed-setting processes. Soyabean shows moderate resilience, with slight yield declines from 5.14 per cent to 1.31 per cent, suggesting a potential plateau under future climates. In contrast, safflower exhibits strong positive projections, with yield gains rising from 9.64 per cent to 36.47 per cent, pointing to its robust adaptability to higher temperatures.
Cash crops exhibit diverse patterns. Cotton demonstrates substantial yield improvements, increasing from 21.58 per cent to 31.78 per cent, consistent with Jans et al. (2021), who projected a 50 per cent rise in cotton production under RCP8.5 by 2,099. The heat tolerance and lower water demands of cotton likely explain this positive trend. In contrast, sugarcane faces yield reductions ranging from −4.37 per cent to −15.48 per cent, reflecting its vulnerability to increased temperature and precipitation variability.
Overall, future projections indicate that climate change will significantly reshape Indian agriculture. Wheat, pigeon pea and sugarcane emerge as among the most adversely affected crops, while groundnut, safflower and cotton show notable yield gains; these trends are further elucidated in Table 1, which highlights the specific thermal thresholds that dictate these varying outcomes. These outcomes underscore the urgent need for crop-specific adaptive strategies, focusing on the promotion of climate-resilient varieties and the development of targeted interventions to support the most vulnerable crops.
Crops with the Highest Tolerance (Turning Point)
4.3 Robustness Check
A comparative analysis between FGLS and PCSE reveals substantial alignment in the direction of effects but divergence in coefficient magnitudes and levels of statistical significance, highlighting PCSE’s relative robustness in estimating climate–crop yield relationships. Both estimators confirm the positive and significant impact of temperature on heat-tolerant crops, with PCSE yielding coefficients of 0.4479 for wheat, 0.9869 for sorghum and 1.3850 for pigeon pea, all significant. FGLS likewise reports significant effects, with somewhat smaller magnitudes: 0.6050, 0.6173 and 1.2514, respectively. For the concave temperature–yield relationship in wheat, both models detect significant negative effects, with PCSE estimating −0.0127 and FGLS being slightly steeper at −0.0159. However, differences emerge in the influence of precipitation; PCSE reports a significant negative effect on rice (−0.0037), capturing waterlogging sensitivity, while FGLS presents a non-significant, near-zero estimate (0.0006). As regards wind speed, PCSE indicates a strongly significant and large negative impact on rapeseed (−1.5024), whereas FGLS also identifies a significant but slightly lower effect (−1.2339). Turning to agricultural inputs like fertiliser, FGLS reports highly significant but inflated impacts, such as 1.1118 for rice and 1.6905 for maize, whereas PCSE provides more measured yet significant estimates: 0.1765 and 1.9184, respectively. Mechanisation proxies like electric pumps and tractors are also significant in PCSE but again appear exaggerated under FGLS. Overall, PCSE not only retains statistical significance across critical variables but does so with coefficient sizes that align more closely with field-level dynamics, making it the more dependable estimator for climate impact studies and policy-oriented agricultural analysis.
5. Conclusion and Policy Implication
This study offers a novel and rigorous assessment of climate change impacts on Indian agriculture by examining crop-specific yield responses across a 50-year panel (1966–2016), capturing a period of significant agricultural transformation driven by policy shifts and technological change. Unlike prior research by Gupta et al. (2014), Burke and Emerick (2016) and Chen and Gong (2021), which focused on aggregated crop groups or shorter time frames, this analysis disaggregates yields across four major crop categories such as cereals, pulses, oilseeds and cash crops, allowing for a more comprehensive and differentiated understanding of climatic and non-climatic influences on productivity. The extended temporal coverage strengthens the reliability of long-run crop–climate interaction estimates, while the broader crop selection moves beyond conventional staples, addressing gaps in existing literature. A key novelty lies in the integration of forward-looking projections using Shared Socioeconomic Pathways (SSP 1–2.6 for precipitation and SSP 5–8.5 for temperature), which reflect both mitigation-aligned and high-impact climate trajectories. These projections are combined with non-linear estimation techniques to identify crop-specific yield inflection points for temperature and precipitation, offering insights into critical thresholds where climate stress may trigger disproportionate yield losses. The adoption of advanced econometric methods further ensures robust estimation by addressing issues of heteroscedasticity, autocorrelation and cross-sectional dependence. Overall, the methodological design and empirical depth of the study signify an original contribution to the literature, generating policy-relevant evidence to inform adaptive strategies tailored to crop-specific climate sensitivities.
The research identifies key vulnerabilities and adaptive capacities across major crop groups in response to climate variability, noting that staple cereals like rice, wheat and millet face substantial risks from temperature fluctuations and shifting precipitation patterns, which threaten India’s food security. Water-intensive crops like rice, central to Indian diets, are especially susceptible to yield reductions in water-stressed areas, whereas drought-tolerant cereals like pearl millet and finger millet demonstrate resilience, making them viable substitutes in drought-prone regions. The study recommends the implementation of improved irrigation methods, such as drip and micro-irrigation, and rainwater harvesting to stabilise cereal yields in vulnerable areas (Di Falco & Chavas, 2006). Pulses, particularly chickpeas and groundnuts, exhibit heightened sensitivity to rising temperatures and rainfall variability, emphasising the need for diversification, wherein heat-tolerant crops like chickpea could be alternatives in warmer regions (Jaworski et al., 2023). Soil conservation practices like mulching, no-till farming and bio-fertiliser integration are recommended to bolster the resilience of pulse crops by enhancing soil moisture retention and quality.
As regards oilseeds, soyabean and sunflower emerge as relatively resilient under moderate climate fluctuations, making them suitable for regions experiencing increased variability. The study suggests enhancing irrigation efficiency through the adoption of sub-surface methods and promoting rainwater harvesting infrastructure to support yield stability for these crops even during dry spells. Cash crops like sugarcane and cotton, with high water demands, face intensified risks under climate stress; thus, introducing drought-tolerant varieties and practices like crop rotation and intercropping can help mitigate these vulnerabilities, while enhancing long-term productivity (Da Cruz & Machado, 2023; Okemo et al., 2024). To mitigate food security risks and advance agricultural sustainability, policies should focus on adopting climate-smart agricultural technologies, including advanced breeding, improved irrigation and soil management techniques to increase crop resilience against extreme weather. By diversifying crop production and emphasising sustainable water use, India can reinforce agricultural stability amidst climate change, safeguarding both rural livelihoods and national food security.
Despite its strengths, this study has certain limitations. The CMIE and IndiaStat databases, while comprehensive, have limitations in crop-specific data for variables such as agricultural labour, mechanisation and fertiliser usage, which required us to rely on weighted estimates based on cropping area. While this approach provides a broad perspective, it may not fully capture resource allocation nuances across individual crops, introducing potential inaccuracies. Additionally, data on agricultural equipment use were available only at five-year intervals, necessitating extrapolation that may impact the analysis’s precision. Furthermore, we did not incorporate climate variables like soil moisture and radiation, given their limited availability and inconsistent measurements across states. While these variables are undoubtedly important, they were omitted to ensure robustness in longitudinal data.
Future research should take the next steps towards developing refined adaptation strategies, particularly by testing additional climate variables such as soil moisture, solar radiation and extreme weather events, which were limited in the current study due to data constraints. Comparative analyses of adaptation strategies between India and other developing nations with similar agro-ecological climates would also help identify globally applicable solutions. An interdisciplinary approach combining insights from agricultural science, environmental economics and climate modelling could improve our understanding of effective adaptation techniques. Such an approach, like evaluating the significance of crop genetic variety and innovative irrigation systems, could provide useful insights on enhancing resilience across various agro-ecological zones.
Footnotes
Acknowledgements
We thank the handling editor and anonymous reviewers for their constructive comments and suggestions that significantly improved the manuscript. An earlier version of this paper was presented at the 8th SANEM Annual Economists’ Conference (21–23 February 2025, Dhaka). We are grateful to the participants at these events for their valuable feedback. Any remaining errors are our own.
Data Availability Statement
The data that support the findings of this study are available upon reasonable request from the corresponding author.
Declaration of Conflicting Interests
The authors declared no potential conflicts of interest with respect to the research, authorship and/or publication of this article.
Funding
The authors received no financial support for the research, authorship and/or publication of this article.
Permission to Reproduce Material from Other Sources
All materials included in this manuscript are original or used with appropriate permission and citation.
Appendix A: Notes
A1. For two models, say one is linear and the other in semi-log form, then Sagan’s criteria would prompt towards the best-fit model using the following expression: ∀ = [φ1/g2 φ2] n , where φ1 is the standard deviation of the combined residual (sum of random or fixed error component and overall stochastic error component) of Model 1; φ2 is the standard deviation of the combined residual of Model 2; g2 is the geometric mean of the dependent variable of Model 2, and n is the number of observations. If ∀ < 1, then Model 1 is favoured, and if ∀ > 1, then Model 2 is favoured.
A2. Assuming a linear model, where the focus is on the k × 1 vector β, which represents both climate and non-climate coefficients, if Cov (x, c) ≠ 0, then including ‘γ’ as the error term is inappropriate, and β cannot be estimated without additional information. The model with the error term is fit.
The OLS βˆ estimator becomes biased and inconsistent because f (x, c) ≠ 0, as ‘γ’ is correlated with ‘x’, making pooled OLS inappropriate. Before applying the error correction panel model, it is important to discuss the nature of the unobserved effect so that the model is
A3. The null hypothesis of the Hausman test assumes random effect estimates are consistent and efficient. The Hausman test shows that ‘γ’ has a random effect across states, so a pooled OLS estimation of the unobserved effect model is used as
For pooled OLS to be consistent, the estimation of this equation requires as follows:
The expression for the test is given as
where
The initial model with yit violates certain assumptions because yi, t – 1 and γi are correlated, while ϑit is serially correlated with γi across time periods (see Equation (A3)). This correlation highlights the need for more robust tests and a random effects analysis, as the model’s assumptions are no longer met. Specifically, we must impose stricter assumptions to ensure the robustness of the estimates, including the conditions of strict exogeneity and orthogonality between ϑit and xit.
The assumptions required are:
The random effects approach accounts for serial correlation in the composite error term, where ϑit = γi + Uit thus capturing spatial correlation as well. The suitable regression model for this set-up is expressed as
Here, yit is the observed crop yield for state ith in time period tth, while xit is a vector of K × 1 observation on non-stochastic regressors and Uit is the regression disturbance, defined as
In this case, μ = (μ1,…, μN) represents the random state effects vector, assumed to be independently and identically distributed (IIN) with a variance of μ2. The autoregressive coefficient t measures the degree of correlation, while W is an N × N weight matrix. Equation (A6) indicates that part of the current error depends on past errors, with ϑit representing the new disturbance or shock for the current period. This autoregressive nature of the error term is crucial in the model’s design. The error term εit is thus represented as
Hence, the full disturbance term U is expressed as
Now, under the assumption of normality, the likelihood function is
A4. Marginal effect of temperature:
Marginal effect of precipitation:
Here, y denotes crop yield (tonne/ha), φ s are the regression coefficients associated with the weather variable, its quadratic terms, and φ7 is the regression coefficients of the intersection terms between temperature and precipitation, respectively. Using y = 100 as the baseline, we compute the effects of marginal increases in weather parameters.
The total effect on crop yield is the sum of the marginal effects of temperature and precipitation. Additionally, we calculate the elasticity of crop yield in response to weather variables using the following equation:
Appendix B: Tables
Panel Regression (FGLS) Estimates for Oilseeds and Other Seeds and Cash Crops
| Variables | Oilseeds and Other Seeds | Cash Crops | ||||||||
| Sesamum | Rapeseed and Mustard | Castor | Linseed | Sunflower | Safflower | Soyabean | Oilseeds | Sugarcane | Cotton | |
| Average temperature | 0.6192*** | 0.7804*** | 1.1476*** | 0.1201 | 0.1289 | −1.7207*** | 0.2967 | −0.0313 | 0.7467*** | −0.6357*** |
| (−0.1092) | (−0.1284) | (−0.2778) | (−0.1880) | (−0.4611) | (−0.5160) | (−0.2949) | (−0.0696) | (−0.1088) | (−0.1421) | |
| Temperature square | −0.0114*** | −0.0197*** | −0.0182*** | −0.0008 | −0.0043 | 0.0327*** | −0.0056 | 0.0013 | −0.01469*** | 0.0164*** |
| (−0.0027) | (−0.0032) | (−0.0068) | (−0.0047) | (−0.0102) | (−0.0127) | (−0.0072) | (−0.0017) | (−0.0027) | (−0.0036) | |
| Temperature deviation | 0.1274* | 0.0108 | −0.0809 | 0.0931 | −0.2897 | −0.5962* | −0.0935 | −0.0182 | 0.0550 | 0.0726 |
| (−0.0731) | (−0.0872) | (−0.1783) | (−0.1278) | (−0.2237) | (−0.3289) | (−0.2388) | (−0.0441) | (−0.0774) | (−0.0943) | |
| Wind speed | −0.7964** | −1.2339*** | −1.2690* | −1.0351* | −4.6779*** | −0.6558 | 3.6612*** | 0.0731 | 0.3541 | 0.2432 |
| (−0.3230) | (−0.3532) | (−0.6913) | (−0.6019) | (−1.1410) | (−1.3035) | (−1.3222) | (−0.1830) | (−0.3353) | (−0.3380) | |
| Wind speed square | 0.1531* | −0.0127 | 0.3595** | 0.1020 | 1.3277*** | 0.5786 | −1.0161*** | −0.0187 | −0.0193 | −0.0727 |
| (−0.0814) | (−0.0936) | (−0.1768) | (−0.1682) | (−0.3177) | (−0.3555) | (−0.3571) | (−0.0480) | (−0.0855) | (−0.0893) | |
| Wind speed deviation | 0.0051 | 1.3084*** | −0.4395 | 0.4786* | −0.0211 | −1.7806*** | 0.3853 | 0.0859 | −0.1272 | 0.3074* |
| (−0.1487) | (−0.1671) | (*−0.3176) | (−0.2715) | (−0.4251) | (−0.6479) | (−0.4827) | (−0.0873) | (−0.1563) | (−0.1788) | |
| Evapotranspiration | 0.0012 | −0.0088*** | −0.0023 | 0.0011 | 0.0088** | −0.0332*** | −0.0012 | 0.0017** | 0.0057*** | 0.0014 |
| (−0.0010) | (−0.0016) | −(0.0029) | (−0.0024) | (−0.0043) | (−0.0105) | (−0.0053) | (−0.0008) | (−0.0011) | (−0.0013) | |
| Evapotranspiration square | −4.80E−07 | 4.41E−6*** | 1.35E−06 | −1.85E−06 | −4.99E−5* | 2.033E−5*** | −6.00E−07 | −9.4E−7* | −3.97E−7*** | −1.52E−7* |
| (−6.5E−07) | (−1.0E−06) | −(1.9E−06) | (−1.6E−06) | (−2.7E−06) | (−6.3E−06) | (−3.4E−06) | (−5.4E−07) | (−7.1E−07) | (−8.0E−07) | |
| Evapotranspiration deviation | −0.0009* | −0.0015** | 0.0004 | 0.0000 | 0.0013 | −0.0010 | −0.0014 | 0.0001 | 0.0012** | −0.0006 |
| (−0.0005) | (−0.0006) | −(0.0012) | (−0.0010) | (−0.0015) | 9−0.0025) | (−0.0018) | (−0.0003) | (−0.0005) | (−0.0007) | |
| Precipitation | 0.0004* | 0.0012*** | 0.0015** | 0.0012** | −0.0006 | −0.0008 | 0.0009 | −0.0004** | 0.0000 | −0.0003 |
| (−0.0002) | (−0.00030 | −(0.0007) | (−0.0005) | (−0.0013) | (−0.0013) | (−0.0016) | (−0.0002) | (−0.0002) | (−0.0003) | |
| Precipitation square | −8.00E−08 | −3.5E−7*** | −3.7E−7* | −4.0E−7*** | 1.20E−07 | 8.0E−7** | −1.20E−07 | 9.0E−8* | -1.3E−7** | 2.00E−08 |
| (−5.0E−08) | (−1.1E−07) | −(2.1E−07) | (−1.6E−07) | (−4.6E−07) | (−3.9E−07) | (−6.0E−07) | (−6.0E−08) | (−6.0E−08) | (−7.0E−08) | |
| Precipitation deviation | −1.43E−06 | 1.57E−05 | −0.0004 | −0.0002 | −0.0004 | −0.0001 | 0.0001 | −0.0001 | −0.0001 | −0.0001 |
| (−0.00010 | (−0.0002) | −(0.0004) | (−0.0003) | (−0.0005) | (−0.0006) | (−0.0005) | (−0.0001) | (−0.0002) | (−0.0002) | |
| Rural literacy | 8.24E−6*** | 1.504E−5*** | 9.70E−07 | 4.03E−06 | 0.0005*** | 0.0009*** | 0.0001*** | 2.157E−5*** | 1.319E−6*** | −1.038E−6** |
| (−3.1E−06) | (−4.2E−06) | −(7.3E−06) | (−4.5E−06) | (−7.3E−06) | (−1.4E−05) | (−7.9E−06) | (−2.3E−06) | (−4.1E−06) | (−4.4E−06) | |
| Agricultural labour | −0.4534** | 0.0520** | 0.6092*** | −0.0663** | 0.0042 | 0.5372** | 0.4041*** | 0.0903*** | 0.0268 | 0.3809*** |
| (−0.1918) | (−0.0251) | −(0.0510) | (−0.0320) | (−0.0866) | (−0.2705) | (−0.0507) | (−0.0138) | (−0.0225) | (−0.0745) | |
| Irrigation | 0.0208 | −0.3555*** | −1.1455*** | 0.3321 | 0.3782 | −1.0894* | −1.2499*** | 0.0326 | 0.3322*** | −0.2553* |
| (−0.0468) | (−0.1377) | −(0.2510) | (−0.2066) | (−0.2954) | (−0.6524) | (−0.2617) | (−0.0718) | (−0.1139) | (−0.1519) | |
| Fertiliser | 10.1453 | 0.4301 | 1.3072 | −0.0236 | 2.3845** | 2.1062* | 0.6381 | 0.3745* | −0.2751** | −0.2009 |
| (−10.3604) | (−0.4029) | −(0.8284) | (−0.2823) | (−1.0021) | (−1.2383) | (−0.7599) | (−0.2247) | (−0.1198) | (−0.1392) | |
| Electric pump set | −0.0052 | −2.9032*** | −4.3498*** | −4.4655*** | −1.4920 | 1.7349 | 0.7857 | −0.5896 | −0.1375 | 0.2931 |
| (−0.032) | (−0.770) | −(1.183) | (−1.089) | (−1.167) | (−2.112) | (−1.223) | (−0.364) | (−0.379) | (−0.457) | |
| Tractor | −0.3124 | 10.5697*** | 6.3944 | −0.0519 | 25.0391*** | 84.0748*** | 54.5729*** | 6.9228*** | −12.0456*** | 4.8570 |
| (−0.2221) | (−2.6675) | −(8.5425) | (−0.4496) | (−5.7797) | (−23.0214) | (−12.0269) | (−1.6753) | (−2.5177) | (−3.4133) | |
