Response of bread wheat to sulfur and phosphorus fertilizers in the north central Ethiopia

Emerging research evidences since few years back are indicating that sulfur (S) is becoming a limiting nutrient in some Ethiopian soils. However, these evidences are not sufficient to make a solid conclusion that some soils of Ethiopia have became deficient in S. There is also limited information on the interaction effects of S and phosphorus (P) on bread wheat in Ethiopia. Therefore, an experiment was conducted to evaluate the effects of S and P fertilizers; and their interaction on yield components and yield of bread wheat grown at Gerba and Deneba locations, Northern central Ethiopia. A factorial experiment consisting of three levels (0, 15 and 30 kg ha−1) of S and four levels (0, 11, 22 and 44 kg ha−1) of P was laid out in RCB design with three replications. Results revealed that yield components and yield of wheat were significantly affected by both main and interaction effects of S and P fertilizers in both locations. At Gerba, S applied at 15 and 30 kg ha−1 increased grain yield (GY) of wheat by 32 and 44% over untreated control treatment, respectively. The corresponding increases at Deneba were 29 and 37% over untreated control treatment, respectively. However, significantly higher GY was obtained with treatments involving combined application of S and P than that obtained with single application of S or P. GY gains due to S and P interaction effect ranged from 46–65 to 52–75% over untreated control treatment at Gerba and Deneba, respectively. Optimum GY of wheat was obtained with treatment involving 15 kg S ha−1 + 22 kg P ha−1. The results of this experiment revealed that application of P and S fertilizer has significantly increased yield component, grain and straw yield of wheat compared to unfertilized control plot, indicating insufficient soil P and S content for optimum production of wheat and this was confirmed by very low soil test values of S in both sites. Combined application of S and P produced significantly higher yield of wheat than that obtained with single application of S or P indicating synergistic interaction between these nutrients. In all cases, optimum grain and straw yield of wheat was obtained with treatment involving at 22 P + 15 S kg ha−1, while a partial budget analysis result revealed that a combination of 22P and 15S kg ha−1 produced the highest MMR (54.9%). Thus this treatment is found to be economically feasible treatment for bread wheat production in study area of the district.


Background
Bread wheat (Triticum aestivum L.) is one of the most important cereal crops produced and consumed in Ethiopia. It stands third among cereals in terms of area cultivated annually and grain production next to barley and tef [1]. Traditionally, wheat grains are used to prepare household bread, beverage and pancake.

