Impacts on Soil and Cowpea Plants Fertigated with Sanitary Sewage through Subsurface Drip Irrigation

: The application of sanitary sewage using subsurface drip irrigation can mitigate current challenges related to water availability and food production. However, before implementing these techniques, it is crucial to conduct studies to assess their impacts on soil and plants. The objective of this study was to evaluate chemical changes in the soil and the development of cowpeas subjected to sanitary sewage applied by drippers with different flow rates and installation depths. Drippers were positioned at various depths (0 to 30 cm) and operated with flow rates of 1.6 and 3.8 L h − 1 . Cowpeas were cultivated in pots with clayey soil, using synthetic sanitary sewage based on the maximum limit of nitrogen fertilization. Irrigation management was controlled in terms of soil moisture, which was monitored using TDR probes. The results indicated that reducing the depth of the drippers positively affected grain production and the development of cowpeas. Fertigation with sanitary sewage at greater depths increased soil phosphorus concentrations and base saturation. Dripper depth also influenced soil concentrations of phosphorus, potassium, calcium, and magnesium, while sodium concentrations decreased with greater depth. It is concluded that dripper flow rates did not impact soil chemical parameters or the agronomic characteristics of cowpeas. However, despite nutrient supply at greater depths, the subsurface drip irrigation system proved unsuitable for cowpea production in clayey soil.


Introduction
The increasing scarcity of water for irrigation has become a critical issue globally, particularly in arid and semi-arid regions.The United Nations has predicted that global demand for food will increase by 70%, and agricultural water consumption will rise by 19% by 2050 [1,2].Therefore, it is crucial to conserve water for agriculture to ensure food security.
The disposal of sanitary sewage in agriculture is a sustainable way to provide water and nutrients to plants while recycling nutrients in the soil [3].However, this disposal of sanitary sewage into the soil must be carefully managed, as uncontrolled application can lead to soil salinization, clay dispersion [4,5], and contamination of plants, humans, and the environment, especially when applied through methods such as sprinkler irrigation.
Sprinkler irrigation methods typically increase crop yields by applying large quantities of water and fertilizers; however, they do not promote water conservation or the preservation of the soil's ecological environment [6,7].Nevertheless, water and fertilizer conservation are crucial for agricultural development.In this context, surface drip irrigation (SDI) and subsurface drip irrigation (SSDI) are more effective alternatives.SDI and SSDI systems fall under the category of localized irrigation methods.These systems uniformly apply water directly to the plant root zone through emitters at frequent intervals [8].With Water 2024, 16, 1194 2 of 19 these systems, it is possible to provide water and nutrients close to the plant root system, leading to water savings and increased crop yields [9,10].
Moreover, SDI and SSDI facilitate the use of lower-quality water [11,12] since these systems exhibit high irrigation efficiency, minimizing water losses due to drift.This characteristic is crucial for the disposal of sanitary sewage in agriculture as it reduces the dispersion of pathogens, odors, and animal and human contact.Advantages related to cultivation and cultural practices include increased growth, yield, and production quality due to higher irrigation frequency and the deposition of less mobile nutrients near the root zone [13].
An important aspect to consider in these irrigation systems is the depth of dripper installation, as improper installation can affect the transport and distribution of moisture in the soil [14,15].In general, soil nutrients and the roots of major crops are distributed primarily in the 0-20 cm soil layer.Therefore, a lack of water at this depth can impact mineralization and nutrient availability in the soil [16], especially for phosphorus (P) and potassium (K), whose movements in the soil primarily occur through diffusion [17].In this context, changes in soil moisture dynamics, temperature, pH, and nutrients serve as effective indicators of soil quality [18].
Considering the importance of the safe disposal of sanitary sewage through drip irrigation, it is interesting to associate its study with crops of economic and social significance, such as cowpeas (Vigna unguiculate), primarily produced in the Northeast region of Brazil.This legume constitutes one of the main sources of plant protein, plays a crucial role in human nutrition and contributes significantly to employment and income generation.Currently, due to supply and demand issues exacerbated by the COVID-19 pandemic, the cost of cowpeas is significantly high in the market [19].Furthermore, cowpeas are a short-cycle crop with low water and soil fertility requirements [20], making them a suitable alternative for areas with sanitary sewage disposal.
Given the current challenges of water availability versus food production, the study of alternative and more sustainable agricultural production techniques becomes essential.One such approach is the disposal of sanitary sewage through drip irrigation installed at different depths in cowpea cultivation.Thus, the objective of the present study was to assess the soil chemical changes and the development of cowpeas, specifically the BRS Tumucumaque cultivar, fertigated with synthetic sanitary sewage through drippers with different flow rates and installed at varying depths.

Local Experiment Conditions
The experiment was conducted in a greenhouse located at the Laboratory of Hydraulic Engineering of the Department of Agricultural Engineering at the Federal University of Viçosa (UFV), situated in the municipality of Viçosa, MG, Brazil.The geographical coordinates are 20 • 46 ′ 18.97 ′′ S, 42 • 52 ′ 28.19 ′′ W, WGS-84 Datum, with an altitude of 651 m above sea level.The study period spanned from August 2019 to February 2020.Figure 1 depicts the variations in meteorological elements that occurred inside the greenhouse during the experimental period.

