Kura Clover Living Mulch: Spring Management Effects on Nitrogen

Kura clover living mulch (KCLM) systems have the potential to provide ecosystem services in intensively managed cropping systems while supplying soil mineral nitrogen (N) to the growing cash crop. Living mulch management relies on strong spring suppression to reduce competition between vigorous kura clover and emerging row crop seedlings, but standard suppression management practices utilize widely different modes of action. The objective of this research was to gain insight into the impact of common KCLM management practices on early season N dynamics. Kura clover was mowed, and residue was either harvested or returned before rows were established via strip tillage or banded herbicide. Soil and gaseous N pools were monitored for 12 weeks post initial application of suppression management treatments. An enrichment factor (EF) approach was utilized to compare N pools under managed treatments relative to an unmanaged clover control. Strip tillage increased soil N by 300%, while banded herbicide row establishment increased soil N by 220% relative to the unmanaged control. Pre-plant clover harvest reduced short term soil NO3–N, but during later time intervals there was no relationship between residue management and soil N. We conclude that, for the dual goals of maintaining clover perenniality while providing greater soil N enrichment, strip tillage is superior to band herbicide for row establishment. Additionally, pre-plant clover harvest may open opportunities for dual harvests in a single growing season, increasing economic return while maintaining in-season N contributions from the living mulch.


Introduction
Kura clover (Trifolium ambiguum M. bieb.) is a rhizatomous persistent perennial legume forage crop native to the mountainous regions of eastern Europe [1].Its dense root and rhizome system allows for frequent defoliation and vegetative repropagation in harsh environments [2].This extreme persistence, along with its shade, drought, flood, and cold tolerance, low-growing habit, and nitrogen-rich biomass has prompted researchers to investigate its use as a living mulch in cropping systems of the Upper Midwestern U.S. [3,4].
Kura clover living mulch (KCLM) systems have been recognized for their soil and water conservation benefits, including reduced soil erosion [5], increased water infiltration (Baker, unpublished data), and reduced nitrate leaching [6] and residual N [7], compared to conventionally managed crop rotations.These environmental benefits, however, are often accompanied by reduced yields of the main cash crop [4][5][6][7][8][9][10].Competition between the living mulch and emerging row crop Agronomy 2019, 9, 69 2 of 14 seedlings is a primary factor for delayed crop development and yield loss [4,10,11].This factor has led to the development of more aggressive suppression management practices during the critical spring establishment period [4,8].
Vigorous spring clover regrowth stores large amounts of organic N in protein-rich biomass early in the growing season.When spring growing conditions allow, clover may be harvested as forage [9], but most often, residue is returned via suppression management, and low C/N clover tissues are readily mineralized upon senescence and incorporation [8,12].Clover suppression techniques often utilize combinations of pre-plant mowing [8,9], broadcast chemical suppression [5,6,10,11], strip-tillage [12,13], and/or chemical banding [4,7,14] before and after the planting date.These management techniques vary in the amount of clover disturbance, the degree of incorporation, and the spatial organization of disturbed residues.
Methods of clover suppression and residue incorporation impact the resulting soil environment, which plays an important role in the spatiotemporal mineral N supply from biomass decomposition [15,16].Kura clover living mulch systems could be better utilized if management techniques were designed and chosen based on organic N cycling dynamics.Enhanced understanding of KCLM suppression management techniques may reveal agronomic benefits, such as a reduction in the fertilizer N requirement for the cash crop, while achieving soil and water conservation benefits.
This experiment branches from a two-year N management study for continuous corn in KCLM.Preliminary data from the first year of these experiments suggested increased soil N after herbicide suppression management, leading to our hypothesis that clover disturbance and suppression increases the reactive N supply in KCLM systems [17].The objective of this study was to compare N availability and loss pathways from commonly used KCLM suppression management practices.Understanding nitrogen dynamics and clover recovery after spring clover suppression management could facilitate the design of more stable, resilient, and beneficial companion cropping systems for the Upper Midwestern corn belt.