Open Access
Agriculture & Food Security *Correspondence: wassiehaile@yahoo.co.uk 2 College of Agriculture, Hawassa University, P. O. Box 05, Hawassa, Ethiopia Full list of author information is available at the end of the article Despite, its importance and growing demand for bread wheat in Ethiopia, its production and productivity are desperately very low, about 2.5 t ha −1 which is below the world's average of 3.4 t ha −1 [2]. Consequently, the country imports large amounts (30-50% of total annual demand) wheat grains every year from abroad to meet domestic demand [3] which has grave consequence on foreign currency reserve of the country.
Declining soil fertility is one of the major factors that account for low productivity of wheat in Ethiopia [4,5]. Declining soil fertility which is caused by soil erosion, continuous cropping of same land year after year, deforestation, inadequate replenishment of nutrients lost through different pathways through application of organic and inorganic fertilizers are the main challenge not only to wheat production but also to the production of all crops [4,6]. It is an issue of great concern in Ethiopia as soil nutrient depletion is becoming severe and severe with time since little efforts are being made to reverse the problem. Food security is at the national and global level tends to focus on the supply side of the food equation. The question raised is there enough food available, where food is usually interpreted to mean dietary energy. But availability does not assure access, and enough calories do not assure a healthy and nutritional diet [7]. Among several composite indexes for food security is a large set of proxies for food security [8].
In an effort to overcome declining soil fertility challenge and thereby improve crop productivity including wheat, application of inorganic fertilizer containing nitrogen (N) and phosphorus (P) has long been started in Ethiopia. This is because based on the results of extensive field and laboratory studies some 50 years ago, N and P were identified as the only nutrients deficient in almost all arable lands of Ethiopia despite variations from place to place [9]. Appreciable increases in the yields of several crops used to be obtained due to these inputs [10][11][12]. However, crop yield gain due to N and P fertilizer application is declining over time despite steady increases in fertilizer consumption in Ethiopia [13]. Declining crop yield responses to N and P fertilizers is attributed to decreasing soil organic matter (SOM) content [13]. Moreover, depletion of other nutrients in addition to N and P could be additional factor for decreasing response of crops to N and P fertilizers [9]. Sulfur (S) could be one of the most likely limiting nutrients in Ethiopia soils. It is a building block of protein and a key ingredient in the formation of chlorophyll [14]. It has been reported that S-deficient plants exhibits reduced plant height and stunted growth, reduced tillers, spikelets and delayed maturity [15,16]. Sulfur-deficient plants are shown to be less resistance under stress conditions [17].
In this regards, emerging research evidences are showing that sulfur (S) is one of the nutrients becoming deficient in some Ethiopian soils limiting crop production. For instance, Assefa [18] studied the response of wheat to S application on 18 sites located in Arsi, East Shewa and West Shewa and reported that wheat significantly responded to S fertilizer application in 72% of experimental sites. He further reported that soils of responding sites had S content below critical level (11-13 mg kg −1 SO 4 −2 -S) for optimum production of wheat. Similarly, Kiros and Singh [19] studied the effect of sulfur on two varieties of wheat (Enkoy and Shehan) grown in Andisol and Cambisols in semi-arid region of Ethiopia and found that S application significantly increased grain yield, Nuse efficiency and S uptake of both varieties of wheat.
There is little or no information so far on the response of wheat to S application in Moretena-Jiru district of north Shewa, Ethiopia. Therefore, the designed objective of this experiment was to determine effects of S and P application on the growth, yield components and yield of wheat and to determine the economic feasibility of S and P fertilizer application for wheat production in the study district.

Description of the study locations
The experiment was conducted in 2017/2018 cropping season on two farmers' fields at Gerba and Deneba locations which are 17 km apart each other and they are found in Moretena-Jiru district of north Shewa, about 195 km northwest of Addis Ababa. Geographically, the study locations lie between 9°46′53.3''-9°52′07.5'' N and 39 0 10′38.1''-039°11′28.3'' E. The altitude of the study locations ranges from 2600 to 2700 m.a.s.l with mean annual minimum and maximum temperature ranging from 5.9-11.6 ℃ to 17.4-31.2 ℃, respectively. The study locations and the district as a whole are characterized by having a uni-modal rainfall pattern and receives an average annual rainfall of 921.2 mm. Vertisols is the dominant soil type occurring in the experimental locations and the district. Major crops grown in the Gerba and Deneba locations; and in Moretena-Jiru district as well include wheat, tef, lentil, faba bean, and chickpea, field pea and grass pea in decreasing orders of area cultivated under these crops.

Soil sampling and analyses
After selecting the experimental sites, pre-planting soil samples were collected from each site for the analyses of selected physicochemical properties. Composite soil samples were taken from each site from a depth of 0-20 cm using augur randomly from 15 spots by walking in a zigzag pattern. After thoroughly mixing the composite samples, 1 kg of sub-sample was taken and brought to Debre Birhan agricultural research centre soil laboratory where it was air dried and grounded to pass 2 mm mesh sized sieve.
The processed samples were analyzed for texture following by Bouyoucos hydrometer method [20]. The pH of the soil was measured using pH-water method by making soil-to-water suspension of 1:2.5 ratio and was measured using a pH meter. The soil OC content was determined by wet digestion method [21]. Total nitrogen (TN) was determined by using the modified micro-Kjeldhal method [22], available P (avaP) was analyzed by using Olsen's calorimetric method as described by Olsen et al. [23], and exchangeable K in the soil was extracted with 1 N NH 4 OAc and the amount was estimated using a flame photometer [24]. Soil available sulfur (avaS) was determined turbidimetrically using a spectrophotometer [25]. CEC was measured after saturating the soil with 1 N NH 4 OAc and displacing it with 1 N NaOAc [26].