Treatments and Experimental Design
The experimental design was a randomized complete block design, with three replications, and the treatments were arranged in a split-plot scheme.The main plots consisted of emitters with different flow rates, and the subplots were composed of the depths of installation of the drip tubes.
The experimental plots consisted of emitters operating with flow rates of 1.6 L h −1 and 3.8 L h −1 .The subplots were composed of the following depths of drip tube installation: 0 (surface drip irrigation), 5, 10, 15, 20, 25, and 30 cm (subsurface drip irrigation).The drippers used were the AmnonDrip PC AS model (self-compensating and anti-siphon) from the NaanDanJain ® brand.The irrigation system was evaluated before the start of the experiments, and a Christiansen Uniformity Coefficient (CU) of 100% was found for both flow rates used [21].

Treatments and Experimental Design
The experimental design was a randomized complete block design, with three replications, and the treatments were arranged in a split-plot scheme.The main plots consisted of emitters with different flow rates, and the subplots were composed of the depths of installation of the drip tubes.
The experimental plots consisted of emitters operating with flow rates of 1.6 L h −1 and 3.8 L h −1 .The subplots were composed of the following depths of drip tube installation: 0 (surface drip irrigation), 5, 10, 15, 20, 25, and 30 cm (subsurface drip irrigation).The drippers used were the AmnonDrip PC AS model (self-compensating and anti-siphon) from the NaanDanJain ® brand.The irrigation system was evaluated before the start of the experiments, and a Christiansen Uniformity Coefficient (CU) of 100% was found for both flow rates used [21].
The experimental units consisted of pots cultivated with cowpeas of the BRS Tumucumaque cultivar sown with uncertified commercial seeds of the first generation (S1).The pots had volumes of 65 L, with a height of 60.5 cm and bottom and top diameters of 30 and 44 cm, respectively.The pots were filled with soil to a height of 55 cm, leaving a 5.5 cm rim.The useful soil volume was 58.9 L, and the useful area at the soil level was 0.1434 m 2 .

Experimental Procedure
The soil used in the experiments was collected from the Experimental Area of Irrigation and Drainage of the Department of Agricultural Engineering at the Federal University of Viçosa (UFV).The soil was classified as dystrophic Red-Yellow Latosol [22], and its chemical and physico-hydric characteristics are presented in Table 1.The experimental units consisted of pots cultivated with cowpeas of the BRS Tumucumaque cultivar sown with uncertified commercial seeds of the first generation (S1).The pots had volumes of 65 L, with a height of 60.5 cm and bottom and top diameters of 30 and 44 cm, respectively.The pots were filled with soil to a height of 55 cm, leaving a 5.5 cm rim.The useful soil volume was 58.9 L, and the useful area at the soil level was 0.1434 m 2 .

Experimental Procedure
The soil used in the experiments was collected from the Experimental Area of Irrigation and Drainage of the Department of Agricultural Engineering at the Federal University of Viçosa (UFV).The soil was classified as dystrophic Red-Yellow Latosol [22], and its chemical and physico-hydric characteristics are presented in Table 1.
The collected soil was sieved through a 1-inch screen to prevent coarse materials from entering the pots, thereby minimizing the formation of preferential flow paths.Hoses were placed inside the pots, and they were filled with soil without altering the depth of the drip tube installation.Finally, before conducting the experiment with the disposal of sanitary sewage at different drip depths, other irrigation tests were conducted five months prior to the start of the experiment.This allowed the soil to settle in the pots, facilitating the initiation of the study by providing a well-structured soil environment.
To simulate sanitary sewage, a synthetic solution was used following the method described by Nopens et al. [23].For the preparation of this solution, salts and ingredients were diluted in tap water.The compounds used for the preparation of synthetic sanitary sewage (SSS) are presented in Table 2.  Notes: 1 Mass quantity of salts and ingredients for the production of 1 L of synthetic sanitary sewage; * Hydrated sodium acetate. 2 Chemical oxygen demand.
The fertilization of cowpeas was calculated based on the results of the soil chemical analysis (Table 1).Initially, soil liming was performed with a dose of 100 g m −3 of limestone, equivalent to 320 kg ha −1 [24].Soil fertility correction followed the recommendations of Melo et al. [25] to meet the nutritional demand of cowpeas: 6.7 g m −3 of N (20 kg ha −1 ), 20 g m −3 of P 2 O 5 (60 kg ha −1 ), and 13.3 g m −3 of K 2 O (40 kg ha −1 ).All nitrogen and part of the potassium and phosphorus were applied via SSS, with nitrogen considered the limiting element [26].The remaining phosphorus and potassium were applied during cowpea sowing (Table 3) in the form of single superphosphate and potassium chloride, respectively.
Sowing took place on 26 August 2019, with six cowpea seeds per pot planted at a depth of three centimeters.Later, only two plants per pot were retained.Applications of SSS began 24 days after sowing, with one application every 4 days, totaling 7 SSS applications or 8.26 L of SSS per pot (Table 3).Irrigation management was carried out by considering the soil water balance.Daily monitoring of moisture was conducted by inserting time domain reflectometry (TDR) probes at five depths (7.5, 15, 22.5, 30, and 37.5 cm) on the sides of the pots.The surface drip irrigation treatment was used as a standard to calculate the water demand of the other treatments, thus standardizing the volume of water applied to all plots.The entire volume of applied SSS, calculated by the limiting element, was insufficient to meet the water demand throughout the cowpea cultivation cycle.Therefore, supplementation with good-quality water was carried out to replenish the water lost through evapotranspiration in the pots.

