Spatial Distribution of Nutrient Loads Based on Mineral Fertilizers Applied to Crops: Case Study of the Lobo Basin in C ô te d’Ivoire (West Africa)

: Eutrophication in the Lobo watershed remains a major problem. The work carried out has focused on chemical and biological analyses in the lake or in its immediate environment: they did not sufﬁciently take into account the diffuse transfer of nutrients over the entire watershed. This study aims to assess the nutrient (N and P) loads in the Lobo watershed, an agricultural area, to understand the spatio-temporal impacts of land management practices on eutrophication. The methodology uses two steps: streamﬂow calibration and nutrient (N and P) estimation using the Soil and Water Assessment Tool (SWAT) watershed model. Thus, the nutrient inputs were estimated based on the levels of N and P in every kilogram of Nitrogen-phosphorus-Potassium (NPK) type fertilizers applied by farmers. The average quantities of N and P applied to the crops were 47.24 kg ha − 1 and 21.25 kg ha − 1 . Results show a good performance on ﬂow calibration as evidenced using evaluation criteria R 2 , Nash–Sutcliffe Efﬁciency (NSE), and Percent Bias (PBIAS) of 0.63, 0.62, and − 8.1, respectively. The yields of inorganic N and soluble P varied from 0 to 0.049 kg ha − 1 and from 0 to 0.31 kg ha − 1 . These results show that the crops’ inorganic nitrogen requirements were higher than the demands for soluble phosphorus. Simulations relating to the organic N transfer revealed values ranging from 0.2 to 5 kg ha − 1 , while the transport of organic phosphorus was estimated to vary from 0.3 to 1.3 kg ha − 1 .


Introduction
FAO statistics show that agriculture, hunting, fishing, and forestry supported 2.57 billion people at the start of the new millennium.This figure represents 42 percent of humanity, including those who work within these sectors and their families [1].Therefore, agriculture is essential for the economy in most developing countries and even in industrialized countries like the USA, where agricultural exports reached $177 billion in 2021 [1].Today, few countries are able to experience significant economic growth without agricultural development.
In Côte d'Ivoire, one estimates about 47.5% of the population lives in rural areas, and the agricultural sector occupies 30.7% of the Gross Domestic Product (GDP) [2].Thus, the State of Côte d'Ivoire has set up several Government agencies structures such as National Rural Development Support Agency (ANADER in French), the National Center for Agronomic Research of Cote d'Ivoire (CNRA in French), the Marketing of Food Products Office (OCPV in French) and the Development of Rice Cropping National Office (ONDR in French) to oversee and assist with agricultural management and development for sustainable agricultural growth and self-sufficiency of foods.Faced with the rising demand linked to population growth, the Government of Côte d'Ivoire uses these agencies to encourage traditional farmers and industrialists to use chemicals to increase their agricultural productivity [3].Industrial (or chemical) agricultural fertilizers have emerged as an important feature to increase agricultural production and compensate for the lack of nutrients in soils subjected to intensive cultivation throughout the year or during successive years.In order to obtain the desired benefit from these fertilizers, they are brought to the plant via the soil, depending on the quantity, quality, and time of addition required, according to its growth phases and needs [4].However, excessive and indiscriminate use of these chemical fertilizers is seen due to a lack of training on the sustainable use of chemicals and lax policy on environmental protection and quality standards.This irresponsible use has a major impact on the environment.The most recurrent phenomenon today in Côte d'Ivoire is the eutrophication of water bodies.
The last decades have been marked by the eutrophication issue; this phenomenon is a form of indigestion of aquatic ecosystems stuffed with excessive amounts of nitrogen and phosphorus.In the wake of many human activities (industrial, agricultural, or domestic), these nutrients, used in particular as fertilizers for crops, are indeed dumped into waterways and groundwater [5].Within lakes, this mass influx of nutrients results in the development of plants, such as macroalgae of the green algae type or microalgae of the phytoplankton type, which can be harmful or toxic [6].This all-out plant proliferation can, in particular, cause a decrease in the oxygen concentration in the water and changes in biodiversity, thus leading to a degraded ecological state with a change in the structure and functioning of the ecosystems concerned [6].
The advanced eutrophication of lakes leads to water quality deterioration.Today, this phenomenon has become a major environmental problem in urban areas and fertilized agricultural basins in developing countries around the world [7].Nutrient pollution accelerates the eutrophication process and affects lake ecosystems and human well-being by undermining ecosystem services [8].Therefore, the eutrophication phenomenon and nutrient transfer mobilize a large scientific community.In fact, Xu et al. [9] analyzed the effects of seasonal variation on water quality parameters related to eutrophication in Lake Yangzong.Lin et al. [10] studied the assessment and management of lake eutrophication in Lake Erhai in China.In the same direction, an evaluation of the long-term effectiveness of the management of the Lake Rauwbraken internal load made it possible to control its eutrophication [11].Suresh et al. [12] also analyzed the recent advancement in water quality indicators for eutrophication in global freshwater lakes.In the same vein, Cheng et al. (2024) [13] studied nutrient transport following water transfer in China.
Several studies published previously [14][15][16][17][18][19] have highlighted the Lobo water body's advanced eutrophication state.These studies were based on the analysis of biological and physicochemical parameters from point samples taken from the water body and did not allow an in-depth understanding of the eutrophication phenomenon, in particular, the nutrient transfer in the basin.To remedy this lack in this context, methodologies based on the combination of hydrological models and geographic information systems (GIS) are believed to be more effective.Hydrological models of watersheds are driven using mathematical formulas linked to physical processes and conceptual models.The basin hydrologic functioning includes an analogy and a concept.In the analogy, soils and groundwater are considered reservoirs whose rate of emptying depends on the filling.The physically based models are directly linked to physical processes [20,21].The principle of calculation is based on precise physical and empirical representations of natural processes.Empirical models [22] such as the USLE (Universal Soil Loss Equation) erosion model [23], derived from laboratory or field experiments.These models are global: this the case of the GR4 model which can represent the watershed in a single entity for global assessment [24].However, global-scale models are not suitable for studying hydrological processes of watersheds, because they do not take into account the spatial variability of phenomena such as runoff or erosion directly related to land use.Thus, semi-distributed models such as the Soil and Water Assessment Tool (SWAT), come in addition to overcome this lack by now, considering topography, land cover, soil and soil management practices in the basin.
This study aims to assess the nutrient loads from the amounts of fertilizers applied to different crops in the Lobo watershed in order to understand the spatio-temporal evolution of nutrient transport processes using the SWAT model.

