Effect of Mineral Nitrogen, Phosphorus, and Potassium Fertilizers on the Productivity of Faba Bean ( Vicia faba L.) in Acidic Soils of Wolaita Zone, Southern Ethiopia

Poor soil fertility is among the major factors that limit faba bean production in Wolaita Zone in southern Ethiopia. Therefore, a ﬁeld experiment was conducted in the Kokate Marachare subdistrict of Sodo Zuria District of the zone during the 2019 and 2020 cropping seasons to determine the response of faba bean to diﬀerent rates of nitrogen (N), phosphorus (P), and potassium (K) fertilizers under lime-treated soil conditions by using Tumsa faba bean variety. The treatments consisted of three rates of nitrogen fertilizer (0, 23, and 46kg Nha − 1 ), three rates of phosphorus fertilizer (0, 46, and 92kg P 2 O 5 ha − 1 ), and three rates of potassium fertilizer (0, 30, and 60kg K 2 O ha − 1 ) that were laid out as RCBD and replicated three times per treatment . The results indicated that the N, P, and K fertilizer combination at 23:92:60kgha − 1 increased plant height, the number of branches per plant, and stem girth by 18%, 62.6%, and 55.6%, respectively, compared with the control treatment. A signiﬁcantly high aboveground dry biomass (11.8tha − 1 ), number of pods per plant (17 pods), number of seeds per pod (4 pods), stover yield (6.83t ha − 1 ), and hundred seed weight (88g) were obtained from the N, P, and K fertilizer combination at 23:92:60kg ha –1 . The highest grain yield (4.97 t ha − 1 ) was obtained from the N, P, and K combination at 23:92:60kg ha − 1 , which was 360% higher than the yield obtained from the control treatment. Moreover, the highest mean net beneﬁt (USD 4,109.33ha − 1 ) with an acceptable marginal rate of return of 1,340% was obtained from the N, P, and K fertilizer combination at 23 kgN, 92kg P 2 O 5 , and 60kg K 2 Oha –1 , respectively. Thus, these rates are suggested for faba bean production in the acidic soils of the Wolaita Zone.


Introduction
Faba bean (Vicia faba L.) is one of the major pulse crops grown in the highlands of Ethiopia [1]. Currently, it occupies 31% of the area cultivated for pulses (1,863,445 ha) in the country [2]. e crop plays a significant role in human and livestock feed and in improving soil fertility [3]. However, the productivity of the crop in the country is low (2.12 t ha −1 ) compared with the average yield (3.7 t ha −1 ) obtained in major faba-bean-producing countries of the world [2,4]. Particularly, in Wolaita Zone, faba bean occupies 3.6% of the area of land cultivated for total grain production [5].
However, farmers in the study area harvest a lower average yield (1.2 t ha −1 ) than the national average yield obtained [5].
Poor soil fertility and soil acidity are serious problems constraining faba bean productivity in Ethiopia [6]. In most cases, soils with a soil pH value of less than 5.5 are deficient in macronutrients such as nitrogen (N), phosphorus (P), potassium (K), calcium (Ca), and magnesium (Mg) [7]. In addition, the toxicity of aluminum (Al), iron (Fe), and manganese (Mn) is among the reasons for lower faba bean yield in acidic soils [8].
ese reasons could affect the availability of applied fertilizer and crop nutrient uptake. Shanka et al. [9] and Mesfin et al. [10] indicated that soil acidity and related soil fertility problems are major crop yield limiting factors in the Wolaita Zone. Fanuel et al. [11] also reported that the deficiency of total N, available P, and K constrains crop production in Wolaita Zone. e deficiency of these nutrients requires immediate attention for sustaining production and improving the productivity of the smallholder faba bean farmers in the study area.
Nitrogen is the major essential nutrient applied to the crop for higher vegetative growth, productivity, and quality [12]. In addition, the application of starter N fertilization at the rate of 20 kg N ha −1 enhances the nodulation process in faba bean plants [13]. However, reports in Ethiopia have shown that about 50% of the applied N fertilizer remains unavailable to the crop. is is due to temporary immobilization in soil organic matter or due to losses by leaching, erosion, nitrification, or volatilization [14]. In this regard, the studied districts have been exposed to erosion, soil acidity, and a high loss of N due to leaching [10]. Confirming this problem, [8] reported that N deficiency is a major factor constraining crop production in Wolaita Zone. Shanka et al. [9] also reported that P deficiency in the soils is a constraint to common bean production in Wolaita Zone. Laekemariam et al. [8] indicated that potassium deficiency is also a constraint to crop production in Wolaita Zone. us, the survival, nutrient use efficiency, growth, and grain yield of faba bean are highly affected in Lolita Zone due to a lack of balanced nutrient application [15]. In general, continuous cultivation, together with inappropriate soil management, such as limited soil conservation, inadequate use of bio-fertilizers, and lack of balanced inorganic fertilizer application, has aggravated soil fertility problems in Wolaita Zone [10].
About 8.3% more yield response of faba bean was obtained from the combined application of N, P, and K at 36 kg N, 54 kg P 2 O 5 , and 114 K 2 O kg ha −1 compared with P alone applied at 54 kg P 2 O 5 ha −1 [16]. Abebe et al., 2014 [17] indicated significant contributions of N, P, and K fertilizer application to improved growth and grain yield productivity of faba bean. In many parts of Ethiopia, a significant improvement in the grain yield of faba bean was obtained in response to applying inorganic fertilizers [18]. For instance, Bezabih et al. [19] showed that the combined application of 30 kg P 2 O 5 ha −1 and 60 kg K 2 O ha −1 in rhizobium-inoculated faba bean contributed to 30% more grain yield than the control treatment in Alicho Wuriro Highlands in southern Ethiopia. In agreement with this, Abebe et al. [20] reported that the combined application of 46 kg P 2 O 5 ha −1 and 24 kg K 2 O ha −1 contributed a 38% more grain yield improvement in faba bean than the control treatment in Sekela District in northern Ethiopia. Similarly, a 31% mean improvement in the grain yield of faba bean was reported due to the application of phosphorus fertilizer at 15, 30, 45, and 60 kg ha −1 compared with the control treatment in Boloso Sore District in Wolaita Zone in southern Ethiopia [15].
To improve the grain yield of faba bean, a balanced application of all major nutrients (N, P, and K) is vital [21]. However, the optimum requirements of combined N, P, and K fertilizers for faba bean production are not well known in Wolaita Zone in general and the study area in particular. erefore, this research was conducted to determine the effect of the application of N, P, and K fertilizers on the productivity of faba bean in the Sodo Zuria District in Wolaita Zone. It was hypothesized that the application of balanced N, P, and K fertilizers improves the growth and yield of faba bean.