Experimental Evaluations 2.4.1. Productive Components and CO 2 and Water Emissions
Net CO 2 exchange rate soil (NCER) and water vapor flux assessments were conducted when all treatments were in the reproductive phase using an infrared gas analyzer (IRGA, model LC-Pro+, ADC BioScientific Ltd., Hoddesdon, UK).Measurements were taken directly in the soil with the help of a metal cylinder inserted into the soil profile, under natural temperature conditions and atmospheric CO 2 concentration.The soil temperature was measured with a digital infrared thermometer.
For each experimental unit, grain production per plant was determined and corrected to 13% moisture content.After the end of the productive cycle, cowpea plants were collected, and the leaves and petioles were removed.This procedure was carried out because some treatments were more developed and had already lost part of the leaves (senescence).Subsequently, the stems were dried in an oven at a temperature of 70 ± 3 • C until the dry stem mass per plant (DM-stem) could be obtained.During the same period, the aboveground parts of weeds were removed from each treatment, and the dry mass (DM-weed) was also determined.

Soil Chemical Analyses
At the end of the experiment, the soil from each treatment was divided into 5 layers: 0-10; 10-20; 20-30; 30-40; and 40-55 cm.For each layer, a soil sample was taken, and chemical analysis was carried out.A summary of the analyzed variables and their respective quantification methods is presented in Table 4.

Statistical Analyses
The collected data underwent analysis of variance, and when the F-test was significant at a 5% probability, regression analysis was performed for the drip emitter installation depth factor.The choice of regression models was based on the significance of the regression parameters, verified using the t-test at a 5% probability, and the model's ability to explain the behavior of the studied variable.To conduct an integrated analysis of the various measured variables, two multivariate analyses were employed: principal component analysis (PCA) and k-means clustering analysis.The choice of the number of k-means groups was made based on the Silhouette coefficient [27] and the behavior of treatments in the PCA.The k-means analysis was conducted with all variables under study, and the Euclidean distance measure was used.

Productive Components and CO 2 and Water Emissions
The cowpea crop was evaluated to understand the effects of subsurface drip irrigation when using sanitary sewage.The net CO 2 exchange rate soil (NCER), water vapor flux (WFlux), and soil surface temperature (Ts) provide information on the water status of the soil; that is, they enable an assessment of whether the water distributed in depth was being absorbed by the roots of the crop.The quantification of the dry matter of weeds (DM-weed) provides an indication of whether the applied water reaches the soil surface.The other parameters were evaluated to determine if the crop was taking advantage of the nutrients from the application of synthetic sanitary sewage.
It can be seen in Table 5 that there was no effect in the interaction between dripper flow and the depth of drip emitter installation for any parameter evaluated.The simple effect of dripper flow rates was also not observed for any parameter.In relation to the simple effect of the installation depth of the drippers, the variables WFlux, Ts, DM-weed, grain production per plant (Prod.), and dry matter of the cowpea stem per plant (DM-stem) were significantly influenced.However, the NCER was not affected by the depth of drip emitter installation.
Table 5.The table shows a summary of the analysis of variance for net CO 2 exchange rate soil (NCER), water vapor flux (WFlux), and soil surface temperature (Ts), grain production per plant (Prod.),dry matter of the cowpea stem per plant (DM-stem), dry matter of weeds (DM-weed) in treatments fertigated with synthetic sewage applied by drip emitters operating at different flow rates (flow) and installation depths (depths).Water vapor flux, cowpea grain production, and the dry matter of weeds decreased with the increasing depth of drip emitter installation, with quadratic models being adjusted to explain the behavior of these variables.Cowpea stem dry matter decreased linearly, and soil surface temperature increased linearly with increasing drip emitter installation depth (Figure 2).

Source of Variation
The water vapor flux was 0.199 mmol m −2 s −1 for surface drip irrigation (SDI) and decreased by 20.5% and 57.3% for drip emitter installation depths of 5 and 30 cm, respec-tively.As a result, the soil surface temperature was lower for the SDI treatment, measuring 28.2 • C (Figure 2).

DM-weed
1.06 × 10 Water vapor flux, cowpea grain production, and the dry matter of weeds decreased with the increasing depth of drip emitter installation, with quadratic models being adjusted to explain the behavior of these variables.Cowpea stem dry matter decreased linearly, and soil surface temperature increased linearly with increasing drip emitter installation depth (Figure 2).The water vapor flux was 0.199 mmol m −2 s −1 for surface drip irrigation (SDI) and decreased by 20.5% and 57.3% for drip emitter installation depths of 5 and 30 cm, respectively.As a result, the soil surface temperature was lower for the SDI treatment, measuring 28.2 °C (Figure 2).
Irrigation and fertigation with synthetic sewage through surface drip irrigation (SDI) provided a grain production per plant of 50.65 g.For drip emitter installations at 5 and 30 cm depths, grain production decreased by 18.6% and 56.1% compared to SDI, Irrigation and fertigation with synthetic sewage through surface drip irrigation (SDI) provided a grain production per plant of 50.65 g.For drip emitter installations at 5 and 30 cm depths, grain production decreased by 18.6% and 56.1% compared to SDI, respectively.The dry matter of the cowpea stem was 21.90 g per plant for the SDI treatment, decreasing by 8.4% at the 5 cm depth and 50.1% at the 30 cm depth (Figure 2).SDI favored the development of weeds, resulting in 6.18 g of dry matter of weeds per experimental plot.However, at a drip emitter installation depth of 5 cm, there was a 51% reduction in weed development (Figure 2).

