Productivity and Soil Health Evaluation of Cropping Systems with Hey Won NPK Fertilizers: A Predictive Modeling and Performance Analysis Approach Through Python

R
Rajeev Sikka1,*
S
Sat Pal Sharma1
N
Navjot Singh Brar2
R
Radha Ahuja2
S
Shephali Sachan2
M
Mandeep Kaur2
A
Annu Singh2
1Department of Vegetable Science, Punjab Agricultural University, Ludhiana-141 004, Punjab, India.
2Department of Soil Science, Punjab Agricultural University, Ludhiana-141 004, Punjab, India.

Background: With the intensification of agriculture, optimizing fertilizer use is critical for enhancing crop yields while maintaining soil health. This study evaluates the impact of three NPK fertilizer grades-Hey Won NPK 14-14-14, 14-5-21 and 20-8-12-on the productivity of maize-based cropping systems and their influence on soil properties.

Methods: A two-year field experiment (2022-2024) was conducted at the research farm of the Department of Soil Science, Punjab Agricultural University (PAU), Ludhiana, India, using a split-plot design. Four maize-based cropping systems were assigned to the main plots and five fertilizer treatments to the sub-plots. Key parameters such as grain yield, plant nutrient levels (N, P, K) and soil nutrient status (pre- and post-fertilization) were assessed.

Result: Preliminary results demonstrated significant differences in crop yields, system productivity and equivalent yields due to varying fertilizer treatments. The findings highlight the benefits of targeted NPK applications in boosting crop performance and improving soil nutrient profiles, supporting sustainable agriculture. Additionally, predictive modeling and performance analysis using Python identified key nutrient-yield relationships and the pivotal role of specific crops in enhancing nutrient-use efficiency. The study underscores the value of data-driven strategies for optimizing fertilizer management, offering actionable insights for farmers, agronomists and policymakers to improve productivity while preserving soil health across diverse cropping systems.

The rapidly growing population across the globe is estimated to reach over 9 billion by 2050, affecting sustainable agriculture development and food security. According to United Nations Food and Agricultural Organization, 99.7% of human food is fulfilled by terrestrial environment (Singh and Ryan, 2015). Therefore, it is essential to increase global agricultural production by 60%-110% for the coming 20 years (Grafton et al., 2015). Currently, the overall cropland area covers 11% out of 13 billion ha total land area which delivers about 78% of average per capita calorie consumption (Brevik, 2013).
       
Hence, increasing crop yields is the most suitable way rather than extending arable lands, for ensuring food security. In addition, utilizing balanced and complete fertilizers provides more than 50% increase in crop yields. Thus, the use of fertilizers has been become very evident (Selassie, 2020). There is a well-established significant relationship between food grain production and fertilizer consumption (Sunil et al., 2014). The cereals account for approximately 50.8% of the total world NPK requirement (Heffer, 2013).  The present fertilizer strategy is to increase agricultural production through its efficient and balanced use (Nadeem et al., 2022).
       
The nutrients play an important role in crop growth and development. Crop nutrients are elements which are essential in providing healthy and vigorous plants ultimately leading to higher yields. The secret to good nutrient management is to ensure that crops are receiving the right quantity of nutrients exactly when they need them most for sustainable growth and productivity of different cropping systems (Toor, 2021). As nutrients can be supplied by the soil, they should also be supplemented with organic manures and fertilizers. The nutrients that are required by crops in the largest amounts are N, P and K (Zewide and Melash, 2021). The use of fertilizers efficiently is an important strategy for increasing crop production. The organic matter application significantly and strongly affects the soil health even in small concentration. This keeps the optimum soil microbial population and also balances the nutrient cycling (Singh and Ryan, 2015). In the present study, Hey Won NPK fertilizers have been introduced by Taiwan Fertilizer Company (TFC) with the addition of humic acid, to explore the beneficial effects on crop production and soil health.
       
Complete NPK fertilizer out of the total nutrient composition supply is a major constrain in the country. However, there is limited information available on the effects of complex fertilizers on different crops. The present study investigated the effects of various NPK-grade (Hey Won) fertilizers on maize based cropping system and soil health along with identifying boosting cropping system yield through following Python based prediction and performance evaluation approach.
Site description
 
A field experiment was conducted for two years (2022-2024) at the research farm (Latitude- 30o56′ N, Longitude- 75o52′ E, Climate- semi-arid subtropical and Annual Rainfall- 700-800 mm), Department of Soil Science, Punjab Agricultural University (PAU), Ludhiana, India.
 
Experimental plan
 
The experiment followed a split plot design (plot size: 27 m²) with four maize-based cropping systems as main plots: Maize-potato-summer moong (CS1), Maize-cabbage-spring potato (CS2), Maize-potato-onion (CS3) and Maize-peas-eggplant (CS4). Sub-plots included five NPK treatments: control (T1), 14-14-14 (T2), 14-5-21 (T3), 20-8-12 (T4) and PAU conventional fertilizer (T5). Each treatment was replicated thrice in a fixed layout, totaling 15 treatments (Fig 1). Irrigation was applied as per crop and weather conditions. Baseline soil physico-chemical properties like soil texture (sandy-loamy), pH (7.14), electrical conductivity (0.147 dS/m), organic carbon (0.65%), nitrogen (150.5 kg/ha), phosphorous (29.6 kg/ha) and potassium (161.25 kg/ha) were measured.

Fig 1: Split plot experimental layout.


 
Material and crop description
 
Taiwan Fertilizer Company supplied Hey Won NPK fertilizers, transformed nitrophosphate blends enriched with humic acid to enhance soil properties and nutrient uptake (Bhatt et al., 2022). The tested grades included 14-14-14 (T2), 14-5-21 (T3) and 20-8-12 (T4), each supplying nitrogen (ammonical and nitrate), available phosphorus (P2O5) and water-soluble potassium (K2O). The PAU conventional fertilizer (T5) used urea (N source), SSP - single super phosphate (P source) and MOP- muriate of potash (K source) based on crop-specific recommended doses (RDF). Control (T1) plots received no fertilizer. With maize (PMH2) other crops such as summer moong (SML 1827), onion (Punjab Naroya), peas (Punjab 89), cabbage (Indian), eggplant (Punjab Bharpoor), potato and spring potato (Kufri Pukhraj) were cultivated as cropping system. Details of each crop variety specific spacing and RDF are provided in Thind and Mahal (2021).
 
Crop yield and cropping system use efficiency
 
Grain yields (q/ha) from maize-based cropping systems were recorded over two years and converted to maize equivalent yield (MEY) using minimum support prices (Uddin et al., 2009). System productivity (SP) was calculated as MEY (kg/ha) divided by the system duration (days) and expressed as kg/ha/day (Tomar and Tiwari, 1990). Land use efficiency (LUE) was derived by dividing total crop duration by 365 and expressed as a percentage (Kumar et al., 2014).
 
Plant macronutrient analysis
 
At harvest, grains and fruits were separated, dried and 0.5 g samples were ground. Nitrogen was determined via the Kjeldahl method (Jackson, 1967) after digestion with selenium dioxide-sulfuric acid. Phosphorus and potassium were analyzed from nitric-perchloric acid (3:1) digests using colorimetric and flame photometric methods, respectively (Jackson, 1967). Crop NPK concentrations (%) were summed by system to represent total NPK content per cropping system.
 
Soil analysis
 
Surface soil was taken from 0-15 cm depth and sieved (2 mm) for physico-chemical analysis. Soil pH and EC (1:2 soil:water) were measured using a pH meter and conductivity meter (dS/m). Available nitrogen, organic carbon, phosphorous and potassium were determined using standard methods (Sikka et al., 2024).
 
Statistical analysis
 
ANOVA was performed (GenStat 10th ed.) to assess plant and soil data, with significance tested via LSD (p≤0.05) and Duncan’s post hoc test for pairwise yield comparisons. Polynomial modeling was applied to evaluate cropping system yield dynamics under Hey Won fertilizers. Correlation analysis (SPSS 30.0) identified key plant and soil factors influencing yields, with multicollinearity excluded before final regression modeling. Python streamlined feature selection, regression diagnostics and visualizations (actual vs. predicted and residual plots) to ensure model robustness.
       
Regression equations were developed for each system using significant (p≤0.05 and 0.01) variables. Model performance was assessed using R2, RMSE and MAE. All prediction modeling and evaluation were conducted in Google Colab using Python libraries (scikit-learn, Pandas, NumPy, Matplotlib, Seaborn). Fig 2 outlines the complete workflow from data preparation to model refinement. Future studies may apply advanced models (e.g., decision trees, neural networks) for improved predictions.

