Optimum Rate and Deep Placement of Nitrogen Fertilizer Improves Nitrogen Use Efficiency and Tomato Yield in Nepal

: In Nepal, blanket fertilizer recommendations without considering diverse soil types, nutrient status, climate and crop management practices along with imbalanced fertilization practices by farmers, mainly “urea fertilizer,” have resulted in reduced nitrogen use efﬁciency (NUE) and productivity in tomato production. Optimizing the rate of nitrogen (N) fertilizer, application time and improved application methods could increase crop yields and NUE and reduce environmental costs. This study was conducted to identify the optimum N rate and application method for increased tomato yield and NUE. Multilocation trials ( n = 28) conducted in a randomized complete block design with nine treatments across ﬁve districts included the omission of N, P and K (N0, P0, K0), variable N rates of 100, 150, 200 and 250 kg ha − 1 (N-100, N-150, N-200 and N-250), use of urea briquettes (UB) with deep placement (UBN-150) and a control (CK). N input in UB was reduced by 25% from the recommended N rate of 200 kg ha − 1 considering its expected higher NUE. Yield responses from an NPK omission plot revealed N as the most limiting plant nutrient. Applications of fertilizer at N-100, N-150, N-200 and N-250 increased tomato yield by 27%, 35%, 43% and 27%, respectively, over N0. Tomato yields responded quadratically to the added N fertilizers with optimum rates ranging from 150 to 200 kg ha − 1 across districts. UBN-150 signiﬁcantly increased tomato yield by 12% over N-150 and produced a similar yield to N-200 (the recommended rate). The highest partial factor productivity of nitrogen (PFPN) was observed at N-100 and the highest agronomic efﬁciency of N (AEN) was at N-200. Deep placement of UB at-150 increased PFPN by 8% and 21% and AEN by 27% and 21% compared with N-150 and N-200, respectively. These results have positive implications for developing efﬁcient N fertilization strategies to increase tomato yields and reduce environmental impacts in Nepal.


Introduction
Tomato is grown as commercial produce by Nepalese famers year-round across different agro-climatic zones under both polyhouse and open-field conditions [1,2]. Tomato ranks third after cauliflower and cabbage in area used (20,000 ha) and production (0.3 million metric tons, t) [1]. Despite the large potential for domestic production, the current production does not satisfy national demand and a large amount of tomatoes are imported from India every year [2,3]. Poor irrigation and inefficient nutrient management practices are among the major reasons for lower productivity in Nepal [3][4][5]. Efficient management of soil nutrients such as nitrogen (N), phosphorous (P), potassium (K) and micronutrients such as zinc (Zn), boron (B) and optimum soil moisture providing timely irrigation play a significant role in improving tomato yield.
The major N fertilizer for tomato is prilled or granular urea, which is highly watersoluble and less efficient, particularly when broadcast-the common practice in Nepal. In the broadcasting method, farmers spread urea manually on the surface of the entire field uniformly during planting time as a basal application and top dressing, resulting in a nitrogen use efficiency (NUE) of 20-35%. This means that more than 50% of the applied N is lost to the environment through leaching, surface run-off, volatilization, nitrification and denitrification [6][7][8]. The official recommendation of the Nepal Agricultural Research Council (NARC) for fertilizer use is a blanket type; that is, a uniform rate of N, P and K across diverse agroecological zones that encompasses varied soil types, climates (rainfall and temperature) and crop management practices. This blanket recommendation is due to the limited number of on-farm fertilizer trials and crop management information practiced by farmers in different agricultural domains in Nepal. The fertilizer rate recommendations were also based on the results of on-station trials, with very few trials conducted on farms, and yield responses and NUE from on-station trials will differ from on-farm results obtained under diverse cropping systems and crop management practices [9]. On-station trials are managed by scientist following best management practices (BMPs): tillage practices using machine (soil tilth improvement), line planting methods to maintain uniform plant density, accurate amount of agri-inputs application (fertilizer, irrigation, pesticides, etc.) and close supervision of the field, whereas on-farm trials are managed by farmers with a lack of proper agronomic practices: broadcasting seed and fertilizers, the use of inaccurate agri-inputs and less supervision of the farm, thus lacking BMPs. These improved crop management practices in on-station trials produced higher crop yields compared with on-farm trials [9]. It is thus critical to measure fertilizer responses in farmers' fields and recommend site-specific fertilizer rates and BMPs at the local farm level [10].
Nitrogen is one of the most limiting nutrients, affecting plant growth and yield worldwide [11,12] and in Nepal [13]. Optimum N fertilization is integral to meeting plant nitrogen demand and increasing yields. However, rather than increasing yields, N supplies that are beyond plant demands lead to reduced NUE, higher production costs and large N losses to the environment via leaching and emissions of nitric (NO) and nitrous oxides (N 2 O), the latter a potent greenhouse gas [14][15][16]. Optimizing N fertilization is vital to increase yields and address environmental concerns [17]. Application of N fertilizer at the right rate with multiple splits based on plant nutrient demand has been shown to increase NUE and tomato yields [14,18]. Previous studies have reported that fertilizer management using the 4Rs (right source, right rate, right time and right place) of nutrient stewardship enhances crop productivity and NUE [19].
Recently, improved application methods such as urea deep placement (UDP) or root zone application has been found to increase NUE and yields and reduce N losses to the environment [20,21]. For convenient applications, normal urea is compressed into 1-3 g pellets, commonly called a "urea briquette (UB)," and these are applied in the root zone 7-10 cm below soil surface. Deep placement of UBs retains N longer and releases it slowly over the cropping period, based on the plant nutrient demand [22,23]. UBs can be placed in a single application and at a 20-30% lower N rate than the recommended rate, given its slow release and higher NUE [20,24]. Studies have reported that a single UB application with less N has produced similar or even higher vegetable (tomato, cauliflower, bitter gourd) yields than with the recommended N rate [21,24,25].
Nitrogen is a crucial nutrient for the physiological and metabolic process in a tomato; adequate N availability increases marketable yield [25]. So far, few on-farm studies have explored optimal N rates and practices for tomato production in Nepal. In this study, the agronomic effects of variable N rates and UB with deep placement were tested on tomatoes in farmers fields under polyhouse conditions with drip irrigation across five midhill (subtropical climate) districts to identify the most limiting plant nutrient, determine optimal N rates and assess the effects of UB deep placement on increasing tomato yields and NUE across sites.