Phosphorus, Potassium, and Sodium
To adequately evaluate fertigation with domestic sewage applied via subsurface drip, it is essential to carry out a comprehensive analysis of the different soil parameters.These parameters include an assessment of essential nutrients to ensure that plants receive adequate nutrition during the irrigation process.Furthermore, it is important to monitor soil quality to ensure an environment conducive to healthy plant growth.The analysis of these parameters allows a comprehensive understanding of the effects of fertigation with domestic sewage, assisting in decision making to optimize agricultural production and ensure the sustainability of irrigation systems.
The variable sodium concentration in the 0-10 cm soil layer (Na L0-10), supplied through fertigation with SSS, was significantly influenced by the interaction of different flow rates and drip emitter installation depths.After unfolding the interaction, it is observed that in treatments with drip emitters installed at 5 cm, Na L0-10 was higher for the flow rate of 1.6 L h −1 compared to the flow rate of 3.8 L h −1 .However, for the drip emitter at a 15 cm depth, the concentration of Na L0-10 was higher for the flow rate of 3.8 L h −1 (Table 6).
The concentration of potassium in the 20-30 cm, 30-40 cm, and 40-55 cm soil layers (K L20-30, K L30-40, and K L40-55) received a unique effect from the drip emitter flow rates.The flow rate of 3.8 L h −1 promoted a higher potassium concentration compared to the flow rate of 1.6 L h −1 (Table 6).It can also be observed in this table that the concentration of phosphorus in the 0-10 and 30-40 cm layers, potassium in all layers, and sodium in the 0-40 cm layers were influenced by different drip emitter installation depths in terms of fertigation with SSS and mineral fertilization.Table 6.The table shows a summary of the analysis of variance for phosphorus (P), potassium (K), and sodium (Na) at different soil layers, L0-10; L10-20; L20-30; L30-40; L40-55 cm, and comparison of mean sodium and potassium supplied through fertigation with synthetic sewage applied by drip emitters with different flow rates (flow) and installation depths (depths).It can be observed in Figure 3 that the increase in drip emitter installation depth led to a linear increment in phosphorus concentration in the 0-10 and 30-40 cm layers.The concentration of P L0-10 at the end of the experiment was 16.37 mg dm −3 for SDI and 27.60 mg dm −3 for drip irrigation at a 30 cm depth.Considering the phosphorus inputs through mineral fertilization that the 0-10 cm layer received at sowing (Table 3), an expected increment in phosphorus concentration of 20.43 mg dm −3 was anticipated.In fertigation with SSS, the expected phosphorus increment was 5.96 mg dm −3 .The increase in P L0-10 concentration in the treatment with drip irrigation at a 30 cm depth indicates low phosphorus extraction in the 0-10 cm soil layer, directly influencing grain production and vegetative development of the aboveground cowpea parts (Figures 2 and 3).

Variables
The concentration of P L30-40 at the end of the experiment was 2.27 mg dm −3 for SDI and 2.92 mg dm −3 for drip irrigation at a 30 cm depth (Figure 3).Considering buried drip irrigation at 30 cm, the expected increment in P concentration in the 20-30 cm and 30-40 cm layers through fertigation with SSS was 2.98 mg dm −3 (Table 3).Therefore, it was observed that fertigation with SSS at a 30 cm depth promoted a higher final concentration of phosphorus in the 20-30 cm and 30-40 cm layers compared to fertigation at other depths (Figure 3).The sodium concentration in the soil from the L0-10 to L30-40 layers decreased with increasing drip emitter depth (Figure 3).According to the regression equation, the highest sodium concentration in the 0-10 cm layer and the lowest flow rate (1.6 L h −1 ) was 16.70 mg dm −3 for SDI, and the lowest concentration was 7.25 mg dm −3 for drip irrigation installed at a depth of 23 cm.Considering the same soil layer, for the highest flow rate (3.8 L h −1 ), the highest sodium concentration was 15.00 mg dm −3 for SDI, and the lowest concentration was 7.32 mg dm −3 for drip irrigation at a depth of 30 cm (Figure 3).Increasing the installation depth of the drip emitters led to a linear increment in potassium concentration in the 10-20 cm and 40-55 cm layers.For the potassium concentration in the 0-10 cm, 20-30 cm, and 30-40 cm layers, there is initially an increase in potassium concentration with the increase in the installation depths of the drip emitters, followed by stabilization and a slight decrease in potassium concentration for the greater depths of the drip emitters (Figure 3).
The lowest potassium concentration in all soil layers, after mineral fertilization and fertigation with SSS, occurred via the surface drip irrigation system, with values of 27.10, Water 2024, 16, 1194 10 of 19 11.40, 5.85, 5.62, and 5.72 mg dm −3 for the 0-10 cm, 10-20 cm, 20-30 cm, 30-40 cm, and 40-55 cm layers, respectively.According to the regression equations (Figure 3), the highest potassium concentrations in the 0-10 cm layer were 84.61 mg dm −3 for the emitter installed at a depth of 22 cm; in the 10-20 cm layer, it was 29.73 mg dm −3 for the emitter at a 30 cm depth; in the 20-30 cm layer, it was 14.39 mg dm −3 for the emitter at a 25 cm depth; in the 30-40 cm layer, it was 13.18 mg dm −3 with the emitter at a 20 cm depth; and in the 40-55 cm layer, it was 12.47 mg dm −3 with the emitter at a 30 cm depth.
The sodium concentration in the soil from the L0-10 to L30-40 layers decreased with increasing drip emitter depth (Figure 3).According to the regression equation, the highest sodium concentration in the 0-10 cm layer and the lowest flow rate (1.6 L h −1 ) was 16.70 mg dm −3 for SDI, and the lowest concentration was 7.25 mg dm −3 for drip irrigation installed at a depth of 23 cm.Considering the same soil layer, for the highest flow rate (3.8 L h −1 ), the highest sodium concentration was 15.00 mg dm −3 for SDI, and the lowest concentration was 7.32 mg dm −3 for drip irrigation at a depth of 30 cm (Figure 3).
The lowest sodium concentration in the soil layers between 10 and 40 cm after fertigation with SSS occurred via the subsurface drip irrigation system at a depth of 30 cm, with values of 0.13, 0.91, and 3.98 mg dm −3 , respectively.The highest sodium concentration in the 10-20 cm, 20-30 cm, and 30-40 cm layers occurred with the SDI treatment, with values of 6.29, 7.62, and 9.26 mg dm −3 , respectively (Figure 3).

