Land-surface models (LSMs) typically simulate a single crop per year in a field or location. However, actual cropping systems are characterized by a succession of distinct crop cycles that are sometimes interspersed with long periods of bare soil. Sequential cropping (also known as multiple or double cropping) is particularly common in tropical regions, where the crop seasons are largely dictated by the main wet season. In this paper, we implement sequential cropping in a branch of the Joint UK Land Environment Simulator (JULES) and demonstrate its use at sites in France and India. We simulate all the crops grown within a year in a field or location in a seamless way to understand how sequential cropping influences the surface fluxes of a land-surface model. We evaluate JULES with sequential cropping in Avignon, France, providing over 15 years of continuous flux observations (a point simulation). We apply JULES with sequential cropping to simulate the rice–wheat rotation in a regional 25 km resolution gridded simulation for the northern Indian states of Uttar Pradesh and Bihar and four single-grid-box simulations across these states, where each simulation is a 25 km grid box. The inclusion of a secondary crop in JULES using the sequential cropping method presented does not change the crop growth or development of the primary crop. During the secondary crop growing period, the carbon and energy fluxes for Avignon and India are modified; they are largely unchanged for the primary crop growing period. For India, the inclusion of a secondary crop using this sequential cropping method affects the available soil moisture in the top 1.0 m throughout the year, with larger fluctuations in sequential crops compared with single-crop simulations even outside the secondary crop growing period. JULES simulates sequential cropping in Avignon, the four India locations and the regional run, representing both crops within one growing season in each of the crop rotations presented. This development is a step forward in the ability of JULES to simulate crops in tropical regions where this cropping system is already prevalent. It also provides the opportunity to assess the potential for other regions to implement sequential cropping as an adaptation to climate change.
Introduction
Climate change is likely to impact all aspects of crop production, affecting plant growth, development and crop yield as well as cropping area and cropping intensity . The impact of climate change on agriculture has been the focus of several large collaborative projects such as the Agricultural Model Intercomparison and Improvement Project (AgMIP; ) and the Inter-Sectoral Impact Model Intercomparison Project ISIMIP;. These projects have highlighted the likelihood of competition between crops grown for food and those grown for bio-energy in order to mitigate climate change . discuss how the use of sequential cropping systems may have made it possible for populations in some areas to adapt to large changes in monsoon rainfall between 2200 and 2100 BC. These ancient agricultural practices are common today across most tropical countries but may also be a useful adaptation in areas where traditionally mono-crop systems are used, in order to meet a future rising demand for food or the demand for bio-fuels. This sort of adaptation is already happening in some locations. show that longer growing seasons in the extratropics have made the cultivation of multiple crops in a year at northern latitudes more viable. Warmer spring temperatures in the Brahmaputra catchment have allowed earlier planting of a winter crop, leaving time for a second crop .
The economy of South Asia is highly dependent on the agricultural industry and other industries also with a high demand for water . The most important source of water for this part of the world is the Asian summer monsoon (ASM), which typically occurs between June and September ; this phenomenon provides most of the water resource for any given year. The South Asian crop calendar is defined by the ASM, which has an important influence on the productivity across the whole year , thereby affecting crop production outside the monsoon period.
Intercropping or sequential cropping allow farmers to make the most efficient use of limited resources and space in order to maximize yield potential and lower the risk of complete crop failure. These techniques also influence ground cover, soil erosion and chemical properties, albedo and pest infestation . Intercropping is the simultaneous cultivation of multiple crop species in a single field , while sequential cropping (also called multiple or double cropping) involves growing two or more crops on the same field in a given year . We use the term sequential cropping from here on to avoid confusion with other cropping systems. Sequential cropping systems are common in Brazil where the soybean–maize or soybean–cotton rotations are used and for South Asia where the rice–wheat systems are the most extensive, dominating in many Indian states , across the Indo-Gangetic Plain (IGP) and Pakistan . States such as Punjab, Haryana, Bihar, Uttar Pradesh and Madhya Pradesh account for approximately 75 % of national food grain production for India. Rice–rice rotations are the second most prevalent crop rotation to rice–wheat rotations; these are typically found in the northeastern regions of India and Bangladesh , with some regions cultivating as many as three rice crops per year.
Modelling sequential cropping in land-surface models
The modelling of crop rotations is a regular feature of soil carbon simulations . found that the rice–wheat rotation, common across the IGP, has helped maintain carbon stocks. However, in recent years, the yields of rice and wheat have plateaued, leading farmers to diversify and include other additional crops in the rotation, potentially depleting carbon stocks. The modelling of crop rotations has also been represented in the field of agricultural economics with work regarding sequential cropping being mainly to understand influences on decision-making, therefore focusing on short timescales and at the farm management level .
Many dynamic global vegetation models (DGVMs), used to study the effects of climate change, simulate a single crop in a field per year, both for individual sites and gridded simulations. This may be due in part to some global observation datasets such as reporting only one growing period per year at a given location for most crops . Where different crop calendars are available for different regions, e.g. MIRCA2000 , rice and wheat are divided equally between the kharif (i.e. sown during the monsoon and harvested during the autumn) and rabi seasons (i.e. the drier winter/spring growing season), when in reality wheat is only grown during the rabi season .
The Lund–Potsdam–Jena managed land model (LPJml; ) is one of the few models that is able to simulate sequential cropping. extend LPJml to consider sequential cropping in Africa for two different crops on the same field within a year. specify different growing periods for each crop in the rotation, where the growing period is calculated from the sum of the daily temperatures above a crop-specific temperature threshold. use the method to specify the onset of the main rainy season as the start of the growing season, where growing season is defined as the period of time in which temperature and moisture conditions are suitable for crop growth. The growing period of the first crop in the rotation begins on the first wet day of the growing season, with the second crop assumed to start immediately after harvest of the first crop. find that when considering the impact of climate change, the type of cropping system is important because yields differ between crops and cropping systems. also use a version of LPJml refined for South Asia to estimate water demand and crop production for South Asia. combine the output from two separate simulations, each with different kharif and rabi land-use maps and zonal sowing and harvest dates based on observed monsoon patterns. find that accounting for multiple different crops being grown on the same area at different times of the year improves the simulations of demand for water for irrigation, particularly the timing of the demand. and simulate more than one crop growing on the same area using very different methods; both have highlighted the importance of representing this type of cropping system.
demonstrate that the Interactions between the Soil, Biosphere and Atmosphere (ISBA) land-surface model (LSM) , specifically ISBA-A-gs , is able to represent a 12-year succession of arable Mediterranean crops for a site in Avignon, France . This type of cropping system is not typically represented in LSMs; however, this study showed that the implementation of crop successions in an LSM leads to a more accurate representation of cumulative evapotranspiration over the 12-year period. It would be beneficial for more land-surface models to develop the capability to simulate different cropping systems and link crop production with irrigation both to improve the representation of the land surface in coupled models and to improve climate impacts assessments.
In this paper, we describe and implement sequential cropping in the Joint UK Land Environment Simulator (JULES). We simulate the same Avignon site described in and two states in India to illustrate and evaluate the method implemented in the JULES stand-alone model at version 5.2 for simulating sequential-crop rotations. By using Avignon and India, we simulate two types of crop rotation. We define the Avignon crop rotation as an irregular crop rotation due to the occurrence of long fallow periods, with extended periods of bare or almost bare soil. There are no long fallow periods in the rice–wheat rotation for the northern Indian states of Uttar Pradesh and Bihar, so we refer to this as a regular crop rotation. The rice–wheat rotation is the dominant cropping system for Uttar Pradesh and Bihar, with these states being key producers of these crops. Sequential cropping is used in a series of single-grid-box simulations (where each grid box is 25 km) across these two states for comparison with a point simulation in Avignon. We also use sequential cropping in a regional 25 km resolution gridded simulation to demonstrate that this method can be applied at larger scales. In the India simulations (both the regional and single grid box), the model uses standard soil parameters, with meteorological information provided from a regional climate simulation (Sect. ). However, for the point simulation in Avignon, we use the site characteristics and meteorological information from the site itself (Sect. ). This paper is structured as follows: Sect. describes the JULES-crop model, the rationale and the method for implementing sequential cropping in JULES. The simulations are described in Sect. . In Sect. , we present our hypotheses for assessing the impact of sequential cropping in JULES, discuss the climate of the simulated regions and describe the observations we use to evaluate the model. The results from the evaluation of the simulations and the assessment of the impact of sequential cropping are provided in Sect. and discussed in Sect. , with conclusions in Sect. .
JULES flags used that are new or different from those in .
FlagJULESAvignonIndiaEffectnotationsettingssettingsof switchCanopy radiationschemecan_rad_mod66Selects the canopy radiation scheme.Irrigation demandl_irrig_dmdFTSwitches on irrigation demand.Irrigation schemeirr_crop–2Irrigation occurs when the DVI of the crop is greater than 0.Physiologyl_trait_physFFSwitches on trait-based physiology when true.Sowingl_prescsowTTSelects prescribed sowing.Plant maintenance respirationl_scale_resp_pmFFSwitch to scale respiration by water stress factor. If false, this is leaf respiration only, but if true, it includes all plant maintenance respiration.Crop rotationl_croprotateTTA new switch to use the sequential cropping capability.Irrigation on tilesfrac_irrig_all_tiles–FSwitch to allow irrigation on all or specific tiles.Irrigation on specific tilesset_irrfrac_on_irrtiles–TA new switch to set irrigation to only occur on a specific tile.Specify irrigated tile(s)irrigtiles–6Setting to set the value(s) of the specific tile(s) to be irrigated.Number of tiles irrigatednirrtile–1Setting to set how many tile(s) are to be irrigated.Set a constant irrigation fractionconst_irrfrac_irrtiles–1.0A new setting to set the value(s) of the irrigation fraction for specific tile(s) to be irrigated in the absence of a file of irrigation fractions.
JULES plant functional type (PFT) parameters and values modified for use in this study. We include only the values that have been changed or are new in JULES since .
ParameterJULES notationDescription (units)Winter wheatSorghumSpring wheatRiceTlowt_low_ioLower temperature for photosynthesis (∘C).518515Tuppt_upp_ioUpper temperature for photosynthesis (∘C).30533040neffneff_ioScale factor relating Vcmax with leaf nitrogen concentration.0.8e-30.75e-30.8e-30.95e-3nl(0)nl0_ioTop leaf nitrogen concentration (kgN(kgC)-1.0.0730.070.0730.073fsmc methodfsmc_mod_ioWhen equal to 0, we assume an exponential root distribution with depth.00When equal to 1, the soil moisture availability factor, fsmc, is calculated using average properties for the root zone.11drrootd_ft_ioIf fsmc_mod_io = 0, dr is the e-folding depth (m).0.50.5If fsmc_mod_io = 1, dr is the total depth of the root zone (m).1.51.5p0fsmc_p0_ioParameter governing the threshold at which the plant starts to experience water stress due to lack of water in the soil.0.50.50.50.5μrlnr_nl_ioRatio of root nitrogen concentration to leaf nitrogen concentration.0.390.390.390.39μslns_nl_ioRatio of stem nitrogen concentration to leaf nitrogen concentration.0.430.430.430.43Q10,leafq10_leaf_ioQ10 factor in the Vcmax calculation.1.01.01.01.0
JULES crop parameters used in this study. The sorghum cardinal temperatures are from with the other parameters being those used for maize in . We include only the values that have been changed or added since . Table 3 of provides the original PFT parameters and Table 4 of provides the original crop parameters.
ParameterJULESDescription (units)WinterSorghumSpringRicenotationwheatwheatTbt_bse_ioBase temperature (K).273.15284.15273.15278.15Tmt_max_ioMax temperature (K).303.15317.15308.15315.15Tot_opt_ioOptimum temperature (K).293.15305.15293.15303.15TTemrtt_emr_ioThermal time between sowing and emergence (degree days).35803560TTvegtt_veg_ioThermal time between emergence and flowering (degree days).Table Table Table Table TTreptt_rep_ioThermal time between flowering and maturity (degree days).Table Table Table Table Tmortt_mort_ioSoil temperature (second level) at which to kill crop if DVI>1 (K).273.15281.15273.15281.15fyieldyield_frac_ioFraction of the harvest carbon pool converted to yield carbon.1.01.01.01.0DVIinitinitial_c_dvi_ioDVI at which the crop carbon is set to Cinit.0.00.00.00.0DVIsensen_dvi_ioDVI at which leaf senescence begins.1.51.51.51.5Cinitinitial_carbon_ioCarbon in crop at emergence in kgCm-2.0.010.010.010.01The JULES-crop modelModel description
The JULES model is the land-surface scheme used by the UK Met Office for both weather and climate applications. It is also a community model and can be used in stand-alone mode, which is how it is used in the work presented here. The parameterization of crops in JULES (JULES-crop) is described in and . JULES-crop is a dual-purpose crop model intended for use within stand-alone JULES, enabling a focus on food production and water availability applications, and as the land-surface scheme within climate and Earth system models. JULES-crop has been used in stand-alone mode in recent studies such as and . The aim is that these studies and this one will lead to using JULES in these larger models to allow the feedbacks from regions with extensive croplands and irrigation systems, like South Asia, to have an effect on the atmosphere, e.g. via methane emissions from rice paddies or evaporation from irrigated fields .
JULES is a process-based model that simulates the fluxes of carbon, water, energy and momentum between the land surface and the atmosphere. JULES represents both vegetation (including natural vegetation and crops) and non-vegetation surface types including; urban areas, bare soil, lakes and ice. With the exception of the ice tile, all these tiles can co-exist within a grid box so that a fraction of the surface within each grid box is allocated between surface types. For the ice tile, a grid box must be either completely covered in ice or not . JULES treats each vegetation type as a separate tile within a grid box, with each one represented individually with its own set of parameters and properties, such that each tile has a separate energy balance. The model and the equations it is based on are described in detail in and . Prognostics such as leaf area index (LAI) and canopy height are therefore available for each tile. The forcing air temperature, humidity and wind speed are prescribed for the grid box as a whole for a given height. Below the surface, the soil type is also uniform across each grid box (where the number of soil tiles is set to 1). We use JULES-crop to simulate the crops in this study. The main aim of JULES-crop is to improve the simulation of land–atmosphere interactions where crops are a major feature of the land surface .
Photosynthesis in JULES-crop uses the same parameters and code as the natural plant functional types (PFTs). There are two temperature parameters: Tlow and Tupp; these define the upper and lower temperature parameters for leaf biochemistry and photosynthesis within JULES and are used to calculate the maximum rate of carboxylation of RuBisCO (unstressed by water availability and ozone effects – Vcmax, with units of molCO2m-2s-1), as defined in and reproduced here in Eq. (). Equation () demonstrates the Vcmax at any desired temperature.
1Vcmax=neffnl(0)fT(Tc)1+e0.3(Tc-Tupp)1+e0.3(Tlow-Tc)2fT(Tc)=Q10leaf0.1(Tc-25),
where fT is the standard Q10 temperature dependence (given in Eq. ) and Tc is the canopy temperature. neff represents the scale factor in the Vcmax calculation (in units of molCO2m-2s-1kgC(kgN)-1) and nl(0) the top leaf nitrogen concentration (in units of kgN(kgC)-1). More details regarding the calculation of Vcmax are provided in and . Vcmax is an important component in two limiting factors for photosynthesis: the RuBisCO-limited rate and the rate of transport of photosynthetic products; Eq. () shows the relationship between Vcmax and temperature. Gross primary productivity (GPP) is used to describe the total productivity of a plant; this defines the gross carbon assimilation in a given time. Net primary productivity (NPP) is GPP minus plant respiration; NPP is used in the crop partitioning code and subsequently in the calculation of the yield in JULES. The nitrogen cycle in JULES cannot yet be used with the crop model, so in this study the same assumption is made as in : that crops are not nitrogen limited.
The effective temperature (Eq. ) is the function that the model uses to relate air or leaf temperature to the cardinal temperatures that define a plant's development; these are the base temperature (Tb), maximum temperature (Tm) and optimum temperature (To) and are specific for each crop. Different models define their effective temperature function in different ways; for example, Fig. 1 of provides a number of different possible definitions. The JULES definition described by Eq. () is most similar to type 4 given in . Type 4 increases gradually towards the optimum temperature with a steeper decline from the optimum to the maximum. Other functions have no decline or a flatter top, which can have different effects on the development of the crop. In JULES, the cardinal temperatures and the 1.5 m tile (i.e. air) temperature (T) are used to calculate the thermal time, i.e. the accumulated effective temperature (Teff) to which a crop is exposed . Table summarizes the settings for these temperatures used in this analysis. The crop model integrates an effective temperature over time as the crop develops through these stages, with the carbon partitioned according to the development index (DVI).
Teff=0forT<TbT-TbforTb≤T≤To(To-Tb)(1-T-ToTm-To)forTo<T<Tm0forT≥Tm
The DVI is a function of the thermal time since emergence; therefore, DVI of -1 is sowing, 0 is emergence, and 1 is flowering. Maturity and therefore harvest occur at a DVI of 2 under standard growth conditions but may be harvested earlier in other situations in the model . In reality, the maturity date and the harvest dates are not usually the same date. The integrated effective temperature in each development stage is referred to as the thermal time of that development stage (Eq. and ).
Crop development can also be affected by the length of the day. However, in these simulations, as in , this effect is not included. The thermal time is then used to calculate the rate of crop development or rate of increase of the development index, described by Eq. ():
dDVIdt=TeffTTemrfor-1≤DVI<0(TeffTTveg)for0≤DVI<1TeffTTrepfor1≤DVI<2,
where TTemr is the thermal time between sowing and emergence, TTveg and TTrep are the thermal time between emergence and flowering and between flowering and maturity, respectively. These are calculated either using a temperature climatology from the driving data and sowing dates from observations or using the method presented in to create a reliable sowing and harvest dataset. The advantage of using the method is that there are no missing data, which is often the case when using observed data. Whichever source of sowing and harvest dates are used, the aim is for the crop to reach maturity, on average by the harvest date. The sowing and harvest dates used in the simulations in this analysis are described in Sect. .
Implementing sequential cropping in JULES-cropJULES-crop: rationale for implementing sequential cropping in JULES-crop
JULES-crop is typically run as a single-crop model, represented by the red curve in Fig. , where a primary crop is simulated but no second crop is possible, and the land is left fallow with a minimum surface cover. In many regions, sequential cropping is the main cropping system used, with several crops cultivated one after another. JULES-crop has been developed for implementation in Earth system and climate models for application in adaptation and mitigation studies. Only being able to simulate one crop per year is therefore limiting application in many parts of the world. In the changes to JULES described in Sect. , new controls are implemented to allow the current JULES-crop code to be run more than once in a year at a particular location, so that sequential cropping systems can be represented in JULES. Sequential cropping is available from version 5.7 of JULES; this option is represented by the black curve in Fig. .
The implementation of sequential crops in JULES is part of a project to develop simulations for South Asia to understand the impacts of climate change on both agriculture and water sectors using existing regional climate model (RCM) projections . This will improve understanding of the impacts of climate change and how they affect each other. Sequential cropping provides clear added benefits for the following reasons:
providing a more realistic representation of the observed surface land cover;
allowing the continuous simulation of a location where different crops are grown within the same area, thereby simulating water resource demand from crops;
allowing the climate to affect both the water and crops, while simultaneously allowing interactions between water and crops throughout the year makes it possible to simulate the integrated impacts of climate change on these two sectors; and
providing the opportunity to investigate the impact of adopting sequential cropping for regions where it is not currently used.
The sequential cropping method and the modifications made to JULES-crop
The sequential cropping method implemented into JULES as part of this study is illustrated by the flowchart in Fig. and described here using the Avignon site simulation. The Avignon site is a point run which is assumed to be entirely used to grow sorghum (from spring to late summer) and winter wheat (from winter to early summer). JULES updates the fraction of the site that is allocated to sorghum (winter wheat) just before the sowing date so that the appropriate crop occupies the whole of the site. The fraction of the site that is sorghum (winter wheat) is prescribed in the Avignon case using observed sowing and harvest dates. Once the fraction is updated, the crop is sown; it then develops between the stages of sowing and emergence, emergence and flowering, and flowering and maturity.
In order to simulate the characteristics of a typical sequential cropping location using JULES, we made modifications to both JULES-crop and the irrigation code. To simulate crops in sequence on the same grid box, each crop must be completed cleanly so the second one can be sown accordingly. The specification of a latest harvest date (latestharvestdate) forces the harvest of the first crop regardless of whether it has reached maturity or not. The latestharvestdate is a safeguard built into the model, usually set to a date well after the expected harvest date. If the latestharvestdate is used, an alert is triggered, which provides some initial information to aid the investigation of the problem. In this study, the latestharvestdate is set but never actually required for any of the simulations, which is the ideal scenario. The latestharvestdate safeguard is preferable to the simulation of a crop growing for an unrealistically long time, i.e. developing too slowly and overlapping the next growing season. This is essential for the implementation of sequential cropping at a global or regional scale, where the model is forced to grow crops that are potentially unsuitable for a particular grid box. This is more likely for global simulations, which typically simulate a restricted set of crop types and varieties. These modifications are controlled using the l_croprotate switch (Table ). Therefore, l_croprotate ensures the following:
All crops are initialized at the start of a simulation so that they can be used later when they are needed within the crop rotation being modelled.
If JULES is simulating a crop rotation, the user must supply a latestharvestdate so that the first crop is harvested before the second crop is sown (a latestharvestdate can also be specified without using l_croprotate).
The current JULES default for irrigation allows individual tiles to be specified (when frac_irrig_all_tiles is set to false) but the irrigation is applied as an average across a grid box and therefore actually occurs across tiles. The flag set_irrfrac_on_irrtiles restricts the irrigation to the tiles specified by irrigtiles only (Table ). This new functionality is needed because many locations that include crop rotations include crops that both do and do not require irrigation.
The flowchart shown in Fig. is equally applicable to the India simulations. Rice is therefore represented by the summer crop (green boxes) and wheat is represented by the winter crop (purple boxes). This method could be extended to include as many crops as those occurring in a rotation at a particular location.
Model simulations
The description of the simulations is divided into two sections. Section describes the Avignon point simulations; this is a well-observed site used to describe and demonstrate the sequential cropping method and evaluate it against observations at this location. Avignon is a typical Mediterranean crop succession characterized by a succession of winter and summer crops and in between a period of bare soil. When a summer crop follows a winter crop, the period of bare soil can last up to 9 months. Winter crops are generally seeded October–November with harvest towards the end of June–July. Summer crops are seeded in late April–May and harvested at the end of August. Section describes the simulations of northern India where a more traditional sequential cropping system is commonly used, with a regular rotation between rice during the wetter kharif season and wheat during the drier rabi season. The parameter settings and switches used in JULES for the simulations in this study are provided in Tables , and . The Avignon and India simulations use the same settings wherever possible; these are provided in Table (see Avignon settings and India settings columns).
The PFT parameter settings are also broadly the same between simulations, with the majority of these from and therefore based on natural grasses. The crops are different between the two sets of simulations with winter wheat and sorghum at the Avignon site and spring wheat and rice at the India locations. The PFT parameters used in this study that govern Vcmax: including the lower (Tlow) and upper (Tupp) temperatures for photosynthesis, neff and nl(0) are tuned to the maximum leaf assimilation expression from for each crop (Table ). These values are consistent with the wider literature . The parameters μrl and μsl are the ratios of root-to-leaf and stem-to-leaf nitrogen concentrations, respectively; these are tuned to those given in to lower the plant maintenance respiration, which was high in some of the initial simulations. The crop parameters are mainly from , with maize parameters used for sorghum (Sect. ), except for the cardinal temperatures (Table ) which are from .
The calculation of the soil moisture availability factor (beta, β, Table ) is different between the Avignon and India simulations. β in each layer is 0 below the wilting soil moisture and 1 above a threshold; this is shown in Fig. 1 of . In the Avignon simulations, we assume a rectangular root distribution and the total depth of the root-zone dr to be 1.5 m, equivalent to the observed average maximum root depth over all of the years at the Avignon site. β is then calculated using this maximum root depth together with the average properties of the soil. The India single-grid-box simulations assume an exponential root distribution with an e-folding depth dr of 0.5 m because we do not have an observed root depth for these locations. In all simulations in this study, we adjust the parameters that affect the use of water by the plant so that the plants experience less water stress (this parameter is P0 and is set to 0.5; , Table ). This is because water stress is not the main focus of this analysis, but the representation of soil moisture stress on vegetation is a known issue in JULES; this is the subject of a large international collaborative effort . The individual simulations are described in more detail in Sect. and for the Avignon and India simulations, respectively.
The purpose of including Avignon is because it provides a wealth of observations for evaluating land-surface models, where there is no equivalent site for South Asia. Observations of these fluxes show if the model is correctly representing the fluxes and coverage of the land surface. The purpose of including a simulation that does not use the crop model but approximates crops using grasses is to show how the model performs with the correct LAI and height; i.e. it is a clean test of the representation of leaf photosynthesis, stomatal conductance, water stress and leaf-to-canopy scaling within the model (these parts of the code are shared by both natural vegetation and crops).
Avignon, France, simulations
The Avignon “remote sensing and flux site” of the National Research Institute for Agriculture, Food and Environment (INRAE), described in , provides a well-studied location (Avignon, France 43.917∘ N 4.878∘ E) with several years of crop rotation data. The Avignon simulations focus on the period between 2005 and 2013 with a rotation of just two crops: winter wheat and sorghum. JULES already contains parameterizations for wheat and maize. The wheat in JULES is the spring variety which is similar to the winter wheat crop that is grown in Avignon. Spring wheat does not require a vernalization period, which is a process usually needed for winter wheat varieties to achieve optimum yields . Vernalization is not explicitly implemented in JULES; therefore, spring and winter wheat can be simulated interchangeably. The maize crop is a C4 crop that is similar to sorghum. Therefore, we use these existing parameterizations rather than develop new ones. During this period, two varieties of sorghum were cultivated. In 2009, a fodder crop variety was grown, with a shorter growing season and a larger LAI than the variety for the other two years (2007 and 2011). Therefore, the 2009 sorghum crop is planted much later in the year compared to the other two sorghum seasons (2007 and 2011) but harvested at a similar time. The aim of simulating the crops at this site is to demonstrate the new sequential cropping functionality in JULES and show how the implementation of sequential cropping affects the JULES crops simulated.
The length and detail of the observation record at the Avignon site indicate that it is an ideal site to demonstrate the method being implemented in JULES for simulating sequential cropping. High-resolution meteorological data, important for the practicalities of running the JULES model, are used to run the model using a half-hourly time step; this includes air temperature, humidity, wind speed, rainfall, radiation measurements and atmospheric pressure at a height of 2 m above the surface . Irrigation in Avignon is only applied to the summer crops, i.e. sorghum. The observed irrigation amounts are added to the precipitation driving data at the exact day and time they were applied to the crops . The irrigation and other settings governing irrigation are therefore not switched on in JULES for the Avignon site simulations (Table , column “Avignon settings”). We include simulations for the Avignon site where the crops are represented by grasses (Avi-grass) for comparison with the simulations that use the JULES-crop model. In the Avi-grass simulations, the LAI and the canopy height are prescribed from observations in order to capture the growing seasons correctly without the crop model, and the PFT parameters are adjusted to be the same as the crops. These Avi-grass simulations use the same photosynthesis and respiration calculation as JULES-crop, but this is not allowed to influence LAI as they do in the crop model. This allows the evaluation of the photosynthesis and respiration parts of the model, together with the water and energy fluxes, when the observed LAI and canopy height are used. In the Avi-grass simulations, JULES is not modelling the crops as grasses but fixing some parts of the crops (LAI and canopy height) straight to observations. We also run two simulations that use the crop model: a single-crop (Avi-single) and a sequential-crop simulation (Avi-sequential). In both the Avi-single and Avi-sequential simulations, the LAI and the canopy height are calculated by the model. The JULES total aboveground biomass is calculated from the sum of the stem, leaf and harvest carbon pools for each crop. Observed sowing and harvest dates from are used to calculate the thermal time requirements for each crop represented in the simulations; these are provided in Table . During the periods between each crop, the ground is mostly bare . The only difference between the Avi-sequential and Avi-single simulations is that Avi-single only simulates wheat; therefore, no sowing dates are provided for sorghum.
A schematic describing the single-crop model (red curve) that is part of the standard JULES-crop and the new option for including sequential cropping (black curve). This schematic represents a generic crop at a single location.
A flowchart showing the sequence followed to carry out the crop rotation in JULES. The first step (top green box) in the sequence is to update the first crop fraction; this occurs as or just before the first crop is sown.
A map showing the location of the single-grid-box simulations in the wider context of India on a map of the surface altitude (a) from the regional climate model that is used in the JULES simulations. The same locations are shown in three smaller maps (b–d) that zoom in on the states of Uttar Pradesh and Bihar. The map in panel (b) shows the total monsoon precipitation, (c) shows the minimum temperature, and (d) shows the maximum temperature averaged for the period 1991–2007.
Thermal times in degree days used in this study for the Avignon site; these are based on the observed sowing and harvest dates from .
YearCropSowing dateHarvest dateEmergence – floweringFlowering – maturitySowing DOY2005Winter wheat27 Oct 20051301.3867.5300200627 Jun 20062007Sorghum10 May 200716 Oct 2007647.6791.51302007Winter wheat13 Nov 20071401.0934.031720081 Jul 20082009Sorghum25 Jun 200922 Sep 2009462.5565.31762009Winter wheat19 Nov 20091308.6872.4323201013 Jul 20102011Sorghum22 Apr 201122 Sep 2011679.5830.51122011Winter wheat19 Oct 20111559.61039.7292201225 Jun 2012Uttar Pradesh and Bihar, India, simulations
The India simulations focus on the northern Indian states of Uttar Pradesh and Bihar. These states are key producers of rice and wheat in South Asia and use a regular rice–wheat rotation that is prevalent in this part of India . We include single-grid-box simulations and a regional simulation. The single-grid-box simulations are a selection of four locations from across these two states. For each grid box, both single-crop (referred to as India-single) and sequential-crop (referred to as India-sequential) simulations are run. India-single and India-sequential are set up in the same way, with the only difference being that sowing dates are provided for just one crop. For consistency with the rest of the simulations, only wheat is simulated in India-single. The single-grid-box simulations enable a similar analysis to that described for Avignon (Sect. ), while the regional simulation (this is only a sequential-crop run) is a demonstration of the sequential cropping method being used at larger scales. For the regional simulation, we assume that wheat and rice are grown in every grid box across the two states and the crops are not limited by nutrient availability. The sequential cropping system in this region involves growing rice during the wet monsoon months and an irrigated wheat crop during the dry winter. In these simulations (both single grid box and regional), wheat is only irrigated during its growing period and without applying limits due to water availability (this is referred to as unlimited irrigation). The wheat variety grown in India is spring wheat, which is the standard variety represented by JULES (Sect. ).
The locations of the selected grid boxes are shown on a map of the surface altitude for South Asia in Fig. a. The driving data used for these four simulations are from an RCM simulation run for South Asia for the period 1991–2007 as described below. Figure b, c and d show a close-up view of the locations selected. The map in Fig. b shows the average total monsoon precipitation for the 1991–2007 period, while Fig. c and d show the average minimum and maximum temperatures, respectively, to illustrate that these four grid boxes are representative of the climate of the wider Uttar Pradesh/Bihar region.
In both the single-grid-box and regional India simulations, JULES is run using a 3-hourly time step using driving data from ERA-Interim downscaled to 25 km using the HadRM3 RCM . This RCM simulation is one of an ensemble of simulations produced for the EU-HighNoon FP7 project for the entire Indian subcontinent (25–32∘ N, 79–88∘ E).
The HighNoon simulations are described in detail in previous publications such as and . HadRM3 provides more regional details in addition to the global data with updated 3-hourly lateral atmospheric boundary conditions interpolated to a 150 s time step.
These simulations include a detailed representation of the land surface in the form of version 2.2 of the Met Office Surface Exchange Scheme MOSESv2.2;. JULES has been developed from the MOSESv2.2 land-surface scheme, and therefore the treatment of different surface types is consistent between the RCM and JULES . In the India single-grid-box simulations, sowing dates are prescribed using climatologies calculated from the observed dataset from the government of India, Ministry of Agriculture and Farmers Welfare. Thermal times are calculated using these climatological sowing and harvest dates from and a thermal climatology from the model simulation as described in ; the values used in the simulations here are provided in Table . In the regional simulation, the thermal time requirements are estimated from the sowing and harvest dates provided by the method to avoid problems with missing observed data. The settings used for the India simulations are provided in Table (column “India settings”). Plots of the regional ancillaries for each of rice and wheat are provided in Appendix C.
The sowing day of year (sowing DOY) and thermal times in degree days used in this study for the locations in Uttar Pradesh and Bihar, India (Fig. for a map of the locations), the values given here are based on the observed sowing and harvest dates from .
Table of statistics comparing the Avignon simulations with observations for each type of run: Avi-single (single), Avi-sequential (sequential) and Avi-grass (without the crop model).
The objectives of this study are to evaluate the model against observations where they are available, for the Avignon point and the India single-grid-box simulations, and to test the following hypotheses with regard to the implementation of the presented sequential cropping method in JULES:
Null hypothesis: the inclusion of a secondary crop on the same field does not change the growth and development of the primary crop in an irregular sequential cropping rotation with long fallow periods. Alternative hypothesis: the inclusion of a secondary crop on the same field modifies the growth and development of the primary crop in an irregular sequential cropping rotation with long fallow periods.
Null hypothesis: the inclusion of a secondary crop on the same field does not change the energy and carbon fluxes in an irregular sequential cropping rotation with long fallow periods. Alternative hypothesis: the inclusion of a secondary crop on the same field modifies the energy and carbon fluxes in an irregular sequential cropping rotation with long fallow periods.
Null hypothesis: in a regular rotation without long fallow periods, the inclusion of a secondary crop on the same field does not change the crop development of the primary crop or the grid box energy and carbon fluxes and soil conditions. Alternative hypothesis: in a regular rotation without long fallow periods, the inclusion of a secondary crop on the same field modifies the crop development of the primary crop, the grid box energy and carbon fluxes and soil conditions.
The first set of hypotheses will be assessed by comparison of the observed LAI, canopy height and total aboveground biomass in Avignon with single-crop and sequential-crop simulations. The second set of hypotheses will be assessed by comparison of observed fluxes, GPP, latent heat (LE) and sensible heat (H) in Avignon with single-crop and sequential-crop simulations. For the third set of hypotheses, we compare JULES yields with observed yields for single- and sequential-crop simulations and analyse the same variables as for the first and second sets of hypotheses for four locations across the northern Indian states of Uttar Pradesh and Bihar. We also assess if the implementation of sequential crops affects the soil moisture in a regular sequential-crop system, which does not have long periods of bare soil. For a regular sequential-crop system without long fallow periods, changes in soil moisture are more likely to be due to the effects of sequential cropping and are less likely to be affected by evaporation from bare soil.
More information is provided about the Avignon site in Sect. and the climate across the four India grid box locations in Sect. , including the observations used to evaluate the simulations at these locations.
Avignon
Avignon is characterized by a Mediterranean climate with a mean annual temperature of 287.15 K (14 ∘C) and most rainfall falling in autumn (with an annual average of 687 mm). The Avignon time series of temperature (with a 10 d smoothing applied) is provided in Appendix A in Fig. a and precipitation (10 d totals, which include actual irrigation amounts) in Fig. b . In general, Avignon experiences a fairly regular distribution of rainfall throughout the year and the annual temperature range for Avignon (26 ∘C) is relatively consistent, with only a brief cold snap in early 2012 having a much lower minimum. The Avignon site represents the irregular cropping rotation, chosen because it has been observed and documented over several years (2001 to 2014), growing a range of crops throughout this period. No equivalent site to Avignon has been found for South Asia.
The observations for evaluating the model include canopy height (measured every 10 d), aboveground dry weight biomass (taken at four field locations) and LAI; biomass and LAI are destructive measurements repeated up to six times per crop cycle . H and LE flux measurements are available for several years, enabling the evaluation of the JULES fluxes. Cumulative evapotranspiration (ET) are derived from the half-hourly LE measurements. More information is documented in regarding the site and the observations available. These continuous measurements of surface fluxes provided by the Avignon dataset are a unique resource for evaluating LSMs and for testing and implementing more irregular crop rotations in LSMs.
India
The four India locations (grid boxes) selected for analysis in this study are shown on a map of South Asia in Fig. a with smaller inset plots (Fig. b, c and d) focusing on the sequential cropping region being considered across the states of Uttar Pradesh and Bihar. The climate of this region is summarized in Appendix B (Fig. ), which shows the time series of the average precipitation (Fig. a), temperatures (Fig. b) and vapour pressure deficit (VPD) (Fig. c) at each of these four grid boxes. The different crop seasons are emphasized by the different colour shading, with yellow for wheat and pink for rice. The temperatures (Fig. b) rarely reach the lower cardinal temperatures set in the model (Tb) shown for rice (green) or wheat (orange); however, the high temperatures do exceed the maximum cardinal temperatures (Tm) for these crops, especially those set for wheat. In general, EastBi is cooler than the other locations in more of the years, with the two locations in Uttar Pradesh often being the warmest. The precipitation at each location is variable (Fig. a), with variation in the distribution of precipitation through the monsoon period, which could be important for crop yields. , for example, found that in two seasons with similar rainfall totals, the distribution of the rainfall during the growing season strongly affected groundnut crop yield. There is also a clear seasonal cycle in the VPD, increasing toward the end of the wheat season and decreasing into the rice season. EastBi generally has the lowest VPD, with WestUP and EastUP usually the highest throughout the time series shown (Fig. ). These time series show that there is a gradual change in conditions from west to east across Uttar Pradesh and Bihar with increasing humidity and rainfall and decreasing maximum temperatures from west to east.
District-level area and production data from the International Crops Research Institute for the Semi-Arid Tropics are used to calculate district level yields. These are then gridded at the resolution of the ERA-Interim data (0.25∘) to ensure that the scale of simulated and observed yields matched. We also show average crop yield observations for three 5-year periods between 1993 and 2007 (1993–1997, 1997–2003, 2003–2007). Data from are made available via . are based on previous publications . All the observations used include the period of the single-grid-box simulations, which are from the period of 1991–2007. We show both of these datasets to highlight that there is a range in the estimates of yield for this region.
ResultsEvaluation against observationsAvignon
Figure shows the time series of total aboveground biomass (a), LAI (b) and canopy height (c) for Avi-sequential and Avi-single compared with observations. The crops in JULES (both single and sequential) are developing throughout the crop seasons, with maxima occurring at approximately the correct time for the crops being simulated. Therefore, the lack of vernalization in the model does not affect the simulation of winter wheat in Avignon. Avi-grass simulations are not shown, as these follow the observed canopy height and LAI exactly as these values are prescribed in the simulations without crops. The total aboveground biomass from JULES is plotted as a time series (dashed lines) for comparison with observations, which are provided as a single time series with the crop type confirmed from the timing of the observations (purple asterisks in Fig. a). The increase in biomass for both crops through the start of the season follows the observations quite closely but in most years, especially for wheat, JULES-crop (using either single or sequential crops) does not accumulate enough biomass later in the crop season to reach the observed maxima. The wheat canopy height is very close to observations for all four seasons; however, the wheat LAI is overestimated. The two wheat seasons of 2006 and 2010 are closer to the LAI observations than 2008 and 2012, but the underestimation of the biomass is greater for these seasons. highlight that 2006 and 2008 have atypical rainfall during the wheat season, with 2006 being very dry (256 mm of rain during the wheat season) and 2008 being very wet (500 mm during the wheat season). Therefore, in 2008, Avignon received 73 % of its annual average (Sect. ) during the wheat season alone; these differing conditions could explain the large differences in observed LAI and biomass between the two years .
The shorter observed growing season for the taller fodder variety of sorghum grown in 2009 can be compared with the longer season, smaller variety grown in 2007 and 2011 (solid red line Fig. b and c). No sorghum is simulated in Avi-single, so the dotted red line is equal to zero and not visible in Fig. . Avi-sequential (dashed red) simulates the 2009 sorghum season well, in terms of biomass (Fig. a), LAI (Fig. b) and canopy height (Fig. c); with relatively small differences between the simulations and observations maximum values (LAI of 1 m2m-2 and canopy height of 0.1 m). In the 2007 sorghum season, Avi-sequential (dashed red) overestimates the maximum LAI and canopy height by approximately 2 times the observations (Fig. b and c) and underestimates the total biomass (Fig. a) by about 30 %. For the 2011 season, the Avi-sequential sorghum biomass is equal to the magnitude of the observations; however, the maximum LAI is overestimated by a factor of 4 in the model (similar to 2007) and the maximum canopy height is approximately 2 times the observed maximum.
The peaks in productivity shown in the LAI in Fig. b are consistent with the two years (2006 and 2007) of GPP observations (black line Fig. a). The 2006 wheat crop is represented in the GPP of all three simulations, although it is underestimated in all of them (Fig. a). The GPP in Avi-single is lower than both Avi-sequential and Avi-grass during the second half of the wheat growing period. The decline in GPP at the end of the 2006 wheat season is quite close to the observations for the three simulations, with Avi-grass (red line) being slightly early and both the crop simulations; Avi-sequential (blue line) and Avi-single (cyan line) being slightly late. For the sorghum growing period, the magnitude and timing of the maximum GPP for Avi-sequential (blue line) are a good fit to observations. However, the increase in GPP begins slightly too early for Avi-sequential and slightly late for Avi-grass. The Avi-grass simulations slightly underestimate the maximum GPP during the sorghum season and it occurs a little later than observed (Fig. a). The decline in GPP at the end of the sorghum season occurs at the same time as the observations for both Avi-grass and Avi-sequential. These results are quantified in Fig. , with both Avi-grass (Fig. a) and Avi-sequential (Fig. b) showing a strong linear correlation, with r values greater than 0.8. The values for Avi-single (Fig. c) are lower, with an r value of approximately 0.5; this is because the observations contain seasons with both sorghum and wheat, which is not possible in Avi-single. The statistics discussed here are summarized in the GPP row of Table ). It is likely that the representation of both crops in Avignon in terms of the crop growth and development would be improved by developing specific crop parameterizations for winter wheat and sorghum in JULES.
Comparisons between the simulated energy fluxes (H and LE) and observations are shown in Figs. and d–i. The RMSE and bias values for H and LE are given in Table ; these are generally comparable to those from Table 5 in , which are LE: RMSE of 52.4 Wm-2, bias of -11.8 Wm-2; and H: RMSE of 56.2 Wm-2, bias of 17.6 Wm-2. The linear correlations for H are strong for all three simulations, with r values above 0.7 (Avi-grass, Fig. d; Avi-sequential, Fig. e; and Avi-single, Fig. f). The linear correlations for LE are more variable between the three simulations, with Avi-grass (Fig. g) having the strongest correlation (r value of 0.81), Avi-single (Fig. i) having the weakest correlation (r value of 0.64) and Avi-sequential (Fig. h) in between the two (r value of 0.73). The H values for all three simulations and the LE values for Avi-grass and Avi-sequential are comparable to but lower than those from Table 5 in , which provides correlation values of 0.8 for LE and 0.85 for H. The annual cycles of LE and H are shown in Fig. a and b, respectively. Generally, the seasonal cycles of H and LE are captured well in JULES (Fig. and the time series in Fig. b and c). The annual cycle for LE is close to observations in the first half of the year but too high in the second half for Avi-grass and Avi-sequential. Avi-single is much too low, which explains its lower r value. Overall, H is closer to observations for all three simulations; however, the annual cycles show that both Avi-grass and Avi-single are a little too high and similarly Avi-sequential a little too low, explaining why the r values for this variable were much closer to each other.
India
The water resource requirement of a crop is affected by the size of the crop and its leaf area with overestimation (underestimation) of these characteristics potentially skewing the results towards a higher (lower) resource requirement. Canopy heights for both the regional simulation and the single-grid-box simulations for both rice and wheat across the region are between 0.5 and 0.7 m (not shown), which is an expected value for a typical crop, as described in . Similarly, LAIs in both the regional and single-grid-box simulations are consistent with each other. Figure shows the LAI for each of the four single-grid-box simulations. This shows that wheat LAI from JULES is between 5 and 7 m2m-2 across the locations; this is also an expected value for a crop according to . Rice LAI is lower (between 2 and 4 m2m-2), with the lowest values for WestUP, slightly increasing from west to east locations. For WestUP particularly, rice (solid red line) has a small LAI (Fig. ) but it generates a yield (red asterisks in Fig. ) that falls within the range of the observations for each year. However, wheat (solid black line) generates an LAI that is closer to expected values but a smaller yield compared with observations (Fig. , black asterisks). The inclusion of both observation datasets in Fig. highlights the spread between yield estimates for this region.
Figure shows single-grid-box yields for both the India-single (asterisks) and India-sequential (circles), illustrating that the wheat crops simulated in both India-single and India-sequential are very similar to each other. Figure compares the JULES yields with observed yields from , indicating that the yields change at each of the locations from west to east. However, comparisons with observed yields should not be overinterpreted because the combined stresses that will be implicit in the observations, such as nutrient and water stress are not accounted for in these simulations. We also assume that rice and wheat are grown in rotation everywhere, which will be an overestimation. The observed yields, particularly for wheat are larger to the west and reduce to the east (shown for the whole region for rice – Fig. a; and wheat – Fig. c). In both the India-single and India-sequential crop simulations, JULES simulates a reduction in wheat yields from west to east; however, this is not as extreme in the simulations as it is in the observations because JULES tends to underestimate the wheat yield especially for the western locations. Across both observation datasets, the bias is smaller for the other locations, with rice and wheat yields within the range of the observations for most years (Fig. ). For EastBi, the rice yields are often toward the top of the range provided by the two observed datasets but still within the range of the observations (Fig. ). The spatial distributions of JULES rice and wheat yields from the regional simulation are shown in Fig. and reflect the results shown in the single-grid-box simulations.
Crop growth and developmentAvignon
Figure shows that there are only small differences between Avi-single (dotted line) and Avi-sequential (dashed), particularly for total aboveground biomass (Fig. a) and canopy height (Fig. c). The LAI in JULES is typically more sensitive to changes in conditions than the biomass or canopy height, but even these differences between Avi-single and Avi-sequential LAI are small. The observations show that there was a sorghum crop during the summer immediately before the 2008 and 2012 wheat crop; during these years, the LAI of Avi-single is slightly larger than that of the Avi-sequential and observations. In Avi-single, this sorghum crop is not present, which could affect the condition of the soils in the model at the time the wheat is sown. Overall, Fig. shows that a similar wheat crop is simulated in both the single- and sequential-crop simulations, which is expected because both use the same parameterization to calculate crop development. Including a secondary crop on the same field in Avi-sequential does not really affect the primary crop growth and development, which means that Avi-single and Avi-sequential simulate a very similar primary (wheat) crop.
India single-grid-box simulations
The LAI in Fig. together with the yields in Fig. show that the wheat crop in India-sequential is indistinguishable from the wheat crop in India-single. The canopy height is very similar to the LAI; therefore, it is not shown here. In each of the figures showing LAI and yield, it is only occasionally possible to see the two simulations, highlighting that the crop growth and development of the primary crop is not really changed by the implementation of a secondary crop at the same location in JULES.
Time series of total aboveground biomass (a), LAI (b) and canopy height (c) for the Avignon site for wheat (black) and sorghum (red) for observations (solid lines) and simulations using the observed sowing and harvest dates: Avi-sequential (dashed) and the wheat-only Avi-single (dotted) simulations for the period between 2005 and 2013 using observed sowing and harvest dates. Simulations with prescribed LAI and canopy height are not shown here as these follow the observed LAI and canopy height. The dotted red sorghum line in the Avi-single simulation is not visible because this is at zero for the wheat-only simulation. Observed aboveground biomass in panel (a) is shown by purple asterisks. The standard deviation of the measurements is shown to represent the uncertainty in the observations.
Time series of GPP (a), H(b) and LE (c) for the Avignon site compared with observations (black lines). For H(b) and LE (c), the whole period from 2005 to 2012 is shown, while GPP (a) shows the period 2005–2007 due to availability of observations. In the GPP plot, only one complete winter wheat (yellow) and one complete sorghum season (pink) are highlighted. The following model simulations are also shown: Avi-grass (red), Avi-sequential (blue) and a wheat-only Avi-single (cyan). In each plot, a 10 d smoothing has been applied to the daily data.
Scatter plot comparing the observed rice and wheat yields against JULES simulations at each of the India locations shown in Fig. , with India-sequential shown by circles (rice in red and wheat in black). India-single shown by asterisks (rice shown in red and wheat shown in blue). Note that the simulations assume that rice and wheat are grown everywhere and the crops are not limited by pests and diseases or nutrient availability, which will affect the comparison with the observations that implicitly include these aspects.
Annual climatology of fluxes (in day of year) for WestUP for India-sequential (black) and India-single (red): carbon fluxes: NPP (a) and GPP (b). Heat fluxes: sensible heat H(c) and latent heat LE (d). Soil moisture variables: β(e) and soil moisture availability in the top 1 m of soil (f). Moisture fluxes: evapotranspiration (g) and non-evapotranspiration moisture fluxes (h).
Energy and carbon fluxesAvignon
The time series of GPP (Fig. a) shows that all three simulations largely follow each other closely (Avi-sequential, blue line; Avi-single, cyan line; Avi-grass, red line), except during the sorghum growing period, which is not represented using the single-crop version of JULES. During the sorghum season, the land cover in JULES is therefore represented by bare soil, so the GPP goes to zero (dashed cyan line in Fig. a). The wheat season fluxes of H (Fig. b) and LE (Fig. c) also follow each other closely for all three simulations. However, during the sorghum growing period, the effect on the energy fluxes is more complex, with Avi-single having a lower LE and higher H than Avi-sequential and Avi-grass. In Avi-grass and Avi-sequential, it is possible to include representation of the two observed crops and therefore represent the actual land cover at the Avignon site; this is not possible in Avi-single. The inclusion of this secondary crop on the same field in Avi-sequential modifies the energy and carbon fluxes for the part of the year the secondary crop is being represented but the primary crop fluxes remain similar to those in Avi-single.
India single-grid-box simulations
Wheat is simulated in both India-single and India-sequential; it is grown during the last 25 d of one year and the first 140 d of the following year. The carbon fluxes during the wheat season are very similar in both India-sequential and India-single; these are shown for WestUP GPP in Fig. a and NPP in Fig. b and similarly for the other India locations in Appendix B, Figs. a–b, a–b and a–b. The peak for wheat is therefore the first peak shown in NPP (Fig. a) and GPP (Fig. b). Wheat NPP begins to decline at around day 41 of the year in both India-sequential and India-single, which is quite early in the season and may have a direct impact on the yield. The decline in NPP could be related to the way the carbon is partitioned to different parts of the plant. A short time series of India-sequential showing how carbon is partitioned to the different parts of the plant for wheat (black) and rice (red) is shown in Appendix B (Fig. ). The allometric coefficients that control the partitioning of carbon to the different parts of the crop in JULES are currently those from ; it is possible that the results could be improved for South Asia if these were tuned to more appropriate values for the crops there. Rice is planted between day 141 and day 330 of a typical year; this is not represented in the India-single simulation; therefore, the NPP and GPP go to zero during the rice season.
Similar to the carbon fluxes, Fig. shows the fluxes of H (Fig. c) and LE (Fig. d) for WestUP and the other India locations in Appendix B, Figs. c–d, c–d and c–d. At each location, the energy fluxes during the wheat season for both India-sequential and India-single are largely indistinguishable from each other. The main differences in the energy fluxes occur during the rice season, where India-single has a lower LE and higher H than India-sequential. We investigate the changes in LE further, using the evapotranspiration (for WestUP shown in Fig. g) and non-evapotranspiration moisture fluxes (for WestUP shown in Fig. h). The sequential-crop method affects the individual components of the moisture flux, resulting in a larger component from the evapotranspiration flux and a lower component from non-evapotranspiration fluxes than India-single; this effect is visible across all four India locations (Appendix B, Figs. , and g and h).
Soil moisture
The crop rotation in Avignon has long fallow periods between crops, which means that the effect of sequential crops on soil moisture will also be influenced by soil evaporation. The crop rotation in India does not have long fallow periods, which allows analysis of the effects of sequential crops on soil moisture for the India single-grid-box simulations.
India single-grid-box simulations
The effect of sequential crops on the soil moisture availability factor (β, Sect. ) is similar to the carbon and heat fluxes, with India-sequential and India-single being consistent with each other for the wheat season. Figure e shows β for WestUP (β at the other India locations are shown in Appendix B, Figs. e, e and e). β is mainly affected during the secondary growing period (rice) with India-sequential showing a larger decrease in soil moisture availability over a shorter part of the year than India-single (Fig. e). Unlike β and the fluxes of carbon and heat, the available soil moisture in the top 1.0 m is affected by sequential cropping throughout the year; this is shown for WestUP in Fig. f. The available soil moisture in the top 1.0 m has much larger variations in India-sequential than in India-single even outside the rice season across the four locations, but particularly for the westernmost locations of WestUP (Fig. f) and EastUP (Fig. f). The effect of sequential crops on the available soil moisture in the top 1.0 m reduces for the more eastern locations of WestBi (Fig. f) and EastBi (Fig. f).
These moisture fields show that WestUP has the lowest available soil moisture and therefore β value, suggesting this location is likely to be the most water stressed. It is also the western locations that are more affected by the implementation of sequential crops. WestBi, on the other hand, is least affected by the implementation of sequential crops, often with the highest β and the most consistent available soil moisture in the top 1.0 m across the year of the four locations. This is consistent with the temperature and precipitation time series shown in Fig. , where the locations to the east are wetter and cooler than those to the west. This means there is more available soil moisture in the top 1.0 m for the eastern locations compared with the western locations.
Discussion
The sequential cropping system is used around the world especially in the tropics but is relatively underrepresented in land-surface models and consequently climate simulations. Most crop models are only able to simulate one crop per year, which limits their use in many parts of the world. Often this means that only one of the crops in a rotation are modelled explicitly, with the rest of the year represented by bare or almost bare soil. In a single-crop model, the correct land cover is not represented for the season that is not being simulated, which means that water, carbon and energy fluxes are incorrect for this season. An alternative approach used by ISIMIP simulates each crop for a small fraction of each grid box, thereby producing a yield for each crop for each grid box; this is then post-processed to give larger-scale yields. The ISIMIP approach is effective if yields are the primary focus, but it is not applicable for understanding the use of water resources across different seasons. The real benefit of the sequential cropping capability is that it can represent the observed land cover in a single simulation, which improves the simulation of water, energy and carbon fluxes. Including sequential crops enables a continuous simulation for multiple years and seasons, which provides an understanding of the evolving demand for resources.
We present a method for simulating sequential crops in JULES with a view to enabling JULES-crop to represent these cropping systems in Earth system and climate models for application in adaptation and mitigation studies. We demonstrate its use for a site in Avignon, France, which uses the sorghum–winter wheat rotation and the Indian states of Uttar Pradesh and Bihar, where the rice–spring wheat rotation is prevalent. The wheat season is represented in all sequential- and single-crop simulations, producing a similar wheat crop in terms of LAI, canopy height and yield (or as observed in Avignon, total aboveground biomass). This is expected because the sequential-crop method does not change the parameterization of crops in JULES but changes the control of the crop code to allow it to be used more than once per year. This means that the crop growth and development are not affected by the inclusion of sequential cropping in JULES for neither Avignon nor India. For Avignon, there are small differences in the LAI of the single-crop simulation but only for wheat seasons which occur immediately after a sorghum season that is not being represented. These small differences in LAI in Avignon could be due to differences in soil conditions, where sorghum is not being simulated but the irrigation water is still being applied via the precipitation.
The representation of crops, either using the crop model (Avi-sequential) or using grasses (Avi-grass), has a similar effect on the surface fluxes, which shows that the code used by the wider vegetation in JULES as well as the crop model is approximating leaf-level photosynthesis, stomatal conductance, water stress and leaf-to-canopy scaling correctly when compared with observations. For simulations that include both crops (Avi-sequential and Avi-grass), fluxes of carbon and heat are generally captured well for Avignon compared with observations. There are regular differences in the fluxes during the wheat-only single-crop simulations at all locations. The single-crop simulation cannot simulate the correct land cover during the second crop of the rotation; for Avignon, this is the sorghum season, and for India, the rice season. In Avignon, compared with observations, the LE is too low and H is too high during the sorghum season not being explicitly represented in JULES. The surface fluxes for the India grid box simulations are consistent with those in Avignon, with a higher H and lower LE in India-single than in India-sequential. In the sequential simulations, the representation of the land surface is closer to the observed land cover because bare soil and the relevant crops can be reproduced within the simulation when they are typically observed. This analysis shows that the main heat and carbon fluxes are changed by the implementation of sequential cropping, with the main differences occurring during the period where a crop is observed but is not usually represented in the simulation.
The regular rice–wheat rotation across India does not have long periods of bare soil, allowing the effect of sequential crops on soil moisture to be considered. The climate of Uttar Pradesh and Bihar, two of the main producers of rice and wheat in India, is highly variable, both in terms of temperature (ranging from 7 to 52 ∘C) and rainfall (between 0 and 15 mmd-1). Across these states, there is a cooling moistening trend from west to east, making conditions for growing crops very different across a relatively limited area. In the India-single simulation, the available soil moisture is less variable than the India-sequential simulation; this is not limited to an individual crop season but is evident throughout the year. In general, this increased variability is more pronounced for the more arid western locations, with the wetter eastern locations being less affected. Sequential cropping is therefore modifying the soil moisture availability across all the crop seasons. The aridity of these more westerly locations could be exacerbated by the lack of irrigation in the RCM driving data that prescribe the JULES climate in these simulations. This region is intensively irrigated , which means that there is a significant contribution to atmospheric humidity from evaporation due to irrigation and the recycling of water into precipitation that cannot be accounted for here. estimate that as much as 35 % of the evaporation moisture from the Ganges Basin is recycling within the river basin. Therefore, in these simulations, we are missing the part of the water cycle that allows evaporation from the surface to affect the humidity.
WestUP has the least favourable conditions for crops, with the least available soil moisture and higher temperatures than the other locations, yet the observed yields are higher than the eastern locations. The observed yields at EastBi are the lowest of the four locations, where the cooler wetter conditions should be more conducive to achieving higher yields; these are neither observed nor modelled. This is likely to be due to the differing management practices between the two states of Uttar Pradesh and Bihar. Uttar Pradesh is characterized by high agricultural productivity with effective irrigation systems and early adoption of new management practices . Bihar, on the other hand, has lower agricultural productivity, farms tend to be smaller and more fragmented, irrigation systems are less effective , and adoption of new technology is also slower due to the lack of available machinery . Yield gap parameters are included in many crop models in order to account for the impact of differing nutrient levels, pests, diseases and non-optimal management , thus explaining the difference between potential and actual yield under the same environment . This is not included in these simulations.
The regional sequential simulation shows that the sequential cropping method presented here can be used at larger scales. From the analysis for Avignon and the four India grid boxes, we anticipate that the inclusion of sequential crops will have a significant impact on the surface fluxes of energy, carbon and soil moisture at regional scales. However, in the simulation presented, it is assumed that crops are grown everywhere across Uttar Pradesh and Bihar. For a fair and comprehensive evaluation of the impacts of sequential crops at regional scales, this simulation requires further development to represent the types of crops, together with where and when they are grown. Future work will compare a regional sequential simulation (based on the larger-scale simulation presented here), with an equivalent single-crop simulation that includes realistic crop fractions, to understand the impact of sequential crops on a regional scale. Irrigation is also important for this region; while this is included in these simulations, it is only applied during the dry wheat season with no irrigation during the monsoon rice season. Some irrigation may also occur during monsoon breaks if there is a significant drop in soil moisture; this would be a useful future JULES development. Future developments to include the effects of heat stress, soil nutrient limitations, pests and diseases on crops in JULES would also be useful. The simulations shown here consider a small number of rotations, crops and regions. However, different varieties and types of crops, the timing of sowing and harvesting, together with many possible irrigation options, can have a large impact on the model results. This is an important consideration for future work and should be investigated fully when applying this method to new areas.
Conclusions
In this paper, we describe and implement a new sequential cropping capability in JULES enabling more than one crop to be simulated at a given location during a growing season (available in JULES from version 5.7). This is another step towards being able to include JULES-crop in Earth system and climate models, allowing more accurate representation of land use and surface coverage in regions where two or more crops are grown in rotation. This is important for accurate representation of the fluxes of carbon and energy in climate simulations in regions with more than one crop grown in one location within a year. Sequential crops reduce H and increase LE, thereby modifying the contribution from the different components of the latent heat flux. For crop rotations without long fallow periods, including sequential crops also affects the availability of soil moisture throughout the year. More realistic fluxes of carbon, water and energy are important for understanding the impacts of climate change. The continuous simulation of all crops throughout the year also provides a more complete picture of the total demand for water resources, which is important for climate impact assessments.
We show that JULES simulates two realistic crops in a growing season both in Avignon (wheat and sorghum) and across Uttar Pradesh and Bihar (rice and wheat). For Avignon and India, the maxima of LAI, canopy height and biomass (yield) occur at approximately the correct times of the year. For Avignon, the observed GPP and energy (H and LE) fluxes are reproduced, correlating well with observations. There are only small differences between using the crop model and using grasses to represent the crops at this site, indicating that JULES-crop can reproduce the LAI and canopy height well enough to compare well with the observed surface fluxes. The representation of crops in Avignon could be improved by including crop-specific parameterizations of winter wheat and sorghum in the model, although sorghum would probably require two different sets of parameters for a significant improvement because the two varieties grown at the site are so different.
A regional simulation provides the basis for future sequential-crop simulations for a larger area. This cropping system is likely to be a feature of the future land surface, not just in the tropics but globally as an adaptation to climate change. Therefore, we encourage other modelling communities to include sequential crops in their models, so that it can benefit from being part of future model intercomparison projects such as AgMIP and ISIMIP and be developed further.
Avignon comparison
Time series of temperature (a) and precipitation which includes the observed irrigation amounts added at the exact day and time they were applied to the crops (b) in Avignon for the time period analysed (2005–2012).
Annual cycle of H(a) and LE (b) compared with observations (black line) at the Avignon site between 2005 and 2013. Annual cycles for the simulations are also shown: Avi-grass (red line), Avi-sequential (blue line) and Avi-single (cyan line).
Comparison of observations and model fluxes at the Avignon site: GPP between 2005 and 2008 for Avi-grass (a), Avi-sequential (b) and Avi-single (c). H between 2005 and 2013 for Avi-grass (d) and Avi-sequential (e), Avi-single (f). LE between 2005 and 2013 for Avi-grass (g), Avi-sequential (h) and Avi-single (i).
India single-grid-box comparison
Time series of monthly precipitation (a), temperature (b) and vapour pressure deficit (c) at each of the India locations shown by the solid lines (WestUP – black, EastUP – red, WestBi – blue and EastBi – cyan). Panel (b) also shows the minimum (“×”) and maximum (“+”) temperatures for each of the locations for each month together with the JULES cardinal temperatures (horizontal lines) for rice (green) and wheat (orange): max temperatures (dotted line), optimum temperatures (solid line) and base temperatures (dashed line).
Time series of crop harvest pool (solid lines) with the JULES yield from the sequential-crop run at the time it is output by the model (asterisks) for rice (red) and wheat (black) at each of the India locations shown in Fig. . Also shown are two sets of observations: annual yields from shown by the filled circles and 5-year averages from shown by the filled triangles (following the same colours with rice shown in red and wheat in black).
Time series of the LAI of rice (red) and wheat (black) at each of the India locations shown in Fig. . India-sequential is shown by solid lines and India-single is shown by dashed lines although they are indistinguishable from each other.
Annual climatology of fluxes (in day of year) for EastUP for India-sequential (black) and India-single (red): carbon fluxes: NPP (a) and GPP (b). Heat fluxes: sensible heat H(c) and latent heat LE (d). Soil moisture variables: β(e) and soil moisture availability in the top 1 m of soil (f). Moisture fluxes: evapotranspiration (g) and non-evapotranspiration moisture fluxes (h).
Annual climatology of fluxes (in day of year) for WestBi for India-sequential (black) and India-single (red): carbon fluxes: NPP (a) and GPP (b). Heat fluxes: sensible heat H(c) and latent heat LE (d). Soil moisture variables: β(e) and soil moisture availability in the top 1 m of soil (f). Moisture fluxes: evapotranspiration (g) and non-evapotranspiration moisture fluxes (h).
Annual climatology of fluxes (in day of year) for EastBi for India-sequential (black) and India-single (red): carbon fluxes: NPP (a) and GPP (b). Heat fluxes: sensible heat H(c) and latent heat LE (d). Soil moisture variables: β(e) and soil moisture availability in the top 1 m of soil (f). Moisture fluxes: evapotranspiration (g) and non-evapotranspiration moisture fluxes (h).
Time series of each crop carbon pool: leaf (solid lines), root (dashed), stem (dotted) and harvest (dash–dot) with the JULES yield at the time it is output by the model (asterisks) for rice (red) and wheat (black) at each of the India locations shown in Fig. for a subset of years of the simulation between 1998 and 2001.
India regional simulation
The values used in the regional JULES ancillary for rice. Sowing date (a) and latest possible harvest date (b), both in units of day of year. Thermal time for the vegetative stage (c) and thermal time for the reproductive stage (d), both in units of degree days
.
The values used in the regional JULES ancillary for wheat. Sowing date (a) and latest possible harvest date (b), both in units of day of year. Thermal time for the vegetative stage (c) and thermal time for the reproductive stage (d), both in units of degree days.
A comparison of observed rice yields (a) with JULES rice yields (b) and observed wheat yields (c) with JULES wheat yields (d) for the period 1991–2007 across Uttar Pradesh and Bihar. Note that the simulations assume that rice and wheat are grown everywhere and the crops are not limited by nutrient availability.
Code and data availability
To use the JULES code, registration is required and code is freely available subject to completion of a software
licence .
The version of the model used in this analysis is an enhanced JULESvn5.2; this branch is available at https://code.metoffice.gov.uk/trac/jules/browser/main/branches/dev/camillamathison/vn5.2_croprotate_irrigtiles.
The developments contained in this branch are now implemented into the trunk of JULES from version 5.7. The regional climate model datasets used will hopefully be available via the Centre for Environmental Data Analysis (CEDA) catalogue. The Avignon rose suite that uses the vn5.2 branch is available from this link: https://code.metoffice.gov.uk/trac/roses-u/browser/a/y/2/2/5/u-ay225-cropmodel in order to aid future model development. Please seek authorization from andre.chanzy@inrae.fr and albert.olioso@inrae.fr or contact Sebastien Garrigues (sgarrig@gmail.com) to use the Avignon dataset, which is published in .
Author contributions
AJC, PF and AW provided general scientific guidance throughout the paper and helped prioritize the initial requirement for including sequential cropping in JULES, toward the aim of modelling integrated impacts for South Asia. CD provided gridded observations of crop yield courtesy of ICRISAT. SG and SM provided Avignon data. KW provided the code for generating the crop ancillary files and an initial suite for Avignon without crops.
Competing interests
The authors declare that they have no conflict of interest.
Acknowledgements
The authors thank the INRAE EMMAH laboratory in Avignon (France) in charge of the “flux and remote sensing site” for the supply of the data. The authors thank the reviewers and the editor for their comments and suggestions. The authors thank Kate Halladay for helping Camilla Mathison understand the soil moisture and irrigation implications.
Financial support
The research leading to these results has received funding from the European Union Seventh Framework Programme FP7/2007–2013 under grant agreement no. 603864. Camilla Mathison, Pete Falloon, Andy Wiltshire and Karina Williams were supported by the Met Office Hadley Centre Climate Programme funded by BEIS and Defra.
Review statement
This paper was edited by Wolfgang Kurtz and reviewed by Katharina Waha and three anonymous referees.
ReferencesAllen, R. G., Pereira, L. S., Raes, D., and Smith, M.: Crop evapotranspiration
– Guidelines for computing crop water requirements – FAO Irrigation and
drainage paper 56, Food and Agriculture Organization of the United Nations,
available at: http://www.fao.org/3/X0490E/x0490e00.htm (last access: June 2019), 1998.Best, M. J., Pryor, M., Clark, D. B., Rooney, G. G., Essery, R. L. H., Ménard, C. B., Edwards, J. M., Hendry, M. A., Porson, A., Gedney, N., Mercado, L. M., Sitch, S., Blyth, E., Boucher, O., Cox, P. M., Grimmond, C. S. B., and Harding, R. J.: The Joint UK Land Environment Simulator (JULES), model description Part 1: Energy and water fluxes, Geosci. Model Dev., 4, 677–699, 10.5194/gmd-4-677-2011, 2011.Betts, R. A.: Integrated approaches to climate-crop modelling: needs and
challenges, Philos. T. R. Soc. B, 360, 2049–2065, 10.1098/rstb.2005.1739, 2005.Bhattacharyya, T., Pal, D., Easter, M., Batjes, N., Milne, E., Gajbhiye, K.,
Chandran, P., Ray, S., Mandal, C., Paustian, K., Williams, S., Killian, K.,
Coleman, K., Falloon, P., and Powlson, D.: Modelled soil organic carbon
stocks and changes in the Indo-Gangetic Plains, India from 1980 to 2030,
Agr. Ecosyst. Environ., 122, 84–94,
10.1016/j.agee.2007.01.010,
soil carbon stocks at regional scales, 2007.Biemans, H., Speelman, L., Ludwig, F., Moors, E., Wiltshire, A., Kumar, P.,
Gerten, D., and Kabat, P.: Future water resources for food production in five South Asian river basins and potential for adaptation – A modeling study,
Sci. Total Environ., 468–469, Supplement, S117–S131,
10.1016/j.scitotenv.2013.05.092, 2013.Biemans, H., Siderius, C., Mishra, A., and Ahmad, B.: Crop-specific seasonal estimates of irrigation-water demand in South Asia, Hydrol. Earth Syst. Sci., 20, 1971–1982, 10.5194/hess-20-1971-2016, 2016.Bodh, S. P. C., Rai, S. J. P., Sharma, S. A., Gajria, S. P., Yadav, S. M.,
Virmani, S. S., and Pandey, S. R.: Agricultural Statistics at a Glance 2015,
Ministry of Agriculture and Farmers welfare, Directorate of Economics and
Statistics, available at: http://eands.dacnet.nic.in (last access: June 2019), 2015.Bondeau, A., Smith, P. C., Zaehle, S., Schaphoff, S., Lucht, W., Cramer, W.,
Gerten, D., Lotze-Campen, H., Müller, C., Reichstein, M., and Smith, B.:
Modelling the role of agriculture for the 20th century global terrestrial
carbon balance, Glob. Change Biol., 13, 679–706,
10.1111/j.1365-2486.2006.01305.x,
2007.Caldwell, R. M. and Hansen, J. W.: Simulation of multiple cropping systems with
CropSys, pp. 397–412, Springer Netherlands, Dordrecht,
10.1007/978-94-011-2842-1_24, 1993.Calvet, J.-C., Noilhan, J., Roujean, J.-L., Bessemoulin, P., Cabelguenne, M.,
Olioso, A., and Wigneron, J.-P.: An interactive vegetation SVAT model tested
against data from six contrasting sites, Agr. Forest Meteorol.,
92, 73–95, 10.1016/S0168-1923(98)00091-4, 1998.Challinor, A., Wheeler, T., Craufurd, P., Slingo, J., and Grimes, D.: Design
and optimisation of a large-area process-based model for annual crops,
Agr. Forest Meteorol., 124, 99–120,
10.1016/j.agrformet.2004.01.002,
2004.Clark, D. B., Mercado, L. M., Sitch, S., Jones, C. D., Gedney, N., Best, M. J., Pryor, M., Rooney, G. G., Essery, R. L. H., Blyth, E., Boucher, O., Harding, R. J., Huntingford, C., and Cox, P. M.: The Joint UK Land Environment Simulator (JULES), model description Part 2: Carbon fluxes and vegetation dynamics, Geosci. Model Dev., 4, 701–722, 10.5194/gmd-4-701-2011, 2011.Cong, W.-F., Hoffland, E., Li, L., Six, J., Sun, J.-H., Bao, X.-G., Zhang, F.-S., and Werf, W. V. D.: Intercropping enhances soil carbon and nitrogen,
Glob. Change Biol., 21, 1715–1726, 10.1111/gcb.12738,
2015.Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M. A., P., G. B., Bauer, Bechtold, P., Beljaars, A. C. M., van de Berg, L., Bidlot, J., Bormann, N., Delsol, C., Dragani, R.,
Fuentes, M., Geer, A. J., Haimberger, L., Healy, S. B., Hersbach, H.,
Hólm, E. V., Isaksen, L., Kållberg, P., Köhler, M., Matricardi, M., McNally, A. P., Monge-Sanz, B. M., Morcrette, J.-J., Park, B.-K., Peubey, C., de Rosnay, P., Tavolato, C., Thépaut, J.-N., and Vitart, F.: The
ERA-Interim reanalysis: configuration and performance of the data
assimilation system, Q. J. Roy. Meteor. Soc.,
137, 553–597, 10.1002/qj.828,
10.1002/qj.828, 2011.Dury, J., Schaller, N., Garcia, F., Reynaud, A., and Bergez, J. E.: Models to
support cropping plan and crop rotation decisions. A review, Agron. Sustain. Dev., 32, 567–580, 10.1007/s13593-011-0037-x, 2012.Erenstein, O. and Laxmi, V.: Zero tillage impacts in India's rice-wheat
systems: A review, Soil Till. Res., 100, 1–14,
10.1016/j.still.2008.05.001,
2008.Erenstein, O., Farooq, U., Malik, R., and Sharif, M.: On-farm impacts of zero
tillage wheat in South Asia's rice-wheat systems, Field Crop. Res., 105,
240–252, 10.1016/j.fcr.2007.10.010,
2008.Essery, R. L. H., Best, M. J., and Cox, P. M.: MOSES 2.2 technical
documentation, Hadley Centre Technical Note, 30, available at:
https://library.metoffice.gov.uk/Portal/DownloadImageFile.ashx?objectId=1038 (last access: January 2021),
2001.Fischer, R.: Definitions and determination of crop yield, yield gaps, and of
rates of change, Field Crop. Res., 182, 9–18,
10.1016/j.fcr.2014.12.006, 2015.Frieler, K., Levermann, A., Elliott, J., Heinke, J., Arneth, A., Bierkens, M. F. P., Ciais, P., Clark, D. B., Deryng, D., Döll, P., Falloon, P., Fekete, B., Folberth, C., Friend, A. D., Gellhorn, C., Gosling, S. N., Haddeland, I., Khabarov, N., Lomas, M., Masaki, Y., Nishina, K., Neumann, K., Oki, T., Pavlick, R., Ruane, A. C., Schmid, E., Schmitz, C., Stacke, T., Stehfest, E., Tang, Q., Wisser, D., Huber, V., Piontek, F., Warszawski, L., Schewe, J., Lotze-Campen, H., and Schellnhuber, H. J.: A framework for the cross-sectoral integration of multi-model impact projections: land use decisions under climate impacts uncertainties, Earth Syst. Dynam., 6, 447–460, 10.5194/esd-6-447-2015, 2015.Garrigues, S., Olioso, A., Calvet, J. C., Martin, E., Lafont, S., Moulin, S., Chanzy, A., Marloie, O., Buis, S., Desfonds, V., Bertrand, N., and Renard, D.: Evaluation of land surface model simulations of evapotranspiration over a 12-year crop succession: impact of soil hydraulic and vegetation properties, Hydrol. Earth Syst. Sci., 19, 3109–3131, 10.5194/hess-19-3109-2015, 2015a.Garrigues, S., Olioso, A., Carrer, D., Decharme, B., Calvet, J.-C., Martin, E., Moulin, S., and Marloie, O.: Impact of climate, vegetation, soil and crop management variables on multi-year ISBA-A-gs simulations of evapotranspiration over a Mediterranean crop site, Geosci. Model Dev., 8, 3033–3053, 10.5194/gmd-8-3033-2015, 2015b.Garrigues, S., Boone, A., Decharme, B., Olioso, A., Albergel, C., Calvet, J.-C., Moulin, S., Buis, S., and Martin, E.: Impacts of the Soil Water
Transfer Parameterization on the Simulation of Evapotranspiration over a
14-Year Mediterranean Crop Succession, J. Hydrometeorol., 19,
3–25, 10.1175/JHM-D-17-0058.1, 2018.Goswami, B. and Xavier, P. K.: Dynamics of “interna” interannual variability
of the Indian summer monsoon in a GCM, J. Geophys. Res.-Atmos., 110, D24104, 10.1029/2005JD006042, 2005.Griffiths, F. E. W., Lyndon, R., and Bennett, M.: The Effects of Vernalization
on the Growth of the Wheat Shoot Apex, Ann. Bot.-London, 56, 501–511,
10.1093/oxfordjournals.aob.a087035, 1985.Harding, R., Blyth, E., Tuinenburg, O., and Wiltshire, A.: Land atmosphere
feedbacks and their role in the water resources of the Ganges basin, Sci. Total Environ., 468–469, Supplement, S85–S92,
10.1016/j.scitotenv.2013.03.016, 2013.Harper, A. B., Williams, K. E., McGuire, P. C., Duran Rojas, M. C., Hemming, D., Verhoef, A., Huntingford, C., Rowland, L., Marthews, T., Breder Eller, C., Mathison, C., Nobrega, R. L. B., Gedney, N., Vidale, P. L., Otu-Larbi, F., Pandey, D., Garrigues, S., Wright, A., Slevin, D., De Kauwe, M. G., Blyth, E., Ärdo, J., Black, A., Bonal, D., Buchmann, N., Burban, B., Fuchs, K., de Grandcourt, A., Mammarella, I., Merbold, L., Montagnani, L., Nouvellon, Y., Restrepo-Coupe, N., and Wohlfahrt, G.: Improvement of modelling plant responses to low soil moisture in JULESvn4.9 and evaluation against flux tower measurements, Geosci. Model Dev. Discuss. [preprint], 10.5194/gmd-2020-273, in review, 2020.Hatfield, J. L. and Prueger, J. H.: Temperature extremes: Effect on plant
growth and development, Weather Clim. Extremes, 10, 4–10,
10.1016/j.wace.2015.08.001, 2015.Hu, S., Mo, X., and Lin, Z.: Optimizing the photosynthetic parameter Vcmax by
assimilating MODIS-fPAR and MODIS-NDVI with a process-based ecosystem model,
Agr. Forest Meteorol., 198–199, 320–334,
10.1016/j.agrformet.2014.09.002,
2014.Hudson, R.: Management of Agricultural, Forestry, Fisheries and Rural
Enterprise – Volume I, EOLSS Publications, available at: https://books.google.co.uk/books?id=-eGvCwAAQBAJ (last access: June 2019), 2009.
ICRISAT: District Level Database Documentation, Tech. rep., International
Crops Research Institute for the Semi-Arid Tropics, Hyderabad, 2015.Iizumi, T. and Ramankutty, N.: How do weather and climate influence cropping
area and intensity?, Glob. Food Secur.-Agr., 4, 46–50,
10.1016/j.gfs.2014.11.003,
2015.Jones, R. G., Noguer, M., Hassell, D. C., Hudson, D., Wilson, S. S., Jenkins, G. J., and Mitchell, J. F.: Generating high resolution climate change
scenarios using PRECIS, Met Office Hadley Centre, Exeter, UK, 40 pp.,
available at: https://www.metoffice.gov.uk/binaries/content/assets/metofficegovuk/pdf/research/applied-science/precis/tech_man_v2.pdf (last access: January 2021), 2004.JULES Collaboration: Rose suite access page, available at: https://code.metoffice.gov.uk/trac/jules/browser/main/branches/dev/camillamathison/vn5.2_croprotate_irrigtiles,
last access: January 2021.Kumar, P., Wiltshire, A., Mathison, C., Asharaf, S., Ahrens, B., Lucas-Picher, P., Christensen, J. H., Gobiet, A., Saeed, F., Hagemann, S., and Jacob, D.:
Downscaled climate change projections with uncertainty assessment over India using a high resolution multi-model approach, Sci. Total Environ., 468–469, Supplement, S18–S30,
10.1016/j.scitotenv.2013.01.051,
2013.
Kumar, R., Singh, R., and Sharma, K.: Water resources of India, Current
Sci. India, 89, 794–811, 2005.Laik, R., Sharma, S., Idris, M., Singh, A., Singh, S., Bhatt, B., Saharawat, Y., Humphreys, E., and Ladha, J.: Integration of conservation agriculture
with best management practices for improving system performance of the
rice–wheat rotation in the Eastern Indo-Gangetic Plains of India,
Agr. Ecosyst. Environ., 195, 68–82,
10.1016/j.agee.2014.06.001,
2014.Liu, L., Xu, X., Zhuang, D., Chen, X., and Li, S.: Changes in the Potential
Multiple Cropping System in Response to Climate Change in China from
1960–2010, PLoS ONE, 8, e80990, 10.1371/journal.pone.0080990,
2013.Mahajan, A. and Gupta, R. D. (Eds.): The Rice–Wheat Cropping System,
pp. 109–117, Springer Netherlands, Dordrecht, 10.1007/978-1-4020-9875-8_7, 2009.Makino, A.: Rubisco and nitrogen relationships in rice: Leaf photosynthesis and
plant growth, Soil Sci. Plant Nutr., 49, 319–327,
10.1080/00380768.2003.10410016, 2003.Mathison, C., Wiltshire, A., Dimri, A., Falloon, P., Jacob, D., Kumar, P.,
Moors, E., Ridley, J., Siderius, C., Stoffel, M., and Yasunari, T.: Regional
projections of North Indian climate for adaptation studies, Sci. Total Environ., 468–469, Supplement, S4–S17,
10.1016/j.scitotenv.2012.04.066,
2013.Mathison, C., Wiltshire, A. J., Falloon, P., and Challinor, A. J.: South Asia river-flow projections and their implications for water resources, Hydrol. Earth Syst. Sci., 19, 4783–4810, 10.5194/hess-19-4783-2015, 2015.Mathison, C., Deva, C., Falloon, P., and Challinor, A. J.: Estimating sowing and harvest dates based on the Asian summer monsoon, Earth Syst. Dynam., 9, 563–592, 10.5194/esd-9-563-2018, 2018.Met Office Science Repository Service: Scientific Collaboration Trac page, available at: https://code.metoffice.gov.uk/trac, last access: January 2021.Monfreda, C., Ramankutty, N., and Foley, J. A.: Farming the planet: 2.
Geographic distribution of crop areas, yields, physiological types, and net
primary production in the year 2000, Global Biogeochem. Cy., 22, GB1022, 10.1029/2007GB002947,
2008.Mueller, B., Hauser, M., Iles, C., Rimi, R. H., Zwiers, F. W., and Wan, H.:
Lengthening of the growing season in wheat and maize producing regions,
Weather Clim. Extremes, 9, 47–56,
10.1016/j.wace.2015.04.001, 2015.
Nicklin, K.: Seasonal crop yield forecasting in semi-arid West Africa, chapt. 4, PhD
thesis, University of Leeds, 1, 2012.Noilhan, J. and Planton, S.: A Simple Parameterization of Land Surface
Processes for Meteorological Models, Mon. Weather Rev., 117, 536–549,
10.1175/1520-0493(1989)117<0536:ASPOLS>2.0.CO;2, 1989.Ogbaga, C.: Regulation of Photosynthesis in Sorghum in response to drought, PhD
Thesis, University of Manchester, 1, 1–186, 10.13140/RG.2.1.4756.5208,
2014.Olsovska, K., Kovar, M., Brestic, M., Zivcak, M., Slamka, P., and Shao, H. B.:
Genotypically Identifying Wheat Mesophyll Conductance Regulation under
Progressive Drought Stress, Front. Plant Sci., 7, 1111,
10.3389/fpls.2016.01111,
2016.Osborne, T., Gornall, J., Hooker, J., Williams, K., Wiltshire, A., Betts, R., and Wheeler, T.: JULES-crop: a parametrisation of crops in the Joint UK Land Environment Simulator, Geosci. Model Dev., 8, 1139–1155, 10.5194/gmd-8-1139-2015, 2015.Penning de Vries, F. W. T., Jansen, D., ten Berge, H., and Bakema, A.: Simulation of Ecophysiological Processes of Growth in Several Annual Crops,
Simulation monographs, Pudoc,
available at: https://books.google.co.uk/books?id=G-cjxJlr71wC (last access: August 2019), 1989.Petrie, C. A., Singh, R. N., Bates, J., Dixit, Y., French, C. A. I., Hodell, D. A., Jones, P. J., Lancelotti, C., Lynam, F., Neogi, S., Pandey, A. K.,
Parikh, D., Pawar, V., Redhouse, D. I., and Singh, D. P.: Adaptation to
Variable Environments, Resilience to Climate Change: Investigating Land,
Water and Settlement in Indus Northwest India, Curr. Anthropol., 58,
1–30, 10.1086/690112,
2017.Pires, G. F., Abrahão, G. M., Brumatti, L. M., Oliveira, L. J., Costa, M. H.,
Liddicoat, S., Kato, E., and Ladle, R. J.: Increased climate risk in
Brazilian double cropping agriculture systems: Implications for land use in
Northern Brazil, Agr. Forest Meteorol., 228–229, 286–298,
10.1016/j.agrformet.2016.07.005,
2016.Portmann, F. T., Siebert, S., and Döll, P.: MIRCA2000-Global monthly
irrigated and rainfed crop areas around the year 2000: A new high-resolution
data set for agricultural and hydrological modeling, Global Biogeochem. Cy., 24, GB1011,
10.1029/2008GB003435, gB1011, 2010.Ramankutty, N., Evan, A. T., Monfreda, C., and Foley, J. A.: Farming the
planet: 1. Geographic distribution of global agricultural lands in the year
2000, Global Biogeochem. Cy., 22, GB1003, 10.1029/2007GB002952, 2008.Ray, D. K., Ramankutty, N., and, N. D. M.: Recent patterns of crop yield
growth and stagnation, Nat. Commun., 3, 1293–1300,
10.1038/ncomms2296, 2012a.Ray, D. K., Ramankutty, N., Mueller, N. D., West, P. C., and Foley, J. A.:
Harvested Area and Yield for 4 Crops (1995–2005), available at:
http://www.earthstat.org/harvested-area-yield-4-crops-1995-2005/ (last access: August 2019),
2012b.Rivington, M. and Koo, J.: Report on the Meta-Analysis of Crop Modelling for
Climate Change and Food Security Survey, Climate Change, Agriculture and Food
Security Challenge Program of the CGIAR, available at:
https://cgspace.cgiar.org/rest/bitstreams/9114/retrieve (last access: July 2019),
2010.Robertson, M., Brooking, I., and Ritchie, J.: Temperature Response of
Vernalization in Wheat: Modelling the Effect on the Final Number of Mainstem
Leaves, Ann. Bot.-London, 78, 371–381, 10.1006/anbo.1996.0132, 1996.Rosenzweig, C., Jones, J., Hatfield, J., Ruane, A., Boote, K., Thorburn, P.,
Antle, J., Nelson, G., Porter, C., Janssen, S., Asseng, S., Basso, B., Ewert, F., Wallach, D., Baigorria, G., and Winter, J.: The Agricultural Model
Intercomparison and Improvement Project (AgMIP): Protocols and pilot studies,
Agr. Forest Meteorol., 170, 166–182,
10.1016/j.agrformet.2012.09.011, 2013.Rosenzweig, C., Elliott, J., Deryng, D., Ruane, A. C., Müller, C., Arneth, A., Boote, K. J., Folberth, C., Glotter, M., Khabarov, N., Neumann, K.,
Piontek, F., Pugh, T. A. M., Schmid, E., Stehfest, E., Yang, H., and Jones, J. W.: Assessing agricultural risks of climate change in the 21st century in
a global gridded crop model intercomparison, P. Natl. Acad. Sci. USA, 111, 3268–3273, 10.1073/pnas.1222463110, 2014.ROSE Suite JULES Collaboration: Sequential Crop rose suite access page, available at:
https://code.metoffice.gov.uk/trac/roses-u/browser/main/branches/dev/camillamathison/vn5.2_croprotate_irrigtiles, last access: January 2021.Schaphoff, S., von Bloh, W., Rammig, A., Thonicke, K., Biemans, H., Forkel, M., Gerten, D., Heinke, J., Jägermeyr, J., Knauer, J., Langerwisch, F., Lucht, W., Müller, C., Rolinski, S., and Waha, K.: LPJmL4 a dynamic global vegetation model with managed land Part 1: Model description, Geosci. Model Dev., 11, 1343–1375, 10.5194/gmd-11-1343-2018, 2018.
Sacks, W. J., Deryng, D., Foley, J. A., and Ramankutty, N.: Crop planting
dates: an analysis of global patterns, Global Ecol. Biogeogr., 19,
607–620, 2010.Shannon, S., Smith, R., Wiltshire, A., Payne, T., Huss, M., Betts, R., Caesar,
J., Koutroulis, A., Jones, D., and Harrison, S.: Global glacier volume
projections under high-end climate change scenarios, The Cryosphere, 13,
325–350, 10.5194/tc-13-325-2019, 2019.Sharma, B. and Sharma, H.: Status of Rice Production in Assam, India, J.
Rice Res., 3, e121, 10.4172/2375-4338.1000e121, 2015.
Simmons, A., Uppala, S., Dee, D., and Kobayashi, S.: ERA-Interim: New ECMWF
reanalysis products from 1989 onwards, ECMWF Newsletter – Winter 2006/07,
110, 25–35, 2007.Sinclair, T., Jr, P. P., Kimball, B., Adamsen, F., LaMorte, R., Wall, G.,
Hunsaker, D., Adam, N., Brooks, T., Garcia, R., Thompson, T., Leavitt, S.,
and Matthias, A.: Leaf nitrogen concentration of wheat subjected to elevated
[CO2] and either water or N deficits, Agr. Ecosyst. Environ.,
79, 53–60, 10.1016/S0167-8809(99)00146-2,
2000.Tuinenburg, O. A., Hutjes, R. W. A., Stacke, T., Wiltshire, A., and
Lucas-Picher, P.: Effects of irrigation in india on the atmospheric water
budget, J. Hydrometeorol., 15, 1028–1050,
10.1175/JHM-D-13-078.1, 2014.Waha, K., van Bussel, L. G. J., Müller, C., and Bondeau, A.: Climate-driven
simulation of global crop sowing dates, Global Ecol. Biogeogr., 21,
247–259, 10.1111/j.1466-8238.2011.00678.x, 2012.Waha, K., Müller, C., Bondeau, A., Dietrich, J., Kurukulasuriya, P.,
Heinke, J., and Lotze-Campen, H.: Adaptation to climate change through the
choice of cropping system and sowing date in sub-Saharan Africa, Global Environ. Chang., 23, 130–143,
10.1016/j.gloenvcha.2012.11.001,
2013.Wang, E., Martre, P., Zhao, Z., Ewert, F., Maiorano, A., Rötter, R. P.,
Kimball, B. A., Ottman, M. J., Wall, G. W., White, J. W., Reynolds, M. P.,
Alderman, P. D., Aggarwal, P. K., Anothai, J., Basso, B., Biernath, C.,
Cammarano, D., Challinor, A. J., Sanctis, G. D., Doltra, J., Dumont, B.,
Fereres, E., Garcia-Vila, M., Gayler, S., Hoogenboom, G., Hunt, L. A.,
Izaurralde, R. C., Jabloun, M., Jones, C. D., Kersebaum, K. C., Koehler, A.-K., Liu, L., Müller, C., Kumar, S. N., Nendel, C., O'Leary, G.,
Olesen, J. E., Palosuo, T., Priesack, E., Rezaei, E. E., Ripoche, D., Ruane, A. C., Semenov, M. A., Shcherbak, I., Stöckle, C., Stratonovitch, P.,
Streck, T., Supit, I., Tao, F., Thorburn, P., Waha, K., Wallach, D., Wang, Z., Wolf, J., Zhu, Y., and Asseng, S.: The uncertainty of crop yield
projections is reduced by improved temperature response functions, Nat. Plants, 3, 17102, 10.1038/nplants.2017.102, 2017.Warszawski, L., Friend, A., Ostberg, S., Frieler, K., Lucht, W., Schaphoff, S.,
Beerling, D., Cadule, P., Ciais, P., Clark, D. B., Kahana, R., Ito, A.,
Keribin, R., Kleidon, A., Lomas, M., Nishina, K., Pavlick, R., Rademacher, T. T., Buechner, M., Piontek, F., Schewe, J., Serdeczny, O., and
Schellnhuber, H. J.: A multi-model analysis of risk of ecosystem shifts under
climate change, Environ. Res. Lett., 8, 044018, 10.1088/1748-9326/8/4/044018,
2013.Warszawski, L., Frieler, K., Huber, V., Piontek, F., Serdeczny, O., and Schewe, J.: The Inter-Sectoral Impact Model Intercomparison Project (ISI-MIP):
Project framework, P. Natl. Acad. Sci. USA, 111,
3228–3232, 10.1073/pnas.1312330110, 2014.Williams, K., Gornall, J., Harper, A., Wiltshire, A., Hemming, D., Quaife, T., Arkebauer, T., and Scoby, D.: Evaluation of JULES-crop performance against site observations of irrigated maize from Mead, Nebraska, Geosci. Model Dev., 10, 1291–1320, 10.5194/gmd-10-1291-2017, 2017.Williams, K. E. and Falloon, P. D.: Sources of interannual yield variability
in JULES-crop and implications for forcing with seasonal weather forecasts,
Geosci. Model Dev., 8, 3987–3997, 10.5194/gmd-8-3987-2015,
2015.Williams, K. E., Harper, A. B., Huntingford, C., Mercado, L. M., Mathison, C. T., Falloon, P. D., Cox, P. M., and Kim, J.: How can the First ISLSCP Field Experiment contribute to present-day efforts to evaluate water stress in JULESv5.0?, Geosci. Model Dev., 12, 3207–3240, 10.5194/gmd-12-3207-2019, 2019.Xue, W.: Evaluation of biophysical factors driving temporal variations in C
gain, water use and yield production in Rice, PhD thesis, Department of Plant
Ecology, University of Bayreuth, 1, 1–230, available at:
https://www.researchgate.net/publication/295562013_Evaluation_of_biophysical_factors_driving_temporal_variations_in_carbon_gain_water_use_and_yield_production_in_rice (last access: March 2019),
2015.Zhang, G., Dong, J., Zhou, C., Xu, X., Wang, M., Ouyang, H., and Xiao, X.:
Increasing cropping intensity in response to climate warming in Tibetan
Plateau, China, Field Crop. Res., 142, 36–46,
10.1016/j.fcr.2012.11.021,
2013.