Site Description and Experimental Design
A field study was conducted from 29 May through 22 August, 2018 to investigate spatiotemporal N dynamics in KCLM systems after spring agronomic management.Plots were located at the University of Minnesota Research and Outreach Center in Rosemount, MN (44.73 • N, 93.09 • W) on a Waukegan silt loam (fine-silty over sandy or sandy-skeletal, mixed, superactive, mesic Typic Hapludolls).Soils from an adjacent experiment were grid sampled to 0.3 m (n = 64) and contained an average of 20.5 g kg −1 organic carbon and 5.7 pH in KCl [17].Endura kura clover was seeded at 11 kg ha −1 in 2006 and used as a living mulch for row crop production from 2008-2009, rhizomes were dug up with a potato digger for vegetative repropagation in 2010 [18], row crop production commenced in 2011-2014, and three hay cuttings and one seed harvest were taken from 2015-2017.In 2015, P and K were applied according to soil test values.
Two main clover management factors were examined: (i) seed-row establishment, with the row prepared either mechanically by rotary strip-tillage ('T'), or chemically using a banded herbicide or 'kill' ('K') application, and (ii) residue management, with mowed residue removed from ('Rm') or returned to ('Rt') the plots, resulting in four residue-row treatment combinations, i.e., T/Rm, T/Rt, K/Rm, and K/Rt (Table 1).An additional unmanaged control treatment, which was not mowed and had no row establishment, was also examined.Plots were arranged in a randomized complete block design with three replications of the five treatments and were 4.7 m (6 rows) wide by 7.6 m.In addition, the four plots within each block receiving row and residue management were split by 'row' and 'interrow' zones after row management treatments were applied (Figure A1).Clover was mowed to 50 mm on 29 May, and on 31 May cut clover residue was raked and removed.Rows were established on 4 June.Kill treatments received N-(phosphonomethyl) glycine (glyphosate) at 9.35 L a.e.ha −1 applied with a walk-behind sprayer unit in 0.3-m bands spaced every 0.76 m.Strip-till treatments were tilled using a rotary zone tillage tool (Northwest Tillers, Yakima, WA, USA) which created 0.3-m wide strips on 0.76-m intervals.An additional 2.3 L a.e.ha −1 of glyphosate was broadcast on managed plots (all treatments except the control) on 22 June to inhibit clover regrowth [5,6,10].Row-crops were not seeded into plots after row establishment so that the measured soil N pools could be isolated from soil N uptake by the main cash crop.

Soil N
Soils were collected bi-weekly from 17 May to 5 July and on 17 July, 6 August, and 22 August.Samples were taken at three depths (0-50, 50-150, and 150-300 mm) with a 20-mm i.d.coring device.Sample location was random in the unmanaged control and at the center of a randomly selected pair of row and interrow zones in managed plots.Samples were weighed and homogenized before 38 mL of 2 M KCl solution was added to a 10 g subsample of wet soil and shaken at 120 rpm for 1 h.Soil slurry was filtered through 11 µm filter paper and the extract was analyzed for NO 3 -N and NH 4 -N concentrations with the Greiss-Ilosvay with cadmium reduction and the sodium salicylate-nitroprusside methods, respectively, modified for flow-through injection analysis (Lachat, Loveland, CO, USA) [19].Cadmium reduces NO 3 -N to NO 2 -N, which is detected by the Greiss-Ilosvay method, therefore reported values for NO 3 -N are the sum of NO 2 -N and NO 3 -N.A second subsample of 5 g wet soil was oven dried for at least 24 h and weighed to determine gravimetric water content.Core volume and adjusted dry sample mass were used to calculate bulk density which was used to convert soil N concentrations on a soil mass basis (e.g., mg N g −1 ) to a per hectare basis (e.g., kg N ha −1 ) by depth interval.

