Skip to main content
Have a personal or library account? Click to login
Harvested Area, Climate Variability, and Cereal Production Dynamics in Algeria: Var Evidence, 1990–2024 Cover

Harvested Area, Climate Variability, and Cereal Production Dynamics in Algeria: Var Evidence, 1990–2024

Open Access
|Sep 2026

Figures & Tables

Table 1.

Variables and data sources

SymbolVariable and definitionTransformationUnitSource
LPRODFAOSTAT “Cereals, primary” – Production; annual national cereal production.Natural logarithmMetric tonnesFAOSTAT
LHAFAOSTAT “Cereals, primary” – Area harvested; annual area that reaches harvest.Natural logarithmHectaresFAOSTAT
LPRECCCKP national (ADM0) precipitation indicator (pr); annual accumulated precipitation.Natural logarithmMillimetresWorld Bank CCKP
LTEMPCCKP national (ADM0) average surface-air temperature indicator (tas).Natural logarithm; level-form robustness checkDegrees CelsiusWorld Bank CCKP
D1Low-precipitation anomaly indicator: 1 if PREC < 67.72 mm; 0 otherwise.Deterministic binary controlBinaryAuthors’ calculation
D2High-precipitation anomaly indicator: 1 if PREC > 95.83 mm; 0 otherwise.Deterministic binary controlBinaryAuthors’ calculation

[i] Source: own elaboration based on FAOSTAT and World Bank Climate Change Knowledge Portal data.

Table 2.

Descriptive statistics

StatisticLPRODLHALPRECLTEMP
Mean14.932214.73354.38883.1617
Median15.076114.80364.43433.1633
Maximum15.618215.11404.67743.1855
Minimum13.676213.87204.04863.1280
Std. Dev.0.51740.31850.17930.0125
Skewness–1.0809–1.2807–0.4489–0.8860
Kurtosis3.25174.01771.92464.1667
Jarque-Bera6.908911.07882.86186.5651
Probability (Jarque-Bera)0.03160.00390.23900.0375
Observations35353535

[i] Note: Statistics refer to the logged variables over 1990–2024.

[ii] Source: own elaboration.

Fig. 1.

Time-series plots of the log-transformed variables, 1990–2024

Source: own elaboration based on FAOSTAT and World Bank Climate Change Knowledge Portal data.

Table 3.

Unit-root and stationarity tests

Panel A. Conventional ADF and PP unit-root tests
VariableADF: CADF: C+TADF: nonePP: CPP: C+TPP: noneDecision
LPROD–4.4949***–5.7846***0.1515–4.6596***–5.7846***0.5140I(0)
LHA–5.7619***–3.9701**0.4212–5.9625***–6.3377***0.0718I(0)
LPREC–3.3528**–3.8020**–0.6048–3.2801**–3.7718**–1.4537I(0)
LTEMP–2.1619–4.1701**0.4554–1.9491–4.2693***1.4484Stationary under trend specification
Panel B. Zivot-Andrews unit-root test with one endogenous structural break
VariableZA statisticSelected lagBreak year5% critical valueDecision
LPROD–7.106402002–5.0733Reject unit root
LHA–4.692732016–5.0733Do not reject
LPREC–5.667602002–5.0733Reject unit root
LTEMP–5.878101996–5.0733Reject unit root

[i] Panel A: C = constant; C+T = constant and deterministic trend. The null hypothesis is a unit root. LTEMP rejects the unit-root null only when the trend is included.

[ii] Panel B: Model C allows one endogenous break in intercept and trend; maximum lag = 4; lag selected by AIC. Critical values: 1% = −5.5756, 5% = −5.0733, 10% = −4.8267.

[iii] Source: own elaboration.

Fig. 2.

VAR Stability Diagnostic Based on Inverse AR Roots

Source: own elaboration based on EViews stability diagnostics.

Table 4.

VAR lag order selection

LagLogLLRFPEAICSCHQ
0142.6939NA3.34e-09–8.168369–7.618718–7.986175
1173.445048.04858*1.36e-09*–9.090312*–7.807793*–8.665193*
2181.396310.436092.46e-09–8.587268–6.571881–7.919225
3196.385415.925913.21e-09–8.524087–5.775832–7.613118

* indicates the lag selected by the criterion.

Source: own elaboration.

Table 5.

Quandt-Andrews unknown-breakpoint tests for VAR parameter stability

VAR equationMaximum LR FCandidate yearp-valueDecision at 5%
LPROD2.478720010.3049Do not reject stability
LHA4.150919970.0192Reject stability
LPREC2.100020030.4865Do not reject stability
LTEMP3.188519970.1052Do not reject stability