Study Area
The Lobo watershed at Nibéhibé is located in west-central Côte d'Ivoire between 6 • 0 0 and 7 • 0 0 West longitude and 6 • 54 and 8 • 0 0 North latitude.It covers about 6442.66 km 2 in area (Figure 1).The plateaus about 200 to 400 m altitude compose the relief.This makes it a little bit contrast [25].The underlying geology is primarily composed of granites, while the peaks of precipitation are observed in the month of September with 275 mm and in May with 150 mm.This makes it possible to define two seasons: a dry season and a rainy season.The dry season is from November to February, while the rainy season is observed from March to October.The average annual rainfall recorded in the study area over the period 1979-2014 was 1335 mm.This rainfall value remains lower than that reported by Yao et al. [26], which was 1437.4 mm from 1943 to 2010 in the same basin using the GR4J model.The estimated annual recharge is 84 mm in 2019 to 66.4 mm in 2020.As for the average direct recharge, it is estimated at 44 mm and 57.3 mm, respectively, in 2018 and 2019, represents about 4% and 5% of the precipitation [27].
The study area's vegetation is composed of dense, moist, semi-deciduous forest and cleared mesophilic forest.Ferralitic soil, strongly or moderately desaturated, governs this area.The main activity is agriculture based on coffee, cocoa, rubber, oil palm, and food crops [28].The drinking water supply comes from the Lobo reservoir exploited by SODECI.
Appl.Sci.2023, 13, x FOR PEER REVIEW 3 of 22 of watersheds, because they do not take into account the spatial variability of phenomena such as runoff or erosion directly related to land use.Thus, semi-distributed models such as the Soil and Water Assessment Tool (SWAT), come in addition to overcome this lack by now, considering topography, land cover, soil and soil management practices in the basin.This study aims to assess the nutrient loads from the amounts of fertilizers applied to different crops in the Lobo watershed in order to understand the spatio-temporal evolution of nutrient transport processes using the SWAT model.

Study Area
The Lobo watershed at Nibéhibé is located in west-central Côte d'Ivoire between 6°0'0'' and 7°0'0'' West longitude and 6°54' and 8°0'0'' North latitude.It covers about 6442.66 km 2 in area (Figure 1).The plateaus about 200 to 400 m altitude compose the relief.This makes it a little bit contrast [25].The underlying geology is primarily composed of granites, while the peaks of precipitation are observed in the month of September with 275 mm and in May with 150 mm.This makes it possible to define two seasons: a dry season and a rainy season.The dry season is from November to February, while the rainy season is observed from March to October.The average annual rainfall recorded in the study area over the period 1979-2014 was 1335 mm.This rainfall value remains lower than that reported by Yao et al. [26], which was 1437.4 mm from 1943 to 2010 in the same basin using the GR4J model.The estimated annual recharge is 84 mm in 2019 to 66.4 mm in 2020.As for the average direct recharge, it is estimated at 44 mm and 57.3 mm, respectively, in 2018 and 2019, represents about 4% and 5% of the precipitation [27].
The study area's vegetation is composed of dense, moist, semi-deciduous forest and cleared mesophilic forest.Ferralitic soil, strongly or moderately desaturated, governs this area.The main activity is agriculture based on coffee, cocoa, rubber, oil palm, and food crops [28].The drinking water supply comes from the Lobo reservoir exploited by SO-DECI.

Data
Data used in this study were a Digital Elevation Model (DEM), land use map, soil data (soil physicochemical and biological properties, soil map), hydro-climate, and agronomic data.The Figure 2 shows essential input data for the SWAT model.

Data
Data used in this study were a Digital Elevation Model (DEM), land use map, soil data (soil physicochemical and biological properties, soil map), hydro-climate, and agronomic data.The Figure 2 shows essential input data for the SWAT model.

Digital Elevation Model (DEM)
The Digital Elevation Model (DEM) of 30 m × 30 m resolution, in the Northern Hemisphere projection UTM Zone 30, downloaded from the website https://search.earthdata.nasa.gov/search,accessed on 1 February 2018, has been used to extract the hydrographic network and delineate the watershed.

Soil Data
To run the SWAT model, soil physicochemical and biological properties (texture, available water content, hydraulic conductivity, bulk density, and organic carbon content of different soil layers) and soil map are required.These data are available in the FAO database [30] and have been completed using physicochemical properties cited above obtained after field sampling and in-laboratory analyses.The dominant soil of the basin is the sandy-clay-loam, about 38%, followed by the sandy-loam soil (32%) of the basin.

Hydro-Climate Data
The climate data at a daily time step were obtained from the SWAT model website http://globalweather.tamu.edu/,accessed on 1 February 2018.These data include minimum and maximum temperatures, precipitation, relative humidity, solar radiation, and wind speed.The use of these data is justified by the fact that in the watershed, there are few data available that can allow adequate hydrological modeling of the watershed.Some authors, such as Mbungu and Kahaigili [31], believe that global gridded climatological datasets (GGCDs) could be combined with hydrological modeling to understand hydrological processes in data-scarce environments.Therefore, for this study, eleven (11) climatic stations were used over the period 1979-2014.The monthly streamflow data were collected from the National Drinking Water Office (ONEP in French) for the period from 1981 to 1994 at the Nibéhibé hydrometric station.

