Ammonia Emissions, Exposed Surface Area, and Crop and Weed Responses Resulting from Three Post-Emergence Slurry Application Strategies in Cereals

: Ammonia (NH 3 ) emissions resulting from the ﬁeld application of livestock slurry has both negative human health and environmental impacts. However, decreasing the exposed surface area (ESA) of slurry upon application can reduce NH 3 volatilization by limiting its atmospheric exposure. In the present study, three strategies for depositing slurry within a growing crop were studied, including: 1. standard trailing hoses (SA hose ), 2. trailing shoes (SA shoes ), and 3. the combination of rigid tines and trailing shoes (SA tines+shoes ). Application methods interact with the soil to varying degrees and were evaluated within the context of contemporary weed management practices, namely in cereals receiving inter-row hoeing. SA hose , SA shoes , and SA tines+shoes were compared in three coinciding experiments that assessed slurry ESA, NH 3 emissions, and crop and weed effects. SA tines+shoes resulted in smallest ESA, 70–72% and 61–66% less than SA hose and SA shoes , respectively. However, in only one of three site–years did SA shoes and SA tines+shoes reduce NH 3 emissions compared to SA hose , by 46% and 29%, respectively. Crop yields, nitrogen (N) accumulation in crop biomass, and intra-row weed biomass were unaffected by the placement method. In heavily crusted soils, the SA tines+shoes prototype worked well; however, the functional differences among placement strategies were not great enough to detect crop and NH 3 effects.


Introduction
Ammonia (NH 3 ) emissions resulting from the land application of animal manure constitutes a significant agricultural pollutant.The direct and indirect effects of NH 3 pollution are detrimental to both human and environmental health.Once airborne, NH 3 contributes to the formation of fine particulate matter in the atmosphere [1] and the development of disease in humans upon inhalation [2].Wet and dry deposition of atmospheric nitrogen (N) also adds to the acidification and eutrophication of terrestrial and aquatic ecosystems [3,4].Within the EU-28, approximately 93% of NH 3 emissions can be traced to agricultural activities [5], with 26% of total emissions stemming from animal manure application to soils in 2018 [6].The EU has set substantial goals for reducing NH 3 emissions by 19% within the 25-year period from 2005 to 2030 [7].Achieving these objectives requires a multifaceted approach focused on improved animal rearing techniques, on-farm manure handling, storage, and processing, as well as the development of field application strategies that reduce NH 3 volatilization.
Methods known to effectively mitigate ammonia emissions during land application of liquid manure (slurry) include strategies that reduce atmospheric exposure by increasing Agronomy 2022, 12, 2441 2 of 19 infiltration rate and reducing exposed surface area (ESA) [8,9].As slurry infiltrates into the soil, the quantity of aboveground N susceptible to loss decreases, thus reducing emission potential.NH 3 volatilization occurs at the interface between exposed slurry and the atmosphere; therefore, when aboveground ESA is reduced or infiltration rate increases, emissions decrease.
Existing methods for low-emission slurry application include surface application with trailing hoses, shallow application with trailing shoes, and sub-surface application with open and closed slot injection [10].Open slot injection places slurry in a channel that is not covered by soil, while closed slot injection places slurry below the soil's surface, so it is not exposed to the atmosphere.Closed slot injection is considered the most effective strategy for reducing NH 3 emissions [11,12].When implemented before planting a cereal crop, injection improves yields by reducing atmospheric N loss and by placing fertilizer near the crop's rooting zone [13].An associated increase in crop vigor may also aid in the suppression of weeds, further contributing to increased yields [14,15].However, injection also has its drawbacks.Weed response may differ among species with contrasting root morphologies; Melander et al. [16] report that upon fertilizer injection, the biomass of shallow-rooted species are reduced, while deeply rooted species benefit.Closed slot injection is also highly disruptive, rendering its implementation impractical within a growing crop, where significant damages are not offset by improved N conservation and crop acquisition [17].Therefore, closed slot injection is best suited for bare field settings preceding sowing, and open slot injection is most applicable in grasslands.Ample horsepower is also required to accommodate high draught forces associated with pulling injectors; this limits implement width and further contributes to crop damages and soil compaction resulting from additional tractor passes [18].Finally, a tradeoff in volatile organic compound emissions must be considered.Compared to surface application, closed slot injection reduces NH 3 emissions by approximately 61% but increases emissions of the potent greenhouse gas, nitrous oxide (N 2 O), by 196%; this "pollution swapping effect" is described in greater detail by Emmerling et al. [19].
Methods less disruptive than injection also exist, suitable for post-emergence application within a growing cereal crop.Camera and GPS guidance technologies have made the accurate placement of liquid manure between narrowly spaced rows possible, using tools that either deposit at the surface or interact shallowly with the soil.These strategies also demand less power during operation and can therefore be used at increased working rates, covering a wider surface area with each tractor pass [18].
In Northern Europe, trailing hoses represent the standard method for placing liquid manure in a growing cereal crop.Hoses drag across the soil surface, depositing slurry in bands between crop rows, reducing ESA and NH 3 volatilization compared to broadcast, full-width spreading [20].Placing slurry beneath the canopy of a growing crop has also been shown to reduce NH 3 emissions [9,21,22] and is, therefore, a method well suited to springtime application in established winter cereals.However, trailing hoses are prone to side-to-side shifting during operation and may bounce on top of the soil, causing some splashing; thus, hoses are not ideal when targeting fertilizer placement in close proximity to crop rows.In addition, when compared to trailing shoes, hoses may produce higher ESA and NH 3 emissions [23].
A variety of trailing shoe designs are now available on the market.Typically, these shoes function by causing some shallow soil disturbance, creating a lightly cultivated zone or depression ahead of where the slurry is placed.Cultivation immediately preceding slurry deposition reduces NH 3 loss by increasing infiltration rate [24].In addition, some shoes are designed to funnel slurry into a narrow channel, reducing ESA and NH 3 emissions [25,26].Trailing shoes are more rigid in their construction compared to hoses, preferable for targeting placement near crop rows.However, if substantial soil crust is present, trailing shoes may not function as intended, riding on top of the soil without causing any disturbance; the present study investigates whether combining the operation of rigid tines and trailing shoes will help resolve this issue.When operating any slurry placement tool within a growing crop, the risk of inflicting damages that may negatively affect yields exist simultaneously alongside the potential positive effects, including N conservation, improved infiltration, and improved crop uptake.Therefore, when comparing slurry application strategies, it is necessary not only to evaluate resulting ESA and NH 3 emissions; one must also assess crop and weed responses.
In the present study, three slurry application methods are compared.First, standard surface application with trailing hoses.Second, application with trailing shoes designed to create a shallow furrow ahead of where the slurry is placed.Third, a novel prototype applicator composed of rigid tines, followed by trailing shoes.Tines were added to break through any hardened crust that has formed on the soil's surface; field conditions typical for loam-and clay-type soils in spring, when slurry is spread in over-wintered cereals.Application methods are evaluated within the inter-row hoed cereal system, an effective weed management strategy gaining popularity among organic growers in Northern Europe [16,27].Hoeing is a robust and highly effective method for controlling inter-row weeds [28]; however, hoeing has little effect on weeds growing in the intra-row zone, capable of negatively affecting crop yields [29].
This study aims to evaluate how the three slurry application methods described affect crop and intra-row weed response, NH 3 emissions, and slurry ESA.It is hypothesized that similar crop yields, biomasses, and grain quality responses will be observed across slurry application treatments.It is anticipated that damages resulting from more disruptive placement strategies (trailing shoes and tines plus trailing shoes) will be offset by increased N uptake due to improved infiltration and targeted placement.While total intrarow weed biomass may not differ among application strategies, individual weed species are expected to respond differently depending on key morphological traits.Taprooted weeds are hypothesized to experience improved growth and inferior control with targeted placement tactics (trailing shoes and tines plus trailing shoes), while weeds with shallow fibrous roots will experience reduced growth.Finally, greater slurry NH 3 emissions and ESA are expected to result from the standard trailing hoses, and these measures will be reduced with the trailing shoes and further reduced with the tines plus trailing shoes.The present study contributes a comprehensive comparison of three low-draft slurry placement methods well suited for post-emergence use in cereals; one prototype and two pre-existing designs are evaluated.Coinciding experiments allow results to be contextualized within the interdependent dynamics of slurry ESA, NH 3 emissions, plant N uptake, and crop and weed effects.