Fig 2: Prediction modeling and performance evaluation flowchart of maize-based cropping systems performed by PYTHON.

Crop and cropping system yield
 
Hey Won fertilizers (14:14:14, 14:5:21 and 20:8:12) significantly (p<0.05) improved yields in maize-based systems over two years, with the most notable increases in spring potato, potato, cabbage, eggplant and maize (Table 1). Maximum yields were crop-specific: potato, eggplant and onion under T2; spring potato and cabbage under T3; pea, maize and summer moong under T4. This variability reflects crop genetic differences and fertilizer-crop compatibility (Setiawati et al., 2020). Control (T1) consistently yielded the lowest, while RDF (T5) out performed control but underperformed compared to Hey Won treatments.

Table 1: Effect of 14:14:14, 14:5:21 and 20:8:12 application on average grain yield (q/ha) of different crops grown in maize-based cropping system.


       
Yields under T2, T3 and T4 surpassed T5, with eggplant (T2  and T3, +77.01%) and summer moong (T2, +1.27%) showing the highest and lowest percentage increases, respectively. T3 notably boosted yields in four crops, likely due to improved nutrient synergy, root development and humic acid-enhanced soil health (Kolage et al., 2018; Bhatt and Singh, 2022).
       
System-wise, M-C-Sp (570.33 q/ha) and M-P-O (466.42 q/ha) under T3 outperformed M-Pe-E (411.40 q/ha) and M-P-Sm (339.71 q/ha) under T2. The M-C-Sp system’s superior yield suggests enhanced nutrient and resource use efficiency. Regression analysis (R2 = 0.80-0.96) confirms strong model-data fit (Fig 3). This approach allows us to observe how changes in fertilizer types and application rates can interact with the specific components of each cropping system, leading to varying yield outcomes (Shi et al., 2021).

Fig 3: Line graph presenting different Hey Won fertilizer treatment application on maize-based cropping system yield (q/ha) with polynomial regression analysis.


       
Overall, Hey Won fertilizers improved crop and system yields, aligning with Chimonyo et al. (2019) and Abid et al. (2020), who emphasized integrated nutrient management benefits. Vegetable-based systems may also offer better economic returns depending on crop and conditions. Walder et al. (2023) and this study highlight higher yields under organic-based fertilizers versus pure chemical inputs. Intercropping (e.g., maize-potato) further enhances yield and water efficiency, reducing water use by 3-13% (Xie et al., 2021).
 
Equivalent yield and resource use efficiency
 
Table 2 and 3 show significant differences in MEY and SP across cropping systems. The maize-potato-onion (M-P-O) system recorded the highest MEY (324.28 q/ha) and SP (117.49 kg/ha/day), while maize-potato-summer moong (M-P-Sm) had the lowest MEY (225.60 q/ha) and SP (80.00 kg/ha/day), likely due to poor maize-summer moong compatibility (Xie et al., 2021). M-Pe-E and M-C-Sp outperformed M-P-Sm in MEY, with onion inclusion boosting nutrient uptake (Singh and Sikka, 2007).

Table 2: Effect of 14:14:14, 14:5:21 and 20:8:12 application on maize equivalent yield.



Table 3: Effect of 14:14:14, 14:5:21 and 20:8:12 application on system productivity.


       
Fertilizer treatments significantly affected MEY and SP, with T3 (14:5:21) consistently outperforming other treatments, while the control (T1) yielded the least (Xie et al., 2021). Also, the M-Pe-E system showed the highest LUE (79.18%), while M-C-Sp had the lowest (71.51%), Fig 4, suggesting longer land utilization but potentially lower productivity due to factors like soil nutrient depletion (Walder et al., 2023). Similar patterns were noted by Kumar et al. (2014) in jute-based and by Sammauria et al. (2020) in pearl-millet based cropping systems. Ali et al. (2021) described the maximum pearl millet EY and SP in cotton-summer pearl millet cropping system under 75% RDN- inorganic fertilizer + 25% RDN-FYM, an integrated fertilization approach.

Fig 4: MEY, SP and LUE levels on different maize-based cropping system.


       
Other studies confirm the benefits of intercropping: Maize-soybean systems improved yield, water use and land productivity (Liu et al., 2018; Xu et al., 2020). Islam et al. (2020) reported maximum system productivity, profitability, sustainable yield index, production efficiency and relative economic efficiency but least land use efficiency in Maize + Green gram(1:2)- green gram + maize (1:1)-tomato, cropping system indicating that the sustainable production might not linked with the efficiently utilization of land cover. Sikka et al. (2024) also reported significant yield and MEY improvements with Sardar amin granules and Bentonite Sulphur in maize-based systems. The highest production efficiency was observed in soybean-onion cropping system than soybean-potato under 100% RDN (Recommended Dose Nutrient) applied with FYM@ 5 t per ha and biofertilizer suggesting the important role of onion and nutrient fertilizer amendments in improving the cropping system (CS3) performance (Patil et al., 2024).
 
Plant macronutrient status
 
Nitrogen (N), phosphorus (P) and potassium (K) levels varied significantly under different cropping systems and Hey Won fertilizer treatments (Fig 5). The M-P-O system recorded the highest N (11.46%) and K (3.36%) levels, likely due to the synergistic interaction of maize, potato and onion (Xie et al., 2021). Conversely, N (6.05%) and K (1.16%) were lowest in M-C-Sp and M-P-Sm, respectively. P content peaked in M-C-Sp (2.25%) and was lowest in M-P-O (0.93%), suggesting cabbage with spring potato may enhance P uptake under certain conditions, but may limit N and K due to competition (Thummanatsakun and Yampracha, 2018).

Fig 5: (A) Nitrogen (B) Phosphorous (C) Potassium levels under different maize-based cropping system and Hey Won fertilizer treatments.


       
Macronutrients spiked under T2, slightly declined at T3, stabilized at T4 and dropped at T5, with T1 showing the lowest levels. Despite lower nutrient levels under T3, yield and productivity were highest, indicating that efficient nutrient use and agronomic practices also drive yield (Xie et al., 2021). Nitrogen was generally dominant across all systems and treatments. Notably, P exceeded K in M-C-Sp, reflecting the plant’s ability to manage limiting factors.
       
Previous studies (Kumar et al., 2014; Ghosh et al., 2020; Jiang et al., 2024) support these findings, highlighting the importance of balanced NPK fertilization for enhanced nutrient uptake and yield. Our highest yield under T3 and strong M-P-O system performance align with Luitel et al. (2024), who reported optimal onion bulb yield under medium organic matter and high K2O in maize-based systems.
 
Soil health status
 
The root-nutrient uptake association is very strong at 0-15 cm as crops roots are highly concentrated at this depth. No significant changes were observed in soil pH, EC and K after two cropping cycles, indicating soil stability and balanced K dynamics (Table 4). However, fertilizer treatments significantly influenced soil OC, N and P levels. T3 (14% N, 5% P2O5, 21% K2O) showed the highest OC (at par with T2) and N, along with elevated P. Fertilization likely boosted microbial activity, enhancing nutrient cycling and soil fertility (Dinca et al., 2022). Available P peaked under T4, while K remained stable across T2, T3 and T4, suggesting no nutrient imbalance. T5 (100% RDF) improved soil parameters over control but was less effective than Hey Won fertilizers. Cropping systems had no notable effect on soil properties.

Table 4: Effect of 14:14:14, 14:5:21 and 20:8:12 fertilizers on soil health of maize-based cropping system.


       
Organic amendments with fertilizers significantly improve soil health and productivity (Dhaliwal et al., 2019). The humic acid in Hey Won fertilizers appears beneficial for sustaining soil quality (Sikka et al., 2024). Similarly, integrating organic and inorganic sources enhances soil fertility and yield (Bangre et al., 2024). Singh et al. (2021) also reported increased OC and N under high NPK doses in maize-vegetable pea systems, boosting yield and profitability.
 
Correlation, regression analysis (Prediction modeling) and Performance evaluation
 
Pearson correlation analysis was conducted for four cropping systems (CS1-CS4) using key variables: yield, MEY, SP, plant NPK and soil OC, N and P (Tables 5-8). Yield correlated significantly (p≤0.05/0.01) with P.P and P.K in CS1; P.N in CS2; P.N, P.K and S.OC in CS3 and P.K in CS4. This highlights system-specific nutrient-yield dynamics requiring tailored nutrient management (Roy et al., 2006; Yousaf et al., 2016). Yield, MEY and SP showed near-perfect correlations (r≈1.000), confirming their alignment as performance indicators (Bahadur et al., 2024).