Description of the Study Site
Field trials were conducted under plastic tunnels in five districts (Doti, Palpa, Kavre, Dang and Surkhet) in the mid-hill region (subtropical climate) of Nepal ( Figure S1). Average minimum and maximum temperatures during the growing period were 11.9 • C (ranging from 3.6 to 21.5 • C) and 23.3 • C (ranging from 14.4 to 31.0 • C) across districts. During the growing period, temperatures were higher in Dang and Surkhet than the other three districts (Palpa, Doti and Kavre) ( Figure S2). Average total rainfall across districts was 365 ± 82 mm, ranging from 251 to 459 mm ( Figure S2). Elevation across the study sites ranged from 497 to 1828 meters above sea level.

Soil properties
Soil samples from three locations in each experimental site were collected at a vertical depth of 10 cm and pooled into a composite sample for each site across districts (n = 28). Crop residues, leaf litter and pebbles were removed from the soil surface before collecting the samples. Collected soil samples were mixed thoroughly, removing the foreign materials (pebbles, stones, roots and other plant parts) and kept in zip-lock plastics with proper labelling (farm ID, GPS location, farmers name). Soil samples were oven-dried at 40 • C for three days and passed through a 2 mm sieve before analysis. Soil pH, organic matter (OM), texture, total N, available P 2 O 5 and available K 2 O were determined using both wet chemistry and spectroscopy as described by Pandit et al. [24]. Soils were characterized as a moderately acidic silty loam (pH ranging from 5.9 to 6.6). Soil OM ranged from 16.0 to 31.7 g kg −1 , total N from 0.1 to 1.7 g kg −1 , available P 2 O 5 from 21.3 to 64.6 mg kg −1 and available K 2 O from 86.3 to 224.6 mg kg −1 ( Table 1).