Treatments, design and experimental procedure
A factorial experiment consisting of three levels (0, 15 and 30 kg ha −1 ) of S and four levels (0, 11, 22 and 44 kg ha −1 ) of P was laid out in RCB design with three replications. The above S levels of 0, 15 and 30 kg ha −1 were as coded S0, S1 and S2, respectively. Similarly, P levels of 0, 11, 22 and 44 kg ha −1 P were coded as P0, P1, P2 and P3, respectively. K 2 SO 4 and TSP were used as S and P sources, respectively. The test crop, wheat variety, Menzaie was sown in a unit plot size of 3.6 × 3.4 m with row spacing of 20 cm apart at a rate of 150 kg ha −1 . The whole doses of K 2 SO 4 and TSP fertilizers were applied basal applied in both sides of rows just before sowing as per the treatment. To compensate for K applied along with S in S treatment plots, K was applied to all non-S treatment receiving plots in the form of KCl. The recommended dose of N fertilizer for wheat production around the study area which is 167 kg ha −1 [27] was applied uniformly to all plots in the form of urea. The urea-N was split applied in which one half of N was applied at planting and the remaining one half was applied one month after planting and first weeding. All agronomic management of the trials was done as per the recommendation for the crop [28].

Data collection and measurements
Data on plant height (PHT) at harvest, number of fertile tillers per meter square (NFTPM), at harvest, spike length (SPL), number of spikelets per spike (NSPS), number of grains per spike (NGPS), thousand grain weigh (TGW), grain yield (GY), straw yield (SY) and harvest index (HI) data were collected at appropriate stage and time pertinent for each parameter.

Data analysis
The collected data were subjected to ANOVA. After verifying the homogeneity of error variances, combined analysis of variance was done using the procedure of SAS software version 9.4 [29]. Mean comparisons were done by least significant difference (LSD) at the 5% level of significance.

Soil physicochemical properties
Pre-planting soil analyses data of selected physicochemical properties of samples collected from experimental location at Geraba and Deneba are summarized in Table 1. The soils of both sites were belonging to clay textural class. Gerba soil has soil reaction near to neutral, whereas the soil of Deneba is slightly alkaline reaction. The SOM and TN content of soils of both locations are in low categories according to Jones [24]. The available P content of both locations are slightly higher than 6 mg kg −1 , lower than the critical soil available P value established for some Ethiopian soils which is 10 mgkg −1 [30]. The S content of both soils were also found to be below the critical value of S required for optimum production of wheat which is 10-13 mg kg −1 SO 4 −2 -S [31]. The CEC of both soils was also in high category according to rating established by Landon [32].

Effects of sulfur and phosphorus fertilizers on growth of bread wheat
Main effects of S and P on PHT and NFTPM of wheat in Gerba and Deneba locations are presented in Table 2. S1 and S2 treatments significantly increased PHT of wheat in both locations compared with control (S0), but they were statistically at par with each other. This is in agreement with the findings [33] who reported that application of S significantly increased PHT of wheat. PHT was also significantly increased by main effects of P application. The highest PHT of was produced by P3 which was statistically at par equal to P2 followed by P1 and the least PHT was produced by P0 in both locations. NFTPM which are one of the major yield contributing factors [34] were significantly affected by main effects of both S and P fertilizers (Table 2). S1 and S2 increased NFTPM by 8.5 and 4% over S0 in Gerba, respectively. Similarly, in Deneba, S1 and S2 increased NFTPM by 8 and 4.5% over S0, respectively. In line with this, Assefa [18] reported that application of S significantly increased number of tillers per plant of wheat grown in Keteba and Bekejo sites of east Shewa zone of Ethiopia. According to him, S applied at 10 kg ha −1 produced the highest number of tillers per plant relative to the control.
Phosphorus has also significantly increased NFTPM in both locations. In Gerba, P2 produced the highest NFTPM followed by P3, P1 and P0 in that order whereas in Deneba, all P rates have significantly produced higher NFTPM relative to P0, but they were statistically at par among each other. The apparent significant increase in NFTPM of wheat in both locations is in agreement with Ibrahim et al. [35] who reported increased in number of fertile tillers in wheat due to increase in P levels over control.
There was also significant interaction effect between S and P on PHT in both locations (Table 3). Significantly higher PHT of wheat was produced treatment involving combined applications of S and P than that obtained with control treatment and treatments with single application of S or P in both locations. The present result is in agreement with Kumar et al. [36] who reported that combined application of P and S produced significantly higher PHT of wheat than that produced by single application either S or P. But S by P interaction was not significant for NFTPM of wheat in both locations indicating that change in the levels of S doesn't change the effect of P levels on NFTPM and vice versa.