Table 5: Correlation matrix among several plant and soil variables for M-P-Sm (CS1) cropping system.



Table 6: Correlation matrix among several plant and soil variables for M-C-Sp (CS2) cropping system.



Table 7: Correlation matrix among several plant and soil variables for M-P-O (CS3) cropping system.



Table 8: Correlation matrix among several plant and soil variables for M-Pe-E (CS4) cropping system.


       
Post-multicollinearity removal (Hair et al., 2014), regression models were developed for each system using Python (Table 9) to predict yield, with independent predictors showing minimum residuals (Kutner et al., 2005). Model performance was assessed via R², RMSE and MAE. CS3 outperformed others (R² = 0.975, RMSE = 15.34, MAE = 12.45), followed by CS1, CS2 and CS4. A composite score, integrating these metrics, confirmed CS3 as the optimal system (score = 0.52) (Crookston et al., 2021).

Table 9: Python based regression and composite score analysis for maize-based cropping systems.


       
Plots of actual vs. predicted yield and residuals vs. predicted yield (Fig 6-9) validated CS3’s high explanatory power and prediction accuracy (Chatterjee and Hadi, 2015). This suggests CS3 as the best maize-based system for enhancing yield. These findings align with prior studies linking nutrient dynamics and yield prediction using correlation and regression tools (Chimonyo et al., 2019; Macholdt et al., 2020). Similar modeling and composite scoring approaches have been applied by Di Paola et al. (2023) and Murthy et al. (2022) in other cropping systems.

Fig 6: (A) Actual vs. Predicted (B) Residual vs. Predicted yield plots of CS1 maize cropping system.



Fig 7: (A) Actual vs. Predicted (B) Residual vs. Predicted yield plots of CS2 maize cropping system.



Fig 8: (A) Actual vs. Predicted (B) Residual vs. Predicted yield plots of CS3 maize cropping system.



Fig 9: (A) Actual vs. Predicted (B) Residual vs. Predicted yield plots of CS4 maize cropping system.

In a present study, the maize-potato-onion (M-P-O) system proved most productive, with Hey Won fertilizers-enriched with humic acid-outperforming standard recommended fertilizer doses (RDF) by enhancing nutrient uptake, yield and soil health. Yield improvements were notable across maize, potato, spring potato, summer moong and pea crops, with minimal adverse effects on soil properties.
       
Correlation and regression analyses identified plant N, P, K and soil organic carbon as key drivers of yield. The M-P-O system showed superior regression model performance (high R², low RMSE  and MAE), confirming its efficiency in maximizing productivity. Polynomial modeling further highlighted the advantage of crop combinations with complementary growth patterns. This study emphasizes the importance of data-driven strategies, such as regression modeling, for optimizing fertilizer management and improving cropping system sustainability. These findings provide valuable tools for farmers and researchers to enhance productivity while safeguarding soil health.
The present study was supported by Taiwan Company Fertilizers and Punjab Agricultural University, Ludhiana, India.
 