Calcium, Magnesium, and Sodium Saturation Index
It can be seen in Table 7 that there was no effect of the interaction of dripper flow and depth of drip emitter installation for any parameter evaluated.However, isolated effects of dripper flow and depth of drip emitter installation were verified.The variables calcium (Ca) concentration in layers 0 to 40 cm, magnesium (Mg) concentration in layers 0 to 20 cm, and sodium saturation index (SSI) in layers 0 to 40 cm were significantly influenced by the installation depths of the drip emitters.
Table 7.The table shows a summary of the analysis of variance for calcium (Ca), magnesium (Mg), and sodium saturation index (SSI) in different soil layers: L0-10, L10-20, L20-30, L30-40, and L40-55 cm, fertigated with synthetic sewage applied by drip emitters operating at different flow rates (flow) and installation depths (depths).In Figure 4, it can be observed that the concentration of Ca in the 20-30 cm layer increased linearly with greater depths of drip emitter installation.The lowest concentration was 1.24 cmol c dm −3 , which was found in the SDI system.In the 30-40 cm layer, the adjusted model was quadratic, with the lowest concentration of Ca being 1.22 cmol c dm −3 for the SDI treatment, and the highest concentration was 1.40 cmol c dm −3 with the emitter installed at 20 cm depth.Analyzing the magnesium concentration in the first soil layer, a quadratic model was fitted as a function of the depths of the drip emitter installation.Considering this model, the lowest Mg concentration was 0.65 cmolc dm −3 for SDI, and the highest concentration was 0.93 cmolc dm −3 for drip irrigation at 24 cm depth (Figure 4).

Source of Variation
In the 10-20 cm layer, the Mg concentration increased linearly with drip emitter depths.The lowest concentration was 0.46 cmolc dm −3 for SDI, and the highest concentration was 0.61 cmolc dm −3 for drip irrigation at a 30 cm depth (Figure 4).
The values of SSI in the 0 to 40 cm layers decreased with increasing drip emitter depth (Figure 4).According to the regression equations, the highest SSI in the 0-10 cm layer was Analyzing the magnesium concentration in the first soil layer, a quadratic model was fitted as a function of the depths of the drip emitter installation.Considering this model, the lowest Mg concentration was 0.65 cmol c dm −3 for SDI, and the highest concentration was 0.93 cmol c dm −3 for drip irrigation at 24 cm depth (Figure 4).
In the 10-20 cm layer, the Mg concentration increased linearly with drip emitter depths.The lowest concentration was 0.46 cmol c dm −3 for SDI, and the highest concentration was 0.61 cmol c dm −3 for drip irrigation at a 30 cm depth (Figure 4).
The values of SSI in the 0 to 40 cm layers decreased with increasing drip emitter depth (Figure 4).According to the regression equations, the highest SSI in the 0-10 cm layer was 1.25% for SDI, and the lowest SSI was 0.53% for drip irrigation at 24 cm depth.The highest SSI values in the 10-20 cm, 20-30 cm, and 30-40 cm soil layers, after fertigation with synthetic sanitary sewage, occurred via the SDI system, resulting in 0.63%, 0.77%, and 0.92%, respectively.The lowest SSI values in the 10-20 cm, 20-30 cm, and 30-40 cm layers occurred with the treatment via subsurface drip irrigation at a 30 cm depth, with values of 0.01%, 0.09%, and 0.39%, respectively (Figure 4).