Effects of sulfur and phosphorus on yield components of wheat
The main effect results revealed that both S and P fertilizers have significant effect on yield components of wheat grown in Gerba and Deneba locations (Table 4). Both S1 and S2 have significantly increased SPL, and NGPS of wheat relative to S0 except SPL at Deneba. However, there was no significant difference between S1 and S2 rates with regards to these parameters of wheat in both sites. Significantly higher NSPS was recorded on S2 on Gerba and S1 on Deneba as compared to S0 but there was no significant difference between S1 and S2. Significantly higher TGW was recorded on S1 than S0 and S2 on both locations. In line with this finding, Togay et al. [37] studied the effects of different levels of S on barley in Turkey and found that application of S has significantly increased the above parameters relative to the control. Similarly, all P rates have also significantly increased SPL, NGPS and TGW of wheat in both locations relative to the control (P0) except NSPS at Gerba. But there were no significant differences among P rates with respect to these parameters of wheat indicating that optimum levels of yield components of wheat can be achieved at P1  rate in both sites. In soils with low available P contents such as at Gerba and Deneba locations in this case, application of P significantly increases yields of crops including wheat through its positive impacts on the growth and yield components of crops [38]. The interaction effects of S and P were non-significant on SPL, NSPS, NGPS in both locations suggesting that the change in the levels of S was not affected by the change in the levels of P with respect to these parameters and vice versa. However, there was a significant interaction effect between S and P on TGW in both locations (Fig. 1). The result implies that the effect of changing levels of S on TGW was significantly changed with a change in the levels of P and vice versa. In both locations significantly higher TGW of wheat was obtained with treatment involving combined applications of S and P and the highest TGW was obtained with treatment involving S + P applied at 15 and 22 kg ha −1 , respectively, in both locations.