Soil-Atmosphere Gas Exchange
Nitrous oxide emissions were measured bi-weekly from 24 May-5 July using non-steady-state chambers [20].Chamber bases measuring 0.7 m long × 0.36 m wide × 0.1 m deep were installed by trenching base dimensions with an electric chainsaw and pressing the acrylic base at least 50 mm into the ground.Bases were placed randomly within unmanaged control plots and randomly in managed treatments with the condition that the base spanned one row-width, containing equivalent ratios of row and interrow zones to field scale proportions.Bases in non-tilled treatments were installed 48 h before initial sampling to mitigate high gas flux from disturbance while bases in tilled treatments were installed immediately after tillage and hours before sampling to capture the effect of soil mixing on soil-to-atmosphere gas emissions, as done previously [21,22].Chamber displacement by clover biomass was calculated using fresh clover moisture content determined via oven-drying and density obtained via water displacement so that fresh clover volume could be calculated from in-season dry biomass measurements.The chamber volume used in the flux calculations was adjusted by subtracting the wet clover volume from the total above ground chamber volume [23].Biomass that exceeded chamber top height was folded into each top upon chamber placement [23].Atmospheric N 2 O-N concentration was sampled from the top of the clover canopy from each control plot with a 12 mL polypropylene syringe.These samples represented the initial (time 0) measurement for all chambers within the corresponding block.Insulated and vented chamber tops were then placed and secured onto bases with binder clips.Gas samples were collected from each chamber at 20, 40, and 60 min after chamber top placement and the 12 mL samples were immediately transferred to glass vials sealed with butyl rubber septa.Samples were handled and analyzed according to [20] with a 5890A Gas Chromatography analyzer (Hewlett-Packard, Palo Alto, CA, USA) in conjunction with a 7000 Headspace Autosampler (Teledyne Tekmar, Mason, OH, USA).Gas fluxes were calculated using the restricted quadratic method, where quadratic regression is used unless the second derivative of the resulting quadratic regression function is greater than 0, in which case linear regression is used [24].
Ammonia (NH 3 -N) emissions were measured bi-weekly from 24 May to 5 July.Transparent chambers were assembled and modified based on a semi-open chamber design as described previously [25].Twenty milliliters of 0.5 M H 2 SO 4 , containing 2% v/v glycerol solution and a 25 × 3 × 200 mm polyurethane foam strip were placed in a 125 mL bottle so that the foam was saturated in the acid solution.The bottle was then placed at the base of each chamber and the foam was suspended vertically with the bottom of the strip placed in the excess acid solution.After 3-4 d, the foam strips were removed from the chambers and placed into their respective bottles for transport to the lab.Fresh strips and solution were then installed into each chamber.Acid glycerol solution was added to the used strip and excess solution to reach the initial 20 mL volume before dilution with 30 mL of H 2 O.The total 50 mL solution was shaken for 0.5 h before filtration through 11 µm filter paper and the extract was analyzed for NH 4 -N concentration using the sodium salicylate-nitroprusside method modified for flow-through injection analysis [26] (Lachat, Loveland, CO, USA).

Clover Sampling
Aboveground biomass was collected on 17 May, 31 May, 13 June, 21 June, 28 June, and 26 July.A 0.5 m 2 quadrat was placed in each plot, with the condition that it contained row and interrow zones in proportion to the field scale and was not previously sampled.Clover biomass was cut to 10 mm and collected from within each quadrat.Clover samples were dried at 60 • C for at least 3 days before dry mass was taken.Subsamples were pulverized and analyzed for total C and N concentrations using the Dumas dry combustion method with an elemental analyzer [27] (VarioMax, Elementar, Langenselbold, Germany).