Materials and Methods
Separate experiments were required to evaluate the effect of liquid manure application strategy on crop and intra-row weed response (EXP CROP+WEED ), NH 3 emissions (EXP NH3 ), and slurry ESA (EXP ESA ).In 2019, one replication of EXP CROP+WEED (EXP CROP+WEED.A ) and EXP NH3 (EXP NH3.A ) was performed simultaneously at the Foulum Research Center, Denmark (56.49• N, 9.58 • E) in a field possessing a loamy-sand-soil texture and using spring barley (Hordeum vulgare L., cv.RGT Planet) as the test crop (SITE A ). From 2019 to 2020, experiments were performed at the Flakkebjerg Research Center, Denmark, using winter wheat (Triticum aestivum L., cv.Cheriff) as the test crop.In fields possessing a sandyloam-soil texture, one replication of EXP CROP+WEED (EXP CROP+WEED.B ), two replications of EXP NH3 (EXP NH3.B1 and EXP NH3.B2 ), and one EXP ESA trial (EXP ESA.B ) were performed (SITE B ; 55.32 • N, 11.39 • E).A second repetition of EXP CROP+WEED (EXP CROP+WEED.C ) was implemented at a nearby site (SITE C ; 55.33 • N, 11.41 • E).Spring barley and winter wheat were sown uniformly across experiments and sites.At SITE A , spring barley was sown at 24 cm row spacing and a target density of 300 plants m −2 .At SITE B and SITE C , winter wheat was sown at a 25 cm row spacing and a target density of 300 plants m −2 .Table 1 provides an overview of experiments performed across site-years.For all EXP CROP+WEED trials, three factors were evaluated within a full factorial, randomized complete block design with four replications.Experimental factors included pig slurry application method, N rate, and weed management treatment.In total, three post-emergence slurry application methods were studied, including (i) standard trailing hoses (SA hose ), (ii) trailing shoes (SA shoes ), and (iii) a prototype model comprising a rigid tine, designed to break through the soil's crust, followed by a trailing shoe (SA tines+shoes ).Two N rates, calculated as kg NH 4 -N ha −1 of pig slurry, were evaluated, (i) 100 kg ha −1 (N 100 ) and (ii) 50 kg ha −1 (N 50 ).Finally, two weed management treatments were evaluated, (i) weedy plots, receiving no herbicide treatment (WMT weedy ), and (ii) weed-free plots, receiving herbicide treatment (WMT weed-free ).All plots received inter-row hoeing with 16-cm shares.Plot dimensions were 3 by 13 m at SITE A and 2.5 by 12.5 m at SITE B and SITE C .A total of 96 plots were included in each experiment; 3 slurry application methods * 2 N rates * 2 weed management treatments * 4 blocks.Table S1 provides a record of the implementation date for each EXP CROP+WEED field operation, as well as crop, surrogate weed, and ambient weed growth stages.
Two surrogate weeds with differing root morphologies were sown together in a strip, horizontally through each plot's center.Surrogate weeds included white mustard (Sinapis alba L., cv.Lotus), possessing a deep taproot, and Italian ryegrass (Lolium multiflorum Lam., cv.Fabio), exhibiting shallow fibrous root growth.The target population of both surrogate weed species S. alba and L. multiflorum was 30 plants m −2 ; therefore, the combined target density of both surrogate weed species was 60 plants m −2 .At SITE A , surrogate weeds were sown in strips 3 m wide within four days of planting the crop and established successfully (Table S1).At SITE B and SITE C , surrogate weed sowing was attempted in late March 2020.The planter was raised during operation to avoid damaging the growing winter wheat crop, dropping seeds onto the soil's surface; however, both S. alba and L. multiflorum failed to establish due to a lack of precipitation following planting.In place of surrogate weeds, two prevalent and naturally occurring ambient weed species were studied instead.Ambient weeds studied at SITE B and SITE C included mayweed (Matricaria L.), possessing a deep taproot, and annual bluegrass (Poa annua L.), exhibiting shallow fibrous root growth.
Slurry was applied within the growing cereal crop in the inter-row approximately 8 cm from crop row centers.Different trailing shoe models were used when implementing the SA shoes and SA tines+shoes treatments at SITE A versus SITE B and SITE C .In addition, differing prototypes combining the operation of a rigid tine and trailing shoe were used in SA tines+shoes plots at SITE A versus SITE B and SITE C .Bomech trailing shoes (Bomech B.V., Albergen, Netherlands) were employed at SITE A .The arm connecting the Bomech shoe to the toolbar is fabricated from steel, as is the V-shaped hull that runs along the bottom of the shoe responsible for creating a shallow soil opening where the slurry is deposited.The Bomech shoe can respond to greater downward pressure owing to its rigidity; therefore, when operating in lighter soils at SITE A (Table 1), the tool functioned as expected, creating a furrow in the soil while placing slurry in SA shoes plots.The prototype used in SA tines+shoes plots at SITE A was designed, so the tine and the Bomech shoe operated together, mounted as a single unit on the same toolbar (Figure S1).It was observed that the tine flexed and bounced abruptly during operation while breaking through deeper soil layers; thus, the Bomech shoe moved with the tine, resulting in some splashing of slurry.Vogelsang trailing shoes (Vogelsang GmbH & Co. KG, Essen/Oldb., Germany) were used at SITE B and SITE C (Video S1).The arm connecting the Vogelsang shoe to the toolbar is composed of rigid plastic and, therefore, more flexible than the Bomech model.When operating at SITE B and SITE C in heavy soils with significant surficial crusting, the Vogelsang shoes could not penetrate the soil crust on their own.In SA shoes plots, the Vogelsang shoes rode evenly on top of the soil, depositing slurry without creating a depression.The prototype used in SA tines+shoes plots at SITE A was altered for use at SITE B and SITE C .The rigid tine and the Vogelsang trailing shoe were mounted on separate toolbars, thus operating independently (Video S2).The tine was able to break through the hardened soil crust, and the Vogelsang shoe rode smoothly behind in the narrow opening created by the tine while depositing slurry.For reference, Video S3 shows slurry application with trailing hoses (SA hose ); this method was used at SITE A , SITE B , and SITE C .
Slurry analyses were performed using standard methods for dry matter content [30], total nitrogen [31], and total ammoniacal N (TAN) [32].Viscosity measurements were performed using a rotational viscometer (DV-II + P Viscometer, Brookfield, AMTAK, Middleboro, MA, USA) with an LV-2 spindle at 10 RPM and slurry temperature of 19 ± 0.8 • C. All analysis results and application rates for EXP NH3 and EXP ESA can be found in Table 2. Achieved slurry application rates in EXP CROP+WEED.A , N 50, and N 100 treatments were 2.3 kg m −2 (3.54 g NH 4 -N m −2 ) and 4.6 kg m −2 (7.08 g NH 4 -N m −2 ), respectively; and in EXP CROP+WEED.B and EXP CROP+WEED.C , N 50 and N 100 treatments were 1.5 kg m −2 (5.29 g NH 4 -N m −2 ) and 3.0 kg m −2 (10.57g NH 4 -N m −2 ), respectively.Following application, the slurry surface pH was measured with a flat-tip pH electrode (Orion TM 8135BN ROSS TM , Combination Flat Surface pH Electrode, Fischer Scientific, Loughborough, UK).Measurements were taken as close to the slurry-air interface as possible.Results can be seen in Figure S2 in Supplementary Materials.