Fertilizer Treatments, Research Design and Crop Management
Multilocation trials (n = 28) were conducted under plastic tunnels across five districts: Doti (n = 6), Palpa (n = 4), Kavre (n = 6), Dang (n = 6) and Surkhet (n = 6) during September-October 2018 (see supplementary materials). Trials were arranged in a randomized complete block design (RCBD) consisting of nine treatments (Table 2) Table 2. The N rate in UB was reduced by 25% from the recommended N (200 kg N ha −1 ) to account for the expected higher N use efficiency with the deep placement of UB [20,24]. During land preparation, a trial plot in a plastic tunnel (12 m × 5 m) was plowed 2-3 times and pulverized well before planting. The plot size was 4.8 m 2 (4 m × 1.2 m). Phosphorous and Potassium fertilizers were applied as a basal dose at a rate of 150 kg ha −1 each (5 cm depth and 5 cm apart from planting line) in the form of DAP and MOP, respectively. Each plot was mulched with UV plastics and holes were made to transplant seedlings (25-30-day-old seedlings at the 3-4 leaf stages) at a spacing of 60 cm x 40 cm. In each plot, 20 plants were grown. The hybrid tomato variety Srijana was used. The Srijana variety is recommended for the tropical (terai) and sub-tropical (mid hills) climates in Nepal with 70-80 days maturity and has a yield potential of 105 to 110 ton ha −1 . Urea-N was top-dressed in two equal splits as per treatment (Table 2) 20 days after transplanting and fruit bearing. In UB treatments, two UBs of 2.7 g (total 5.4 g) were applied to each plant at a depth of 7-10 cm and 5 cm apart during the first top dressing 20 days after transplanting. Drip irrigation was provided immediately after transplanting and sufficient moisture maintained until two weeks after transplanting. Later, irrigation was provided as needed and during critical moisture requirement stages (flowering and fruiting). Staking (bamboo sticks of 170 cm long) was provided in each row between four plants for support and to promote growth. Pruning was carried out, allowing one stem (single central stem) or two stems (double stem) for fruit production. For single central stems, buds coming from the nodes were continuously removed and, for double stems, all buds below 30-40 cm were removed, allowing one bud to grow at 40-45 cm above the ground in addition to the central steam. Hand weeding, earthing up and pest disease management were conducted as required.

PFPN and AEN
Tomatoes were harvested manually when the fruit turned yellow-red from three quadrants of 1 m 2 (3 m 2 in total) in each treatment plot. Fruit was harvested 5-6 times per plant and the yield measured in each harvest. Total yield was recorded after final harvest in each plot.

Statistical Analysis
Data were analyzed using R software, version 3.6.2 (R Core Team 2019, Vienna, Austria). A linear mixed-effect model was used to test the effect of fertilizer treatments (fixed factor) on tomato yields across districts (random factor) using the 'lmer' function in the 'lme4' package, fitted by the Restricted Maximum Likelihood (REML) model. One-way ANOVA was performed to assess the fixed effect of variable N rates applied in the form of conventional urea and UBs at 150 kg N ha −1 on PFPN and AEN. A post hoc Tukey test (p = 0.05) was performed to identify the least significant differences between treatment means for tomato yield, PFPN and AEN. Model checking was performed through basic diagnostic plots (normal Q-Q and residuals plots). Quadratic regression was performed to identify the relationships between various N rates and tomato yields. Prior, linear regres-sion was also tested to see which model fit the data well, and then quadratic regression was used based on the fitted line with a higher coefficient of determination (R 2 ) value over the linear model to explain the relationship between N fertilizer supply and yields. The N-response curve (N rate vs yield) was used to identify the right rate of N fertilizer to produce the highest yield across districts. A linear regression model was fitted to identify the relationship between various soil chemical parameters (pH, OM, total N, available P and K) and tomato yield across districts.

Tomato Yield
Both fertilizer treatments ( Figure 2) and location (Figure 3) had significant effects on the tomato yield. There was no interaction between the fertilizer practice and district locations.

Effects of Nutrient Omission
The tomato yield was significantly affected by nutrient omission across the districts (Figure 2a). To identify the limiting nutrient for yield, NPK omission was compared with the recommended NPK fertilizer (RP: N200 treatment) and a control (without NPK fertilizer: CK). On average, nitrogen omission (N0) reduced the tomato yield by 43% (41.5 ± 1.9 t ha −1 ) compared with RP (72.8 ± 3.5 t ha −1 ). P omission (P0) and K omission (K0) reduced the yield by 23% and 22% compared with RP. There was no significant difference between N0 and CK, while a significant difference was observed with P0 and K0 over CK, illustrating that N was the most limiting nutrient.

N source and Placement
The agronomic effect of UB deep placement was compared with conventional urea use and practices. N-200 (the recommended rate) and N-150 (the same rate as UB) were compared with UB deep placement applied at 150 kg N ha −1 (UBN-150). N-150, N-200 and UBN-150 increased the yield (p < 0.05) by 43%, 50% and 49% over CK and by 35%, 43 %, and 42% over N0 (Figure 2c). The application of UB and CU at 150 kg N ha −1 showed a significant effect on the tomato yield with an average increase of 12% (71.9 ± 3.4 t ha −1 ) over CU (63.6 ± 3.6 t ha −1 ) (Figure 2c).

