Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

The purpose of this document is to inform potential end-users about the functionality of the about the functionality of the APSIM Wheat Yield Data Yield Data Cube. Here, we describe the methodology, including the source data, APSIM crop model configurations and outputs available configurations and outputs available to end-users via an APIan API.  See Agriyieldz API Guide for details on how to access the product via the API.

Table of Contents

Background Background 

CSIRO have created a unique data product about wheat crop yields for Australia. This data product can provide insight into the wheat production potential, the variability of production and the riskiness of production potential. Data are available at the field scale, where fields, farms, or regions can be compared for productivity and variability in crop yield. Data have been created using the APSIM wheat crop model, a soil grid, and climate data from the SILO patched point data base. Data can be exportedexported, via an API, to end-users for reprocessing and integrating into other data platforms. 

Methodology Methodology 

The APSIM Next Gen Crop model (Holzworth et al., 2018) was used to generate simulations of wheat crop yield across the grain growing belt of Australia in every year from 1990 to 2020. Simulations were created for each possible wheat crop field in Australia. CSIRO have mapped 1.7 million fields across the continent, and these simulations are designed to augment those field boundaries. These data allow for example, an individual or organisation to make field to field, farm to farm and region to region comparisons relating to production potential, and the variation in production.    

Locations

At each location At each location (i.e. the centroid of a  the centroid of a field), the functional soil  the functional soil was selected from a the ASRIS soil map. The soils were selected from 18 possible candidates, that represent the functional properties of soils across Australia (Searle et Australia (Searle et al., 2021 2021). These soils are defined defined with sufficient complexity to provide insight into how complexity to provide insight into how crops and pastures will grow under grow under different climates. They  They are not a literal translation of soil type, as defined by the Australian Soil Classification system. Rather, they are  Rather, they are a subset of the soils from the APSOIL database. The  The complete list of soils used used to create the simulations are provided in Table 2.  Simulations in Table 2.  Simulations were only executed for the for the soil at the exact centroid exact centroid location of the field, and simulations were only executed for sensible candidate soils for the location. For example, there are no peaty soils (51) in Western Australia, and simulations were not run for that soil in WA.


Table 1. Number of simulations performed across Australia to create the APSIM wheat data cube


Unique Locations

Total Simulations

NSW

1488

6792

QLD

1012

5385

SA

665

3442

VIC

1129

4041

WA

795

2915


Table 2. List of soils used to create the APSIM wheat data cube (Location of original soil profile in parentheses).

Apsoil Number

Apsoil Description

179

Brown Chromosol (Temora)

175

Red Chromosol (Coolamon)

1194

Silty Clay Loam over Light Clay (Goonumbla)

1197

Sandy Clay Loam over Sandy Clay over Light Medium Clay (Tullibigeal)

878

Hydrosol (Babinda)

1211

Sandy earths (Yandanooka )

192

Sandy Loam over Clay Loam (Parkes)

650

Redoxic Hydrosol (Fairymead Mill )

424

Yellow Deep Sand (Buntine)

1209

Coloured Sand (Allanooka )

419

Shallow Gravel (Buntine)

1108

Shallow loamy duplex (Salmon Gums)

300

Sandy loam over light-medium clays (Roseworthy)

1014

Grey Vertosol (Pilliga)

1272

Black Vertosol (Capella)

1267

Black Vertosol (Fernlees)

1168

Brown Vertosol (Bundella)

1255

Shallow dark clay loam on limestone (Padthaway)

Management

To create the data cube, general crop management rules were created to grow a crop. Regardless of location, or season, sowing took place on the 15th Maythe 15th May, with the crop emerging 5 days later. The broadly adapted wheat cultivar, Wyalkatchem, was selected. The objective was to grow a crop to its water limited yield potential, where yields were not limited by a lack of nutrients. To achieve yield potential, initial total soil nitrogen was 300 kg N/ha. Initial water, at January 1, was 10% of capacity, and initial residues were 1000 kg/ha wheat stubble. 

Data

For each location, in each year from 1990 to 2020, data are available. The primary output is wheat yield, but additional variables relating to soil water, growing season rainfall and rainfall and occurrence of extreme events such as frost, and heat stress are stress are available (Table 3). The outputs from  The outputs from the 30-year simulations of wheat are stored in a database. Geographic co-ordinates are stored, and it is possible to interrogate the data for a single year and for multiple years. Comparisons  Comparisons are possible between fields, regions and regions and states. 

Table 3. Outputs, available from 1990 to 2020 for every field.    

Wheat Grain Weight (g/m2)

  

  

Grain

Yield 

Yield 

Wheat Above Ground Weight (g/m2)

   

   

Total Wheat

Biomass  

Biomass  

Plant available water at Sowing (mm)

  

  

Starting soil water status, at

sowing   

sowing   

Plant available water at Harvest (mm) 

Soil water at the end of the season, at

harvest 

harvest 

In Crop Rain (mm) 

Total amount of rainfall between sowing and

harvest  

harvest  

Plant available water capacity (mm) 

Amount of water that can be held in the

soil 

soil 

Mild Frost Count (days) 

Number of days where the temperature falls

between 0◦ C and – 2◦ C during flowering 

between 0◦ C and – 2◦ C during flowering 

Moderate Frost Count (days) 

Number of days where the temperature falls between -

2◦ C

2◦ C and –

4◦ C during flowering 

4◦ C during flowering 

Severe Frost Count (days) 

Number of days where the temperature drops below -

4◦ C during flowering 

4◦ C during flowering 

Mild Heat Count (days) 

Number of days where the temperature rises

between 32◦ C and 34◦ C

between 32◦ C and 34◦ C, from

flowering  

flowering  

Moderate Heat Count (days) 

Number of days where the temperature rises between

34◦ C and 36◦ C

34◦ C and 36◦ C, from

flowering 

flowering 

Severe Heat Count (days) 

Number of days where the temperature exceeds

36◦ C

36◦ C, from

flowering  

flowering  

Sense Testing and Data screening

Outputs from the model were compared to a simple French- Schultz modelSchultz model, that predicts yield in response to growing season rainfall and starting soil moisture.  Additional outputs about temperature extremes for heat and frost were extracted from the data. The goal was to create a simple evaluation that allowed the user to quickly gauge if the model is performing as expected for expected for the desired location. APSIM is more sophisticated than the than the French-Schultz modelSchultz model, so perfect correlations are not expected. However However, over a 30 year simulationa 30 year simulation, at one location, the mean error the mean error should be less than 20%.  Maps   Maps of the mean error, between the French-Schultz model Schultz model and APSIM are created.    

Results

APSIM output, for WA, SA, Vic, NSW and Qld are presented in Figures one to five. The figures illustrate the spatial distribution of long term (30 year) potential yield provided by APSIM.  The mean deviation at a single location, from a French-Schultz equation for WA, SA and Vic are generally low, and average less than 20%. Within a year, the deviation can be higher, as the APSIM model is more complex, and accounts for a greater range of abiotic stresses than the French-Shultz equation. For example, In NSW and Qld, the mean deviation is higher, in part because the French-Shultz model does not summarise potential yield as well, where in-season rainfalls are often low.

...

Output from the APSIM wheat data cube can be used to compare locations. In Figure 11,  the yields over 30 years, from two paddocks in Western Australia are compared. The paddock with low yield originates from the low rainfall zone (-31.5936 ,  118.2347). This particular field had a long term mean of 2.45 t/ha, with a standard deviation of 0.62 t/ha and a coefficient of variation of 26%. The other paddock originates from the high rainfall zone (-34.5714, 117.01). The long-term mean, standard deviation and coefficient of variation for this field was 3.52t/ha, .84 t/ha and 24%.  Data can be extracted from as many fields as needed, and where necessary aggregated to create regional assessments.Image Removed


...

Figure 11. Crop yields, generated by the APSIM crop model from 1990 to 2020 for fields located in the high rainfall zone and low rainfall zone of Western Australia.

...