Introduction
Ecological Agriculture (EA) aims for a natural and friendly relationship with the environment to promote plant and animal biodiversity (de los Ríos et al., 2016), unlike the Green Revolution, which is associated with the extensive use of chemical fertilizers and heavy-duty machinery (Ceccon, 2008). This has led agriculture to several negative impacts, creating simplified and unstable ecosystems and causing environmental deterioration (Rubio-Sánchez, 2016).
It is in this context that EA uses techniques which are favorable to the environment with changes in crop rotation and the use of green manures, vermicomposts (worm humus), poultry manure, biofertilizers, and efficient microorganisms (EM) in order to maintain soil productivity (Guiberteau & Labrador, 1991).
Another productive and environmentally friendly strategy is the implementation of agroforestry systems, also called SAF, a sustainable alternative that helps increase soil productivity (Casanova-Lugo et al., 2011). In addition, these systems contribute to multiple environmental services such as mitigating the conflict of land use in Colombia, especially considering that our country has a great diversity of species. This allows for the diversification of said systems, of which coffee production is one of the predominant modalities (Henao, 2020).
These traditional techniques use the resources that are available, namely organic fertilizers, in this case worm humus, which can be an efficient solution to replace mineral fertilizers, as they help maintain balance in the soil (Milanés et al., 2005). Another technique that has been implemented is the use of biofertilizers that improve and conserve the physical, chemical, and biological conditions of the soil, given that a large part of the organic matter comes from enzymatic digestion by microorganisms, and it is from this organic matter that we obtain the nutrients that are easily assimilated by the soil (Trejo-Escareño et al., 2013).
Moreover, EMs are also used as organic fertilizers in agriculture, with the great benefit of increased production, in addition to helping prevent and reduce the attack of pests and diseases. They are produced as a mixture of natural resources that is highly efficient and beneficial (Callisaya and Fernández, 2017). Finally, poultry manure provides benefits such as increased moisture retention, increased organic matter, and improved soil fertility through increased amounts of macro and micronutrients (González et al., 2013).
It has been suggested that the implementation of these EA techniques, as well as the different covering options, provides several benefits to the soil. According to the National Commission for Biodiversity Knowledge and Use (Comisión nacional para el conocimiento y uso de la biodiversidad) (CONABIO, 2015), shade-grown coffee is particularly promising in terms of conservation, with benefits such as greater biomass, more soil nutrients, and microbial biodiversity, fewer weeds and harmful insects, and a better water and microclimate balance. Trees like orange, banana, and guava, among others, generate benefits to the soil by providing organic matter produced by leaf litter and atmospheric nitrogen fixation (Ruelas-Mojardín, 2014). In Colombia, coffee is a product with special altitude, latitude, and climate conditions. This means that Colombian coffee culture has one fundamental requirement: the quality of the soil (FNC, 2012). This is why the agroecological conditions of San José de Pare, Boyacá, are suitable for coffee plantations; it has a cultivated area of 617 ha, with plants cultivated under the shade or coverings due to the high radiation that is emitted in this area (Alcaldía Municipal San José de Pare Boyacá, 2020).
Soil is one of the most important resources for life on the planet (García et al., 2012). Therefore, physicochemical indicators such as gravimetric moisture, pH, CEC, phosphorus, organic carbon, and nitrogen are so important, since they allow determining the degree of vulnerability of a given soil depending on the activities it has been exposed to. This includes the buffering capacity of the soil and the availability of water and nutrients for plants and microorganisms. Biological indicators (phosphate solubilizing bacteria, phosphate solubilizing fungi, cellulolytic bacteria, cellulolytic fungi, and nitrogen fixing bacteria) show the diversity and abundance in soil components based on factors such as the presence of micro and macroorganisms, bacteria, earthworms, and soil fungi (Bautista et al., 2004).
Another important indicator is soil enzymes, which are involved in microorganism activity. Together with other soil organisms, they are linked to energy transfer, biochemical processes, and environmental quality. Additionally, they are considered to be an important indicator of soil fertility (providing adequate amounts of nutrients) and agronomic, ecological, and anthropogenic influence (Doran & Parkin, 1994).
Finally, despite the importance of EA techniques for soils, there are still not enough studies showing possible changes in soil quality indicators as a result of the application of said techniques over time, specifically in coffee plantations under different coverings. For this reason, the main question of this research is: what changes occur in the soil quality indicators of coffee plantations under different coverings, which are associated with the implementation of agriculture ecological techniques in the municipality of San José de Pare, Boyacá?
Materials and methods
Description of the study area
This research was carried out in a rural area located in the Colombian municipality of San José de Pare, in the department of Boyacá, at an altitude of 1545 m.a.s.l. and average temperatures of 13 and 24 ºC, with the following geographical coordinates: 6º01.84' N longitude, 73º31' W latitude. This municipality is located in the Republic of Colombia, 217 km north of Bogotá DC, in the northeastern region of the department of Boyacá, and it shares borders with the municipalities of Chitaraque, Moniquirá, San José de Pare, Santana, and Togüi.
Experimental design and soil sampling
As shown in Figure 1, sampling was carried out by selecting three coffee plots (twelve coffee plants in each plot) under three different coverings: banana trees (CP), mandarin trees (CM), and three citrus fruit trees (CC) (mandarin, lemon, and orange). The EA techniques applied were: 1) worm humus (Californian earthworm) (Milanés et al., 2005); 2) biofertilizer, obtained from the mixture of legumes, molasses, yeast, chicha, poultry manure, rabbit manure, wood ash, and milk (Trejo-Escareño et al., 2013); 3) efficient microorganisms (EM) from a mix of leaf litter, rice grits, molasses, and milk (Callisaya & Fernández, 2017); 4) poultry manure (González et al., 2013); 5) Red clover green manures (Trifolium pratense) (Fonseca-López et al., 2020); and 6) cover crops (Geissert et al., 2017). As the EA techniques were being applied, soil samples were collected on a bimonthly basis (months 0, 2, 4, 6, 8, 10) for ten months, for a total of 6 sampling moments. Each plot provided three samples composed of 10 different subsamples, which were taken in zigzag from the upper 20 cm of the soil. The soil samples were then homogenized, labeled, and stored in bags at different temperatures depending on their purpose: the bags for physicochemical analysis at room temperature, those for microorganism abundance analysis at 4 °C, and those for enzymatic activities at -20 °C. The samples were processed in the shortest possible time after sampling (maximum one week after collection).
Physicochemical parameters of the soil
All samples underwent physicochemical analysis according to the following methods:
pH (potentiometer): 1:1, soil:water. (Bernal et al., 2014).
Moisture (gravimetric method): it consists of taking a soil sample, weighing it before and after drying, and calculating its moisture content. The soil sample is considered to be dry when its weight remains constant at a temperature of 105 ºC (IDEAM, 1994).
Phosphorus (Bray II): the soil phosphorus is extracted with 25 ml of the extractor solution of HCl 0.1N and NH4F 0.03N, shaken, filtered, and the available phosphorus is colorimetrically quantified by reaction with chloromolybdic acid and stannous chloride to form a blue complex whose intensity is proportional to the amount of the element present (García & Ballesteros, 2006).
Organic carbon (Walkley-Black): this method is based on the oxidation of the organic carbon in the soil by means of a solution of potassium dichromate (K2Cr2O4) and the heat of the reaction generated by mixing with concentra-ted sulfuric acid (H2SO4). After a waiting time, the mixture is diluted, phosphoric acid is added to avoid Fe3+ interferen-ces, and the residual potassium dichromate is titrated with ferrous sulfate (FeSO4 1N. pH 7) (Vela et al., 2012).
Nitrogen (Kjendahl): which consists of three stages digestion, distillation and titration, where the sample is left in digestion with pure sulfuric acid, then an absorbent solution such as boric acid is used for the distillation stage and finally a titration with 1N sulfuric acid is performed (Jarquín et al., 2011).
Abundance of functional groups of microorganisms
To determine the abundance of functional microorganism groups, we performed plating and plate counting in different media, to later quantify the colony forming units (CFU) per gram of dry soil. For the seeding and counting of phosphate solubilizing fungi and bacteria, we used the modified Sundara and Sinha medium (Sundara & Sinha, 1963). The procedure was carried out by incubating the bacteria at 28 ºC for 40 h and the fungi at 20 ºC for six days; for cellulolytic bacteria and fungi, we used a medium whose only source of carbon was 1% carboxymethyl cellulose, where bacteria were incubated at 28 ºC for 48 h and fungi at 20 ºC for six days (Avellaneda-Torres et al., 2014); and, finally, for the nitrogen-fixing bacteria, we used the N-free Rennie medium (Rennie, 1981) and CFU plate count by maintaining the temperature at 28 ºC for 48 h (Avellaneda-Torres et al., 2020).
Enzymatic activity
With the purpose of evaluating the enzymatic activities, we considered the association with phosphorus, carbon, and nitrogen cycles. For the phosphorus cycle, the study determined acid and alkaline phosphatase activities, which were quantified by determining the p-nitrophenol that was released after the incubation with p-nitrophenyl phosphate solution for 1 h at 37 ºC. The p-nitrophenol was photometrically determined at 400 nm. For the carbon cycle, the β-glucosidase activity was determined by quantifying the p-nitrophenol released after the incubation with a p-nitrophenyl glucoside solution for 1 h, at 37 ºC. The p-nitrophenol was photometrically determined at 400 nm. As for the nitrogen cycle, the urease activity was determined, which is based on the colorimetric determination at 690 nm of the released ammonia after incubation with a solution of urea for 2 h, at 37 ºC (Avellaneda-Torres et al., 2018). Finally, the study determined the protease activity by incubating the soil samples for 2 h at 50 ºC using casein as substrate (Caballero-Vanegas et al., 2018) after incubation for 2 h at 50 ºC and a pH 8.1. In this method, the released amino acids reacted with the Folin-Ciocalteu reagent in an alkaline solution to form a blue complex that was colorimetrically evaluated (700 nm) (Avellaneda-Torres et al., 2018).
Statistical analysis
To determine if there were significant differences in the physicochemical parameters, the abundance of functional groups of microorganisms, and the enzymatic activities throughout the months of application of EA techniques in the soils with coffee plantations under three different coverings (CP, CM, and CC), we carried out univariate and multivariate data analyses. To this effect, the assumptions of normality and homogeneity of variances were verified using the Kolmogorov-Smirnov and Barlett tests, respectively, followed by the Kruskal-Wallis tests with a post-hoc Mann-Whitney analysis to determine significant differences. The tests were performed with the Stats package (R Core Team, 2019). A principal component analysis (PCA) was also performed using the factoextra package (Kassambara & Mundt, 2020). The analyses were performed with the R statistical software, version 3.5.3. (R Core Team, 2019).
Physicochemical parameters
According to the results presented in Table 1, it appears that, on average, we can find a 29.3% for CP in the Ɵg, which shows that it is a soil with good moisture content, as well as an average of 37.4 and 48.7% for CM and CC, respectively. This means that, according to its water content, it is a usable soil. In general, the pH in registers strongly acid values between (5.0 and 5.7), thus indicating that some of the elements that are available in the soil do not have adequate levels. On the other hand, we found an average CEC of 26.4 meq.100g-1 for CP and 38.7 meq.100g-1 for CM, which indicates that it is a medium CEC soil, whereas 41.4 meq.100g-1 for CC imply a high CEC soil. Regarding the phosphorus content, an average of 30.6% was obtained for the three plots, indicating that phosphorus in the soil under study is at high levels. It was observed that carbon levels were high, as they ranged between 4.3 and 8.7%. Finally, nitrogen content values were high (between 0.4 and 0.9%).
According to Table 1, the Ɵg shows fluctuations throughout the months with each of the coverings. However, they are not statistically significant. Moreover, the pH of the soils did not change significantly for the three evaluated coverings as the months EA technique application elapsed. In addition, there were no significant changes in the CEC during this period. For organic carbon, a significant increase was observed in the three coverings. The same behavior was seen in the case of phosphorus, where a significant increase was evident in each of the coverings, with the highest statistical increase in month 10 for both phosphorus and carbon. There was no significant change for nitrogen over the months of application of the EA techniques.
Functional groups of microorganisms
The results for the abundance of functional groups of microorganisms show a statistically significant increase as the months of application of EA techniques elapsed in the case of phosphate solubilizing bacteria (Figure 2a). Phosphate solubilizing fungi (Figure 2b) show a significant increase, whereas cellulolytic bacteria (Figure 2c) did not show statistically significant changes during this period. Cellulolytic fungi (Figure 2d) increased statistically during the last months of study, and, finally, nitrogen-fixing bacteria (Figure 2e) did not show any significant changes.
Enzymatic activity
Figure 3 shows that there is no significant change in the enzymatic activity corresponding to alkaline phosphatase (Figure 3a), acid phosphatase (Figure 3b), β-glucosidase (Figure 3c), urease (Figure 3d), and protease (Figure 3e) during the months of application of the EA techniques. On the other hand, the statistical analyses for enzymes in each crop (banana, mandarin, and citrus) accompanying the coffee plantations showed the same results, which is why they are not presented in this paper.
Principal component analysis of different variables in response to EA techniques
The cumulative variance of the PCA in (Figure 4) was 60.2%, which constitutes an indicator of good quality in the analysis of the results. It can be observed that, towards the last two months of sampling and application of EA techniques, there is a correlation between the variables of phosphorus, organic carbon, phosphate solubilizing fungi, phosphate solubilizing bacteria, and cellulolytic fungi. pH, nitrogen, and enzymatic activities such as β-glucosidase and urease appear in a lower proportion. The behavior of these variables is opposite to that observed for gravimetric moisture, showing that the lower the soil moisture with respect to the months of application, the higher the percentages of the different variables opposite to it. Likewise, another inversely proportional relationship takes place between pH and phosphate solubilizing bacteria, which indicates that the count of these bacteria depends on the change in soil pH.
Component 1 has a variance of 60%, and, together with component 2 of variance 40, they represent 100% of the total variation.
Figure 5 shows that the accumulated variance of the PCA was 100%, which indicates a low stress analysis for the understanding of the variables. It is observed that CP is correlated with the protease and phosphate solubilizing bacteria. CM shows a correlation with variables such as organic carbon, CEC, cellulolytic fungi, nitrogen-fixing bacteria, acid and alkaline phosphatase, and β-glucosidase. Finally, there is a higher correlation for CC among variables such as gravimetric moisture, pH, phosphorus, nitrogen, microorganisms, phosphate solubilizing fungi, and cellulolytic carbon bacteria, as well as for enzymes such as urease.
Component 1 has a variance of 9.6 and groups 60% of the variation, and, together with component 2 of variance 6.4, which groups 40% of the variation, they represent 100% of the total variation.
Discussion
This research focused on evaluating whether the application of EA techniques generates significant changes in soil quality indicators in coffee plantations with different covering crops (CP, CM, and CC). It was found that variables such as phosphorus, organic carbon, phosphate solubilizing bacteria, phosphate solubilizing fungi, and soil cellulolytic fungi increased significantly during the months of application of the EA techniques, namely, vermicomposting, liquid humus, biofertilizer, EM, green manures, strengthening of polycultures, and composted poultry manure.
Said applications of EA techniques, together with agroforestry systems, integrate practices in which a woody component (tree or shrub) is combined with an agricultural one, thus allowing an improvement in the soils with respect to the production per unit area. In other words, it enhances productivity, thus obtaining more than a single product from the same soil unit (Henao, 2020) and promoting polyculture which is one of the practices of applied EA techniques.
The increase in phosphorus may be due to the enrichment of organic compounds from biofertilizers, which can improve phosphorus dynamics, resulting in increased nutrient availability (Rosatto et al., 2014).
In contrast, the increase in soil organic carbon may be due to the application of green manure and humus techniques, since they are carbon-rich, which allows a greater biointensity of the soil (Martínez et al., 2008). Likewise, this increase is linked to an important component such as humus, which can represent up to 90% of the total carbon in the soil (Jaramillo, 2011), meaning that adding organic matter to the soil makes the incorporation of vermicompost stand out due to its high stability, high content of bacterial fiber, and high content of nutrients that can be assimilated by plants (Salinas-Vásquez et al., 2014).
Furthermore, the use of EA techniques promotes the processes of some microorganisms that can synthesize organic acids, releasing phosphate ions from the exchange surface, thus increasing the concentration of these ions in the soil solution and therefore increasing their availability (Tapia & García, 2013). Additionally, the organic amendments used could cause an increase in soil phosphorus due to their derived composition (water, legumes, yeast, cow, and rabbit manure), given that they are associated with an increase in organic matter, which stores soil nutrients that improve its structure, as well as with the concentration of phosphorus in the soil, as plant roots absorb phosphorus (Múnera, 2014). Regarding biological indicators, the results obtained showed statistical increases in phosphate solubilizing bacteria and fungi, as well as in cellulolytic fungi. These results coincide with those found by other studies that evaluated the effect of the application of mycorrhizae in coffee growing systems (Avila et al., 2020).
Many soil microorganisms can transform insoluble phosphorus into forms that can be assimilated by plants, such as phosphate solubilizing bacteria. The increase in such bacteria may be related to the ionic forms of Pi (primary soil minerals) such as the H2PO4- ion, which can be absorbed by plants due to its solubility in water. These are dependent on pH, which ranges between 4.0 and 6.0 (Restrepo et al., 2015). Considering that the pH of this study was between 5.1 and 5.6, carrying out these processes is considered feasible. Likewise, this is correlated with the increase in inorganic phosphorus in the soil, which may be associated with a greater number of colony-forming units of phosphate solubilizing microorganisms as the application of EA techniques for the purposes of this study progressed.
As for phosphate solubilizing fungi, a significant increase was evidenced before the application of the EA techniques, mainly with the biofertilizer, which yielded positive responses for crops (Patiño & Sanclemente, 2014). These microorganisms transform insoluble mineral compounds into soluble minerals that are absorbed by plants (phosphorus solubilization) (Alejo-Iturvide et al., 2016).
Additionally, we observed a significant increase in cellulolytic fungi, possibly because of the increase in organic matter related to the application of the EA techniques, which promote an increase in cellulolytic microorganisms, and the presence of cellulose as a source for carbon, since this is the major component of plant biomass found in organic matter (Monserrate et al., 2014). These results match those found by other authors when comparing organic farming systems in coffee with conventional ones (Garzón et al., 2020). Additionally, efficient microorganisms (EM) secrete beneficial substances such as vitamins, organic acids, antioxidants, among others. These antioxidants promote the decomposition of organic matter and increase the content of humus in the soil (Liriano et al., 2015).
As for soil enzymatic activities, in general, it can be observed that there was no significant increase related to the application of EA techniques. In the case of acid and alkaline phosphatase, the absence of significant changes during the develop- ment of this research may be due to the fact that a high P content can inhibit the action mechanism of the enzymes by having a greater availability of nutrients in the soil in labile form (Narváez et al., 2010). The β-glucosidase enzyme, in turn, is one of the most abundant and important enzymes for soil quality. However, no significant change was found, which may be due to the fact that, in order to demonstrate its changes, longer times are required for the application of EA techniques (Gutiérrez et al., 2008).
On the other hand, regarding the covering crops that accompany the coffee crop (CP, CM, and CC), and as it can be observed in Figure 5 for ACP, the nitrogen-fixing bacteria, acid phosphatase, alkaline phosphatase, and β-glucosidase increased for CP. For CC, the pH, phosphorus, nitrogen, phosphate solubilizing fungi, cellulolytic bacteria, and urease also increased. The same behavior was observed for organic carbon, cellulolytic fungi, and acid phosphatase with both CC and CM. This has to do with the relationships between the increase in organic matter and the application of techniques such as vermicomposts, biofertilizers, efficient microorganisms, and green manures. This coincides with reports that soil coverings contribute to the stimulation and optimization of soil biological processes, improving soil fertility, counteracting erosion processes, and favoring the presence of populations of beneficial organisms (Puertas et al., 2008).
Moreover, it can be seen in Figure 4 (PCA) that there is a direct correlation for months 8 and 10 with nine of the variables that were analyzed, namely phosphorus, organic carbon, nitrogen, cellulolytic fungi, nitrogen-fixing bacteria, phosphate solubilizing fungi, and β-glucosidase, which are variables that increased due to the application of EA techniques because they generated a high content of organic matter and high biological activity in soils, thus highlighting the importance of applying these methods over time (Gutiérrez V. et al., 2008).
Conclusions
According to the univariate data analysis of the six variables evaluated on physicochemical properties, a significant increase in the phosphorus and organic carbon was reported in the coffee soils under study. It was also observed that, due to the application of the EA techniques, the five functional groups of quantified microorganisms showed a significant increase in the abundance of phosphate solubilizing bacteria and fungi, as well as cellulolytic fungi.
As for the multivariate data analysis, the PCA determined that, during months 8 and 10 of application of the EA techniques, there was a significant increase in the content of phosphorus, organic carbon, nitrogen, the abundance of cellulolytic fungi, nitrogen-fixing bacteria, phosphate solubilizing fungi, and β-glucosidase.