Cowpea Crop Response to Mineral and Organic Fertilization in SE Spain

: Mineral fertilization is considered to be useful for improving soil fertility and yields. However, its use is linked to global warming and soil and water pollution by its rapid mobilization. On the other hand, organic fertilization is recommended to maintain or improve soil organic carbon and total nitrogen stocks while contributing to climate change mitigation. The main goal of this study was to assess the effect of two different fertilizer types, mineral and organic, during three cowpea crop cycles on the soil’s physicochemical properties, enzyme activities, crop yield, crop quality and nutritional composition when considering two cowpea cultivars (Feij ã o frade de ﬁo preto (FP) and Feij ã o frade de ﬁo claro (FC)). The use of mineral fertilizers was seen to contribute to improved soil fertility due to the increase in soil properties, such as recalcitrant carbon, total nitrogen, ammonium content, available K and available Mg. On the other hand, organic fertilizers only increased the nitrate content in the soil. There were no differences in terms of cowpea crop yield, quality and nutritional composition by fertilizer type. Thus, both fertilizer types contributed to the same crop yield and quality, and thus the use of organic fertilizers can result in a sustainable alternative for maintaining cowpea crop yield and quality.


Introduction
Cowpea (Vigna unguiculata (L.) Walp.) was first cultivated in Africa, and its cultivation has extended to tropical and semi-arid regions of the world [1]. Its seeds, leaves and green pods are valued for their nutritional properties, since its grains are rich in protein (about 25%), vitamins, fibers, minerals and other nutrients [2,3]. This legume is a drought-tolerant and warm weather crop which can be used not only in human food, but also for animal feed and the recovery of soil fertility through its use as green manure [4]. Cowpea as a legume crop needs a symbiotic association with effective nitrogen-fixing bacteria to attain maximal benefits from symbiotic N 2 fixation and thus save costly N fertilizers. Cowpea establishes symbiotic relationships with several species of the phyla Alpha-and Betaproteobacteria [5].
The use of mineral or organic fertilizers is directly related to soil fertility, soil carbon sequestration, greenhouse gas (GHG) emissions and crop yield. Mineral fertilizers play an important role in terms of the rapid increase of soil fertility and crop yield due to their high nutrient content and ease of availability [6]. However, this type of fertilization could not be the most effective or adequate, considering the maintenance of soil fertility over time, since an excessive application of mineral fertilizers can lead to low nutrient efficiency and deterioration of the soil and environment [7,8]. In terms of carbon sequestration, the use of mineral fertilizers has shown to have a rapid influence on the soil, since it can promote mineralization of a stable fraction of organic matter due mainly to their high N rates, or