Hydrogenionic Potential, Electrical Conductivity, Potential Acidity, and Base Saturation
It can be seen in the summary table of the analysis of variance provided in Table 8, that only the potential acidity (H + Al) in the 30-40 cm soil layer was affected by the interaction between dripper flow and depth of drip emitter installation.In treatments with drip emitters installed at 5, 15, 25, and 30 cm, H + Al L30-40 was higher for the flow rate of 1.6 L h −1 compared to the flow rate of 3.8 L h −1 .The potential acidity (H + Al) in the 10-20 cm layer of the soil was the only parameter present in Table 8 that suffered an isolated effect from the different flow rates of the drippers.H + Al L10-20 for the flow rate of 1.6 L h −1 was 2.56 cmol c dm −3 , and H + Al L10-20 for the flow rate of 3.8 L h −1 was 2.44 cmol c dm −3 (Table 8).It can also be observed from this table that the soil pH in the 30-40 cm layer, electrical conductivity (EC) in the 0-10 cm layer, potential acidity (H + Al) in the 20-30 cm and 30-40 cm layers, and base saturation (V) in the 0-10 cm, 10-20 cm, 30-40 cm, and 40-55 cm soil layers were significantly influenced by the depths of drip emitter installation.
Figure 5 shows that the soil pH in the 30-40 cm layer exhibited quadratic behavior as a function of the different depths of drip emitter installation.The lowest pH found was 5.31 for the treatment with SDI, and according to the regression equation, the highest pH was 5.47 for the drip emitter installed at a 17 cm depth.The potential acidity (H + Al) in the 10-20 cm layer of the soil was the only parameter present in Table 8 that suffered an isolated effect from the different flow rates of the drippers.H + Al L10-20 for the flow rate of 1.6 L h −1 was 2.56 cmolc dm −3 , and H + Al L10-20 for the flow rate of 3.8 L h −1 was 2.44 cmolc dm −3 (Table 8).It can also be observed from this table that the soil pH in the 30-40 cm layer, electrical conductivity (EC) in the 0-10 cm layer, potential acidity (H + Al) in the 20-30 cm and 30-40 cm layers, and base saturation (V) in the 0-10 cm, 10-20 cm, 30-40 cm, and 40-55 cm soil layers were significantly influenced by the depths of drip emitter installation.
Figure 5 shows that the soil pH in the 30-40 cm layer exhibited quadratic behavior as a function of the different depths of drip emitter installation.The lowest pH found was 5.31 for the treatment with SDI, and according to the regression equation, the highest pH was 5.47 for the drip emitter installed at a 17 cm depth.The lowest concentrations of phosphorus, potassium, and magnesium in the 0-10 cm soil layer for the SDI system (Figures 3 and 4) contributed to the lowest electrical conduc-tivity (EC) for this system, which was 263 µS cm −1 .The highest EC was 554.80 µS cm −1 for buried drip irrigation at a 20 cm depth (Figure 5).In the same figure and, according to the regression equation, it can be observed that the highest values of H + Al in the 20-30 cm and 30-40 cm soil layers were found in the surface drip system, at 2.59 and 2.73 cmol c dm −3 , respectively.According to the regression equation, the lowest values of H + Al in the 20-30 cm and 30-40 cm layers occurred in treatments with drip emitters installed at 15 cm depth, being 2.43 and 2.47 cmol c dm −3 , respectively (Figure 5).
Base saturation (V) is an indicator of the overall soil fertility conditions.Due to higher productivity and nutrient extraction in the 0-10 cm soil layer, the SDI system exhibited the lowest V value, which corresponded to 56.60%.The highest V in the 0-10 cm layer was 64.94% for the treatment with the drip emitter at 20 cm depth (Figure 5).For the 30-40 cm layer and considering the regression equation adjusted to the data, the lowest V was 40.3% for SDI, and the highest V was 44.91% for the system at a 17 cm depth.According to the regression equation, V in the 40-55 cm layer increased linearly with increasing depths of drip emitters, with the lowest V being 39.6% for SDI and the highest V being 45% for drip irrigation at 30 cm depth.

Principal Component Analysis and Cluster Analysis
For a better understanding of the influence of drip emitter installation depths and emitter flow rates on the sets of variables studied, principal component analysis (PCA), along with k-means cluster analysis, was performed for all variables together (Figure 6).The lowest concentrations of phosphorus, potassium, and magnesium in the 0-10 cm soil layer for the SDI system (Figures 3 and 4) contributed to the lowest electrical conductivity (EC) for this system, which was 263 µS cm −1 .The highest EC was 554.80 µS cm −1 for buried drip irrigation at a 20 cm depth (Figure 5).In the same figure and, according to the regression equation, it can be observed that the highest values of H + Al in the 20-30 cm and 30-40 cm soil layers were found in the surface drip system, at 2.59 and 2.73 cmolc dm −3 , respectively.According to the regression equation, the lowest values of H + Al in the 20-30 cm and 30-40 cm layers occurred in treatments with drip emitters installed at 15 cm depth, being 2.43 and 2.47 cmolc dm −3 , respectively (Figure 5).
Base saturation (V) is an indicator of the overall soil fertility conditions.Due to higher productivity and nutrient extraction in the 0-10 cm soil layer, the SDI system exhibited the lowest V value, which corresponded to 56.60%.The highest V in the 0-10 cm layer was 64.94% for the treatment with the drip emitter at 20 cm depth (Figure 5).For the 30-40 cm layer and considering the regression equation adjusted to the data, the lowest V was 40.3% for SDI, and the highest V was 44.91% for the system at a 17 cm depth.According to the regression equation, V in the 40-55 cm layer increased linearly with increasing depths of drip emitters, with the lowest V being 39.6% for SDI and the highest V being 45% for drip irrigation at 30 cm depth.

Principal Component Analysis and Cluster Analysis
For a better understanding of the influence of drip emitter installation depths and emitter flow rates on the sets of variables studied, principal component analysis (PCA), along with k-means cluster analysis, was performed for all variables together (Figure 6).Overall, it can be observed that the points corresponding to treatments with shallower depths are separated from points of greater depth.However, for flow rates, there was no well-defined pattern, corroborating with the previously presented results, which indicated the small influence of flow rate on the studied variables.Based on the clustering, it can be observed that, generally, the treatments with SDI were grouped together, while the other treatments, corresponding to the SSDI, were divided into two other groups.

