Soybean and lentils are major legume crops grown in Saskatchewan, with an average yield of 1398 kg·ha−1 and 1453 kg·ha−1, respectively . Soybean is well suited to more moist regions of the prairies such as the Dark Brown and Black soil zones in Saskatchewan  and is reported to be responsive to starter N in other jurisdictions . Lentil, which is better adapted to drier soil conditions, is mainly grown in southern and west-central portions of the province in the Brown and Dark Brown soil zones .
Growers use MAP and urea as major N and P sources in the northern Great Plains of North America, including western Canada . Further, adding fertilizer in the seed-row at planting/seeding is an effective strategy for fertilizer placement, especially to provide a starter source of nutrient for early crop nutrition and growth. Seed-row placement is identified as especially effective for immobile nutrients like P in cold, prairie soils. This is to ensure the supply of the nutrient early on to the roots of seedlings when it is needed for cell division early in the plant growth cycle .
Even though soybean and lentil can supply the majority of their N requirement through BNF via the Rhizobium symbiosis , it takes one to two weeks for nodules to become established. In a controlled environment (24˚C day, 18˚C night) study of nitrogenase enzyme activity of soybean (Glycine max var. Kent), nitrogenase activity was nil until 12 days of age of soybean and multiplication of bacteroids per vesicle only began after 20 days . Cooler soil conditions at the time of planting such as occur in northern prairie soils in early spring will further delay biological activities such as atmospheric nitrogen (N2) fixation. Starter N in the seed-row of a legume can help to carry the crop by supplying N until the nodules get established and begin fixation . Further, fertilization is deemed more important under stress conditions such as when root development is hindered by root disease, as well as when available N is limited by residue immobilization or any other environmental factors. Therefore, a small amount of starter N-P fertilizer may enhance the early season growth and nutrient uptake that could translate into a final grain yield benefit.
A concern surrounding seed-row placement of fertilizer with crops is that too much fertilizer in close proximity to seed can damage the seedling through desiccation from an osmotic effect. As well many N fertilizer sources can also cause significant damage through ammonia toxicity. Excessive fertilizer in the seed-row can to serious damage or death of seedlings and reduced stand count. Significant reductions in emergence of soybean in a Saskatchewan soil were observed when rates of P applied as 11-52-0 fertilizer exceeded 20 kg P2O5 ha−1 . Currently, guidelines are available for maximum safe rates of P applied as MAP and also combinations of MAP with potash for common pulse crops . However, there is no information on the tolerance of soybean and lentil to seed-placed fertilizer blend that contain equal proportions of N and P from a MAP and urea fertilizer blend. Therefore, a field study was conducted to assess the effect of varying rates of seed placed fertilizer blend: 50% Urea + 50% MAP applied at 0, 10, 20 and 30 kg N ha−1 on soybean and lentil emergence, grain and straw production, and nutrient uptake. Given the limited information on response of soybean to fertilization under northern prairie conditions, the influence of fertilization on the Ndfa via biological fixation was determined using an 15N label application to the soybean and use of non-fixing wheat as a reference crop.
2. Materials and Methods
2.1. Site Description
The field study was initiated on an annually cultivated field located at SE36-20-04-W3 (50.733107N; 106.424019W). This study site was located approximately 7 km south and 7 km east of the town of Central Butte, within the Rural Municipality of Enfield, No. 194, Saskatchewan, Canada. The site is located within the Brown soil zone but is very close (~2 km) to the boundary between the Brown and Dark Brown soil zone. The soils transition from loam to clay loam in texture and are classified as a Brown Chernozem . The soil map unit in which the site is situated describes the site as consisting predominantly of soils of the Haverhill association. Haverhill association can be described as Brown soils that formed in loamy glacial till. It occurs on undulating and undulating dissected landscapes with gentle slopes. Surface texture is commonly described as loam . Prior to the field trial in 2018, wheat (Triticum aestivum) was grown on the site in 2017, and canola in 2016.
2.2. Initial and Post-Harvest Soil Sampling
To characterize baseline soil properties at the field trial site, including nutrients before the experiment, pre-seeding composite samples of soil were taken from the 0 - 15, 15 - 30 and 30 - 60 cm depths on May 2018. Ten soil cores were taken in a transect across the site area with a Dutch auger. Soil samples were placed in a cooler for transport back to Saskatoon and then stored in a refrigerator at 5˚C until preparation for chemical analysis. Each soil sample collected was air-dried, sieved and the <2 mm fraction was retained and analyzed for various extractable nutrients (extractable P, potassium (K), nitrate ( )-N, sulfate ( )-S) and chemical properties (pH, electrical conductivity, % organic carbon) (Table 1).
Post-harvest soil sampling took place in September, 2018. A composite sample was obtained for each plot by collecting two cores within each plot (in-row and between-row) from all treatments and combining the sub-samples according to depth. A hydraulic punch truck fitted with a 5 cm diameter barrel was used to remove samples at three depth increments (0 - 15, 15 - 30 and 30 - 60 cm). Samples were placed in a refrigerator at 5˚C until preparation for chemical analysis, which consisted of air-drying samples at 30˚C followed by grinding with a Wiley mill to pass a 2 mm sieve. After processing, the samples were stored at room temperature until chemical analysis.
Table 1. Summary of baseline soil properties in soil cores collected from the field site locations in spring 2018. Values are means from analysis of ten individual soil cores collected in May across the field site before any treatments or field operations were conducted.
†P and K = Modified Kelowna extractable PO4-P and K (Qian et al., 1994). ‡N = CaCl2 extractable nitrate, NO3-N (Houba et al., 2000). §S = CaCl2 extractable sulphate, SO4-S (Houba et al., 2000). ¶pH measured in a 1:2 soil: water suspension (Hendershot et al., 2008). #EC measured in a 1:2 soil: water suspension (Miller and Curtin, 2008).
2.3. Experiment Design and Outline of Treatments
The experiment was set up as a randomized complete block design (RCBD) with four replicates of each treatment. Four rates of starter N: 0, 10, 20 and 30 kg N ha−1 were used in the form of a blend of 50% urea (46-0-0) and 50% MAP (11-52-0) with a blend analysis of 28-26-0. In addition to N, this supplied 0, 10, 20 and 30 kg P2O5 ha−1. The fertilizer products were obtained from a local agricultural fertilizer retailer (Nutrien Inc). Fertilizer was placed in the seed-row with the lentil and soybean under ~15% seed bed utilization (opener spread/row spacing) using a seeding tool designed and utilized for seeding small scale field plots. Modern varieties were selected of soybean (Glycine max L. cv. NSC Watson RR2Y), and small red lentil (Lens culinaris L. cv. CDC Maxim). Soybean was selected as the crop for the N fixation study as there is limited information on how N fixation in short season soybean responds to starter fertilizer in the prairies. Hard red spring wheat (Triticum aestivum L. cv Brandon) was seeded as the reference crop for soybean N uptake and BNF.
The area of each plot was 3.0 m × 1.0 m with 3 seed-rows per plot. Row length was 3.0 m and row spacing was 25.4 cm (10 inch spacing) as this is a common spacing configuration for these crops in Saskatchewan. At the first-trifoliate leaf stage of the soybean, a 1 m × 1 m subplot was established in each of the soybean and reference crop plots. A 10 atom % excess 15N-(15NH4) (15NO3) fertilizer was uniformly applied to each subplot at the rate of 5 kg N ha−1 in a liquid form dissolved in deionized water.
2.4. Field Operations
2.4.1. Site Preparation
The whole plot area was lightly rototilled and harrowed to distribute residue and smooth the seed bed. On May 9, 2018, prior to seeding, herbicide was applied for pre-plant control of weeds. Assigned plots for soybean and lentils were sprayed with the liquid herbicide product Vector® (540 g.a.i glyphosphate) at 0.8 l/ac and Aim® (carfentrazone) at 0.017 l/ac.
On the May 22, 2018, prior to seeding, soybean was inoculated with Bradyrhizobium japonicum in peat form (Nodulator® SCG peat-BASF). At the time of seeding, a granular form of this inoculant was also added to the seed row to provide a double inoculation to ensure adequate populations of the B. japonicum in the soil to enable infection and onset of N fixation in the soybean. The lentils were inoculated with Rhizobium leguminosarum (Nodulator® XL peat-BASF) peat-based inoculant. Seeding operations took place on May 23, 2018, with the date chosen to confirm that soil temperature at the depth of seeding was 10˚C . Prior to seeding, a packer was run across the plot area to firm the seed bed. The entire study area was seeded with a small plot disk seeder at 25.4 cm spacing, for a total of three rows within each 3 m × 1 m plot and at a rate of 79 kg soybean seeds·ha−1 and 56 kg lentil seed·ha−1 as per recommended rates . A packer wheeler, attached behind the disk seeder, helped ensure good seed to soil contact. Care was taken to ensure a suitable seeding depth of 2.5 cm, as recommended .
2.4.3. Weed Control
Application of herbicide after seeding was performed for the first time on June 15, 2018 for weed clean-up with a hand sprayer. Soybean and lentil plots were sprayed with imazamox at 11.7 g·ac−1 and wheat plots were sprayed with fluroxapyr at 0.125 L·ac−1 and 2,4-D at 0.340 L·ac−1.
A second application of herbicides was made on June 27, 2018. Imazamox was applied at 11.7 g·ac−1 rate with a hand sprayer between the rows of lentils and soybean plots to remove any persisting weeds from the first application, along with some hand hoeing. Fluoxapyr was applied between the rows of the wheat to control weeds.
Fungicide was applied on June 27, 2018 using a field sprayer. Soybean and lentil plots were sprayed with prothioconazole and trifloxystrobin as Delaro® at 0.350 L·ac−1.
2.4.4. Climate Data
Climate data (Table 2) was obtained from a weather station located about 500 meters from the field site.
2.4.5. 15N Application
The 15N application treatments were made on June 20, 2018 for quantifying N fixation using the 15N isotope dilution technique for soybean, and hard red spring wheat (Triticum aestivum L. cv Brandon) was used as the non-fixing reference crop . Prior to application of 15N, a stock solution was made by dissolving 58.8 g of 10 atom% excess 15N-(15NH4) (15NO3) fertilizer in 400 ml deionized water to provide the desired rate of application in the field of 5 kg N ha−1. When soybean reached the first trifoliate leaf, a 1 m × 1 m subplot on the north
Table 2. Comparison of mean monthly precipitation (mm) and temperature (˚C) during 2018 growing season at the field site to the 25-year (1992-2017) average.
HM = Historical mean (1992-2017) .
end of each plot was marked in each of the soybean plots and wheat reference crop plots for application of 15N. A plastic frame was placed on the ground to confine the fertilizer solution within the plot area. For field application of 15N, 10 ml of the stock solution was dissolved in 4 L of deionized water for application, followed by another 4 L of water application over the plot area to rinse off the residual fertilizer contained on the leaves into the soil.
2.4.6. Crop Harvesting and Processing
The crops were harvested at physiological maturity. The harvesting operation of lentil and non-15N labelled wheat was conducted on August 14, 2018 and soybean (15N applied) and 15N applied wheat harvesting operations were conducted on September 04, 2018. Two sub-samples of 1 m length from each plot were hand cut approximately 2.5 cm above ground. Crop samples from each plot were placed in tagged cloth bags and allowed to air dry for 1 - 2 weeks prior to threshing operations.
Soybean, lentil and wheat were threshed, and grain was cleaned in October 2018. Prior to threshing, whole biomass (grain + straw) harvested yields were measured and reported on a kg·ha−1 basis. Threshing, de-awning and cleaning was accomplished using a Wintersteiger® LD 350 thresher. After threshing, the grain was weighed. To obtain above ground straw yield, the grain weight was subtracted from the whole biomass (grain + straw) yields. A sub sample of straw and grain of each crop (soybean, lentil and wheat) was collected in an 8-dram vial. Finally, grain was finely ground using a NutriBullet Balance® grinder/mixer for determination of seed nutrient content by acid digestion.
For 15N analysis, samples (soybean and 15N wheat) were further oven-dried over night at 60˚C to a stable weight, and subsequently ground again to a fine powder using an 8000D Mixer/Mill® ball mill with dual clamps and stored in 8-dram vials for mass spectrometer analysis.
2.5. Laboratory Analyses
2.5.1. Acid Digestion for N and P Uptake
An acid digest of ground grain and straw was conducted . Briefly, 0.25 g (±0.001 g) of finely ground plant grain or straw was weighed into glass digestion tubes and 5 mL of concentrated (conc.) sulfuric acid (H2SO4) was added. Samples were placed on a digestion block at 360˚C for 30 min. Following this, samples were removed from the digestion block, allowed to cool, and 0.5 mL H2O2 was added. Samples were then placed on the digestion block an additional three times for 30 min, adding H2O2 after each heating period. Finally, samples were placed on the digestion block for 1 h. After samples were allowed to cool, distilled water was added to dilute the final volume of the sample to 75 mL to achieve a final concentration within the detection limit of the instrumentation. Samples were placed in a refrigerator until analysis for N and P by Technicon automated colorimetry.
2.5.2. Plant Available Soil and
Plant available soil and were extracted according to the methods described in . 20g of air-dried soil was weighed into an extraction bottle and 40 ml of pre-prepared 0.01 M CaCl2 solution (1.11 g CaCl2 into 1 L of distilled water) was added into it. The resulting solutions were shaken for 30 min at 142 RPM on a rotary shaker. After shaking, the suspension was filtered through Whatman No. 42 filter paper into 8-dram vials. Vials were refrigerated until analysis for soil and by automated colorimetry and plasma emission spectroscopy, respectively.
2.5.3. Plant Available Soil P and K
The Modified Kelowna extraction procedure  was used to extract plant available orthophosphate and potassium (K+). In this method, 3 g of air-dried soil was weighed into an extraction bottle and 30 ml of pre-prepared Kelowna solution (0.25 M acetic acid + 0.015 M NH4F + 0.25 M NH4OAc) was added into the bottle. Samples were then shaken horizontally on a rotary shaker for 5 min at 142 RPM. After shaking, samples were filtered through VWR 454 filter paper into 8 1/2 vials, and the filtrate was stored at 5˚C until analysis by AA-FE spectrometry for K and automated colorimetry for P.
2.5.4. 15N Analysis
The finely ground soybean and wheat grain and straw sub-samples were weighed (grains ~2.0 mg and straw ~4.0 mg) using a Sartorius Microbalance® micro-balance and encapsulated using 6 × 4 mm tin capsules (Elemental Microanalysis®) into an approximately spherical shape with the air pressed out of the encapsulated sample. The encapsulated grain and straw samples were then analyzed for percent N and atom % 15N using a Costech® ECS4010 elemental analyzer coupled to a Delta V Advantage® Mass Spectrometer, with a standard for the spectrometry measurement.
The nitrogen isotope dilution method provides direct evidence for N2 fixation since the 15N concentration in the plants exposed to 15N2 is greater than the 0.366% natural abundance but is diluted relative to the non-fixing reference crop if BNF occurs . The isotope dilution method introduced by McAuliffe  was used to determine N2 fixed (kg·ha−1) by calculating percentage of N derived from atmosphere (%Ndfa) and multiplying by total plant N calculated for the grain and straw separately:
Further, Hardarson and Danso  described that in the absence of any supply of N other than soil and 15N-labelled fertilizer, both plants will contain the same ratio of 15N/14N, since they are taking N of similar 15N/14N composition, but not necessarily the same total quantity of N. Therefore, the total amount of N derived from atmosphere (Ndfa, kg N ha−1) in plant biomass was calculated by % 15N atom excess (Ndff) of fixing and non-fixing crop as:
Above equation can be further classified as,
The total N in the fixing crop can be obtained by multiplying the dry matter by the % N as;
2.7. Statistical Analysis
Statistical analyses were conducted using PROC GLIMMIX of SAS, Version 9.4 (SAS Institute Inc, 2017). This study was conducted as an RCBD with a total of four treatments with four replicates. An ANOVA was conducted, with treatments as a fixed effect and block as a random effect. Where applicable, mean comparisons were performed using least significant differences (LSD; equivalent to Fisher’s protected LSD) at a significance level of 0.05. PROC UNIVARIATE was used to determine if the measured values of a treatment were normally distributed.
The mean soybean emergence measured at 19 days after seeding (Table 3) shows that at the 20 kg N, 20 kg P2O5 ha−1 seed placed fertilizer rate treatment, the emergence was significantly reduced compared to the 10 kg·ha−1 and control (no fertilizer treatment) (P = 0.0270). Similarly, above 10 kg·ha−1, lentil emergence was significantly lower (P = 0.0272).
Table 3. Mean soybean and lentil plant emergence counts (number of plants emerged in 1 m row length) at 19 days after seeding in 2018 at Central Butte site.
†Reported values are the means of four replicates (n = 4) of each treatment. Within a column, means followed by the same letter are not significantly different at P > 0.05 according to Fisher's protected LSD method. SEM = standard error of mean. CV = Coefficient of variation.
3.2. Grain and Straw Yield
Mean soybean and lentil grain and straw yield for the starter N, P fertilizer applied at different rates is shown in Figure 1 and Figure 2, respectively. For soybean yield (Figure 1), there was a trend for fertilization with starter fertilizer to result in higher soybean yield compared to the unfertilized control, with increasing rate resulting in slight decrease in mean grain yield. However, the effect was not statistically significant (P = 0.3994) among the treatments.
For lentil (Figure 2), mean grain and straw yield followed a similar trend to soybean, with slightly higher mean yield at the 10 and 20 kg N, P2O5 ha−1 rates, and showed decrease with the treatments from 10 to 30 kg·ha−1. However, even more so than for the soybeans, the yield differences among treatments were small and not statistically significant (for lentil grain and straw, P = 0.5339 and P = 0.1458, respectively. Overall response to fertilization appeared less in lentil compared to soybean.
3.3. Grain and Straw Nutrient Uptake
Mean soybean and lentil N and P uptake in grain and straw is shown in Table 4. In soybean, the N uptake in grain and straw was not significantly increased by N and P fertilization (P = 0.5227). There was a trend for about 10% - 20% higher mean N uptake in soybean yield components with fertilization, compared to control (no fertilizer). Increasing fertilizer rate had little evident impact on soybean grain or straw yield. For P uptake in soybean grain and straw, mean P uptake in soybean grain for all the rates (10, 20, 30 kg N, P2O5 ha−1) was significantly higher (P = 0.0425) than the unfertilized control, but among the fertilizer rates there were no differences. Straw N and P uptake by soybean showed a similar pattern to grain N and P uptake, but straw N uptake of fertilized treatments compared to unfertilized control have increased N uptake by 15% - 20%.
In lentils, both N and P uptake were quite similar among treatments (Table 4) and there was no significant difference in N and P uptake of grain in the different treatments, but a trend to higher N and P uptake in the grain. Straw N uptake was significantly higher in fertilized lentil straw. Compared to N and P uptake increases in soybean, the increases in uptake for lentil with fertilization
Figure 1. Mean soybean grain (yellow) and straw (green) yield measured in fall 2018 at Central Butte field study site. Fertilizer is blend of 50% urea and 50% MAP by weight (analysis: 28-26-0). Reported values are the means of four replicates (n = 4) of each treatment. For the same yield component (grain or straw yield), means followed by the same letter are not significantly different at P > 0.05. The multi-treatment comparisons were made using Fisher's protected LSD method.
Figure 2. Mean lentil grain (yellow) and straw (green) yield measured in fall 2018. Reported values are the means of four replicates (n = 4) of each treatment. Fertilizer is blend of 50% urea and 50% MAP by weight (analysis: 28-26-0). For the same yield component (grain or straw yield), means followed by the same letter are not significantly different at P > 0.05. The multi-treatment comparisons using Fisher’s protected LSD method.
Table 4. Mean soybean and lentil N and P uptake (kg·ha−1 uptake in grain and straw) measured in 2018 at Central Butte field study site.
† Reported values are the means of four replicates (n = 4) of each treatment. Within a column, means followed by the same letter are not significantly different at P > 0.05 according to Fisher's protected LSD method. SEM = standard error of mean.
versus without were smaller, only about 4% - 7% increase. Overall, soybeans were more responsive to starter N-P fertilizer in both increased yield and nutrient uptake compared to lentil, especially for P.
3.4. Nitrogen Fixation
The mean % ndfa in soybean grain and straw is provided in Table 5. The % ndfa in soybean grain was significantly lower with higher rates of the N, P2O5 urea-MAP fertilizer blend added. A reduction in % ndfa in grain and straw of about 10% can be observed at the 30 kg N, P2O5 ha−1 rate compared to the unfertilized control that can be mainly attributed to the effect of the added N. This is expected as N fertilizer addition at higher rates is known to significantly decrease the proportion of N in the legume plant derived from biological fixation of atmospheric N.
4.1. Emergence and Yield
In this study, the blend of 50% urea and 50% MAP giving a product blend analysis of 28-26-0, supplied an available N to P2O5 ratio of 1:1. The mean soybean and lentil plant emergence counts at 19 days after seeding (Table 3) showed the emergence was reduced by 20% - 25% at 20 and 30 kg N, P2O5 ha−1 compared to the 10 kg·ha−1 and control (no fertilizer) treatment. In a previous controlled environment study by Weiseth , MAP was placed in a seed-row with soybean, and only rates above 20 kg P2O5 ha−1 (40, 60 and 80 kg P2O5 ha−1) showed a significant reduction in emergence using a soil similar to that of the current study. Greater apparent injury per unit of P2O5 observed in the current study compared
Table 5. Mean N derived from atmosphere through biological fixation in soybean grain and straw measured in 2018 field trial at Central Butte.
† Reported values are the means of four replicates (n = 4) of each treatment. Within a column, means followed by the same letter are not significantly different at P ≤ 0.05. The multi-treatment comparisons using Fisher's protected LSD method. SEM = standard error of mean. CV = Coefficient of variation. ‡ % Ndfa denotes percentage of N derived from the atmosphere.
to Weiseth  is explained by use of a urea-MAP blend rather than MAP alone and also drier conditions experienced in the field, especially in May and June (Table 2), compared to optimum moisture conditions in a growth chamber. “Fertilizer burn” is related to salt effect of the fertilizer and ammonia/ammonium toxicity . Free NH3 is toxic and can move freely through the cell wall, whereas NH4+ cannot, which explains why N fertilizers that produce large amounts of free ammonia in the vicinity of the granule, like urea, contribute more to germination and seedling damage. Therefore, obviously a blend with urea will cause more damage than MAP alone due to hydrolysis of urea into NH3 and CO2. It is therefore suggested for both soybean and lentil, that as the blend approaches 50% urea and 50% MAP, the starter N, P2O5 rate is reduced to 10 kg N, P2O5 ha−1.
Low soil moisture and temperature can also adversely impact emergence along with too much fertilizer placed in the seed-row as discussed above. The current recommended minimum average soil temperature at seeding according to Saskatchewan Pulse Growers  is 5˚C and 10˚C for lentils and soybean, respectively. Mean air temperature at site at the time of seeding was ~15˚C and this was above the historical mean for the month but the average precipitation was below the average throughout the growing season (from May to August), about 20% and 4% lower in May and August, respectively (Table 2). This low moisture condition in spring seeding may also have aggravated the salt and ammonia injury potential of the seed-row fertilizer, particularly for soybean. Low precipitation along with high sand content and low organic matter content in the site will lead to greater concentration of the fertilizer salt in the soil solution, reduced ammonia absorption by water and colloids, and thereby lead to more injury at high N rates for both soybean and lentil. This result agrees with the controlled environment study done using soil from same field.
In 2017, the average grain yield per unit area for soybean in the area the study was conducted was reported as 1398 kg·ha−1 . The mean soybean grain yield obtained in this field trial in 2018 is higher than the average, with yield of 2003, 1926, and 1862 kg·ha−1 for the rates of 10, 20, and 30 kg N, P2O5 ha−1, respectively. The trend of slight yield reduction with higher rates of seed-row urea-MAP blend may be explained from the emergence reduction and reduced plant stand with the higher rates. The lack of large positive yield responses to fertilization reflects the good ability of the legumes, especially lentil, to access existing P reserves in the soil.
In a meta-analysis of soybean grain yield responses to N fertilization , there was a positive response to fertilizer N in about half of the published studies, with the average yield increase from N fertilizer addition in the N-fertilizer responsive studies of 520 kg·ha−1 (n = 154). The magnitude of the response did not significantly differ among N rate categories of 0 - 50, 50 - 100 and >100 kg N ha−1. The results of this 2018 field study indicate that 10 kg N, P2O5 ha−1 as starter seed-row placed urea-MAP blend would be optimal for the soybean in terms of maximizing plant stand and yield on this soil under the conditions of the study. This is supported by the results of Gai  that used urea as a starter fertilizer at four different N rates (0, 25, 50 and 75 kg N ha−1), with an average yield soybean yield around 3100 kg·ha−1 for all the treatments but no significant difference among the treatments.
As for soybean, the mean lentil grain and straw yield was slightly higher in 10 and 20 kg N, P2O5 ha−1 treatments than control, and yield showed a decrease from 10 to 30 kg N, P2O5 ha−1, but there was no significant difference among treatments. The mean lentil grain yield in this study: 2610, 2523, and 2334 kg·ha−1 at the rate of 10, 20 and 30 kg N, P2O5 ha−1, respectively, was higher than the provincial average yield. Slightly lower mean yield with seed-row rates of N, P2O5 above 10 kg·ha−1 for lentils observed in this study may be attributed to fertilizer injury reducing emergence and stand count. Therefore, similar to soybean, 10 kg N, P2O5 ha−1 as seed row placed starter 28-26-0 fertilizer for lentil can be considered best agronomic performer.
4.2. Crop Nutrition and Biological Nitrogen Fixation
In soybean, grain N uptake was not significantly increased by fertilization, but there was a trend for increased N uptake by about 10% - 20% with N, P fertilization, compared to the control (no fertilizer). The P uptake in soybean grain at all the rates of urea-MAP blend fertilization (10, 20, 30 kg N, P2O5 ha−1) was significantly higher (P = 0.0334) than the unfertilized control but among the fertilizer rates the effect was similar. The increase in P uptake can be explained by the P in the fertilizer contributing to soil available P supplies and additional P uptake by roots that was eventually translocated to the grain. For N, increases in soil available N from fertilization may be countered by reduced biological nitrogen fixation. The trend in increased N uptake by soybean with starter N, P2O5 fertilization may reflect additional uptake of soil and fertilizer N without equivalent reduction in N derived from biological nitrogen fixation, as well as enhancement of N uptake associated with P fertilization. There is typically a correlation of nutrient uptake with final crop yield for most plant nutrients . This relationship may exist even for legumes that can fix N biologically if it is not present in the soil. The linear relationship between soybean grain yield and total N uptake with above-ground biomass had a slope of 12.7 kg grain per kg N . Further, it was mentioned that on average, a soybean crop yielding 5000 kg·ha−1 accumulates about 400 kg N ha−1 in its aboveground biomass, and that N must be provided from indigenous soil resources, the biological fixation process, and/or fertilizer.
In lentils, both N and P uptake in yield components were similar among treatments and there was no significant difference in N and P uptake in grain and straw in the different treatments. Compared to lentil, soybeans were more responsive to starter N-P fertilizer in increased nutrient uptake, especially for P. Uptake amounts for lentil similar to that obtained in the current study were observed in a study conducted in the Moist Dark Brown soil zone at Scott, SK in 2014 with grain N and P uptakes of 90.1 kg N ha−1 and 7.8 kg P ha−1, respectively . Corresponding to the lack of large yield response of the soybean and lentil to starter N and P fertilization observed at the field site in this thesis, there were also no visual symptoms of N and P deficiency in any of the soybean and lentil treatment plots throughout the growing season. In this field study, for both soybean and lentil, highest grain N and P uptake was recorded at 10 kg N, P2O5 ha−1, which had highest grain yield. This is also the fertilizer rate that produced minimal reduction in emergence.
It has been noted by other researchers that little starter N is typically needed for good pulse yield because they can fix their own N from the atmosphere in the symbiotic relationship with Rhizobium once they establish nodules, and too much available N can reduce the contribution from symbiotic fixation . Concomitant with this concept, the % ndfa in soybean grain (Table 5) was lower with higher rates of N fertilizer, with about 10% reduction in the 30 kg N, P2O5 ha−1 rate compared to the unfertilized control. A similar effect was observed for % ndfa in straw but was not significantly different among the treatments. In this work, the highest soybean grain % ndfa with no fertilizer was 81%, with 79%, 72%, and 71% ndfa for 10, 20, and 30 kg N ha−1, respectively. The % ndfa values reported for soybean in the current study tend to be slightly higher than other studies. However, nodulation of soybean on 0 N treatment in the field was visually noted to be good. The % ndfa reported for soybean was 60% in , 57% in  and 50% in  and comparatively lower than the results obtained in this work. Nevertheless, a similar 15N dilution technique in a Black Chernozmem soil near Rosthern, SK has shown similar results with 74% ndfa for N in soybean grain .
A negative exponential relationship was typically observed between N fertilizer rate and N2 fixation when N was applied in the top 0 - 20 cm of soil or on the soil surface . The % ndfa results gained from the field study in this work supports this, as highest % ndfa was in the unfertilized treatment and lowest in 30 kg N, P2O5 ha−1. Therefore, the 10 kg N, P2O5 ha−1 rate of 28-26-0, which gave a % ndfa that was not significantly different from the control and produced the highest mean yield and N uptake, may be the best treatment for maximizing yield and contribution of biological fixation from soybean.
The field study conducted with soybean and lentil in 2018 in south-central Saskatchewan showed a rate of 10 kg N, P2O5 ha−1 as seed-row placed urea-MAP blend (28-26-0) appeared to be rate that did not significantly reduce emergence, stand count or biological nitrogen fixation and was sufficient to maximize yield, N and P uptake. The 2018 field site moisture conditions were low compared to the historical mean and overall conditions were very dry. Further, the current field study was conducted over only one season for lentil and soybean. Therefore, expanding field sites and soils would be recommended to encompass and match the environments where the different pulse crops are grown.
This work was supported by the Agriculture Development Fund-Government of Saskatchewan and Saskatchewan Pulse Growers. Further, the authors would like to thank Cory Fatteicher, Ryan Hangs, Gravel Wang, Ranjan Kar for their assistance during the experiments.
 Gai, Z., Zhang, J. and Li., C. (2017) Effects of Starter Nitrogen Fertilizer on Soybean Root Activity, Leaf Photosynthesis and Grain Yield. PLoS ONE, 12, e0174841.
 Schoenau, J.J., Qian, P. and King, T. (2007) Strategies for Improving the Efficiency and Crop Safety of Starter Fertilizer Phosphorus and Potassium. Agriculture Development Fund Project 20050725 Final Report.
 Environment and Climate Change Canada (2018) Historical Climate Data. http://climate.weather.gc.ca/historical_data/search_historic_data_e.html
 Hardarson, G. and Danso, S.K.A. (1990) Use of 15N Methodology to Assess Biological Nitrogen Fixation. In: Hardanson, G., Ed., Use of Nuclear Techniques in Studies of Soil-Plant Relationships, International Atomic Energy Agency, Vienna, 129-160.
 Thomas, R.L., Sheard, R.W. and Moyer, J.R. (1967) Comparison of Conventional and Automated Procedures for Nitrogen, Phosphorus, and Potassium Analysis of Plant Material Using a Single Digestion. Agronomy Journal, 59, 240-243.
 Houba, V.J.G., Temminghoff, E.J.M., Gaikhorst, G.A. and Van Vark, W. (2000) Soil Analysis Procedures Using 0.01 M Calcium Chloride as Extraction Reagent. Communications in Soil Science and Plant Analysis, 31, 1299-1396.
 Qian, P., Schoenau, J.J. and Karamanos, R.E. (1994) Simultaneous Extraction of Available Phosphorus and Potassium with a New Soil Test: A Modification of Kelowna Extraction. Communications in Soil Science and Plant Analysis, 25, 627-635.
 McAuliffe, C., Chamblee, D.S., Uribe-Arango, H. and Woodhouse, W.W. (1958) Influence of Inorganic Nitrogen on Nitrogen Fixation by Legumes as Revealed by 15N. Agronomy Journal, 50, 334-337.
 Salvagiotti, F., Cassman, K.G., Specht, J.E., Walters, D.T., Weiss, A. and Dobermann, A. (2008) Nitrogen Uptake, Fixation and Response to fertilizer N in Soybeans: A Review. Field Crops Research, 108, 1-13.
 Xie, J., Schoenau, J. and Warkentin, T.D. (2018) Yield and Uptake of Nitrogen and Phosphorus in Soybean, Pea, and Lentil and Effects on Soil Nutrient Supply and Crop Yield in the Succeeding Year in Saskatchewan. Canadian Journal of Plant Science, 98, 5-16.
 Yang, J.Y., Drury, C.F., Yang, X.M., Jong, R. De, Huffman, E.C., Campbell, C.A. and Kirkwood, V. (2010) Agriculture, Ecosystems and Environment Estimating Biological N2 Fixation in Canadian Agricultural Land Using Legume Yields. Agriculture, Ecosystems & Environment, 137, 192-201.
 Peoples, M.B., Brockwell, J., Herridge, D.F., Rochester, I.J., Alves, B.J.R., Urquiaga, S., Boddey, R.M., Dakora, F.D., Bhattarai, S., Maskey, S.L., Sampet, C., Rerkasem, B., Khan, D.F., Hauggaard-Nielsen, H. and Jensen, E.S. (2009) The Contributions of Nitrogen-Fixing Crop Legumes to the Productivity of Agricultural Systems. Symbiosis, 48, 1-17.