Role of Modelling in International Crop Research : Overview and Some Case Studies

Crop modelling has the potential to contribute to global food and nutrition security. This paper briefly examines the history of crop modelling by international crop research centres of the CGIAR (formerly Consultative Group on International Agricultural Research but now known simply as CGIAR), whose primary focus is on less developed countries. Basic principles of crop modelling building up to a Genotype × Environment × Management × Socioeconomic (G × E × M × S) paradigm, are explained. Modelling has contributed to better understanding of crop performance and yield gaps, better prediction of pest and insect outbreaks, and improving the efficiency of crop management including irrigation systems and optimization of planting dates. New developments include, for example, use of remote sensed data and mobile phone technology linked to crop management decision support models, data sharing in the new era of big data, and the use of genomic selection and crop simulation models linked to environmental data to help make crop breeding decisions. Socio-economic applications include foresight analysis of agricultural systems under global change scenarios, and the consequences of potential food system shocks are also described. These approaches are discussed in this paper which also calls for closer collaboration among disciplines in order to better serve the crop research and development communities by providing model based recommendations ranging from policy development at the level of governmental agencies to direct crop management support for resource poor farmers.


Introduction
Current rates of yield growth and overall production are not sufficient to satisfy future demand [1] and climate change is predicted to add at least another 60 million hungry people by 2050 [2].Rosegrant et al. [3] predicted that-in the absence of unprecedented, coordinated measures to raise productivity-consumers will pay double for staple foods by 2050 in real terms.While crop-focused research still achieves impressive productivity gains [4], multiple approaches must be employed to improve leverage of knowledge, expertise, and data to achieve greater returns on investment in agriculture and maintain global food security under the added challenges of climate change, resource scarcity, growing populations, and changing demand.
Crop modelling has the potential to significantly contribute to global food and nutrition security.New technologies and conceptual breakthroughs have contributed to a better understanding of crop performance and yield gaps, genetic gains, better prediction of pest and insect outbreaks, more efficient irrigation systems, and optimized planting dates.Newer developments include the use of remote-sensed data and mobile phone technology linked to simulation models [5]; improved geographic information system (GIS) techniques [6]; strategic crossing and selection models making use of high-throughput genotyping and phenotyping to increase genetic gains [7]; as well as data sharing and standardization in the new era of big data [8].Renewed interest in crop modelling that addresses the challenges of climate change has already led to collaborative initiatives such as the Agricultural Model Intercomparison and Improvement Project (AgMIP) [9][10][11], the CGIAR Platform for Big Data in Agriculture (launched in 2017) [8], as well as other initiatives for dataset standardization [12].Additional investments and large-scale collaborations are needed to improve the collection of open access, easy to use data for crop modelling purposes [13].
This paper briefly outlines the history of crop modelling, including major theoretical advances and their practical applications by international crop research centres of the CGIAR.It shows that theory and field-level implementation evolved in tandem, with models and their applications making mutual improvements and refinements over time.Several successful applications of models are presented, highlighting their contributions in food security at multiple levels.Basic principles of crop modelling, building up to a new Genotype × Environment × Management × Socioeconomic (G × E × M × S) paradigm, are explained.The application of these principles to environmental modelling, as well as crop management and breeding, is explained using multiple examples from CGIAR Centre's work throughout developing regions.Socioeconomic applications include foresight analysis of agricultural systems under global change scenarios, and the consequences of potential food system shocks.Using crop modelling to achieve strong global development impacts depends on the contribution of a wide range of stakeholders-from academia to private companies-to large, multi-location datasets.The need for greater collaboration between modelling disciplines is a driver of this paper.
Crop modelling in the target environments of the CGIAR centres has its own unique challenges in addition to the fact that many of the regions face direct climate change and other resource limitations, e.g., related soil erosion and unavailability of good quality water.For example, crop models require comprehensive weather data for simulations to run and to explain G × E data, which is often not available or insufficient.Other aspects of environmental characterization of target environments may also be lacking, including data on soils types, fertility regimes and even water application.While the expansion of crop modelling approaches is still somewhat limited by the availability of quality data for calibration, validation, and evaluation, recent progress in the development of large-scale agro-ecological datasets, facilitated by global initiatives such as AgMIP [10] and the Big Data Platform in Agriculture [8], are helping fill the data gap.
Modelling innovations can address increasing concerns on sustainable food production, nutrition, and natural resource management challenges, and their consequences on closely associated socioeconomic issues such as conflict, migration, human health, and gender inequality.While this Agronomy 2018, 8, 291 3 of 46 review takes stock of the state-of-the-art of crop related modelling within CGIAR, it also highlights where improvements are desirable and feasible.
This review considers crop productivity models within the framework of the widely used conceptual model G × E × M (GEM), representing the three most important biophysical variables that directly determine crop growth, and their interactions.Namely, expression of the genetic potential of a cultivar (G), the physical environment of the crop (E), and crop management inputs (M).However, these three elements are also influenced indirectly by socioeconomic factors (S), such as supply and demand of inputs/outputs, finance and credit, agricultural policies, and the adaptive practices of farmers and consumers, thus extending the paradigm to G × E × M × S.

Environment
In crop productivity models, E must be defined in terms of the variables that most strongly affect growth.This presents some challenges: firstly, E has almost infinite dimensions.The crop inhabits the soil and the atmosphere, where it interacts with physical, chemical, and biotic factors.Light and temperature regimes are in constant flux and essential molecules exist in solid, liquid, and gas phases-often interchangeably.When defining E, macroscopic, microscopic, and biotic domains must also be considered.Which variables most strongly affect growth is still a subject for research.Secondly, data of many of these factors are not readily available-hardly surprising considering that crops are grown on over 1 billion ha of land worldwide.Thirdly, while some of these dimensions are predictable, the key drivers of productivity-temperature, radiation, and rainfall-are highly variable on most relevant time scales and subject to climate change.Long term detailed weather datasets are crucial, preferably at a high density.It is for these reasons that GIS (Geographic Information System) has emerged as one of the most important and necessary disciplines in crop research and modelling in recent decades [6].

Management
M is essentially a subset of E. Yet it is given special status since, unlike E, the factors associated with M are controlled.The most widely adopted crop models are those that demonstrate how yield gaps may be closed using M interventions, namely fertilizer, irrigation, planting methods/dates, tillage, and other crop husbandry approaches [214].The most advanced models evaluate a whole package of M interventions within a farmer-friendly decision support tool [215].Taking a long-term view, such models can also be used to assess the broader environmental consequences of M.

Genotype
G is the most challenging aspect of crop modelling, partly because of the complexity and diversity of biological organisms, and partly because-unlike livestock, for example-plants do not systematically control their temperature, making their interaction with E a constantly moving target.While models of G may be robust for 'model species' grown in controlled environments, and conceptual models of G may provide hypotheses about candidate traits for future target environments (e.g., [216]), the greatest challenge for crop modelling is to simulate G × E interaction to refine crop breeding strategies and deploy new cultivars.Crop models are typically driven by the effect of E on growth, especially water, nitrogen (N), radiation and temperature.However, genotypic components are less accurately simulated, except for development rate, i.e., the duration to flowering and grain filling period (see [217,218] and references therein).

Socioeconomics
While S is potentially quantifiable and can inform the parameterization of crop models and their scenario developments, this is only possible when detailed household-level survey data is available.For example, World Bank's Living Standards Measurement Study-Integrated Surveys on Agriculture provides a nationally-representative, comprehensive snapshot of farmers' livelihoods and plot-level farm management practices, yet its coverage is currently limited to eight countries in sub-Saharan Africa.When available, this type of micro-level data provides context for on-farm productivity and the off-farm socioeconomic environments that impact farming households, such as their proximity to markets, access to extension and financial services, and capacity to manage production risks.This information can help define the crop management baseline, parameterize model input data, calibrate the model, evaluate model performance, and develop practical what-if scenarios that are affordable for farmers.Socioeconomic inputs (S) are discussed in Section 6, in the context of calculating food systems risks.At the more aggregate level, market interactions combined with policies, institutional arrangements, and governance define the socioeconomic circumstances in which crop management decisions (both choice of genotypes and management practices) are made at the farm level.The combination of the bio-physical environment and socioeconomic circumstances define the development domains where specific technologies can or cannot prevail [219][220][221].

Interactions
The full G × E × M × S paradigm can be overly complex, and most practical applications of crop modelling to date relate to M recommendations [215].Thus, it is practical to consider two main classes of interaction: those where G is considered as a fixed/generic effect (i.e., E × M and M × S interactions), or where G is a variable itself (G × E, G × M and G × E × M) (Figure 1, Sections 4-6).
of interaction: those where G is considered as a fixed/generic effect (i.e., E × M and M × S interactions), or where G is a variable itself (G × E, G × M and G × E × M) (Figure 1, Sections 4-6).
The modelling and research approaches outlined above are developed further in dedicated sections below.However, since models are only as good as their inputs, this review also addresses data standardization and sharing, as well as establishing and using phenotyping networksactivities that will hugely strengthen the scope and reliability of crop modelling.