Description of the Experimental Site.
e study was conducted for two consecutive years during the 2019 and 2020 main cropping seasons (June to November) on farmers' fields in the Kokate Marachere subdistrict of Sodo Zuria District in Wolaita Zone, southern Ethiopia ( Figure 1). e study site was selected purposively based on its high faba bean production potential. e experimental site is located at 7°25′21″ N latitude and 37°46′52″ E longitude and an altitude of 2,156 m above sea level [9]. In the 2019 and 2020 cropping seasons, the mean monthly temperature of Sodo Zuria District ranged from 15°C to 23.8°C and from 14.5°C to 25.8°C, respectively (Table 1; Figure 2). In addition, in the 2019 and 2020 cropping seasons, the total monthly precipitation of the study district was 1,187 mm and 1,376 mm, respectively [22]. e dominant soil type of the study area is silty clay, which is deep and clay in texture [9].

Planting Material.
A faba bean variety "Tumsa," which was released in 2010 by the Holetta Agricultural Research Center in Ethiopia, was used. It needs 700-1,000 mm of rainfall for high yield and grows at an altitude of 900-2,800 m above sea level. e variety needs 120-130 days to reach maturity [23].

Treatments and Experimental Design.
e treatments consisted of three rates of N (0, 23, and 46 kg N ha −1 ), P (0, 46, and 92 kg P 2 O 5 ha −1 ), and K (0, 30, and 60 kg K 2 O ha −1 ). e middle rates of N and P were fixed based on the Ethiopian Institute of Agricultural Research's recommendation for faba bean production [22]. Although there is no recommended K rate for faba bean by the Ethiopian Institute of Agricultural Research, the rate of K was fixed based on research work in Ethiopia [13]. e experiment was laid out as a randomized complete block design in factorial arrangement and replicated three times per treatment. e randomization in the 2019 cropping was used in 2020 without any change in location in order to avoid biases. e plot size was 2.8 m × 2.1 m (5.8 m 2 ) with 0.5 m spacing between plots and 1 m between blocks. e inter-and intrarow spacings were 40 and 10 cm, respectively, and each row and plot consisted of 21 and 168 plants, respectively [23]. 2 International Journal of Agronomy 2.4. Soil Sampling and Analysis. Before planting, surface soil samples (0-30 cm depth) were collected in a zigzag pattern from 16 spots from the experimental field by using augur. en the soil samples were mixed to form a 1 kg composite sample. Afterward, the soil samples were air-dried, ground, and sieved through a 2 mm sieve to determine their pH and available phosphorus. e soil samples were analyzed at the Hawassa Soil Laboratory for soil pH, CEC, soil organic carbon (OC), total N, available P, exchangeable cations, and soil texture. Soil pH (1:2.5 soil-to-water ratio) was measured using a glass electrode pH meter, as described by [24]. e CEC of the soil was determined from NH 4 OAcsaturated samples, which was measured through distillation using the micro-Kjeldahl procedure. e soil OC was determined by the chromate acid oxidation method [25]. Total nitrogen was analyzed using the macro-Kjeldahl digestion method, followed by the ammonium distillation and titration method [26]. Soil available P was analyzed using the Olsen method [27]. Exchangeable K was extracted by the ammonium acetate (1 M NH 4 OAc at pH 7) extraction method, as described by Rowell [28], and determined by flame photometry. e particle size distribution was determined following the Bouyoucos hydrometer method [29], and the textural class was determined based on the soil textural triangle using the International Soil Science Society (ISSS) system [28].  e experimental field was ploughed three times with local oxendrawn maresha, followed by manual seed-bed preparation, and laid out according to the experimental scheme. Two seeds per hill were sown by hand at a seeding rate of 200 kg ha −1 with an intra-and inter-row spacing of 10 cm and 40 cm, respectively [23], on 8 June 2019 and 9 June 2020. Later, the seedlings were thinned to have one plant per hill. Weeds were removed by hand, weeding two times at 25 and 50 days after crop emergence. Harvesting was performed manually using hand sickles. e crop was harvested on 4 November 2019 and 17 November 2020 in the two cropping seasons.

Fertilizer Application.
Phosphorus fertilizer was applied in the form of TSP (46% P 2 O 5 ) at the time of planting at the specified rate. Nitrogen and K were applied at the specified rate in the form of urea (46% N) and KCl (62% K 2 O), respectively, by splitting into 1/3 at planting and the remaining 2/3 at the active vegetative stage before flowering [30]. e liming rate was fixed based on greenhouse findings as 2 t ha −1 and applied in broadcast 2 months before planting [9]. e Tumsa faba bean variety best performing in a pot experiment was used for the study.

