In the current study, we combined additional in situ 15N2O measurements together with high temporal resolution measurements of N2O presented in Takeda et al. (2021a) and the recovery of 15N-labelled fertiliser presented in Takeda et al. (2021b) conducted on a commercial sugarcane farm in Burdekin, QLD (19° 37' 4'' S, 147° 20' 4'' E) from October 2018 to August 2019 to attribute N2O losses to fertiliser and soil (non-fertiliser) N sources. We conducted a complementary field experiment on a commercial sugarcane farm in Mackay, QLD (21° 14' 4'' S, 149° 04' 6'' E) from October 2019 to August 2020 and present the results together.
Study site
The climate is tropical in both Burdekin and Mackay and the average annual rainfall is 945 mm (Bureau of Meteorology, Site 033002, Ayr) and 1598 mm (Bureau of Meteorology, Site 033119, Mackay), respectively. The annual average daily minimum and maximum temperatures are 18.0 and 29.2°C in Burdekin and 19.1 and 26.5°C in Mackay. The soil is classified as Brown Dermosol and Brown Kandosol in the Australian Soil Classification (Isbell 2016), or Luvisol and Fluvisol in the World Reference Base (WRB) Classification (IUSS Working Group 2014) at the Burdekin and Mackay sites, respectively. Sugarcane varieties Q240 and Q208 were planted in 2015 and 2016 and the crop was the 3rd ratoon during the experiment at the Burdekin and Mackay sites, respectively. Irrigation was applied by furrow irrigation at the Burdekin site and overhead sprinkler at the Mackay site. Sugarcane was burnt before harvest at the Burdekin site. ‘Green cane trash blanketing (GCTB)’, a practice where the cane is harvested green and the trash (crop residues) is spread over the ground, is practised at the Mackay site. Site conditions are summarised in Table 1.
Table 1
Climate conditions, soil properties at 0-0.2 m depth and crop management practices at the Burdekin and Mackay sites
Category | Variable | Burdekin | Mackay |
Experiment | Coordinates | 19° 37' 4'' S, 147° 20' 4'' E | 21° 14' 4'' S, 149° 04' 6'' E |
| Period | Oct. 2018 to Aug. 2019 | Oct. 2019 to Aug. 2020 |
Climate | Annual mean rainfall (mm) | 945 | 1598 |
| Cumulative rainfall over the experimental period (mm) | 1180 | 1418 |
| Annual mean daily minimum temperature (°C) | 18.0 | 19.1 |
| Annual mean daily maximum temperature (°C) | 29.2 | 26.5 |
Soil (0.0-0.2 m) | BD (g cm−3) | 1.3 | 1.1 |
| pH (H2O) | 6.92 | 4.13 |
| Total C (%) | 1.60 | 1.35 |
| Total N (%) | 0.08 | 0.09 |
| Clay (%) | 10.5 | 22.2 |
| Silt (%) | 48.5 | 15.9 |
| Sand (%) | 41.0 | 61.9 |
| Mineral N (kg N ha−1) | 37.0 | 31.8 |
Crop management | Cultivar | Q240 | Q208 |
| Crop | 3rd ratoon | 3rd ratoon |
| N rate treatments (kg N ha−1) | 0, 150, 200 and 250 | 0, 100, 150, 200 and 250 |
| SIX EASY STEPS N rate (kg N ha−1) | 200 | 150 |
| N fertiliser product | Urea | Urea |
| Fertiliser application | Two-sided banding | Stool splitting |
| Irrigation management | Furrow | Overhead |
| Estimated total irrigation amount (mm) | 600 | 180 |
| Trash management | Burnt | GCTB |
Experimental design
The experimental design of the Burdekin site was a randomised strip design with four plots across two strips for each N treatment, detailed in (Takeda et al. 2021a). The experiment at the Mackay site had a completely randomised block design with three replicates per treatment (around 10 m in length) and three planting rows per plot with a row distance of 1.56 m, accompanied by an unfertilised control (0N) plot with three subplots (20 m of length). Fertiliser N rate treatments included 0N, 150 kg N ha-1 (150N), 200 kg N ha-1 (200N) and 250 kg N ha-1 (250N), plus 100 kg N ha-1 (100N) at the Mackay site only. The recommended N rate in SIX EASY STEPS best management practice of the Australian sugar industry (Schroeder et al. 2010) was 150N at the Mackay site and 200N at the Burdekin site. Urea was applied by banding the fertiliser 10 cm deep and 30 cm from the bed centre on both sides of the cane row at the Burdekin site and by stool splitting 10 cm deep at the bed centre of the cane row at the Mackay site. For the 15N recovery in the soil and the plant, a 2.0 m section was exempted from the application of N fertiliser in each plot and 15N enriched urea fertiliser (5 atom%) in solution was manually applied at the corresponding rate (i.e. 100N, 150N, 200N and 250N) matching the N fertiliser placement at the respective site.
Measurement of N 2 O and CO 2 emissions using an automated chamber system
Soil-borne N2O and CO2 emissions were measured at a high temporal resolution using an automated chamber system (Grace et al. 2020) from 17 October 2018 to 15 August 2019 at the Burdekin site and from 3 October 2019 to 24 August 2020 at the Mackay site. Details of the automated chamber system are given in Supplementary materials S1.1. Manual gas sampling was conducted for the control plots of the Mackay site by the static closed chamber method (Friedl et al. 2017), detailed in Supplementary materials S1.2. The placement of the chambers accounted for N fertiliser placement and irrigation practice at each site: At the Burdekin site, chambers were installed covering the area from (a) the fertiliser band to the centre of the bed (bed chamber) and (b) the fertiliser band the centre of the furrow (furrow chamber). At the Mackay site, bed chambers (a) were placed at the centre of the bed (i.e., on the fertiliser band) and furrow chamber measurements (b) were substituted with those from the control plots.
Measurement of fertiliser and soil-derived N 2 O emissions using a 15N method
Emissions of fertiliser-derived N2O were measured as described by (Friedl et al. 2017). Micro plots were established alongside the main plots with N fertiliser rates of 150, 200 and 250 kg N ha-1 at the Burdekin site and with 100, 150, 200 and 250 kg N ha-1 at the Mackay site. The micro plots had a completely randomised block design with four replicates. A steel base (0.22 m × 0.22 m at the Burdekin site and 0.2 m × 0.4 m at the Mackay site) was installed in each micro plot and 15N enriched urea fertiliser (70 atom%) was applied inside the base at the corresponding rates. Gas sampling was conducted with static closed chambers at the Burdekin site from November 2018 to February 2019 and with semi-automated chambers at the Mackay site from October 2019 to January 2020, detailed in Supplementary materials S1.3. The gas samples were analysed for the concentration of N2O using a Shimadzu GC-2014 Gas Chromatograph (Shimadzu, Kyoto, Japan) and for 15N2O using an Isotope Ratio Mass Spectrometer (IRMS) (20–22 Sercon Limited, UK).
Plant and soil sampling and analyses
Plant and soil samples were taken from each of the 2.0 m sections just before the farmer’s harvest (on 27–28 August 2019 at the Burdekin site and 25–26 August 2020 at the Mackay site). The procedure of plant and soil sampling and analyses at the Burdekin site are detailed in Takeda et al. (2021b) and briefly described here. Plant sampling procedures are detailed in Supplementary materials S1.4. Soil samples were taken at three to four points between the bed and furrow centres using a soil corer and a post-hole driver down to 1.0 m. At the Burdekin site, the sampling points were 0, 0.25, 0.50, 0.75 m away from the bed centre. At the Mackay site, those were 0, 0.12, 0.40, 0.80 m away from the bed centre in 100N and 250N and the three points except for 0.12 m away from the bed centre in 150N and 200N. Each soil core sample of 1.0 m was separated into 0–0.2, 0.2–0.4, 0.4–0.7, 0.7–1.0 m soil depths and subsamples (about 100 g) were taken. The dried plant and soil samples were then finely ground and analysed for N and 15N content via IRMS analysis (20–22 Sercon Limited, UK). The stalk samples were further dried with a vacuum oven at 40°C for 48 hours before fine grinding to avoid aggregation due to sugar.
15 N calculations
Fertiliser 15N recovery (%) in the plant and soil was then calculated by 15N mass balance (Friedl et al. 2017; Rowlings et al. 2016). The percentage of N in individual soil and plant material samples (‘sinks’) derived from 15N-labelled fertiliser (Ndff plant or soil) was calculated from
\(Ndff \left(\%\right)=\frac{{\%}^{15}N excess of sink}{{\%}^{15}N excess of fertiliser}\times 100\) [1]
where the %15N excess used for all sources and sinks was the 15N abundance less an adjustment of %15N measured for the corresponding samples in the 0N plots for background enrichment. Calculations for fertiliser 15N recovery in plant and soil are detailed in Takeda et al. (2021b) and given in Supplementary materials S1.5. For the gas samples from the 15N micro plots, the measured N2O %15N was corrected to account for the N2O existing in the headspace at 0 min after chamber closure, detailed in Supplementary materials S1.6. After the correction, the proportion of N2O emissions derived from fertiliser (Ndff N2O) at the fertiliser band was calculated by Equation [1]. Then, Ndff N2O was gap-filled per N rate over the crop growing season on a daily basis at each site. The gap-filled daily Ndff N2O per N rate from the micro plots was then applied to daily N2O emissions per plot measured by the automated chamber system to calculate fertiliser-derived N2O emissions as follows:
\({Fertiliser derived N}_{2}O {emissions}_{i, j}= {area weighted N}_{2}O {emissions}_{i, j}\times \frac{{Ndff {N}_{2}O}_{i, j}}{100}\) [2]
where i and j indicate days after fertilisation and chamber position (i.e. bed or furrow). At the Burdekin site, the adjusted Ndff N2O at the fertiliser band in micro plots was used for both bed and furrow chambers because both chambers covered the fertiliser band and N2O emissions did not differ between the positions. At the Mackay site, adjusted Ndff N2O at the furrow was assumed to be zero because the furrow chamber did not cover the fertiliser band. The rationale behind these assumptions is detailed in Supplementary materials S1.7. Soil-derived N2O emissions (i.e. non-fertiliser-derived N2O emissions) was calculated by the difference between total N2O emissions and fertiliser-derived N2O emissions per plot on a daily basis.
Auxiliary measurements
For soil ammonium (NH4+) and nitrate (NO3−) measurements, soil samples (0–20 cm depth) were taken near the fertiliser band in each plot one day after fertilisation, every 3–7 days for the first three months and then monthly thereafter. Soil NH4+ and NO3− were extracted by adding 100 mL of 2 M KCl to 20 g of air-dried soil and shaking the solution for one hour, followed by NH4+ and NO3− content measurements using a Gallery™ Discrete Analyzer (Thermo Fisher Scientific, USA).
Statistical analysis
Statistical analyses and graphical presentations in this study were conducted using R statistical software version 3.5.2 (R Core Team 2018) with a significant level set at p < 0.05. Daily N2O emissions were calculated by averaging the measured hourly fluxes over a 24-h period from each chamber and multiplying by 24. Missing daily N2O emissions between measurements were imputed by linear interpolation, using a package “imputeTS” (Moritz and Bartz-Beielstein 2017). Area-weighted N2O emissions were calculated by the ratio bed:furrow = 1:1 at the Burdekin site, while bed:furrow = 1:2 at the Mackay site. Cumulative N2O emissions were then calculated by summing the daily average of each chamber over each growing period. The 2-D interpolation of the fertiliser 15N remaining in the soil from the bed centre to the furrow centre down to 1.0 m depth at harvest was conducted for each plot by fitting a thin-plate spline using a package “mgcv” (Wood 2011). Gap-filling of Ndff N2O values between measurements over the crop growing season was conducted by fitting generalised additive mixed models (GAMMs), using a package “mgcv” (Wood 2011) and detailed in Supplementary materials S1.8. Briefly, GAMMs can quantify non-linear relationships without specifying the functional forms (De Rosa et al. 2020; Dorich et al. 2020), which were used to analyse Ndff N2O over time in response to days after fertilisation (DAF) and N rates. Two-way analysis of variance (ANOVA) was performed to assess significant differences in cumulative (fertiliser and soil-derived) N2O emissions, plant N uptake and fertiliser 15N recoveries between the N rate treatments and the sites. Subsequently, the Tukey-Kramer test was performed when the ANOVA yielded P values < 0.05, using a package “agricolae” (de Mendiburu 2020). To evaluate the response of fertiliser 15N loss to N rates across sites, (generalised) linear mixed models were fitted with the site as a random factor, using a package “lme4” (Bates et al. 2015).