Agronomic Data
The types of fertilizers, the doses applied to each type of crop, the types of crop, agricultural practices, and agricultural calendars on the watershed were used in this study.Agronomic data come from the processing of the National Center for Agronomic Research technical sheets.They are available online at https://cnra.ci/nos-fiches-techniques/,accessed on 1 January 2020.

Computer Software
The computer tools used in this study consist of QSWAT 1.9 in the QGIS 2.6.1 interface, SWATCUP 2019, ENVI 4.7, and the WGN program.QSWAT 1.9 made it possible to delimit the watershed into sub-watersheds to extract the hydrographic network and to simulate the water flow and nutrient transfer.SWATCUP 2019 was used to calibrate the water flow.ENVI 4.7 has been useful in the processing of satellite images for obtaining land use.As for the WGN program, it made it possible to calculate the statistical parameters from the climatic data.

Methods
The methodology based upon the SWAT (Soil and Water Assessment Tool) model was structured into two main stages: flow calibration and nutrient estimation.

•
SWAT Model description SWAT is a semi-distributed agro-hydrological model [32] with a daily time step.SWAT is used for evaluating water quality and quantity, agricultural practices management, surface and subsurface water management, and sediment, nutrient, and pesticide transfer.SWAT is used worldwide by a large scientific community.
This model is based on the water balance equation, which is as follows [32]: where SWt is the final soil water content (mm) of the day, SWo is the in-soil initial water content (mm), Rday is the daily rainfall (mm), IRR is irrigation volume added to the soil (mm), Qsurf is the surface runoff volume (mm), Ea is the actual evapotranspiration (mm), and Wseep is percolation loss from the soil profile into the shallow aquifer (mm).

• SWAT-CUP and SUFI-2 algorithm
Uncertainties may occur while setting up and running the model.Uncertainties come either from measurement errors or from the model itself [33][34][35].To be used efficiently, the model's performance must be tested.This study tested the correlation between observed data and simulated outputs statistically and graphically.Thus, Swat-cup (SWAT Calibration Uncertainty Program) was used to optimize the flow model.The program is the most widely used tool by the SWAT community.SWAT-CUP allows to perform calibration, uncertainty, and global sensitivity analysis automatically (Table 1).Several methods are integrated into the SWAT-CUP program.We have Sequential Uncertainty Fitting version 2 (SUFI-2) [36,37], Generalized Likelihood Uncertainty Estimation (GLUE) [38], Particle Swarm Optimization (PSO) [39], Parameter Solution (PARASOL) [40] and Markov Chain Monte Carlo (MCMC) [41].

SUFI-2
In SUFI 2, it is considered that the uncertainty in the simulations is observed in a uniform way.The sources of uncertainties are the driving variables, the conceptual model, parameters, and measured data.

GLUE
In this method, once the general probability has been defined, all the parameters are randomly sampled from the previous distribution.The parameters are thus grouped either into a behavioral set or into a non-behavioral set by comparing them to a given threshold probability.The parameters are then weighted according to their behavior.Finally, the uncertainty is predicted.

PSO
Here, the uncertainty prediction method is based on stochastic population optimization.The optimization is performed from a random sampling of parameters.

PARASOL
During the PARASOL method, a global optimization criterion (GOC) is first fixed.The method seeks to minimize the objective functions (OF) or GOC from the Shuffle Complex algorithm (SCE-UA).
MCMC MCMC proceeds with a random sampling, which adapts to the posterior distribution.
For this study, the SUFI 2 method was used because it permits concomitantly to perform global parameter sensitivity analysis, uncertainty analysis, and calibration [42].

Global sensitivity analysis
The analysis of global sensitivity begins with the definition of the objective function.Objective functions within SUFI2 include R 2 , Chi2, NS, and R 2 multiplied by the line regression coefficient b, bR 2 , and SSQR (sum of squared errors) coefficients.However, for this study, NS and R 2 coefficients were the objective functions used.NS is used to analyze the strength of the model predictions and while R 2 coefficient is related to the correlation between simulations and observed values.Then, the values of the parameters to be optimized, are chosen based on Equation (2) [43]: where b j is the jth parameter and m is the number of parameters to be estimated.
The parameter global sensitivity in each simulation is factored in using the sensitivity matrix J, of g (b j ) [43]: where b j is the parameter, g is global sensitivity, i is all possible combinations of two simulations, and j is the number of columns (number of parameters).
Latin Hypercube sampling is the sampling method used in this study.It makes it possible to calculate the global sensitivity by calculating a multiple regression system of the parameters considered [43]: where . σ s and σ m are considered to be the respective standard deviations of the simulated and observed data; µ s and µ m are the averages of the simulated and observed data.
The relative impact of each parameter on the flow is identified from a statistical test.This impact is noted bi.The global sensitivity is related to the variation of the objective function, which depends on the observed variation of the parameter considered when all parameters change.

•
Uncertainty analysis P-factor and R-factor are used to calculate the uncertainties at 2.5% (Xl) and 97.5% (Xu) percentiles of the cumulative distribution at each point.P-factor indicates the percentage of observed data bracketed using 95PPU.R-factor is the average thickness of the 95PPU band obtained from the standard deviation of the observed data [44].
with Q i obs observed flow, Q i sim simulated flow, Q obs average observed flow and Q sim average simulated flow.

