Annual canarygrass (Phalaris canariensis L.) is a specialty crop grown in Canada and the harvested grain is primarily used to feed wild and domesticated bird species. A field experiment was conducted at 5 locations in both 2012 and 2013 to study the response of annual canarygrass development and grain yield to the combined effects of fungicide (propiconazole + trifloxystrobin) and nitrogen (N) fertilizer, and to determine the minimum number of site years required to detect these effects. The experimental design was a split plot with fungicide application (none or fungicide) as the main plot and N application rate as the sub plot (10, 20, 30, 50, 70, 90 kg N·ha−1). There was a linear increase of 14.5% in grain yield as the rate of N fertilizer increased. The fungicide application increased the grain yield 8.5% by increasing kernel density from 15 197 kernels m−2 to 16 288 kernels m−2. There was no interaction between the N rate and fungicide application. The application of a fungicide did not increase the responsiveness of annual canarygrass to N fertilizer. The lack of an interaction between N and fungicide application indicates that these two practices can be managed independently by annual canarygrass producers. To optimize grain yield producers should apply 50 kg N·ha−1 and apply a fungicide to increase grain yield especially in regions where septoria leaf mottle is prevalent. The number of site years needed to consistently detect the impact of N and fungicide on the grain yield were 4 and 5 site years, respectively.
Introduction
Annual canarygrass (Phalaris canariensis L.), also called canary seed, is a cereal crop grown in the prairie region of Canada, which is the largest producer and exporter of annual canarygrass worldwide (Tridge 2020). In 2018, 100 000 hectares were harvested in Canada, producing around 146 000 tonnes of grain worth $73 730 000 (CAD) (Agriculture and Agri-Food Canada 2020; Tridge 2020). Currently, the major use of annual canarygrass is as a feed for both caged and wild bird species. The grain of annual canarygrass has higher concentrations of the eight essential amino acids than wheat (Triticum aestivum L.) or corn (Zea mays L.) and a higher level of protein then most comparable grains (Putnam et al. 1996). Due to its unique nutrient profile, it is starting to enter the food market for human consumption, as a gluten-free cereal. In 2016, Health Canada (2016) approved glabrous (hairless) annual canarygrass, a gluten-free grain for human consumption. Currently, small amounts are used in health food drinks produced and sold in Latin America.
Annual canarygrass originates from the Mediterranean region and is a grain crop that is most closely related to oat (Avena sativa L.) (Putnam et al. 1996). While it is currently grown in Europe, Asia, Australia, and the Americas, the biggest producer is Canada, specifically Saskatchewan (Tridge 2020). It is an annual crop that can reach between 50 cm and 100 cm tall with a compressed oval shaped panicle. This crop prefers clay soils with good moisture retention due to its shallow rooting depth (Putnam et al. 1996). Annual canarygrass is usually planted in May and takes 95–110 d to mature (May et al. 2012a).
Annual canarygrass growth, development, and grain yield can benefit from the application of nutrients, especially N and chloride (Cl−) (May et al. 2012a, 2012b). The positive impact of N fertilizer on annual canarygrass grain yield has been observed in several studies. A field trial at Indian Head from 1978–1983 found a curvilinear increase in grain yield as the rate of fertilizer N increased, with no effect on kernel weight or maturity (Holt 1988). The suggested rate of application fell between 50–75 kg N·ha−1 (Holt 1988). May et al. (2012a) reported that an increase in N fertilizer up to 78 kg·ha−1 resulted in an increase in grain yield. Both studies found that grain yield increased as N rate increased, but overall yield and the yield increase was lower (710 to 2200 kg·ha−1) compared with other cereal crops. May et al. (2004) reported that oat yields can increase to 4100 kg·ha−1 with 80 kg N·ha−1, a similar maximum N rate compared with annual canarygrass, but with a much higher grain yield. O’Donovan et al. (2015), found that barley (Hordeum vulgare L.) had a higher grain yield (4600 kg·ha−1) than annual canarygrass but required a higher N rate (110 kg N·ha−1) to maximize grain yield and Biswas and Ma (2016) reported that with 180 kg N·ha−1, corn yields can reach up to 14 190 kg·ha−1. When evaluating an appropriate N rate to recommend for annual canarygrass production it is important to consider the sensitivity of the N rate recommendation to changes in the price of crop and N fertilizer. This becomes even more important if in the future N fertilizer prices are elevated to reduce excess N fertilization and green house gas production (Meyer-Aurich et al. 2020).
As the supply of N to the crop is increased there are physiological structural changes that can occur in a cereal plant to increase its susceptibility to lodging including basal internodes with smaller diameters, fewer support roots per stem, and a lower root dry weight per stem and as lodging increases grain yield often decreases (Berry et al. 2004; Easson et al. 1993). May et al. (2012a) found that increased lodging as the N rate increased occurred in a few site years, but not consistently across years and locations. Also, a microclimate can be created by lodging that fosters disease development that further reduces leaf area as well as grain yield (Nakajima et al. 2008).
Annual canarygrass faces a number of biotic stresses during the growing season including weeds, aphids, and disease. These pests can lower yield and increase the production costs incurred to control these pests, thereby lowering the profitability of annual canarygrass. For annual canarygrass, the most relevant disease is septoria leaf mottle caused by Septoria triseti Speg. (Kaminski 1994). The fungus produces pale tan to gray lesions that appear on the leaf blades and sheaths (Berkenkamp 1989). First observed on annual canarygrass in 1988, septoria leaf mottle restricts photosynthetic capacity, which causes the plants to produce smaller kernels, lowering grain yield, and causing an economic loss to growers (Kaminski 1994). Field scale trials conducted between 2008 and 2013 at Indian Head found that annual canarygrass had a large and consistent grain yield response in every year to a fungicide application to control septoria leaf mottle (Holzapfel 2014). Holzapfel (2014) found that the increase in grain yield from a fungicide ranged from 4% to 67%, and the average response over the 6 years was 23%. In addition, he found that annual canarygrass was more responsive to a fungicide application than the other cereals tested including, oat, wheat and barley.
Previous research has demonstrated that annual canarygrass responds to both the application of a fungicide and the application rate of N fertilizer; however, the combined effect of a fungicide application and N fertilizer has not been studied (Holzapfel 2014; May et al. 2012a). The objectives of this study were to determine if a fungicide application and nitrogen management have any synergistic effects on the growth, development, grain yield and gross return of annual canarygrass and to determine the minimum number of site years required to detect these effects using the data from this study.
Materials and Methods
This study was conducted in 2012 and 2013 across 5 locations in Saskatchewan: Melfort (52°49′6.29″ N, 104°36′11.77″ W and elevation 478 m), Indian Head (50°31′58.48″ N, 103°39′2.91″ W and elevation 582 m), Swift Current (50°16′49.56″ N, 107°47′58.39″ W and elevation 748 m), Redvers (49°34′21.82″ N, 101°41′53.54″ W and elevation 594 m) and Scott (52°21′40.73″ N, 108°50′15.85″ W and elevation 655 m). The experimental design was a split plot with fungicide application (none or Stratego®, propicoazole 62 g a.i. ha−1 and trifloxystrobin 62 g a.i. ha−1) as the main plot and N application rate as the sub plot (10, 20, 30, 50, 70, 90 kg N·ha−1). A gap of at least one sub plot was placed between the unsprayed and sprayed main plots to ensure that spray drift did not reach the unsprayed main plots during fungicide application and that disease develop in the unsprayed main plots did not quickly spread back into the sprayed main plots. The applied N rates were not adjusted for the residual N level in the soil. There were four replicates at each location. Each site used their own research seeder to seed the experiment. Plot size ranged from 1.4 m × 5.5 m in Swift Current and Melfort to 4 m × 10 m in Scott. (Table 1). The experiments were seeded between 15 May and 2 June. CDC Bastia was the annual canarygrass cultivar used and was seeded at a rate of 35 kg·ha−1. All seed lots of annual canarygrass used in this study had a germination percentage above 85%. The source of N was urea with a N concentration of 46%. At each site, phosphorus (P), potassium (K), Cl and sulphur (S) were side banded during seeding to each plot to provide 20 kg P2O5·ha−1, 10 kg K2O·ha−1, 6.5 kg Cl·ha−1 and 10 kg S·ha−1. Plot drills using hoe type openers were used in a no-till cropping system. All fertilizer was side banded at seeding, placed 3 cm to 4 cm to the side and 7 cm below the seed. Glyphosate was applied before seeding at Indian Head at 360 g a.e.·ha−1. In crop broadleaf herbicides were determined for each seeding date and site according to weed species density. The fungicide application occurred after the flag leaf was fully emerged with a plot sprayer.
Table 1.
Agronomic information for each of the individual site years.
Data collection
Soil was collected before seeding at depths of 0–15 cm and 15–60 cm for determining the availability of nutrients N, P, and K. A sodium bicarbonate (NaHCO3) extraction procedure (Hamm et al. 1970) was used to estimate residual soil N (nitrate; NO3), P, and K. Plant density was measured approximately three weeks after seeding. The number of plants were counted at two separate locations in the plot and at each location plants were counted in 1 m of one row. Panicle density, plant height, and lodging were each measured during grain filling (Zadocks 73 to 83). For each plot, the number of panicles was measure at two separate locations in the plot and at each location panicles were counted in 1 m of one row. Plant height was measured at two locations in each plot prior to harvest. Lodging of plants was a visual rating on the entire plot prior to harvest on a scale of 1–10, where 1 represented full standing plants and 10 represented plants completely flat on the ground. Leaf disease severity caused by septoria leaf mottle was assessed using the McFadden scale (0 to 11), a quantitative ordinal scale, on the entire plant canopy before plant maturity was reached (Zadocks 80 to 83) (McFadden 1991). Grain yield was determined with a self-propelled plot combine harvesting a between 7.7 to 22.5 m2 of each plot depending on the site and year (Table 1). The harvest occurred between 15 Sept. and 16 Oct. (Table 1). Grain yield was corrected to a kernel moisture level of 13% (Canaryseed Development Commission of Saskatchewan 2021). Kernel weight was calculated by approximately taking a 40 mL volume of seed using a small measuring cylinder and placing it in an automated seed counter, which resulted in counting between 700 and 1000 kernels and expressed as grams per 1000 kernels. Kernel density was calculated by dividing the grain yield per square metre by the kernel weight of a single seed. Test weight was measured as specified by the Canadian Grain Commissions Official Grain Grading Guide (2016) using a cox funnel to fill a 0.5 L cylindrical cup and then weighting the cup (g 0.5 L−1) after a hardwood striker is used to scalp off the excess grain from the top of the cup.
Change in gross return
The gross return for each combination of N fertilizer cost and seed value was calculated by using the following equation:
Three N prices were chosen: 0.50, 0.65 and 1.00 $·kg−1. The 0.50 and 0.75 represents the current range in the cost of N fertilizer while the $1.00 kg−1 price represents an elevated price that might be established by government policy to reduce N fertilizer use to combat global warming (Meyer-Aurich et al. 2020). Three grain prices were chosen: 0.44, 0.66 and 0.88 $·kg−1. These prices cover the usual price range of the annual canarygrass crop price from low to high (Government of Saskatchewan 2020a). Together, the three N fertilizer costs and three crop prices created nine combinations. The change in gross return was calculated by taking the gross return of each N rate and subtracting the gross return at the lowest N rate of 10 kg·ha−1. This provided a measure of how gross return changed over the range of N rates used in this study.
Statistical analysis
Mixed model analysis of data was conducted using the PROC MIXED model procedure in SAS software (Littell et al. 2006). Nitrogen application rate and fungicide application were assigned as fixed rates while the effect of the replicate and site year was considered random. Linear and quadratic contrasts were used to determine the effect of N on the measured variable. Unequal spacing between N treatments were use, 10, 20 30, 50, 70, 90 kg N·ha−1, to ensure that changes at low N rates, 10 to 30 N·ha−1 were captured and then the interval was widened to increase the probability of exceeding the N required to maximize grain yield. The contrast coefficients were calculated in SAS using PROC IML (SAS Institute 2013). Statistical significance was declared at P < 0.05 for all analysis. The random effects were only examined when they had a significant effect and account for more than 10% of the total variation associated with the random effects.
A second analysis was conducted to assess the statistical power from this existing data set when a different number of site years and replicates were used for the analysis. More specifically, the number of site years examined was 4 to 9 (maximum number of site years from the existing dataset) (Stroup 1999; W.W. Stroup personal communication). Four to 10 replicates were considered in the analysis. An iterative process was set up to run a set of SAS code for a progressively greater number of site years. The first step included PROC SURVEYSLECT (SAS Institute 2013) to randomly select the site years from the existing data set. A simple random sampling method was used. The portion of the existing dataset associated with these sites was filtered according to the site years selected. The filtered dataset was analyzed using the PROC MIXED procedure of SAS (Littell et al. 2006; SAS Institute 2013) where site year, fungicide, N fertilizer rate were designated as fixed effects, and replicate was considered as a random effect. The covariance parameters from this analysis were saved into a dataset for future use. The filtered dataset was then used as a base dataset which included four to 10 replicates (Stroup 1999; W.W. Stroup, personal communication). A series of SAS data steps were then used to create these 7 separate datasets with varied replication. This dataset with varied replication was analyzed using the PROC MIXED procedure of SAS (Littell et al. 2006; SAS Institute 2013) where site year, fungicide, N fertilizer rate were designated as fixed effects, and replicate was considered as a random effect. This analysis was repeated for each level of replication and was seeded with the covariance parameters saved from the first PROC MIXED analysis. The parms statement, pdata option was used to seed these covariance parameters. The final step with each site year iteration was to calculate power (Littell et al. 2006). After the power calculation, the process began again with an analysis with one more site year. The iterations continued until the analysis for the final number of site years (nine) was completed.
Results and Discussion
Temperature and precipitation
Monthly temperature averages and total monthly precipitation during the growing season are shown in Table 2 (Environment and Climate Change Canada 2020). September was included in the growing season due to the ability of annual canarygrass to keep growing, tillering and increasing grain yield if there is sufficient precipitation and soil moisture. In 2012, early season precipitation, Apr. and May, was above the 30-year average at all 5 locations. At three locations, Indian Head, Melfort and Scott precipitation for the rest of the growing season (Apr. to Sept.) stayed above the 30-year average while at Swift Current and Redvers, precipitation during the rest of the growing season was low enough to result in the total precipitation for the growing season being below the 30-year average. In 2013, precipitation was lower than normal for most locations over both the early and entire growing season, with Redvers having the lowest precipitation at only 67% of normal. Swift Current was the exception with 142% of the 30-year average for growing season precipitation. All test sites had higher than normal temperature in 2012, ranging from 102% of the 30-year normal in Indian Head and Melfort, to 111% of the normal at Swift Current. Temperature was below normal for all test sites in 2013, ranging from 90% of the 30-year normal at Indian Head to 99% at Swift Current.
Table 2.
Summary of climatic conditions at sites located near Indian Head, Melfort, Scott, Swift Current, and Redvers, Saskatchewan in 2012 and 2013.
Plant and panicle density
Plant density was not affected by N, fungicide application or the N × fungicide application interaction (Table 3). This indicated that increasing the N rate had no impact on germination or emergence of the annual canarygrass. Plant density ranged from 203 to 211 plants m−2 (Table 4). Since the fungicide was applied after the plant density was measured, it would not be expected to affect plant density. Panicle density responded to N rate but not to fungicide application or the interaction of fungicide × N (Table 3). As the N rate increased, there was a linear increase in the panicle density from 408 to 462 panicles m−2. Increasing the panicle density increases the yield potential of the annual canarygrass. Similar trends were reported by May et al. (2012a) where N fertilizer rate did not impact plant density while causing a curvi-linear increase in panicle density from 465 to 507 panicles m−2 as the N rate increased from 20 to 70 kg N·ha−1. In this study the increase in panicle density started to increase at a slightly higher N rate, 30 kg N−ha−1 compared with 20 kg N·ha−1 in the previous study and in this study it panicle density continued to increase to 90 kg N·ha−1 while in the previous study panicle density peaked at 70 kg N·ha−1. This difference may be due to higher panicle density obtained by the previous study at a low N rate 465 panicles m−2 compared with the 407 panicles m−2 at a low N rate in this study so less N was required in the previous study to maximize panicle density. The increase in panicle or head density as the N rate increased has been observed in other small grains under similar environmental conditions including oat, wheat and barley (May et al. 2020a, 2020b; O’Donovan et al. 2015).
Table 3.
Analysis of variance for the effect of fungicide and nitrogen on the development and grain yield of annual canarygrass over all sites and years.
Table 4.
The effect of fungicide and nitrogen (N) on the development and grain yield of annual canarygrass.
Kernel density
Kernel density was affected by N and fungicide application, but not by the N × fungicide interaction (Table 3). The kernel density increased linearly as the N rate increased from 14 240 kernels m−2 with 10 kg N·m−2 to 16 648 kernels m−2 with 90 kg N·m−2 (Table 4). Fungicide increased kernel density from 15 197 kernels m−2 to 16 288 kernels m−2. As a yield component, this increase in the kernel density increased the yield potential of the annual canarygrass.
Kernel weight
Fungicide application and the N rate × fungicide application interaction had no effect on kernel weight (Table 3). There was a small linear decrease of 2.6% in kernel weight as N rate increased from 7.59 g·1000 kernels−1 at 10 kg N·ha−1 to 7.4 g·1000 kernels−1 at 90 kg N·ha−1 (Tables 3 and 4). Kernel weight is the last yield component that can adjust to changes in environmental conditions during seed filling. This decrease in kernel weight as N rate increases indicates that the yield potential created by the increasing N rate became limited by other factors including environmental conditions. In contrast to this study, both May et al. (2012a) and Holt (1988) found that kernel weight was not affected by N rate. Holt (1988) had precipitation at or above the 30-year average during the years he conducted his study. May et al. (2012a) had precipitation that was closer to the current study with site years with precipitation at or above the 30-year average and other site years with a low level of precipitation during the growing season below the 30-year average. May et al. (2012a) did have a small numerical decrease in kernel weight as the N rate increased, but it was not statistically significant. It is important to note that the increase in yield potential as the N rate increased occurred in the yield components prior to kernel weight.
In this study the application of a fungicide affected kernel density but not kernel weight. Similar results in which kernel weight was the yield component affected by the fungicide application where found in oat (May et al. 2014a) and wheat (Kelley 2001); however, other studies found that both kernel density and kernel weight were increased by a fungicide application in other cereal crops including oat (Dietz et al. 2019), durum wheat [Triticum turgidum L. ssp. durum (Desf.) Husn] (May et al. 2014b), and barley (Safavi et al. 2012).
Height and lodging
Plant height was affected by N rate but not by an application of fungicide or the N rate × fungicide application interaction (Table 3). There was a linear increase in plant height from 91 cm to 95 cm as the N rate increased from 10 to 90 kg·ha−1 (Table 4). This small increase in plant height would not be expected to have a big impact on the development or grain yield of the annual canarygrass. May et al. (2012a) reached the same outcome, with N having a small effect on crop height, peaking at 101.7 cm with 70 kg N·ha−1. Holt (1988) reported a much bigger increase in plant height from 68 cm to 88 cm as the N rate increased from 0 to 100 kg·ha−1.
Lodging was not affected by a fungicide application. Holzapfel (2014) noted that on a field scale a fungicide application tended to reduced lodging improving the harvestability of the annual canarygrass crop. There was a small linear increase in lodging from 1.0 to 2.4 as the N rate increased from 10 to 90 kg·ha−1 (Table 4); however, this increase in lodging would not be expected to result in a yield decrease given that overall lodging especially at the higher N rates remained relatively minor. May et al. (2012a) also found that increasing the N rate resulted in a small increase in lodging, from 2.0 to 2.7 as the N rate increased from 20 to 100 kg N·ha−1 with little or any impact on grain yield. In this study, no significant interaction was found between fungicide and N for lodging. There was a significant site year × N interaction for lodging with two of the nine sites-years having a moderate increase in lodging from 1 to 5 as the N rate increased, and four of the nine site years having no change in lodging as the N rate increased (data not shown). Therefore, in this study, increasing the N rate occasionally induce moderate lodging at specific site-years but at most site years the impact of N rates on lodging in annual canarygrass was low.
Leaf disease
Leaf disease severity was not affected by N rate, fungicide application or the N rate × fungicide application interaction and was fairly low throughout the trial (Tables 3 and 4). Relatively low levels of disease, 2.2 to 3.6, occurred when averaged across the site years and this may be why the overall fungicide was not significant. However, there was a significant site year ×fungicide interaction with 32% of the variation due to the random effects for leaf disease, indicating that at some locations, fungicide application had an effect on disease severity. The no fungicide treatment had elevated levels of disease at 4 out of the 9 site years with the application of a fungicide reducing the disease rating by at least 1 on the rating scale. The low level of disease at the other 5 sites years resulted in the overall effect of fungicide on disease being non-significant. The elevated disease severity was recorded at two locations, Indian Head with McFadden ratings of 5.7 in 2012 and 8.2 in 2013 and Scott in with McFadden ratings of 7.8 in 2012 and 6.1 in 2013. These four site years do not have environmental conditions in common with growing season precipitation ranging from 141% to 73% of the 30 y average, but what they do have in common are that large commercial fields of annual canarygrass are grown near Indian Head and Scott. This could allow inoculum of septoria leaf mottle to enter the plots at Indian Head and Scott earlier in the growing season and be detected by the disease ratings. A lack of consistent effect of fungicide on disease severity ratings combined with a consistent beneficial impact on grain yield may be due to measuring disease severity too early in the growing season. It should be noted that from personal observation that the severity of septoria leaf mottle does not increase until canopy closer occurs at anthesis. In addition annual canarygrass will initiate tiller later in the growing season if precipitation occurs during anthesis and early seed filling resulting in later cohorts of tillers. In the current study, the only yield component that was increased by a fungicide application was kernel density, which suggests that the panicles could set and support more seeds during reproductive development. Unfortunately, we did not take measurements to determine if this increase was on all panicles or mostly from a specific cohort of tillers (first, second or third) and if a late cohort was filling kernels in late Aug. and Sept. It is interesting to note that increasing the N rate had no impact on the disease severity observed in this study and this lack of response was consistent among site years.
Days to maturity
Days to maturity was affected by N rate but not by fungicide or the N rate × fungicide application interaction (Table 3). As the N rate increased annual canarygrass maturity was delayed from 95.9 to 96.6 d (Table 4). Overall the days to maturity observed in the current study was slightly shorter than the 100–110 d reported by Putnam et al. (1996); however, it is representative of the earlier maturity found in cultivars developed at the University of Saskatchewan breeding program, which typically reach maturity in 98 d (Government of Saskatchewan 2020b). In contrast to the current study, May et al. (2012a) reported that maturity did not respond to different rates of N fertilizer.
Grain yield
Fungicide application and N rate both affected grain yield while the interaction of N rate × fungicide application had no affect on grain yield (Table 3). The application of a fungicide increased grain yield by 8.5% from 1139 to 1236 kg ha−1 (Table 4). The yield component that was responsible for increase in grain yield was kernel density. This indicates that the fungicide application either prevented seed abortion or increased initial seed set at anthesis. This suggests that at least a portion of this increase is due to improved control of septoria leaf mottle preventing a reduce supply of photoassimilates. The beneficial impact of a fungicide application on the grain yield of annual canarygrass was supported by Holzapfel (2014). In field scale trials near Indian Head, SK from 2004–2013, Holzapfel (2014) reported that, overall grain yield increased 23.2% following fungicide application versus the 8.5% increase observed in the current study.
As the N rate increased from 10 to 90 kg N·ha−1, there was a 14.4% linear increase in grain yield from 1087 kg·ha−1 to 1243 kg·ha−1 (Fig. 1). In this study, panicle density and kernel density were the yield components that responded to the applied N and contributed to the increase in grain yield. May et al. (2012a) found that panicle density was the yield component in annual canarygrass that responded to N fertilizer and contributed to grain yield. The lack of an interaction between N and fungicide application on grain yield indicates that growers can manage N fertilizer and disease control using foliar fungicides separately. This also indicates management of septoria using foliar fungicides in annual canarygrass does not impact the responsiveness of annual canarygrass to N fertilizer. In contrast to the current study, Holt (1988) reported a curvilinear increase of 79% in grain yield as the fertilizer N rate increased from 0 to 100 kg N·ha−1, whereby grain yield increased from 710 to 1270 kg·ha−1. May et al. (2012a) also reported a small curvilinear increase of 7% in grain yield by increasing N rate up to 78 kg·ha−1 with a maximum yield of 1216 kg·ha−1. Most of the increase occurred between the 20–40 kg N·ha−1 range, and using the coefficient of variation of each N rate across site years, they found that N fertilizer rates above 40 kg N·ha−1 had greater variation around the mean than the 20 and 40 kg N·ha−1 rates. The larger increase from applied N fertilizer observed by Holt (1988) may be due to the fact that, the research conducted by Holt was in a conventional tillage system verses a no-till cropping system in this study and in the May et al. (2012a) study. Long term no-till tends to have a higher level of residual N in the soil than a conventional cropping system (Jat et al. 2018; Lafond et al. 2011). In addition, the larger increase may be partially due to the fact that Holt (1988) started at a lower N rate, 0 kg N·ha−1, than this study with an initial N rate of 10 kg N·ha−1 and May et al. (2012a) with an initial N rate of 20 kg N·ha−1. Holt found that grain yield increased by 220 kg·ha−1 as the N rate increased from 0 to 25 and by 340 kg·ha−1 as the N rate increased from 25 to 100 kg N·ha−1.
Test weight
Test weight was not affected by fungicide application but had a significant linear decrease with increasing N application (Tables 3 and 4). The test weight decreased from 345.3 g 0.5 L·ha−1 at 10 kg N·ha−1 to 339.9 g 0.5 L−1 at 90 kg N·ha−1. This is an important factor as test weight may indicate the physical quality of the grain, which helps to determine the price at which a farmer can sell the grain. Holt (1988) also found that the test weight decreased by a small amount from 308 to 304 g 0.5 L−1 as the N rate increased from 25 to 100 kg·ha−1. It is interesting that in a genetically related cereal with a panicle type of inflorescence, oat, test weight declines as the fertilizer N rate is increased (Lafond et al. 2013; May et al. 2004; 2020b). Currently there is no clear explanation for this relationship.
Change in gross return
The change in gross return as the N rate increased was investigated using three N fertilizer prices and three grain prices. Under most of the nine combinations, the change in gross return as the N rate increased was not significant and the contrasts are included in Table 5 and Supplementary Table S1 (cjps-2020-0221suppla.docx)11. Since both the N rate and the contrasts were not significant no N rate can be recommended based on gross return; however, based on yield, N rates that provided the highest yields were at or above 50 kg·ha−1 N. It is interesting to observe that as the price of the annual canarygrass increased to 0.88 $·kg−1 and the cost of the N decreased to 0.50 $·kg−1, the change in gross return had its greatest range of 0 to 112 $·ha−1 and a linear contrast with a P value of 0.0685 (Table 5); however, as the cost of nitrogen was increase to 0.75 and 1.0 $·N kg−1 the range in gross return was numerically lower making it more difficult to detect differences as can be seen from the P values of the linear contrasts 0.1132 at 0.75 $·N kg−1 and 0.2730 at 1.0 $·kg−1. The variation in the data is too great to further refine the changes in gross return due to fertilizer and grain price.
Table 5.
The change in gross returns compared with 10 kg N ha−1 as the nitrogen (N) rate increased at three grain prices and three nitrogen prices over all sites and years.
Power of detection related to replication and location
When a study is conducted on a crop with limited research dollars, a better understanding of the number of replications and locations required to detect statistical differences is important when planning for future research. Statistical power is 1 – the type II error (not rejecting the null hypothesis when the null hypothesis is incorrect); therefore, as the statistical power increases the probability of a type II error decreases. Power analysis can be used to calculate the minimum sample size required to detect an effect of a given size. The approach used in this study and outlined by Stroup (1999) has been used in a wide range of research including Wireworms, [Limonius californicus (Mannerheim) (Coleoptera: Elateridae)], a balsam fir forest, dairy nutrition and fish community composition (Bilotta et al. 2016; D’Orangeville et al. 2013; Kononoff and Hanford 2006; Loiewski and Wenninger 2019). When the data from this study was used to estimate the number of replications within a location and the number of locations required to detect significant treatment differences in grain yield the results varied between fungicide and N rate (Table 6). In this study we wanted to consistently obtain a power of 0.80 (Di Stefano 2003) or more before recommending a set number of replicates or locations for future studies. The analysis randomly selected individual site years to include and since the severity of disease varies among locations, the power to detect differences fluctuates as we add site years that differ in disease severity (Table 6). To consistently detect the overall effect of a fungicide application on grain yield over the region in which the study was conducted, 4 site years are recommended (Table 6). Increasing the replicates at each location did not have as great an impact as the number of site years. To differentiate the effects of a fungicide application on grain yield at different locations (site year × fungicide) the number of replicates at the four locations would have to be increased from 4 to 6 or the number of locations would have to be increased from 4 to 5 to obtain a power greater than 0.80.
Table 6.
The change in the statistical power to detect significant differences as the number of replicates and site years change for the variables investigated in this study.
To consistently detect the impact of N fertilizer on grain yield over the region in which the study was conducted, the analysis indicates that 3 site years would be required. To differentiate between grain yield responses to N at different locations, 4 site years with 4 replicates or 3 site years with 8 replicates at each location would be required (Table 6). The need for fewer site years to detect the impact of N fertilizer than a fungicide application is probably due to the impact of N fertilizer varying on a range from small to large among the different site years while fungicides will have a large affect at the locations with elevated levels of disease and low to none at locations with low levels of disease. In this study we could not detect a interaction between N fertilizer and a fungicide using the entire data set (Table 3) so as expected an interaction was not detected using a subset of the entire data set (Table 6).
Conclusion
This study examines the relationship of N fertilizer and fungicides on crop development, septoria leaf mottle, grain yield and physical grain quality of annual canarygrass. Both N fertilizer and fungicide application increased grain yield. The response curve for N is slight with a rate of 50 kg N·ha−1 being sufficient to maximize yield and economic returns. In annual canarygrass it does not appear that the N rate recommended from this study would change if the cost of N was significantly increased. The results of this study indicate that growers cannot use leaf disease management with fungicides as a means to increase annual canarygrass’s responsiveness to N fertilizer. The lower impact of disease on grain yield observed in this study compared with field scale research (Holzapfel 2014) indicates that more research is required into the effects of septoria leaf mottle on annual canarygrass and how to best measure the impact of disease on annual canarygrass. When measuring the impact of a fungicide on grain yield, a power analysis indicates that at least 5 site years will be required for future research. Finally, both N fertilizer and fungicide applications are important management practices in annual canarygrass. The lack of an interaction between N and fungicide application indicates that these two practices can be managed independently by annual canarygrass producers. The results from this study indicate that producers should use 50 kg N·ha−1 at seeding and apply a fungicide to control septoria leaf mottle after full emergence of the flag leaf.
Acknowledgements
The technical support of Rebecca Davies, Randy Shiplack, Kevin Willoughby, Kathy Ringdal, and Orla Willoughby at Indian Head, Lana Shaw at South East Research Farm, Mike Hall at the East Central Research Farm and Darwin Leach at the Melfort Research Station is recognized and appreciated. Statistical advice provided by Dr. Stevenson is appreciated. Funding was provided by Saskatchewan Ministry of Agriculture, the Canaryseed Development Commission of Saskatchewan and Agriculture and Agri-Food Canada.