Data Collection
To avoid edge effects, samples were not collected from the outermost rows or from the first 50 cm at the top and bottom of each plot.Achieved in-field crop, surrogate weed, and ambient weed densities (plants m −2 ) were measured before herbicide application, inter-row hoeing, or slurry application.Crop density (plants m −2 ) was calculated from four counts of 1-m length row per plot.Intra-row surrogate and ambient weed counts were made within 8 by 100 cm quadrats (0.08 m 2 ), centered on random sections of crop row.Intra-row ambient weed density (plants m −2 ) was recorded individually for each species observed within four quadrats per plot.At SITE A , counts of intra-row surrogate weed density (plants m −2 ) were made within two quadrats per plot, positioned in the surrogate weed strip.
After barley surpassed principal growth stage six (flowering and anthesis), intrarow crop, surrogate weed, and ambient weed biomass samples were cut from six 8 by 100 cm quadrats per plot.Quadrats were centered on crop rows, encompassing only the intra-row zone.Having been treated with herbicide, surrogate weed strips were absent in WMT weed-free plots; thus, the distribution of the six biomass samples taken within each plot was slightly different for WMT weed-free and WMT weedy .In WMT weed-free plots, all six biomass cuts were made randomly throughout the plot and combined into a single sample.In WMT weedy plots, four biomass cuts were made outside the surrogate weed strip, and two cuts were made inside the surrogate weed strip; biomass samples taken from these zone were evaluated separately.All samples were sorted into four categories.At SITE A, plant biomass was divided into (i) spring barley, (ii) S. alba, (iii) L. multiflorum, and (iv) ambient weeds.At SITE B and SITE C , biomass was divided into (i) winter wheat, (ii) Matricaria, (iii) P. annua, and (iv) other ambient weeds.All biomass samples were dried for at least 24 h (h) at a temperature of 80 • C and then weighed.
Following the collection of plant biomass dry weights, biomass samples from WMT weedy plots were pulverized.Sub-samples of ground plant tissue were analyzed to determine percent N using the Dumas combustion procedure [33,34].Information on N concentration was employed to calculate N accumulation (g N m −2 ) among the four previously described plant categories.
Before harvesting the crop, the surrogate weed strip in the plot's center was cut away.Final plot length and the number of harvested rows were recorded, and crop yields (kg ha −1 ) were calculated based on the final plot area.Harvested grain was cleaned, and final yield weights were standardized to 15% moisture content.
Grain quality measures were obtained after harvest.Moisture content (%), protein content (%), and bulk density (kg hL −1 ) were measured using a near-infrared spectroscopy analyzer (Infratec TM 1241 Grain Analyzer, Foss A/S; Buchmann et al., 2001).Thousand kernel weight (g) was acquired by weighing four samples of 200 grain kernels per plot.

Data Analysis
JMP ® software, version 14.3.0(SAS Institute Inc., Cary, NC, USA) was used to analyze EXP CROP+WEED data.Each site (SITE A , SITE B , and SITE C ) was evaluated separately.All crop effects were analyzed using a three-way analysis of variance (ANOVA); subsequent means comparisons were performed using Tukey's HSD tests.In the three-way ANOVA, fixed variables consisted of main effects, weed management treatment (WMT), nitrogen rate (N), slurry application method (SA) and interaction effects, WMT*N, WMT*SA, N*SA, and WMT*N*SA; block served as a random variable.
Data from WMT weedy plots were analyzed to evaluate intra-row weed effects using a two-way ANCOVA.Subsequent means comparisons were performed using Tukey's HSD tests.In the two-way ANCOVA, N, SA, and N*SA served as fixed variables, and block served as a random variable.The corresponding measure of intra-row weed density (plants m −2 ) served as the continuous explanatory variable when analyzing each class of weed biomass data.For example, when analyzing P. annua biomass (g m −2 ) data, P. annua density (plants m −2 ) serves as the continuous explanatory variable.It is important to note that intra-row weed density measures were recorded before implementing inter-row hoeing, herbicide application, or slurry application; therefore, variables are not confounded and meet the assumption of independence.
Anderson-Darling tests were applied to residuals to ensure a normal distribution, and Levene's tests were performed to ensure equality of variance (α = 0.05).Where required, log 10 , log 10 (x + 1), square root, and square root (x + 1) transformations were performed to meet model assumptions.

