Nitrogen, deep phosphorus and potassium – a picture says a thousand words
Nitrogen, deep phosphorus and potassium – a picture says a thousand words
Author: Jayne Gentry (DPI QLD), Andrew Erbacher (DPI QLD), Henry Baskerville (DPI QLD), Megan Hunter (DPI QLD), Cameron Silburn (DPI QLD), David Lester (DPI QLD), Julie Boddington (DPI QLD), Isabella Macpherson (DPI QLD), Meg Tate (DPI QLD) | Date: 05 Mar 2025
Take home message
- Queensland grain growers identified managing nitrogen and declining soil fertility levels as their ‘trickiest’ (riskiest) decisions
- Crop responses to the application of nitrogen and deep phosphorus are variable
- Trials showed an increase in crop growth early in the season under deep phosphorus with increases in nutrient uptake
- Yield responses varied. An increase in yield of ~500 kg/ha was recorded on a soil that had very low levels of subsurface P when both additional phosphorus and nitrogen was applied
- Paddock variability drives crop responses. Understanding paddock variability can help inform research outcomes and may assist in improved management
- Large strip trials utilising proximal sensing and precision technology are helping understand crop responses encompassing paddock variability.
Background
RiskWi$e (the National Risk Management Initiative) is a 5-year national initiative aimed at helping growers understand and improve the risk-reward outcomes when making decisions in the face of uncertainty. To deliver RiskWi$e, a participatory action research methodology (an approach to research that pro-actively involves members of communities affected by that research in the research itself) has been employed. In Queensland the project is centred around Capella, Gindie, Roma, Brigalow, Meandarra, St George and Goondiwindi. Issue identification sessions were conducted for all groups, during which we explored participants ‘trickiest/riskiest’ decisions. Managing nitrogen (N) and declining soil fertility levels (N, phosphorus (P) and potassium (K)) were identified as being the riskiest. As a result, 30+ growers, with support from the project team have implemented on-farm trials investigating nitrogen management, declining soil fertility and managing soil constraints. This paper will explore three of these trials, situated at Wallumbilla, Surat and St Geroge.
Trial design and implementation
Treatments, trial designs and logistics were negotiated one-on-one with each of the growers and their agronomist. This type of research is not a one size fits all, rather purpose built for each grower, their research question, paddock, machinery and skill set. All paddocks were mapped using the online tool ConstraintID to quickly, easily and cheaply identify paddock variability to inform trial location. Paddocks were soil sampled, and trials implemented in 2023. The N treatment strips were implemented using grower machinery, with the design allowing the grower to plant and harvest. Deep P (& K) was applied at right angles to the N strips in two of the trials (Wallumbilla and Surat) in a strip plot design. At Wallumbilla deep P was applied using experimental machinery due to limited capacity of grower machinery, whereas it was applied by the grower at Surat. In contrast, the trial at St George was a split plot design with all treatments running parallel down the paddock (12 m wide “deep” treatments in 36 m “N” strips), all implemented by the grower. “Deep” treatments were ripped in at ~20 cm. There were no “rip only” control treatments as deep placed fertiliser logistically cannot be applied unless ripped in so any potential response to ripping is inherently part of the overall “deep” treatment response (Table 1). However, this interaction is being investigated in the small plot research of the GRDC NGN Deep P in SW Queensland (DAQ2404-008BGX). All treatments were replicated twice. It is envisaged that these trials will have N treatments reapplied each season for the remainder of the project, approximately four years.
Research questions
These trials are investigating the following questions:
- What are the comparative efficiencies and legacies of different nitrogen strategies (including rate and timing of application)? What are the risks and rewards of these approaches? What are the impacts on N losses, movement of N in the soil profile, efficiency of crop N use, yield and gross margin?
- Is the application of deep P (and K) effective and viable for increasing yield and profit compared across a range of nitrogen strategies? What are the risks and rewards and impacts on yield and cumulative gross margin?
Treatments
Table 1 outlines the treatment list for each site.
Table 1. Outline of nitrogen (N), deep phosphorus (P) and potassium (K) treatments at each of the three Queensland sites.
Yield potential (YP) is calculated annually on modelled yield.
Treatment name | Treatment description | Wallumbilla | Surat | St George |
---|---|---|---|---|
Nitrogen treatments | ||||
Nil | No fertiliser | ![]() | ![]() | ![]() |
Early N 50% YP | N budget targeting 50% yield potential, applied early in the fallow | ![]() | ![]() | ![]() |
Early N 75% YP | N budget targeting 75% yield potential applied early in the fallow | ![]() | ![]() | ![]() |
Early N 100% YP | N budget targeting 100% yield potential applied in the fallow | ![]() | ![]() | |
Late N 50% YP | N budget 50% yield potential applied at planting | ![]() | ||
Late N 75% YP | N budget 75% yield potential applied at planting | ![]() | ||
Deep treatments | ||||
Deep P | Deep P (150 kg/ha Granulock® Z; 30 kg P/ha) ripped in at a depth of ~20 cm | ![]() | ![]() | ![]() |
Deep P + K | Deep P + K (150 kg/ha Granulock® Z; 30 kg P/ha) plus (100 kg/ha Muriate of Potash; 50 kg K/ha) ripped in at a depth of ~20 cm | ![]() | ![]() |
Rates of nitrogen fertiliser are calculated annually via N budgeting based on modelled yield potentials in line with treatment rules and agreed to by the grower (Table 2). Due to sufficient soil N, additional N was not required on some treatments (note: only treatments with additional applied N will be presented). Nitrogen is not applied to pulse crops. All trials were planted to wheat in 2024. Approximate costs were calculated for each treatment (Tables 2 & 3).
Table 2. Target yield potential (YP) and quantity and cost of N fertiliser applied to wheat in 2024.
Site | Treatment | Yield potential wheat (t/ha) | N required* (kg/ha) | Avail. soil N** (kg/ha) | Applied N rate (kg/ha) | Applied urea rate (kg/ha) | Cost of urea applied*** ($/ha) |
---|---|---|---|---|---|---|---|
Wallumbilla | 50% YP | 2.3 | 82 | 60 | 23 | 50 | 40 |
75% YP | 3.2 | 115 | 60 | 55 | 120 | 96 | |
Surat
| 50% YP | 2.3 | 82 | 130 | 0 | 0 | 0 |
75% YP | 3.2 | 115 | 130 | 0 | 0 | 0 | |
100% YP | 6.0 | 215 | 130 | 85 | 184 | 147 | |
St George
| 50% YP | 2.2 | 93 | 100.5 | 0 | 0 | 0 |
75% YP | 3.0 | 126 | 100.5 | 23 | 50 | 40 | |
100% YP | 5.5 | 232 | 100.5 | 92 | 200 | 160 | |
Assumptions: *wheat @ 12% protein; **20–30 kg N/ha mineralisation over fallow in 0–90 cm soil profile; ***urea cost $800/t |
Table 3. Cost of Deep P and K application and yield improvement to break even per hectare for wheat grown in 2024.
Treatment | Ripper cost ($) | Granulock Z cost ($) | Muriate of Potash Cost ($) | Total cost ($/ha) | Breakeven yield increase required (t/ha) |
---|---|---|---|---|---|
Deep P | 75 | 180 | - | 255 | +0.85 |
Deep P & K | 75 | 180 | 70 | 325 | +1.08 |
Assumptions: Ripper cost $75/ha; Granulock® Z cost $1200/t; Muriate of Potash cost $700/t; wheat price $300/t |
Data collection
The following data is collected for specific treatments for growers’ crop rotations over the life of the trial:
- Soil mineral N and water at planting and harvest
- Biomass cuts during vegetative growth, at anthesis and physiological maturity
- Maturity biomass cuts are processed to determine grain yield, protein, quality and nutrient concentration
- Grain yield and protein from grower harvester yield & protein monitor
- Rainfall data
- Satellite imagery (NDVI).
Satellite imagery from the PlanetScope constellation was used to monitor crop growth and responses throughout the growing season. Vegetation indices, typically the normalised vegetation index (NDVI), measuring crop greenness, was used as an indicator of crop biomass and its distribution and variation across the trial. A key benefit remote sensing is the ability to capture transient patterns in the crop which may not be present when physical measurements are collected.
Electromagnetic (EM) surveys were used determine the apparent soil electrical conductivity (ECa). Several soil properties contribute to the ECa of a soil, including soil water, texture (especially the clay fraction) and ions in the soil (mainly sodium). Mapping ECa is a useful tool to identify changes in these soil properties and thus identify differing soil types within a paddock. This allows better interrogation of treatment responses, which may be affected by soil type.
Collection and analysis of harvester yield monitor data was used to assess treatment responses in this project and trial design was used to enable this methodology. Treatment strips run in line with traffic were applied at least one header width wide. Treatment strips running perpendicular to the traffic direction were wide enough to allow sufficient grain to allow yield monitoring to register discrete samples for each treatment, working within the limitations of the harvester size, processing and yield monitor systems. One header was used to harvest the trial area to avoid calibration issues between different headers. Raw data from the yield map was cleaned and analysed using statistical methods supported by AAGI, the GRDC’s Analytics for the Australian Grains Industry- however statistical results were not finalised at the time of writing this paper and as such all ‘results’ should be treated with caution. Headers allowed us to measure the entire plot in the trial (some are over 2 km long), as opposed to small individual locations when taking hand samples in the crop. These datasets provide insight into crop responses across the entire site and other factors which may be influencing plant growth.
Soil characterisation
The Wallumbilla, Surat and St. George sites were surveyed and characterised prior to implementing treatments (Table 4). The St George site had two distinct soil types in terms of clay content; hence they were soil sampled separately (Figure 1). Wallumbilla and Surat had very low levels of subsurface P while the St George site had acceptable levels. Current knowledge of critical limits suggests none of the sites are critically low in subsurface K.
Table 4. Soil characterisation for trial sites in 2023.
Site | Depth (cm) | Phosphorus Colwell (mg/kg) | Phosphorus BSES | Effective cation exchange capacity (ECEC (meq/100g) | Exc. potassium (meq/100g) |
---|---|---|---|---|---|
Wallumbilla | 0-10 | 17 | 24 | 23.3 | 0.63 |
10-30 | 5 | 5 | 28.9 | 0.34 | |
Surat | 0-10 | 12 | 26 | 41.4 | 0.67 |
10-30 | 4 | 18 | 42.9 | 0.41 | |
St. George - Clay | 0-10 | 28 | 204 | 38.1 | 1.49 |
10-30 | 14 | 215 | 38.7 | 1.12 | |
St. George - Sand | 0-10 | 52 | 82 | 7.5 | 0.55 |
10-30 | 57 | 168 | 9.7 | 0.48 |
Figure 1. St George trial site paddock (left) with soil types and treatments overlayed onto trial area (right).
Biomass – weight, nutrient content and NDVI
A combination of methods were deployed to monitor responses across all trials during the growing season. This was particularly useful, as satellite imagery assisted tracking transient crop responses whilst hand cuts provided quantification of biomass production and nutrient uptake.
Vegetative biomass and nutrient uptake
Biomass was sampled at the Wallumbilla and Surat trials at the vegetative stage (~8 weeks after planting) to assess the impact of N status, with and without deep P and K. At the Wallumbilla site, no responses to either N timing or application rate were apparent (Table 5). At the Surat site, there appeared to be more biomass, N and P uptake in the 100% YP treatment compared to the Nil treatment (Table 6). At this stage of growth, responses to deep P were evident at both sites, producing more biomass and increasing P uptake relative to the Nil treatment (Tables 5 and 6). The application of deep K also increased K uptake at the Surat site (Table 6).
Table 5. Wheat biomass and nutrient uptake at vegetative growth state at Wallumbilla in 2024.
Shoot biomass (t/ha) | N uptake (kg N/ha) | P uptake (kg P/ha) | ||||
---|---|---|---|---|---|---|
Deep P treatment | Nil | Deep P | Nil | Deep P | Nil | Deep P |
Nil | 1.5 | 1.8 | 53.4 | 67.6 | 3.5 | 4.6 |
Early N 50% YP | 1.1 | 1.4 | 45.2 | 62.1 | 2.1 | 3.6 |
Early N 75% YP | 1.3 | 1.4 | 52.3 | 59.8 | 2.6 | 4.0 |
Late N 50% YP | 1.3 | 1.3 | 53.3 | 51.4 | 3.0 | 3.3 |
Late N 75% YP | 1.2 | 1.5 | 49.3 | 58.8 | 2.5 | 4.0 |
*These data have not been statistically analysed and observations or trends should be treated with caution. |
Anthesis biomass and nutrient uptake
Biomass was sampled again at the Wallumbilla and Surat trials at anthesis to assess the impact of N status, with and without deep P and K.
Responses to N at the Wallumbilla site were variable, with no clear trends apparent regarding effects of either N application rate, or timing on biomass or N uptake (Table 7). At the Surat site, biomass responses to N application were less apparent at anthesis compared to vegetative growth stage, but the 100%YP treatment appeared to have a higher N uptake than the Nil. (Table 8).
At the Wallumbilla site, a response to deep P appears to be present at anthesis with respect to both biomass and P uptake (Table 7). Across all N treatments there is a trend observed around increased shoot biomass and increased N uptake with the application of deep P. In contrast, at the Surat site, the P response which appeared to be evident early in the season had diminished by anthesis (Table 8).
Maturity biomass and nutrient uptake
Biomass was sampled at maturity at all sites to assess the impact of nitrogen status, with and without deep P and K.
Finalising data from the laboratory is still in progress. However, preliminary data suggests that the trends apparent at anthesis followed through to maturity at all sites.
Table 6. Wheat biomass and nutrient uptake at vegetative growth state at Surat in 2024.
Shoot biomass (t/ha) | N uptake (kg N/ha) | P uptake (kg P/ha) | K uptake (kg K/ha) | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
Nil | Deep P | Deep P + K | Nil | Deep P | Deep P + K | Nil | Deep P | Deep P + K | Nil | Deep P | Deep P + K | |
Nil | 4.1 | 4.6 | 4.3 | 132.4 | 145.4 | 136.7 | 7.2 | 7.9 | 8.3 | 130.8 | 128.3 | 140.3 |
Early N 100 % YP | 4.3 | 4.6 | 4.8 | 143.3 | 153.1 | 167.8 | 7.9 | 9.5 | 9.7 | 148.6 | 138.2 | 162.7 |
Table 7. Wheat biomass and nutrient uptake at anthesis at Wallumbilla in 2024.
Shoot biomass (t/ha) | N uptake (kg N/ha) | P uptake (kg P/ha) | ||||
---|---|---|---|---|---|---|
Nil | Deep P | Nil | Deep P | Nil | Deep P | |
Nil N | 7.3 | 8.1 | 100.5 | 111.9 | 13.2 | 14.6 |
Early N 50% YP | 6.9 | 8.4 | 83.9 | 120.8 | 10.4 | 11.7 |
Early N 75% YP | 7.6 | 8.5 | 112.4 | 120.7 | 11.4 | 13.5 |
Late N 50% YP | 7.7 | 9.2 | 99.6 | 131.4 | 10.8 | 13.8 |
Late N 75% YP | 7.3 | 8.3 | 96.0 | 116.2 | 10.9 | 12.8 |
Table 8. Wheat biomass and nutrient uptake at anthesis at Surat in 2024.
Shoot biomass (t/ha) | N uptake (kg N/ha) | P uptake (kg P/ha) | K uptake (kg K/ha) | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
Nil | Deep P | Deep P + K | Nil | Deep P | Deep P + K | Nil | Deep P | Deep P + K | Nil | Deep P | Deep P + K | |
Nil | 10.9 | 10.8 | 10.3 | 158.8 | 158.0 | 156.2 | 15.3 | 13.5 | 14.5 | 212.4 | 167.7 | 179.2 |
Early N 100 % YP | 11.1 | 10.5 | 10.8 | 195.8 | 173.2 | 171.2 | 17.7 | 16.8 | 14.6 | 242.8 | 185.1 | 179.3 |
*Data has not been statistically analysed and observations or trends should be treated with caution. |
NDVI
NDVI from satellite imagery was captured at regular intervals during the growing season. The Surat trial was interesting to follow (data not shown). Large differences were seen in the deep P treatment strips in the first months of the trial, supported by biomass data (Table 6). However, these distinct differences slowly decreased as the crop matured, with no yield responses evident by harvest (Figure 6). This suggests that although the deep P supported additional early crop growth, this benefit was not translated into yield. This response was most likely driven by timing of rainfall.
Yield and protein
Yield and protein results were also collected using a variety of methods. Header yield monitor data was captured for all three sites and in addition, protein mapping data from Surat, with the raw data cleaned, ready for statistical analysis. Hand harvest data was also collected to ground truth yield maps (data not shown). Grain from these samples was analysed for protein.
Wallumbilla
This trial began showing quite severe crown rot symptoms which would have impacted yield responses. The harvester yield map (Figure 2) suggests an increase in yield can be observed the deep P treatments, pending statistical analysis. However, there did not appear to be any response to N treatments in the absence of P. The addition of deep P appeared to increase yields across all N treatments (Figure 3). This site had the lowest levels of P in the subsurface (similar Colwell and lower BSES P) of the three sites (Table 4). There appears to be a positive yield response to Deep P alone (Deep P Nil N treatment) of ~200 kg/ha, however this response appeared to benefit from the application of N with higher responses to P in the presence of early and late applied N at the 75% YP level (Early N 75% YP & Late N 75% YP) of ~500 kg/ha. This response supports the theory that sufficient N is required to maximise a response to deep P. The cost of this treatment was ~$350/ha and is in a basic gross margin deficit of ~$175/ha after year 1. Where there was no deep P applied there appeared to be no yield response to N. Nitrogen increased grain protein levels, with a trend observed around further increases from the application of deep P (Figure 4). It will be interesting to continue to follow crop responses over time as it is anticipated that yield benefits from deep P may be observed for many years to come. N treatments have been designed as strategies as opposed to set quantities, and annual budgets will be calculated prior to each crop, hence carry over N will be incorporated into future N budget requirements. This paddock has been double cropped into sorghum and will results will continue to be collected.
Figure 2. Yield map – Wallumbilla.
Figure 3. Grain yield - Wallumbilla (yield map).
*Data has not been statistically analysed and observations or trends should be treated with caution.
Figure 4. Grain protein - Wallumbilla (NIR data).
*Data has not been statistically analysed and observations or trends should be treated with caution.
Surat
Additional N was only applied to the Early N 100% YP treatment. Surat yield responses were interesting. This site had an obvious response to deep P early in the season, however, as the NDVI showed, this visual response was lost as the crop matured. The yield map (Figure 5) shows the variability within the paddock. Imagery indicates that the deep P and deep P + K strips at the top of the paddock (rows B and C, bottom of page) were lower yielding than the strips at the other end of the paddock (E and F, top of page). Since these replicates were averaged, we suggest that paddock variability is potentially skewing results. The data extracted from this map indicates that the Nil treatment was the highest yielding (Figure 6). Paddock variability is likely the main driver of differences in the data however, we are interrogating this data further with assistance from the AAGI team. The soil test results indicate that this soil may not be limited in K in the subsurface (Table 4). The application of N increased grain protein (Figure 6). Again, this is only one season of data from a long-term trial, results will be collected in subsequent years.
Figure 5. Protein (left) and yield (right) maps – Surat.
Note: gaps in the protein map are a result of missing data.
Figure 6. Grain yield and protein – Surat (yield monitor).
*Data has not been statistically analysed and observations or trends should be treated with caution.
St George
Additional N was required to be applied to meet N treatment level requirements for the Early N 75% and 100% YP treatments whilst soil N was sufficient in other treatments. The St George site has two distinct soil types within the trial (Figure 1). When considering the differences in cation exchange capacity and nutrient availability of different soil textures and types, it was expected that observations of treatment effects would differ for soil types. This is an example of how averaging across this paddock variability rather than taking variable soils into account can dramatically impact results. Consequently, when yield monitor data (Figure 7) was interrogated with clay and sandy soil separated very different trends were observed. The clay soil, on average, yielded ~0.7 t/ha more than the sand (not statistically analysed). The N 100% YP treatment appeared to increase yield by ~200 kg/ha in the clay above the nil, but decreased yield on the sandier soil by approximately the same amount (Figure 8). There was no obvious response to deep P and K in the clay. However, there was an apparent increase in yield from deep P in the sand, the highest being ~300 kg/ha in the N 100% YP + deep P treatment. This was interesting as soil test results indicated the sand had higher levels of P in the subsurface compared to the clay (Table 4). This trial has now been planted to mungbeans 2025 and will continue to be monitored.
Figure 7. Yield map – St George.
Figure 8. Grain yield, clay and sandy soils - St. George (yield monitor).
Conclusion
Preliminary results show a range of responses both to nitrogen and deep P with statistical analysis yet to be completed. None of the trials indicated a response to deep P + K compared with deep P alone, noting results are yet to be analysed statistically. The Wallumbilla site, with the lowest levels of subsurface P (BSES), appeared to show the greatest response to deep P. Increases were seen in biomass production, nutrient uptake and yield, with the greatest yield observed under the highest N treatment (75% YP) with deep P of ~500 kg/ha. What was interesting was the general lack of response to N. In almost all trials there was no economic response to N, noting data is still being collated, and statistical analysis is yet to be completed.
The value of these trials is twofold. Firstly, the growers and their agronomists have been integral in identifying the research questions and determining treatments so that each site is customised to suit their situation. Secondly, is their scale. The paddock scale has allowed the project team to deploy new research methodologies using satellite imagery and EM mapping combined with commercial scale equipment including yield and protein monitors. As a result, we have been able to easily visualise and account for paddock variation. Data will continue to be collated which will then be shared with growers as part of the RiskWi$e initiative. Part of this process will be determining the risks and rewards of these various nutrition strategies. These trials will continue to be monitored over the coming years, which will provide more informed insights.
Acknowledgements
The research undertaken as part of this project is made possible by the significant contributions of growers through both trial cooperation and the support of the GRDC, the author would like to thank them for their continued support.
RiskWi$e is conducted in partnership with grower groups, 6 action research groups, research/extension partners, CSIRO and the Grains Research and Development Corporation (GRDC). The Queensland Department of Primary Industries is leading the RiskWi$e initiative on various topics across Queensland in collaboration with JB Ag Services and Iker Ag Consulting.
Contact details
Jayne Gentry
Queensland Department of Primary Industries
Toowoomba, Queensland
jayne.gentry@daf.qld.gov.au
Date published
March 2025
GRDC Project Code: DAQ2303-011BGX, DAQ2404-008BGX,