Effects of sulfur and phosphorus on yield of wheat
Both grain yield (GY) and straw yield (SY) of wheat were significantly affected by main effects of S application in both Gerba and Deneba locations (Table 5). At Gerba, S1 and S2 treatments have significantly increased GY by 13 and 12% over the control (S0), respectively, and SY by 39 and 38% over S0, respectively. Similarly, at Deneba, S1 and S2 rates have significantly increased GY by 19 and 20% over S0, respectively and same treatments increased SY by 24 and 26% over S0, respectively. The apparent significant response of wheat to S application in both locations found in this study indicates low level soil S for optimum production of wheat in the study areas. This is substantiated by low level or below critical soil test S values (Table 1) in both soils for optimum crop production. These results are in agreement with the finding of Assefa [18] Who studied the response of wheat to S application on 18 sites located in Arsi, east Shewa and west Shewa and reported that wheat significantly responded to S fertilizer application in 72% of experimental sites. He further reported that soils of responding sites had S content below critical level (10-13 mg kg −1 SO 4 −2 -S) for optimum production of wheat. In another study by Khan et al. [39] reported that S applied at 20 kg ha −1 at stem elongation stage significantly increased GY of wheat in Peshawar, Pakistan by 28.5% over untreated control. According to DeRuiter and Martin [40], wheat yield can be increased by up to 42% due to S fertilizer depending on the inherent S level in a particular soil. The further noted that usually best crop response to S is obtained at application rates between 10 and 20 kg ha −1 , even though much higher application rates S (49-60 kg ha −1 ) have been reported in highly calcareous soils in India [41]. The main effect of S application on HI was significant at Gerba but not at Deneba locations (Table 5). At Gerba, both levels of S have significantly increased HI relative to the control but they were statistically at par with each other. This is in line with the finding of Togay et al. [37] who reported that application of S significantly increased HI in barley.
GY and SY of wheat were also significantly affected by main effect of P application in both locations ( Table 5). All P levels produced significantly higher GY and SY compared to control in both sites. The GY and SY increasing effects of P levels were in the order of P3 = P2 > P1 > P0. This suggests that the soil of the study sites is low in their available P content for optimum growth and productivity of wheat. This was substantiated by low level soil test P values of both sites shown in Table 1. This is in line with reports of several researchers that P is the second most limiting nutrients in almost all arable soils of Ethiopian [4,9,30].
The interaction effects of S and P on GY and SY of wheat are shown in Fig. 2. The grain yields as influenced by different levels of P and S are presented in Fig. 2. The addition of 22 kg P ha −1 with 15 kg S ha −1 produced maximum grains yield (3526.8 kg ha −1 ) which was significantly different from other treatment combination except 22P + 30S, 44P + 15S and 44P + 30S kg ha −1 . This treatment has a yield advantage of 40.1% over control treatment. However, grain yield become decreased numerically but not significantly as interactive application rate of P and S beyond 22P and 15S kg ha −1 . The minimum grain yield (2076.7 kg ha −1 ) was recorded in control plot (untreated with P and S) shown in Fig. 2. Even though, the overall grain yield of the wheat was lower compared to potential yield of wheat, the treatment combinations showed that P and S at the rate of 22 and Table 5 Main effects of S and P fertilizers on yield (tha −1 ) of wheat at Gerba and Deneba locations a Means within column followed by same letter(s) are not statistically different from each other at 0.05 probability level, **significant at 0.01 probability level, ns, nonsignificant 15 kg ha −1 , respectively, was an advantageous. The straw yield of wheat was also influenced significantly (P ≤ 0.01) by different levels of P and S application interactively. The higher straw yield (5120.0 kg ha −1 ) was produced when 22 kg P ha −1 and 15 kg S ha −1 were applied interactively (Fig. 2) which have straw yield boosted by 53.4% over control treatment. Similar trends were also observed like grain yield due to the influence of different P and S application rates. The minimum straw yield (2385.6 kg ha −1 ) was recorded in control plot.
In agreement with the current finding, both synergistic and antagonistic interactions between S and P have been reported in literature [42][43][44]. For instance, Uddin et al. [45] reported that there was a synergistic interaction between S and P combined applied at 20 and 40 kg ha −1 , respectively, on yield of rice grown in Joydebpur, Bangladesh. According to these authors, the above treatment increased GY and SY of rice by 61 and 33% over the control, respectively. But combined application of S and P above these rates resulted in antagonistic interaction effect on rice as was revealed by declining yield response. Similarly, Islam et al. [46] reported synergistic effect of these nutrients on chickpea when 15 kg ha −1 S was combined applied with 40 kg ha −1 P. But they found antagonistic effect when 30 kg ha −1 S was combined applied with 80 kg ha −1 P.
The synergistic interaction effect of S and P that occurs at specific combinations levels of these nutrients on crops including wheat could be attributed to increased uptake and utilizations of soil nutrients due to welldeveloped roots and shoots [47]. On the contrary, antagonistic interaction effect between S and P that may occur at higher dose of S when it is combined applied with P or vice versa could be attributed to competition between the two nutrients for uptake sites on plant roots as both nutrients are taken by roots as negatively charred ions (H 2 PO 4 1− /HPO 4 2− -P and SO 4 2− -S). This is substantiated by the finding of Muhammad et al. [47] who reported that application of high dose of S (75 kg ha −1 ) resulted in decreased uptake of P by maize plant and vice versa ultimately resulted in decrease grain yield. In general, the interaction between S and P could be synergistic or antagonistic depending on the initial soil fertility status, the amounts of nutrient applied and types of crop [48]. The current result has similar trend with Islam [35] reported that application of S and P at higher rate (40 kg S and 35.2 kg P ha −1 ) showed reduction in yields. Other findings [49] revealed that combined application of P and S at 50 kg P 2 O 5 and 20 kg S ha −1 improve yield of wheat over control. Similarly, the interaction of P and S levels increased their beneficial synergistic effect and applying of the individual Sulfur levels with P level produced increases in the grain yield of wheat over the control [50].