PFPN and AEN
Significant differences were observed among the N fertilizer treatments (N-100, N-150, N-200, N-250 and UBN-150) in the PFPN and AEN ( Figure 5). PFPN and AEN were calculated to identify the efficiency of the added fertilizer in raising tomato yields.    Figure 5. Average partial factor productivity of N (PFPN (a)) and agronomic N use efficiency (AEN (b)) of various N rates and sources across districts; mean ± SE. Different letters above an error bar represent significant differences between different N treatments (post hoc Tukey test, p < 0.05).

Discussion
Nitrogen fertilizer strongly influenced the tomato yields across the district sites ( Figure 2). The omission of N fertilizer (N0) reduced the yield by twice as much as omission of P (P0) and K (K0) (Figure 2a), illustrating that N was the most limiting plant nutrient for tomato growth and development in those study districts, corroborating the results of previous studies in Nepal [13,23]. The fact that we found no significant differences between N0 (receiving a full dose of P and K) and CK indicates that there were no significant effects of P and K addition on the tomato yield. However, when compared with full NPK (recommended rate), on an average, there were significant differences in the tomato yield with P0 and K0 (Figure 2a), possibly due to synergistic effects among N, P and K [26,27]. Positive interactions (synergistic) among applied fertilizers increase nutrient use efficiency and yields [27]. The results from the NPK omission plot indicate that farmers should focus on optimizing N fertilizer management and adopt a balanced fertilization to have a synergistic effect and obtain higher tomato yields across the study sites.
After identifying N as the most important yield-limiting nutrient, the next step is to find strategies to increase NUE and achieve a better yield. Optimizing N fertilizer is an integral contributor to crop yield [15,25]. Our study identified the optimal N rate based on the yield response of tomato, which varied across the districts (Figure 4). Tomato yield increased with increasing N rate up to 200 kg N ha −1 [28] and began to decline beyond that rate (Figure 4a). When the rate was increased from 200 to 250 kg N ha −1 , the yield was reduced from 11 to 16% in Palpa, Dang and Surkhet, whereas negligible yield increases (parabolic) were observed in Doti and Kavre with an average of 2.5% (Figure 4b-f). Earlier studies have reported parabolic quadratic responses to different N rates, but the optimum rate varied in different studies based on different agro-ecological domains [11,12,14,15,29]. Li et al. [12] reported 300 kg N ha −1 to be the optimum for tomato yield (110 t ha −1 ), beyond which there was no yield increase and the yield declined when N rates were raised to 450 kg ha −1 . Similarly, another study [30] showed optimum tomato yields at 180 kg N ha −1 , beyond which the yield increment was negligible up to 300 kg N ha −1 . In our study, the yield increased by 43% at the optimum N rate of 200 kg ha −1 (Figure 4a). This is in close agreement with the meta-analysis conducted by Cheng et al. [25], where the optimal N rate ranged from 230 to 345 kg ha −1 , with an average increase of 59% in tomato yields (Figure 4a). In addition to its agronomic benefits, the identification of an optimum N rate is of high economic value in terms of reducing fertilizer costs and increasing profits, as well as environmental benefits in the form of greater NUE and reduced N leaching and nitrous oxide emissions [7,30].
Improved N application methods also increase NUE and help to better meet crop nutrient demand [28,31]. In our study, single UB applications of 150 kg N ha −1 with deep placement increased yields by 12% over conventional urea applied at the same N rate. UB with 25% less N produced similar tomato yields as conventional urea applied at the recommended rate (200 kg N ha −1 ). Few studies have assessed the agronomic effect of UB on vegetables, and markedly not on tomato. Thus, our results are discussed in reference to the agronomic effect observed in other vegetable and crops. In line with our results, Akter et al. [31] reported a 7.4% increase in bitter gourd yields with use of UB instead of the recommended conventional urea. In our previous studies, the application of UB increased cauliflower yield by 24% over the recommended CU [24]. The benefits of UB for cereal crops such as maize and rice are well documented [20,32,33].
The characteristics of the experimental sites, including the varied soil characteristics and biophysical factors such as temperature and rainfall, exerted significant and varied effects on the tomato yields, which were higher in Dang, Surkhet and Palpa and lower in Doti and Kavre ( Figure 3) [34][35][36]. However, in our study, moderate variation in soil characteristics across districts (Table 1) did not exert a significant effect on tomato yield ( Figure S3). Across the districts, soil pH, OC, total N, and available P and K were within the optimum range required for better tomato growth and development. Temperature variations might have affected yields, as temperature significantly influences leaf growth, truss appearance, pollination and fruiting [36]. Higher temperatures also favor early growth and enhanced photosynthesis under regularly irrigated conditions [35,36]. Higher tomato yields in Dang and Surkhet could thus be due to the higher temperatures at the experimental sites than at those in Doti and Kavre ( Figure S2). In Palpa, higher yields may owe to a greater amount of available P (Table 1), which plays a significant role in root and shoot growth and fruit development [37].
PFPN and AEN are important indices that help farmers and commercial growers to gauge the efficiency of added N fertilizer and its contributions to marketable yield [14]. PFPN decreased drastically with an increasing N rate (Figure 5a). AEN was found consistent while adding N from 100 kg ha −1 to 200 kg ha −1 , after which it decreased drastically with further N up to 250 kg ha −1 (Figure 5b). The decreased PFPN and AEN with excessive N fertilizer was in line with earlier studies [12,14]. UBN-150 increased PFPN by 8% and 21% and AEN by 22% and 27% over conventional urea applied at same rate (N-150) and the recommended N practice (N-200) ( Figure 5). In accordance with this, previous studies have reported significant benefits of UB with deep placement in increasing NUE [13,38] and reducing N losses to the environment [6,20].
Considering the associated tomato yield increases (Figure 2b,c and Figure 4a) and AEN ( Figure 5), the application of 200 kg N ha −1 in the form of CU and UB with 25% less N input (UBN-150) seems promising for achieving optimum tomato yields. Based on the yield response across the districts (Figure 4b-f), N could be applied at a rate of 150 kg ha −1 in Palpa and Kavre, and 200 kg ha −1 in Dang, Surkhet and Doti. The application of N fertilizer beyond this rate would not increase yield, and would reduce NUE and increase N losses through leaching and nitrous oxide emissions [6,19].
The average N use by farmers across the districts was 70 kg ha −1 (ranging from 27 to 286 kg ha −1 ) with an average tomato yield of 37 t ha −1 (ranging from 16 to 64 t ha −1 ) (Figure 1). In our result, N applied at a rate of 200 kg N ha −1 (recommended N) following improved management practices and UB with deep placement produced an average yield of 72.4 t ha −1 (Figure 2c), an increase of around 35.4 t ha −1 . This indicates that the use of the optimum N rate and an efficient N source with improved technology along with balanced fertilization (P and K addition) can increase the tomato yield by 49% over farmers' current fertilizer management practices. Moreover, considering the cost of urea (0.18 US$ kg −1 ) and UB (0.22 US$ kg −1 ) and keeping other input and management (seed and labor) costs constant, the use of 200 kg N in the form of CU and 150 kg N in the form of UB can provide a savings of around US $15,500 per hectare over farmers' current nutrient management practices.