The Physical Environment
Environmental data are a crucial input for crop, environmental, and systems modelling.Their availability in terms of time lines, recording frequency, and location limits the ability and the representativeness of models.Ideally, sub-daily or daily data for at least a minimum set of environmental variables, particularly weather data (Table 2), are required [23,24,222,223].Due to the increasingly global nature of cropping systems, the density of automatic weather stations is not sufficient, public access to data is limited, and/or the parameters needed for crop models are not sufficient (especially for larger areas in developing countries).In such cases, interpolated outputs from large climate datasets such as CRU (Climate Research Unit) or Worldclim are commonly used to provide monthly long-term average or long-term timeline environmental data for the models [224][225][226][227].However, datasets like CRU have low resolution, which may cause problems with loss of spatial detail, especially in areas with high spatial heterogeneity [228].
Weather generators integrated into crop models or standalone models such as WGEN, SIMMETEO, or MARKSIM are commonly used to generate daily time series data [229][230][231][232].The same applies to climate change applications to incorporate future climate scenarios into crop models, where all data utilized must be downscaled using dynamic or statistical approaches from larger global circulation model outputs and then run through a weather generator [230,[233][234][235][236][237].The modelling and research approaches outlined above are developed further in dedicated sections below.However, since models are only as good as their inputs, this review also addresses data standardization and sharing, as well as establishing and using phenotyping networks-activities that will hugely strengthen the scope and reliability of crop modelling.

The Physical Environment
Environmental data are a crucial input for crop, environmental, and systems modelling.Their availability in terms of time lines, recording frequency, and location limits the ability and the representativeness of models.Ideally, sub-daily or daily data for at least a minimum set of environmental variables, particularly weather data (Table 2), are required [23,24,222,223].Due to the increasingly global nature of cropping systems, the density of automatic weather stations is not sufficient, public access to data is limited, and/or the parameters needed for crop models are not sufficient (especially for larger areas in developing countries).In such cases, interpolated outputs from large climate datasets such as CRU (Climate Research Unit) or Worldclim are commonly used to provide monthly long-term average or long-term timeline environmental data for the models [224][225][226][227].However, datasets like CRU have low resolution, which may cause problems with loss of spatial detail, especially in areas with high spatial heterogeneity [228].
Weather generators integrated into crop models or standalone models such as WGEN, SIMMETEO, or MARKSIM are commonly used to generate daily time series data [229][230][231][232].The same applies to climate change applications to incorporate future climate scenarios into crop models, where all data utilized must be downscaled using dynamic or statistical approaches from larger global circulation model outputs and then run through a weather generator [230,[233][234][235][236][237].
Soil data are equally important, and similar problems to those outlined above limit their use for larger areas.Specific sites (e.g., research stations) often have detailed soil profile data available, which can be used to create high quality soil input files, but the outputs of models for larger areas (country, regional, or global level) are limited in their reliability.Recent efforts have attempted to add detail to global soil databases and make them available for the modelling community [238][239][240][241].
Satellite-derived datasets, such as NASA's Tropical Rainfall Monitoring Mission (TRMM) and the Global Precipitation Mission (GPM), indicate that it might soon be possible to use a combination of remotely-sensed and ground-based data to collect precise, near real-time environmental data to run models at any given location-even in environments where data are scarce [242][243][244][245]. Remote sensing with a variety of sensors already allows monitoring and acquisition of surface soil moisture and solar radiation [246][247][248], as well as soil chemical and physical properties [249] (see Section 4.4).
Table 2. Generic standards for field phenotyping platforms to support simulation modelling (adapted from [250]). Treatments:

•
Managed stress environments and precision treatments with adequate irrigation options for the range of target environments to be adequately simulated.
Another aspect of modelling E includes that of predicting climate change, for which a vast body of literature exists.Given the timeframe for crop genetic improvement, the outputs of these models can be used by crop breeders and policymakers to plan long-term strategic investments.Pest and disease pressure is another key factor in modelling crop response to E, which is covered in Section 3.2.

Tools for Monitoring and Managing Pests and Diseases
Pests and diseases are an ever-present threat to crops worldwide, since they continually evolve to new pathotypes.As climate change and weather fluxes make it harder to predict their spread and virulence, multi-variable predictive models are becoming even more valuable as tools for crop protection.For pests, models can link organisms' physiological and phenological behaviour to environmental factors.For disease modelling, the life cycle is often represented by an infection chain, which corresponds to functional traits that can be converted to quantifiable processes and their responses linked to biological and environmental factors [252].Within the CGIAR, these tools are being used to answer several relevant questions [121], such as: (1) how can we predict pest/disease vector population and disease transmission dynamics in the presence of multiple factors?(2) What impact do control agents have within an integrated pest management system?(3) How can we better understand and predict pest/disease propagation patterns over time and space?and (4) What are the priority locations/regions for conducting risk assessments and limiting potential invasions of pests/diseases [253]?
Modelling for pest and disease management generally utilizes an inductive or deductive approach.In the inductive ("top-down") approach, pest/disease occurrence is linked to the prevailing bio-climatic variable in a specific area, and then used to predict the likelihood of that organism occurring in another area [121].Algorithms such as Maximum Entropy [254] and Genetic Algorithm for Rule Set Production [255] utilize the inductive approach.The deductive ("bottom-up") approach starts by modelling a species' response to climate variables and applies the model to predict the climate suitability of the organism.Inductive and deductive approaches can be combined using Climate and Population Modelling [256], which compares location function, and the NCSU (North Carolina State University)/APHIS (Animal and Plant Health Inspection Service) Plant Pest-Forecasting system (NAPPFAST) [257].
Case Study: Insect Life Cycle Modelling Software for Pest and Disease Management One example of a deductive approach is applied in the Insect Life Cycle Modelling (ILCYM) open-source software package, developed by the International Potato Centre (CIP) [168,258].It represents a breakthrough in conceptualization, development, and implementation of a tool for pest and disease management within the CGIAR.It uses non-linear functions [259,260] and includes stochastic functions for variability in individual development times within a population [260].ILCYM supports the analysis of life table data of insect species obtained from laboratory studies at constant temperatures and estimates, step-by-step, a mathematical expression to represent the species' development time, development rate, mortality, senescence, survival, and reproduction [121].Functions derived for each stage are combined to give a process-oriented, temperature-driven, and age-stage structured insect phenology/population model.Linkages between the phenology model and the landscape are based on geo-reference daily/monthly minimum and maximum temperature values using a cosine or sinus function to mimic patterns under natural conditions.Further, three risk indices (establishment, generation, and activity) are derived to assess the potential distribution and abundance of the pest in a given location [121,167,261].
Most of the tools developed to enhance management of pests/diseases require an estimate of model parameters [121].Efforts have been made to address this, but challenges can remain when selecting a mathematical equation for a modelling approach.Paradoxically, obtaining good fit does not automatically indicate that the selected mathematical equation or method accurately describes the dynamics and interactions between the crop and the pest/disease [121,262,263].Different model formulations may give excellent and equal fit to observed data, while their arithmetical characteristics produce dissimilar performance under field conditions [121].
Despite advances in modelling, the development and use of pest and disease models in isolation, as presented here, is not holistic and powerful enough.A reductionist approach, concentrating on a component rather than on the whole crop system [252] is not sufficient because, in reality, complex interactions exist between agricultural ecosystem components.Comprehensive modelling of crop development and yields requires multitrophic interactions, pest and disease impacts, ecosystem functions and processes to be combined, socio-ecological and socio-economic challenges, as well as respective feedback loops, biophysical flows, time delays, and nonlinearity within functional fluctuating bounds (Figure 2) [121].Overall, as a remedy to current approaches and methods, we propose to put the study of the whole crop system before that of the components to achieve sustainable management in crop production [121,252].

Figure 2.
The figure displays a system thinking approach in which the study of the whole crop system is placed before that of the components such as pest or disease.The approach links all system components to provide a comprehensive assessment of crop production within the ecosystem with emphasis on multitrophic interactions including delays and side effect taking into account socioeconomical and socio-ecological challenges.Figure was adapted with permission from H. Tonnang [121].

Crop Models for Management Practices
Crop simulation models have been successfully deployed in developed countries to manage farm resources including crop choice, crop variety, sowing date, plant density, irrigation, and inputs of nitrogen and phosphorus, but the quality of input data remains a challenge.Further, current problems related to sustainable intensification of agricultural systems cannot be solved by only optimizing these management practices, and many biophysical models lack components that computer metrics need to assess sustainable intensification related to soil health, for example [212].Most current crop/cropping system models lack the capability to simulate more sophisticated management aspects that may co-limit crop productivity.These include management options such as the use of disease-resistant varieties, modified tillage, soil additives, micronutrients, changing row geometries, and more precise timing and placement of fertilizers, drip irrigation, use of surface mulches for modification of crop micro-environments, and erosion.Some of these limitations have been addressed in specific model modifications, but are not generally included in the widely distributed versions of the models.Some of these model applications are discussed in detail in [116], and efforts to address these limitations are ongoing [252].

Case study: Exploring Sustainable Crop Management Options to Reduce Groundwater Table Decline in Northwest India
Irrigated rice (Oryza sativa L.)-wheat (Triticum aestivum L.) in northwest India is entirely dependent upon groundwater pumped from the underlying aquifer, resulting in alarming rates of

Biological and Agro-Ecological processes
Figure 2. The figure displays a system thinking approach in which the study of the whole crop system is placed before that of the components such as pest or disease.The approach links all system components to provide a comprehensive assessment of crop production within the ecosystem with emphasis on multitrophic interactions including delays and side effect taking into account socio-economical and socio-ecological challenges.Figure was adapted with permission from H. Tonnang [121].

Crop Models for Management Practices
Crop simulation models have been successfully deployed in developed countries to manage farm resources including crop choice, crop variety, sowing date, plant density, irrigation, and inputs of nitrogen and phosphorus, but the quality of input data remains a challenge.Further, current problems related to sustainable intensification of agricultural systems cannot be solved by only optimizing these management practices, and many biophysical models lack components that computer metrics need to assess sustainable intensification related to soil health, for example [212].Most current crop/cropping system models lack the capability to simulate more sophisticated management aspects that may co-limit crop productivity.These include management options such as the use of disease-resistant varieties, modified tillage, soil additives, micronutrients, changing row geometries, and more precise timing and placement of fertilizers, drip irrigation, use of surface mulches for modification of crop micro-environments, and erosion.Some of these limitations have been addressed in specific model modifications, but are not generally included in the widely distributed versions of the models.Some of these model applications are discussed in detail in [116], and efforts to address these limitations are ongoing [252].

Case Study: Exploring Sustainable Crop Management Options to Reduce Groundwater Table Decline in Northwest India
Irrigated rice (Oryza sativa L.)-wheat (Triticum aestivum L.) in northwest India is entirely dependent upon groundwater pumped from the underlying aquifer, resulting in alarming rates of reductions of the groundwater table and the potential for saline groundwater intrusion into fresh groundwater [264].There is therefore great interest in finding sustainable solutions to arrest the decline in groundwater depth.
The Agricultural Production Systems sIMulator (APSIM) crop model (defined in Section 4.2) was applied to evaluate the effects of different crop management practices (rice variety duration, sowing date, replacing rice with other crops, and the use of conservation agriculture) on the land and water productivity of the total system, and on the components of water balance (Table 3) [265,266].
The simulations gave the following results: • If minimizing permanent loss of water from the system is the objective, the best method is a partial conservation agriculture rice-wheat system with a short duration rice variety.

•
If maximizing productivity is the sole objective, a full conservation agriculture system with a medium duration variety is the best option.

•
If the objective is to produce the maximum output per unit of evapotranspiration, this can be achieved using a short duration variety planted in mid-July under full conservation agriculture [265,266].

•
Replacing rice with rainy season maize (Zea mays L.) in the system can maintain yields while reducing the total irrigation amount by ~80%, which can contribute to huge reductions in pumping costs and energy use.Maize-based systems also reduce permanent water loss from the system by 200 mm (Figure 3).
This approach can be extended to rainfed rice systems worldwide as a risk reducing strategy in a variable monsoon.given that testing all options in field is not feasible and are too expensive.Shaded bars are 25th-75th percentiles; whisker caps are 10th and 90th percentile and black dots 5th and 95th percentile, horizontal line is median.

E × M Interactions
Crop simulation models can effectively account for the yield and growth effects of the interaction of environment with crop management (E × M), assuming fixed G effects.One example is the AgroClimate program run by the University of Florida [267], which provides farmers with strategic advice regarding which crop to plant as a function of the El Niño-Southern Oscillation patterns.Probabilistic distribution of: (a) rice equivalent yield (REY), (b) irrigation, and (c) evapotranspiration of a rice-wheat system and maize-wheat system in central Punjab, India over 40 years .An example of how model outcome can be used to guide the field experimentation, given that testing all options in field is not feasible and are too expensive.Shaded bars are 25th-75th percentiles; whisker caps are 10th and 90th percentile and black dots 5th and 95th percentile, horizontal line is median.

E × M Interactions
Crop simulation models can effectively account for the yield and growth effects of the interaction of environment with crop management (E × M), assuming fixed G effects.One example is the AgroClimate program run by the University of Florida [267], which provides farmers with strategic advice regarding which crop to plant as a function of the El Niño-Southern Oscillation patterns.Another example is Yield Prophet ® , a web-based service that provides field-specific yield forecasts to optimize the investments of approximately 900 farmers, based on the Agricultural Production Systems sIMulator (APSIM) model, one of the major used crop models in agricultural systems [268].APSIM was originally developed to model biophysical processes of crop production in farming systems with the purpose of evaluating options of management under climate risk conditions, it has been broadly used in research and more applied questions, decision making on farm level, crop supply chain planning, evaluation of policy options and others.The recent introduction of Smartphone applications aims to reach even more farmers than conventional tools, such as the Yield Prophet Lite application [269] and PANI (Program for Advanced Numerical Irrigation) [270].However, the use of smartphone applications presents some challenges for some small farmers, such as the lack or unreliable internet connection, difficulties to recharge the mobile devices, illiteracy or poor education, high costs, among others.The RiceAdvice application, developed by Dr. Saito from AfricaRice for target yield, nutrient management, crop calendar and good agricultural practices, has partially solved some of these issues.The application runs offline, an important fact given poor connectivity and high internet prices in some areas in Africa, and an extensive training for service providers and users has been provided, thus overcoming some issues of illiteracy.
Case Study: Irrigation Scheduling in Bangladesh's Ganges Delta Region In Bangladesh's southern delta, farmers mostly grow grass pea (Lathyrus sativus L.) and mungbean (Vigna radiata L.) during the Rabi (dry) season (December-April) [271].These crops are generally not irrigated, instead drawing on residual soil moisture and the fact that the water table is close to the surface, fluctuating between 1-2.5 m depth in the dry season [272].In this environment, tactical irrigation is a useful and cheap option to increase water use efficiency, especially as these crops can be a major source of cash income for smallholder farmers, while significant yield losses are often caused by sub-optimal plant establishment due to dry conditions around sowing.
The International Maize and Wheat Improvement Center (CIMMYT) conducted several experiments on maize, mungbean, and wheat in the delta region to better understand irrigation needs and develop a decision support app for farmers and irrigation service providers.Figure 4 shows the simulated soil water content, with and without capillary upflow (a function of soil type, rooting depth, and distance to the water table), based on work by Talsma [273] and an algorithm developed by W. Meyer for the farm-level water balance model BASINMAN [274].The soil-water balance used in the smart phone application is based on the widely used CERES model, of which several specific crop variants are included in DSSAT (Decision Support System for Agrotechnology Transfer) [275], and can therefore be used under most conditions.Crop water use was estimated using percent ground cover, which can be estimated from optical remote sensing data [276,277].In Bangladesh, field sizes tend to be too small to make use of the free Sentinel-2 satellite data, given its limited resolution.An application is being developed that estimates percent ground cover from RGB photos taken with a smartphone.
addresses the needs of irrigation service provider and farmers, as both receive a weekly SMS informing them whether a field needs to be irrigated or not.PANI demonstrates that it is possible to bring crop simulation models, or components thereof, to the fingertips of smallholder farmers.It is possible to automatically run the model with representative soil and weather information based on the GPS location; the only inputs required from the user are sowing and irrigation dates and ground cover photos taken every 10-15 days.The green line represents the drained upper limit, while the red indicates the target stress level (50% of drained upper limit until anthesis and 30% thereafter).With a water table present, irrigations were required 7 and 32 days after planting, whereas without a water table, two more irrigations were required at 77 and 92 days after planting.The green line represents the drained upper limit, while the red indicates the target stress level (50% of drained upper limit until anthesis and 30% thereafter).With a water table present, irrigations were required 7 and 32 days after planting, whereas without a water table, two more irrigations were required at 77 and 92 days after planting.
The simulation model has been integrated into a smartphone application called PANI (Program for Advanced Numerical Irrigation), which runs on a daily time step and uses forecasted daily maximum and minimum temperatures to predict irrigation needs one week in advance.PANI addresses the needs of irrigation service provider and farmers, as both receive a weekly SMS informing them whether a field needs to be irrigated or not.PANI demonstrates that it is possible to bring crop simulation models, or components thereof, to the fingertips of smallholder farmers.It is possible to automatically run the model with representative soil and weather information based on the GPS location; the only inputs required from the user are sowing and irrigation dates and ground cover photos taken every 10-15 days.
Bio-physical models can provide M options at both macro-and micro-environmental levels [215].

M × S Interactions
Yield gap is defined as the difference between potential yield and the actual yield obtained in farmers' fields [278].Potential yield assumes that neither water, nutrients, nor biotic stress are limiting.Water-limited potential yield (instead of potential yield) is commonly used as a reference under rainfed conditions, and a further distinction is made for the yield potential of partially irrigated crops.Potential yield, water-limited potential yield, and partially irrigated potential yield can all be simulated in crop simulation models using optimal agronomic management (i.e., phenotype, sowing data, and planting density) as inputs [278].Neither nutrient nor abiotic stresses are considered.
Many models, such as the Global Yield Gap Atlas, base optimum sowing dates, plant density, and cultivar on the dominant practices used by farmers.Crop simulation models have the advantage of being able to assess the year-to-year variability of potential yield (e.g., [18,59]).Often, actual yield data are obtained for just one or a limited number of years, thus crop simulation models can be used to assess whether the year for which the yield gap was analysed was representative or not [279].More sophisticated analyses of e.g., nutrient-or pest-limited yield gaps, require specific inputs for the models, which often are not available in sufficient detail.
Yield gap models enable production constraints on farmers' fields to be identified [279].Region-or field-specific optimized crop management recommendations can then be generated using the example of farmers whose fields achieve the highest economic return.
Yield gap analysis has traditionally focused on agronomic crop production capacity and actual yields, but socioeconomic factors also play an important role.While the representation of yield gaps in Figure 5 ultimately involves G × E × M × S interactions, it highlights the interface between the biophysical and human or socioeconomic aspects of productivity.For any given E, M interventions are determined largely by S, since they must be both feasible and cost-effective for the farmer.Potential yield and water-limited yield represent G × E potential, depending on the water regime.Nutrient-limited yield and part of the environment-reduced yields (due to biotic and abiotic pressures) can often be overcome by M (through proper fertilization, integrated pest management, etc.), assuming such interventions are economically feasible and socially acceptable.
Many models, such as the Global Yield Gap Atlas, base optimum sowing dates, plant density, and cultivar on the dominant practices used by farmers.Crop simulation models have the advantage of being able to assess the year-to-year variability of potential yield (e.g., [18,59]).Often, actual yield data are obtained for just one or a limited number of years, thus crop simulation models can be used to assess whether the year for which the yield gap was analysed was representative or not [279].More sophisticated analyses of e.g., nutrient-or pest-limited yield gaps, require specific inputs for the models, which often are not available in sufficient detail.
Yield gap models enable production constraints on farmers' fields to be identified [279].Regionor field-specific optimized crop management recommendations can then be generated using the example of farmers whose fields achieve the highest economic return.
Yield gap analysis has traditionally focused on agronomic crop production capacity and actual yields, but socioeconomic factors also play an important role.While the representation of yield gaps in Figure 5   Actual yield is often limited by many interacting factors.Attempts have been made to incorporate this by introducing the concept of attainable yield, and the attainable yield gap.This is defined in terms of the best local [280] or global [281] practices, under similar biophysical and climatological circumstances.
Figure 5. Graphical presentation of the various factors contributing to the yield gap. Actual yield is often limited by many interacting factors.Attempts have been made to incorporate this by introducing the concept of attainable yield, and the attainable yield gap.This is defined in terms of the best local [280] or global [281] practices, under similar biophysical and climatological circumstances.
Attainable yield levels are estimated using bio-economic models that incorporate household and community characteristics, and calibrated household-level objective functions, including risk management and imperfect information.In practice, farmers' actual yields may diverge from predicted yields due to idiosyncratic circumstances not accounted for in the model, such as the occurrence of household-level catastrophes.Models can be calibrated to take farm-specific circumstances into account and-if the level of yield reporting is sufficient and readily quantifiable G × E × M × S interactions are included-the difference between actual farmers' yields and modelled yields will be small.Bio-economic modelling can provide a set of M options relevant under different S scenarios [196,282].Bio-economic modelling covers the ex-ante impact assessment of new technologies (including new varieties), combining socio-economic insights with explicit biological and biophysical process models.Extensive reviews of bio-economic models have been published (see [196,[283][284][285][286][287][288][289][290]) and more detailed information about these models is described in Section 6.1.Bio-economic modelling at the farm level can take one of several approaches, some of which make assumptions about farmer behaviour and use these assumptions to analyse trade-offs (e.g., [5,[291][292][293][294][295][296][297]).
Recent research by van Dijk et al. [298] uses econometric analyses at the household level to demonstrate how some of these issues might be resolved by integrating M and S approaches to yield gap measurement.The authors deconstruct the conventional yield gap indicator into four components to provide a better understanding of why actual farm yields fall below their potential: (1) the technical efficiency yield gap; (2) the economic yield gap; (3) the feasible yield gap; and (4) the technology yield gap ( [299,300]).These components are an alternative way of identifying the various components contributing to the yield gap as highlighted in Figure 5.
Building on these insights, and using ex-ante bio-economic modelling approaches, M × S interaction models can explore the cost-benefit of new technologies under different socioeconomic and climate scenarios, if interoperable agronomic, genetic, and socioeconomic datasets are available.Further issues related to the data availability are discussed in Section 6.The approach also requires more widely accessible bio-economic models at various levels of aggregation, which combine the biophysical process models with behaviour-based economic models, preferably as open-source models that allow communities of practice to be established [5].

Technological Innovations Aiding Crop Management Models
Remote sensing has the potential to characterize growth at high-throughput over space and time, and is thus a valuable tool for extending the inputs and calibration of crop simulation models.Optical data provide good estimates of crop light interception [276,277] and have been successfully used to drive crop simulation models for yield and biomass predictions.Remote sensing can also estimate stress parameters; the widely used Normalized Difference Vegetation Index (NDVI) can indicate the cumulative effects of stress throughout the crop growth cycle (especially when referenced against non-stressed control plots).Other indices, such as the photochemical reflectance index [301] or chlorophyll fluorescence indices [302], can estimate the instantaneous effect of stress on radiation use efficiency.Satellite data can also detect key phenological stages, such as crop emergence, end of leaf growth in cereals, and maturity.This information reduces the burden of having to gather accurate crop management information from farmers, which is otherwise a major bottleneck for the widespread use of models.
Case Study: Anomalies in Irrigation Experiments in Maize Farmer's Fields Detected with Remote Sensors As an example, remote sensing helped better understand perceived anomalies in an irrigation experiment conducted in maize farmers' fields in Bangladesh as part of a project developing PANI (see Section 4.2).In-season measurements of canopy temperature and ground cover showed that there were larger differences in crop growth within irrigation treatments than among them (Figure 6).This is because farmers intentionally create a relatively water-impermeable soil layer (a "hard pan") at a depth of 20-30 cm when cultivating fields for summer-season rice planting which affect root penetration and therefore the expression of canopy temperature.The land preparation for the maize irrigation experiment consisted of one pass with a power tiller and subsequent shaping of the beds.This was evidently not sufficient to brake the hard-pans in some of the previous rice paddy fields, as can be seen for the central plot of the irrigation experiment in Figure 6.The upper part suffered from dry conditions, which resulted in lower ground cover and higher canopy temperature as compared to the lower part of the plot.Soil measurements indicated the presence of hard pans related to the previous cropping history, and showed higher moisture availability below the hard pan in experimental plots with severe drought symptoms, indicating that plants in these plots were not able to penetrate the hard pan (Schulthess et al., pers. comm.).
fields, as can be seen for the central plot of the irrigation experiment in Figure 6.The upper part suffered from dry conditions, which resulted in lower ground cover and higher canopy temperature as compared to the lower part of the plot.Soil measurements indicated the presence of hard pans related to the previous cropping history, and showed higher moisture availability below the hard pan in experimental plots with severe drought symptoms, indicating that plants in these plots were not able to penetrate the hard pan (Schulthess et al., pers.comm.). Figure 6.Plots of irrigation experiment near Patuakhali, Bangladesh, in 2015.The nine maize plots had been imposed over several farmers' fields.The hard pans of some of the fields, created for the previous paddy rice crop, were not broken up when preparing the land for the maize experiment.This can be seen for the central plot, where plants suffered from dry conditions in the upper part, resulting in lower ground cover and higher canopy temperature as compared to the lower part of the plot.
The availability of remote sensing data has been greatly improved by recent launches of new public and private remote sensing satellite systems.The combined use of remote sensing and crop simulation models will make management recommendations more dynamic, with the ability to account for actual and forecasted weather conditions (including the prediction of outcomes in different weather years).Smart phones can also help gather information on crop management [5], which then can be used to parameterize model runs.It is already theoretically possible to do this, but practical hurdles remain.The stakes are high, as the input for all model components needs to be accurate and the simulation results reliable.
There is the potential to use specific model components to generate recommendations regarding the timing of crop management activities, i.e., whether nitrogen top dressing or irrigation needs to be conducted earlier, at the same time, or later than usual, as function of actual weather conditions.Similarly, the abovementioned PANI recommendation can be turned into irrigation scheduling recommendations, informing the farmer whether a crop needs to be irrigated in shorter or longer intervals than usual.These kinds of recommendations can be generated without any input from the farmers.The only requirements are a properly calibrated and validated model, and historic and current weather data.
Janssen et al. [303] discuss integrating emerging information and communication technologies with big data sources, and how this can be incorporated into farmer or industry decision making using crop models.The increased access of rural areas to Information Communication Technologies (ICT), reduced costs in satellite imagery acquisition, and cloud-or supercomputer-based analysis for Figure 6.Plots of irrigation experiment near Patuakhali, Bangladesh, in 2015.The nine maize plots had been imposed over several farmers' fields.The hard pans of some of the fields, created for the previous paddy rice crop, were not broken up when preparing the land for the maize experiment.This can be seen for the central plot, where plants suffered from dry conditions in the upper part, resulting in lower ground cover and higher canopy temperature as compared to the lower part of the plot.
The availability of remote sensing data has been greatly improved by recent launches of new public and private remote sensing satellite systems.The combined use of remote sensing and crop simulation models will make management recommendations more dynamic, with the ability to account for actual and forecasted weather conditions (including the prediction of outcomes in different weather years).Smart phones can also help gather information on crop management [5], which then can be used to parameterize model runs.It is already theoretically possible to do this, but practical hurdles remain.The stakes are high, as the input for all model components needs to be accurate and the simulation results reliable.
There is the potential to use specific model components to generate recommendations regarding the timing of crop management activities, i.e., whether nitrogen top dressing or irrigation needs to be conducted earlier, at the same time, or later than usual, as function of actual weather conditions.Similarly, the abovementioned PANI recommendation can be turned into irrigation scheduling recommendations, informing the farmer whether a crop needs to be irrigated in shorter or longer intervals than usual.These kinds of recommendations can be generated without any input from the farmers.The only requirements are a properly calibrated and validated model, and historic and current weather data.
Janssen et al. [303] discuss integrating emerging information and communication technologies with big data sources, and how this can be incorporated into farmer or industry decision making using crop models.The increased access of rural areas to Information Communication Technologies (ICT), reduced costs in satellite imagery acquisition, and cloud-or supercomputer-based analysis for faster processing and delivery (e.g., Descartes Labs daily image analysis), improves the ability of farmers and technicians to use crop modelling tools for day-to-day management decisions.

Crop Simulation Models for Genetic Improvement
Analysis of interactions between G and E, and between G and E and M, requires a deep understanding of their biological basis.Mechanistic crop models such as APSIM [304] or Simple Simulation Modelling (SSM) [305] can predict crop performance, based on equations characterizing the responses of crop growth and development to environmental cues.These equations are derived from experimental data, making these models both process-and statistics-driven.Therefore, APSIM and SSM provide accurate predictions of crop production in relation to genotype (G), environment (E), soil and management factor (M), while addressing the long-term resource management issues.
Models like APSIM and SSM have successfully predicted growth and yield in chickpea (Cicer arietinum L.) [143,145,306], lentil (Lens culinaris L.) [205,206], bean (Phaseolus vulgaris L.) [307], soybean (Glicine max (Linn.)Merr.) [152], peanut (Arachis hypogaea L.) (Vadez et al., pers.comm.), sorghum (Sorghum bicolor L.) [308], wheat [309], and maize [310][311][312].Models can be used to assess the intensity and frequency of water stress patterns, enabling more targeted breeding and agronomic management interventions for specific environments [313,314].For example, Kholova et al. [151] used APSIM to demonstrate that several genetic traits underlying a staygreen phenotype in sorghum increased grain yield across sorghum production regions in India (Figure 7).Some of these trait modifications resulted in trade-offs between grain and stover productivity.For example, breeding for a smaller leaf canopy reduced pre-anthesis water use, which improved water availability for grain filling and ultimately increased yields, but was at the expense of stover productivity, which is unfavourable for farmers interested in crop residues for fodder [151].

Crop Simulation Models for Genetic Improvement
Analysis of interactions between G and E, and between G and E and M, requires a deep understanding of their biological basis.Mechanistic crop models such as APSIM [304] or Simple Simulation Modelling (SSM) [305] can predict crop performance, based on equations characterizing the responses of crop growth and development to environmental cues.These equations are derived from experimental data, making these models both process-and statistics-driven.Therefore, APSIM and SSM provide accurate predictions of crop production in relation to genotype (G), environment (E), soil and management factor (M), while addressing the long-term resource management issues.
Models like APSIM and SSM have successfully predicted growth and yield in chickpea (Cicer arietinum L.) [143,145,306], lentil (Lens culinaris L.) [205,206], bean (Phaseolus vulgaris L.) [307], soybean (Glicine max (Linn.)Merr.) [152], peanut (Arachis hypogaea L.) (Vadez et al., pers.comm.), sorghum (Sorghum bicolor L.) [308], wheat [309], and maize [310][311][312].Models can be used to assess the intensity and frequency of water stress patterns, enabling more targeted breeding and agronomic management interventions for specific environments [313,314].For example, Kholova et al. [151] used APSIM to demonstrate that several genetic traits underlying a staygreen phenotype in sorghum increased grain yield across sorghum production regions in India (Figure 7).Some of these trait modifications resulted in trade-offs between grain and stover productivity.For example, breeding for a smaller leaf canopy reduced pre-anthesis water use, which improved water availability for grain filling and ultimately increased yields, but was at the expense of stover productivity, which is unfavourable for farmers interested in crop residues for fodder [151].i.e., more in the northern sub-region where the prevalence of severe stress is higher, less so in the two southern sub-regions where the prevalence of mild stress is higher.Adapted from Kholova et al. [151]., more in the northern sub-region where the prevalence of severe stress is higher, less so in the two southern sub-regions where the prevalence of mild stress is higher.Adapted from Kholova et al. [151].
In another modelling study, Vadez et al. [143] compared the effects of altering root-related genetic traits on the grain yield of chickpea and found that, while the best genetic trait increased yields by 10%, a modest irrigation of 30 mm at the beginning of pod filling increased yields by 40%.This illustrates the potential of mechanistic crop models to compare genetic and agronomic options and their importance as decision-making tools in research investments.
Finally, crop models have been used to undertake yield gap analysis and lay out potential options for improving crop productivity at a country level.For instance, Soltani et al. [157] showed that chickpea production in Iran could be tripled by using short duration cultivars and applying a single irrigation of 60 mm at the beginning of pod filling.

Genomic Selection Models Incorporating Physiological Traits and Insights
Genomic selection exploits the within family variance by means of genome-wide molecular markers/sequence data without the need to phenotype the candidates for selection; this saves time and enhances genetic gains and the response to selection per unit of time.While genomic selection models do not require a physiological understanding of the traits of interest, understanding the physiological underpinnings of grain yield can enable the development of new genomic selection approaches with improved accuracies, both within and across environments.A relatively simple way to utilize data on physiological traits is to use multivariate mixed models for prediction.These models enable information to be "borrowed" from relatives and secondary traits that are genetically correlated with the trait of interest.Multivariate prediction models have outperformed univariate prediction models in both animal and plant species [315,316], especially when secondary traits are observed on both the model training and validation sets [317].
The use of multivariate models with aerial measurements of canopy temperature and NDVI as secondary traits on both training and validation sets increased pedigree and genomic selection accuracy for grain yield breeding in optimal, severe drought, and late heat environments in Obregon, Mexico, for both replicated and unreplicated secondary traits (Figure 8) [115].Across five environments, multivariate models outperformed univariate models by up to 70% when secondary trait data was not replicated.
In another modelling study, Vadez et al. [143] compared the effects of altering root-related genetic traits on the grain yield of chickpea and found that, while the best genetic trait increased yields by 10%, a modest irrigation of 30 mm at the beginning of pod filling increased yields by 40%.This illustrates the potential of mechanistic crop models to compare genetic and agronomic options and their importance as decision-making tools in research investments.
Finally, crop models have been used to undertake yield gap analysis and lay out potential options for improving crop productivity at a country level.For instance, Soltani et al. [157] showed that chickpea production in Iran could be tripled by using short duration cultivars and applying a single irrigation of 60 mm at the beginning of pod filling.

Genomic Selection Models Incorporating Physiological Traits and Insights
Genomic selection exploits the within family variance by means of genome-wide molecular markers/sequence data without the need to phenotype the candidates for selection; this saves time and enhances genetic gains and the response to selection per unit of time.While genomic selection models do not require a physiological understanding of the traits of interest, understanding the physiological underpinnings of grain yield can enable the development of new genomic selection approaches with improved accuracies, both within and across environments.A relatively simple way to utilize data on physiological traits is to use multivariate mixed models for prediction.These models enable information to be "borrowed" from relatives and secondary traits that are genetically correlated with the trait of interest.Multivariate prediction models have outperformed univariate prediction models in both animal and plant species [315,316], especially when secondary traits are observed on both the model training and validation sets [317].
The use of multivariate models with aerial measurements of canopy temperature and NDVI as secondary traits on both training and validation sets increased pedigree and genomic selection accuracy for grain yield breeding in optimal, severe drought, and late heat environments in Obregon, Mexico, for both replicated and unreplicated secondary traits (Figure 8) [115].Across five environments, multivariate models outperformed univariate models by up to 70% when secondary trait data was not replicated.As described above, physiological trait data could also be potentially used in prediction models that incorporate crop growth models and genome-wide markers to predict breeding candidate performance in untested environments.Physiological trait data has not yet been utilized in such As described above, physiological trait data could also be potentially used in prediction models that incorporate crop growth models and genome-wide markers to predict breeding candidate performance in untested environments.Physiological trait data has not yet been utilized in such models because the traits utilized in crop growth models are relatively difficult to measure on many individual plants.
Rather than combine crop growth models and genome-wide markers in one, unified model, Heslot et al. [318] used a crop model and weather data to derive stress covariates by growth stage.They then used these stress covariates, along with genome-wide markers, to explain the residuals from genomic prediction using a machine-learning approach.If sufficient phenotypic data were available such that crop models could be parameterized for every individual, it would be possible to compute growth stages for each individual in each environment, perhaps improving predictive abilities.

G × E Interactions
Physiological breeding crosses parents with different complex traits to accumulate alleles associated with yield increases [7].The strategy is based on testing models of theoretically complementary traits and evaluating selected progeny in multi-location yield trials to inform G × E models.Analysing the performance of new lines developed via physiological breeding helps inform crop design and refine experimental environments, demonstrating which combinations of traits/alleles improve yield, and in which environments [319].Physiological and genetic dissection of complex traits provides valuable insights; for example, by establishing the genetic bases for trait synergies and trade-offs (e.g., [320,321]).Such outputs can inform simulation models that test hypotheses regarding new levels of trait expression and trait combinations across a range of environments [145,322].
Many statistical models have been developed to study G × E for plant breeding.Unlike simulation or process driven models-that are driven by environmental inputs and modified by genotype information-, statistical models are entirely stochastic.Outputs are based on associations among genetic and environmental variables with probability largely determining the significance of conclusions.Statistical models are very flexible and can be designed for almost any data set assuming that error (and other terms as in regression or multivariate analyses) can be modelled effectively.
These models have helped plant breeders assess the stability of economically important traits and predict the performance of newly developed genotypes under various environmental conditions.The development of methods for modelling G × E precedes the development of analyses of variance, as it could be better modelled by fitting a multiplicative (product) operator than by a sum formula [323].A simple regression of a line's performance on the environmental mean was first proposed and used by Yates and Cochran [324].Other models to study G × E are based on the principal component decomposition of matrix G × E [325][326][327][328][329]. Researchers employed the principal component operator for modelling G × E in the context of linear mixed-effect models [330][331][332], while others evaluated using structured co-variance matrices to model G × E in the context of pedigree-based mixed models for estimating the breeding values [333][334][335].
The models mentioned above do not attempt to identify genetic causes of G × E in the chromosome by decomposing total G × E into chromosome regions or specific genes.However, differences between individuals can be explained by regressing phenotypes on molecular markers, and differences between environments can be evaluated by regressing phenotypes on a reduced number of environmental covariates (e.g., temperature, precipitation, soil moisture, solar radiation).By combining information on individual genomes (e.g., molecular markers) with environmental covariables, it is possible to account for G × E.An example of such an approach is the ordinary least squared factorial regression model [336][337][338][339].
One advance in using statistical approaches in breeding models was developed by Yin et al. [340].Their concept was that variety-specific parameters could be linked to genetic information through QTL (quantitative trait loci) and eventually a specific part of the genome.In breeding programs, varieties can be characterized using QTLs, enabling G × E to be predicted through crop models.Working with the common bean, Hwang et al. [341], established dynamic QTL effect models for node addition rate, daily rate of progress from emergence toward flowering, and daily maximum main stem node number.
As modern genotyping and sequencing technologies develop, molecular marker information has become highly dimensional.The number of markers now greatly exceeds the number of data-points available for model fitting, posing important statistical and computational challenges.However, whole-genome regression methods allow information about hundreds of thousands of markers to be jointly considered [342].Such models capture major-effect genes, as well as genomic regions with small effects.Environmental information is also becoming highly dimensional due to information systems that allow continuous monitoring of environmental climatic conditions.Previous models cannot cope with the highly dimensional nature of genomic and E information.
Standard genomic selection models were recently extended to multi-environment settings.For instance, Burgueño et al. [343] used a multi-environment version of the Genomic BLUP (GBLUP), modelling G × E using genetic correlations, and found that the multi-environment GBLUP had a much higher prediction accuracy than the single-trait GBLUP.Meanwhile Jarquin et al. [344] proposed a model where the main effects of, and interactions between, markers and environmental covariates are introduced.Recently, Lopez-Cruz et al. [345] proposed a marker × environment interaction model where marker effect and genomic values are partitioned into components that are stable across environments (main effects), and others that are environment-specific (interactions).This interaction model is useful when selecting for stability and adaptation to target environments.Genomic prediction accuracy consistently and substantially increased when incorporating G × E and marker × environment interactions.

G × E × M Interactions
G × E × M interactions are difficult to model with precision, particularly when extrapolating to new environments and genotypes.The challenges have been highlighted by the work of the AgMIP Wheat group, who showed the uncertainty encountered when comparing 30 different models; nevertheless, multi-model ensembles with more robust predictions have been developed [346].A powerful innovation will be to encompass the following types of approaches into crop simulation modelling for more robust outputs:

•
Statistical modelling of G × E × M (and lower order interactions) as a stochastic tool for testing specific elements of simulation models and to parameterize or calibrate input variables; • Genomic selection models using high density molecular markers to parameterize G and thus test the value of specific genomic regions and make predictions in different target environments [115]; • Physiological breeding models to test different combinations of traits and alleles through strategic crossing and progeny testing in well-defined target environments [7]; Testing of G × M interactions, which lend well to reductionist research approaches because-unlike E, which is infinitely variable in the field-M can be varied for one or more factor within a single E.
Montesinos-Lopez et al. [347] recently developed the theory of a Bayesian three-way, multi-trait, multi-environment genomic model.The proposed Bayesian model (with unstructured variance-covariance) was evaluated on two datasets and improved prediction accuracy, compared to the model with diagonal and identity variance-covariance structures [347].The three-way Bayesian genomic interaction model uses genomic and/or pedigree information and could be extended to include more factors, potentially making it applicable for assessing G × E × M and G × E × M × S in the context of crop modelling.
Intercropping is a well-known practice in Africa that is reported to have advantages over mono cropping through several mechanisms.To analyse the performance of intercropping systems, crop simulation models can be helpful to identify those mechanisms in such a complex cropping system.The first ecophysiological model for interplant competition INTERCOM was developed by Kropff and van Laar [17] at IRRI and Wageningen University and was well documented with first applications in crop weed systems.The model has subsequently been used in several systems for intercropping in Africa, Switzerland, China and the Netherlands.In Africa, with AfricaRice and Wageningen University, the model was used to identify the optimum relay crop in rice systems with a highly competitive (Cajanus c.) and a not competitive species (Stylosanthus spp.)where it was found that relay cropping with the competitive species was highly beneficial when late sown.This was a very interesting counter intuitive result [60,62].Intercropping model versions have been developed by the APSIM teams as well and have been used for agroforestry systems [348].However intercropping modelling has not moved to the level of intensively used decision support systems yet.As solutions are being sought for many issues in relation to sustainable agricultural systems in Africa in relation to climate change and new upcoming pests and diseases, the intensification of the use of intercrops and research to improve them through crop models such as INTERCOM may be very useful to enhance the agronomic research effort.

Bio-Economic Modelling
Bio-economic modelling dates to the 1990s and includes mathematical programming models [284,294,295,349,350], multi-agent models [351][352][353], and other simulation models [291,354].These models are sometimes combined with econometric analyses that fine-tune the results.Linking the results of crop models to these lower aggregation level models requires the determination of relevant adjustment factors, depending on the issues being addressed.For instance, if biotic stress in the analysis is linked to certain weather patterns and stochastic probabilities of occurrence, the crop growth model output needs to be combined with disease progress curves and damage functions to adequately account for biotic stress.While many of these bio-economic models focus on crop production, a number also take into account crop-livestock interactions [195,197,294,355].
Household survey data on actual yields and management practices can be used to develop specific, reduced form models that estimate yield reducing factors for non-optimal input levels.However, the use of reduced form models in econometric analyses of overall farm household behaviour poses difficulties for soft-coupling with crop model outcomes, and may be better analysed using structural models.The difference between econometrically estimated models of farm household behavior and structural models is that the structural models such as the bio-economic models highlighted in the previous paragraph, describe the behavior as structural equations in mathematical terms, whereas the reduced form models provide only the statistical associations between measured inputs into the decision making process and the measured outputs in terms of key indicators [282].
At higher levels of aggregation, integrated assessment models can illuminate interactions between climate and socioeconomic global drivers.Models such as IMPACT [86,356], the International Agricultural Prospects Model [357], and MAGNET [358,359] have demonstrated their usefulness in exploring these interactions.Global economic models have been assessed within AgMIP [2,90,360,361] and this approach will soon be applied to micro-level models.This allows models to be improved and fine-tuned to better answer questions related to the socioeconomic factors associated with technology change in agriculture.Within the CGIAR for example, a joint study by CIMMYT and IFPRI [89] implemented a spatially-explicit bio-economic modelling framework, incorporating the CERES-Wheat model in DSSAT and a farm-level profitability analysis of local wheat production.This integrated framework assessed the economic viability of potential wheat expansion in 12 countries across sub-Saharan Africa, and found that they are using less than 10% of their potential area for profitable local production.This type of approach can be utilised on a regional or global scale by agricultural ministries and funding agencies to help determine potential returns on investments.

Foresight Analysis
Foresight analysis is an important instrument for monitoring and enhancing our understanding of the evolving context in which agricultural systems operate.Given the research-to-adoption lag, the needs of future beneficiaries should inform current priorities, while also considering likely biophysical, socioeconomic, and politico-institutional contexts at the time when technologies (including varietal improvement) come to fruition.Both beneficiaries' needs and the context in which they operate are shaped by global drivers of change, pressures, and events.The primary global drivers of change are climate change, demographic developments, and global economic development.Secondary drivers of change include migration, urbanization, shifting gender roles, resource depletion, technological convergence, and changing political environments.Pressures resulting from these global drivers include pollution; decreases in natural ecologies; land degradation (and the resulting changes in environmental productivity); health related issues, from chronic malnutrition to obesity; and unequal wealth and income distribution.Events are major occurrences with global or regional effects.These events are difficult to accurately predict, but their occurrence is both probable and plausible at any given time.They can include weather extremes (e.g., droughts, floods, and the El Niño phenomenon); geophysical disasters (earthquakes and volcanic eruptions); outbreaks of new biotic stresses (such as recent outbreaks of maize lethal necrosis and rust in east Africa, and wheat blast in Bangladesh); economic shocks (such as steeply rising food prices); and political and social unrest.One positive factor that holds the middle ground between a driver and an event is disruptive technology.A clear example of a disruptive technology is the Green Revolution, which lies at the foundation of the CGIAR.
These drivers, pressures, and events are forcing farmers, civil societies, and countries to reassess priorities and interventions.The primary effect is the change in agro-ecological production potential (related to climate change) and the comparative advantage of different crops in different locations (linked to long-term trends in commodity prices).
Using a wide range of models and approaches across different scientific disciplines when analysing of the effects of global drivers, pressures, and events, improves our ability to realize potential opportunities, while responding in a timely fashion to avoid major food security crises.For example, to assess the potential role of agricultural technologies and innovations in foresight modelling studies, IFPRI's IMPACT model incorporates the DSSAT Cropping System Model to simulate grid cell-specific crop yield responses to scenarios of changes in technology adoption and environmental conditions, and uses this information to simulate their long-term implications on international food trade, food security, and markets [95].
Rosegrant et al. [92] used this analytical framework to assess the potential role of 12 promising technologies on international food trade, food security, and natural resources in 2050.The authors developed an integrated economic modelling framework with DSSAT, where the models of maize, rice, and wheat were modified to simulate technologies not previously implemented in the model (e.g., effects of anthesis-silking interval to simulate drought tolerance trait of maize).Another example is the analysis of climate change on future maize production.Tesfaye et al. [112] showed that climate change will affect maize yields across sub-Saharan Africa in the future (2050 and 2080), and that the extent of the impact at a given period will vary considerably between input levels, regions, and maize mega-environments.These effects may have a negative impact on maize consumption and daily calorie intake, though Rosegrant et al. [356] showed that international trade has the potential to offset some of the negative impacts.These analyses that build on the IMPACT model and use gridded crop models [93,362], in combination with socioeconomic models, have helped focus breeding efforts on specific abiotic stresses.

Case Study: Calculating Food Systems Risks
Global insurance company Lloyd's of London released a peer-reviewed, industry approved 2015 study on the impacts of acute disruption to the global food supply [363].Experts considered how weather-related events can disrupt the food supply, and can be exacerbated by additional political and economic risks in a chain-like reaction.Their scenario was based on a strengthened warm phase of the El Niño climate pattern, which would generate flooding and drought in different locations across the globe, as well as winds that spread the Ug99 wheat stem rust across a wide territory (Figure 9).Shock scenarios were developed and explored using the IFPRI IMPACT model, which gives a rough interpretation of weather changes leading to drought at a large grid level.Lloyd's model was constructed by de-trending globally aggregated FAO data from 1961 to 2013 to identify probable changes in crop area and yield.It should be noted that the estimated production response to the shock was deemed drastic but plausible by the authors (Further detail on assumptions about responses at the country level is included in the report, including production loss, trade restriction mechanisms, social instability, and currency fluctuations.The report authors note that "the probability of occurrence is estimated as significantly higher than the benchmark return period of 1:200 years applied for assessing insurers' ability to pay claims against extreme events", ( [363], p. 19)).The effects emulate recent country-level effects on major crops, such as maize in southern Africa, where major yield losses have been observed.An improved modelling system with more detailed grid data and more sophisticated crop models could deliver better insight, though more data on soils and weather conditions are essential to make progress in the explorative capacity of the approaches.constructed by de-trending globally aggregated FAO data from 1961 to 2013 to identify probable changes in crop area and yield.It should be noted that the estimated production response to the shock was deemed drastic but plausible by the authors (Further detail on assumptions about responses at the country level is included in the report, including production loss, trade restriction mechanisms, social instability, and currency fluctuations.The report authors note that "the probability of occurrence is estimated as significantly higher than the benchmark return period of 1:200 years applied for assessing insurers' ability to pay claims against extreme events," ([363], p. 19).).The effects emulate recent country-level effects on major crops, such as maize in southern Africa, where major yield losses have been observed.An improved modelling system with more detailed grid data and more sophisticated crop models could deliver better insight, though more data on soils and weather conditions are essential to make progress in the explorative capacity of the approaches.

Figure 9.
Plausible scenario of a food systems shock.The Lloyd's report aids crop modellers by clarifying the likely impacts of a plausible sudden food supply shock due to simultaneous occurrence of both biotic and abiotic stress on supply chains and societies [363].Adapted from Lloyd's report with the permission from the authors.
According to the IMPACT bio-economic model, production shocks would generate global losses of 7% for wheat, 10% for maize, 11% for soybeans, and 7% for rice (Figure 10).This would generate a quadrupling of prices for the first three crops, and a five-fold increase for rice.The impacts would extend deep into society, creating humanitarian crises where food shortages arise.Moreover, stock markets would be de-stabilized, with their values estimated to fall by 10% in Europe and 5% in the USA.
While Lloyd's report is a forecasting exercise, not a prediction, it underscores the fragility of the food system, its global interconnectedness, and its central role in maintaining social stability.Moreover, the complexity of socio-ecological systems confirms that crop modellers have a key role to play in understanding and preparing for major food supply disruptions.According to the IMPACT bio-economic model, production shocks would generate global losses of 7% for wheat, 10% for maize, 11% for soybeans, and 7% for rice (Figure 10).This would generate a quadrupling of prices for the first three crops, and a five-fold increase for rice.The impacts would extend deep into society, creating humanitarian crises where food shortages arise.Moreover, stock markets would be de-stabilized, with their values estimated to fall by 10% in Europe and 5% in the USA.While Lloyd's report is a forecasting exercise, not a prediction, it underscores the fragility of the food system, its global interconnectedness, and its central role in maintaining social stability.Moreover, the complexity of socio-ecological systems confirms that crop modellers have a key role to play in understanding and preparing for major food supply disruptions.
Figure 10.Supply chain and societal impacts of a plausible food systems shock.The Lloyd's study [363] found that the impact of such a weather event would be serious and widespread.Adapted from Lloyd's report with the permission from the authors.

Global Phenotyping Networks
To be credible and useful, crop simulations must hold true across a wide range of environments.They would therefore benefit immeasurably from global phenotyping networks that characterize breeding materials in key ecologies using harmonized protocols [13].Such networks can capitalize on existing infrastructure, but there is often a need for investment to ensure greater control of experimental treatments (Table 2) and of trait measurement (e.g., [364]).A networking approach has already been applied to initiatives such as Drought Tolerant Maize for Africa and Improved Maize for African Soils, and in the use of field phenotyping tools provided by various organizations and developed for ground and aerial sensing (e.g., [365]).Various phenotyping initiatives have emerged recently (e.g., International Plant Phenotyping Network, European Plant Phenotyping Network, North American Plant Phenotyping Network, Wheat Initiative, etc.) that could help promote common standards, at least within crops.The added value of high-throughput phenotyping methods (e.g., [366]) is recognized, as is their generic applicability for measuring many high priority traits across major cereal crops [250].As a relatively new field, high-throughput phenotyping lends itself well to the introduction of common standards.Some generic standards for field phenotyping platforms to support simulation modelling are presented in Table 2.

Data Sharing and Standardization
All the modelling approaches outlined here could be strengthened by validation from independent datasets.Many institutions-including CGIAR-are increasingly obliged to openly share research data.Aside from helping avoid redundancy and adding robustness to datasets, more widespread sharing of data enables modelling to take place on a larger scale (e.g., [367]).Broad agreement on core datasets, measurement protocols, and standardization of experimental treatments and definitions could exponentially increase the power of shared data [12,250].To simulate plant growth in models, the minimum information required includes daily weather data, soil characterization, basic agronomic inputs, and at least two measures of biomass, while more detailed and time-consuming measurements vary from model to model [250].Better environmental data is required from many more countries if the effects of climate change are to be predicted [368].
The limited availability of large-scale, interoperable datasets is recognized as bottleneck to the use of crop modelling and, more broadly, agricultural systems modelling [5].CGIAR centres manage Figure 10.Supply chain and societal impacts of a plausible food systems shock.The Lloyd's study [363] found that the impact of such a weather event would be serious and widespread.Adapted from Lloyd's report with the permission from the authors.

Global Phenotyping Networks
To be credible and useful, crop simulations must hold true across a wide range of environments.They would therefore benefit immeasurably from global phenotyping networks that characterize breeding materials in key ecologies using harmonized protocols [13].Such networks can capitalize on existing infrastructure, but there is often a need for investment to ensure greater control of experimental treatments (Table 2) and of trait measurement (e.g., [364]).A networking approach has already been applied to initiatives such as Drought Tolerant Maize for Africa and Improved Maize for African Soils, and in the use of field phenotyping tools provided by various organizations and developed for ground and aerial sensing (e.g., [365]).Various phenotyping initiatives have emerged recently (e.g., International Plant Phenotyping Network, European Plant Phenotyping Network, North American Plant Phenotyping Network, Wheat Initiative, etc.) that could help promote common standards, at least within crops.The added value of high-throughput phenotyping methods (e.g., [366]) is recognized, as is their generic applicability for measuring many high priority traits across major cereal crops [250].As a relatively new field, high-throughput phenotyping lends itself well to the introduction of common standards.Some generic standards for field phenotyping platforms to support simulation modelling are presented in Table 2.

Data Sharing and Standardization
All the modelling approaches outlined here could be strengthened by validation from independent datasets.Many institutions-including CGIAR-are increasingly obliged to openly share research data.Aside from helping avoid redundancy and adding robustness to datasets, more widespread sharing of data enables modelling to take place on a larger scale (e.g., [367]).Broad agreement on core datasets, measurement protocols, and standardization of experimental treatments and definitions could exponentially increase the power of shared data [12,250].To simulate plant growth in models, the minimum information required includes daily weather data, soil characterization, basic agronomic inputs, and at least two measures of biomass, while more detailed and time-consuming measurements vary from model to model [250].Better environmental data is required from many more countries if the effects of climate change are to be predicted [368].
The limited availability of large-scale, interoperable datasets is recognized as bottleneck to the use of crop modelling and, more broadly, agricultural systems modelling [5].CGIAR centres manage data from multiple crops and locations, much of which is collated within the AgTrials database.CGIAR collaborations result in a vast quantity of data; the International Wheat Improvement Network alone has collected over 20 million raw phenotypic data points [369].Such datasets, strengthened by the ideas on platforms and standardization described above, present enormous opportunities for crop modelling, but greater investment is required to effectively automate data extraction from different databases.The CGIAR Platform for Big Data in Agriculture was launched in 2017 with the explicit purpose of making CGIAR datasets findable, accessible, interoperable and reusable (FAIR).The platform brings together CGIAR, academia, and other actors from public and private sector and civil society to harness the capabilities of big data to enhance the impact of international agricultural research.One recent progress from the platform to help filling the data gap has been the development of the open access, searchable data harvester GARDIAN (Global Agriculture Research Data Innovation and Acceleration Network), which enables the discovery of agricultural datasets and publications across the CGIAR system and beyond.
A global "science commons" attitude should also be encouraged (for both public entities and private organizations), at least for non-propriety items like E data and generic responses of crops to E and M factors.Funding bodies could facilitate timely sharing of data by investing more explicitly in research publication, with the condition of public access [13].Such interventions could also result in more thorough reporting of experimental treatments and conditions by researchers, as well as greater availability of datasets that are normally not written up.

Conclusions
Moving forward, we propose the following developments for the crop modelling field: 1.
New insights into the biology of plant development and response to environmental cues will naturally enable model development, as it becomes possible to incorporate additional functions into the main algorithms of models.While this will contribute to a continuous improvement, it contains an inherent risk of making models so complex that they lose their initial purpose of being a simpler representation of the reality, but not the reality itself.A recent paper by Soltani and Sinclair [309] illustrates the value of keeping models simple.While continuous improvements to the science of crop models are critical, there should be a balance between the level of complexity and approximation of the processes they model.

2.
While crop models are primarily being developed and used as scientific research tools, their stakeholders extend far beyond research institutions to include farmers, policymakers, development practitioners, and the private sector.Maintaining close links with this international community is a crucial part of developing and applying models, and one where CGIAR centres can play a major role through their network of public and private sector partners at all levels.Coordinated linkages between stakeholders will help define the crop modelling community's priorities and stimulate further development in the modelling science, as well as transparently communicate the capacity and caveats of crop modelling approaches.For example, crop models are increasingly equipped to provide stochastic estimates of crop growth and yields (e.g., Yield Prophet ® ) and potential impacts of altering specific genetic or agronomic traits (e.g., [370]).Meanwhile, foresight studies increasingly incorporate crop models and couple them with large-scale economic models to address questions from policymakers on the potential impacts of global change scenarios.

3.
Scientists, including many working in CGIAR Centres, are progressively picking robust and complex methods from advanced mathematics, computer science, and physics principles to address the challenges outlined in this review.Artificial neural networks and cellular automata, coupled with fuzzy logic percolation and individual/agent-based approaches, are among the modelling methods being applied to improving crop management (Figure 2).Artificial neural networks are a system inspired by the functioning of neural networks of animal brains [371].Its main advance is that it permits the use of data collected at any suitable scale, bypassing the ambiguities that can occur when fitting equations to estimated parameters [371].Meanwhile cellular automate is a spatially and temporally discrete modelling concept that describes the dynamics through interactions and synchronous evolution using elements such as grid of cells, states, neighbourhood, transition rules, and time step [372].This approach, when coupled with fuzzy logic (i.e., representation of the level of truth with scaled values from 0 and 1), yields a robust model capable of incorporating vague and imprecise knowledge to guide and improve pest management [372].4.
As great strides are made in harnessing the genotypic information of more and more traits, the potential exists to interface a layer of genetic information within crop coefficients.Crop coefficient(s), leading to a specific response, could be driven by QTL with a value of putative effect.This kind of progress would theoretically enable the prediction of ideotype performance, based on marker content [373][374][375].However, more knowledge is needed to better understand epistatic interactions and their role in determining phenotypes.

5.
Genomic selection is increasingly being used to select promising recombinants, based on their marker content.Genomic selection has shown that breeding can be done saving time and money and thus improving the genetic gains.Crop models could therefore be interfaced to deal with the portion of the phenotypic variance that is accounted for by these interactions.6.
The G × E × M × S paradigm offers scope to address the complexities of agri-food systems, using data and models, without falling into the trap of over-simplification, inherent in modelling, by harnessing data and tools across multiple scientific disciplines and thus, enhancing the potential for creating impact for the target beneficiaries of the CGIAR.
Designing intelligent agri-food systems is central to meeting global development challenges.The renewed focus on the systems-level has created significant opportunities for modellers to participate in enhancing the impact of science on development.However, a coherent approach-based on principles of transparency, cooperation, and innovation-is essential to achieve this.Modelling is the best tool we have for grasping the big picture and truly comprehending how the whole is greater than the sum of its parts.Millions of farmers, and the societies that depend on their production, are relying on us to step up to the plate.

Figure 1 .
Figure 1.Crop modelling factors and relevant interactions.Four-way modelling captures the complex interactions of genotype, environment, management, and socioeconomics.

Figure 1 .
Figure 1.Crop modelling factors and relevant interactions.Four-way modelling captures the complex interactions of genotype, environment, management, and socioeconomics.

Figure 3 .
Figure 3. Probabilistic distribution of: (a) rice equivalent yield (REY), (b) irrigation,and (c) evapotranspiration of a rice-wheat system and maize-wheat system in central Punjab, India over 40 years.An example of how model outcome can be used to guide the field experimentation, given that testing all options in field is not feasible and are too expensive.Shaded bars are 25th-75th percentiles; whisker caps are 10th and 90th percentile and black dots 5th and 95th percentile, horizontal line is median.
Figure 3.Probabilistic distribution of: (a) rice equivalent yield (REY), (b) irrigation, and (c) evapotranspiration of a rice-wheat system and maize-wheat system in central Punjab, India over 40 years.An example of how model outcome can be used to guide the field experimentation, given that testing all options in field is not feasible and are too expensive.Shaded bars are 25th-75th percentiles; whisker caps are 10th and 90th percentile and black dots 5th and 95th percentile, horizontal line is median.

Figure 4 .
Figure 4. Simulated soil water content with and without capillary upflow.Bars at the bottom indicate daily upflow from the water table.The green line represents the drained upper limit, while the red indicates the target stress level (50% of drained upper limit until anthesis and 30% thereafter).With a water table present, irrigations were required 7 and 32 days after planting, whereas without a water table, two more irrigations were required at 77 and 92 days after planting.

Figure 4 .
Figure 4. Simulated soil water content with and without capillary upflow.Bars at the bottom indicate daily upflow from the water table.The green line represents the drained upper limit, while the red indicates the target stress level (50% of drained upper limit until anthesis and 30% thereafter).With a water table present, irrigations were required 7 and 32 days after planting, whereas without a water table, two more irrigations were required at 77 and 92 days after planting.
ultimately involves G × E × M × S interactions, it highlights the interface between the biophysical and human or socioeconomic aspects of productivity.For any given E, M interventions are determined largely by S, since they must be both feasible and cost-effective for the farmer.Potential yield and water-limited yield represent G × E potential, depending on the water regime.Nutrient-limited yield and part of the environment-reduced yields (due to biotic and abiotic pressures) can often be overcome by M (through proper fertilization, integrated pest management, etc.), assuming such interventions are economically feasible and socially acceptable.

Figure 5 .
Figure 5. Graphical presentation of the various factors contributing to the yield gap. Actual yield is often limited by many interacting factors.Attempts have been made to incorporate this by introducing the concept of attainable yield, and the attainable yield gap.This is defined in terms of the best local[280] or global[281] practices, under similar biophysical and climatological circumstances.

Figure 7 .
Figure 7. Crop modelling predictions of trait alterations on yield (top chart) in different stress scenarios of a large sorghum growing area in India (bottom map).The top graph shows APSIM predictions of grain yield across five stressed environments in a standard genotype that is insensitive to vapour pressure deficit (VPD; orange bars) and in the same genotype when made sensitive to VPD (green bars).The bottom map shows the areas where introgression of the VPD-sensitivity trait would have the strongest effect (based on the frequency of different stress scenarios in the growing regions), i.e., more in the northern sub-region where the prevalence of severe stress is higher, less so in the two southern sub-regions where the prevalence of mild stress is higher.Adapted from Kholova et al.[151].

Figure 7 .
Figure 7. Crop modelling predictions of trait alterations on yield (top chart) in different stress scenarios of a large sorghum growing area in India (bottom map).The top graph shows APSIM predictions of grain yield across five stressed environments in a standard genotype that is insensitive to vapour pressure deficit (VPD; orange bars) and in the same genotype when made sensitive to VPD (green bars).The bottom map shows the areas where introgression of the VPD-sensitivity trait would have the strongest effect (based on the frequency of different stress scenarios in the growing regions), i.e., more in the northern sub-region where the prevalence of severe stress is higher, less so in the two southern sub-regions where the prevalence of mild stress is higher.Adapted from Kholova et al.[151].

Figure 8 .
Figure 8. Accuracy of univariate and multivariate genomic prediction models for yield.Comparison of ordinary univariate pedigree and genomic selection model accuracies with those of multivariate pedigree and genomic selection models incorporating either replicated or unreplicated secondary trait data in optimal, severe drought, and late heat environments.(UV = univariate; MV = multivariate).Adapted from Rutkoski et al., 2016 [115].

Figure 8 .
Figure 8. Accuracy of univariate and multivariate genomic prediction models for yield.Comparison of ordinary univariate pedigree and genomic selection model accuracies with those of multivariate pedigree and genomic selection models incorporating either replicated or unreplicated secondary trait data in optimal, severe drought, and late heat environments.(UV = univariate; MV = multivariate).Adapted from Rutkoski et al., 2016 [115].

Figure 9 .
Figure 9. Plausible scenario of a food systems shock.The Lloyd's report aids crop modellers by clarifying the likely impacts of a plausible sudden food supply shock due to simultaneous occurrence of both biotic and abiotic stress on supply chains and societies[363].Adapted from Lloyd's report with the permission from the authors.

Table 1 .
Modelling developments by CGIAR centres, 1970s to present.
1990s FURTHER WORK ON RICE, BEANS, MAIZE International Rice Research Institute (IRRI): ORYZA model family developed: IRRI, Wageningen UR & National Agricultural Research Systems of 9 Asian countries, Systems Analysis and Simulation in Rice Production (SARP) project.Yield gap analysis, identifying management constraints, climate change analysis, declining yields, Modelling interplant competition (INTERCOM models), Water management optimization, Linking crop and pest models, Plant type design and breeding strategies, Comparing yields in long term trials across years to yield potential, Optimize grain yields with minimum methane emissions, Predict impact of climate change on rice production in Asia using all

Table 3 .
[265,266] analysis between different system components under different objectives for various crop establishment methods, rice sowing dates, and variety duration[265,266].FP-Farmer's practice, CA-Conservation Agriculture, REY-Rice equivalent yield, ET-Evapotranspiration, WP ET -Water productivity based on evapotranspiration, WP I -Water productivity based on irrigation.