Model Setup
The model implemented follows the methodology adopted by Koua et al. [60].The Lobo watershed in Nibéhibé covers an area of 6442.66 km 2 .The delimitation of the watershed and the extraction of the hydrographic network were possible thanks to the processing of the DEM under QSWAT.The D8 algorithm [61] integrated into QGIS makes it easy to extract watershed boundaries and river networks using QSWAT.This step is followed by the integration of meteorological data, land use, pedology, and their physico-chemical properties.Then, the hydrological response units (HRU) combining land use type, soil type, and slope in a sub-basin were calculated.HRUs are the basis for calculating the SWAT model.They highlight the differences in evapotranspiration and other hydrological conditions between different types of land cover.The multiple HRU method of the SWAT model made it possible to constitute 163 HRUs.After that, QSWAT established the various input data tables.This made it possible to integrate information relating to the watershed, pedology, hydrometeorology, the HRUs themselves, water resources and their use, watershed management, wetlands, septic tanks, different reservoirs on the watershed, and the watershed master file.At this stage, the hydrological simulation of the Lobo watershed can be done.

Streamflow Calibration Process
The flow calibration at a monthly time step was possible by using observed data at the Nibéhibé hydrometric station using the SUFI-2 method [37,38] in the SWAT-CUP program [62].For this study, ten (10) flow parameters have been used.These parameters were the CN2 expressing the humidity conditions II, the quantity of water available in the soil SOL_AWC, the ESCO factor which translates the compensation coefficient of the soil evaporation, RCHRG_DP which is the fraction of percolation water from the deep aquifer, the re-evaporation coefficient of groundwater GW_REVAP, ALPHA_BF (Baseflow alpha factor), the water depth threshold in the shallow aquifer depth GWQMN required for return flow to occur, the water depth threshold in the shallow aquifer REVAPMN for revaporation to occur, the groundwater delay time GW_DELAY, the CANMX (Maximum storage of canopy).Flow calibration was performed over the period 1981 to 1994.A total of 500 iterations were performed during the monthly water flow simulation process.Then, the Nash-Sutcliffe coefficient (NSE), the percentage of bias (PBIAS), and the R 2 determination coefficient were calculated in order to judge the robustness of the SWAT model on this watershed.

Nutrient Loads Estimation
While reservoir eutrophication processes are significantly influenced by nutrient contents in the water, there is no information on the reservoir water quality.Because of the lack of observation data on nitrogen and phosphorus, the estimation of nutrient inputs to the reservoir was made on the basis of the quantity of NPK fertilizers applied to crops in the watershed because these nutrients come from NPK chemical fertilizers [17,18].Data relating to fertilizers, such as the amount, duration, and frequency of application, were integrated into SWAT according to crops by sub-basin and HRU.The mineral fertilizers are of the xN-yP-zK formula, x, y, and z representing, respectively, the proportions in percentage (%) of N, P, and K (Potassium) in 1 kg of fertilizer used.If the fertilizer formula is known, the quantities of N, P, and K are determined (Equations ( 8) and ( 9)) [33]: P(kg) = %P × 1kg(of the fertilizer)/100 (9) For this study, only N and P quantities were estimated.
From the soil, the crops take nitrogen in the form of nitrates [33].The principle of simulating the transfer of nitrogen in SWAT takes into consideration the natural and artificial processes of its introduction into the soil.Figure 3 shows the nitrogen cycle in the soil.Mineral nitrogen from mineral fertilizers is applied to the soil, either in the form of ammoniacal nitrogen or in the form of nitrates.The ammoniacal form is unstable and volatile.Nitrate remains the most stable form in the soil.It is in this form that it is taken up by plants.The organic form of nitrogen contains three pools (fresh, stable and active)."Fresh" organic nitrogen is associated with crop residues and microbial biomass, while "active and stable" organic nitrogen is associated with humus.The two forms of nitrogen (organic and mineral) are readjusted every day at the URH scale according to the inputs and transformation processes that cause the nitrogen to change form.For a large quantity applied, the remainder not used by plants is transferred by runoff water, especially since it is very soluble in water [53].Nitrate transportation within the watershed is performed via surface runoff, lateral flow, or percolation.The amount of nitrate contained in runoff is assessed from the nitrate concentration in moving water.This concentration is a proportion of the nitrate mass in the runoff volume [33].It is calculated according to Equation (10): where Conc NO3,mobile is the concentration of nitrate in mobile water for a given soil layer (kg mm −1 ), NO3 ly expresses the quantity of nitrate (kg ha −1 ), w mobile is the amount of mobile water in the layer (mm), θ e is the porosity fraction on which the anion exclusion depends, and SAT ly is the water content when the soil layer is saturated (mm).Surface runoff transports nutrients to the top 10 mm above the soil.Equation ( 13) was used to calculate the amount of nitrate in runoff.