Conclusions
Our study suggests that nitrogen is the most important factor among the primary nutrients affecting tomato yields. Optimizing N fertilization appears to be crucial for increasing yield, NUE and tomato farmers' profits. However, optimum N rates can vary across soil types and agro-ecologies and ranged from 150 to 200 kg N ha −1 . An N rate beyond 200 kg N ha −1 decreased yields and NUE. This study is representative of a subtropical climate; thus, future studies should be focused in tropical, temperate and other agroecological zones to recommend N fertilization across Nepal. Moreover, the right timing of N fertilization can substantially increase NUE and yields and reduce environmental pollution, thus future studies should be focused on investigating the right time for N fertilization to achieve higher yields and a healthy ecosystem.
The use of UB allowed the use of up to 25% less N (150 kg N ha −1 ) than the conventional urea dose of 200 kg N ha −1 without reducing yields. UB showed the highest AEN among the various N treatments, illustrating that deep placement or root zone applications are efficient in increasing yields and reducing environmental costs. Our results suggest the advisability of government policy support for large-scale production and distribution of UB in Nepal and the need to scale out the optimization of N fertilization and UB with deep placement at the farm level for sustainable agronomic and economic development in tomato farming.

Supplementary Materials:
The following supporting information can be downloaded at: https://www. mdpi.com/article/10.3390/soilsystems6030072/s1, Figure S1: Field trial locations across five districts in Nepal (Source: Department of Survey, Nepal); Figure S2: Annual minimum, maximum temperatures, and precipitation during growing period of tomato across five districts in the year 2018. Temperature and rainfall data were extracted from National Oceanic and Atmospheric Administration (NOAA; https://psl.noaa.gov) and rainfall estimates from Rain Gauze and Satellite Observations (CHIRPS; ftp://ftp.chg.ucsb.edu/pub/ org/chg/products/CHIRPS-2.0) respectively; Figure S3: Relationship between soil chemical properties and yield to identify indigenous soil nutrient supply across the districts (n = 28).