All authors declared that there is no conflict of interest.

  1. Abid, M., Batool, T., Siddique, G., Ali, S., Binyamin, R., Shahid, M.J., Rizwan, M., Alsahli, A.A., Alyemeni, M.N. (2020). Integrated nutrient management enhances soil quality and crop productivity in maize-based cropping system. Sustainability12(23): 10214. DOI: 10.3390/su122310214.

  2. Ali, S., Patel, A.M. and Sharma, S. (2021). Impact of cropping systems and resource conservation techniques on productivity and profitability of systems. Indian Journal of Agricultural Research. 55(2): 175-180. doi: 10.18805/ IJARe.A-5403.

  3. Bahadur, R., Kumar, N., Tripathi, U., Singh, A.K., Kumari, G., Alagesan, A., Senthilkumar, T., Malathi, G., Prabhavathi, J. (2024). Enhancing productivity and sustainability through integrated farming systems: A review. International Journal of Environment and Climate Change. 14(1): 719-725. DOI: 10.9734/ijecc/2024/v14i13887.

  4. Bangre, J., Sinha, N.K., Mohanty, M., Jayaraman, S., Sahoo, R.N., Dwivedi, A.K., Singh, D., Mishra, A., Kumar, S., Wanjari, R., Jha, P., Kumar, D., Mishra, R., Mandloi, S., Singh, U.K., Chaudhary, R., Reddy, K.S., Prabhakar, M., Singh, V.K., Rao, C.S. (2024). Long-term impact of inorganic fertilizers and farmyard manure on soil quality and productivity in subtropical Vertisols under a soybean-wheat cropping system. Land Degradation  and Development. 35(14): 4257-4270. DOI: 10.1002/ldr.5220.

  5. Bhatt, P., Singh, V.K. (2022). Effect of humic acid on soil properties and crop production-A review. The Indian Journal of Agricultural Sciences. 92(12): 1423-1430. doi: 10.56093/ ijas.v92i12.124948.

  6. Brevik, E.C. (2013). Soils and Human Health: An Overview. In: Brevik, E.C., Burgess, L.C., editors. Soils and Human Health. CRC Press, USA: Boca Raton, FL. p.29-56.

  7. Chatterjee, S., Hadi, A.S. (2015). Regression Analysis by Example. Wiley.

  8. Chimonyo, V.G.P., Snapp, S.S., Chikowo, R. (2019). Grain legumes increase yield stability in maize based cropping systems. Crop Science. 59: 1-14. DOI:10.2135/cropsci2018.09.0532.

  9. Crookston, B., Yost, M., Bowman, M., Veum, K. (2021). Relationships of on farm soil health scores with corn and soybean yield in the midwestern United States. Soil Science Society of America Journal. 86: 91-105. doi: 10.1002/saj2.20355.

  10. Dhaliwal, S.S., Naresh, R.K., Mandal, A., Walia, M.K., Gupta, R.K., Singh, R., Dhaliwal, M. K. (2019). Effect of manures and fertilizers on soil physical properties, build-up of macro and micronutrients and uptake in soil under different cropping systems: A review. Journal of Plant Nutrition. 42(20): 2873-2900. DOI: 10.1080/01904167.2019.1659337.

  11. Di Paola, A., Di Giuseppe, E., Gutierrez, A. P., Ponti, L., Pasqui, M. (2023). Climate stressors modulate interannual olive yield at province level in Italy: A composite index approach to support crop management. Journal of Agronomy and Crop Science. 209: 475-488. DOI: 10.1111/jac.12636.

  12. Dinca, L.C., Grenni, P., Onet, C., Onet, A. (2022). Fertilization and Soil Microbial Community: A Review. Applied Sciences. 12(3): 1198. doi: 10.3390/app12031198. 

  13. Ghosh, D., Brahmachari, K., Brestic, M., Ondrisik, P., Hossain, A., Skalicky, M., Sarkar, S., Moulick, D., Dinda, N.K., Das, A. et al. (2020). Integrated weed and nutrient management improve yield, nutrient uptake and economics of maize in the rice-maize cropping system of Eastern India. Agronomy. 10(12): 1906. DOI: 10.3390/agronomy10121906.

  14. Grafton, R.Q., Williams, J., Jiang, Q. (2015). Food and water gaps to 2050: Preliminary results from the global food and water system (GFWS) platform. Food Security. 7: 209-220.

  15. Hair, J.F., Black, W.C., Babin, B.J. anderson, R.E. (2014). Multivariate Data Analysis. Pearson Education.

  16. Heffer, P. (2013). Assessment of fertilizer use by crop at the global level, 2010-2010/11. International Fertilizer Industry Association (IFA), Paris, France.

  17. Islam, M., Nath, L.K., Samajdar, T. (2020). Sustainable diversification of maize (Zea mays L.)-legumes cropping systems for productivity, profitability and resource use efficiency in West Garo Hills of Meghalaya, India. Legume Research- An International Journal. 43(3): 427-431. doi: 10.18805/ LR-3970.

  18. Jackson, L.W.R. (1967). Effect of shade on leaf structure of deciduous tree species. Ecology. 48(3): 498-499. doi: 10. 2307/1932686.

  19. Jiang, M., Dong, C., Bian, W. et al. (2024). Effects of different fertilization practices on maize yield, soil nutrients, soil moisture and water use efficiency in northern China based on a meta-analysis. Science Reports. 14: 6480. doi: 10.1038/s41598-024-57031-z.

  20. Kolage, A.K., Thawal, D.W., Mane, R.B., Pawar, P.P. (2018). Production potential and economic analysis of maize based vegetable cropping sequences. J. Agric. Res. Technol. 43(1): 30-35.

  21. Kumar, M., Singh, S.R., Jha, S.K., A, S., Mazumdar, S.P., Singh, A., Kundu, D.K., Mahapatra, B.S. (2014). System productivity, profitability and resource use efficiency of jute (Corchorus olitorius) based cropping systems in the eastern Indo- gangetic plain. The Indian Journal of Agricultural Sciences. 84(2): 209-213.

  22. Kutner, M.H., Nachtsheim, C.J., Neter, J. and Li, W. (2005). Applied Linear Statistical Models. McGraw-Hill.

  23. Liu, X., Rahman, T., Song, C., Yang, F., Su, B., Cui, L., et al. (2018). Relationships among light distribution, radiation use efficiency and land equivalent ratio in maize-soybean strip inter- cropping. Field Crop. Res. 224: 91-101. doi: 10.1016/ j.fcr.2012.09.019.

  24. Luitel, B.P., Bhusal, Y., Bhandari, B.B. (2024). Plant characteristics and bulb yield of onion as influenced by farmyard manure and potassium application under maize-based cropping system at Dailekh. Agriculture Development Journal. 17(1): 78-90. DOI: 10.3126/adj.v17i1.67874.

  25. Macholdt, J., Styczen, M.E., Macdonald, A., Piepho, H.-P., Honermeier, B. (2020). Long-term analysis from a cropping system perspective: Yield stability, environmental adaptability and production risk of winter barley. European Journal of Agronomy. 117: 126056. DOI: 10.1016/j.eja.2020.126056.

  26. Murthy, C.S., Poddar, M.K., Choudhary, K.K., Pandey, V., Srikanth, P., Ramasubramanian, S., Senthil Kumar, G. (2022). Paddy crop insurance using satellite-based composite index of crop performance. Geomatics, Natural Hazards and Risk. 13(1): 310-336. doi: 10.1080/19475705.2021.2025155.

  27. Nadeem, M.Y., Zhang, J., Zhou, Y., Ahmad, S., Ding, Y., Li, G. (2022). Quantifying the impact of reduced nitrogen rates on grain yield and nitrogen use efficiency in the wheat and rice rotation  system of the Yangtze river region. Agronomy. 12(4): 920.

  28. Patil, S.S., Hiwale, S.D., Shinde, R.S., Yadav, K.R. and Patil, S.R. (2024). Effect of integrated nutrient management levels on soybean and its residual effect on production efficiency and economic efficiency under soybean-onion and soybean- potato cropping system. Indian Journal of Agricultural Research. 58(5): 816-819. doi: 10.18805/IJARe.A-6059.

  29. Roy, R.N., Finck, A., Blair, G.J., Tandon, H.L.S. (2006). Plant nutrition for food security: A guide for integrated nutrient management. FAO Fertilizer and Plant Nutrition Bulletin.

  30. Sammauria, R., Meena, O.P., Yadav, M.R., Gupta, K.C., Gupta, A.K. and Yadav, H.L. (2020). Crop diversification of predominant pearl millet based cropping system for higher productivity, resource-use efficiency and profitability in semi-arid Rajasthan. Legume Research-An International Journal. 43(4): 568-572. doi: 10.18805/LR-4077.

  31. Selassie, Y.G., Molla, E., Muhabie, D., Manaye, F., Dessie, D. (2020). Response of crops to fertilizer application in volcanic soils. Heliyon. 6(12).

  32. Setiawati, M., Prayoga, M., Stöber, S., Adinata, K., Simarmata, T. (2020). Performance of rice paddy varieties under various organic soil fertility strategies. Open Agriculture. 5(1): 509-515. doi: 10.1515/opag-2020-0050.

  33. Shi, Z., Liu, D., Liu, M., Hafeez, M.B., Wen, P., Wang, X., et al. (2021). Optimized fertilizer recommendation method for nitrate residue control in a wheat-maize double cropping system in dryland farming. Field Crops Research. 271: 108258. DOI: 10.1016/j.fcr.2021.108258. 

  34. Sikka, R., Sekhon, B.S., Simranpreet, K., Sahni, T., Radha, A., Puttaswamy, C. (2024). Effect of sardar amin granules and bentonite sulfur on the productivity of maize-based cropping system in a sandy-loam soil. Environment and Ecology. 42(1): 27-33. doi: 10.60151/envec/QIBL9811.

  35. Singh, A., Pathania, P., Sharma, T., Sharma, S. (2021). Effect of different crop sequences on soil nutrient status, nutrient uptake and crop yield in Western Himalayas of India. Indian Journal of Agricultural Research. 55(6): 733-738. doi: 10.18805/IJARe.A-5774.

  36. Singh, B., Ryan, J. (2015). Managing fertilizers to enhance soil health. International Fertilizer Industry Association, Paris, France. 1: 1-24.

  37. Singh, T., Sikka, R. (2007). Performance of basmati rice-based cropping systems for productivity, profitability and soil health. J. Res Punjab Agricult Univ. 44(3): 177-180.

  38. Sunil, C.N., Kumar, V.V.N., Sanilkumar, M.G. (2014). A new species of garnotia (poaceae) from Kerala, India. Journal of the Botanical Research Institute of Texas8(2): 517-521. https: //www.jstor.org/stable/26549398.

  39. Thind, S.K., Mahal, J.S. (2021). Package of practices for cultivation of vegetables. Additional Director of Communication for Punjab Agricultural University, Ludhiana, 39.

  40. Thummanatsakun, V., Sukunya, Y. (2018). Effects of interaction between nitrogen and potassium on the growth and yield of cassava. International Journal of Agricultural Technology. 14: 2137-2150.

  41. Tomar, S.S., Tiwari, A.S. (1990). Production potential and economics of different crop sequences. Ind J. Agron. 35(1-2): 30-35. 

  42. Toor, M.D., Adnan, M., Rehman, F.U., Tahir, R., Saeed, M.S., Khan, A.U. and Pareek, V. (2021). Nutrients and their importance in agriculture crop production; A review. Ind. J. Pure App. Biosci. 9(1): 1-6.

  43. Uddin, M.J., Quayyum, M.A., Salahuddin, K.M. (2009). Intercropping of hybrid maize with short duration vegetables at hill vall eys of Bandarban. Bangladesh J. Agricul Res. 34(1): 51-57. 

  44. Walder, F., Büchi, L., Wagg, C., Colombi, T., Banerjee, S.,  Hirte, J., Mayer, J., Six, J., Kelle, T., Charles, R., van der Heijden, M.G.A. (2023). Synergism between production and soil health through crop diversification, organic amendments and crop protection in wheat-based systems. Journal of Applied Ecology. 60: 2091-2104. DOI: 10.1111/1365- 2664.14484.

  45. Xie, J., Wang, L., Li, L., Anwar, S., Luo, Z., Zechariah, E., Kwami Fudjoe, S. (2021). Yield, economic benefit, soil water balance and water use efficiency of intercropped maize/ potato in responses to mulching practices on the semiarid loess plateau. Agriculture. 11(11): 1100. Doi: 10.3390/ agriculture11111100.

  46. Xie, K., Cakmak, I., Wang, S., Zhang, F., Guo, S. (2021). Synergistic and antagonistic interactions between potassium and magnesium in higher plants. The Crop Journal. 9(2): 249-256. doi: 10.1016/j.cj.2020.10.005.

  47. Xu, Z., Li, C., Zhang, C., Yu, Y., van der Werf, W., Zhang, F. (2020). Intercropping maize and soybean increases efficiency of land and fertilizer nitrogen use; A meta-analysis. Field Crops Research. 246: 107661. doi: 10.1016/j.fcr.2019. 107661.

  48. Yousaf, M., Li, X., Zhang, Z., Ren, T., Cong, R., Ata-Ul-Karim, S.T., Fahad, S., Shah, A.N., Lu, J. (2016). Nitrogen fertilizer management for enhancing crop productivity and nitrogen use efficiency in a rice-oilseed rape rotation system in China. Front. Plant Sci. 7: 1496. DOI: 10.3389/fpls.2016. 01496.

  49. Zewide, I., W., M. (2021). Review on macronutrient in agronomy crops. Journal of  Nutrition and Food Processing. 4(6): 1-7.  doi: 10.31579/2637-8914/062.

