Greenhouse Gas Emissions and Crop Yields From Winter Oilseed Rape Cropping Systems are Unaffected by Management Practices

Winter oilseed rape is traditionally established via plough-based soil cultivation and conventional sowing methods. Whilst there is potential to adopt lower cost, and less intensive establishment systems, the impact of these on greenhouse gas emissions have not been evaluated. To address this, field experiments were conducted in 2014/2015 and 2015/2016 to investigate the effects of 1) crop establishment method and 2) sowing method on soil greenhouse gas emissions from a winter oilseed rape crop grown in Ireland. Soil carbon dioxide, nitrous oxide and methane emission measurements were carried out using the static chamber method. Yield (t seed ha−1) and the yield-scaled global warming potential (kg CO2-eq. kg−1 seed) were also determined for each management practice. During crop establishment, conventional tillage induced an initially rapid loss of carbon dioxide (2.34 g C m−2 hr−1) compared to strip tillage (0.94 g C m−2 hr−1) or minimum tillage (0.16 g C m−2 hr−1) (p < 0.05), although this decreased to background values within a few hours. In the crop establishment trial, the cumulative greenhouse gas emissions were, apart from methane, unaffected by tillage management when sown at a conventional (125 mm) or wide (600 mm) row spacing. In the sowing method trial, cumulative carbon dioxide emissions were also 21% higher when plants were sown at 10 seeds m−2 compared to 60 seeds m−2 (p < 0.05). Row spacing width (125 and 750 mm) and variety (conventional and semi-dwarf) were found to have little effect on greenhouse gas emissions and differences in seed yield between the sowing treatments were small. Overall, management practices had no consistent effect on soil greenhouse gas emissions and modifications in seed yield per plant countered differences in planting density.