NO3
is the quantity of nitrate retained during surface runoff (kg ha −1 ), β represents the nitrate percolation coefficient, Conc , expresses the concentration of nitrate in the runoff water for the top 10 mm of soil (kg mm −1 ), Q is the runoff height of the given day (mm).
Surface runoff allows organic nitrogen, mostly attached to sediments, to transfer to the main channel in the watershed.The amount of organic nitrogen transported with the sediments to the watercourse is calculated using Equation ( 14) [63,64]: where orgN is the amount of organic nitrogen contained in surface runoff transported in kg ha −1 , ConC is the concentration in g t −1 of soil organic nitrogen in the first 10 mm above soil (g t −1 ), sed is the amount of sediment produced on the given day (metric tons), the area is the HRU area (ha), and ε is the nitrogen enrichment rate.The nitrogen enrichment rate is estimated according to Menzel [65].Thus, Equation The amount of mobile water in the soil layer is assumed to be equal to that lost by runoff, lateral flow, or percolation and can be estimated by Equations ( 11) and ( 12).w mobile = Q surf + Q lat,ly + w perc,ly for top 10 mm (11) w mobile = Q lat,ly + w perc,ly for lower soil layers, Q surf is the depth of surface water runoff on the given day (mm), Q lat,ly is the height of water escaping from the layer laterally (mm), w perc,ly is the amount of percolation water in the underlying soil layer on the given day (mm).
Surface runoff transports nutrients to the top 10 mm above the soil.Equation ( 13) was used to calculate the amount of nitrate in runoff.
NO3 surf is the quantity of nitrate retained during surface runoff (kg ha −1 ), β NO3 represents the nitrate percolation coefficient, Conc NO3,mobile expresses the concentration of nitrate in the runoff water for the top 10 mm of soil (kg mm −1 ), Q surf is the runoff height of the given day (mm).
Surface runoff allows organic nitrogen, mostly attached to sediments, to transfer to the main channel in the watershed.The amount of organic nitrogen transported with the sediments to the watercourse is calculated using Equation ( 14) [63,64]: where orgN surf is the amount of organic nitrogen contained in surface runoff transported in kg ha −1 , ConC orgN is the concentration in g t −1 of soil organic nitrogen in the first 10 mm above soil (g t −1 ), sed is the amount of sediment produced on the given day (metric tons), the area is the HRU area (ha), and ε Nsed is the nitrogen enrichment rate.
The nitrogen enrichment rate is estimated according to Menzel [65].Thus, Equation ( 15) was used to calculate the nitrogen enrichment rate: where Conc sed,surq is the in-runoff sediment concentration (t m −3 ).The concentration of sediment in surface runoff is calculated according to Equation ( 16): Conc sed,surq = sed 10.area hru .Q surf (16) where sed is the amount of sediment on a given day (metric tons), area hru is the HRU area (ha), and Q surf is the amount of runoff on a given day (mm).Like nitrogen, SWAT also divides the phosphorus terrestrial cycle into six different interconnected pools [33].These six pools are grouped according to their characteristics into two forms.Three phosphorus pools are inorganic form, and the other three are organic form."Fresh" organic phosphorus is associated with crop residues and microbial biomass, and mineral phosphorus is associated with soil humus.The mineralization of organic phosphorus only takes place from the active pool containing the humic substances and from the residues' pool.The mineral forms of phosphorus are further subdivided into three "pools": solution, active and stable.The equilibrium between the pool in solution and the active pool is established quickly, from a few days to a few weeks, whereas the dynamic between the active and stable pool evolves much more slowly (Figure 4).and the active pool is established quickly, from a few days to a few weeks, whereas the dynamic between the active and stable pool evolves much more slowly (Figure 4).The movement of soluble phosphorus in the soil is by diffusion.Due to its low mobility, transfer by surface runoff only concerns soluble phosphorus stored in the top 10 mm of soil [64].The amount of soluble P transferred is calculated from Equation ( 17): . The movement of soluble phosphorus in the soil is by diffusion.Due to its low mobility, transfer by surface runoff only concerns soluble phosphorus stored in the top 10 mm of soil [64].The amount of soluble P transferred is calculated from Equation ( 17): where P surf is the quantity of soluble phosphorus retained during runoff, expressed in kg ha −1 , P soluble,surf is the quantity of soluble phosphorus in solution in the first 10 mm (kg ha −1 ), Q surf is the quantity of surface water runoff on the day in question (mm), ρ b is the soil density over the first 10 mm expressed in t m −3 , depth surf is the depth of the surface layer (10 mm), and k d,surf is the partition coefficient of phosphorus in the soil (m 3 t −1 ).Surface runoff to the main channel allows the organic and mineral phosphorus attached to soil particles to transfer.All the nutrient parameters were integrated into the SWAT model, and nutrient flows were simulated running one time only the model with stream flow calibrated parameters.

Streamflow Parameter Global Sensitivity
The parameter global sensitivity analysis has shown that GW_REVAP, GWQMN, RCHRG_DP, ESCO, and CN2 are the most sensitive parameters to the flow (Tables 2 and 3).It is noted that the most sensitive among these parameters mentioned above are those related to the groundwater flow (GW_REVAP, GWQMN, and RCHRG_DP).This sensitivity analysis showed a possible connection between groundwater and surface water within the watershed.

Streamflow Calibration and Uncertainty
The analysis of the results showed that the SWAT model is suitable for flow evaluation in the Lobo watershed at Nibéhibé.Indeed, R 2 , NSE, and PBIAS had values of 0.63, 0.62, and −8.1, respectively.The values of P_factor and R_factor were 0.48 and 0.52 (Table 4).The graph in Figure 5 shows the results of the water flow calibration.SWAT model predictions generally follow the trend of observed flow.However, we note a shift in the flow peaks at the end of the simulation period (1981)(1982)(1983)(1984)(1985)(1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994).
From 1 January 1981 to 1 January 1990, the model has globally reproduced the observed flow correctly.However, from 1 January 1990 to 31 December 1994, the opposite is observed.Indeed, the SWAT model could not reproduce the peak flows and overestimated the base flow.
The uncertainties in the predictions are highlighted by P-factor and R-factor.The P-factor and R-factor values obtained were 0.48 and 0.52, respectively.This situation is explained by the fact that the SWAT model could not faithfully restore the peakflow from 1 January 1983 to 1 January 1985.The uncertainties in the predictions are highlighted by P-factor and R-factor.The Pfactor and R-factor values obtained were 0.48 and 0.52, respectively.This situation is explained by the fact that the SWAT model could not faithfully restore the peakflow from 1 January 1983 to 1 January 1985.