Productivity and Soil Health Evaluation of Cropping Systems with Hey Won NPK Fertilizers: A Predictive Modeling and Performance Analysis Approach Through Python

R
Rajeev Sikka1,*
S
Sat Pal Sharma1
N
Navjot Singh Brar2
R
Radha Ahuja2
S
Shephali Sachan2
M
Mandeep Kaur2
A
Annu Singh2
1Department of Vegetable Science, Punjab Agricultural University, Ludhiana-141 004, Punjab, India.
2Department of Soil Science, Punjab Agricultural University, Ludhiana-141 004, Punjab, India.

Background: With the intensification of agriculture, optimizing fertilizer use is critical for enhancing crop yields while maintaining soil health. This study evaluates the impact of three NPK fertilizer grades-Hey Won NPK 14-14-14, 14-5-21 and 20-8-12-on the productivity of maize-based cropping systems and their influence on soil properties.

Methods: A two-year field experiment (2022-2024) was conducted at the research farm of the Department of Soil Science, Punjab Agricultural University (PAU), Ludhiana, India, using a split-plot design. Four maize-based cropping systems were assigned to the main plots and five fertilizer treatments to the sub-plots. Key parameters such as grain yield, plant nutrient levels (N, P, K) and soil nutrient status (pre- and post-fertilization) were assessed.

Result: Preliminary results demonstrated significant differences in crop yields, system productivity and equivalent yields due to varying fertilizer treatments. The findings highlight the benefits of targeted NPK applications in boosting crop performance and improving soil nutrient profiles, supporting sustainable agriculture. Additionally, predictive modeling and performance analysis using Python identified key nutrient-yield relationships and the pivotal role of specific crops in enhancing nutrient-use efficiency. The study underscores the value of data-driven strategies for optimizing fertilizer management, offering actionable insights for farmers, agronomists and policymakers to improve productivity while preserving soil health across diverse cropping systems.

The rapidly growing population across the globe is estimated to reach over 9 billion by 2050, affecting sustainable agriculture development and food security. According to United Nations Food and Agricultural Organization, 99.7% of human food is fulfilled by terrestrial environment (Singh and Ryan, 2015). Therefore, it is essential to increase global agricultural production by 60%-110% for the coming 20 years (Grafton et al., 2015). Currently, the overall cropland area covers 11% out of 13 billion ha total land area which delivers about 78% of average per capita calorie consumption (Brevik, 2013).
       
Hence, increasing crop yields is the most suitable way rather than extending arable lands, for ensuring food security. In addition, utilizing balanced and complete fertilizers provides more than 50% increase in crop yields. Thus, the use of fertilizers has been become very evident (Selassie, 2020). There is a well-established significant relationship between food grain production and fertilizer consumption (Sunil et al., 2014). The cereals account for approximately 50.8% of the total world NPK requirement (Heffer, 2013).  The present fertilizer strategy is to increase agricultural production through its efficient and balanced use (Nadeem et al., 2022).
       
The nutrients play an important role in crop growth and development. Crop nutrients are elements which are essential in providing healthy and vigorous plants ultimately leading to higher yields. The secret to good nutrient management is to ensure that crops are receiving the right quantity of nutrients exactly when they need them most for sustainable growth and productivity of different cropping systems (Toor, 2021). As nutrients can be supplied by the soil, they should also be supplemented with organic manures and fertilizers. The nutrients that are required by crops in the largest amounts are N, P and K (Zewide and Melash, 2021). The use of fertilizers efficiently is an important strategy for increasing crop production. The organic matter application significantly and strongly affects the soil health even in small concentration. This keeps the optimum soil microbial population and also balances the nutrient cycling (Singh and Ryan, 2015). In the present study, Hey Won NPK fertilizers have been introduced by Taiwan Fertilizer Company (TFC) with the addition of humic acid, to explore the beneficial effects on crop production and soil health.
       
Complete NPK fertilizer out of the total nutrient composition supply is a major constrain in the country. However, there is limited information available on the effects of complex fertilizers on different crops. The present study investigated the effects of various NPK-grade (Hey Won) fertilizers on maize based cropping system and soil health along with identifying boosting cropping system yield through following Python based prediction and performance evaluation approach.
Site description
 
A field experiment was conducted for two years (2022-2024) at the research farm (Latitude- 30o56′ N, Longitude- 75o52′ E, Climate- semi-arid subtropical and Annual Rainfall- 700-800 mm), Department of Soil Science, Punjab Agricultural University (PAU), Ludhiana, India.
 
Experimental plan
 
The experiment followed a split plot design (plot size: 27 m²) with four maize-based cropping systems as main plots: Maize-potato-summer moong (CS1), Maize-cabbage-spring potato (CS2), Maize-potato-onion (CS3) and Maize-peas-eggplant (CS4). Sub-plots included five NPK treatments: control (T1), 14-14-14 (T2), 14-5-21 (T3), 20-8-12 (T4) and PAU conventional fertilizer (T5). Each treatment was replicated thrice in a fixed layout, totaling 15 treatments (Fig 1). Irrigation was applied as per crop and weather conditions. Baseline soil physico-chemical properties like soil texture (sandy-loamy), pH (7.14), electrical conductivity (0.147 dS/m), organic carbon (0.65%), nitrogen (150.5 kg/ha), phosphorous (29.6 kg/ha) and potassium (161.25 kg/ha) were measured.

Fig 1: Split plot experimental layout.


 
Material and crop description
 
Taiwan Fertilizer Company supplied Hey Won NPK fertilizers, transformed nitrophosphate blends enriched with humic acid to enhance soil properties and nutrient uptake (Bhatt et al., 2022). The tested grades included 14-14-14 (T2), 14-5-21 (T3) and 20-8-12 (T4), each supplying nitrogen (ammonical and nitrate), available phosphorus (P2O5) and water-soluble potassium (K2O). The PAU conventional fertilizer (T5) used urea (N source), SSP - single super phosphate (P source) and MOP- muriate of potash (K source) based on crop-specific recommended doses (RDF). Control (T1) plots received no fertilizer. With maize (PMH2) other crops such as summer moong (SML 1827), onion (Punjab Naroya), peas (Punjab 89), cabbage (Indian), eggplant (Punjab Bharpoor), potato and spring potato (Kufri Pukhraj) were cultivated as cropping system. Details of each crop variety specific spacing and RDF are provided in Thind and Mahal (2021).
 
Crop yield and cropping system use efficiency
 
Grain yields (q/ha) from maize-based cropping systems were recorded over two years and converted to maize equivalent yield (MEY) using minimum support prices (Uddin et al., 2009). System productivity (SP) was calculated as MEY (kg/ha) divided by the system duration (days) and expressed as kg/ha/day (Tomar and Tiwari, 1990). Land use efficiency (LUE) was derived by dividing total crop duration by 365 and expressed as a percentage (Kumar et al., 2014).
 
Plant macronutrient analysis
 
At harvest, grains and fruits were separated, dried and 0.5 g samples were ground. Nitrogen was determined via the Kjeldahl method (Jackson, 1967) after digestion with selenium dioxide-sulfuric acid. Phosphorus and potassium were analyzed from nitric-perchloric acid (3:1) digests using colorimetric and flame photometric methods, respectively (Jackson, 1967). Crop NPK concentrations (%) were summed by system to represent total NPK content per cropping system.
 
Soil analysis
 
Surface soil was taken from 0-15 cm depth and sieved (2 mm) for physico-chemical analysis. Soil pH and EC (1:2 soil:water) were measured using a pH meter and conductivity meter (dS/m). Available nitrogen, organic carbon, phosphorous and potassium were determined using standard methods (Sikka et al., 2024).
 
Statistical analysis
 
ANOVA was performed (GenStat 10th ed.) to assess plant and soil data, with significance tested via LSD (p≤0.05) and Duncan’s post hoc test for pairwise yield comparisons. Polynomial modeling was applied to evaluate cropping system yield dynamics under Hey Won fertilizers. Correlation analysis (SPSS 30.0) identified key plant and soil factors influencing yields, with multicollinearity excluded before final regression modeling. Python streamlined feature selection, regression diagnostics and visualizations (actual vs. predicted and residual plots) to ensure model robustness.
       
Regression equations were developed for each system using significant (p≤0.05 and 0.01) variables. Model performance was assessed using R2, RMSE and MAE. All prediction modeling and evaluation were conducted in Google Colab using Python libraries (scikit-learn, Pandas, NumPy, Matplotlib, Seaborn). Fig 2 outlines the complete workflow from data preparation to model refinement. Future studies may apply advanced models (e.g., decision trees, neural networks) for improved predictions.