INTRODUCTION
Atmospheric concentrations of the greenhouse gases (GHG) carbon dioxide (CO 2 ), nitrous oxide (N 2 O) and methane (CH 4 ) continue to rise globally due to anthropogenic activities. Agriculture is responsible for approximately 12% of global GHG emissions with livestock systems, soil cultivation, rice production and crop residue management making the more significant contributions (Ciais et al., 2014). In terms of land use impacts, croplands make one of the major contributions to agricultural GHG emissions through various farming and management activities. Field operations such as soil tillage, sowing, fertilizer addition and chemical treatment that are normally required to maximize plant productivity and yield can have significant impacts on GHG emissions through perturbations in the carbon (C), nitrogen (N) and water dynamics of these ecosystem (Bondeau et al., 2007;Osborne et al., 2010). Mitigation of GHG emissions from agricultural sources involves measures that aim to increase soil organic carbon (SOC) sequestration, reduce GHG emission rates, or both, through improved management practices Minasny et al., 2017;Ogle et al., 2019). Cropland management practices that involve the adoption of less intensive soil cultivation and agronomic interventions that improve agricultural productivity may can lead to reductions in GHG emissions and determine whether these ecosystems function as sinks or sources of C (Ceschia et al., 2010).
Soil CO 2 is produced by autotrophic and heterotrophic processes. Autotrophic respiration is derived from roots and/ or the metabolism of photosynthetic substrates released by roots (Högberg and Read, 2006) whereas heterotrophic respiration is associated with the microbial decomposition of SOC or root exudates (Trumbore, 2000). Agriculture accounts for ∼60% of the global anthropogenic N 2 O emissions primarily due to increased N-fertilizer use (Smith, 2017), with croplands accounting for at least 80% these emissions (Tian et al., 2019). Nitrous oxide production in soils arise from nitrification; the conversion of ammonium (NH 4 + ) to nitrite (NO 2 − ) and subsequently nitrate (NO 3 − ) in aerobic soils, and denitrification, the sequential reduction of NO 3 − to gaseous N 2 O or N 2 in anaerobic soils (Butterbach-Bahl et al., 2013). The soil-atmosphere exchange of CH 4 is largely governed by the balance between net CH 4 production by methanogens and net CH 4 oxidation/uptake by methanotrophs (Serrano-Silva et al., 2014). Although upland arable soils are generally aerobic and considered CH 4 sinks, intensive soil management practices have reduced the capacity of many soils to oxidise CH 4 (Suwanwaree and Robertson, 2005). Land management practices thus often govern whether cropland soils are net sinks or sources of GHGs (Ceschia et al., 2010).
Tillage and sowing are management practices that influence crop establishment, plant growth, nutrient uptake, canopy/soil microclimate and yield (Sharratt and McWilliams, 2005;Malhi et al., 2006;Soane et al., 2012), and subsequently impact on C and N and the GHG balance in croplands (Ceschia et al., 2010;Moors et al., 2010). Conventional tillage (CT) encompasses soil inversion, and crop residue incorporation, facilitating seedbed preparation for the succeeding crop. However, CT practices result in the mechanical disruption of soil aggregates and the release of protected organic C from soil organic matter (SOM) resulting in enhanced CO 2 emissions. Altering the turnover rate of SOM can directly impact C sequestration and the emissions of CO 2 , N 2 O and CH 4 (Six et al., 2004;Abdalla et al., 2013;Abdalla et al., 2016;Shakoor et al., 2021). Agricultural practices that minimize or reduce tillage operations can increase soil aggregate formation potentially retaining ∼30% of crop residues on the soil surface increasing nutrient availability and reducing soil erosion (CTIC, 2004). These practices include minimum tillage (MT), involving shallow cultivation to a depth of 5-10 cm, and strip tillage (ST), combining cultivated strips (25 cm depth) with direct drilling whilst leaving the interrow spaces unaffected, have been promoted as alternative management practices to CT (Davies and Finney, 2002;Morris et al., 2010) that could reduce GHG emissions and increase C sequestration. In combination with these approaches crop residues may also be retained on the soil surface to protect the C in soil aggregates. The argument being that these non-inversion tillage systems will reduce the exposure and subsequent oxidation of SOC and lower GHG emissions.
In addition to CO 2 , soil tillage can also affect the emissions of N 2 O and CH 4 . Both N 2 O and CH 4 have global warming potentials (GWP) that are 265 and 28 times higher than CO 2 , respectively (Myhre et al., 2013). In some cases, reductions in CO 2 emissions that have been attributed to a particular tillage management may be off-set by increases in the emissions of N 2 O (Six et al., 2004). The rate of N 2 O production is controlled by factors that are affected by tillage intensity, such as the soil waterfilled pore space (WFPS), soil organic C availability and temperature (Butterbach-Bahl et al., 2013). The reported effects of RT operations on N 2 O emissions are equivocal, showing either increases (Ball et al., 1999;Shakoor et al., 2021), decreases (Kessavalou et al., 1998;Chatskikh and Olesen, 2007;Ussiri et al., 2009) or similar N 2 O emissions (Abdalla et al., 2010;O'Neill et al., 2020).
Methane fluxes are regulated by factors such as soil moisture, temperature, oxygen availability, SOM content and C and N availability (Jacinthe and Lal, 2005). Therefore, tillage practices, through their effects on soil physicochemical properties, can influence net CH 4 emissions (Hütsch, 1998). Well-aerated soils with higher oxygen availability that facilitate CH 4 diffusion to methanotrophic sites tend to have higher rates of CH 4 uptake than poorly drained soils that restrict CH 4 oxidation (Prajapati and Jacinthe, 2014). Soils may also differ in their capacity to oxidise methane, however, the mechanism(s) involved is not fully understood (Lang et al., 2020).
Oilseed rape (Brassica napus L.; OSR) is the third most important oil crop in the world after oil palm and soybean (FAOSTAT, 2021). The winter OSR variety is predominantly cultivated in temperate climatic regions of Europe, where the mean yield of winter and spring OSR varieties (2.7 t ha −1 ) consistently exceed the global average (2.1 t ha −1 ) (FAOSTAT, 2021). In Ireland, WOSR is cultivated as a break crop in cereal production systems on ∼10,000 ha annually (Zahoor et al., 2015;CSO, 2020) and traditionally established via CT sowing operations (Forristal and Murphy, 2010). Compared to CT practices, the crop area under reduced tillage in Ireland is small at 40,000 ha (Meade and Mullins, 2005) although there is increasing interest by growers in the adoption of RT approaches for crop establishment for environmental and economic reasons.
Several studies have quantified GHG emissions from WOSR and/or canola systems. Management practices varied considerably between these experiments spanning 3 decades of research where the soil was cultivated by CT, RT or NT practices or the crops were established using a wide range of row spacing's and seed rates ( Table 1; Walter et al., 2015 and references therein). Another confounding factor in previous studies is that the chambers used for GHG emission measurements were either placed in the inter-row space between plant rows or enclosed both the plant row and inter-row areas. This could complicate an assessment of the impact of different management practices as the presence of plants could directly or indirectly have an impact on GHG emissions. Given the broad range of management practices used it is also difficult to generalize about the reason(s) for any observed differences in GHG emissions.
Other agronomic management practices, such as row spacing and seed rate, are often optimized by growers to encourage rapid plant emergence, canopy closure and development from sowing to harvest. Crop row spacing and plant density influence canopy architecture, which affects solar radiation interception, and the utilization of water and nutrients. Increasing row width and/or reducing the seed rate may expose surface soils to increases in solar radiation, precipitation and wind, resulting in microclimaterelated modifications in soil temperature and moisture that influence GHG emissions.
Although there are some studies that have examined the link between WOSR yields, tillage practice, row spacing or seed rate (Christian and Bacon, 1990;Vann et al., 2016) the link between cultivation practices, GHG emissions and yield has not yet been examined. The objectives of this study were to examine the effects of: 1) non-inversion tillage and row spacing, 2) row spacing and seed rate, using two WOSR cultivars, on CO 2 , N 2 O and CH 4 emissions, compared to the conventional management practice typically used on Irish farms (CT; 125 mm row spacing, 60 seeds m −2 ). Given the importance of agronomic practices that can contribute to GHG mitigation whilst having no yield penalty we also assess the effects on GWP per unit of yield (i.e., the yieldscaled GWP) in different WOSR cultivation systems.