Crop Data Collection
2.6.1. Growth Data. All growth data were recorded from five randomly selected plants from the central rows at the flowering stage. Plant height was measured from the base to the tip of the plant by using a measuring tape. e stem girth was recorded by using a digital caliper that was attached to the center of a stem. Lodging was assessed based on a scale of 1-5, where 1 (0-15°) indicates no lodging, 2 (15-30°) indicates 25% lodging, 3 (30-45°) indicates 50% lodging, 4 (45-60°) indicates 75% lodging, and 5 (60-100°) indicates 100% lodging [31]. e scale was determined by the angle of inclination of the main stem from the vertical line to the base of the stem by visual observation [32]. e score was then calculated as described by [33]. e number of branches was determined by counting all branches originating from the main stem, and secondary growth branches were excluded.

Yield and Yield Component Data.
e total aboveground sundry biomass was determined by weighing the shoots along with the seeds using a sensitive balance. e number of pods per plant and the number of seeds per pod were determined from five randomly selected plants in the net plot area at harvest. e grain yield (kg ha −1 ) was determined after threshing the sun-dried plants harvested from each net plot area, and the yield was adjusted at 10% moisture content [34] by using Dicky Johns' hand moisture tester (M20P, V-Tech, Tokyo, Japan). e moisture correction factor was obtained by the following formula: where Y is the actual moisture content measured and is the standard moisture content for pulse crops, that is, 10%. erefore, according to Birru [34], the 10% adjusted grain yield was calculated as follows: adjusted grain yield � moisture correction factor grain yield, obtained from each plot. e hundred seed weight was determined by counting the number of seeds randomly  taken from each plot, and the weight was adjusted to 10% moisture content. e harvest index (HI) was calculated as the ratio of grain yield to the aboveground dry biomass yield. e stover yield was calculated as the difference between the total aboveground biomass and grain yield.

Data Analysis.
e homogeneity of variances was evaluated using the F-test, as described by Gomez and Gomez [35]. Since the F-test showed homogeneity of the variances of the 2 years for most of the agronomic parameters, a combined analysis of variance (ANOVA) was performed using SAS version 9.4 [36] following the procedure described by Gomez and Gomez [35]. e treatment means were separated using the least significant difference test at a 5% level of significance.

Partial Budget Analysis.
A partial budget analysis was performed, as described by [37]. Since both grain and straw are important for farmers, the partial budget analysis considered the mean grain and straw yields of each treatment in 2019 and 2020. Economic analysis was performed using the prevailing average market prices for inputs at planting and outputs at the time the crop was harvested. All costs and benefits were calculated on a hectare basis in Ethiopian Birr. Total costs that varied (TCV) included the sum of the costs incurred on the purchase of NPK fertilizer, fertilizer transport cost, and the cost of application of the fertilizers. e transport costs of NPK fertilizer were estimated from market and farm gate prices. e costs of N, P, and K fertilizers were fixed based on the market price of the agricultural input supply enterprise of Ethiopia. Actual grain and straw yields were adjusted downward by 10% to reflect the difference between the experimental yield and the yield the farmers would expect to get from the same treatment [37]. us, the gross field benefit, total variable costs, net benefit, and marginal rate of return were calculated.

Physical and Chemical Properties of the Experimental Soil.
e results of the soil's physical and chemical analysis are shown in Table 2. e texture of the soil at the experimental site is silty clay. is shows that the soil has limitations in soil pH and nutrient availability for crop production. According to the soil pH rating by Murphy [38], the soil at the experimental site is strong to moderately acidic. e value is far below the range suitable for most crops and optimum for the availability of nutrients, that is, 6.5-7.5 [42]. Moreover, faba bean grows best in soils with pH values ranging from 6.5 to 9.0 [43]. us, the strong soil acidity could be one of the major factors responsible for reducing the grain yield of the crop, which is lower than the national average yield by about 100%. Hence, the application of lime, which increases soil pH, is necessary for improving faba bean yield in the experimental area. is is because lime application increases soil pH, in which Ca reacts with H + at the exchange site and neutralizes it, thereby increasing the pH of the soil [44]. e CEC value of the soil is moderate according to the rating by Hazelton and Murphy [39], which indicates its retention of adequate cations. e soil has moderate organic carbon and total N contents according to the rating by Berhanu [40].
us, this indicates the ability of the soil to supply organic carbon and mineralizable nitrogen for the proliferation of soil biota, which is important for soil biochemical processes that increase the mobility of nutrients, such as P and others, for plant uptake [45]. e available P content of the soil is low according to the rating by Cottenie [41].
us, the application of P fertilizer is required for improving faba bean yield in the study area [42]. e exchangeable K in the soil is low according to the rating by Berhanu [40], which indicates the level is not sufficient for plant growth. Since critical values for K that begin to limit plant growth are around 0.2-0.5 cmol (+) kg ha −1 according to Gourley [46], the content of exchangeable K is low, and the soil requires external application of the nutrient. us, the external application of N, P, and K fertilizers is in order.