Fig 2: Prediction modeling and performance evaluation flowchart of maize-based cropping systems performed by PYTHON.

Crop and cropping system yield
 
Hey Won fertilizers (14:14:14, 14:5:21 and 20:8:12) significantly (p<0.05) improved yields in maize-based systems over two years, with the most notable increases in spring potato, potato, cabbage, eggplant and maize (Table 1). Maximum yields were crop-specific: potato, eggplant and onion under T2; spring potato and cabbage under T3; pea, maize and summer moong under T4. This variability reflects crop genetic differences and fertilizer-crop compatibility (Setiawati et al., 2020). Control (T1) consistently yielded the lowest, while RDF (T5) out performed control but underperformed compared to Hey Won treatments.

Table 1: Effect of 14:14:14, 14:5:21 and 20:8:12 application on average grain yield (q/ha) of different crops grown in maize-based cropping system.


       
Yields under T2, T3 and T4 surpassed T5, with eggplant (T2  and T3, +77.01%) and summer moong (T2, +1.27%) showing the highest and lowest percentage increases, respectively. T3 notably boosted yields in four crops, likely due to improved nutrient synergy, root development and humic acid-enhanced soil health (Kolage et al., 2018; Bhatt and Singh, 2022).
       
System-wise, M-C-Sp (570.33 q/ha) and M-P-O (466.42 q/ha) under T3 outperformed M-Pe-E (411.40 q/ha) and M-P-Sm (339.71 q/ha) under T2. The M-C-Sp system’s superior yield suggests enhanced nutrient and resource use efficiency. Regression analysis (R2 = 0.80-0.96) confirms strong model-data fit (Fig 3). This approach allows us to observe how changes in fertilizer types and application rates can interact with the specific components of each cropping system, leading to varying yield outcomes (Shi et al., 2021).

Fig 3: Line graph presenting different Hey Won fertilizer treatment application on maize-based cropping system yield (q/ha) with polynomial regression analysis.


       
Overall, Hey Won fertilizers improved crop and system yields, aligning with Chimonyo et al. (2019) and Abid et al. (2020), who emphasized integrated nutrient management benefits. Vegetable-based systems may also offer better economic returns depending on crop and conditions. Walder et al. (2023) and this study highlight higher yields under organic-based fertilizers versus pure chemical inputs. Intercropping (e.g., maize-potato) further enhances yield and water efficiency, reducing water use by 3-13% (Xie et al., 2021).
 
Equivalent yield and resource use efficiency
 
Table 2 and 3 show significant differences in MEY and SP across cropping systems. The maize-potato-onion (M-P-O) system recorded the highest MEY (324.28 q/ha) and SP (117.49 kg/ha/day), while maize-potato-summer moong (M-P-Sm) had the lowest MEY (225.60 q/ha) and SP (80.00 kg/ha/day), likely due to poor maize-summer moong compatibility (Xie et al., 2021). M-Pe-E and M-C-Sp outperformed M-P-Sm in MEY, with onion inclusion boosting nutrient uptake (Singh and Sikka, 2007).

Table 2: Effect of 14:14:14, 14:5:21 and 20:8:12 application on maize equivalent yield.



Table 3: Effect of 14:14:14, 14:5:21 and 20:8:12 application on system productivity.


       
Fertilizer treatments significantly affected MEY and SP, with T3 (14:5:21) consistently outperforming other treatments, while the control (T1) yielded the least (Xie et al., 2021). Also, the M-Pe-E system showed the highest LUE (79.18%), while M-C-Sp had the lowest (71.51%), Fig 4, suggesting longer land utilization but potentially lower productivity due to factors like soil nutrient depletion (Walder et al., 2023). Similar patterns were noted by Kumar et al. (2014) in jute-based and by Sammauria et al. (2020) in pearl-millet based cropping systems. Ali et al. (2021) described the maximum pearl millet EY and SP in cotton-summer pearl millet cropping system under 75% RDN- inorganic fertilizer + 25% RDN-FYM, an integrated fertilization approach.

Fig 4: MEY, SP and LUE levels on different maize-based cropping system.


       
Other studies confirm the benefits of intercropping: Maize-soybean systems improved yield, water use and land productivity (Liu et al., 2018; Xu et al., 2020). Islam et al. (2020) reported maximum system productivity, profitability, sustainable yield index, production efficiency and relative economic efficiency but least land use efficiency in Maize + Green gram(1:2)- green gram + maize (1:1)-tomato, cropping system indicating that the sustainable production might not linked with the efficiently utilization of land cover. Sikka et al. (2024) also reported significant yield and MEY improvements with Sardar amin granules and Bentonite Sulphur in maize-based systems. The highest production efficiency was observed in soybean-onion cropping system than soybean-potato under 100% RDN (Recommended Dose Nutrient) applied with FYM@ 5 t per ha and biofertilizer suggesting the important role of onion and nutrient fertilizer amendments in improving the cropping system (CS3) performance (Patil et al., 2024).
 
Plant macronutrient status
 
Nitrogen (N), phosphorus (P) and potassium (K) levels varied significantly under different cropping systems and Hey Won fertilizer treatments (Fig 5). The M-P-O system recorded the highest N (11.46%) and K (3.36%) levels, likely due to the synergistic interaction of maize, potato and onion (Xie et al., 2021). Conversely, N (6.05%) and K (1.16%) were lowest in M-C-Sp and M-P-Sm, respectively. P content peaked in M-C-Sp (2.25%) and was lowest in M-P-O (0.93%), suggesting cabbage with spring potato may enhance P uptake under certain conditions, but may limit N and K due to competition (Thummanatsakun and Yampracha, 2018).

Fig 5: (A) Nitrogen (B) Phosphorous (C) Potassium levels under different maize-based cropping system and Hey Won fertilizer treatments.


       
Macronutrients spiked under T2, slightly declined at T3, stabilized at T4 and dropped at T5, with T1 showing the lowest levels. Despite lower nutrient levels under T3, yield and productivity were highest, indicating that efficient nutrient use and agronomic practices also drive yield (Xie et al., 2021). Nitrogen was generally dominant across all systems and treatments. Notably, P exceeded K in M-C-Sp, reflecting the plant’s ability to manage limiting factors.
       
Previous studies (Kumar et al., 2014; Ghosh et al., 2020; Jiang et al., 2024) support these findings, highlighting the importance of balanced NPK fertilization for enhanced nutrient uptake and yield. Our highest yield under T3 and strong M-P-O system performance align with Luitel et al. (2024), who reported optimal onion bulb yield under medium organic matter and high K2O in maize-based systems.
 
Soil health status
 
The root-nutrient uptake association is very strong at 0-15 cm as crops roots are highly concentrated at this depth. No significant changes were observed in soil pH, EC and K after two cropping cycles, indicating soil stability and balanced K dynamics (Table 4). However, fertilizer treatments significantly influenced soil OC, N and P levels. T3 (14% N, 5% P2O5, 21% K2O) showed the highest OC (at par with T2) and N, along with elevated P. Fertilization likely boosted microbial activity, enhancing nutrient cycling and soil fertility (Dinca et al., 2022). Available P peaked under T4, while K remained stable across T2, T3 and T4, suggesting no nutrient imbalance. T5 (100% RDF) improved soil parameters over control but was less effective than Hey Won fertilizers. Cropping systems had no notable effect on soil properties.

Table 4: Effect of 14:14:14, 14:5:21 and 20:8:12 fertilizers on soil health of maize-based cropping system.


       
Organic amendments with fertilizers significantly improve soil health and productivity (Dhaliwal et al., 2019). The humic acid in Hey Won fertilizers appears beneficial for sustaining soil quality (Sikka et al., 2024). Similarly, integrating organic and inorganic sources enhances soil fertility and yield (Bangre et al., 2024). Singh et al. (2021) also reported increased OC and N under high NPK doses in maize-vegetable pea systems, boosting yield and profitability.
 
Correlation, regression analysis (Prediction modeling) and Performance evaluation
 
Pearson correlation analysis was conducted for four cropping systems (CS1-CS4) using key variables: yield, MEY, SP, plant NPK and soil OC, N and P (Tables 5-8). Yield correlated significantly (p≤0.05/0.01) with P.P and P.K in CS1; P.N in CS2; P.N, P.K and S.OC in CS3 and P.K in CS4. This highlights system-specific nutrient-yield dynamics requiring tailored nutrient management (Roy et al., 2006; Yousaf et al., 2016). Yield, MEY and SP showed near-perfect correlations (r≈1.000), confirming their alignment as performance indicators (Bahadur et al., 2024).