Experiment 1: Crop Cultivation
The different cultivation treatments ("Exp. 1") were established in the Hockey field at Knockbeg (52°51′42″N, 6°56′28″W) around 5 km from the Teagasc Crops Research Facility in 2014/2015 ( Figure 1B). The soil at this site is a sandy loam-to-loam texture (Conry, 1987) and before cultivation had been under a permanent pasture for at least 10 years. Since the 1990s the site has been under CT with continuous crop cultivation (van Groenigen et al., 2011).
Winter oilseed rape (Brassica napus L. cv. Compass) was cultivated using CT, strip tillage (ST) and MT. The CT treatment consisted of a primary cultivation, soil inversion to a depth of 200-250 mm, with a 5-fin mouldboard plough followed by one pass of a roller to consolidate the soil after ploughing. Secondary cultivation (power harrowing) was performed immediately prior to sowing the seed. The ST treatment is a modified non-inversion tillage technique where strips of soil (200 mm) are cultivated to the conventional plough depth (200-250 mm) whilst the inter-row spacing between plants is left uncultivated with the previous crop residue retained on the soil surface (He-va Sub-Tiller, Denmark). The MT treatment consisted of non-inversion soil cultivation to a depth of 100-125 mm using a stubble cultivator with tines spaced 300 mm apart followed by leveling discs and a cage roller (Horsch, Terrano FX3, United States). A dynamic 3 m cultivator drill (Vaderstad, Rapid 300S, Sweden) capable of working in a range of seedbed types delivered seed through individual hydraulic metering units to a disc coulter.
Plants were sown at a rate of 60 seeds m −2 at two row spacing's, 125 and 600 mm, giving a total of five treatments as there was no 125 mm spacing for the ST treatment. The treatments (25 m × 5 m) were laid out in a randomised block design with four replications. Calcium ammonium nitrate (CAN) was applied

Soil CO 2 , N 2 O and CH 4 Emission Measurements
Circular stainless-steel collars (225 mm ∅, area 0.03 m 2 ) with rubber gaskets were inserted into the inter-row spacing (100 mm depth) of all cultivation treatments. Soil CO 2 emission measurements were carried out using an infrared gas analyzer (EGM-4, PP Systems, United Kingdom) by slowly sealing an unvented chamber (0.0034 m 3 ) onto a rubber O-ring gasket lining the collar. The enclosure time was 120 s with readings taken every 20 s. After sampling, the chamber was removed, and the CO 2 concentration was allowed to equilibrate before placement on the next collar. Gas accumulation was linear within the headspace and measurements with R 2 > 0.8 were retained for analysis (Widén and Lindroth, 2003). Additionally, in Exp. 2, four PVC Collars (160 mm ∅) were inserted into bare soil which had plants removed adjacent to the main experimental plots to estimate soil basal/heterotrophic respiration (R h ).
Emissions of N 2 O and CH 4 were measured using the manual closed chamber method (Chadwick et al., 2014). Accumulation of N 2 O within the chamber headspace, measured at four time points over 1 hour (T 0 , T 20 , T 40 , T 60 ), was determined to be linear on 84% of occasions (O'Neill et al., 2020). The unvented chambers (above) were sealed onto the collars for an enclosure period of 40 min (T 40 ) and sampling carried out between 9.00 and 13.00 h (Barton et al., 2015). The chambers were covered with aluminium foil to prevent solar radiation induced temperature changes inside the headspace during the enclosure period.
For gas sampling, headspace air was withdrawn through a stopcock fitted to the chamber vent using a 20 ml polypropylene syringe (BD Plastipak, Spain). The chamber headspace was mixed by flushing air slowly with the syringe plunger twice prior to the withdrawal of the gas sample. Using a hypodermic needle, samples were immediately transferred into pre-evacuated 7 ml glass exetainers (Sigma-Aldrich, United Kingdom) fitted with double wadded septa (Labco, High Wycombe, United Kingdom). The exetainers were injected with a 12 ml sample to create an overpressure and prevent back diffusion of ambient air during storage. Four ambient air samples were also taken near ground level before and after each sampling occasion to obtain a surrogate time zero (T 0 ) sample for each chamber (Chadwick et al., 2014;Charteris et al., 2020). Sampling frequency was increased during the period of fertilizer application: four times per week for 2 weeks, then twice a week for 2 weeks, then once a week until harvest. Additional samples were taken before or after precipitation events, which are known to stimulate denitrification (Butterbach-Bahl et al., 2013).
Analysis of N 2 O and CH 4 concentrations were carried out with a Gas Chromatograph with a 63Ni electron capture detector (ECD) at 60°C with a flame ionization detector (FID) at 300°C (Bruker Scion 456, Germany). Samples were injected into the GC using a Combi-PAL auto-sampler (CTC Analytics AG, Switzerland). Results were expressed in parts per million by volume (ppmv).
Daily GHG emissions were calculated using the Eq. 1: where ΔC/Δt is the rate of change of CO 2 (T o to T 2 ; IRGA) N 2 O and CH 4 concentration (T 0 to T 40 ; GC), where ΔC is the change in concentration of the gas in the headspace volume (ppmv or ppbv), Δt is the enclosure time period (minutes), MW is the molar mass of CO 2 -C (12 g), N 2 O-N (28 g) or CH 4 -C (12 g), P is the atmospheric pressure at the time of sampling (Pa), R is the gas constant (8.314 J mol −1 K −1 ), T is the air temperature at the time of sampling (K), V is the headspace volume within the chamber (m 3 ) and A is the area covered by the base (m 2 ). Cumulative GHG emissions (±SE) were calculated by trapezoidal integration of the daily means. The global warming potential (GWP) was calculated for each management practice by converting N 2 O and CH 4 emissions to CO 2 equivalent emissions (CO 2 -eq.) for a 100 years time horizon. The radiative forcing potential used relative to CO 2 was 265 for N 2 O and 28 for CH 4 (Myhre et al., 2013). Yield-scaled GWP was calculated by dividing the CO 2 -equivalent emissions by the seed yield (harvested at 9% moisture) and expressed in units of kg CO 2 -eq. kg −1 . The contribution of CO 2 to GWP was excluded based on: 1) the assumption that the soil CO 2 -C efflux was largely off-set by high rates of net primary productivity (C input) and biomass removal at harvest (C export)  and 2) the absence of accurate crop residue (straw and root C) data to quantify the annual change in SOC (Mosier et al., 2006).