Plant Height.
e height of faba bean plants varied significantly (p ≤ 0.05) in response to the main effects of the growing year and the rates of mineral N, P, and K fertilizers and their interactions. However, the interaction effect of the growing year with nitrogen, phosphorus, and potassium fertilizer rates was nonsignificant (Table 3). e interaction of applied nitrogen, phosphorus, and potassium fertilizers significantly (p ≤ 0.05) influenced faba bean plant height. us, the N, P, and K combination at 46: 92:30 kg ha −1 produced 23.3% more height compared with the control treatment (Table 4). e highest plant height obtained by the previously mentioned treatment was statistically similar for all N, P, and K combinations, except for a single application of these nutrients (Table 4). e result signified the vital role of combinations of the three nutrients to improve the plant height of faba bean. e possible reason for the largest plant height at the N, P, and K combination of 46:92:30 kg ha −1 could be associated with the synergic effect of N, P, and K nutrients, which enhances plant height. Particularly, the applied P and K fertilizers might facilitate more N fixation, which provides better vegetative growth and plant height of faba bean [47]. Because of the cumulative role of the nutrients in cell division, cell expansion, and enlargement, this might have ultimately contributed to the increase in the plant height of faba bean. Consistent with this result, Abou-Amer et al. [48] reported that the application of mineral N, P, and K fertilizers at 60 kg N ha −1 , 60 kg P 2 O 5 ha −1 , and 80 kg K 2 O ha −1 leads to the growth of significantly taller faba plants than plants grown with the control treatment. Similarly, Bezabih et al. [19] reported that the combined effect of rhizobium inoculation and 60 kg ha −1 K fertilization results in taller plants (55.16 cm) compared to the controls (39.83 cm). e height of faba bean also significantly (p ≤ 0.001) varied in response to the growing year (Table 3). us, a greater plant height was obtained during 2019 (123.69 cm) compared with 2020 (118.34 cm; Table 4). e increased International Journal of Agronomy plant height in the 2019 cropping year could be attributed due to the relatively suitable precipitation (1,187 mm) availability as compared to 1,378 mm in 2020 (Table 1). Consequently, more nutrient uptake from the soil may have led to the more vigorous growth of the plants.

Stem
Girth. e stem girth was significantly (p ≤ 0.001) affected by the rates of mineral N, P, and K fertilizers and their interaction. However, the stem girth did not significantly respond to the effect of the growing year and its interaction with mineral N, P, and K fertilizer rates (Table 3).
Increasing the rate of nitrogen fertilizer with the increasing rates of both phosphorus and potassium significantly increased the stem girth (Table 4). In contrast, the lowest stem girth was recorded with either nil or N, P, or K fertilizer treatment alone (Table 4). us, the N, P, and K combination at 46:92:30 kg ha −1 produced the highest stem girth (9.58 mm), and it was statistically at par with the stem girth obtained in response to the N, P, and K combinations of 46:92:60 kg ha −1 (Table 4). us, the interaction of N, P, and K at 46:92:30 kg ha −1 resulted in about 56.2% higher stem girth compared with the application of no fertilizer (Table 4). e result might be associated with the combined role of the three fertilizers in improving soil fertility and overall crop growth. e results are consistent with Ali et al., who reported that the combined application of N, P, and K of 36:54:114 kg ha −1 results in a 23% improvement in stem girth of faba bean than the controls. Bezabih et al. [19] also reported that the P and K fertilizer combination of 30: 60 kg ha −1 in rhizobium-inoculated faba bean produces a considerably wider stem girth (7.4 mm) than the controls (3.8 mm). e results imply that the synergistic effect of N, P, and K fertilizers improves the soil fertility status and enhances the stem girth [49].   International Journal of Agronomy 7

Lodging Percentage.
e lodging percentage of faba bean varied significantly (p ≤ 0.001) in response to the effect of the growing year and rates of mineral N, P, and K fertilizers and their interactions, whereas the interaction effect of the growing year with N, P, and K fertilizer rates was nonsignificant (Table 3). e minimum (0%) and maximum (39.1%) faba bean lodging percentages were recorded for potassium applied at 30 and 60 kg ha −1 and nitrogen applied at 23 kg ha −1 , respectively ( Table 4). Application of N and P without K resulted in the highest lodging percentage, whereas the increasing rate of K with N and P significantly reduced faba bean lodging. is might be due to the role of the applied K fertilizer in increasing the root growth, which improves water and nutrient uptake [50].
us, more uptakes of applied nutrients from the soil resulted in vigorous faba bean plants, which might have reduced the lodging percentage. In addition, potassium has an essential function in reducing the disease incidence of stem rot [51], which minimizes the lodging of the crop. e lodging percentage of faba bean also significantly (p ≤ 0.01) varied in response to growing seasons. Higher lodging (21.6%) was observed during 2020 compared with 2019 (12.3%; Table 4). e higher lodging in the 2020 cropping season might be linked to relatively higher precipitation (1,378 mm) than in the 2019 cropping season (1,187 mm; Table 1), which might have exposed the soil to more erosion, which might have exposed the soil to more erosion, which further affected the soil density, aggravating lodging. In addition, a higher lodging percentage in the 2020 cropping season might be associated with rainstorm damage at the vegetative stage, given the high flooding in planting in most farmers at experimental subdistricts.