Ammonia Emissions (EXP NH3 ) 2.2.1. Experimental Design, Treatment Specifications, and Data Collection
In plots adjacent to the EXP CROP+WEED trials, relative NH 3 emissions were compared among the three slurry application strategies in a randomized complete block design with three replications.The experiment was repeated a total of three times, once at SITE A (EXP NH3.A ) and twice, in separate locations, at SITE B (EXP NH3.B1 and EXP NH3.B2 ).EXP NH3.A was initiated on 15 May 2019, the same day that slurry was applied in EXP CROP+WEED.A .Experiments at SITE B were started on 17 April 2020, the day that slurry was applied in EXP CROP+WEED.B (EXP NH3.B1 ), and one week later on 24 April 2020 (EXP NH3.B2 ).A target application rate of 100 kg NH 4 -N ha −1 was used; however, a computational error resulted in a low application rate of 35.4 kg NH 4 -N ha −1 at EXP NH3.A (Table S2).
Ammonia emissions were measured using a system of nine dynamic chambers (wind tunnels) and continuous online measurements with a cavity ring down spectroscopy (CRDS) instrument (G2103 NH 3 Concentration Analyzer, Picarro, CA, USA).See Pedersen et al. [23] for a detailed description of the wind tunnel system and instrumentation used to measure NH 3 emissions in the present study.To assess SA hose , SA shoes , and SA tines+shoes treatments, equipment described in Section 2.1.1 was operated within plots without applying slurry, creating the soil disturbance associated with each application strategy; liquid slurry was then applied manually using a watering can, and measurements commenced immediately thereafter.Within each tunnel, NH 3 concentration in the air was logged every 1 to 2 s.In addition, background concentrations of NH 3 were measured at the air intake of each tunnel.The total retention time for each tunnel was 104 min and the total measuring time for each tunnel was 113 h.Before beginning experiments, recovery throughout the system was tested with a standard NH 3 gas and was found to be 96% at minimum.

Data Analysis
An average of the last 30 s of measurements for each measuring cycle was used for emission calculations.An average of the three backgrounds was subtracted for each measurement cycle.Concentrations lower than the detection limit (three times the standard deviation of the background measurements) were set to zero.The NH 3 flux (F NH3 [g min −1 m −2 ]) was calculated from the NH 3 concentration (C), the volumetric airflow rate (q) and the emitting area (A) (Equation ( 1)).
The cumulative emissions were calculated using the trapezoid rule [35].Differences among treatments were analyzed for each experiment individually using a one-way ANOVA with treatment as the independent variable.In each experiment, the nine wind tunnels were divided into three blocks, each containing one tunnel with each treatment, with a randomized block design.A single wind tunnel was used as an observational unit, and the cumulative NH 3 was used as the response variable.Tukey's HSD test (confidence interval of 95%) was used to investigate differences among treatments.Homogeneity of soil and crops within the experimental area was assumed.experiment, the achieved slurry ESA among SA hose , SA shoes , and SA tines+shoes treatments were studied.The experiment, EXP ESA.B , possessed randomized complete block designs with three replications.Exposed surface area of the slurry following application was measured by employing methods developed and presented by Pedersen et al. [36]; fluorescent dye was used to detect and quantify slurry on top of the soil, and images collected within dark chambers were analyzed to determine ESA (m 2 m −2 ).Exposed surface area measurements were taken every 10 min for a period of 60 min, followed by every 20 min for a period of 80 min, resulting in a total of 10 images per dark chamber.To mimic SA hose , SA shoes , and SA tines+shoes treatments, the same methods described in Section 2.2.1 were used; soil disturbance was replicated by operating equipment within treatment plots, and then slurry was applied by hand using a watering can.

Data Analysis
The time decay of slurry ESA was analyzed using a generalized linear mixed model defined with the Gamma distribution and the logarithmic link function.Responses were observed as change in slurry ESA per square meter (m 2 m −2 ) over time for each dark chamber; the model contained two fixed effects: one factor representing the dark chamber and a continuous explanatory variable (regression) representing the logarithm of the time shifted.
According to the model, ESA at time t (with t = 0, 0.015, 0.334, 0.484, 0.667, 0.834, 0.984, 1.35, 1.65, 1.984, 2.317, and 2.6 h) in the b th block, is represented by the random variable A bt , and is conditional on the residual random component U bt , which is Gamma distributed with conditional expectation given by: which is equivalent to: Model specification is completed by stating that the random components U 10 , . . .,U 2.66 are identically distributed with U 10 ∼ N 0, σ 2 r , and for b = 1, 2, and 3, Moreover, the residual random component corresponding to each dark chamber is independent per construction.In this way, we account for dependencies between the observations made at the same dark chamber at different observation times.
The time shift parameter γ was estimated by maximizing the likelihood profile, obtained with marginal models, and defined with fixed values of γ.Additionally, initial ESA (obtained at time t = 0) and final ESA (obtained at the last observation time, t = 2.66) were estimated as exp(α)γ β and exp(α)(γ + 2.66) β , respectively.We used parametric bootstrap techniques, with 1000 bootstrap replicates, to construct confidence intervals and treatment comparison tests, for the parameter and the quantities calculated with the parameters referred above.