Soil and Climatic Measurements
Soil mineral N (NH 4 + and NO 3 − ) concentrations were determined during the spring-summer growth period. Soil cores (0-10 cm depth) were taken weekly after fertilization for 1 month and every 3-4 weeks thereafter until harvest. Soils were Frontiers in Environmental Science | www.frontiersin.org September 2021 | Volume 9 | Article 716636 initially stored at 4°C and extracted either on the same day of collection, or within 24 h, using 2 M KCL at a ratio of 5:1 (v:w) water: soil (Maynard et al., 1993). The NH 4 + and NO 3 − concentrations of the extract were analyzed with an Aquakem 600 discrete analyzer (Thermo Fisher Scientific, United States). Soil volumetric moisture content (%) (GS3, Decagon Devices, United States) and soil temperature ( o C; Exp. 2 only) (ELE International, Bedfordshire, United Kingdom) were measured (50-100 mm depth) in the inter-row spacing adjacent to the collars. Water-filled pore space (WFPS) was calculated using Eq. 2: where Bd is the bulk density measured at the soil surface (70 mm depth) and Pd is particle density estimated at 2.65 g cm −3 (Linn and Doran, 1984). Climatic measurements of mean air temperature ( o C), atmospheric pressure (Pa), rainfall (mm) and soil temperature ( o C) were taken from the Met Eireann automated weather station at Oakpark, Co. Carlow, Ireland, which was located 5 km from the Exp. 1 and 30 km from the Exp. 2 sites.

Statistical Analysis
Statistical analysis was conducted using SAS 9.4 (Cary, NY, United States) and R software (R 3.6.1, R Core Team, 2019). Normality and homogeneity of variance were checked using histograms and residual graphs. Where necessary, log or square root transformations (y or y + constant) were applied to data to achieve homogeneity of variance. For Exp. 1, the PROC GLIMMIX procedure in SAS was used to test for differences between the treatments. Tillage (CT, MT, and ST) and row spacing (125 and 600 mm) were the main effects. To account for the missing treatment (i.e., ST125) in the model, additional parameters were created, each with two levels: a "Nested Group" (Factors CT + MT; Single ST) and a "Variable Group" (1 CT125, 2 all others). The _residual_ option was used in the RANDOM statement. Significant pairwise differences were determined according to the simulate posthoc test (p < 0.05). For Exp. 2, a linear-mixed effects model (lmer) in the "lme4" (Bates et al., 2012) and "lmerTest" packages (Kuznetsova et al., 2017). Variety (cvs. Compass and Troy), row spacing (125 and 750 mm) and seed rate (10 m −2 and 60 m −2 ) were set as the main effects with block and block x variety set as the random effects. Model effects were examined using the emmeans function (Lenth, 2019) where ANOVA and pairwise comparisons (normal and backtransformed data) were examined by the Tukey method at a significance level of p < 0.05. Plots were made using the "ggplot2" package in R (Wickham, 2016).