Environmental Conditions
Soil moisture and temperature data were collected at the center of each zone in managed plots and randomly within the control.Sensors were installed vertically at the 76-mm depth and measurements were taken at 0.5-hr intervals.Single replicates from each zone in each treatment were monitored from 28 May-10 July.Daily minimum and maximum temperatures from each zone were used to calculate cumulative soil heat units with a minimum threshold value at 10 • C to account for limited microbial activity below this value [28,29].
Daily precipitation, minimum air temperature, and maximum air temperature were obtained from the National Weather Service Cooperative Observer Station no.217107 for the time period beginning on 20 April and ending on 22 August.Daily maximum and minimum air temperatures were averaged to obtain single daily average air temperature values and cumulative precipitation was calculated beginning on 20 April.Daily average air temperature and cumulative precipitation in 2018 were compared graphically to the 1981-2010 historical average.

Data Analysis
Soil NO 3 -N, NH 4 -N, and the sum of NO 3 -N and NH 4 -N (total inorganic N; TIN) concentrations in soil from sampled depth intervals were summed across depth intervals.The summed concentrations were plotted against time and trapezoidally-integrated to represent cumulative soil N availability across the entire analysis period [30][31][32][33].Individual N 2 O flux measurements were similarly time-integrated to determine cumulative emissions over the sampling period.Cumulative NH 3 -N emissions were determined by the summation of individual flux measurements, since these values represented the total cumulative flux between sampling dates.
Values of time-integrated NO 3 -N, NH 4 -N, TIN, and N 2 O-N and NH 3 -N emissions were analyzed using an enrichment factor (EF), which was calculated based on Equation (1): where C represents the measured variable and b represents the experimental block.The EF approach was based on calculations commonly used in 15 N isotope, contaminant, and mineral-ore analyses, where measured concentrations are compared to baseline or reference isotope, elemental, or mineral concentrations present in earth's atmosphere, soils, or crust [34][35][36].For this study, the EF represents the magnitude of the difference between N pools under KCLM management treatments relative to baseline values under unmanaged clover.Categories were defined for enrichment factor values, where, <100% represents depletion, 100 ≤ EF C < 200% represents slight enrichment, 200 ≤ EF C < 400% represents moderate enrichment, and EF C ≥ 400% represents high enrichment of the N pool variable in the managed treatment relative to the unmanaged clover.Statistical analysis utilized the MIXED procedure of SAS 9.4 at p ≤ 0.05 (SAS Institute, Cary, NC, USA).Scatter plots of predicted and residual values were evaluated for homogeneity of variance and normality [37] with the UNIVARIATE procedure of SAS; these requirements were met for all dependent variables.The data were organized into four groups that were subjected to separate statistical analyses, each during different time periods corresponding with management type, as follows (Supplementary Tables S1-S3).
Short-term soil N response to residue management: This analysis focused on the 6-d period between mowing (29 May) and application of row treatments (4 June), and aimed to evaluate the effects of the management 'system' (i.e., Rm, Rt, or control) on time-integrated NO 3 -N, NH 4 -N and TIN, and the effects of residue management (i.e., Rt or Rm) on the EFs corresponding to these soil N variables.
Zone-differentiated N response to residue and row management: This analysis examined the EF variables NO 3 -N, NH 4 -N, TIN, and NH 3 -N following row establishment (4 June) through the end of the experiment (22 August), and their relationship with residue management, row establishment (i.e., K or T), and zone (i.e., row or interrow).
Zone-weighted N response to residue and row management: This analysis addressed the effect of residue management and row establishment on the EF variables NO 3 -N, NH 4 -N, TIN, NH 3 -N, and N 2 O-N over the entire sampling period of the corresponding N variable.For this analysis, N variables observed in the row and interrow zones following row establishment were weighted by relative zone area before time-integration and calculation of the EFs.This analysis also examined clover biomass and biomass-N variables sampled on 28 June as affected by residue management and row establishment, since the majority of clover regrowth was expected to occur prior to row crop canopy closure.
Cumulative N gas emissions: This analysis evaluated the effect of the management system (i.e., K/Rt, K/Rm, T/Rt, T/Rm, or control) on N 2 O-N and NH 3 -N emissions occurring over the entire study period (24 May-5 July).For this analysis, NH 3 -N fluxes observed in the row and interrow zones in the managed plots were combined after weighting by relative zone area.This weighting procedure was not necessary for N 2 O-N because measurements were made across both zones with a single chamber.