Crop and Weed Responses
Achieved crop densities at SITE A , SITE B , and SITE C were 240, 242, and 363 plants m −2 , respectively.At SITE A , surrogate weeds S. alba and L. multiflorum were established at densities of 20 and 28 plants m −2 , respectively.At SITE B and SITE C , surrogate weeds failed to establish due to lack of moisture; therefore, ambient weeds, P. annua and Matricaria were studied, each occurring naturally at an average density of 8 and 11 plants m −2 at SITE A , respectively, and 64 and 17 plants m −2 at SITE B .Ambient weed densities at SITE A SITE B and SITE C were 207, 103, and 33 plants m −2 , respectively.Ambient weed populations at each site, characterized by the three most predominant species were, at SITE A , 46% Common lambsquarters (Chenopodium album), 23% European field pansy (Viola arvensis), and 12% P. annua, at SITE B , 33% Matricaria, 23% P. annua, and 19% Stickywilly (Galium aparine), and at SITE C , 63% P. annua, 17% Matricaria, and 8% Quackgrass (Elymus remens).
Winter wheat yields (kg ha −1 ) were on average 11% and 7% greater in WMT weed-free plots compared to WMT weedy , and 35% and 54% greater in N 100 plots compared to N 50 at SITE B and SITE C , respectively (Table 3).Across treatments at SITE C , SA hose resulted in winter wheat yields 11% greater than SA tines+shoes .However, differences were not significantly different (p > 0.05) at the other sites.Crop biomass was consistently affected by nitrogen rate across site-years (Table S2).On average, spring barley biomass (g m −2 ) was 12% greater at N 100 compared to N 50 at SITE A , and winter wheat biomass (g m −2 ) was 27% and 38% greater at N 100 compared to N 50 at SITE B and SITE C , respectively.At SITE B , winter wheat biomass was also 7% greater in WMT weed-free plots than WMT weedy plots, and SA tines+shoes produced 13% less biomass than SA hose and SA shoes .
Weed biomass (g m −2 ) was largely unaffected by the treatment variables studied.Nitrogen rate and slurry application strategy did not affect the biomass of taprooted surrogate weeds, S. alba at SITE A , or Matricaria at SITE B and SITE C (Table 4).The biomass of P. annua, possessing shallow fibrous roots, was 66% greater at N 100 compared to N 50 across slurry application treatments at SITE B .However, no significant effects were observed for P. annua at SITEc or L. multiflorum at SITE A (Table 5).The intra-row biomass of naturally occurring ambient weeds was also unaffected by treatments across sites (Table 6).
Nitrogen accumulation (g N m −2 ) among crop and weed plant biomass components is presented in Table S3.At SITE A , N capture in spring barley was on average 20% greater in N 100 plots compared to N 50 .However, statistically significant differences were not present among other plant biomass components samples from SITE A (i.e., L. multiflorum, S.alba, and ambient weeds).At SITE B , winter wheat within the SA shoes treatment had 38% more g N m −2 than SA tines+shoes .Across slurry application methods, winter wheat N accumulation was 65% and 68% greater in N 100 plots compared to N 50 at SITE B and SITE C , respectively.A greater degree of difference in N accumulations for N 100 versus N 50 among winter wheat compared to spring barley may be linked to the 49% increase in achieved NH 4 -N application rates at SITE B and SITE C versus SITE A .At SITE B , N in P. annua and ambient weed components were 90% and 335% greater at N 100 than N 50 , respectively, while no differences among these groups were present at SITE C .
Grain quality measures were somewhat affected by treatment variables across siteyears.At SITE A , grain protein (%) was on average 13% greater in N 100 than N 50 plots; also, grain protein was on average 7% and 8% greater in SA shoes compared to SA tines+shoes and SA hose , respectively (Table S4).At SITE B and SITE C , on average in N 100 plots, grain bulk density was 1% and 2% greater, grain protein was 7% and 11% greater, and thousand kernel weight was 3% and 5% greater than N 50 at SITE B and SITE C , respectively.Table 4.The effect of nitrogen rate (N 100 and N 50 ) and slurry application method (SA hose , SA shoes , and SA tines+shoes ) on intra-row taprooted Sinapis alba (SITE A ) and Matricaria biomass (SITE B and SITE C ). Data from each site were analyzed separately.A connecting letters report displays differences among treatments, whereby means not sharing a similar letter are significantly different (p < 0.05).Significant main effects and interactions are highlighted in bold text (p < 0.05).Standard errors are displayed in parentheses.

Ammonia Emissions
In experiment EXP NH3.A , there were significantly lower cumulative emissions following slurry application by SA tines+shoes and SA shoes compared to SA hose (p = 0.0007).However, in experiments EXP NH3.B1 and EXP NH3.B2 , there were no significant differences between the treatments.Cumulative emissions can be seen in Table 7, and flux over time can be seen in Figure S3.
Table 7. Cumulative NH 3 emission as a percentage of applied total ammoniacal N (TAN) and g N m −2 113 h after field application of 22.8 tonne ha −1 pig slurry to a spring barley crop on loamy sand (EXP NH3.A at SITE A ) and 30 tonne ha −1 pig slurry to a winter wheat crop on sandy loam (EXP NH3.B1 and EXP NH3.B2 at SITE B ).A connecting letters report displays differences among treatments, whereby means not sharing a similar letter are significantly different (p < 0.05).Standard errors are displayed in parentheses.Following application, slurry surface pH showed the same trends for all three application methods (Figure S2).Within the first two hours, pH increased approximately one unit; thereafter, it dropped below the initial value.At the time of final pH measurements (3.5-4.5 h after application), no liquid was observed at the soil surface.

Exposed Surface Area
Slurry ESA immediately after application, and 2.66 h after application, differed among all treatments.Both initial ESA (t = 0 h) and final ESA (t = 2.66 h) were greatest in the SA hose treatment, second largest in SA shoes , and were lowest in SA tines+shoes (Table 8; Figure 1).Initial ESA at t = 0 h was 3.3 and 2.6 times greater in SA hose and SA shoes compared to SA tines+shoes , respectively.Final ESA at t = 2.66 h was 3.5 and 3.0 times greater in SA hose and SA shoes compared to SA tines+shoes , respectively.2) and ( 3) to analyze EXP ESA.B data.Average slurry exposed surface area (ESA, m 2 m −2 ) immediately after (t = 0 h) and 2.66 h after (t = 2.66 h) the application of 30 tonne ha −1 of pig slurry to winter wheat in sandy loam soil are also estimated and presented for three application methods (SA hose , SA shoes , and SA tines+shoes ).Lower and upper 95% confidence interval limits are shown in parentheses.A connecting letters report displays differences among parameters and estimates, whereby means not sharing a similar letter are significantly different (p < 0.05).

Discussion
Achieved winter wheat yields were low compared to the regional average.Within N100 plots, yields were 29% and 25% less than the 2020 mean reported for Zealand, Denmark, at SITEB and SITEC, respectively [37].Spring barley yields, however, were comparable in 2019; across N100 plots, achieved yields were only 6% lower than the reported average for West Jutland, Denmark [37].While placing slurry within a growing spring barley crop is not a recommended practice, the application of 100 kg N ha −1 is representative of a typical rate for organic farms, with conventional rates ranging from 140 to 150 kg N ha −1 .Reduced winter wheat yields are partly explained by having employed low N rates throughout the experiment; recommended rates for winter wheat range between 170 and 200 kg N ha −1 [38].In addition, splitting fertilizer placement into two springtime events, typically applying a mineral fertilizer first and slurry second, has shown to improve yields.Compared to a single application event, Olesen et al. [39] report