Table 5: Correlation matrix among several plant and soil variables for M-P-Sm (CS1) cropping system.



Table 6: Correlation matrix among several plant and soil variables for M-C-Sp (CS2) cropping system.



Table 7: Correlation matrix among several plant and soil variables for M-P-O (CS3) cropping system.



Table 8: Correlation matrix among several plant and soil variables for M-Pe-E (CS4) cropping system.


       
Post-multicollinearity removal (Hair et al., 2014), regression models were developed for each system using Python (Table 9) to predict yield, with independent predictors showing minimum residuals (Kutner et al., 2005). Model performance was assessed via R², RMSE and MAE. CS3 outperformed others (R² = 0.975, RMSE = 15.34, MAE = 12.45), followed by CS1, CS2 and CS4. A composite score, integrating these metrics, confirmed CS3 as the optimal system (score = 0.52) (Crookston et al., 2021).

Table 9: Python based regression and composite score analysis for maize-based cropping systems.


       
Plots of actual vs. predicted yield and residuals vs. predicted yield (Fig 6-9) validated CS3’s high explanatory power and prediction accuracy (Chatterjee and Hadi, 2015). This suggests CS3 as the best maize-based system for enhancing yield. These findings align with prior studies linking nutrient dynamics and yield prediction using correlation and regression tools (Chimonyo et al., 2019; Macholdt et al., 2020). Similar modeling and composite scoring approaches have been applied by Di Paola et al. (2023) and Murthy et al. (2022) in other cropping systems.

Fig 6: (A) Actual vs. Predicted (B) Residual vs. Predicted yield plots of CS1 maize cropping system.



Fig 7: (A) Actual vs. Predicted (B) Residual vs. Predicted yield plots of CS2 maize cropping system.



Fig 8: (A) Actual vs. Predicted (B) Residual vs. Predicted yield plots of CS3 maize cropping system.



Fig 9: (A) Actual vs. Predicted (B) Residual vs. Predicted yield plots of CS4 maize cropping system.

In a present study, the maize-potato-onion (M-P-O) system proved most productive, with Hey Won fertilizers-enriched with humic acid-outperforming standard recommended fertilizer doses (RDF) by enhancing nutrient uptake, yield and soil health. Yield improvements were notable across maize, potato, spring potato, summer moong and pea crops, with minimal adverse effects on soil properties.
       
Correlation and regression analyses identified plant N, P, K and soil organic carbon as key drivers of yield. The M-P-O system showed superior regression model performance (high R², low RMSE  and MAE), confirming its efficiency in maximizing productivity. Polynomial modeling further highlighted the advantage of crop combinations with complementary growth patterns. This study emphasizes the importance of data-driven strategies, such as regression modeling, for optimizing fertilizer management and improving cropping system sustainability. These findings provide valuable tools for farmers and researchers to enhance productivity while safeguarding soil health.
The present study was supported by Taiwan Company Fertilizers and Punjab Agricultural University, Ludhiana, India.
 