Environmental Conditions
Environmental conditions during the study period were near the 1981-2010 climate average (Figure 1).The average daily temperature was 19.9 • C in 2018 and 17.6 • C from 1981-2010.Cumulative precipitation was 436 and 457 mm in 2018 and 1981-2010, respectively (Figure 1).chamber.Environmental conditions during the study period were near the 1981-2010 climate average (Figure 1).The average daily temperature was 19.9 °C in 2018 and 17.6 °C from 1981-2010.Cumulative precipitation was 436 and 457 mm in 2018 and 1981-2010, respectively (Figure 1).

Short-Term Soil N Response to Residue Management
Nitrate and TIN were significantly affected by the management system (Table 2).Timeintegrated soil NO3-N concentration following mowing and residue management was 105% greater in the residue returned (Rt) system than the unmanaged control while the residue removed (Rm) system was not different from the Rt or control systems (Figure 2a).Soil TIN concentrations in the Rt and Rm treatments were 87 and 49% greater than the control, respectively (Figure 2b).The EFs for all three soil N variables did not differ between the Rt and Rm systems in the period following mowing and preceding row treatment application.Soil TIN concentrations over this sampling period can be found in the Appendix (Figure A2).-----------EF, %-------

Short-Term Soil N Response to Residue Management
Nitrate and TIN were significantly affected by the management system (Table 2).Time-integrated soil NO 3 -N concentration following mowing and residue management was 105% greater in the residue returned (Rt) system than the unmanaged control while the residue removed (Rm) system was not different from the Rt or control systems (Figure 2a).Soil TIN concentrations in the Rt and Rm treatments were 87 and 49% greater than the control, respectively (Figure 2b).The EFs for all three soil N variables did not differ between the Rt and Rm systems in the period following mowing and preceding row treatment application.Soil TIN concentrations over this sampling period can be found in the Appendix A (Figure A2).-----------EF, %-------

Zone Differentiated N Response to Residue Management and Row Establishment
The row establishment by zone interaction was significant for NO3-N and TIN EFs following row establishment through the end of the study (Table 3).Nitrate and TIN in the row of the tilled (T) treatment were highly enriched compared to the unmanaged control, 675 and 479%, respectively, and were statistically greater than the row of the killed (K) treatment and the K and T interrow zones.The row in the K treatment was moderately enriched in NO3-N (333%) and was statistically greater than the interrow in the K and T treatments.Nitrate concentrations in the interrow of the K and T treatments were moderately (202%) and slightly enriched (183%), respectively, and were not different from each other.Soil TIN enrichment in the K interrow (188%) was not different from the K row (247%) or the T interrow (160%), but enrichment in the K row was greater than the T interrow (Figure 3).Soil TIN concentration over this sampling period can be found in the Appendix (Figure A3).

Zone Differentiated N Response to Residue Management and Row Establishment
The row establishment by zone interaction was significant for NO 3 -N and TIN EFs following row establishment through the end of the study (Table 3).Nitrate and TIN in the row of the tilled (T) treatment were highly enriched compared to the unmanaged control, 675 and 479%, respectively, and were statistically greater than the row of the killed (K) treatment and the K and T interrow zones.The row in the K treatment was moderately enriched in NO 3 -N (333%) and was statistically greater than the interrow in the K and T treatments.Nitrate concentrations in the interrow of the K and T treatments were moderately (202%) and slightly enriched (183%), respectively, and were not different from each other.Soil TIN enrichment in the K interrow (188%) was not different from the K row (247%) or the T interrow (160%), but enrichment in the K row was greater than the T interrow (Figure 3).Soil TIN concentration over this sampling period can be found in the Appendix A (Figure A3).N  --------* EF, %-------