Discussion
Achieved winter wheat yields were low compared to the regional average.Within N 100 plots, yields were 29% and 25% less than the 2020 mean reported for Zealand, Denmark, at SITE B and SITE C , respectively [37].Spring barley yields, however, were comparable in 2019; across N 100 plots, achieved yields were only 6% lower than the reported average for West Jutland, Denmark [37].While placing slurry within a growing spring barley crop is not a recommended practice, the application of 100 kg N ha −1 is representative of a typical rate for organic farms, with conventional rates ranging from 140 to 150 kg N ha −1 .Reduced winter wheat yields are partly explained by having employed low N rates throughout the experiment; recommended rates for winter wheat range between 170 and 200 kg N ha −1 [38].In addition, splitting fertilizer placement into two springtime events, typically applying a mineral fertilizer first and slurry second, has shown to improve yields.Compared to a single application event, Olesen et al. [39] report the three-year yield average of winter wheat increasing by 2 to 4% when an equivalent N rate is split.Reduced N rates of 50 and 100 kg N ha −1 and one application event were chosen to enhance resource competition among crop and weeds, and to simplify the interpretation of treatment outcomes.
Nitrogen rate had the most consistent effect on crop growth across site-years.Crop biomass increased with elevated N input at all sites, with winter wheat demonstrating a stronger response compared to spring barley.This result is logical, given that postemergence fertilizer placement is ideal for meeting N demand of winter wheat [40], while pre-emergence fertilizer placement is best suited for spring barley [41].At SITE B and SITE C , the accumulation of N in winter wheat biomass was consistently greater in N 100 plots compared to N 50 , as were crop yield, grain protein, bulk density and TKW.However, a yield response was not present in spring barley, and among quality parameters, only grain protein was greater at N 100 .The absence of yield response could be due to the late application of slurry.Improved N capture in spring barley helps to explain the increase in crop biomass and grain protein at SITE A .By comparing the difference in N uptake between N 50 and N 100 , apparent N recovery (ANR) of NH 4 -N in applied slurry could be estimated.When Pedersen et al. [8] applied a mineral fertilizer on the same soil type ANR was 80% for spring barley.In spring barley at SITE A , average ANR was 41%, and in winter wheat at SITE B , average ANR was only 47%, indicating a significant loss of NH 3 .Corresponding NH 3 emission rates observed at SITE A and SITE B ranged from 9 to 17% and 25 to 42%, respectively.At SITE C , ANR was 73% indicating only small losses of NH 3 ; high utilization of NH 4 -N, despite some NH 3 loss, could be attributed to the mineralization of organic N in the pig slurry [8].
It is important to note that the field conditions at SITE A differed considerably from SITE B and SITE C .Soil texture at SITE A was lighter, possessing a higher percentage of sand (Table 1), and because spring barley was employed as the test crop, field cultivation preceding sowing further loosened the soil shortly before slurry was applied.As was hypothesized, standard trailing hoses were found to produce the greatest total NH 3 emissions at SITE A ; however, the effect of combining tine and trailing shoe did not reduce emissions compared to trailing shoes alone.This lack of difference is partially explained by the absence of a surface crust on light loose soils, limiting the functional difference between SA shoes and SA tines+shoes treatments (see Section 2.1.1).Soil texture is known to affect NH 3 emissions; sandy soils possess a greater slurry infiltration capacity than heavier clay-type soils, leading to decreased emissions.Both sandy soils [25,42] and cultivation preceding slurry application are known to increase slurry infiltration capacity and decrease NH 3 emissions [9,43]; in the case of SITE A , reducing the functional difference between SA hose versus SA shoes and SA tines+shoes .
At SITE A , a decrease in atmospheric N loss associated with reduced NH 3 emissions in combination with targeted placement in SA shoes and SA tines+shoes treatments did not translate into increased N uptake in spring barley as expected.More invasive slurry placement strategies may inflict yield penalties by causing root damage; however, among the tactics causing some shallow soil disturbance (SA shoes and SA tines+shoes ), no yield losses were observed.It can be interpreted that substantial spring barley damages were not imparted by SA shoes or SA tines+shoes strategies.If significant damages had been acquired, yield loss would have to be offset by improved N capture.Instead, crop yield, N accumulation, biomass, and grain quality measures did not differ among slurry application treatments, with the only exception being that SA shoes produced higher grain protein compared to other application strategies.
Soils at SITE B and SITE C had a much higher clay content than SITE A (Table 1).Because winter wheat was employed as the test crop, a lack of cultivation in the spring and low precipitation preceding slurry placement resulted in heavily crusted soils (Figure S4).Trailing shoes did not operate as expected, riding on top of the hardened surface crust without creating a depression, therefore, decreasing the functional difference between SA hose and SA shoes treatments.Results from Pedersen et al. [23] are supportive of the notion that soil type can impact a slurry placement method's ability to reduce NH 3 emissions; improved abatement by trailing shoes compared to trailing hoses were observed in coarse sand only, whereas no effects of placement method were observed in heavier soils, including a loamy sand and sandy loam.Meanwhile, in the present study, the prototype combing the operation of tines and trailing shoes operated as intended, with the straight tine breaking through the soil crust and the shoe riding behind depositing slurry into the soil break created (see Section 2.1.1 and Video S2).Despite SA tines+shoes functioning as designed and in opposition to the expected outcome, no significant differences in cumulative NH 3 emissions were observed among slurry application strategies.Although, when ordered from highest to lowest emissions, slurry application treatments performed as expected in EXP NH3.B2 , with SA hose > SA shoes > SA tines+shoes , and deviated slightly from our hypothesis, in EXP NH3.A and EXP NH3.B1 , SA hose > SA tines+shoes > SA shoes .Accordingly, there were no significant differences in crop N uptake between application methods.
SA tines+shoes reduced initial (t = 0 h) and final (t = 2.66 h) slurry ESA in EXP NH3.B .Slurry was spread at SITE B and SITE C following a long period without precipitation; therefore, cracks that developed under dry conditions in heavy clay soils provided a pathway for liquid slurry to funnel into when deposited in a band at the soil's surface (Figure S4).While SA tines+shoe functioned as intended, producing the lowest initial and final ESA, it was not enough to reduce total NH 3 emissions (Table 7).Under field conditions where hardened crusted soils are present, SA tines+shoes shows promise for reducing ESA and limiting NH 3 volatilization compared to SA shoes and SA hose , if large cracks have not developed providing an alternate route for rapid infiltration using simpler tools.However, the formation of large fissures in clay soils is common especially in overwintered crops following the freeze-thaw cycles of winter [44].When using surface application under these conditions, a larger slurry ESA may, in fact, increase the infiltration rate by broadening the surface area whereby liquid slurry is able to access and infiltrate into surficial cracks.Thus, an important subject of future research would be to investigate the use of slurry placement strategies under a variety of soil conditions.
It was hypothesized that crop damages caused by the shallow disturbance of SA shoes and SA tines+shoes treatments would be offset by improved N uptake, resulting in similar yields across slurry application methods.At SITE B , greater N uptake was observed for SA shoes compared to SA tines+shoes , and while this translated into increased biomass in SA shoes , grain yield effects were not observed among treatments.At SITE C , SA hose produced greater yields than SA tines+shoes , but differences in N uptake and winter wheat biomass were not observed among slurry application strategies.Again, there is little evidence to suggest that the SA shoes or SA tines+shoes treatments inflicted significant crop damages, and contrary to the expected outcome, at SITE B , SA tine+shoe reduced N uptake in winter wheat.As expected, slurry placement strategies did not have an effect on grain quality parameters at SITE B or SITE C , and ambient weed biomass did not differ among slurry application treatments.
Weed effects were observed at SITE B only; increasing N rate resulted in greater N accumulation in P. annua and other ambient weeds; P. annua biomass was also greater in N 100 plots.The observation of both positive and null effects across site-years may be accounted for by differences in weed species' response to N rate [45].The dominant ambient weed at SITE B was G. aparine, a species whose growth is consistently enhanced by increasing N [46]; while ambient weed biomass at N 100 was 2.9 times greater than N 50 at SITE B , significant differences were not observed due to high variability.As N rate increases, a positive growth response for P. annua has also been documented [47,48].However, it is surprising that no P. annua response was observed at SITE C, where P. annua accounted for 63% of the ambient weed population.
Weed species were expected to respond differently to slurry placement treatments tested; however, this effect was never observed.In studies where species-specific effects are present, differences in slurry placement depth among treatments are greater; for example, when depositing slurry at the soil's surface is compared to deep injection.A possible explanation for our having observed no response difference among shallow-rooted and taprooted weeds is that those treatments evaluated do not differ enough in slurry placement depth, either interacting shallowly or on top of the soil's surface.Upon injecting liquid fertilizer to a depth of 5 cm, Melander et al. [16] observed improved growth in deeply rooted ambient weed, Tripleurospermum indorum, and surrogate weed, Sinapis alba, as well as suppressed growth in shallowly rooted weed, Stellaria media.Reduced growth in S. media is likely due partly to having placed fertilizer below the weed's rooting zone, leading to decreased N capture, while improved growth in T. indorum and S. alba may be linked to increased N capture from having deposited fertilizer within these weeds' rooting zone.Results from Petersen [13] support this notion; when injecting slurry to a depth of 5 to 15 cm was compared to broadcast spreading followed by incorporation, injection resulted in a 58% increase in ambient weed N capture.In the present study, SA shoes and SA tine+shoes were expected to improve slurry infiltration into the soil column, resulting in increased N accumulation by taprooted weeds.The SA hose treatment places slurry at the soil's surface and was expected to have a slower infiltration rate, thus increasing accumulation among shallow rooted weeds.In opposition to our hypotheses, increased N uptake and plant biomass were not observed among shallowly rooted surrogate weeds, L. multiflorum or P. annua, for the SA hose treatment, nor among taprooted surrogate weeds, S. alba or Matricaria, in plots receiving SA shoes or SA tines+shoes treatments.Results from the present study indicate that compared to SA hose , slurry application with SA shoes or SA tines+shoes have neither positive nor negative implications for intra-row weed management within the hoed cereal system.
Inter-row hoeing was performed throughout the present study so that crop and weed response to slurry placement strategies could be assessed under contemporary management practices.Hoeing is an effective management tactic [28], and the advance of precision guidance technologies has eased implementation and increased adoption within the last decade [49].It may be hypothesized, however, that weed response to slurry placement treatments may become less pronounced when inter-row weeds are controlled, and assessments focus solely on the response of intra-row weeds.If inter-row weeds are not controlled, spreading slurry in a wide band between crop rows becomes much less efficient than when slurry is purposefully and precisely placed close to the crop row; in this scenario, weeds growing in the inter-row zone will have greater access to available N by proximity when compared to crop plants growing in the intra-row.Therefore, if weed effects are measured across both the inter-and intra-row zones, we would expect to see significant differences between banded and targeted placement strategies [14,15].In contrast, when hoeing is performed, and only intra-row weeds remain, the response of crop and weeds to slurry placement strategy becomes more reliant on slurry infiltration and targeted placement within the soil column rather than targeted slurry placement on the soil's surface [16].While it was expected that SA shoes and SA tines+shoes would have greater effects on weeds throughout the present study, our implementation of inter-row hoeing and focus on intra-row weeds may also help to explain the general absence of weed effects.