Crop Cultivation
Meteorological data was similar for both sites during the gas measurement periods (Figure 2). The mean air temperatures were 10.9°C (2.9-19.9°C) and 11.0°C (3.6-19.8°C) for Exp. 1 and Exp. 2, respectively. Cumulative rainfall was similar from the period of 27th February to 22nd July each year in Exp. 1 (261 mm) and Exp.
2 (265 mm), representing 28 and 37% of the annual rainfall, respectively. During the growing seasons, the wettest months for Exp. 1 and Exp. 2 were May (90 mm) and April (64 mm) whilst the driest months were April (26 mm) and July (24 mm), respectively.
A rapid loss of CO 2 was observed after the implementation of CT (maximum recorded value of 2.34 g C m −2 hr −1 ) at a rate 2.5 times that of ST (maximum recorded value of 0.94 g C m −2 hr −1 ) and 14.6 times that of MT (maximum recorded value of 0.16 g C m −2 hr −1 ) (Figure 3). These decreased rapidly to 0.16 g C m −2 hr −1 (CT) and 0.07 g C m −2 hr −1 (ST) after the tillage events and then remained largely stable for the rest of the measurement period. Excluding the anomalously high value found in the MT treatment after 3 days, the daily soil CO 2 emissions ranged from 0.04 to 1.38 g C m −2 h −1 across all treatments during the 12 days period.
Daily soil CO 2 emissions where generally higher during the early part of the growing season and showed two peaks, each occurring around 7 days after the first and second N fertilizer applications, respectively ( Figure 4A). Soil CO 2 emissions tended to be greater in the CT125 system with maximum daily emissions of 4.0 and 4.2 g C m −2 d −1 . The emissions of CO 2 converged towards similar values in all systems by mid-April until June, where the mean CO 2 loss was 1.16 ± 0.07 g C m −2 d −1 . Higher CO 2 emissions were generally observed between 32 and 70% WFPS, with a tendency for lower emissions under both drier (25% WFPS) and wetter (85% WFPS) conditions. Cumulative CO 2 emissions were not significantly affected by the crop cultivation treatment, with values ranging from 1,083 to 1,683 kg C ha −1 , although higher CO 2 emissions were generally found in the 125 mm row spacing treatment ( Table 3).
Slight increases in N 2 O were observed in the ST600 system (11.6-27.3 g N ha −1 d −1 ) ca. 10 days after the first fertilization ( Figure 4B). The highest N 2 O emissions were observed in the MT125 (53.8 g N ha −1 d −1 ), MT600 (62.8 g N ha −1 d −1 ) and ST600 (70.0 g N ha −1 d −1 ) treatments 4 days after the second N application and a second peak occurred in the ST600 (69.3 g N ha −1 d −1 ) treatment 2 days later. These were observed at 72-79% WFPS when the soil temperature was >10°C ( Supplementary Information, Supplementary Figure S1), however, the majority of the N 2 O emissions occurred across a wide range of WFPS (25-90%) ( Figure 4B). Maximum N 2 O emissions were lower in the CT125 and CT600 treatments, with values of 30.5 and 18.5 g N ha −1 d −1 , respectively, occurring 58 days after fertilization. There was no significant effect of tillage or row spacing on the cumulative N 2 O emissions, with values ranging from 0.81 to 2.05 kg N ha −1 ( Table 3).
Daily CH 4 emissions displayed similar temporal patterns and ranged from −10.7 to 5.4 g C ha −1 d −1 ( Figure 4C). Net CH 4 emissions preceded by significant uptake of CH 4 were observed after the second N application. With some exceptions, net CH 4 uptake was generally sustained until June. For the cumulative CH 4 uptake values, a significant tillage x nested effect (p < 0.05) was observed with a 55% increase in uptake in the CT (−0.34 ± 0.03 kg C ha −1 ) compared to the MT systems (−0.22 ± 0.03 kg C ha −1 ), with the total CH 4 uptake in the ST treatment midway between the CT and MT treatments ( Table 3).