Study Site and Experimental Design
This study was carried out in the same field in Cartagena in southeastern Spain (37 •  The soil was a Haplic Calcisol [27] with a clay loam texture, the main characteristics of which are shown in Table 1. The mean annual temperature of the study area was 18 • C, and the mean annual precipitation was 275 mm. Annual potential evapotranspiration surpassed 900 mm. Cowpea seeds were sown with a spacing of 100 cm between rows and 20 cm between plants (5 plants m −2 ). No herbicide treatment was provided, and the crops were kept free of weeds through hand-hoeing when necessary.
Every year, a surface application of 16,000 kg ha −1 of goat and sheep manure was carried out in all plots before sowing. This manure had the following characteristics: pH = 8.32 ± 0.07; electrical conductivity = 21.2 ± 1.01 mS cm −1 ; total organic carbon content = 307 ± 11 g kg −1 ; total nitrogen content = 13.4 ± 0.7 g kg −1 ; P 2 O 5 content = 1.8 ± 0.08%; and K 2 O content = 4.6 ± 0.2%. Fertilizer application in the cowpea plots started between two or three weeks after sowing and continued until harvesting. In the cowpea crops, 30 kg ha −1 of N and 2.4 kg ha −1 of P 2 O 5 were applied by mineral fertigation of ammonium nitrate (33.5% N) and monoammonium phosphate (61% P 2 O 5 , 12% N) and a commercial liquid organic fertilizer (Bombardier, Agroquímicos los Triviños, Spain; 10.7% Processes 2021, 9, 822 3 of 14 w/v N, 0.7% w/v P 2 O 5 ) in the organic fertigation. Cowpea residues were removed from the field and not applied in the soil as green manure. The irrigation was established on the basis of the evapotranspiration rate, crop coefficient and climatic conditions of rainfall.

Soil and Plant Sampling
The soil was sampled after harvesting the cowpea crop for three successive years. All plots were sampled at 0-20 cm (A p horizon). Three random soil samples per plot were collected and homogenized to obtain a composite sample, which was air-dried for 7 days, sieved (<2 mm) and stored at room temperature until analyzed. Enzyme activities were also measured in the air-dried samples, since this property is medium-term stable in stored air-dried samples of Mediterranean semiarid soils [28].
The cowpea crop yield was determined by weighing all the pods in each plot harvested when the seeds were dried at the end of the crop cycle. The weight of 100 seeds and protein content in the seeds were recorded as crop quality parameters.

Soil Analyses
The following parameters were measured: bulk density by the cylinder method [29]; soil pH and electrical conductivity (EC) in deionized water (1:2.5 and 1:5 w/v, respectively); soil texture by the Bouyoucos method [30]; equivalent calcium carbonate using the volumetric method (Bernard calcimeter) [31]; SOC by the wet oxidation method using K 2 Cr 2 O 7 [32]; recalcitrant carbon (RC) and labile carbon (LC) by the method of double acid hydrolysis [33]; aggregate stability (AS) by the method proposed by [34], based on the application of simulated rainfall with a known intensity; total nitrogen (N t ) by the Kjeldahl method [35]; cation exchange capacity using BaCl 2 as an exchangeable salt [36]; NO 3 − , extracted with deionized water in a 1:10 soil-extractant ratio [37] and measured by ion chromatography (Metrohm 861); NH 4 + , extracted with 2 M KCl in a 1:10 soil-extractant ratio [37] and colorimetrically measured [38]; and available phosphorus (P), extracted according to the Burriel-Hernando method [39] using a Burriel-Hernando solution (0.2 g CaCO 3 , 0.17 g MgCO 3 , 5 mL glacial acetic acid and 0.2 mL H 2 SO 4 in 2 L deionized water) in a 1:25 soil-extractant ratio. The available boron (B) was determined in the deionized water extract [40]. Exchangeable Ca, Mg, Na and K were determined in the BaCl 2 extract for CEC, and the P, B, Ca, Mg, Na and K concentrations were measured using ICP-MS (Agilent 7500CE). The β-glucosidase activity was based on the determination of p-nitrophenol released after incubation at 37 • C with β-D-glucopyranoside [41], while β-glucosaminidase Processes 2021, 9, 822 4 of 14 activity was based on the determination of p-nitrophenol released after incubation with pnitrofenil-β-D-glucopyranoside at 37 • C [42]. The dehydrogenase activity was determined using p-iodonitrotetrazolium chloride as a substrate and measuring the absorbance of the iodonitrotetrazolium formazam (INTF) produced [43]. The arylesterase activity was based on the determination of p-nitrophenol released after incubation with p-nitrophenyl acetate at 37 • C [44]. The cellulase activity was assessed by the determination of gearbox sugars using amorphous cellulose as a substrate [45,46], and the urease activity was based on the determination of ammonium released after incubation of the soil with urea at 37 • C [47].

Plant Analyses
The plant samples were oven dried and ground (A11 Basic, IKA, Staufen, Germany) before incinerating at 500 • C. The ashes were dissolved in 0.6 N HNO 3 and analyzed for P, Ca, Mg, Na and K by ICP-MS (7500 CE, Agilent, Madrid, Spain). The nitrogen (N) content was determined by the Kjeldahl method [35]. NO 3 − was extracted with deionized water in a 1:50 plant-extractant ratio [37] and measured by ion chromatography (Metrohm 861, Metrohm Hispania, Madrid, Spain). The total organic carbon (TOC) was quantified by the total combustion.

Statistical Analyses
The data were checked to ensure normal distribution using the Kolmogorov-Smirnov test and log-transformed when necessary to ensure normal distribution. The data were submitted to two-way repeated ANOVA measures, with the year (2014, 2015 and 2016) as a within-subject factor and the cowpea cultivar (FP and FC) and fertilizer type (mineral and organic) as between-subject factors. The relationships among the properties were studied using Pearson's correlations. A principal components analysis (PCA) was performed with all data to study the structure of dependence and correlation established among the variables studied with both fertilizer types. Statistical analyses were performed with IBM SPSS version 22 software for Windows (IBM, New Yok, United States).

Soil Physicochemical Properties
As a general pattern, the crop year significantly affected all soil properties, with general variations along the three crop cycles in most properties (Tables 2 and 3). The SOC, Nt, AS, RC, NO 3 − and P were higher in 2016 than in the first two years. The NH 4 + showed significantly higher values during the first year, while LC, Ca, Mg, K and Na showed significantly higher values during the second year.
The cowpea cultivar only significantly influenced the RC and NH 4 + ( Table 2). Both soil properties were higher in the soils cultivated with FC than FP in 2015 and 2016. With regard to fertilizer type, the use of mineral fertilizers significantly contributed to higher values of N t , RC, NH 4 + , Mg and K, mainly in 2015 and 2016 compared with the organic fertilizers. However, the NO 3 − content was higher using organic rather than mineral fertilizers. There was no significant effect of fertilizer type on the SOC, AS, LC, Ca, P and Na. The interaction of the year and fertilizer type was significant with respect to AS, RC and NO 3 − . AS and NO 3 − showed the highest values using organic fertilizers during 2016, while RC showed the highest values when using mineral fertilizers during 2016. The interaction of the year and cowpea cultivar was significant for the RC, which was higher in the soil cultivated with FC during 2016 than during the last two years. The interaction of the cowpea cultivar and fertilizer type did not affect any of the soil properties. Finally, the interaction of the three factors (year, cowpea cultivar and fertilizer type) was significant for NH 4 + , with the highest values in the FP cultivar using mineral fertilizers during 2016. N t was positively correlated with the RC (r = 0.82, p < 0.01). The Ca content was positively correlated with the Mg (r = 0.96, p < 0.01) and K contents (r = 0.85, p < 0.01). The Na content was negatively correlated with the available P (r = −0.83, p < 0.01). Table 2. Soil organic carbon, total nitrogen, aggregates stability, recalcitrant carbon, labile carbon, ammonium and nitrate contents in the soil of cowpea crops using mineral and organic fertilizers over three years. Values are mean ± standard deviation (n = 4).

Cowpea Cultivar a Fertilizer
Type

Soil Enzyme Activities
Enzyme activities were significantly influenced by the crop year ( Table 4). The βglucosidase and cellulase activities showed the highest values during 2014. The dehydrogenase and arylesterase activities showed the highest values during 2015 and 2016, respectively. The activity of β-glucosaminidase was significantly higher during 2015 and 2016. The cowpea cultivar only significantly influenced β-glucosaminidase activity, with higher values in the FP than in the FC cultivar. The fertilizer type only significantly influenced the β-glucosidase activity, which was higher with the use of mineral fertilizer. The only interaction between factors that was significant was that of the year and fertilizer type for arylesterase activity. Arylesterase showed significantly higher activity in the cowpea cultivated with organic fertilizers during the first year. Dehydrogenase was positively correlated with Ca (r = 0.87, p < 0.01), Mg (r = 0.90, p < 0.01) and K (r = 0.88, p < 0.01). The other enzyme activities were not correlated with any other property.

Cowpea Yield, Crop Quality and Nutritional Characteristics
As a general trend, the year had a significant effect on the cowpea yield and quality and on the nutritional characteristics (Tables 5 and 6). The cowpea crop showed the lowest yield in 2015, which was related to a serious virosis episode in this crop during this crop cycle. The cowpea yield and weight of 100 seeds were highest during 2016, while the Ca, Mg, K and P contents showed higher values during 2015. The cowpea cultivar significantly affected the cowpea yield and weight of 100 seeds. The cowpea yield was significantly higher in the FC cultivar, while the weight of 100 seeds was higher in the FP cultivar. The fertilizer type did not significantly affect the cowpea yield, crop quality or nutritional characteristics, indicating that both organic and conventional fertilizers contributed to the same crop yield and quality every year. The interaction of the year and fertilizer was significant for Ca, Mg, K and P, which showed higher values using organic fertilizers during 2015. The other interactions among factors were not significant. The cowpea yield was negatively correlated with Ca (r = −0.83, p < 0.01), Mg (r = −0.81, p < 0.01), K (r = −0.88, p < 0.01) and P (r = −0.86, p < 0.01). Table 4. Dehydrogenase, β-glucosaminidase, β-glucosidase, arylesterase and cellulase activities in the soil of cowpea crops using mineral and organic fertilizers for three successive years. Values are mean ± standard deviation (n = 4).

Interrelationships between Soil and Crop Properties
The PCA performed with the studied soil physicochemical properties, enzyme activities, cowpea crop yield, crop quality and nutritional characteristics using both fertilizer types showed that 83.3% of the total variation could be explained by the first five PCs. PC1 (which explained 33.2% of the variation) and PC2 (which explained 19.1% of the variation) did not separate the mineral and organic fertilizers ( Figure 1A). PC1 was related to soilexchangeable K, Mg and Ca, soil NH 4 + , the nutrient content in the seeds and soil enzyme activities (Table 7). PC2 was related to soil B, soil AS, exchangeable Na, bioavailable P and enzyme activities. PC3, which explained 15.5% of the variation, slightly separated the mineral (positive scores) and organic fertilizers (negative scores) ( Figure 1B). PC3 was related with higher values of EC, N t and SOC and lower values of pH (Table 7). PC4 (which explained 9.5% of the variation) and PC5 (which explained 5.9% of the variation) did not separate the mineral and organic fertilizers ( Figure 1B,C) and were related to the protein content in the seeds, the weight of the seeds (PC4) and the soil LC (PC5). Therefore, PCA confirmed the lack of differences between both types of fertilizers in terms of crop properties and most soil properties during the three crop cycles. P and enzyme activities. PC3, which explained 15.5% of the variation, slightly separated the mineral (positive scores) and organic fertilizers (negative scores) ( Figure 1B). PC3 was related with higher values of EC, Nt and SOC and lower values of pH (Table 7). PC4 (which explained 9.5% of the variation) and PC5 (which explained 5.9% of the variation) did not separate the mineral and organic fertilizers ( Figure 1B,C) and were related to the protein content in the seeds, the weight of the seeds (PC4) and the soil LC (PC5). Therefore, PCA confirmed the lack of differences between both types of fertilizers in terms of crop properties and most soil properties during the three crop cycles.

Discussion
The inclusion of a cowpea crop for three years increased the soil's properties, accompanied by an enhance in the crop yield and quality over time. Cultivation of the cowpea acts as a source of different types of C through rhizodeposition, which stimulates the growth of rhizosphere microorganisms [48] and facilitates nutrient availability for the plant [49]. In addition to this, we observed the effects of the specific cowpea cultivar on the soil properties, soil enzyme activities and crop yield and quality. This is not surprising, as previous reports have observed that the release of root exudates is controlled by the plant genotype [50,51]; thus, the plant genotype could have different effects on the soil properties and soil enzyme activities. On the other hand, this effect would be even larger if cowpea residues were used as green manure, increasing the nutrients' availability through mineralization of its residues [52][53][54]. In this regard, in this study, we observed a negative relationship between the exchangeable cations and crop yield. In fact, the lowest cowpea crop yield was observed during 2015, which coincided with the highest nutrient availability in the soil (Tables 3 and 5). This could be explained by the fact that soil sampling was carried out after two weeks of harvesting, and thus the higher the cowpea production and the higher the quantity of nutrients absorbed from the soil by the plants to sustain that production, reducing their concentrations in the soil at the end of the crop cycle. These results indicate that there was no nutrient limitation under either of the fertilized systems, conventional or organic. Since soil nutrients are the main factors limiting the crop yield, a high crop yield mainly depends on nutrient availability through fertilizer application [55], also being highly efficient with the use of organic fertilizers.
As far as the fertilizer type is concerned, both had a similar effect on the microbial activity, with similar soil enzyme activities. In turn, this study showed through PCA that the application of mineral and organic fertilizers had the same effect on the soil's properties ( Figure 1). However, mineral fertilizers contributed to higher RC and N t contents, which was not associated with a higher crop yield. This may be explained by the fact that the application of mineral fertilizers provides a source of N directly available for crop growth, and this N input acts as a driver for the C cycle in soils since it increases biomass production. Thus, it maintains and builds up the SOC stocks [56][57][58]. On the other hand, organic fertilizers contributed to an increase in NO 3 − content in the soil. This result may be explained by the fact that the application of fertilizers in the cowpea crop concluded three weeks before soil sampling, and the use of organic fertilizers produced a greater residual effect than mineral fertilizers by continuous mineralization and nitrification in the soil [59]. In fact, some authors have reported that the residual effect from mineral N fertilizers after a long-term fertilization experiment was low, or even that no effect was observed [60][61][62].
The lack of the expected results with improved soil physicochemical and biological properties with the use of organic fertilizers could be explained by the fact that these fertilizers release nutrients slowly, and only a fraction of the N and other nutrients become available to plants in the first year after application [63,64]. However, the nutrient availability was enough to provide high crop yields, with similar values to those with mineral fertilizers. In terms of crop yield, the literature reviewed showed that mineral fertilization of the soil is useful for improving yields through the increase of SOC and N stocks [65,66]. However, this type of fertilization could not be the most effective in light of a high environmental cost [67]. In addition to this, a higher soil N concentration following mineral fertilization may induce microbe-mediated N transformation processes, leading to high N 2 O emissions [12] or the lack of synchrony between the N release and crop demand, which can stimulate N 2 O emissions [68]. As a consequence, the use of organic fertilizers has been strongly recommended as a replacement for some or all mineral fertilizers, since this kind of fertilization not only includes mitigation of global warming [69], but their use also results in the reduction of NO 3 − and can improve crop yields. Despite this, the results of our study indicated, on one hand, that the use of organic fertilizers led to a higher NO 3 − content in the soil due to a greater residual effect, and on the other hand, there were no significant differences in the crop yield, crop quality or nutritional characteristics due to the fertilizer type ( Figure 1). In this context, the use of organic fertilizers may be beneficial to crop yields in the long term [70]. However, there are a number of scientific papers which conclude that the use of organic fertilizers alone as substitutes for mineral fertilizers is not enough to maintain the productivity of high-yielding crops [71,72], while an effective solution could be the partial substitution of mineral fertilizers with organic fertilizers [73]. Irrigation and fertilization are two important factors for obtaining higher productivity [74]. Nevertheless, some other factors such as the cultivar, climate and soil type are also responsible for changing the yield. In this study, we confirmed that the crop yield was indeed influenced by the cowpea cultivar.

Conclusions
In conclusion, our results show that the cowpea cultivar influences the soil's physicochemical properties and crop yield and quality. Mineral fertilizers increased the RC, nitrogen forms (N t and NH 4 + ) and nutrient availability (Mg and K), but the concentration of soil nutrients with organic fertilizers was also suitable for providing high crop yield and quality. As a consequence, there were no significant differences in the cowpea yield, crop quality or nutritional characteristics by fertilizer type. Thus, although mineral fertilizers increased the soil organic matter and some nutrients in the short term, the nutrient availability and microbial activity were high with the use of both types of fertilizers, leading to the same crop yield. Funding: This research was funded by the European Commission through the Seventh Framework Programme project n • 613781 "Enhancing of legumes growing in Europe through sustainable cropping for protein supply for food and feed" (EUROLEGUME).