Nutrient Requirements for Crops
After integrating the different crops and fertilizers used in the Lobo watershed, the SWAT model estimated the average nutrient requirement (N and P) of the main crops grown in the watershed (Table 5).We note that the average needs observed by the CNRA is 41 kg ha −1 for nitrogen and 28 kg ha −1 for phosphorus.One notices an overestimation of N of +15% and an underestimation of P of −24% compared to CNRA results.

Nutrients Fluxes 3.3.1. Nutrient Requirements for Crops
After integrating the different crops and fertilizers used in the Lobo watershed, the SWAT model estimated the average nutrient requirement (N and P) of the main crops grown in the watershed (Table 5).We note that the average needs observed by the CNRA is 41 kg ha −1 for nitrogen and 28 kg ha −1 for phosphorus.One notices an overestimation of N of +15% and an underestimation of P of −24% compared to CNRA results.The mineral nitrogen and soluble phosphorus spatial distribution on the basis of fertilizers used by crop type is presented in Figure 6A,B.In Figure 6A, the variations in mineral nitrogen vary from 0 to 0.049 kg ha −1 .The highest mineral nitrogen fluxes are observed overall in the northern part of the watershed.The sub-basins concerned are 1, 2, 3, 4 and 5. Sub-basin 3 presents the largest contribution (0.049 kg ha −1 ) of the entire Lobo catchment area taken in Nibéhibé.The south and southwest parts also show high flows but a little less compared to the north part.These quantities of nitrogen per hectare vary from 0 to 0.037 kg ha −1 .The sub-basins with large quantities are sub-basins 28, 29 (containing the Lobo water reservoir), 30, and 31.It is generally observed that the mineral nitrogen inputs over the entire watershed are weak.Figure 6B shows the spatial variation in soluble phosphorus levels by sub-watershed.The soluble phosphorus flux varies from 0 to 0.31 kg ha −1 .The higher values are observed in the central, eastern, and western parts of the basin with rates of up to 0.31 kg ha −1 .It is also noted that the flows of mineral nitrogen are lower than those of phosphorus.This difference stems from the NPK fertilizer formula used.The sub-basins with high inputs are 1, 2, 3, 4, 5, 6, and 17, with rates varying between 0.8 and 1.3 kg ha −1 .The results of this study show that the contributions of nutrients of exogenous and organic origin remain higher than those of mineral origin.

Nitrates and Soluble Phosphorus Concentrations in Streams
The monthly mean nitrate concentrations in mg L −1 were also estimated by the SWAT model.The predictions showed values ranging from 14 to 30 mg L −1 (Figure 7A).The result showed that the left bank and downstream part of the Lobo River are more vulnerable to an excessive intake of nitrates.This excessive inflow is extended to the outlet of the

Organic Nitrogen and Organic Phosphorus Transferred per Sub-Basin
The simulation of organic nitrogen transport by sub-basin revealed inputs varying from 0.2 to 5 kg ha −1 (Figure 6C).One notes that the sub-basins in the southern part, except sub-basin 16, present a greater contribution of organic nitrogen and organic phosphorus than those of sub-basins located in the central and northern areas.The sub-basins with higher rates of organic nitrogen and organic phosphorus are 12, 13, 14, 24, and 29, with values reaching 3.7 to 5 kg ha −1 .As for organic phosphorus, it is shown in Figure 6D.Unlike the distribution of organic nitrogen, we observe higher amounts only in the northern part of the watershed.
The sub-basins with high inputs are 1, 2, 3, 4, 5, 6, and 17, with rates varying between 0.8 and 1.3 kg ha −1 .The results of this study show that the contributions of nutrients of exogenous and organic origin remain higher than those of mineral origin.

Nitrates and Soluble Phosphorus Concentrations in Streams
The monthly mean nitrate concentrations in mg L −1 were also estimated by the SWAT model.The predictions showed values ranging from 14 to 30 mg L −1 (Figure 7A).The result showed that the left bank and downstream part of the Lobo River are more vulnerable to an excessive intake of nitrates.This excessive inflow is extended to the outlet of the watershed (Nibéhibé) through sub-basin 29, where the Lobo water reservoir is located.Sub-basin 15 is the most vulnerable, with concentrations up to 30 mg L −1 .Figure 7B presents the monthly average concentrations of soluble phosphorus in mg L −1 .The values vary from 0.1 to 1.8 mg L −1 (Figure 7B).The rivers in sub-watersheds 30 and 31 have the highest soluble phosphorus concentrations (1.3-1.8 mg L −1 ).Sub-basins 21, 26, 28 and 29 also recorded high concentrations (0.9-1.3 mg L −1 ).One noted that both for nitrates and soluble phosphorus, the sub-basin 29 containing the Lobo water reservoir recorded a significant contribution.A sampling campaign carried out made it possible to measure nitrate and phosphorus concentrations in surface water in the watershed.The measurement points have been represented on the maps of Figure 7A,B.These measurements revealed nitrate concentrations ranging from 7.4 to 74 mg L −1 and phosphorus ranging from 2.3 to 66.6 mg L −1 .In the absence of historical data that can be used to calibrate nutrient inputs (N and P), these measurement points were superimposed on the spatial distribution map in order to verify the robustness of the SWAT model predictions.It appears that high concentration points (dark blue for nitrates and red for phosphorus) coincide with areas of streams where predictions show a high value (orange and red).The points of high concentrations may remain the same over the years.However, the amount of nitrates and phosphorus accumulated can increase significantly over the years.These measurements revealed nitrate concentrations ranging from 7.4 to 74 mg L −1 and phosphorus ranging from 2.3 to 66.6 mg L −1 .In the absence of historical data that can be used to calibrate nutrient inputs (N and P), these measurement points were superimposed on the spatial distribution map in order to verify the robustness of the SWAT model predictions.It appears that high concentration points (dark blue for nitrates and red for phosphorus) coincide with areas of streams where predictions show a high value (orange and red).The points of high concentrations may remain the same over the years.However, the amount of nitrates and phosphorus accumulated can increase significantly over the years.