[i] Note: The test searches the 15% trimmed sample (1997–2019). Candidate years correspond to the maximum LR F-statistic; a candidate year is not treated as a confirmed break when the stability null is not rejected.

[ii] Source: own elaboration.

Table 6.

VAR residual diagnostic tests

DiagnosticStatisticdfp-valueAdditional statisticAdditional pDecision
LM, lag 1LRE = 7.6842160.9577Rao F = 0.45640.9584No serial correlation
LM, lag 2LRE = 16.8861160.3930Rao F = 1.07450.3978No serial correlation
LM, lags 1–2LRE = 32.4358320.4453Rao F = 1.01470.4686No serial correlation
Joint skewness–40.3859––Normality not rejected
Joint kurtosis–40.9030––Normality not rejected
Joint Jarque-BeraJB = 5.195780.7365––Normality not rejected
White heteroskedasticity–1000.1359––Homoskedasticity not rejected

[i] Note: The respective null hypotheses are no residual serial correlation, multivariate normal residuals and homoskedasticity.

[ii] Source: own elaboration.

Table 7.

VAR Granger causality/block-exogeneity Wald tests for LPROD

Excluded variabledfChi-squarep-value
LHA11.35870.2438
LPREC10.11660.7327
LTEMP12.98040.0843
All three38.64130.0345

[i] Note: Included observations = 34 after lag adjustment. The test is interpreted as predictive content, not structural causality.

[ii] Source: own elaboration.

Fig. 3.

Bootstrap Impulse Response Functions of LPROD

Note: Shaded bands are 95% standard percentile-bootstrap confidence intervals based on 2,000 EViews repetitions. Benchmark Cholesky ordering: LPROD, LHA, LPREC, LTEMP.

Table 8.

Variance decomposition of total cereal production (LPROD)

PeriodS.E.LPROD (%)LHA (%)LPREC (%)LTEMP (%)
10.4647100.000.000.000.00
20.482793.072.920.123.87
40.507884.358.801.295.55
60.519480.8211.221.696.24
80.524979.2212.321.886.56
100.527678.4712.851.976.71

[i] Values are descriptive point estimates of LPROD forecast-error variance shares under the benchmark Cholesky ordering. They are conditional on the model, horizon and identification scheme; uncertainty intervals are not reported, so no inferential significance is attached to differences among the shares.

[ii] Source: own elaboration.

Table 9.

Summary of robustness checks

Robustness checkBenchmarkRobustness resultInterpretation
Trend-augmented VARJoint Granger p = 0.0345With trend: joint p = 0.6883; LHA = 0.9278, LPREC = 0.2841, LTEMP = 0.5202Benchmark joint predictive result is trend-sensitive
Alternative Cholesky orderingH10 FEVD: 78.47, 12.85, 1.97, 6.71%Alternative: LPROD 78.47%, LHA 14.48%, LPREC 0.17%, LTEMP 6.87%Shares depend on recursive identification; none is interpreted causally
Temperature in levelsJoint p with LTEMP = 0.0345With TEMP: joint p = 0.0341; LHA = 0.2424, LPREC = 0.7365, TEMP = 0.0832Log transformation does not drive the benchmark result; trend sensitivity remains stronger
Quandt-Andrews testsBenchmark assumes parameter stabilityStability rejected only for LHA (p = 0.0192); not rejected for LPROD, LPREC or LTEMPInstability is equation-specific; no common break is imposed
Bootstrap IRFsConventional asymptotics are weak in N = 352,000-repetition 95% percentile-bootstrap bands; LTEMP horizon-2 lower bound effectively zeroClimate IRF claims are weakened by bootstrap uncertainty
Without D1 and D2Joint p = 0.0345; H10 LPREC FEVD = 1.97%Joint p ≈ 0.0506; LTEMP p ≈ 0.0357; H10 FEVD: LPROD 75.14%, LHA 10.06%, LPREC 7.40%, LTEMP 7.41%D1/D2 absorb part of LPREC information; precipitation results are specification-sensitive

[i] H10 = horizon 10. FEVD figures are point estimates. All robustness results are interpreted as model-conditional evidence.

[ii] Source: own elaboration.

DOI: https://doi.org/10.17306/j.jard.2026.3.00022r1 | Journal eISSN: 1899-5772 | Journal ISSN: 1899-5241
Language: English
Page range: 348 - 360
Accepted on: Sep 7, 2026
Published on: Sep 30, 2026
Published by: Poznań University of Life Science
In partnership with: Paradigm Publishing Services
Publication frequency: 4 issues per year

© 2026 Siham Boudab, Nesrine Arous, Sami Bendjeddou, Hamza Ramli, published by Poznań University of Life Science
This work is licensed under the Creative Commons Attribution 4.0 License.