Zone Weighted N Response to Residue and Row Management
Row establishment significantly affected EFs for NO3-N and TIN following mowing through the end of the experiment (Table 4).The EF for NO3-N was greater in the T treatment (304%) than the K treatment (220%), corresponding with greater enrichment of TIN in the T treatment (229%) than the K treatment (179%).Biomass and biomass-N were not affected by row or residue management.Soil TIN concentration and clover biomass over this sampling period can be found in the Appendix (Figure A4 and Figure A5).Zone significantly affected NH 3 -N emissions, with greater enrichment in the row (170%) than the interrow (120%), corresponding with cumulative NH 3 -N emissions of 0.68 and 0.45 kg ha −1 , respectively.

Zone Weighted N Response to Residue and Row Management
Row establishment significantly affected EFs for NO 3 -N and TIN following mowing through the end of the experiment (Table 4).The EF for NO 3 -N was greater in the T treatment (304%) than the K treatment (220%), corresponding with greater enrichment of TIN in the T treatment (229%) than the K treatment (179%).Biomass and biomass-N were not affected by row or residue management.Soil TIN concentration and clover biomass over this sampling period can be found in the Appendix A (Figures A4 and A5).-------* EF, %-------

Cumulative N Gas Emissions
Emission of N 2 O-N and NH 3 -N over the total study period (24 May-5 July) was not affected by management system (Table 5).

Discussion
Soil N concentration responded to mowing and residue management in the 6 d period between mowing and row establishment (Figure A2).The control treatment had less soil NO 3 -N than when residue was returned (Rt) and less soil TIN than when residue was returned or removed (Rm).In the context of utilizing kura clover as a living mulch to supply N for a cash crop, an absence of a residue management effect on the measured variables suggests that early surface residue removal does not influence soil N pools.Enrichment of soil NO 3 -N and TIN after row establishment was high in the T row, moderate in the K row, and slight in both T and K interrows.Likewise, zone-weighted EFs for soil NO 3 -N and TIN were significantly affected by row management, with greater enrichment in the T treatment than the K treatment.
Soil N not taken up by plants can be lost to environmentally damaging gaseous N pools; NH 3 -N, NO x -N, and N 2 O-N [38][39][40].Volatilized NH 3 -N was slightly enriched in both the row and interrow after row establishment and significantly greater in the row than the interrow.This result may suggest that elevated soil N in the row increased NH 3 -N volatilization; however, NH 4 -N was not significantly affected by zone, indicating that NH 3 -N was volatilized directly from senescing clover residue.Bursts of NH 3 -N emissions after chemical or mechanical senescence of legume crops have been reported previously [41][42][43].Despite enrichment in NH 3 -N volatilization in clover under living mulch management, the magnitude of emission was low in comparison to soil mineral N pools and annual NH 3 -N emission from local cropland [44].This may be due to soil acidity (5.7 in the surface layer), since ammonia volatilization is a pH-dependent reaction in which NH 4 -N is hydrolyzed when high concentrations of OH − are present [45].
In a previous study, N 2 O-N emission during the spring management period in a corn-soybean (Glycine max L. (Merr.))rotation was 98% greater in corn and 161% greater in soybean in a KCLM system compared to conventional management [12].Findings from the present study further highlight the potential for elevated N 2 O-N emissions during living mulch management.Further research is needed to better understand the mechanisms and factors governing partial denitrification and N 2 O-N emissions in KCLM systems, including the partitioning of organic and inorganic N sources of N 2 O-N, and the feasibility of applying management practices to reduce N 2 O emissions, for example, the use of nitrification inhibitors [46].
Greater enrichment of soil N with strip tillage row establishment suggests that soil N contribution from the living mulch is most influenced by biomass incorporation, where more intensive row establishment methods increase mineralizable biomass supply, soil-biomass contact, soil temperature, organic matter mineralization rate, nitrification rate, or a combination of these factors (Figure A6) [47,48].Despite more intensive management in the tilled treatment, biomass and biomass-N did not differ among the Rt, Rm, K, or T treatments by mid-summer, indicating that clover regrowth is similar among residue and row management treatments (Figure A5).These results suggest that zone tillage is preferable to banded herbicide for KCLM cropping systems to supply in-season N to the cash crop.Zone tillage increases soil temperature, soil contact with senescing biomass, and soil N contributions to growing cash crops [47].Removal of cut kura clover residues did not reduce soil N pools, supporting a spring forage harvest before seeding the primary cash crop when conditions allow.Clover harvest prior to planting the cash crop would probably be most feasible as haylage rather than bales due to the risk of slow field drying prior to baling, along with greater harvest losses and the increased time commitment of haying during the spring planting season [49][50][51], although, harvest method will be most influenced by the equipment and storage capabilities of the individual producer.The option to utilize kura clover as a forage crop may also offer flexibility within crop rotations when unfavorable field conditions, poor market outlooks, or market opportunities arise near the planting date.These results open new possibilities for optimization of KCLM management practices that might produce multiple harvests in a single growing season, maintain clover perenniality, reduce reliance on fertilizer N, increase system resilience, and improve soil health.
Concepts presented here improve the current understanding of legume-row crop intercropping systems, but further research is necessary to fully quantify N dynamics in KCLM production systems.These findings should be useful for optimization of spring agronomic management practices for KCLM production systems.Future research should address variation in system performance as influenced by kura clover stand age, cropping history, and extreme weather conditions for the development of stable and resilient KCLM cropping systems.

