Comparing mechanistic and empirical model projections of crop suitability and productivity: implications for ecological forecasting

• Published in: Global ecology and biogeography 2013-08, Vol.22 (8), p.1007-1018
• Main Authors: Estes, L. D., Bradley, B. A., Beukes, H., Hole, D. G., Lau, M., Oppenheimer, M. G., Schulze, R., Tadross, M. A., Turner, W. R.
• Format: Article
• Language: English
• Subjects: Abundance; Accuracy; Animal and plant ecology; Animal, plant and microbial ecology; Biological and medical sciences; Climate change; Corn; Crop ecology; Crops; Ecological modeling; ecological niche model; Ecosystem models; empirical model; Empirical modeling; Fundamental and applied biological sciences. Psychology; General aspects; habitat suitability model; intercomparison; maize; mechanistic model; Modeling; Productivity; South Africa; Species; species distribution model; Studies; Synecology; Zea mays
• ISSN: 1466-822X; 1466-8238; 1466-822X
• DOI: 10.1111/geb.12034

Aim: Intercomparison of mechanistic and empirical models is an important step towards improving projections of potential species distribution and abundance. We aim to compare suitability and productivity estimates for a well-understood crop species to evaluate the strengths and weaknesses of mechanistic versus empirical modelling. Location: South Africa. Methods: We compared four habitat suitability models for dryland maize based on climate and soil predictors. Two were created using maximum entropy (MAXENT), the first based on national crop distribution points and the second based only on locations with high productivity. The third approach used a generalized additive model (GAM) trained with continuous productivity data derived from the satellite normalized difference vegetation index (NDVI). The fourth model was a mechanistic crop growth model (DSSAT) made spatially explicit. We tested model accuracy by comparing the results with observed productivity derived from MODIS NDVI and with observed suitability based on the current spatial distribution of maize crop fields. Results: The GAM and DSSAT results were linearly correlated to NDVI-measured yield (R² = 0.75 and 0.37, respectively). MAXENT suitability values were not linearly related to yield (R² = 0.08); however, a MAXENT model based on occurrences of high-productivity maize was linearly related to yield (R² = 0.62). All models produced crop suitability maps of similarly good accuracy (Kappa = 0.73-75). Main conclusions: These findings suggest that empirical models can achieve the same or better accuracy as mechanistic models for predicting both suitability (i.e. species range) and productivity (i.e. species abundance). While MAXENT could not predict productivity across the species range when trained on all occurrences, it could when trained with a high-productivity subset, suggesting that ecological niche models can be adjusted to better correlate with species abundance.