Discussion
The observed flow data used for this study are available for the period 1981-1994.It is assumed that there is no major land use and rainfall trend over this period.In ungauged hydrological catchments where there is a scarcity of historical data, this hypothesis makes it possible to simulate flow and nutrient dynamics [66].The study revealed that flow parameters such as GW_REVAP, GWQMN, RCHRG_DP, ESCO, and CN2 impact the water flow the most.The same conclusion was made by Anoh et al. [67] in the Taabo watershed bordering that of Lobo.Along the same lines, Dakhlalla et al. [68] showed using a study on the evaluation of flow, nutrient, and sediment transfer parameters sensitivity and uncertainty that hydrological parameters such as GWQMN, ESCO, and CN2 are sensitive to flow.ESCO and CN2 are also parameters considered to significantly influence streamflow in western Mississippi [69].Moreover, CN2, being a parameter very sensitive to runoff, a decrease in precipitation impacts surface runoff, which will modify the values of CN2 [70][71][72] in [73].However, a parameter can be identified as not sensitive.However, this does not really mean that it does not influence the flow in the watershed [74].
The values of R 2 , NSE, and PBIAS of 0.63, 0.62, and −8.1, respectively, highlighted the capacity of the SWAT model according to the evaluation criteria proposed by Moriasi et al. [50].Upon the Taabo watershed, similar results were obtained by Anoh et al. [75].Indeed, the average NSE obtained by these authors was 0.7 during the calibration period.A PBIAS value of −8.1 obtained in this study showed a very good performance of the model [50].For flow, any P-factor > 0.7 and R-factor < 1.5 are considered acceptable.Therefore, the uncertainties on the flow predictions are considered acceptable [76].P-factor has a value of 0.48 which is therefore less good since the wish is for this value to be close to 1 [44].On the other hand, R-factor of 0.52 is acceptable.A model is said to be perfect when R-factor is close to 0. Very often, these values are very difficult to reach.This is why some authors recommend a balance between these two parameters [44].In our case, P-factor is 0.48 and R-factor 0.52.There is a certain balance between these parameters.The prediction uncertainties are therefore acceptable.Despite these shortcomings of the model, the results obtained in this study constitute a solid basis for hydrological modeling in the region with a view to the fight against the eutrophication of surface water facing the population of this region.Compared with the results of the [77], the predictions of SWAT show an underestimation of the N rates of −9.9% and of −5.55% for P rates.The results show low inorganic nitrogen levels.These low inorganic nitrogen levels in the Lobo watershed have also been highlighted by the work of [17,18].This could be explained by the fact that the recommendations of agricultural advisers to farmers regarding chemical fertilizers are respected.Then, nature, including crops, is able to digest the applied fertilizers, and only a small amount remains in the environment.The flows of mineral nitrogen are lower than those of phosphorus.This difference stems from the NPK fertilizer type used.Indeed, surveys carried out with ANADER (National Agency for Rural Development Support) and CIDT (Ivorian Textile Development Company) revealed that all agricultural activities practiced in the watershed use NPK fertilizer.In fact, NPK-10-18-18 fertilizer is used equally well for Coffee, Cocoa and cotton.However, coffee and cocoa are the main industrial crops in this region of Cote d'Ivoire.This fertilizer formula contains 18% phosphorus against 10% nitrogen in 1 kg of this fertilizer, hence the high soluble phosphorus levels relative to nitrogen.Leakage of nitrogen and phosphorus by the SWAT model takes place outside the root zone, i.e., primarily by infiltration for nitrates and runoff for phosphorus, although the model estimates that a fraction of nitrogen migrates by leaching in runoff [78].The variability of nitrogen and phosphorus leaks between sub-basins is conditioned by the vulnerability of soils to runoff and infiltration.Thus, soils favorable to runoff generate higher transfers of phosphorus and reduce those of nitrogen by reducing infiltration and promoting denitrification [79].In general, the levels of mineral nitrogen are low compared to the levels of phosphorus.This is because nitrogen is almost always one of the deficient nutrients in cultivated soils with fragile organic status and particularly in ferralitic and ferruginous soils [80].The contributions of nutrients of exogenous and organic origin remain higher than those of mineral origin.These results are in accordance with the work of Maïga et al. [17] who showed that the fertilizer used by farmers throughout the Lobo watershed was manure, mainly used in the Séguéla area located to the north of the basin.According to Maïga et al. [17], the nutrient losses are of natural origin.Thus, natural fertilizers could provide 10 to 50 kg ha −1 for nitrogen (N) and 0.15 to 0.75 kg ha −1 for phosphorus (P) per hectare in the Lobo watershed.SWAT model simulations gave values lower than measured values in the field.However, the results have proved that the points of high concentrations fit well in the area of high nutrient concentrations predicted by the SWAT model.However, the amount of nitrates and phosphorus accumulated can increase significantly over the years.This explains the concentrations measured in the field, which reached as high as 74 mg L −1 for nitrates and 66.6 mg L −1 for phosphorus, the level above water quality criteria set by the World Health Organization (WHO) which stipulates that water with more than 50 mg L −1 of nitrates is considered polluted.In addition, most unpolluted freshwater contains between 10 and 50 µg L −1 [81] in [82].This shows the state of pollution (eutrophication) of the Lobo water reservoir located in sub-basin 29, which is used for drinking water supply to Daloa city and around.These results agree with those from Dié [17] who highlighted the eutrophication phenomenon of the Lobo reservoir.These authors have strived to examine the causes and consequences of such a phenomenon.In the same vein, [83] showed that the Lobo undergoes very advanced permanent eutrophication.The study carried out by Lefebvre [3] highlighted phosphorus values that can reach 0.066 mg L −1 .The OECD (Organization for Economic Cooperation and Development) threshold values between 0.035 and 0.1 mg L −1 classify Lobo reservoir in the category of hypereutrophic lakes [84].