Number of Branches per Plant.
e number of branches per plant responded significantly to the main effects of year of planting, mineral N, P, and K fertilizers, and interaction of N, P, and K fertilizers, while year interaction with N, P, and K fertilizers were not significant (Table 3). e number of branches per plant significantly (p ≤ 0.001) responded to the combined application of N, P, and K fertilizers. In this regard, the increasing rate of N application consistently increased the number of branches per plant with the increasing rates of both P and K applications (Table 4). e result means that the highest number of branches per plant was obtained at either medium or high rates of nitrogen (23 or 46 kg N ha −1 ) and potassium (30 or 60 kg K 2 O ha −1 ) fertilizers combined with medium phosphorus (46 kg P 2 O 5 ha −1 ) fertilizer (Table 4). However, increasing the rate of P further than 23 kg N ha −1 was no longer required for enhancing the number of branches per plant. us, the highest number of branches per plant (six) was recorded in response to the N, P, and K combination of 46:46:30 kg ha −1 (  (Table 4).
In contrast, the lowest number of branches per plant (two) was recorded from the nonfertilized treatment, which was statistically similar to plants that received nil K with N and P combinations (Table 4). e results of this study indicated that application of N, P, or K fertilizer alone or planting with nil fertilizer application would not lead to the maximum number of branches per plant. However, the interaction of N, P, and K fertilizers led to a significantly higher number of branches per plant. Phosphorus stimulates root development that is necessary for the plant to get nutrients from the soil, which facilitates legumes in fixing nitrogen in the soil through their roots [52]. Potassium is associated with the movement of water, nutrients, and carbohydrates in plant tissue, which affects ATP production, which in turn regulates the rate of photosynthesis [53]. Overall, plants that receive sufficient nitrogen have high rates of photosynthesis and typically exhibit vigorous growth that is attributed to the number of branches [54]. is implies that balanced nutrition facilitates the healthy growth and good photosynthesis and, consequently, an increased number of branches per plant [48]. In line with this result, Walled et al. [55] indicated that the combined application of N, P, and K fertilizers leads to the highest number of branches per plant of faba bean (4.61) than the controls (4.22). e results of this study also revealed that the 2019 growing season produced more branches per plant (4.25) than the 2020 growing season (3.55; Table 4). e fewer number of branches per plant in the 2020 cropping season might be attributed to higher rainfall (1,378 mm; Table 1), which increased the losses of applied nutrients through runoff and leaching [47]. Consequently, the growth and the number of branches might have been negatively affected.

Number of Pods per Plant.
e main and interaction effects of N, P, and K rates significantly (p ≤ 0.01) influenced the number of pods per plant. However, the main effect of the year, three fertilizers rate interaction with year, did not significantly influence the numbers of pods per plant (Table 5).
Increasing the rate of N application significantly and consistently increased the number of pods per plant with increasing rates of both P and K applications (Table 6). Consequently, the highest number of pods per plant (17) was produced in response to the N, P, and K combination of 46: 92:60 kg ha −1 , which was statistically at par with the number of pods per plant obtained with the N, P, and K combination of 46:92:30 kg ha −1 (Table 6). However, the lowest average number of pods per plant (4) was produced by plants that grew with nil fertilizer (Table 6). For example, the number of pods per plant obtained in response to the N, P, and K combination of 46:92:60 kg ha −1 exceeded the numbers of pods per plant produced in response to nil application of the 8 International Journal of Agronomy three fertilizers (control treatment) by 467%. is signifies that the soil of the study area is, in fact, deficient in the availability of these three major nutrients, and an adequate supply of balanced N, P, and K nutrients is required to enhance the number of pods per plant. Similar results were obtained by Ali et al. [56], who reported higher number of pods (32) in response to the combined application of N, P, and K at 36:54:114 kg ha −1 compared with the controls (5). Furthermore, Bezabih et al. [19] showed that the application of P fertilizer at the rate of 60 kg P 2 O 5 ha −1 on rhizobium-  International Journal of Agronomy

Number of Seeds per
Pod. e number of seeds per pod of faba bean varied significantly (p ≤ 0.01) in response to the main effect of the growing year and the rates of N, P, and K fertilizers applied and their interactions. However, the interaction effects of the growing year with N, P, and K fertilizer rates were nonsignificant (Table 5). e number of seeds per pod increased in response to increasing rates of N, P, and K fertilizers (Table 6). us, the highest number of seeds per pod (4.17) was recorded with the N, P, and K combination of 46:92:30 kg ha −1 , which was 52.7% more than the control treatment (Table 6). e highest number of seeds per pod obtained with the aforesaid treatment was statistically at par with the number of seeds per pod obtained with 46 kg N ha −1 × 92 kg P 2 O 5 ha −1 × 60 kg K 2 O ha −1 . e results of this study indicated that the application of a balanced N, P, and K fertilizer rate leads to a higher number of seeds per pod. In agreement with this result, Abou-Amer et al. [48] reported a significantly higher number of seeds per pod (4.63) in response to the combined application of N, P, and K fertilizers at 60:60:80 kg ha −1 than the number of seeds per pod (2.26) in the controls. e result also revealed that faba bean planted in the 2019 growing season produced more seeds per pod (3.44) than that planted in 2020 (3.19; Table 6). e higher number of seeds per pod in the 2019 cropping year might be linked to less lodging in that cropping season (Table 6), since lodging affects the number of seeds per pod by interfering with crop growth processes, such as water uptake, light interception, and healthy growth [33]. In addition, the lodging percentage was negatively associated (r � −0.26 * * * ) with the number of seeds per pod of faba bean (Table 7).

Aboveground Biomass.
e aboveground biomass yield of faba bean was significantly (p ≤ 0.001) influenced by the effects of the planting year and the interaction of N, P, and K fertilizer rates. However, the interactions of the year of planting with N, P, and K fertilizer rates did not significantly influence the aboveground biomass yield (Table 5).
Increasing the rate of nitrogen application significantly and consistently increased the aboveground biomass yield of the crop plant with increasing rates of both P and K application (Table 6). e plants with sufficient nitrogen experience high rates of photosynthesis and typically exhibit vigorous growth that contributes to higher biomass yield in faba bean [57]. In addition, phosphorus is a vital component of ATP, which is involved in overall growth from the beginning of seedlings to seed formation and maturity [58].
us, phosphorus fertilizer plays a vital role in the biomass yield improvement of the faba bean. In addition, potassium plays an important role in cell growth and is used for higher and good-quality biomass yield [59]. Overall, the enhanced production of biomass in response to the increased rates of the three fertilizers could be attributed to the physiological roles that the nutrients play in plant growth. is implies that the balanced application of the three fertilizers fulfils the growth requirements of the faba bean through synergic effects [56]. us, the highest aboveground biomass yields of the crop were produced in response to the application of 23 kg N ha −1 and 46 kg N ha −1 each combined with 46 kg P 2 O 5 ha −1 or 92 kg P 2 O 5 ha −1 plus 60 kg K 2 O ha −1 .
However, the lowest aboveground biomass yield was produced in response to the nil application of all three fertilizers. e application of N, P, and K fertilizers at the rate of 23:92:60 kg ha −1 resulted in about 4-fold biomass increment compared with the unfertilized plot. ere is a longheld belief that K is not a limiting nutrient in Ethiopian soils, yet the finding verifies that the crops respond well to the application of K with N and P fertilizers. is shows that the biomass production of the plant is enhanced in response to increasing the rates of the three fertilizers, verifying that the soil of the study area is, in fact, deficient in the availability of these three major nutrients. Consistent with this result, Fatima et al. [49] reported that the application of N, P, and K fertilizers at 20:60:40 kg ha −1 increased the biomass yield of soybean by 201% compared with the control treatment. Similarly, Bezabih et al. [19] also showed a 73% improvement in the biomass yield of faba bean due to the combined application of 30 kg P 2 O 5 ha −1 and 60 kg K 2 O ha −1 in rhizobium-inoculated faba bean. e aboveground dry biomass was significantly (p ≤ 0.001) influenced across the two cropping years. Higher biomass (7.81 t ha −1 ) was obtained in the 2019 cropping season compared with the 2020 copping season (7.29 t ha −1 ). In this regard, the growing year explained 7.1% of the total variability for the aboveground biomass (Table 6). e higher biomass in the 2019 cropping season might be linked with the positive association growth parameters such as plant height (r � 0.35 * * * ), stem girth (r � 0.44 * * * ), and the number of branches (r � 0.50 * * * ), which were higher than those in the 2020 cropping season. us, the aforementioned parameters are indirectly attributed to the increase in the aboveground dry biomass in 2019 than in the 2020 cropping season. However, the aboveground dry biomass was negatively associated (r � 0.50 * * * ) with the lodging percentage (Table 7), which was lower in the 2019 cropping season and indirectly reduced the loss in the aboveground biomass.

Grain Yield.
e grain yield of faba bean was significantly affected by the effect of the growing year and the interaction of mineral N, P, and K fertilizers. However, the interaction effect of N, P, and K rates with the growing year did not significantly influence the grain yield (Table 5).
Similar to the aboveground biomass yield, increasing the rate of N application significantly (p ≤ 0.001) and consistently increased the grain yield of the crop plant with increasing rates of both P and K application (Table 8). is means that the highest grain yields were obtained at either medium or high rates of phosphorus (46 or 92 kg P 2 O 5 ha −1 ) and potassium (30 or 60 kg K 2 O ha −1 ) fertilizers combined with the medium level of nitrogen (23 kg N ha −1 ) fertilizer (Table 8). Accordingly, the highest grain yield (4.97 t ha −1 ) was produced in response to the combined application of N, P, and K at the rate of 23:92:60 kg ha −1 . is result indicates that increasing the rate of P and K fertilizers combined with the medium rate of nitrogen improves the grain yield. However, increasing the rate of N further than 23 kg N ha −1 was no longer required for enhancing the grain yield of the crop, indicating an already optimized rate at this level of N. is result also shows that nitrogen deficiency is of a medium level in the study area and is not so severe as the deficiency of phosphorus and potassium. However, the higher grain yield obtained in response to the application of 23 kg N ha −1 compared with 46 kg N ha −1 might be associated with the requirement of legumes for little nitrogen as a starter fertilizer [13]. For example, the maximum grain yield produced by faba bean in response to the interaction effect of 23 kg N ha −1 × 92 kg P 2 O 5 ha −1 × 60 kg K 2 O ha −1 exceeded the grain yield produced in response to nil application of the three fertilizers (control treatment) by 360%, that is, about a 4-fold increment (Table 8). is shows that the grain yield was enhanced in response to increasing rates of the three fertilizers, verifying that the soil of the study area is, in fact, deficient in the availability of these three major nutrients. e maximum grain yield obtained in this study was found to be about 5-fold more than the grain yield obtained, on average, by farmers in Wolaita Zone, which amounts to 1.2 t ha −1 [5]. In this study, the grain yield obtained at nil rates of the three fertilizers was comparable to the average grain yield of farmers in the study area obtained per hectare of land, which is about 1 ton [5]. Furthermore, the faba bean grain yield obtained in this study was about 2-fold higher than the national average grain yield of faba bean, which is about 2.12 t ha −1 [2]. e enhanced yield of the crop in response to the application of the three mineral fertilizers can be attributed to the supply of balanced nutrition to the crop. us, the physiologically synergistic roles played by nitrogen, phosphorus, and potassium may have improved the growth of the plant and the partitioning of photosynthate to seeds [49]. In this connection, the applied K fertilizer may have contributed to increased root growth, thereby improving water and nutrient uptake, as suggested by Marek et al. [50]. is situation creates an ideal condition for the uptake of applied N and P from the soil. Moreover, P has an essential function in the energy metabolism of plants, playing a key role in the energy supply in the N fixation process [12]. Phosphorus and potassium supply increases the supply of N in the plant by facilitating fixation, which maximizes the vegetative growth of the crop [54]. e combined application of P and K increases the level of N in the plants, which maximizes the vegetative growth of the crop [56]. is may enhance photosynthetic efficiency and lead to the vigorous growth of the crop, which may increase the partitioning of photosynthate to the grain [60]. In line with this result, Abou-Amer et al. [48] reported that the combined application of N, P, and K at 60 kg N, 60 kg P 2 O 5 , and 80 kg, K 2 O ha −1 , respectively, produced 67.8% more grain yield than no fertilizer applied treatment. Also, N, P, and K fertilizer applied either a single or two nutrient combinations were less grain yield productivity of faba bean (Table 8). In agreement with this result, Ali et al. [56] reported that the combined application of N, P, and K at 46:92: 31 kg ha −1 improved the grain yield of faba bean by 4.7% and 1.1% than the combined application of P and K at 92: 30 kg ha −1 and P at 92 kg ha −1 , respectively.
A significantly (p ≤ 0.05) higher grain yield (3.38 t ha −1 ) was obtained in 2019 than in the 2020 cropping season. us, faba bean that was planted in 2019 had a 2.4% more grain yield than that planted in the 2020 cropping season (Table 8). Similar to the biomass yield, the grain yield result might be linked to either positive or negative associations with the growth, yield, and yield component parameters. For instance, grain yield was positively associated with plant height (r � 0.34 * * * ), thickness (r � 0.42 * * * ), number of branches (r � 0.44 * * * ), number of nodules per root (r � 0.19 * ), aboveground biomass (r � 0.97 * * * ), number of pods per plant (r � 68 * * * ), number of seeds per pod (r � 0.50 * * * ), hundred seed weight (r � 0.69 * * * ), and harvest index (r � 0.16 * ; Table 8), which were higher in the 2019 cropping season than in the 2020 cropping season. Consequently, the aforementioned parameters indirectly led to the increment in grain yield in the 2019 than in the 2020 cropping season. On the other hand, grain yield was negatively and significantly associated with logging (r � 0.38 * * * ) and stover yield (r � 0.93 * * * ; Table 7), which is lower in the 2019 cropping season and indirectly reduce the loss in the grain yield of faba bean in a cropping year.

Stover
Yield. Stover yield responded significantly to the main effect of the year of planting as well as to the main effects of mineral N, P, and K fertilizers. e applied mineral N, P, and K fertilizers also interacted to influence the variables of the plant. However, the fertilizers and growing year did not interact to influence the variables (Table 5). e increasing rate of N application significantly increased the stover yield with increasing rates of both P and K fertilizers (Table 8). Hence, a higher stover yield was obtained in response to the combination of medium or higher rates of nitrogen (23 or 46 kg N ha −1 ) and phosphorus (46 or 92 kg P 2 O 5 ha −1 ) and a higher rate of potassium (60 kg K 2 O ha −1 ). us, the application of 23 kg N × 92 kg P 2 O 5 × 60 kg K 2 O ha −1 led to the production of 430% the highest stover yields (6.83 t ha −1 ) than the unfertilized plot (Table 8).
is shows that nitrogen and phosphorus play a significant role, whereas potassium plays a more prominent role in enhancing the stover production of faba bean, since 60-70% of potassium taken up by legume plants is found in the stover [61]. However, about 30% of nitrogen and phosphorus taken up by legume crops are found in the stover, and the remaining 70% are found in the seed [62]. e results also confirm the low availability of the three nutrients in the soil of the experimental site, as indicated by the results of the chemical analysis. e result is consistent with the finding of Ali et al. [56], who reported the highest stover yield of faba bean (7.16 t ha −1 ) with the combined application of N, P, and K fertilizers at 15:22.5: 48 kg ha −1 .
e stover yield was also significantly (p ≤ 0.01) different in the experimental years in which a higher stover yield was recorded during 2019, that is, 10.8% higher than the 2020 cropping season. e higher stover yield in the 2019 International Journal of Agronomy cropping season might be due to the negative association with the lodging percentage (r � −0.35 * * * ), which was higher in the 2020 cropping season, indirectly attributed to the stover yield (Table 7), since lodging affects healthy growth [63], which results in a reduction in stover yield.
3.3.6. Hundred Seed Weight. e hundred seed weight of faba bean was significantly affected by the growing year and the interaction of mineral N, P, and K fertilizers. However, the interaction effect of N, P, and K rates with the growing year did not significantly influence the grain yield (Table 5). e increasing rate of N fertilizers significantly (p ≤ 0.05) and consistently increased hundred seed weight across the increasing rates of P. However, increasing N fertilizers across the increasing rate of K was not markedly high (Table 8). Consequently, the highest hundred seed weight (106.67 g) was obtained with the combination of N, P, and K fertilizers at the rate of 23:92:60 kg ha −1 , which was statistically at par with the hundred seed weight obtained with the combinations of N, P, and K fertilizers at the rate of 23:92:30 and 46: 92:30 kg ha −1 (Table 8). However, the lowest hundred seed weight was obtained with no fertilizer application. us, for example, the hundred seed weight obtained with N, P, and K fertilizers applied at 23:92:60 kg ha −1 was 21.2% more compared with the nontreated plot. e results indicate that phosphorus plays a critical role in increasing the hundred seed weight, followed by nitrogen. However, with the increased rate of potassium from 30 to 60 kg K 2 O ha −1 , the role appeared to have lowered in enhancing the hundred seed weight. is is consistent with the postulation that most of the P and N taken up by plants (about 70%) is found in the seeds, whereas most of the potassium is found in the stover [64]. e enhancement of the hundred seed weight in response to the interaction effect of the three fertilizers can be attributed to a balanced supply of the nutrients for uptake by the plants, as suggested by Havlin et al. [52]. e enhancement of the hundred seed weight could be linked to the synergistic roles the nutrients may have played in plant growth, which led to increased hundred seed weight [21]. For instance, K fertilizer is crucial for root growth [50], which may have created an ideal condition for the uptake of the applied N and P nutrients. In addition, phosphorus plays an important role in energy supply in the N fixation process [12] and increases seed weight [21]. Abou-Amer et al. [48] indicated that combined application of mineral fertilizers at 60 kg Nha −1 × 60 kg P 2 O 5 × 80 kg K 2 O ha −1 produces about 68% higher hundred seed weight than planting with no fertilizer application. e faba bean that was planted in the 2019 cropping season had a significantly higher hundred seed weight (103 g), which was about 6.0% more than the hundred seed weight in the 2020 cropping season (Table 8). e higher hundred seed weight in the 2019 cropping season might be associated with relatively lower rainfall (1,378 mm) than in the 2020 cropping season (Table 1). us, this may have helped the crop use nutrients efficiently for better growth, in which the seed weight increased. However, the relatively higher rainfall (1,378 mm) in the 2020 cropping season might have negatively affected the applied nutrients by erosion and waterlogging, and consequently, the assimilated translocation to the grain reduced. Hegab et al. [65] and Ahmed et al. [66] reported significant grain weight variability due to the cropping season of the faba bean. Furthermore, Stelling et al. [67] reported a difference in the hundred seed weight of the large size faba bean genotype from 112 to 178 g due to the prevailing moisture during the growing season.

Harvest Index.
e harvest index of faba bean significantly (p ≤ 0.05) varied in response to the main effects of N and K fertilizers, whereas the main effect of year and P fertilizer, the interaction of year with N, P, and K fertilizers, did not significantly influence the harvest index (Table 5).
Nitrogen applied at the rates of 23-46 kg N ha −1 resulted in a significantly low harvest index, which indicates the negative effect of N on biomass partitioning to the grains under increased N rates. e result is consistent with Khamooshi et al. [57], who reported that increased application of nitrogen from 0 kg N ha −1 to 60 kg N ha −1 significantly reduced the harvest index in faba bean (Table 9). In contrast, Ghafoor [66] reported that the mean harvest index of the faba bean increases with increasing application of combined N and P fertilizers from 0 kg N and P 2 O 5 ha −1 (40.5%) to 98 kg N and 250 kg P 2 O 5 ha −1 (43.9%). However, the increased rate of potassium from 30 to 60 K kg ha −1 was statistically at par, which indicates less photoassimilate production and ultimate partitioning into the stover compared with partitioning into the grain (Table 9). In contrast, Khalil

Partial Budget Analysis.
e adjusted grain and straw yield of faba bean were used to calculate the gross field benefits for the partial budget analysis. All N, P, and K fertilizer rate interactions were considered to analyze the partial budget of the faba bean. Furthermore, treatments having marginal rates of return (MRRs) below 100% were considered low and unacceptable to farmers and were, thus, eliminated [44]. is was because such a return would not balance the cost of the investment; however, MRR must be above 100% to have an attractive profit margin.
e maximum net benefit of USD 4,109.33 ha −1 was recorded at the N, P, and K combination of 23:92:60 kg ha −1 , followed by the N, P, and K combinations of 23:92: 30 kg ha −1 , 23:46:60 kg ha −1 , and 23:46:30 kg ha −1 (Table 10). However, the minimum net benefit of USD 917.69 ha −1 was gained on faba bean planted with no fertilizer (controls). Overall, the economical net gain at the N, P, and K combination of 23:92:60 kg ha −1 was about 4.48-fold higher compared with the control treatment (Table 10). However, the low financial capacity of smallholder farmers and the increasing price of mineral fertilizers are a cause for no-tolow inorganic fertilizer application rates in Ethiopia [64].
Economic analysis revealed that the highest mean net benefit (USD 4,109.33 ha −1 ) with an acceptable marginal rate of return (1,340%) was obtained with the N, P, and K combination of 23:92:60 kg ha −1 (Table 11). However, a   Table 11). However, the N, P, and K combination of 23:92:30 kg ha −1 was not considered due to the marginal rate of return below 100%. erefore, the N, P, and K combination of 23:92: 60 kg ha −1 was more economical than the rest of the treatments.