Conclusions
Observations made in the present study illustrate the considerable impact that soil conditions have on slurry deposition outcomes.In the presence of soils possessing a hard surficial crust, the successful operation of a prototype device was observed, whereby the combined operation of rigid tines and trailing shoes (SA tines+shoes ) reduced the initial and final slurry-exposed surface area when compared to standard trailing hoses (SA hose ) and trailing shoes alone (SA shoes ).Thus, the present study provides proof of concept for the prototype design in terms of developing a low draft, post-emergence slurry application device that is capable of successfully operating as intended in heavily crusted soils.However, the slurry placement methods tested were evaluated under the pretext that the adoption of new technologies must be justified by reduced ammonia (NH 3 ) emissions and improved crop growth.Therefore, the results ultimately do not warrant advocating the use of SA shoes or SA tines+shoes over SA hose across a range of soil textures and conditions.No significant effects on crop or weeds were detected among treatments, and compared to SA hose , NH 3 emissions were reduced by SA shoes and SA tines+shoes in only one of three site-years when operating in loose sandy soil.Therefore, deeper slurry placement within the soil column is likely necessary to consistently mitigate NH 3 emissions and illicit a positive crop response.This information may be of use to equipment manufacturers and researchers who wish to further develop and evaluate post-emergence slurry placement strategies.Future research on this topic should be conducted across a wide range of soil conditions or be considered site-specific.