Conclusions
Considering the criteria for evaluating the flow setting, we note that the SWAT model is capable of simulating transfers in the Lobo watershed at Nibéhibé.The evaluation criteria R2, NSE, and PBIAS had values of 0.63, 0.62, and −8.1, respectively.P-factor P and R-factor gave 0.48 and 0.52.On the basis of these values obtained, it appears that the SWAT model is indicated for the simulation of flows and nutrient transfers in tropical agricultural watersheds, in particular that of the Lobo.The estimation of the average quantities of N and P applied to the crops gave 47.24 kg ha −1 of N and 21.25 kg ha −1 of P. The flows of inorganic N and soluble P vary from 0 to 0.049 kg ha −1 for inorganic N and from 0 to 0.31 kg ha −1 for soluble phosphorus.The northern and southern parts of the basin are more vulnerable to a large input of inorganic N when the central part remains vulnerable to the soluble phosphorus transfer.Simulations relating to the organic N transfer revealed values ranging from 0.2 to 5 kg ha −1 , the southern part showing the highest fluxes.As for the flow of organic phosphorus, they vary from 0.3 to 1.3 kg ha −1 .Unlike organic N, the greatest amounts are observed in the north.By comparing the flows of mineral (inorganic) nutrients with those of organic nutrients, one can say that the eutrophication observed in the Lobo water reservoir (sub-basin 29) would mainly have an exogenous origin (domestic, etc.) and organic.Predictions of nitrates concentrations and soluble phosphorus in the stream made using the SWAT model revealed concentrations of nitrates between 14 and 30 mg L −1 and soluble phosphorus between 0.1 and 1.8 mg L −1 .Field measurement points with high concentrations usually fall in areas of the stream with high concentrations simulated by SWAT.Apart from the limitations of this study related to the quality and quantity of the input data, it allowed us to gain a spatial understanding of the evolution and distribution of nutrients in the watershed.It also made it possible to identify the major inflow areas that would further contribute to eutrophication in the Lobo reservoir.To limit the impact of eutrophication on watercourses, it is necessary to regulate the use of fertilizers and phytosanitary products.In areas of high nutrient supply, it would be essential to implement vegetative filter strips around water bodies in order to capture the nutrients flow coming from upstream.These results constitute a solid basis which can serve as a basis in the fight against the phenomenon of eutrophication of water bodies in the watersheds of Côte d'Ivoire.

Figure 2 .
Figure 2. Essential input data for the SWAT model.Digital Elevation Model (DEM) The Digital Elevation Model (DEM) of 30 m × 30 m resolution, in the Northern Hemisphere projection UTM Zone 30, downloaded from the website https://search.earthdata.nasa.gov/search,accessed on 1 February 2018, has been used to extract the hydrographic network and delineate the watershed.

Figure 2 .
Figure 2. Essential input data for the SWAT model.

Figure 3 .
Figure 3. Nitrogen transformation process in the soil taken into account by the SWAT model [33].

Figure 3 .
Figure 3. Nitrogen transformation process in the soil taken into account by the SWAT model [33].

Figure 4 .
Figure 4. Phosphorus transformation process in the soil is taken into account by the SWAT model [33].

Figure 4 .
Figure 4. Phosphorus transformation process in the soil is taken into account by the SWAT model [33].

Figure 6 .
Figure 6.Annual mean spatial distribution of nutrients per sub-basin over the period 1981-1994 in the basin with 1 to 31 the number of each sub-basin.

Figure 6 .
Figure 6.Annual mean spatial distribution of nutrients per sub-basin over the period 1981-1994 in the basin with 1 to 31 the number of each sub-basin.
Appl.Sci.2023, 13, x FOR PEER REVIEW 16 of 22 watershed (Nibéhibé) through sub-basin 29, where the Lobo water reservoir is located.Sub-basin 15 is the most vulnerable, with concentrations up to 30 mg L −1 .Figure 7B presents the monthly average concentrations of soluble phosphorus in mg L −1 .The values vary from 0.1 to 1.8 mg L −1 (Figure 7B).The rivers in sub-watersheds 30 and 31 have the highest soluble phosphorus concentrations (1.3-1.8 mg L −1 ).Sub-basins 21, 26, 28 and 29 also recorded high concentrations (0.9-1.3 mg L −1 ).One noted that both for nitrates and soluble phosphorus, the sub-basin 29 containing the Lobo water reservoir recorded a significant contribution.A sampling campaign carried out made it possible to measure nitrate and phosphorus concentrations in surface water in the watershed.The measurement points have been represented on the maps of Figure 7A,B.(A).Nitrate concentrations (B).Soluble phosphorus concentrations

Figure 7 .
Figure 7. Spatial distribution of monthly mean nitrate and soluble phosphorus concentrations over the period 1981-1994 in the basin with 1 to 31 the number of each sub-basin.

Figure 7 .
Figure 7. Spatial distribution of monthly mean nitrate and soluble phosphorus concentrations over the period 1981-1994 in the basin with 1 to 31 the number of each sub-basin.

Table 2 .
Summary of selected SWAT parameters' global sensitivity and parameters' fitted values for simulating streamflow.

Table 3 .
Summary of selected SWAT parameters' global sensitivity and parameters' fitted values for simulating streamflow (continued from Table2).
a: Add parameter value to an existing value; v: Replace the given value with another; r: Multiply the existing parameter value by (1 + a given value of the same parameter).

Table 5 .
Nutrient requirements for crops.