Sowing Treatments
Daily GHG emissions for the two WOSR varieties are illustrated in Figure 5. Transient fluctuations in soil CO 2 emissions were found during the season, ranging from 0.11 to 3.20 g C m −2 d −1 ( Figure 5A), with no clear seasonal trend. Nitrogen fertilizer application tended to initially suppress CO 2 emissions, which was followed by increases in CO 2 emissions 1-2 weeks later. Soil CO 2 emissions increased temporarily in the 125/10 and 750/10 plots in June, however, little variation existed between sowing treatments during the main growing season.
Daily soil N 2 O emissions were clearly associated with N fertilizer applications ( Figure 5B). Higher N 2 O emissions were found from the cv. Troy plots particularly the 125/10 treatment, and were also observed for soil temperatures >10°C and WFPS values approaching 80% (Supplementary Information, Supplementary Figure S1). Nitrous oxide emissions fell to background rates around 1 month after the final N application.
Daily soil CH 4 emissions ranged from −38.2 to 26.3 g C ha −1 d −1 across all sowing treatments ( Figure 5C). Low net CH 4 uptake rates occurred after mineral N was applied to the soils. The sharp increase in CH 4 uptake in April coincided with a decrease in soil NH 4 + and NO 3 − concentrations in late-April ( Supplementary  Information, Supplementary Figures S2, S3). Several treatments associated with the Compass and Troy cultivars displayed transient net CH 4 emission peaks but, in general, there was a consistent trend of increasing net CH 4 uptake towards the summer months. The highest net CH 4 uptake rates occurred within a narrow temperature range of 13-15°C but a broader range (26-45%) of WFPS. Combining the data from Exp. 1 and Exp. 2, WFPS explained around one quarter of the variance (R 2 0.29) in the CH 4 emissions (Supplementary Information, Supplementary Figure S4). Neither the cumulative N 2 O nor CH 4 emissions were affected by row spacing, seed rate or variety (Table 4).