Supplementary Materials:
The following are available online at https://www.mdpi.com/article/10.3390/agronomy12102441/s1, Table S1: Summary of dates, crop growth stages, and surrogate weed growth stages at the time of performing EXP CROP+WEED field operations and in-field data collection.Table S2: The effect of weed management treatment (WMT weedy and WMT weed-free ), nitrogen rate (N 100 and N 50 ), and slurry application method (SA hose , SA shoes , and SA tines+shoes ) on spring barley biomass (SITE A ) and winter wheat biomass (SITE B and SITE C ). Data from each site were analyzed separately.A connecting letters report displays differences among treatments, whereby means not sharing a similar letter are significantly different (p < 0.05).Significant main effects and interactions are highlighted in bold text (p < 0.05).Standard errors are displayed in parentheses.Table S3: The effect of nitrogen rate (N 100 and N 50 ) and slurry application method (SA hose , SA shoes , and SA tines+shoes ) on nitrogen accumulation within crop and intra-row weed plant biomass components; Hordeum vulgare, Sinapis alba, Lolium multiflorum (SITE A ) Triticum aestivum, Matricaria, Poa annua (SITE B and SITE C ) and other ambient weeds (SITE A , SITE B , and SITE C ). Data from each site were analyzed separately.A connecting letters report displays differences among treatments, whereby means not sharing a similar letter are significantly different (p < 0.05).Significant main effects and interactions are highlighted in bold text (p < 0.05).Standard errors are displayed in parentheses.Table S4: The effect of weed management treatment (WMT weedy and WMT weed-free ), nitrogen rate (N 100 and N 50 ), and slurry application method (SA hose , SA shoes , and SA tines+shoes ) on spring barley (SITE A ) and winter wheat (SITE B and SITE C ) grain bulk density, protein, and thousand kernel weight.Data from each site were analyzed separately.A connecting letters report displays differences among treatments, whereby means not sharing a similar letter are significantly different (p < 0.05).Significant main effects and interactions are highlighted in bold text (p < 0.05).Standard errors are displayed in parentheses.Figure S1: Prototype combining the operation of a tine and the Bomech trailing shoe, used in SA tines+shoes plots at SITE A . Figure S2: Slurry surface pH following field application of 30 tonne ha −1 pig slurry to a winter wheat crop growing on sandy loam (EXP NH3.B1 and EXP NH3.B2 ).Standard errors are displayed in bands.Figure S3: Ammonia emissions following field application of 22.8 tonne ha −1 pig slurry to a spring barley crop on loamy sand (EXP NH3.A ) and 30 tonne ha −1 pig slurry to a winter wheat crop on sandy loam (EXP NH3.B1 and EXP NH3.B2 ).Standard errors are displayed in bands.

2. 3 .
Exposed Surface Area (EXP ESA ) 2.3.1.Experimental Design, Treatment Specifications, and Data Collection Positioned next to EXP CROP+WEED.B and performed in conjunction with the EXP NH3.B2

Figure 1 .
Figure 1.Fitted exposed surface area (ESA) curves overlaid with ESA raw data following field application of 30 tonne ha -1 pig slurry to a winter wheat crop on sandy loam soil (EXPESA.B).Standard errors for raw data are displayed in bands.

Figure 1 .
Figure 1.Fitted exposed surface area (ESA) curves overlaid with ESA raw data following field application of 30 tonne ha −1 pig slurry to a winter wheat crop on sandy loam soil (EXP ESA.B ).Standard errors for raw data are displayed in bands.
Figure S4: Image of soil conditions at SITE B and SITE C , featuring substantial crusting and cracking.Video S1: Slow motion video of the Vogelsang trailing shoe operating in SA shoes plots at SITE B and SITE C .Video S2: Slow motion video of the prototype combining the operation of a tine and the Vogelsang trailing shoe operating in SA tines+shoes plots at SITE B and SITE C .Video S3: Slow motion video of hoses operating in SA hose plots at SITE A , SITE B , and SITE C .Author Contributions: Conceptualization, M.R.M., J.P. and B.M.; methodology, M.R.M., J.P., B.M. and T.N.; validation, J.P.; formal analysis, M.R.M. and J.P.; investigation, M.R.M. and J.P.; data curation, M.R.M. and J.P.; writing-original draft preparation, M.R.M.; writing-review and editing, M.R.M., J.P., P.S., B.M. and T.N.; visualization, M.R.M. and J.P.; supervision, P.S., B.M. and T.N.; project administration, T.N.; funding acquisition, B.M., T.N.M.R.M. was responsible for implementing and analyzing experiments on crop and weed response.J.P. was responsible for implementing and analyzing experiments on ammonia emissions and slurry exposed surface area.All authors have read and agreed to the published version of the manuscript.Funding: This study was part of the NUGA project No. 3400-16-1112 funded by the Green Growth and Development Programme (GUDP) under the Ministry of Environment and Food of Denmark.

Table 1 .
Overview of experiments performed in 2019 and 2020 across sites.

Table 2 .
Slurry properties and application rate used in EXP NH3 and EXP ESA .Standard errors are displayed in parentheses.

Table 3 .
The effect of weed management treatment (WMT weedy and WMT weed-free ), nitrogen rate (N 100 and N 50 ), and slurry application method (SA hose , SA shoes , and SA tines+shoes ) on spring barley grain yield (SITE A ) and winter wheat grain yield (SITE B and SITE C ), presented at 15% moisture content.Data from each site were analyzed separately.A connecting letters report displays differences among treatments, whereby means not sharing a similar letter are significantly different (p < 0.05).Significant main effects and interactions are highlighted in bold text (p < 0.05).Standard errors are displayed in parentheses.

Table 5 .
The effect of nitrogen rate (N 100 and N 50 ) and slurry application method (SA hose , SA shoes , and SA tines+shoes ) on intra-row fibrous-rooted Lolium multiflorum (SITE A ) and Poa annua biomass (SITE B and SITE C ). Data from each site were analyzed separately.A connecting letters report displays differences among treatments, whereby means not sharing a similar letter are significantly different (p < 0.05).Significant main effects and interactions are highlighted in bold text (p < 0.05).Standard errors are displayed in parentheses.
a Data underwent a square root (x + 1) transformation before analysis; transformed means (Trans.)and standard errors are presented along with back-transformed means (Back-trans.).Treatment comparisons made throughout the text of the present article are calculated from back-transformed means.b Lolium multiflorum density (no.m −2 ) was used as a covariate in the analysis for SITE A , and Poa annua density (no.m −2 ) was used in the analysis for SITE B and SITE C .

Table 6 .
The effect of nitrogen rate (N 100 and N 50 ) and slurry application method (SA hose , SA shoes , and SA tines+shoes ) on intra-row ambient weed biomass (SITE A , SITE B and SITE C ). Data from each site were analyzed separately.A connecting letters report displays differences among treatments, whereby means not sharing a similar letter are significantly different (p < 0.05).Significant main effects and interactions are highlighted in bold text (p < 0.05).Standard errors are displayed in parentheses.aAmbient weed density (no.m −2 ) was used as a covariate in the analysis.Surrogate weeds, Sinapis alba, and Lolium multiflorum, were not included in the ambient weed density measure for SITE A , and ambient weeds, Matricaria and Poa annua, were not included for SITE B and SITE C .

Table 8 .
Parameter estimates used in Equations (