Partial budget analysis
The data related to partial budget analysis are given in Table 6. It is clear from the data that maximum net benefit (45351.7 ET Birr) was obtained from combined Depend on dominancy analysis, treatments of S and P at 15 kg S without P, 44 kg P without S, 30 kg S without P, 30 kg S with 22 kg P, 15 kg S with 44 kg P, 30 kg S with 44 kg P, 30 kg S with 11 kg P ha −1 and control were dominated by the rest four treatments and they were also excluded from further economic analysis. The data regarding the marginal rate of return (MRR) revealed that maximum MRR (5490%) was obtained when S and P were applied at the combined rate of 15 kg S and 22 kg P ha −1 followed by rate of 11 kg P ha −1 without S (2502%). Minimum MRR was (175%) recorded in treatment where application of P at 22 kg ha −1 without S.
Data from Table 6 clearly revealed that the non-dominated treatments associated with MRR are greater than 100%. This implies that the four non-dominated treatments are economically feasible alternative to the other dominated treatments. MRR, 5490% means the producer obtained an additional income of 54.9 Ethiopian birr per a unite cost they have invested. Generally, treatment combination of S and P at 15 kg S and 22 kg P ha −1 gave better MRR value relative to the other four non-dominated treatments and profitability can be optimized by using this treatment.

Conclusion
The results of this experiment revealed that application of P and S fertilizer has significantly increased yield component, grain and straw yield of wheat compared to unfertilized control plot, indicating insufficient soil P and S content for optimum production of wheat and this was confirmed by very low soil test values of S in both sites. Combined application of S and P produced significantly higher yield of wheat than that obtained with single application of S or P indicating synergistic interaction between these nutrients. But this was true for treatments involving combined applications low level of S and increasing levels of P or vice versa. Otherwise, there was decreasing trend of yield gains with treatments involving combined applications at higher levels of both nutrients. In all cases, optimum grain and straw yield of wheat was obtained with treatment involving at 22 P + 15 S kg ha −1 . While, partial budget analysis result revealed that, combination of 22P and 15S kg ha −1 produced the highest MMR (54.9%) and thus, this treatment is found to be economically feasible treatment for bread wheat production in study area of the district.

Acknowledgements
The senior author of this manuscript would like extend his thanks to Ato Lisanu Getaneh and Ato Getaneh Shegaw, technicians in the soil laboratory of Debre Berhan agricultural research center for their great support in analyzing soil sample.
Authors' contributions SAA generated the research idea, drafted into proposal containing details of background information, objectives, treatments, experimental design, etc. After defending his proposal and getting approval, he conducted the field and laboratory experiments, collected all the necessary data, analyzed the data and wrote his thesis. WHW served as major supervisor of SAA and played a role in guiding the research starting from shaping the proposal up to reviewing and edition the draft thesis of SAA. WTB was co-supervisor of SAA and involved in commenting and reviewing the draft thesis of SAA. This manuscript was extracted from the thesis of SAA. All authors read and approved the final manuscript.

Funding
This study was financed by Amhara region agricultural research Institute and we are very grateful to it.

Availability of data and materials
All relevant data generated from the study and reported in the manuscript are included in this article. Further data sets are available from the senior and/or corresponding author upon request.

Treatment
Adj