Conclusions
e results of this study demonstrated that the phenology, growth, yield, yield components, and nutrient use efficiency are significantly influenced by the interaction of N, P, and K fertilizer rates. A maximum faba bean grain yield of 4,970 kg ha −1 was recorded in response to the N, P, and K fertilizer combination of 23:92:60 kg ha −1 . e result was in statistical parity with the grain yield obtained in response to the N, P, and K fertilizer combinations of 23:46:30 kg ha −1 and 46:46:60 kg ha −1 . However, the lowest grain yield (1.08 t ha −1 ) was recorded in faba bean planted with no fertilizer application. Overall, the grain yield produced in response to the N, P, and K fertilizer combination of 23:92: 60 kg ha −1 was about 4-fold higher than that with no fertilizer treatment and the average grain yield obtained by the farmers per hectare in the study area. In addition, the faba bean grain yield obtained in this study was about 2-fold higher than the national average grain yield of faba bean, which is about 2.12 t ha −1 . In general, for a better grain yield (4.97 t ha −1 ) and the highest mean net benefit (USD 4,109.33 ha −1 ) with an acceptable marginal rate of return (1,340%), farmers are advised to use N, P, and K fertilizers in a combination of 23:92:60 kg ha −1 . In Ethiopia, there is a belief that the soils contain enough or sufficient quantity of the K nutrient. Consequently, management has focused on N-and P-containing fertilizers. However, this finding suggests the need for balanced fertilization including K. Meanwhile, further research that evaluates the integration of mineral N, P, and K and organic fertilizers is suggested for soil health and higher productivity.

Data Availability
Data will be made available upon request to the corresponding author.

Conflicts of Interest
e authors declare that there are no conflicts of interest.

Authors' Contributions
B. A. T. contributed to the conceptualization, methodology, software provision, formal analysis, investigation, data curation, visualization, original draft preparation, funding, and administration, and N. D., T. T., F. L., and Y. A. contributed to validation, resource provision, review and editing, and supervision. All authors have read and agreed to the final version of the manuscript.