Crop Management and GHG Emissions
In this study, short-term CO 2 emissions after tillage operations were characterized by an initially rapid increase followed by a fast Frontiers in Environmental Science | www.frontiersin.org September 2021 | Volume 9 | Article 716636 exponential decline (Figure 3). The CT treatment resulted in a maximum CO 2 efflux rate ca. 2.5 times that of ST and 14.6 times that of MT which is consistent with an effect of gaseous diffusion due to soil disturbance (Jackson et al., 2003;La Scala et al., 2006;Reicosky and Archer, 2007;Morell et al., 2010). The large loss of CO 2 from the MT plot (0.34 g C m −2 hr −1 or 8.15 g C m −2 d −1 ) ca. 3.5 days after tillage coincided with a rainfall event, a phenomenon often called the "Birch Effect" (Gebremichael et al., 2019). Increased organic matter mineralization after MT in combination with favorable soil moisture and temperature may have led to the relatively larger efflux of CO 2 compared to CT and ST (Franzluebbers et al., 1995;Jabro et al., 2008). Overall, however tillage had very little impact on short-term CO 2 emissions. Soil CO 2 emissions during spring-summer were unaffected by variations in tillage intensity ( Table 3 and Table 4) in line with earlier studies that reported no significant effects of RT systems on CO 2 emissions (Kessavalou et al., 1998;Chatskikh et al., 2008;Abdalla et al., 2014). Tillage management can affect soil CO 2 emissions through its influence on soil moisture, soil temperature and soil organic C accumulation (Buyanovsky et al., 1986;Hendrix et al., 1988;Franzluebbers et al., 1995;Fortin et al., 1996). In RT systems this is often associated with the retention of crop residues that lower soil temperatures and increase soil moisture content (Chen and McKyes, 1993). However, daily soil temperatures and WFPS values varied little between treatments indicating that the effect of contrasting tillage regimes on residue retention had little effect on CO 2 emissions. In this study, the highest soil CO 2 emissions occurred during the spring period of the growing season suggesting that phenology, through its effects on root respiration processes and/or the rhizodeposition of labile C, contributed to CO 2 production (Rood et al., 1984;Whipps, 1990;Franzluebbers et al., 1995;Rochette and Flanagan, 1997;Jans et al., 2010). In OSR, the highest values for leaf area index were reached close to flowering, which is consistent with the higher soil CO 2 emissions observed in March and April across all experiments (except cv. Troy in Exp. 2). Crops allocate, on average, 21% of their photosynthetically fixed C to their roots of which a smaller proportion of this C is released to the soil in root exudates (Pausch and Kuzykayov, 2018). The absence of variation between tillage management may reflect either a low supply of labile C belowground or labile C substrate was not limiting CO 2 production. Tillage management had little impact on N 2 O emissions in this study ( Table 2 and Table 3) in line with earlier work examining CT and RT systems (Abdalla et al., 2010;Liu et al., 2016). Peak N 2 O emissions that were found in the MT and ST but not in the CT treatments coincided with a higher WFPS (>70%) and higher temperatures (>10°C) as noted in previous studies (Adviento-Borbe et al., 2007;Abdalla et al., 2010;Žurovec et al., 2017;O'Neill et al., 2020). This indicates that denitrification is the main source of N 2 O production and the N 2 O emissions are largely dependent on the extent of soil anaerobiosis. A lower WFPS and an associated increase in soil oxygenation may explain the absence of any impact of fertilization on N 2 O emissions in the CT systems ( Figure 4B).
Cumulative CH 4 uptake was ca. 55% greater in the CT compared to the MT treatments (p < 0.05, Table 2). Plaza-Bonilla et al. (2014) also reported a significantly higher cumulative uptake of CH 4 in CT (2.69 kg C ha −1 ) compared to NT (1.16 kg C ha −1 ) under Mediterranean dryland conditions (p < 0.05). The authors suggested that the greater CH 4 oxidation found under CT might be explained by the short duration of NT (3 years) and the possible lack of differences in soil pore structure and methanotrophic communities between CT/NT plots. Reduced tillage practices can however result in a more porous and stable soil structure that facilitates CH 4 diffusion (Ball et al., 1997;Hütsch, 1998;Ussiri et al., 2009) and the long-term implementation of RT (>40 years) has been shown to restore the CH 4 oxidation capacity of arable soil (Ussiri et al., 2009;Jacinthe et al., 2014). Although the reason(s) for the observed differences between tillage treatments in this study are not clear, it could be explained by both the enhanced diffusion of atmospheric CH 4 into the soil and the greater CH 4 oxidation rates under CT as the soils became drier in summer.
Soil WFPS had the largest effect on CH 4 emissions and explained approximately one quarter of the variability in the daily CH 4 emissions (Supplementary Information, Supplementary Figure S4). Higher soil water content inhibits the diffusive transport of CH 4 and oxygen to active microbial sites, consequently reducing CH 4 uptake in arable soils (Flessa et al., 1995;Drewer et al., 2012). The results of this study are consistent with the recent review by Cowan et al. (2020) who found soil volumetric water content as the strongest predictor of CH 4 emissions across agricultural soils in the United Kingdom and Ireland (R 2 < 0.1). This indicates that soil water content and its impact on diffusion processes may override the effect of tillage management in regulating CH 4 exchange in arable soils. Row spacing was found to have no significant effect on GHG emissions. In OSR, row spacing influences plant height, branching, leaf area, pod number and the overall canopy architecture. Plants sown in narrow rows reach canopy closure quicker than wide row cropping systems and thus, it was hypothesised that row spacing would influence factors such as root distribution, root C input, water use, N uptake, light interception and soil temperature, and directly or indirectly influence GHG emissions (Sharratt and McWilliams, 2005). Zapata et al. (2021) reported higher CO 2 emissions from soybean compared to wheat (wide vs. narrow row spacing) due to prolonged exposure of surface soil to direct sunlight and high intensity rainfall. The absence of a row spacing effect in this study may be due to comparable above-ground vegetative growth rates and rapid canopy closure, irrespective of plant stand structure, thereby limiting the potential microclimatic effects of row spacing on GHG emissions.
Seed rate was a significant factor affecting GHG emissions in Exp. 2 ( Table 3). Mean cumulative CO 2 emissions (Compass and Troy) were 21% higher in the 10 seeds m −2 compared to the 60 seeds m −2 treatments (p < 0.05). Seed rate determines the plant density and the canopy architecture and, like row spacing, will influence light interception, water and nutrient use and soil  (Pribyl, 2010). Frontiers in Environmental Science | www.frontiersin.org September 2021 | Volume 9 | Article 716636 microclimate. For OSR crops, Roques and Berry (2016) found that high seed rate plots accelerated flowering and earlier maturation whereas plots with low seed rates remain greener for longer. The longer a plant retains its leaves, or delays leaf senescence, the greater the photosynthetic capacity of the stand and the amount of labile C that is released via root exudates to the rhizosphere. Soil CO 2 emissions in the bare soil plots (R h ) also tended to be higher during pod development in June (data not shown). This result points to a parallel increase in the rate of OM decomposition in the cropped soils caused by a reduction in shading and increases in soil temperature in the 10 seeds m −2 treatments. An increase in belowground OM inputs during vegetative growth and/or decomposition of native SOM at maturity may have both contributed to the greater soil CO 2 emissions.

Variation in Yield and GWP Across Different Management Practices
Despite the use of a wide range of a management practices (tillage intensity, row spacing, seed rate and variety), seed yield was unaffected ( Table 2 and Table 3). These results support earlier studies that found no effect of tillage, row spacing and seed rate on final seed yields (Degenhardt and Kondra, 1981;Christensen and Drabble, 1984;Christian and Bacon, 1990;Bonari et al., 1995;Vann et al., 2016;Wynne et al., 2020). Although seed yield and yield components are determined by plant density (Leach et al., 1999;Diepenbrock, 2000;Rathke et al., 2006;Kuai et al., 2015), individual OSR plants can counter variations in plant density by increasing yield per plant (Diepenbrock, 2000), so that the seed yield is not compromised by sowing technique.
The two varieties also produced similar seed yields ( Table 3). The shorter height of the cv. Troy plants had no detrimental impact on yield components and many of the commercially available semi-dwarf genotypes produce similar yields to conventional varieties (Sieling and Kage, 2008). Miersch et al. (2016) noted that semi-dwarf varieties yielded higher than the conventional types when N is the limiting resource. The nonlimiting soil mineral N concentrations in this study (0-10 cm)  may have resulted in the similar yields, although semi-dwarf varieties, such as cv Troy, may also perform better in low N input systems compared to conventional varieties. Semi-dwarf varieties may also be favoured by farmers for their higher lodging resistance plus easier harvesting and pesticide applications (Sieling and Kage, 2008). Yield-scaled GWP values were unaffected by crop management practice. In Exp. 2, cv. Troy had a yield scaled GWP that was double that of cv. Compass (p < 0.05). Cumulative N 2 O emissions did not significantly differ between treatments, however, the greater GWP was attributed to N-fertilizer applied in excess in the 10 seeds m −2 cv. Troy plots, as indicated by the large N 2 O emission peaks after the third N application ( Figure 5). Sowing at a higher seed rate may circumvent the high GWP of cv. Troy but puts the variety at an inferior position to cv. Compass insofar as low seed rates can achieve similar yields to conventional sowing without adversely affecting GWP.
Our yield-scaled GWP values 0.05-0.16 kg CO 2 -eq. kg −1 are within the range of those for crops reported in the literature. Linquist et al. (2012) calculated a yield-scaled GWP of 0.166 for wheat and 0.185 kg CO 2 -eq. kg −1 for maize in their metaanalysis of 57 cropland sites which included rice 0.657 kg CO 2 -eq. kg −1 . Rajaniemi et al. (2011) found slightly higher yield-scaled GWP values for barley, oats, wheat and rye (0.57, 0.57, 0.59 and 0.87 kg CO 2 -eq. kg −1 ). Gan et al. (2012) calculated a farm-scale carbon footprint of 0.281-0.317 kg CO 2eq. kg −1 for barley succeeding oilseed rape. These studies, however, included the GWP of additional upstream sources such as fuel emissions, production and use of fertilizers and seed production.
The short measurement campaign (6-months) may pose some limitations for the longer-term interpretation of results. We did not take gas measurements after tillage in Exp. 2 meaning we may have missed any transient emissions of CO 2 , N 2 O or CH 4 , whilst the low temporal frequency of measurements in the winter/spring of Exp. 1 could have introduced large errors in the estimation of the cumulative GHG emissions when determined by linear interpolation. To reduce any bias, only the higher temporal frequency measurements made during spring were considered as representative of yield scaled GWP. Annual GWP values are also required to cover the fallow and post-tillage periods as these can be a large source of GHG emissions from croplands (Linquist et al., 2012). Including all factors leading to CO 2 emissions onfarm (tillage, herbicide, pesticide, fertilization) as well as upstream activities (harvest, fertilizer manufacturing, other chemical inputs, transport) (Sainju, 2016) would provide a more complete assessment of the GHG balance of WOSR production systems.
FIGURE 5 | Soil CO 2 , N 2 O and CH 4 emissions for each sowing method treatment in Exp. 2: the conventional ("cv. Compass") and semi-dwarf ("cv. Troy") biomass varieties sown at 125 mm or 750 mm row spacing with 10 or 60 seeds m −2 . Vertical dashed lines indicate fertilizer application and the vertical lines on each data point represent standard error of the means.
Frontiers in Environmental Science | www.frontiersin.org September 2021 | Volume 9 | Article 716636 CONCLUSION Overall, management practice had a minimal effect on GHG emissions, crop yield and yield scaled GWP. Soil tillage and row spacing had no significant impact on CO 2 and N 2 O emissions. Consequently, the management approach used for WOSR cropping systems can be tailored to the particular conditions and economic factors. The low seeding rate resulted in higher CO 2 emissions irrespective of WOSR variety, whilst the yieldscaled GWP of cv. Troy was, on average, higher than that of the cv. Compass. Sowing the semi-dwarf variety at low seeding rates may therefore increase the risk of higher C and/or N emissions from OSR systems. The overall findings indicate that considerable flexibility is possible in the way that OSR crops are established, grown and managed according to local/regional conditions and economics, without compromising yield or GHG emissions. Longer-term, high temporal resolution measurements are required to examine whether the observed GHG and GWP values can be attributed more to management or climatic effects.

DATA AVAILABILITY STATEMENT
The datasets presented in this article are not readily available because The data analyzed in this study is subject to the following licenses/ restrictions: confidential. Requests to access the datasets should be directed to Macdara O'Neill, macdara.oneill@teagasc.ie.

AUTHOR CONTRIBUTIONS
MO'N, GL, PF, and BO conceived of and designed the experiments. MO performed the experimental work and gaseous emission measurements. MO'N, GL, PF, and BO wrote the manuscript.

FUNDING
This work was supported through the Teagasc Walsh Fellowship Scheme awarded to MO'N (RMIS 6363).