Figure 1 .
Figure 1.Daily average air temperature and cumulative precipitation from 20 April-22 August in 2018 and the 1981-2010 average.

Table 2 .
Treatment means and significance of soil N pools post residue management prior to row establishment.

Figure 1 .
Figure 1.Daily average air temperature and cumulative precipitation from 20 April-22 August in 2018 and the 1981-2010 average.

Table 2 .
Treatment means and significance of soil N pools post residue management prior to row establishment.

Figure 2 .
Figure 2. Time-integration (left axis) and enrichment factor (EF) (right axis) means and standard error for soil nitrogen variables by management system.(a) Soil NO3-N; (b) Soil TIN.Columns with the same letter are not significantly different at p ≤ 0.05.

Figure 2 .
Figure 2. Time-integration (left axis) and enrichment factor (EF) (right axis) means and standard error for soil nitrogen variables by management system.(a) Soil NO 3 -N; (b) Soil TIN.Columns with the same letter are not significantly different at p ≤ 0.05.
† Within a column for a given fixed effect, means followed by the same letter are not significantly different at p ≤ 0.05.* Enrichment factor of measured N variables.‡ Banded herbicide row establishment.§ Zone tillage row establishment.│ Residue returned after mowing.|| Residue removed after mowing.

Figure 3 .
Figure 3. Enrichment factor (EF) and standard error for soil NO3-N and soil inorganic N (TIN) following row establishment as a function of band herbicide (K) or zone till (T) row establishment and zone (row or interrow).Columns with the same letter are not significantly different at p ≤ 0.05.

Table 4 .
Treatment means and significance of soil and gaseous N pools over the entire sampling period.Clover biomass and biomass-N are from the 28 June sampling.

Figure 3 .
Figure 3. Enrichment factor (EF) and standard error for soil NO 3 -N and soil inorganic N (TIN) following row establishment as a function of band herbicide (K) or zone till (T) row establishment and zone (row or interrow).Columns with the same letter are not significantly different at p ≤ 0.05.

Table 4 .
Treatment means and significance of soil and gaseous N pools over the entire sampling period.Clover biomass and biomass-N are from the 28 June sampling.

Figure A1 .
Figure A1.Plot layout and design within a single replication.

Figure A2 .
Figure A2.Soil TIN concentration as a function of management system following residue management, where clover was returned (Rt) or removed (Rm) from the plot and preceding row establishment and with unmanaged clover (Control).

