INTRODUCTION
Sugarcane production continues to be the main agricultural activity in Belize, making a significant contribution to the agricultural sector, national income and foreign exchange earnings: 7.8% of gross domestic product (GDP), 10% of jobs and 6% of foreign exchange earnings (Statistical Institute of Belize, 2015). However, the current sugarcane management system is unsustainable (Chi et al., 2017), due to the inefficient fertilizer and herbicide application as well as inappropriate cultivation and harvesting practices which leads to soil degradation (Geissen et al., 2009), environmental deterioration, detrimental effects on soil and water quality, and nutritional imbalance (De León & Gonzalez, 2011).
Intensive agricultural systems such as sugarcane monoculture has also clearly negative impacts on biodiversity conservation (Malézieux et al., 2009), including soil macroinvertebrates populations (Huerta et al., 2006). Thus, it reduces taxonomic richness, density and biomass, indicating high disturbance level of the edaphic environment due to land use intensity, like tillage or fertilizer (Gardi et al., 2002; Cabrera et al., 2011). For instance, in Cuba, Cabrera et al. (2011) report higher densities and biomass of soil macroinvertebrates for secondary forests in comparison to sugarcane monoculture; and Paoletti (1999) and Siqueira et al. (2016) found that earthworms, myriapods and other invertebrates are sensitive to toxic effects of agro-chemicals. Moreover, soil types can also influence local diversity of edaphic macroinvertebrates (Mathieu et al., 2005). In Belize, the sugarcane industry ignores the impact that agricultural inputs and management could have on the environment and more specifically on soil macroinvertebrates.
During recent years renewed scientific attention has propelled evaluation of soil quality bioindicators due to concerns of soil degradation and the need for a sustainable management of soils (Shukla et al. 2006; Barrios, 2007; López et al., 2011), through evaluation of relation between soil land use and soil macroinvertebrates (Barros et al., 2002; Mathieu et al., 2005; Velasquez et al., 2007), allowing an ecological dynamic balance and biodiversity (Lang-Ovalle et al., 2010). Siqueira et al. (2016) indicates that soil macrofauna may be used as ecological indicators of agriculture practices since they are very sensitive to both chemical and physical soil parameters. Thus, in the search of a sustainable agriculture in sugarcane production soil macroinvertebrates are required to be included as a soil quality indicator.
Diverse processes occurring in the soil are mediated by organisms that lives on it, such as macroinvertebrates (Lavelle et al., 2006), which includes a great variety of edaphic organisms larger than 2 mm in size (Baretta, 2006; Bardgett & van der Putten, 2014). In this regard, macroinvertebrates are important in soil ecosystems for their capacity to decompose organic matter, forming complex food webs essential to a functioning ecosystem (Barnes et al., 2016), as well as supporting most agriculture production systems (production services) through beneficial services that they mediate: soil formation, nutrient cycling and primary production (Lavelle et al., 2006). Moreover, they intervene in soil processes: infiltration, aeration; as well as aggregation, texture (Siqueira et al., 2016), chemical and nutritional properties, and can serve as soil fertility quality bioindicators (Mummey et al., 2002). On the other hand, soil macroinvertebrate communities in each habitat are determined by management practices (physical perturbation, residue distribution, amendments and vegetation) (Bardier, 2005) since their response to changes occurs quickly (Correia & Oliveira, 2000). Hence, sustainable soil land use and management of sugarcane fields can improve or alter the abundance, diversity and richness of these organism.
Filter mud is a by-product from sugarcane industries and characterized as a soft, spongy, amorphous, and dark brown to brownish material (Ghulam et al., 2012). Generally, for every 100 t of sugarcane crushed, about 3 t of filter mud is produced as a by-product (Gupta et al., 2011). During the 2015-2016 crop cycle, the sugar industry in Belize produced 38,775.36 t of filter mud (Sugar Cane Production Committee (SCPC), 2015). Likewise, vinasse, obtained from fermentation and distillation of alcohol, which is acidic in nature, dark brown in color with unpleasant odour, possess a high chemical oxygen demand (COD), and biological oxygen demand (BOD) (Dotaniya et al., 2016). In Belize vinasse is produced by two small distilleries which operate seasonally, depending on consumption, but on average, about 10-18 liters of vinasse is obtained per liter of produced alcohol, giving a total production for the industry of 4,400 m3 per year (personal communication, Belize Sugar Industries, 2017). They can improve soil chemical, physical and biological properties, thus enhancing crop quality and yield (Dotaniya et al., 2016, Lang-Ovalle et al., 2010). Also, they can favor humidity retention, microorganism’s activation (Neugebaver et al., 1992), as bio-fertilizers (Rodríguez, 1994; Vandevivere & Ramírez, 1995), improve soil texture and aggregation (Mojica et al., 1994) and contribute to the conservation of natural resources by recycling carbon and mineral elements (Prado et al., 2013).
Thus, direct effects of organic waste from the sugarcane industry on soil chemical properties are increased nitrogen, phosphorus, and calcium concentrations, increased cation exchange capacity (CEC), hence playing a fundamental role in agricultural production, in the maintenance of soil fertility, and as a soil conditioner (Rossetto et al., 2008). These important aspects have been widely discussed among researchers (Barnes et al., 2016), thus Jones et al. (1994) associated termites to soil carbon cycle and soil fertility level; while Subbler et al. (1998), Lang-Ovalle et al. (2010) and Blouin et al. (2013) indicated that earthworms improve soil structure and porosity of agriculture soils.
Regardless the indispensable role of soil macroinvertebrates in soil management, few studies related to sugarcane production have been conducted (Pasqualin et al., 2012; Benazzi et al., 2013). Therefore, the aim of the present study was to determine the abundance and diversity of soil macroinvertebrates in sugarcane plantations under different managements with chemical fertilizers and organic amendments, under the following hypothesis “the chemical fertilizers used in the production of sugarcane have a negative impact on the abundance and diversity of macroinvertebrates”, and predictions: “the abundance and diversity of soil macroinvertebrates in the soil of sugarcane plantations treated with organic amendments (vinasse and filter mud) will be higher than in the soils treated with agrochemicals”.
MATERIALS AND METHODS
Description of study sites. The study was conducted in northern Belize (Figure 1; 18°00’ and 18°30’ N and 88°34’ and 88°39’ W) with prevailing climate according to Köppen classification system and modified by Garcia (2004) as warm sub-humid (Aw), with rains in summer and a dry season in winter. Specifically, the area falls within the Ax' (w1) climate subtype, intermediate between the wettest and the driest of the sub-humid classes (Garcia, 2004), with a mean annual long-term rainfall and temperature of 1,053.33 mm and 26oC, respectively.
The study lasted two crop cycles (12 months cycle), from May 2013 to May 2015, and on two sites belonging to American Sugar Refinery/Belize Sugar Industries (ASR/BSI), with different soil types; that are superficial red young Cambisol and Fluvisol (WRB, 2014). In addition, one plot of secondary forest, (undisturbed for more than 15 years), on a Cambisol (WRB, 2014) was also selected as an extra control. It is important to mention that the secondary forest was managed in the past under conventional sugarcane, including application of chemical fertilizer.
Fluvisols are very young soils with weak horizon differentiation and may contain high natural fertility which allows cultivation of a wide range of crops from dryland (WRB, 2014); however flood control, drainage or irrigation are required (Driessen & Dudal, 1989); it present a dark coarse-textured topsoil (<15% clay), grey medium or fine textured subsoil over weathering limestone and is locally known in Belize as ‘Pixoy’ (King et al., 1992). On the other hand, Cambisols are young soils, this means, they are not highly weathered, what is verified with beginning horizon differentiation through changes in color, structure or texture. They are derived from material with moderate mineral reserves found over calcareous material, medium OM content (De León & Gonzales, 2011) and are locally known as ‘Lazaro clay’ (King et al., 1992). Cambisols are commonly used in a wide variety of agricultural crops (Driessen & Dudal, 1989).
Establishment and management of experimental fields. The sugarcane crops were established during the spring planting period, that is late May to June, and harvested toward the end of the following two crop cycles at 12 months. Land and soil preparation involved removing vegetative material and preparing a tilth adequate for cane emergence and growth. This activity involved five operations; sub-soiling, ploughing (18, 34" discs), harrowing (24 discs), furrowing and bed formation, after which the cane was planted manually. The present investigation was part of a larger fertilization research with crop-yield purposes, where 11 doses of chemical fertilizer, 11 doses of chemical fertilizer combine with vinasse, and 11 doses of chemical fertilizer combine with filter mud were tested. In the present study, sugarcane plots were selected under three different management systems; application of chemical fertilizer (11 different doses of macro elements: N-P-K: 0-20-0; 0-20-100; 0-0-100; 40-10-150; 40-20-50; 40-30-50; 80-0-100; 80-0-0; 80-20-100; 80-20-0; and 120-10-50), vinasse (10 m3 ha-1), and filter mud (50 m3 ha-1). A control without any chemical fertilizer nor organic amendments was also established. Furthermore, a secondary forest, located close to the sugarcane crops, was also sampled. Both vinasse and filter mud doses used in the study resulted from recommendations of previous works from ASR/BSI consultancies.
In the sugarcane plots a randomized complete block design with 4 repetitions were used giving a total of 56 experimental units (EU) per site and cycle. 44 with chemical fertilizer (11 different doses (see above) X 4 repetitions), 4 EU for vinasse and filter mud respectively (1 dose X 4 repetitions) and 4 EU for control. Each EU consisted of 4 rows of 10 m long and 6 meters wide, being the 2 central rows the sampling unit (SU). The net area per EU was 60 m2; giving a gross area of 3,360 m2 and a net area of 1,680 m2. The distance of ridges between furrows was 2 m, 1.5 m row width and 0.7 m walkway.
Filter mud is rich in N (1.67%), P (2.4%) and OM (52.5%) (Aloma, 1973; Prado & Jaramillo, 2006; Melchor et al., 2008; CATS, 2012; Cifuentes et al., 2011). For application it was transported on a fertilizer spreader hauled by tractor and applied directly on top of the setts of the sugarcane study plots. In our study, the filter applied mud came from ASR/BSI sugar factory.
Vinasse is rich in K, Ca (0.05%) and Mg (0.0129) (Dotaniya et al., 2016). It was transported on cistern type container hauled and operated by a tractor and applied directly on top of the sets.
Macroinvertebrates and soil sampling methodology. In 2013 and 2014 at the end of the rainy season when macroinvertebrates activity is at its highest, since during low humidity earthworms are dormant (Huerta et al., 2008), according to the Soil Biology and Fertility (TSBF) method (Anderson & Ingram, 1993), a monolith of 25 x 25 x 30 cm was dug on the SU in each EU; whereas, in secondary forest three 40 m line transects were established at 60 m from each other for a total of 5 samples (monoliths) per transect (8 m between each sample) giving a total of 15 monoliths per sugarcane cycle (2013-2014). A total of 254 monoliths; 56 monoliths at two sites, two per cycle (56 plots/site x 2 site x 2 cycles) in the sugar cane plantation and 30 monoliths for secondary forest (3 transect x 5 samples per transect x 2 cycle) were performed. From the soil of each monolith soil macroinvertebrates (organisms > 2 mm in size) where collected by hand sorting (meticulously checked on site). They were placed in plastic containers with 96% alcohol for earthworms and 76% alcohol for the rest of macroinvertebrates (insects, arachnids, gastropods) to carry to the laboratory. Each container was labelled to know from which monolith and treatment they come from.
Soil macroinvertebrates handling. The soil macroinvertebrates were characterized by their abundance, biomass and diversity. The counting, weighing and classification were performed at the Sugar Industry Research and Development Institute (SIRDI) laboratory in Belize. The method for quantifying soil macroinvertebrates was through direct observation (Williams, 1987), they were identified with the aid of a stereoscope and weighed in an electronic balance (0.01 g precision). For macroinvertebrates from the Arthropoda and Mollusca Phylum, identification was done up to morphospecies using Gibb and Oseto (2006) key and for earthworms the James (1994) key.
The soil biological properties investigated in this study were: macroinvertebrate morphospecies, macroinvertebrate biomass (g.m2), macroinvertebrate density (ind.m2), ecological index for tropical soil quality based on macrofauna (hereafter soil quality index) (Huerta et al., 2009) and a conversion of the Shannon diversity index (hereafter Shannon_exp index). The presence of soil invertebrates can be used as an indicator of soil quality (Klemens et al., 2003; Ruiz-Camacho, 2004; Lavelle et al., 2006), therefore, under this premises the soil quality index based on macrofauna is an index based on the presence/absence of different groups of soil macrofauna to determine the ecological quality of soils (Huerta et al., 2009). This index considered the number of earthworm species, number of native earthworm species, presence of epigeic earthworms, and presence/absence of litter decomposers from the Order Diplopoda, Isopoda, Coleoptera, Heteroptera and Orthoptera, which integrate the structure of macrodecomposer community. This information is combined with the biomass of earthworms and is passed through a set of rules that determines the soil quality index (Fig. 2) (Huerta et al., 2009). The Shannon diversity index is a commonly used method of measuring species diversity (Stiling, 1999). It assumes that all species from the community are represented in a randomly collected sample. The formula used is:
Soil characterization of plant macronutrients. In each monolith, both on sugarcane fields and secondary forest, 500 g of soil was collected (Bautista et al., 2011) with the purpose of characterize the status of plant macronutrients (N, P, and K) in the soils under study. Soil samples were placed in labelled plastic bags, and dried in the shade, sieved and sent to soil laboratory for analysis. In total, 127 loose samples (56 plots/site x 2 sites) for sugarcane fields, and 15 loose samples for secondary forest (3 transect x 5 samples per transect) were collected and processed. The methods used to determine the plant macronutrients were: Total N by micro-Kjeldahl; K by extraction with ammonium acetate at pH 7 and P by the Olsen method (Olsen & Summers, 1982); according to NOM-021-RECNAT-2000 (Norma Oficial Mexicana (NOM), 2002).
Plant variables. As reference of the response to the different types of fertilizers for crop yield, the following plant variables were measured: yield (t ha-1), number of stools (stool/m2), (number of stems (stem/m), number of active leaves (leaf/plant), stem length (cm) and stem diameter (cm).
Statistical Analysis of data. In order to test differences among the different management systems (chemical fertilizer, vinnase, filter mud, and secondary forest), on soil macroinvertebrates parameters (density, biomass, number of morphospecies), ecological index for tropical soil quality, Shannon_exp index and; plant variables, one-way ANOVA followed by post hoc analysis means test was done (Tukey, p < 0.05) after testing normality. It was as well tested differences between soil type and cycle, on the same parameters, as well on the macronutirients (N, P, K), this only between soil types, using t-test analysis (Minitab® version 17).
To identify relatedness among soil macroinvertebrates variables with the soil quality index, Shannon_exp index and sugarcane yield a principal component analysis (PCA) was performed including treatments (management systems), soil type and sampling cycles.
RESULTS
In the present study a total of 7,634 macroinvertebrate individuals, from the 254 monoliths, were identified and grouped in morphospecies (Table 1). The edaphic macroinvertebrate community was characterized by the presence of different annelids, arthropods and mollusc groups. From the total, annelids were the Phylum with the highest number of individuals with 4,461 (58%).
Phylum | Class | Order | Family | Gender | Species | Sugarcane | Secondary Forest | |||
---|---|---|---|---|---|---|---|---|---|---|
# of msp* | # of ind* | # of msp* | # of ind* | |||||||
Annelida | Clittelata | Oligochaeta | Acanthodrilidae | Dichogaster | D. saliens | NA | 2419 | NA | 9 | |
Dichogaster | D. balaui | NA | 39 | NA | 4 | |||||
Pontoscolecidae | Pontoscolex | P. corethrurus | NA | 1989 | NA | 1 | ||||
Arthropoda | Insecta | Blattodea | NI | NI | NI | 2 | 357 | 0 | 3 | |
Colembola | NI | NI | NI | 2 | 5 | 0 | 2 | |||
Coleoptera | NI | NI | NI | 5 | 499 | 4 | 31 | |||
Dermaptera | NI | NI | NI | 4 | 59 | 2 | 7 | |||
Diptera | NI | NI | NI | 1 | 40 | 0 | 0 | |||
Hemiptera | NI | NI | NI | 1 | 4 | 0 | 1 | |||
Homoptera | NI | NI | NI | 2 | 21 | 0 | 0 | |||
Hymenoptera | Formicidae | NI | NI | 6 | 1155 | 2 | 30 | |||
Isoptera | NI | NI | NI | 1 | 6 | 1 | 16 | |||
Orthoptera | NI | NI | NI | 4 | 117 | 1 | 7 | |||
Thysanoptera | NI | NI | NI | 2 | 6 | 0 | 0 | |||
Chelicerata | Anarea | NI | NI | NI | 5 | 268 | 2 | 9 | ||
Myriapoda | Chilopoda | NI | NI | NI | 4 | 294 | 1 | 6 | ||
Diplopoda | NI | NI | NI | 4 | 146 | 1 | 8 | |||
Mollusca | Gasteropa | NI | NI | NI | NI | 2 | 64 | 1 | 12 |
*per Taxonomic group
msp= morphospecies
The identified morphospecies belonged to 15 different Taxa including Aranea, Blattodea, Chilopoda, Colembola, Coleoptera, Dermaptera, Diptera, Diplopoda, Gasteropoda, Hemiptera, Homoptera, Hymenoptera, Isoptera, Orthoptera, Thysanoptera, and one not identified from the Phylum Mollusca (Table 1).
Arthropoda Phylum showed higher abundance for Hymenoptera in sugarcane fields with 38%, while Coleoptera 16.4%, Blattodea 11.7%, Chilopoda 9.7% and Araneae 8.8%; only the more abundant are enlisted here (n= 576). On the other hand, on secondary forest, Coleoptera showed higher abundance with 23.5%, while Hymenoptera scored 22.7%, Isoptera 12.1%, Gasteropoda 9.1% and Aranea 6.8% (n= 60) (Table 1).
Three species of annelids from the Subclass Oligochaeta were determined Dichogaster saliens, Dichogaster bolaui and Pontoscolex corethrurus; the most abundant was Dichogaster saliens, which along with Dichogaster bolaui are both epigeic and exotic; while Pontoscolex corethrurus endogeic and exotic (Table 1).
Results showed no effect of type of fertilizers on soil macroinvertebrate parameters, nor on plant variables (Table 2, 3). Contrary to what was expected, abundance and diversity of macroinvertebrates were not higher in organic treated soils, neither in soils under secondary forest. However, there were significant differences between soil type in most parameters and in cycles in some parameters (Table 2, 3).
Treatment* (16) | Arthropoda Mollusca Biomass (g.m2) | Earthworm Biomass (g.m2) | Arthropoda Mollusca Density (g.m2) | Earthworm Density (g.m2) | Number Morpho spp | Shannon_ exp exp(H) | Soil quality index | Yield (t ha-1) |
---|---|---|---|---|---|---|---|---|
Control | 1.96±2.44 | 1.68±4.72 | 128.0±118.7 | 69.0±82.4 | 3.00±1.46 | 51.1±100.0 | 0.02±0.075 | 104.1±28.8 |
0-20-0 | 0.70±0.72 | 1.00±1.69 | 97.0±58.6 | 70.0±97.0 | 3.31±1.89 | 230.0±736.0 | 0.10±0.027 | 95.8±25.6 |
0-20-100 | 2.76±3.91 | 1.60±4.07 | 162.0±240.5 | 53.0±59.5 | 3.06±1.57 | 61.4±102.5 | 0.08±0.022 | 99.5±29.8 |
0-0-100 | 1.59±2.38 | 1.86±2.93 | 106.0±81.6 | 76.0±67.3 | 3.19±1.47 | 63.4±102.0 | 0.00±0.000 | 102.7±34.4 |
40-10-150 | 0.91±1.59 | 1.07±2.83 | 68.0±50.8 | 68.0±84.4 | 2.50±1.67 | 87.0±271.5 | 0.01±0.023 | 106.2±33.7 |
40-20-50 | 4.18±7.99 | 0.61±1.08 | 103. ±101.1 | 115.0±150.8 | 3.12±1.82 | 232.0±740.0 | 0.01±0.055 | 97.6±29.8 |
40-30-50 | 3.47±8.61 | 2.76±5.63 | 79.0±52.1 | 94.0±113.4 | 3.06±1.53 | 60.0±102.7 | 0.04±0.175 | 108.8±33.4 |
80-0-100 | 0.74±1.97 | 2.37±5.51 | 78.0±77.9 | 115.0±198.8 | 2.37±1.20 | 23.3±37.3 | 0.02±0.076 | 102.9±20.2 |
80-0-0 | 1.68±2.28 | 3.56±8.18 | 71.0±44.9 | 68.0±70.5 | 2.75±1.39 | 39.1±55.7 | 0.07±0.208 | 113.0±28.2 |
80-20-100 | 1.10±1.68 | 1.09±2.43 | 94.0±84.8 | 81.0±131.4 | 2.81±1.56 | 53.2±100.8 | 0.08±0.125 | 97.1±28.1 |
80-20-0 | 1.11±2.42 | 4.08±9.38 | 68.0±53.7 | 72.0±61.3 | 2.37±1.82 | 49.1±102.0 | 0.00±0.000 | 103.9±41.7 |
120-10-50 | 0.99±1.59 | 1.09±1.97 | 100.0±84.6 | 119.0±283.5 | 2.87±1.20 | 28.4±21.9 | 0.01±0.057 | 104.4±31.4 |
Filter mud | 1.96±2.83 | 1.78±5.15 | 83.0±71.8 | 433.0±1532.0 | 2.44±1.50 | 28.5±38.6 | 0.06±0.140 | 108.5±27.4 |
Vinasse | 3.68±6.42 | 0.87±1.62 | 101.0±93.8 | 229.0±667.0 | 2.87±1.09 | 34.1±47.2 | 0.01±0.040 | 111.4±34.2 |
Secondary Forest (30) | 1.50±3.47 | 0.52±2.63 | 82.1±58.5 | 4.2±12.5 | 3.07±1.76 | 80.5±134.1 | 0.02±0.040 | |
Soil | ||||||||
Fluvisol (112) | 2.7±4.49a | 2.52±5.32a | 129.0±123.0a | 205.0±639.0a | 3.18±1.50a | 82.0±359.0a | 0.05±0.160a | 98.4±20.9b |
Cambisol (142) | 1.21±3.60b | 0.99±3.69b | 66.7±51.7b | 26.1±39.8b | 2.62±1.54b | 67.0±129.0a | 0.02±0.100b | 109.7±36.8a |
Cycle | ||||||||
1 (127) | 1.81±3.24a | 1.83±3.86a | 99.0±115.0a | 172.0±605.0a | 2.70±1.54a | 71.0±287.0a | 0.02±0.100b | 82.5±13.7b |
2 (127) | 1.93±4.78a | 1.50±5.13a | 88.9±70.9b | 38.2±47.5b | 3.03±1.53a | 79.0±276.0a | 0.04±0.160a | 125.5±27.0a |
* No significant differences among treatments (p < 0.05). For Soil and Cycle, in columns means with different lowercase letters are significantly different according to t - test (p < 0.05).
Factor | Number of stool (stool/m2) | Number of stem (stalk/m) | Number of active leaves (leaf/plant) | Stem length (cm) | Stem diameter (cm) |
---|---|---|---|---|---|
Soil | |||||
Fluvisol (112) | 2.89±0.73 a | 14.62±0.11 a | 5.44±0.02 b | 195.20±0.47 a | 7.54±0.01 a |
Cambisol (112) | 1.57±0.58 b | 15.07±0.12 a | 5.82±0.04 a | 171.72±0.48 b | 7.58±0.05 a |
Cycle | |||||
1 (112) | 2.89±0.73 a | 14.82±0.11 a | 5.63±0.03 a | 183.46±0.15 a | 7.56±0.01 a |
2 (112) | 2.89±0.74 a | 14.85±0.12 a | 5.61±0.05 a | 183.03±0.25 a | 7.56±0.05 a |
For Soil and Cycle, in columns, means with different lowercase letters are significantly different according to t - test (p < 0.05).
Soil macroinvertebrate morphospecies, biomass and density. The number of morphospecies among type of management systems range from 2.37 to 3.19, both mean values belonging to chemical fertilizer formulas, while filter mud, vinasse, control and secondary forest had 2.44, 2.87, 3.00 and 3.07 mean values, respectively. The number of morphospecies for Cambisol (2.62) was significantly different to Fluvisol (3.18) and there were no significant differences between cycles; cycle 1, 2.70 and cycle 2, 3.03 (Table 2).
The macroinvertebrates belonging to Arthropoda and Mollusca Phylum (hereafter Arthropoda-Mollusca) biomass ranged from 0.70 to 4.18 g.m2 and from 0.52 to 4.08 g.m2 for earthworm biomass. For Arthropoda-Mollusca, the lowest and highest values belonged to chemical fertilizer formulas, the vinasse mean value was 3.68 g.m2, while for both filter mud and control 1.96 g.m2, and 1.50 g.m2 for secondary forest. Earthworm biomass varied as follow; chemical fertilizer 4.08 g.m2, vinasse 0.87 g.m2, filter mud 1.78 g.m2, secondary forest 0.52 g.m2 and control 1.68 g.m2. Fluvisol had significantly higher biomass than Cambisol in both Arthropoda-Mollusca and earthworm biomass, 2.68-1.22 g.m2 and 2.52-0.99 g.m2, respectively. And, for cycle 1 and 2 the values were, for Arthropoda-Mollusca biomass; 1.81 and 1.93 g.m2 and earthworm biomass 1.83 and 1.50 g.m2, respectively, although no significant differences were found (Table 2).
For Arthropoda-Mollusca density, chemical fertilizer formulas had the lowest and the highest mean values 0.70 and 1.62 ind.m2, while vinasse, filter mud, secondary forest and control had 101.00, 83.00, 81.1 and 128 ind.m2, respectively. Earthworm density ranged from 4.27 to 433.00 ind.m2, the lowest value belonged to secondary forest and the largest to filter mud. For vinasse and control the densities were 229.00 and 69.00 ind.m2, respectively, and for chemical fertilizer varied from 53.00 to 119.00 ind.m2. Fluvisol, with 129.00 and 205.00 ind.m2 was significantly higher than Cambisols 66.70 and 26.10 ind.m2 on Arthropoda-Mollusca and earthworm densities, respectively, while cycle 1 was not significantly different to cycle 2 on Arthropoda-Mollusca density, 99.00-88.90 ind.m2, but significantly different on earthworm density, 172.00-38.20 ind.m2, respectively.
Soil quality index and Shannon_exp diversity index. There were not significant differences among management systems on the soil quality index. In general, all values were extremely low, from 0.001 to 0.10. The largest mean value was for chemical fertilizer formula (0.10), while vinasse, filter mud, secondary forest and control had 0.01, 0.06, 0.02 and 0.02, respectively; Fluvisol (0.05) had significantly higher soil quality index than Cambisol (0.02) (p < 0.001); and no differences were observed between cycles (Table 2). The Shannon_exp index (exp H) ranged from 28.5 to 232.0. The largest values were for chemical fertilizer formulas, while for vinasse, filter mud, secondary forest and control, was 34.1, 28.5, 80.5 and 51.1 respectively. No significant differences between soil types and cycles (Table 2).
In the principal component analysis (Figure 3), some soil macroinvertebrate parameters showed correlation to Shannon_exp index and soil quality index: Shannon_exp index with number of morphospecies (r2= 0.64, p = 0.001), and soil quality index with earthworm biomass (r2= 0.32, p = 0.001). However, when correlations were run separately by management system, it could be identified the high affinity of some management systems with some macroinvertebrate parameters. For example, filter mud with earthworm biomass and density. For filter mud the correlation between soil quality index and earthworm biomass and density was r2= 0.66, p < 0.001 and r2= 0.70, p < 0.00, respectively. For control the correlation between Shannon_exp index with number of morphospecies was r2= 0.86, p < 0.001, and between soil quality index with earthworm biomass r2= 0.98, p < 0.001. And for secondary forest the correlation between soil quality index with earthworm biomass was r2= 0.99, p < 0.001, and between Shannon_exp index with Arthropoda-Mollusca density r2= 0.52, p < 0.001.
Plant macronutrients, yield and plant variables. In general, total nitrogen content was low ranging from 0.13-0.17%, while both phosphorus and potassium content were high with value ranging from 21.04-61.93 mg kg⁻¹ and 0.45-0.68 cmol kg-1, respectively. Cambisol had significantly higher total nitrogen content (0.20%) than the site with Fluvisol soil (0.12%) (p < 0.001). With regards to phosphorous content, Cambisol (54.20 mg kg-1) did not differ to Fluvisol (49.90 mg kg-1). Similarly, to total nitrogen content, Cambisol had higher potassium content than Fluvisol, with 0.99 cmol(+) kg-1 soil compare to 0.35 cmol(+) kg-1 soil, respectively (p < 0.001).
No significant differences were observed among type of fertilizers on yield; yield on chemical fertilizer plots varied between 95.75 and 113.01 t ha-1, while vinasse, filter mud and control had a yield of 111.35, 108.53 and 104.13 t ha-1, respectively (Table 2). Cambisol, as with the N and K macroelements, had significantly higher yield (109.70 t ha-1) than Fluvisol (98.40 t ha-1). And cycle 2 (125.50 t ha-1), had significantly higher yield than cycle 1 (82.50 t ha-1) (Table 3).
Moreover, no differences was observed in all plant response variables with regards to type of fertilizers; however, the plants grown in Fluvisol had higher number of stool (2.89 stool m2) and stem length (195.20 cm) than plants in Cambisol (1.57 stool m2; 171.72 cm) respectively, while plants grown in Cambisol had significantly higher number of active leaves (5.82 leaf/plant) than plants in Fluvisol (5.44 leaf/plant) (Table 3). No significant differences were observed between cycles in all plant response variables (Table 3).
DISCUSSION
Soil macroinvertebrate morphospecies, biomass and density. Prolonged cultivation of sugarcane at the study area using agronomic practices with different conventional tillage, fertilization and chemical application in mono-cropping system (Chi et al., 2017) seems to have dramatically affected soil ecosystem, specifically on soil macroinvertebrates diversity and abundance. Our study area, including secondary forest harboured low levels of morphospecies, macroinvertebrate mass and density (Table 2). The practice of burning sugarcane fields before and after harvesting destroys organic matter, main source of energy for macroinvertebrates. In an investigation conducted by Chi et al. (2017) in the study area, found that the organic matter of the soil decreased from 35 to 70% with years (13 to 25) of sugarcane cultivation. In addition, fire kills soil organisms in the immediate zone, further reducing population biomass and density of soil macroinvertebrates (Ramírez et al., 2004, Ramírez et al., 2012).
Jiang et al. (2015) showed that crop residues and organic manure can substantially increase the abundance of macroinvertebrate individuals and communities, with similar positive impact on yields (Gilbert et al., 2008), organic residues from industry are beneficial since they improve soil water conservation and provide a food source for soil macroinvertebrates survival (Deibert & Utter, 1994). In the present study, filter mud (433.0 ind.m2) and vinasse (229.0 ind.m2) had the largest and the second largest mean values on earthworm density respectively. In addition, vinasse had the second largest mean value on Arthropoda-Mollusca biomass (3.68 g.m2), however, no significant difference was obtained. Thus, the benefits of organic manure in our study cannot been proven, two years of applying soil organic amendments may be too little to offset the impact of long-lasting sugarcane cultivation. Therefore, studies are required to observe the effect of soil organic amendments on soil invertebrates in sugarcane production in longer observation time (i.e 3-5 years) (Hernández-Melchor et al., 2008).
Lavelle et al. (1994) indicated that in the tropics, earthworms are generally predominant in macroinvertebrate biomass and can have important contribution on soils and crop productivity (Brown et al., 1999), notwithstanding in the present study earthworms group represented the highest total number of macroinvertebrates (58%), while Arthropoda-Mollusca (42%), their abundance and diversity were very low. The three species of annelids identified in the present study, Dichogaster saliens, Dichogaster bolaui and Pontoscolex corethrurus are exotic. On the other hand, exotic earthworms have a high tolerance range to extreme management conditions (Fragoso, 2001), since they are influenced by the type of vegetation, management conditions and ecological status of the system. In addition, species richness and presence of certain groups can reveal important information of soil quality (Lavelle, 1996; Klemens, 2003), further indicating perturbation of these ecosystems by conventional management practices (Lavelle & Spain, 2001). Studies with earthworms on intensified agriculture demonstrate that on annual crops earthworm’s specific richness and biomass tend to decrease, with native communities replaced by exotic earthworm (Decaëns & Jiménez, 2002). The presence of these species possibly follows better adaptation process in these managed ecosystems assuring their survival, since some species of soil macroinvertebrates are comparatively less affected and can respond positively (Lagerlof & Andren, 1988). For instance, Pontoscolex corethrurus, a geophagous earthworm, is capable of inducing the formation of dense aggregates through the production of casts (Hallaire et al., 2000). This species tend to proliferate very well in disturbed soils due to different macro factors (precipitation, temperature, solar radiation, high reproductive rate) and micro environmental (pH, Organic matter, N, P, S among other), as well as to the diverse management techniques of anthropogenic environments (Fragoso, 2001; Lavelle & Spain, 2001; Edwards & Bohlen, 1996).
Many biological communities present a high number of few abundant species under natural conditions (De Ruiter et al., 2002). This concept applies to many communities of soil organisms in tropical ecosystems as indicated by studies on macroinvertebrates (Portilho et al., 2011). In this regard, Siqueira et al. (2016) in a similar study to ours on sugarcane monoculture and two natural vegetation types (Sandbank and Atlantic Forest) in Brazil obtained similar diversity, respectively, suggesting that at the beginning of sugarcane growth, soil macrofauna communities undergo an initial selection and only individuals better adapted to sugarcane management practices and climatic conditions persist (Portilho et al., 2011).
According to Lavelle and Spain (2001) the main groups of soil macrofauna in terms of their abundance and importance of their activities in soil are earthworms, termites, ants, millipedes, beetle larvae and adults. In our study, the macroinvertebrate community was found to be almost diverse in terms of orders, more representatives by earthworms, Hymenoptera, Coleoptera, Blattodea, Chilopoda, Araneae and other minors taxa. None the less, ants were the most abundant Arthropoda (1,185 individuals, 37.30%) in all systems studied since they generally present rapid colony restructuring, wide foraging area, and social organization (Santos et al., 2017). Even with the disruption caused by the use of fire in sugarcane harvesting (burning prior and after), populations of these species were not adversely affected; similar findings were also observed by Araújo et al. (2004). Moreover, ants are considered ecosystem engineers by increasing soil water infiltration due to their tunnels and can provide valuable ecosystem services (Evans et al., 2010, Lavelle & Spain, 2001) such as building nests on soil and forming small islands with soil organic matter, N and P, which are important elements for plant growth. The importance of the nests is that they present favorable conditions for soil organic matter (SOM) mineralization since they are more exposed to environmental changes. Similar results were obtained by Portilho et al. (2011) in sugarcane monoculture in Brazil obtaining 36.3% ant abundance, thus corroborating with result obtained by Santos et al. (2008) in plant canes, where it was evident its predominance, as well as by Siqueira et al. (2016) with 55%. Even more, it is important to elucidate its use as soil quality bioindicators in highly disturbed areas as they could act as less frequent taxa predators (Andersen & Majer, 2004; Schmidt et al., 2005; Crepaldi et al., 2014).
Soil quality index and Shannon_exp diversity index. The low number of earthworm species (3), total absence of native earthworm species and reduced presence of litter decomposers, such as Diplopoda, Isopoda, Coleoptera, Heteroptera and Orthoptera in the study area explain the low soil quality index observed in our investigation (0-0.10) which is extremely low compared to diverse agroforestry systems (1.00), maize (Zea mays) (0.75), cassava (Manihot sculenta) (0.75), banana (Musa paradisiaca) (0.25), and forests (secondary forest, 30 years) (0.5-0.75) (Huerta et al., 2009; Ramirez et al., 2012). Secondary forest, in our investigation, was not the exception. The history of this land reveals that 15 years in fallow are not enough to recover soil-environmental conditions for macroinvertebrate development after years of sugarcane intense cultivation. The present index can be taken as a good parameter of the ecological health of the soils under sugarcane cultivation, however more studies under local sugarcane conditions are required.
The dominance of a few morphospecies in our study confirms the low species diversity, indicating the higher dominance that particular species from the Hymenoptera, Coleoptera, Chelicerata, Dermaptera, Orthoptera, Chilopoda and Diplopoda Order have on the study community (Magurran, 2004); due to initial selection and persistence of better adapted species to varying management practices (Portilho et al., 2011; Cabrera-Dávila et al., 2017).
Plant macronutrients, yield and plant variables. P and K presented noticeably high values for all sites and management systems, these result could be attributed to the long-term effect of chemical fertilizer programs used the sugarcane region in Belize (Chi et al., 2017), but not for soil biological properties and N levels which were low in the study area. This residual effect of nutrient through long-term application could be beneficial to macrofauna communities since they can offer stable conditions for these organisms, same as those reported by Jiang et al. (2015), which further demonstrated that balanced fertilization in crops (maize-wheat-cotton crop rotation) tended to increase the abundance of soil fauna compared with no fertilizer and unbalanced fertilization such as N, NK, NP and PK. However, this cannot be seen in the present study yet, as the results showed, macroinvertebrate mass, density and diversity were very low.
Overall, yields were high in all experimental plots in comparison to the country’s average yield (44 t ha-1), noting the good development of plant variables in turn indicating the potential of efficient management of chemical and organic amendments for the sugarcane region in Belize. However, yield will depend on application rates, orographic characteristics, type of soil, chemical composition of fertilizers, crop type and climatic conditions (Moran-Salazar et al., 2016). The higher yields in Cambisol (109.7 t ha-1) and in second cycle (first ratoon) (125.50 t ha-1) compare to Luvisol (98.4 t ha-1) and first cycle (82.5 t ha-1) respectively, could be explain by better soil fertility of Cambisol compare to Luvisol and for cycle 2 compare to cycle 1, better establishment of the plants in the field, i.e. stronger development of roots, better climatic conditions, as well as a possible residual effect of organic amendments.
Finally, the proper management of sugarcane cultivation with the aim of balancing soil quality, plant growth and yield, requires a good understanding of all management practices and soils, due to the adverse effect that sugarcane has on productivity. With increasing pressure on sugar demands in Belize, management practices based on organic principles (filter mud and vinasse) will play a crucial role.
CONCLUSIONS
Abundance and diversity of soil macroinvertebrates were not higher in organic treated soils, neither in soils under secondary forest. No significant differences among type of fertilizers and secondary forest on all studied macroinvertebrate parameters were found. Fluvisol had higher values in all evaluated soil macroinvertebrate parameters, except Shannon_exp index. Contrasting, Cambisol had higher N and K values and higher yield. These results indicate that soil type, in a two years trial, had higher weight on determine the studied parameters than the different management systems. Differences between cycles in sugarcane yield could be attributed to seasonal climatic conditions (yearly variation) and better establishment of sugarcane plants for the first ratoon.
Finally, further studies are recommended in sugarcane based on organic amendments application that could maximize yields, improve soil biological properties and reduce costs in Belize.