Discussion
In the present study, it was observed that increasing the depth of drip emitter installation negatively influenced grain yield and vegetative development of cowpea BRS Tumucumaque (Figure 2).Deeper placement of drip emitters leads to reduced water availability in the soil layer of 0-20 cm.This occurs because subsurface drip irrigation (SSDI) exhibits different soil water dynamics compared to surface drip irrigation (SDI).This dynamic impacts soil water content [15,28] and possibly negatively interferes with the absorption of nutrients applied via surface, such as liming and planting fertilization.Overall, it can be observed that the points corresponding to treatments with shallower depths are separated from points of greater depth.However, for flow rates, there was no well-defined pattern, corroborating with the previously presented results, which indicated the small influence of flow rate on the studied variables.Based on the clustering, it can be observed that, generally, the treatments with SDI were grouped together, while the other treatments, corresponding to the SSDI, were divided into two other groups.

Discussion
In the present study, it was observed that increasing the depth of drip emitter installation negatively influenced grain yield and vegetative development of cowpea BRS Tumucumaque (Figure 2).Deeper placement of drip emitters leads to reduced water availability in the soil layer of 0-20 cm.This occurs because subsurface drip irrigation (SSDI) exhibits different soil water dynamics compared to surface drip irrigation (SDI).This dynamic impacts soil water content [15,28] and possibly negatively interferes with the absorption of nutrients applied via surface, such as liming and planting fertilization.
The reduction in soil surface moisture with increasing depths of drip emitter installation is depicted in Figure 2, represented by the water vapor flux.The decrease in water vapor flux due to the increasing depth of drip emitter installation is attributed to the reduction in water content in the upper soil layers.The uniform distribution of soil moisture exerts a significant influence on environmental factors such as nutrient availability, pH, and the temperature of crop root zones [18], thus affecting root growth and crop productivity.Regarding subsurface drip irrigation systems, Lamm [16] reported that when the soil surface layer remains dry, nutrients applied at planting become unavailable to the crop, resulting in reduced yields.
Another important factor regarding soil surface moisture is thermal regulation (Figure 2).This condition is favorable for the germination and proper establishment of cowpea during the hottest periods of the year because water present on the surface mitigates soil temperature [29].Lower soil temperature reduces hypocotyl constriction and minimizes other physicochemical damage to seedlings, negatively affecting their establishment.The optimal temperature range for cowpea is considered to be between 20 and 30 • C [30].Considering this and analyzing Figure 2, it is observed that only the treatment with SDI falls within this range.
Surface wetting of the soil with surface drip irrigation also has its disadvantages [31].One of them is the greater loss of water to the atmosphere via evaporation.Another disadvantage is the promotion of weed seed banks, thus contributing to their establishment and development (Figure 2).However, this condition was not sufficient to compromise cowpea production compared to the other treatments.
Wu et al. [17] found, in sandy soil, that the accumulation of maize dry matter was higher in treatments with surface drip irrigation compared to subsurface drip irrigation systems.In the same study, conducted in clayey soil, the authors observed that maize production in treatments with SDI and SSDI was superior to non-irrigated cultivation by 17% and 12%, respectively.
The main advantage of fertigation with synthetic sanitary sewage through SSDI is the availability of nutrients at the desired depth, especially of less mobile elements in the soil, such as phosphorus.It was observed in Figure 3 that there was an increase in phosphorus concentration in the 30-40 cm soil layer for the deepest drip emitter installation depths.
Wang et al. [18] found that in treatments with drip emitters installed at a depth of 20 cm, the total nitrogen and total phosphorus contents in tomato root were 1.18 and 1.47 times higher, respectively, compared to surface drip irrigation.Phosphorus is prone to fixation at the point of application; thus, fertigation application through drip emitters installed at 20 cm resulted in more assimilable phosphorus at depth compared to surface drip irrigation treatments [32].
The mineral fertilization applied at sowing and distributed on the soil surface (0-10 cm soil layer) accounted for 77.40% of the total phosphorus and 60.38% of the total potassium applied in the present study (Table 3).Considering that water is the primary means of nutrient transport to plants, it is noted that the extraction of these nutrients in the 0-10 cm soil layer decreased with increasing depths of drip emitter installation (Figure 3), inversely corroborating with the water vapor flux (Figure 2).
Wu et al. [17], studying the effect of fertigation by SDI and SSDI on maize cultivation, found that for clayey soil, the moisture content was lower in the SSDI treatment at a depth of 0-20 cm compared to SDI.Furthermore, the extraction of N, P, and K was higher under conditions with higher available moisture.These same authors also found that phosphorus and potassium extraction were greater in treatments with SDI compared to SSDI.Therefore, the low water availability in the 0-20 cm layer can affect soil mineralization and nutrient availability [16], especially for P and K, whose movements primarily occur through diffusion.
The potassium concentration in the 10-20 cm soil layer increased linearly with the increment of drip emitter installation depths.Considering the treatment averages, for drip irrigation installed at 20 cm depth, the potassium concentration was higher, indicating the ability of synthetic sanitary sewage to provide nutrients at desired soil depths (Figure 3).possibly induces oxygen deficiency in the rhizosphere and reduces root volume, imposing a restriction on the plant's ability to extract water and nutrients [14].
For adequate levels of germination and vigor, cowpea seeds can tolerate an EC of up to around 3.0 dS m −1 [30].Considering that the highest EC found in the present study was 554.80 µS cm −1 (Figure 5), which is approximately 0.55 dS m −1 , fertigation with synthetic sewage, supplemented with mineral fertilization, did not compromise the soil chemical condition for the studied period.
The increase in base saturation in the 30-40 cm and 40-55 cm soil layers coincides with the enhancement of fertility through synthetic sewage.However, in the 30-40 cm layer, the drippers at 25 and 30 cm depth showed a reduction in base saturation, indicating nutrient extraction (Figure 5).Similar to what was suggested for potassium leaching in the 40-55 cm layer (Figure 3), there was also an increase in base saturation at these greater depths of dripper installation in the 40-55 cm layer (Figure 5).
Finally, considering the effect of emitter flow rates, it is observed that they mainly influenced the potassium concentration in soil layers below 20 cm (Table 6).It was found that the final K concentration was higher in the treatment with a flow rate of 3.8 L h −1 compared to emitters with 1.6 L h −1 .Lower flow rate emitters exhibit better water distribution, especially in sandy-clayey soils.Thus, in these layers, there was possibly better water distribution, resulting in a more even distribution of nutrients.This favored greater potassium absorption for treatments with emitters of 1.6 L h −1 .Treatments with this same flow rate showed higher potential acidity in the 30-40 cm layer (Table 8) due to the release of H + for potassium absorption [35].
When the emitter application intensity exceeds the soil water infiltration rate, the pressure load around the emitter becomes positive, generating a backpressure that reduces the hydraulic gradient between the interior of the emitter and the soil, consequently reducing the emitter flow rate.In fine-textured soils, the backpressure effect is more pronounced, resulting in a greater reduction in the emitter flow rate.This is one of the main disadvantages of the SSDI system, as it can affect irrigation performance.When water is applied in different quantities, non-uniformities occur in both time and space, consequently affecting the nutrient application rates as well [23].
With principal component and clustering analyses, it is possible to observe that, in general, SDI stands apart from the other depths of SSDI.However, three groups were formed, composed of SDI, intermediate depths of 5 and 10 cm, and greater installation depths of 15, 20, 25, and 30 cm (Figure 6).
An important decision in the use of SSDI for sewage disposal concerns the depth of emitter installation.The appropriate installation depth varies according to the crop, soil type, water source, climate, cultural practices, and irrigation designer preferences.
It was observed that fertigation with synthetic sewage through SDI enabled a higher grain yield of cowpea, which should be related to higher extraction of P, K, Ca, and Mg due to the lower final concentration of these elements in the soil in each studied layer.Therefore, SSDI was not suitable for cowpea production in sandy-clayey soil and in greenhouse conditions.However, it was effective in providing nutrients at depth via synthetic sewage.The responses to soil chemical alterations were not clear for the two studied flow rates.Generally, a slight advantage was observed for the flow rate of 1.6 L h −1 compared to the flow rate of 3.8 L h −1 .Thus, further studies are needed to investigate the effect of different flow rates in SSDI systems.Studies should also be conducted using original domestic sewage to evaluate the microbiological contamination of soil, plants, and operators.