Figure A2 .
Figure A2.Soil TIN concentration as a function of management system following residue management, where clover was returned (Rt) or removed (Rm) from the plot and preceding row establishment and with unmanaged clover (Control).

Agronomy 2018, 8 , 15 Figure A3 .
Figure A3.Soil TIN concentration following row establishment as a function of banded herbicide (K) and zone tillage (T) row management and with unmanaged clover (Control).

Figure A4 .
Figure A4.Soil TIN concentration over the entire sampling period as a function of banded herbicide (K) and zone tillage (T) row management and with unmanaged clover (Control).

Figure A3 . 15 Figure A3 .
Figure A3.Soil TIN concentration following row establishment as a function of banded herbicide (K) and zone tillage (T) row management and with unmanaged clover (Control).

Figure A4 .
Figure A4.Soil TIN concentration over the entire sampling period as a function of banded herbicide (K) and zone tillage (T) row management and with unmanaged clover (Control).

Figure A4 .
Figure A4.Soil TIN concentration over the entire sampling period as a function of banded herbicide (K) and zone tillage (T) row management and with unmanaged clover (Control).

Figure A5 .
Figure A5.Clover biomass and standard error over the entire sampling period as a function of banded herbicide (K) and zone tillage (T) row management and with unmanaged clover (Control).

Figure A6 .
Figure A6.Non-replicated soil temperature heat units for residue return (Rt) and removal (Rm) after mowing management and banded herbicide (K) or zone tillage (T) row establishment and with unmanaged clover (Control) as a percentage of the T/Rm row.Heat units were the sum of the average maximum and minimum daily temperature measurements minus 10 °C.

Figure A5 .
Figure A5.Clover biomass and standard error over the entire sampling period as a function of banded herbicide (K) and zone tillage (T) row management and with unmanaged clover (Control).

Figure A5 .
Figure A5.Clover biomass and standard error over the entire sampling period as a function of banded herbicide (K) and zone tillage (T) row management and with unmanaged clover (Control).

Figure A6 .
Figure A6.Non-replicated soil temperature heat units for residue return (Rt) and removal (Rm) after mowing management and banded herbicide (K) or zone tillage (T) row establishment and with unmanaged clover (Control) as a percentage of the T/Rm row.Heat units were the sum of the average maximum and minimum daily temperature measurements minus 10 °C.

Figure A6 .
Figure A6.Non-replicated soil temperature heat units for residue return (Rt) and removal (Rm) after mowing management and banded herbicide (K) or zone tillage (T) row establishment and with unmanaged clover (Control) as a percentage of the T/Rm row.Heat units were the sum of the average maximum and minimum daily temperature measurements minus 10 • C.

Table 1 .
Chronological clover suppression management by management system.
-Within a column means followed by the same letter are not significantly different at p ≤ 0.05.* Enrichment factor of measured N variables.‡ Unmanaged clover.§ Residue returned after mowing.|| Residue removed after mowing.† Within a column means followed by the same letter are not significantly different at p ≤ 0.05.* Enrichment factor of measured N variables.‡ Unmanaged clover.§ Residue returned after mowing.|| Residue removed after mowing.

Table 3 .
Treatment means and significance of N pools following row management.
-Within a column for a given fixed effect, means followed by the same letter are not significantly different at p ≤ 0.05.* Enrichment factor of measured N variables.‡ Banded herbicide row establishment.§ Zone tillage row establishment.| Residue returned after mowing.|| Residue removed after mowing.
Within a column for a given fixed effect, means followed by the same letter are not significantly different at p ≤ 0.05.* Enrichment factor of measured N variables.‡ Banded herbicide row establishment.§ Zone tillage row establishment.| Residue returned after mowing.|| Residue removed after mowing.

Table 5 .
Treatment means and significance of N emission over the entire sampling period.Unmanaged clover.‡ Banded herbicide row establishment.§ Zone tillage row establishment.| Residue returned after mowing.|| Residue removed after mowing.