Highlights:
Land use scenarios evaluated: forest use, avocado with conventional use and organic management, and macadamia.
Physical and hydrological properties were evaluated at two depths (0 to 20 cm and 20 to 40 cm).
Physical and hydrological properties varied significantly among the four land use scenarios.
Infiltration decreased 40 to 70 % in agricultural systems with respect to forest use.
Land use and depth interaction was significant for bulk density.
Introduction
Global agriculture is characterized by intensive land use in conventional and organic management modalities, causing variations in physical and hydrological properties, regarded as key for agricultural productivity, thus deterioration involves soil chemical and biological activity (Babin et al., 2019; Di Prima et al., 2018; Sihi et al., 2017). About 37.4 % of the Earth's surface is used for agriculture (Centro de Estudios para el Desarrollo Rural Sustentables y la Soberanía Alimentaria [CEDRSSA, 2019]).
Mexico has 26 of the 32 soil groups recognized by the World Soil Resource Reference Base (IUSS Working Group WRB, 2006). This diversity can be explained by several factors such as geographic location, vegetation, topography and climate, which directly influence land use and management. Andosol soils derived from volcanic ash and are important for the development of agriculture, which has been one of the bases of economic growth, and a factor of social stability (IUSS Working Group WRB, 2016). Andosols have good aggregate stability and high permeability, which makes them resistant to erosion; however, under well hydrated conditions or subjected to drastic change they can be susceptible to damage (Meza-Pérez & Geissert-Kientz, 2006). Andosol cover about 0.7 % of the emerged land surface (IUSS Working Group WRB, 2016) worldwide; in Mexico, these soils cover less than 2 % of the continental surface (Secretaría de Medio Ambiente y Recursos Naturales [SEMARNAT], 2016) and, at the regional level, Michoacán accounts for 15.4 % (Bedolla-Ochoa, Bautista, Ihl, & Dubrovina, 2019).
In Mexico, avocado is grown in 28 states, but Michoacán is the main producer with 84 % of the national production with about 112 000 ha. The economic benefits of this crop make Michoacán a significant contributor to the regional economy, but the production model has a high impact on nature, due to land use change in forest lands and soil deterioration due to inadequate management practices and excessive use of agrochemicals. These factors point to Mexico as one of the countries most exposed to desertification; in particular, Michoacán has a severe problem of degradation and potential desertification due to crop expansion (Chávez-León et al., 2012; Villanueva & Zepeda, 2018).
The objective of this study was to evaluate the effect of the change from forest to agricultural land use on physical and hydrological properties of an Andosol. The hypothesis proposes that the change from a forest ecosystem to agricultural use, and agricultural management in relation to the application of organic and conventional management, causes alterations in soil properties.
Materials and Methods
Study area
The study area is located in the ejido of “Toreo El Alto” in Uruapan, Michoacán, located in the southern part of the Purépecha mountain range, between coordinates 19° 28’ 22.2’’ N and 102° 00’ 19.7’’ W. The area has an altitude of 1 890 m and belongs to the Balsas River region; the climate is temperate-humid with summer rains (Cw) (García, 2004), with an average annual temperature between 10 and 27 °C and average annual rainfall that exceeds 1 500 mm. The soil type is Andosol of volcanic origin (Alcalá, Ortiz, & Gutiérrez, 2001). Three main types of vegetation are identified: pine forest, pine-oak forest and mountain mesophyll forest (Bello & Madrigal, 1996).
In an Andosol soil, four experimental plots were selected: one for forest use and three for agricultural use. According to the classification of Rzedowski (2006), the forest plot has coniferous forest vegetation represented by a mixture of species of the genus Pinus and Quercus, which include: Pinus devoniana Lindley, Pinus pseudostrobus Brongn, P. lawsonii Roezl ex Gordon, P. leiophylla Schl. & Cham, Quercus rugosa Neé, Q. laurina Bonpl., Arbustus xalapensis Kunth and Fraxinus udheii (Wenz) Lingelsh, corresponding to an unmanaged and uneven-aged forest, with an average age of 80 years.
The agricultural plots correspond to two avocado orchards (Persea americana Mill var. Hass) with organic and conventional management and an organic macadamia nut orchard (Macadamia integrifolia Maiden & Betche) that uses the same management as the organic avocado orchard; the trees have an average age of 60 years in the three orchards.
Organic management consists of cattle compost, derived from cattle feces with a dose of 50 Mg∙ha-1 (N-P2O5-K2O [39-37-29] + S [18.4 %] + Mg [13.6] + Ca [74 %] and micronutrients), applied in partial shade for a period of three years. Conventional management includes chemical fertilizers such as copper sulfate pentahydrate (CuSO4·5H2O, 600 mL∙ha-1) used as a fungicide and bactericide for preventive use; CO (NH2)2 foliar fertilizer at a dose of 3 kg∙ha-1 (20-30-10) every two months; and an extended-release granules (15-00-00 + 26.6 [CaO]+ 0.3 [B]).
Sampling and soil analysis
A random sampling was carried out in each of the four experimental plots. Four composite samples were collected in each plot at two depths (0 to 20 and 20 to 40 cm). Samples were taken to the soil laboratory of the Faculty of Forestry Sciences of the Universidad Autónoma de Nuevo León (UANL) and were dried outdoors and sieved on a 2 mm sieve to determine texture and moisture content. Methods used for the analysis of soil physical and hydrological properties are shown in Table 1.
Properties | Method | Unit |
---|---|---|
Physical properties | ||
Texture1,2 | Method AS-09 of NOM-021-RECNAT-2000 (Secretaría de Medio Ambiente y Recursos Naturales [SEMARNAT], 2002). | % |
Bulk density1,2 | Gravimetric method (Woerner, 1989). | g∙cm-3 |
MRP1 | Hardness meter/penetrometer (soil hardness tester Yamanaka) (Medina-Guillén et al., 2017). | Kg∙cm-2 |
Porosity1,2 | Estimation by means of bulk density and particle size. | % |
Hydrological properties | ||
Moisture1 | Thermogravimetric technique, NOM-021-RECNAT-2000. | % |
Ks1 | Permeability analysis specified by Japanese Industrial Standards (Das, 2002). | cm∙s-1 |
Fb 1, F0 1 and Fa 1 | Double ring method (Zhang et al., 2017). | mm∙h-1 |
Available water1,2, PWP1,2 (15 bar) and FC1,2 (0.3 bar) | Pressure plate and membrane method (Soil Moisture Equipment Corp., Santa Barbara, CA) (Klute & Dirksen, 1986). | % |
MRP = mechanical penetration resistance, Ks = hydraulic conductivity, Fb = infiltration capacity, F0 = initial infiltration, Fa = cumulative infiltration, PWP = permanent wilting point and FC = field capacity. 1Depth from 0 to 20 cm and 2depth from 20 to 40 cm.
Statistical analysis
The variables were analyzed in a completely randomized experimental design with factorial arrangement to evaluate the effect of land use system and depth changes, and interaction. All variables were tested for normality and homogeneity (n = 32) and the following analyses were applied: Kruskal-Wallis nonparametric test for hydraulic conductivity (Ks), porosity, moisture and sand; one-factor ANOVA for infiltrations and mechanical resistance to penetration (MRP); and factorial ANOVA for available water (Aw), bulk density (BD), clay, silt, permanent wilting point (PWP) and field capacity (FC). Subsequently, a Tukey test was performed for the comparison of means. All variables were analyzed with Pearson's correlation coefficient at both depths (0 to 20 and 20 to 40 cm). Statistical analyses were performed using SPSS statistical package version 22.0 with a confidence level of P ≤ 0.05 (International Business Machines [IBM, 2013]).
Results
Land use effect on physical and hydrological properties at a depth of 0-20 cm
According to Table 2, ANOVA for surface variables such as infiltration (F0, Fb and Fa) and MRP showed significant differences (P ≤ 0.05) among land use scenarios.
Variables (0-20 cm) | SM | DF | MS | MSE | F | P | R2 | |
---|---|---|---|---|---|---|---|---|
F0(mm∙h-1) | 3 588 960 | 3 | 1 196 320 | 549 100 | 2.17** | 0.005 | 0.52 | |
Fb(mm∙h-1) | 491 703 | 3 | 163 901 | 12 451 | 13.16* | 0.042 | 0.86 | |
Fa(mm) | 2 566 634 | 3 | 855 544 | 100 865 | 8.48* | 0.014 | 0.8 | |
MRP (Kg∙cm-2) | 36.43 | 3 | 12.14 | 1.05 | 11.53** | 0.007 | 0.85 |
SM = sum of squares, DF = degrees of freedom, MS = mean square, MSE = mean square error, and R2 = adjusted r-squared. F0 = initial infiltration, Fb = basic infiltration and Fa = cumulative infiltration. **Highly significant differences (P ≤ 0.01); *significant differences (P ≤ 0.05).
Table 3 indicates that, according to Tukey's test, infiltration variables (F0, Fb and Fa) and MRP of agricultural land uses showed significant decreases (P ≤ 0.05) with respect to the forestry system. MRP of the forestry system had the lowest hardness in contrast to agricultural systems with avocado, which showed mean values of 4.5 Kg∙cm-2. Referring to the initial, basic, and cumulative infiltration variables, values were lower compared to those of the forestry system; on average, the three variables decreased 70, 43 and 40 % in conventional avocado, organic avocado and macadamia, respectively. Organic systems (avocado and macadamia) were similar, differing significantly from conventional management, which indicates that the use of organic amendments can be an alternative to improve these hydrological variables.
On the other hand, the Kruskal-Wallis test (Table 3) indicated that Ks and moisture content were significantly different (P ≤ 0.05) among land use scenarios. Ks showed average values from 0.002 to 0.060 cm∙s-1, the lowest value was for organic avocado and the highest value for macadamia. In general, moisture content ranged from 52.47 to 67.45 % in the macadamia crop and forest plot, respectively.
Land use | F0(mm∙h-1) | Fb(mm∙h-1) | Fa(mm) | Ks (cm∙s-1) | M (%) | MRP (Kg∙cm-2) |
---|---|---|---|---|---|---|
Forest use | 2 320 a | 858.08 a | 2 013.56 a | 0.020* | 67.45* | 2.05 c |
Conventional avocado | 600 c | 253.95 c | 618.50 c | 0.009* | 57.54* | 4.62 b |
Organic avocado | 1 470 b | 448.74 b | 1 103.50 b | 0.002* | 66.86* | 4.37 ab |
Macadamia | 1 560 b | 488.14 b | 1 194.40 b | 0.060* | 52.47* | 3.87 a |
Values with different letter (Tukey) and with asterisk (Kruskal-Wallis) represent significant differences (P ≤ 0.05) among land use. F0, Fb and Fa = initial, basic, and cumulative infiltration, respectively; Ks = hydraulic conductivity; M = moisture; MRP = mechanical resistance to penetration.
Effect of land use and depth changes on physical and hydrological properties
Table 4 shows the factorial ANOVA, where the soil use factor led to highly significant differences (P ≤ 0.01) in the variables BD and silt and clay content, and significant (P ≤ 0.05) for the variables FC and PWP, while Aw showed no differences. In the case of the factor depth only PWP and silt content had highly significant differences, while the rest of the variables showed no differences (BD, Aw, FC and clay). In the case of the factor interaction, only BD showed significant differences.
Variables | Land use | Depth | Interaction |
---|---|---|---|
F(3, 24) | F(1, 24) | F(3, 24) | |
BD (g∙cm-3) | 11.21** | 1.34ns | 0.474** |
0 | -0.289 | 0 | |
Aw (%) | 2.72ns | 2.51ns | 0.727ns |
-0.072 | -0.107 | -0.565 | |
FC (%) | 5.371* | 3.178ns | 1.871ns |
-0.008 | -0.126 | -0.199 | |
PWP (%) | 6.766* | 18.668** | 1.707ns |
-0.002 | 0 | -0.262 | |
Silt (%) | 23.45** | 17.38** | 0.339ns |
0 | -0.001 | -0.59 | |
Clay (%) | 12.13** | 0.286ns | 0.943ns |
0 | -0.601 | -0.453 |
**Highly significant differences (P ≤ 0.01); *significant differences (P ≤ 0.05). P < 0.05 is indicated in parentheses. BD = bulk density, Aw = available water, FC = field capacity, PWP = permanent wilting point.
Table 5 indicates that sand and porosity had significant differences (P ≤ 0.05) among land use, for both depths, based on the Kruskal-Wallis test. Porosity was relatively high in all land uses, the systems with the highest and lowest percentages at both depths were forest use with 78.9 % and organic avocado with 69.2 %. In general, sand contents for both depths were high in all land uses with averages of 55 %, except for conventional use which had 35 %.
The comparison of means of physical and hydrological properties by depth (Table 5) indicates that soil BD had low values, typical of Andosols, although forest use had the lowest values (0.55 and 0.65 g∙cm-3), and the highest values were recorded in organic avocado (0.81 and 0.80 g∙cm-3) at both depths.
FC ranged from 55 (macadamia) to 68 % (conventional avocado) for the first depth (0 to 20 cm) and ranged from 55 (organic avocado) to 62 % (forest use) for the second depth (20 to 40 cm). PWP ranged from 34 % (organic avocado) to 45 % (conventional avocado) for the first 20 cm, and at a depth of 20 to 40 cm, it decreased from 23 (macadamia) to 37 % (forest use). Aw ranged between 15 % for forestry system and 23 % for conventional use for the first depth (0 to 20 cm), which can be associated with high sand contents of forestry system; at a depth of 20 to 40 cm it varied from 16 to 26 % for the macadamia and organic avocado systems, respectively.
Silt values ranged from 40 to 60 % for the first depth (0 to 20 cm) and from 30 to 50 % for the second depth (20 to 40 cm). In the case of clay, the four land uses showed low contents, ranging from 6 to 11 % for 0 to 20 cm depth and from 5 to 9 % for 20 to 40 cm depth. The textural class for forest and organic land uses (organic avocado and macadamia) was sandy loam, and for conventional avocado was silt loam.
Variable | Land use | |||
---|---|---|---|---|
Forest use | Conventional avocado | Organic avocado | Macadamia | |
Depth from 0 to 20 cm | ||||
Sand (%) | 48.60* | 27.96* | 51.10* | 56.46* |
Silt (%) | 45.36 a | 60.81 b | 42.13 a | 37.36 c |
Clay (%) | 6.06 a | 11.23 b | 6.77 a | 6.18 a |
BD (g∙cm-3) | 0.55 a | 0.57 a | 0.81 b | 0.69 b |
Porosity (%) | 78.93* | 78.60* | 69.19* | 73.96* |
FC (%) | 58.14 a | 68.73 b | 56.94 a | 55.95 a |
PWP (%) | 42.78 a | 45.66 a | 34.49 b | 39.35 b |
Aw (%) | 15.35 a | 23.06 a | 22.44 a | 16.59 a |
Depth from 20 to 40 cm | ||||
Sand (%) | 56.28* | 37.05* | 59.14* | 67.82* |
Silt (%) | 37.66 a | 53.22 b | 32.91 a | 27.00 a |
Clay (%) | 6.06 a | 9.73 b | 7.95 a | 5.18 a |
BD (g∙cm-3) | 0.65 a | 0.58 a | 0.80 b | 0.73 b |
Porosity (%) | 75.34* | 77.80* | 69.48* | 72.42* |
FC (%) | 62.02 a | 61.60 b | 55.38 a | 40.12 a |
PWP (%) | 37.80 a | 36.89 a | 28.43 b | 23.27 c |
Aw (%) | 24.22 a | 24.71 a | 26.95 a | 16.84 a |
Values with different letter (Tukey) and with asterisk (Kruskal-Wallis) represent significant differences (P ≤ 0.05) among land uses by depth. BD = bulk density, FC = field capacity, PWP = permanent wilting point and Aw = Available water.
Correlation of physical and hydrological properties
According to Pearson's correlation results, at depth 0 to 20 cm, BD correlated negatively with porosity (-1), sand (-0.60), Ks (-0.82), F0 (-0.65), Fb (-0.83) and Fa (-0.76), and positively with clay (0.74), silt (0.53), FC (0.54) and MRP (0.83). There is a strong positive correlation among the variables Ks, porosity and BD. Infiltrations (F0, Fb, Fa) showed no correlation with any other variable except for the relationships between them (Table 6).
BD | Porosity | Sand | Clay | Silt | PWP | FC | Aw | Ks | MRP | F0 | Fb | Fa | M | |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
BD | 1 | |||||||||||||
Porosity | -1.000** | 1 | ||||||||||||
Sand | -0.609* | 0.608* | 1 | |||||||||||
Clay | 0.749** | -0.747** | -0.794** | 1 | ||||||||||
Silt | 0.535* | -0.535* | -0.988** | 0.692** | 1 | |||||||||
PWP | 0.283 | -0.284 | -0.388 | 0.15 | 0.423 | 1 | ||||||||
FC | 0.547* | -0.547* | -0.602* | 0.600* | 0.565* | 0.720** | 1 | |||||||
Aw | 0.437 | -0.436 | -0.388 | 0.670** | 0.293 | -0.177 | 0.556* | 1 | ||||||
Ks | -0.824** | 0.822** | 0.272 | -0.496 | -0.199 | 0.15 | -0.249 | -0.533* | 1 | |||||
MRP | 0.838** | -0.835** | -0.607* | 0.649** | 0.559* | 0.224 | 0.535* | 0.49 | -0.809** | 1 | ||||
F0 | -0.659* | 0.656* | 0.429 | -0.577* | -0.361 | 0.252 | -0.259 | -0.701* | 0.780** | -0.721** | 1 | |||
Fb | -0.832** | 0.830** | 0.452 | -0.662* | -0.366 | 0.247 | -0.282 | -0.732** | 0.924** | -0.781** | 0.921** | 1 | ||
Fa | -0.760** | 0.758** | 0.405 | -0.614* | -0.321 | 0.325 | -0.208 | -0.704* | 0.883** | -0.742** | 0.967** | 0.982** | 1 | |
M | -0.539 | 0.541 | -0.042 | -0.103 | 0.078 | 0.485 | 0.168 | -0.294 | 0.825** | -0.453 | 0.533 | 0.670* | 0.647* | 1 |
**highly significant correlations (P ≤ 0.01); *significant correlations (P ≤ 0.05). Bulk density (BD), porosity, sand, clay, silt, permanent wilting point (PWP), field capacity (FC), available water (Aw), hydraulic conductivity (Ks) and mechanical resistance to penetration (MRP): n = 16. Infiltrations F0, Fb, Fa and moisture (M): n = 12.
On the other hand, at depth 20 to 40 cm, BD correlated strongly with porosity, sand, clay and silt (-1, -0.90, 0.64 and 0.86, respectively). Similarly, porosity correlated strongly with sand, silt and clay with correlation coefficients of 0.90, -0.65 and -0.87, respectively.
BD | Porosity | Sand | Clay | Silt | PWP | FC | Aw | |
---|---|---|---|---|---|---|---|---|
BD | 1 | |||||||
Porosity | -1.000** | 1 | ||||||
Sand | -0.904** | 0.908** | 1 | |||||
Clay | 0.648** | -0.650** | -0.613* | 1 | ||||
Silt | 0.869** | -0.873** | -0.984** | 0.464 | 1 | |||
PWP | 0.419 | -0.426 | -0.640** | 0.469 | 0.613* | 1 | ||
FC | 0.373 | -0.383 | -0.560* | 0.520* | 0.512* | 0.835** | 1 | |
Aw | 0.19 | -0.2 | -0.275 | 0.387 | 0.222 | 0.362 | 0.815** | 1 |
**Highly significant correlations (P ≤ 0.01); *significant correlations (P ≤ 0.05). Bulk density (DA), permanent wilting point (PWP), field capacity (FC), available water (Aw).
Discussion
On a historical basis, deforestation in the state of Michoacán is associated with changes in land use due to agricultural activities under different organic and conventional management (Chávez-León et al., 2012), causing changes in the physical, hydrological, chemical and biological properties of the land. This was observed in the results of the present research, where the variables showed differences compared to the forestry system. This has a direct impact on soil ecosystem services and, consequently, accelerates degradation processes, where agricultural and livestock activities represent the main causes in Mexico (35 % of the land), followed by loss of vegetation cover, urban and industrial development (Etchevers, Saynes, Sánchez, & Roosevelt, 2016). Krasilnikov, Jiménez, Reyna, and García (2011) indicate Mexico needs to take soil conservation measures, given its dominant role in food production, because 44 % of the area is under agricultural activity and indicate that since 2010 there has been greater pressure on soil resources.
Acosta (2007) indicates that physical properties influence soil water dynamics, mainly texture, BD and porosity. Clay content allows a diagnosis of the physical-hydrological interaction of the soil, since increases in these particles lead to physical degradation, mainly modifying hydrological variables (Jiménez-Heredia, Martínez-Bravo, & Mancera-Rodríguez, 2010). This was the case of the conventional avocado plot, where infiltration variables were modified mainly by texture, because it had greater amount of clay and less infiltration.
According to Neris, Jiménez, Fuentes, Morillas, and Tejedor (2012), under natural or unmanaged conditions, Andosols are soils known for high infiltration rates, due to good structural development and aggregate stability; however, these properties are vulnerable to degradation by land use changes. The results of the present study are consistent with the above, where the highest infiltration rates were found under natural (forest) land use. On the other hand, infiltration differences between organic and conventional agricultural uses are marked by the contents of organic matter, which improves chemical, biological, physical and hydrological properties (Williams, Blanco-Canqui, Francis, & Galusha, 2017), being the application of organic amendments a common practice for organic avocado and macadamia.
Narro (1994) points out that Ks values between 0.003 and 0.005 cm∙s-1 indicate high permeability. According to the results of the present study, only the organic avocado plot is outside the range of this assessment (medium permeability), which may be associated with bulk density and porosity values influenced by the degree of mechanization during crop harvesting. On the other hand, Di Prima et al. (2018) mention that sandy loam-textured soils with agricultural management have positive and highly significant correlation between Ks and porosity, having a similar trend in this research. Larios-González, Salmerón-Miranda, and García-Centeno (2014) and Zúñiga, Buenaño, and Risco (2018) report that in Andosols soils it is common to find low BD values (<0.90 g∙cm-3) and high porosity (>65 %) (Ibáñez & Manríquez, 2011), which also coincides with the values estimated in this study.
MRP results show correlation with BD, porosity and texture, which influence hydrological properties, coinciding with Shah et al. (2017), who indicate that high BD and low soil moisture content have a direct impact on the increase of MRP. Parker (2007) and Vásquez, Gómez, Velázquez, Aldrete, and Fierros-González (2011) mention that Andosols with high MRP values exceed the ideal conditions for root growth and development, making them more vulnerable to physical deterioration.
Two textural classes were found (sandy loam and silt loam) in the granulometric analysis of Andosol, which seems to have influenced the infiltration variables; in the case of conventional avocado these variables could have been negatively affected by the texture, because a greater amount of clay was observed (negatively correlated with the infiltration variables) and lower sand content. Even so, Jiménez-Heredia et al. (2010) state that, in general, this type of soil has good conditions for the development of vegetation, because it has optimal drainage characteristics, water storage capacity and aeration, characteristics of volcanic soils (Paz & Sánchez, 2007).
Moisture, FC, PWP and Aw are in the range reported for Andosols with low BD (0.90 g∙cm-3) (Paz & Sánchez, 2007). Soil moisture is mainly affected by altered texture, BD and porosity (Anaya, Mendoza, Rivera, Páez, & Olivares-Martínez, 2016). Ibáñez and Manríquez (2011) indicate that Andosols show high values of Aw and PWP. In this study, the results showed PWP values between 34 and 45 % and for Aw from 15 to 23 %, showing positive correlation of PWP with silt (0.613) and negative with sand (-0.640) at a depth of 20 to 40 cm. Texture is an important variable in water retention at high tensions, due to the influence of soil particle size on fluid absorption and retention (Chicas, Vanegas, & García, 2014; Saxton & Rawls, 2006).
Conclusions
Land use change to agricultural use caused significant variations in the physical and hydrological properties of Andosol. The infiltration process was the most affected; infiltrations decreased 40 to 70 % in agricultural systems compared to forest use. Bulk density, porosity, mechanical resistance to penetration and clay defined the behavior of hydrological variables. Organic agriculture showed minor negative effects compared to conventional agriculture however, both are important for the economy of Michoacán. Therefore, the results of this research can be used to establish criteria and make sustainable decisions on soil resources. This study evidences the consequences of land use change on the behavior of physical and hydrological variables, showing that the conversion from forest use to agricultural use (organic or conventional) generates a significant negative effect on the behavior of soil properties.