All authors declared that there is no conflict of interest.

  1. Abid, M., Batool, T., Siddique, G., Ali, S., Binyamin, R., Shahid, M.J., Rizwan, M., Alsahli, A.A., Alyemeni, M.N. (2020). Integrated nutrient management enhances soil quality and crop productivity in maize-based cropping system. Sustainability12(23): 10214. DOI: 10.3390/su122310214.

  2. Ali, S., Patel, A.M. and Sharma, S. (2021). Impact of cropping systems and resource conservation techniques on productivity and profitability of systems. Indian Journal of Agricultural Research. 55(2): 175-180. doi: 10.18805/ IJARe.A-5403.

  3. Bahadur, R., Kumar, N., Tripathi, U., Singh, A.K., Kumari, G., Alagesan, A., Senthilkumar, T., Malathi, G., Prabhavathi, J. (2024). Enhancing productivity and sustainability through integrated farming systems: A review. International Journal of Environment and Climate Change. 14(1): 719-725. DOI: 10.9734/ijecc/2024/v14i13887.

  4. Bangre, J., Sinha, N.K., Mohanty, M., Jayaraman, S., Sahoo, R.N., Dwivedi, A.K., Singh, D., Mishra, A., Kumar, S., Wanjari, R., Jha, P., Kumar, D., Mishra, R., Mandloi, S., Singh, U.K., Chaudhary, R., Reddy, K.S., Prabhakar, M., Singh, V.K., Rao, C.S. (2024). Long-term impact of inorganic fertilizers and farmyard manure on soil quality and productivity in subtropical Vertisols under a soybean-wheat cropping system. Land Degradation  and Development. 35(14): 4257-4270. DOI: 10.1002/ldr.5220.

  5. Bhatt, P., Singh, V.K. (2022). Effect of humic acid on soil properties and crop production-A review. The Indian Journal of Agricultural Sciences. 92(12): 1423-1430. doi: 10.56093/ ijas.v92i12.124948.

  6. Brevik, E.C. (2013). Soils and Human Health: An Overview. In: Brevik, E.C., Burgess, L.C., editors. Soils and Human Health. CRC Press, USA: Boca Raton, FL. p.29-56.

  7. Chatterjee, S., Hadi, A.S. (2015). Regression Analysis by Example. Wiley.

  8. Chimonyo, V.G.P., Snapp, S.S., Chikowo, R. (2019). Grain legumes increase yield stability in maize based cropping systems. Crop Science. 59: 1-14. DOI:10.2135/cropsci2018.09.0532.

  9. Crookston, B., Yost, M., Bowman, M., Veum, K. (2021). Relationships of on farm soil health scores with corn and soybean yield in the midwestern United States. Soil Science Society of America Journal. 86: 91-105. doi: 10.1002/saj2.20355.

  10. Dhaliwal, S.S., Naresh, R.K., Mandal, A., Walia, M.K., Gupta, R.K., Singh, R., Dhaliwal, M. K. (2019). Effect of manures and fertilizers on soil physical properties, build-up of macro and micronutrients and uptake in soil under different cropping systems: A review. Journal of Plant Nutrition. 42(20): 2873-2900. DOI: 10.1080/01904167.2019.1659337.

  11. Di Paola, A., Di Giuseppe, E., Gutierrez, A. P., Ponti, L., Pasqui, M. (2023). Climate stressors modulate interannual olive yield at province level in Italy: A composite index approach to support crop management. Journal of Agronomy and Crop Science. 209: 475-488. DOI: 10.1111/jac.12636.

  12. Dinca, L.C., Grenni, P., Onet, C., Onet, A. (2022). Fertilization and Soil Microbial Community: A Review. Applied Sciences. 12(3): 1198. doi: 10.3390/app12031198. 

  13. Ghosh, D., Brahmachari, K., Brestic, M., Ondrisik, P., Hossain, A., Skalicky, M., Sarkar, S., Moulick, D., Dinda, N.K., Das, A. et al. (2020). Integrated weed and nutrient management improve yield, nutrient uptake and economics of maize in the rice-maize cropping system of Eastern India. Agronomy. 10(12): 1906. DOI: 10.3390/agronomy10121906.

  14. Grafton, R.Q., Williams, J., Jiang, Q. (2015). Food and water gaps to 2050: Preliminary results from the global food and water system (GFWS) platform. Food Security. 7: 209-220.

  15. Hair, J.F., Black, W.C., Babin, B.J. anderson, R.E. (2014). Multivariate Data Analysis. Pearson Education.

  16. Heffer, P. (2013). Assessment of fertilizer use by crop at the global level, 2010-2010/11. International Fertilizer Industry Association (IFA), Paris, France.

  17. Islam, M., Nath, L.K., Samajdar, T. (2020). Sustainable diversification of maize (Zea mays L.)-legumes cropping systems for productivity, profitability and resource use efficiency in West Garo Hills of Meghalaya, India. Legume Research- An International Journal. 43(3): 427-431. doi: 10.18805/ LR-3970.

  18. Jackson, L.W.R. (1967). Effect of shade on leaf structure of deciduous tree species. Ecology. 48(3): 498-499. doi: 10. 2307/1932686.

  19. Jiang, M., Dong, C., Bian, W. et al. (2024). Effects of different fertilization practices on maize yield, soil nutrients, soil moisture and water use efficiency in northern China based on a meta-analysis. Science Reports. 14: 6480. doi: 10.1038/s41598-024-57031-z.

  20. Kolage, A.K., Thawal, D.W., Mane, R.B., Pawar, P.P. (2018). Production potential and economic analysis of maize based vegetable cropping sequences. J. Agric. Res. Technol. 43(1): 30-35.

  21. Kumar, M., Singh, S.R., Jha, S.K., A, S., Mazumdar, S.P., Singh, A., Kundu, D.K., Mahapatra, B.S. (2014). System productivity, profitability and resource use efficiency of jute (Corchorus olitorius) based cropping systems in the eastern Indo- gangetic plain. The Indian Journal of Agricultural Sciences. 84(2): 209-213.

  22. Kutner, M.H., Nachtsheim, C.J., Neter, J. and Li, W. (2005). Applied Linear Statistical Models. McGraw-Hill.

  23. Liu, X., Rahman, T., Song, C., Yang, F., Su, B., Cui, L., et al. (2018). Relationships among light distribution, radiation use efficiency and land equivalent ratio in maize-soybean strip inter- cropping. Field Crop. Res. 224: 91-101. doi: 10.1016/ j.fcr.2012.09.019.

  24. Luitel, B.P., Bhusal, Y., Bhandari, B.B. (2024). Plant characteristics and bulb yield of onion as influenced by farmyard manure and potassium application under maize-based cropping system at Dailekh. Agriculture Development Journal. 17(1): 78-90. DOI: 10.3126/adj.v17i1.67874.

  25. Macholdt, J., Styczen, M.E., Macdonald, A., Piepho, H.-P., Honermeier, B. (2020). Long-term analysis from a cropping system perspective: Yield stability, environmental adaptability and production risk of winter barley. European Journal of Agronomy. 117: 126056. DOI: 10.1016/j.eja.2020.126056.

  26. Murthy, C.S., Poddar, M.K., Choudhary, K.K., Pandey, V., Srikanth, P., Ramasubramanian, S., Senthil Kumar, G. (2022). Paddy crop insurance using satellite-based composite index of crop performance. Geomatics, Natural Hazards and Risk. 13(1): 310-336. doi: 10.1080/19475705.2021.2025155.

  27. Nadeem, M.Y., Zhang, J., Zhou, Y., Ahmad, S., Ding, Y., Li, G. (2022). Quantifying the impact of reduced nitrogen rates on grain yield and nitrogen use efficiency in the wheat and rice rotation  system of the Yangtze river region. Agronomy. 12(4): 920.

  28. Patil, S.S., Hiwale, S.D., Shinde, R.S., Yadav, K.R. and Patil, S.R. (2024). Effect of integrated nutrient management levels on soybean and its residual effect on production efficiency and economic efficiency under soybean-onion and soybean- potato cropping system. Indian Journal of Agricultural Research. 58(5): 816-819. doi: 10.18805/IJARe.A-6059.

  29. Roy, R.N., Finck, A., Blair, G.J., Tandon, H.L.S. (2006). Plant nutrition for food security: A guide for integrated nutrient management. FAO Fertilizer and Plant Nutrition Bulletin.

  30. Sammauria, R., Meena, O.P., Yadav, M.R., Gupta, K.C., Gupta, A.K. and Yadav, H.L. (2020). Crop diversification of predominant pearl millet based cropping system for higher productivity, resource-use efficiency and profitability in semi-arid Rajasthan. Legume Research-An International Journal. 43(4): 568-572. doi: 10.18805/LR-4077.

  31. Selassie, Y.G., Molla, E., Muhabie, D., Manaye, F., Dessie, D. (2020). Response of crops to fertilizer application in volcanic soils. Heliyon. 6(12).

  32. Setiawati, M., Prayoga, M., Stöber, S., Adinata, K., Simarmata, T. (2020). Performance of rice paddy varieties under various organic soil fertility strategies. Open Agriculture. 5(1): 509-515. doi: 10.1515/opag-2020-0050.

  33. Shi, Z., Liu, D., Liu, M., Hafeez, M.B., Wen, P., Wang, X., et al. (2021). Optimized fertilizer recommendation method for nitrate residue control in a wheat-maize double cropping system in dryland farming. Field Crops Research. 271: 108258. DOI: 10.1016/j.fcr.2021.108258. 

  34. Sikka, R., Sekhon, B.S., Simranpreet, K., Sahni, T., Radha, A., Puttaswamy, C. (2024). Effect of sardar amin granules and bentonite sulfur on the productivity of maize-based cropping system in a sandy-loam soil. Environment and Ecology. 42(1): 27-33. doi: 10.60151/envec/QIBL9811.

  35. Singh, A., Pathania, P., Sharma, T., Sharma, S. (2021). Effect of different crop sequences on soil nutrient status, nutrient uptake and crop yield in Western Himalayas of India. Indian Journal of Agricultural Research. 55(6): 733-738. doi: 10.18805/IJARe.A-5774.

  36. Singh, B., Ryan, J. (2015). Managing fertilizers to enhance soil health. International Fertilizer Industry Association, Paris, France. 1: 1-24.

  37. Singh, T., Sikka, R. (2007). Performance of basmati rice-based cropping systems for productivity, profitability and soil health. J. Res Punjab Agricult Univ. 44(3): 177-180.

  38. Sunil, C.N., Kumar, V.V.N., Sanilkumar, M.G. (2014). A new species of garnotia (poaceae) from Kerala, India. Journal of the Botanical Research Institute of Texas8(2): 517-521. https: //www.jstor.org/stable/26549398.

  39. Thind, S.K., Mahal, J.S. (2021). Package of practices for cultivation of vegetables. Additional Director of Communication for Punjab Agricultural University, Ludhiana, 39.

  40. Thummanatsakun, V., Sukunya, Y. (2018). Effects of interaction between nitrogen and potassium on the growth and yield of cassava. International Journal of Agricultural Technology. 14: 2137-2150.

  41. Tomar, S.S., Tiwari, A.S. (1990). Production potential and economics of different crop sequences. Ind J. Agron. 35(1-2): 30-35. 

  42. Toor, M.D., Adnan, M., Rehman, F.U., Tahir, R., Saeed, M.S., Khan, A.U. and Pareek, V. (2021). Nutrients and their importance in agriculture crop production; A review. Ind. J. Pure App. Biosci. 9(1): 1-6.

  43. Uddin, M.J., Quayyum, M.A., Salahuddin, K.M. (2009). Intercropping of hybrid maize with short duration vegetables at hill vall eys of Bandarban. Bangladesh J. Agricul Res. 34(1): 51-57. 

  44. Walder, F., Büchi, L., Wagg, C., Colombi, T., Banerjee, S.,  Hirte, J., Mayer, J., Six, J., Kelle, T., Charles, R., van der Heijden, M.G.A. (2023). Synergism between production and soil health through crop diversification, organic amendments and crop protection in wheat-based systems. Journal of Applied Ecology. 60: 2091-2104. DOI: 10.1111/1365- 2664.14484.

  45. Xie, J., Wang, L., Li, L., Anwar, S., Luo, Z., Zechariah, E., Kwami Fudjoe, S. (2021). Yield, economic benefit, soil water balance and water use efficiency of intercropped maize/ potato in responses to mulching practices on the semiarid loess plateau. Agriculture. 11(11): 1100. Doi: 10.3390/ agriculture11111100.

  46. Xie, K., Cakmak, I., Wang, S., Zhang, F., Guo, S. (2021). Synergistic and antagonistic interactions between potassium and magnesium in higher plants. The Crop Journal. 9(2): 249-256. doi: 10.1016/j.cj.2020.10.005.

  47. Xu, Z., Li, C., Zhang, C., Yu, Y., van der Werf, W., Zhang, F. (2020). Intercropping maize and soybean increases efficiency of land and fertilizer nitrogen use; A meta-analysis. Field Crops Research. 246: 107661. doi: 10.1016/j.fcr.2019. 107661.

  48. Yousaf, M., Li, X., Zhang, Z., Ren, T., Cong, R., Ata-Ul-Karim, S.T., Fahad, S., Shah, A.N., Lu, J. (2016). Nitrogen fertilizer management for enhancing crop productivity and nitrogen use efficiency in a rice-oilseed rape rotation system in China. Front. Plant Sci. 7: 1496. DOI: 10.3389/fpls.2016. 01496.

  49. Zewide, I., W., M. (2021). Review on macronutrient in agronomy crops. Journal of  Nutrition and Food Processing. 4(6): 1-7.  doi: 10.31579/2637-8914/062.
In this Article
Published In
Agricultural Science Digest

Editorial Board

View all (0)