Conclusions
Reducing the depth of emitter installation positively influences grain production and vegetative development of cowpea cultivated in sandy-clayey soil and in greenhouse conditions.
Fertigation with synthetic sewage at depth increases the concentration of phosphorus in the soil.

Figure 1 .
Figure 1.Environmental conditions during the experimental period.

Figure 2 .
Figure 2. The figure shows the net CO2 exchange rate soil (NCER), water vapor flux (WFlux), soil surface temperature, grain production per plant, dry matter of cowpea stem per plant, and dry matter of weeds in fertigated treatments with synthetic sewage applied by drip emitter installed at different depths.

Figure 2 .
Figure 2. The figure shows the net CO 2 exchange rate soil (NCER), water vapor flux (WFlux), soil surface temperature, grain production per plant, dry matter of cowpea stem per plant, and dry matter of weeds in fertigated treatments with synthetic sewage applied by drip emitter installed at different depths.

Water 2024 ,
16, x FOR PEER REVIEW 15 of 21

Figure 6 .
Figure 6.The figure shows the principal component analysis and cluster analysis of different emitter flow rates (Fq-1.6L h −1 and FQ-3.8 L h −1 ) and for different drip emitter installation depths (D0 to D30).Points of the same color belong to the same group.

Figure 6 .
Figure 6.The figure shows the principal component analysis and cluster analysis of different emitter flow rates (Fq-1.6L h −1 and FQ-3.8 L h −1 ) and for different drip emitter installation depths (D0 to D30).Points of the same color belong to the same group.

Table 1 .
The table shows the chemical and physico-hydric attributes of the soil before the implementation of the experiments.

Table 2 .
The table shows the theoretical composition and respective concentrations for the production of 1 L of synthetic sanitary sewage.

Table 3 .
The table shows the total volume of synthetic sanitary sewage (SSS) applied in each experimental unit and nutrients supplied by SSS and mineral fertilization. ¹

24 ,
16,x FOR PEER REVIEW 14 of 21 means followed by the same lowercase letter in the column do not differ from each other at 0.05 by the F-test.