Next Article in Journal
Key Competencies of Agricultural Managers in the Acute Stage of the COVID-19 Crisis
Next Article in Special Issue
Richness of Rhizosphere Organisms Affects Plant P Nutrition According to P Source and Mobility
Previous Article in Journal
Effects of Different Carbohydrate Sources on Alfalfa Silage Quality at Different Ensiling Days
Previous Article in Special Issue
The Role of Soil N2O Emissions in Agricultural Green Total Factor Productivity: An Empirical Study from China around 2006 when Agricultural Tax Was Abolished
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Earthworm Inoculation Improves Upland Rice Crop Yield and Other Agrosystem Services in Madagascar

1
Laboratoire des Radio-Isotopes, University of Antananarivo, BP 3383, Route d’Andraisoro, Antananarivo 101, Madagascar
2
Eco&Sols, Institut Agro, University Montpellier, CIRAD, INRAe, IRD, 2 Place Viala, 34060 Montpellier, France
*
Author to whom correspondence should be addressed.
Submission received: 11 December 2020 / Revised: 7 January 2021 / Accepted: 11 January 2021 / Published: 13 January 2021
(This article belongs to the Special Issue Thematic of Soil Ecological Functions in Agriculture)

Abstract

:
The effects of earthworm inoculation and cropping systems on upland rice systems were examined over a four-year period in the Highlands of Madagascar. Each year, endogeic earthworms Pontoscolex corethrurus (Rhinodrilidae) were inoculated (EW+) at a density of 75 ind m−2 or were not inoculated (EW0). Inoculation was tested in three cropping systems: conservation agriculture (CA) and traditional tillage with or without residues restitution. Soil and plant properties were measured during the first three years while soil biological properties were assessed at the fourth year. At the end of the experiment, earthworm density was three-fold higher in EW+ than in EW0, demonstrating the success of the inoculation. Earthworm density was more important in CA than in tillage systems. Earthworm inoculation had higher significant effects on soil and plant properties than cropping systems. Earthworm inoculation had positive effects on soil macroaggregation (+43%), aboveground biomass (+27%), rice grain yield (+45%), and N grain amount (+43%). Intensifying earthworm activity in field conditions to meet the challenge of ecological transition is supported by our study.

1. Introduction

Earthworms are involved in many ecological functions leading to the delivery of ecosystem services. Their intense feeding and moving activities have strong impacts (i) on biogeochemical cycles, i.e., carbon (C) transformations [1], nutrient cycling [2], and silicon (Si) availability [3]; (ii) on soil structure and physical properties through the production of biogenic structures [4,5]; (iii) on the control of pests and disease [6,7]; and (iv) on the activity of other soil biota, bacteria [8], fungi [9], nematodes [10], soil microarthropods [11], and (v) on plant growth [12].
Based on this knowledge, the consideration of earthworms for increasing the delivery of ecosystem services from agrosystems (what we call here ‘agrosystem services’) has been proposed for many years [13]. Following their potentially beneficial effects on functions and services, earthworms are generally perceived as animals to be sustained in agricultural systems. It is very likely that the manipulation of the soil community is key to successful restoration of terrestrial ecosystems [14,15]. The main approach consists of changing practices: reducing practices detrimental to earthworms (tillage, pesticides) and developing practices beneficial to earthworms (liming, organic fertilization) [16]. Another approach is inoculation, defined here as the deliberate introduction of living organisms in a given soil with the aim that this action will result in beneficial changes in the dynamics and equilibrium of the environment (i.e., biofertilization). Inoculation, as an active biostimulation approach, has barely been developed at large scales and over a long period [17]. For earthworms, inoculation can concern either cocoons, juveniles, adults, or micro-environments transplantation [18]. Earthworm inoculation needs to be associated to inputs of organic matter aimed at restoring the habitat and favoring earthworm development [19].
Theoretically, it seems possible to inoculate earthworms in all types of soil, under different climates, and for different land-uses, provided that earthworms exist locally and can adapt to soil disturbance. Different trials have been realized in temperate and tropical regions (e.g., [20] in the UK; [21] in Peru; [22] in India; [23] in the West Indies). However, large-scale experiments focusing on agricultural fields are rare and, to our knowledge, only one such assay has been carried out in tea plantations in India [24]. In this experiment, the inoculated earthworms belong to the peregrine species Pontoscolex corethrurus (Rhinodrilidae), and earthworms were inoculated in trenches where tea prunings were buried. This technology was called bio-organic fertilization and resulted in a 240% increase in tea production. An initial objective was to inoculate earthworms at a biomass higher than 30 g m−2 live weight to positively improve/restore soil and plant functions [25]. In the experiment in India, earthworms were reproduced in large covered beds with a cost of production estimated at 3.6 euros kg−1 (i.e., 1400 euros to produce the equivalent of 400 kg live weight per hectare). At this cost, it appears that earthworm inoculation should only be applied to high value crops; this also explains why this technology did not spread worldwide.
In the Highlands of Madagascar, upland crops face many edaphic, environmental, climatic, and economic constraints. Ferralsols are very infertile with many limiting nutrients, especially P, phosphates being rapidly sorbed on clay minerals and oxides, but also N, Ca, Mg [26]. Cereal crop production on these soils is limited by the poor soil fertility, the strong development of pests and diseases (white grubs, blast disease, striga, etc.), and the very low financial ability of smallholder farmers to buy fertilizers and biocides. Fertility improvement is realized with cow manure and, more rarely, pig manure, and crop residues.
Earthworms from Madagascar are relatively well known, with intense surveys being carried out in recent years [27,28,29,30]. In the Highlands, the main species occurring in anthropized environments are exotic species such as P. corethrurus (Rhinodrilidae), Dichogaster bolaui, Dichogaster saliens (Acanthodrilidae), and Amynthas corticis (Megascolecidae). Earthworm densities are higher in perennial than in annual crops, and higher in non-tilled than in tilled situations [31]. P. corethrurus is an endogeic peregrine species found all over the tropics; it is considered as invasive by some authors and tolerates a wide range of environmental conditions [32]. This species has been shown to deeply modify soil structure and soil C and N cycles [32,33], and to increase P availability [34]. This last function is especially important for P-fixing soils such as the Ferralsols from the Highlands of Madagascar. A recent experiment demonstrated that the activity of P. corethrurus can also promote plant health and can especially protect rice against blast disease in Madagascar [7]. A consequence of soil modifications by this species generally results in improved plant growth [21,35].
We hypothesized that soil habitat restoration associated with P. corethrurus inoculation could improve crop production and other agrosystem services, such as C sequestration, and erosion control through an improvement of soil macroaggregation. In the present study, we followed a four-year experimental upland rice (Oryza sativa) crop where earthworms were inoculated or not in three different agricultural systems. Earthworms were inoculated each year at the beginning of the cropping season. Cropping systems were: traditional tillage without residues, traditional tillage with residues, and conservation agriculture, i.e., without tillage and with rotation of an association of maize with the legume dolichos (Dolichos purpureus L.). We hypothesized that earthworm inoculation and consequences on agrosystem services would be more efficient in conservation systems and less efficient in tilled situations without residues.

2. Materials and Methods

2.1. Study Site

The experiment was conducted at Lazaina (18°46′55.59″ S, 47°32′46.3″ E, 1274 m above sea level), in the Analamanga region of the Malagasy highlands. The climate is a tropical altitude climate with an average annual temperature of 20 °C and mean annual rainfall of 1300 mm. The area is characterized with a warm and humid season from October to April and a cool and dry season from May to September. The upland rice cropping season starts at the end of November and ends in April. The study was conducted over four growing seasons: 2013–2014, 2014–2015, 2015–2016, and 2016–2017. Rainfall was higher in the period 2014–2015 (1232 mm during the cropping season, December–April) than in other periods (888–1082 mm).
Soils are Ferralsols with high contents of kaolinitic clays and Al- and Fe-oxides. Under natural vegetation, mainly composed of grasses Aristida sp., soil has a sandy-clay texture with 33% clay, a pH of 5.5, a total organic C content of 20.8 g kg−1, a total N content of 1.3 g kg−1, a C:N of 16, a total P content of 380 mg kg−1, an Olsen P content of 7.12 mg kg−1, and exchangeable K, Ca, and Mg contents of 30.7, 120.7, and 28.3 mg kg−1, respectively. The cationic exchange capacity (CEC) was 1.34 cmol+ kg−1 and the water-holding capacity (WHC) was 0.49 g g−1 dry soil.

2.2. Experimental Design and Plot Preparation

The experimental plot was under Aristida savanna for more than 30 years before being cropped with maize between 2006 and 2012; this was followed by a one-year fallow period before our experiment started in 2013. Four plots of the previous experiment were thus retained to represent the four blocks of our experiment. On the first year of experiment, all subplots were tilled, with an ‘angady’ (traditional iron spade) at 20–25 cm depth and then carefully homogenized within the block, to homogenize the soil. Two factors were tested in this experiment: (i) the inoculation of earthworms with two modalities: with earthworm inoculation (EW+) or without earthworm inoculation (EW0), and (ii) the cropping system with three modalities: traditional tillage with rice residues restitution each year (TR+), traditional tillage without rice residues restitution (TR0), and a conservation agriculture system (CA) including two rotations: rice in rotation with an association of maize (Zea mays) + dolichos (Dolichos purpureus) and association of maize + dolichos in rotation with rice. The soil was not tilled in the CA systems after the first cropping year. The experiment is thus composed of 32 subplots (2 earthworm treatments, 4 cropping systems, 4 replicates, Figure 1). Each subplot had a 5.6 m2 surface area. In the experimental field, the land presented a slope of 8%. To limit soil disturbance and soil runoff flows, all treatments with earthworms (EW+) in a block were gathered in the upper part of the block, while treatments without earthworms were gathered in the lower part. We are aware that this choice will have consequences on the statistics used; thus, this ‘plot position’ was considered in further statistical analysis. The area gathering the four treatments of EW+ were surrounded by a metal sheet buried vertically in the soil up to 40 cm deep and protruding 5 cm above the ground surface. The aim of this was to prevent earthworms from escaping EW+ treatments.
The sowing took place in mid-December. Rice holes were 20 cm apart with 7–8 seeds in each hole, equivalent to 80 kg seed ha−1. We used the rice cultivar ‘FOFIFA 161’. For maize + dolichos, rows of maize were alternated 45 cm with rows of dolichos. Maize seeds (‘IRAT 200’, 3 per hole) and dolichos seeds (D. purpureus with brown seeds, 3 per hole) were 45 cm apart on a row and 40 cm apart on a row, respectively. Rice and maize + dolichos were both fertilized with manure and NPK fertilizer (11% N, 22% P, 16% K) at doses equivalent to 10 Mg ha−1 and 200 kg ha−1, respectively. At the end of the vegetative stage (stem elongation), urea was supplied at a rate of 50 kg ha−1. Weeding was carried out after one month of sowing and was repeated every 15 days until the flowering stages. The main insect pests (Apoderus humeralis, Coleoptera Curculionidae) were monitored and hand-removed.

2.3. Earthworm Inoculation

For this experiment, we used the species P. corethrurus, largely distributed in the anthropized areas of the Highlands. They were collected from an area close to our experiment (less than 100 m), under eucalyptus trees growing in the ‘Bozaka’ savanna. Earthworms were kept in a basin with wet soil. Each year, we introduced earthworms in subplots EW+ to reach a density of 75 ind m−2, slightly higher than the threshold (60 ind m−2, i.e., 30 g m−2) proposed by [25] to obtain a significant effect of earthworms on plant growth. This means a quantity of 6680 earthworms for the whole experiment, i.e., 350 per subplot, each year. Three days, with five people, were necessary to reach this number each year. After collecting the worms, they were transferred to buckets with water, in order to count them and to remove dead or injured worms and earthworms from other species.

2.4. Soil and Plant Analyses

Soil was carefully characterized at the beginning of the experiment in order to get an initial state of soil properties. Soil (0–10 cm) was removed using an auger, with 6 samples per subplot gathered to make a composite sample. Soil C (Walkley–Black method) and N (Kjeldhal method) contents, and exchangeable (resin) inorganic P content (Pi) were measured in the laboratory. Bulk density (0–10 cm) was also measured with 2 cylinders per subplot.
At the end of each cropping season, four samples per subplot were gathered to make a composite sample for analyses. Bulk density was also measured with 2 cylinders per subplot. Soil macroaggregation was also measured by dry sieving at 2 mm [36]. We calculated the percentage of macroaggregates as the weight of aggregates retained by the sieve, divided by the whole weight of sieved soil.
Rice height was regularly measured, each year, at each important development stage, i.e., beginning (one month after sowing) and end (two months after sowing) of the tillering stages, end of flowering stage (3 months after sowing). At rice harvest (end of April), aboveground biomass (sum of shoot and grain biomasses) and grain yield were measured for all plants except dead plants and the border line. Shoots and grains were kept and air-dried for further N and P analyses. Root biomass was measured by hand-collecting all visible roots from two soil monoliths (25 × 25 × 10 cm) per subplot. Roots were then washed, air-dried and weighted.
One month before the end of the 4th cropping season, before the end of the rains, earthworms and nematodes were sampled. Earthworms were hand-sorted from two soil monoliths (25 × 25 × 10 cm) per subplot; they were counted and identified at a species level in the field and put back in the soil. Regarding soil nematodes, two soil samples (0–10 cm) were taken per subplot. Nematodes were extracted in the laboratory, from 100 g fresh soil, using a modified Baermann funnel method [37]. After 48 h of active filtering, nematodes in water were collected and counted using a stereomicroscope.
Here we focus on plant results obtained after the third cropping season in order to best perceive the combined effects of earthworm inoculation and cropping systems. This is also the year when crops were not disturbed by cattle trampling or climatic events.

2.5. Statistical Analyses

Data analyses were performed using the R statistical software [38] with a p-value threshold set at 10%. Prior to analyzing the effects of earthworm inoculation and cropping system on soil and plant properties, we first tested the soil homogeneity at the beginning of the experiment. Thus, initial soil properties were submitted to a one-way analysis of variance (ANOVA) (plot position as factor). A highly significant (p-value < 0.001) (Table 1) effect of plot position was observed for soil chemical variables (total C, total N and Pi). Due to the soil heterogeneity, we tested the main effects (earthworm inoculation and cropping system) on soil properties at the end of cropping seasons by using analysis of covariance (ANCOVA) in which covariates are initial soil properties. To determine the direction of significant effects of ANCOVA, we used multiple comparison tests based on estimated marginal means (hereafter emmeans, package emmeans) with Tukey adjustment.
Regarding the changes in soil properties with time, the difference between mean final soil properties and mean initial soil properties was calculated. The obtained differences can be positive (meaning a gain) or negative (meaning a loss). For instance, difference in soil carbon content was calculated as: ∆C (in g kg−1) = final C content − initial C content. This was calculated in the same way for P content and C and P stocks.
For agronomic parameters and soil biological properties (earthworm and nematode densities), the absence of significant effect of covariates on physico-chemical soil properties observed with ANCOVA (Table 2) led us to realize a three-way ANOVA (earthworm inoculation, cropping system, block) on the whole data. This analysis was followed by a Tukey HSD (honestly significant difference) post hoc test for multiple mean comparison. The normality of the data was verified for all ANOVA models on residuals using the Shapiro-Wilk test.

3. Results

3.1. Soil Biological Properties

At the rainy period of the 4th cropping season, ANOVA results showed that earthworm inoculation significantly increased the total density of P. corethrurus (sum of adults and juveniles) (p-value = 0.002). On average, it was three-fold higher in EW+ treatments than in EW0 treatments (91 ± 11 (individuals) ind m2 vs. 32 ± 12 ind m2, respectively). We noticed that the proportion of the density of juveniles to the density of total earthworms was higher in EW+ (25%) than in EW0 (7%) (Figure 2). Regarding nematode density, we found no significant effect of earthworm inoculation (p-value = 0.106) despite higher values in EW+ (4.07 ± 0.71 ind g1 soil) than in EW0 (2.84 ± 0.45 ind g1 soil).
For the cropping system factor, we observed that the density of earthworms tended to be significantly increased (p-value = 0.064) in CA treatments (89 ± 23 ind m2) compared to tillage systems (49 ± 12 ind m2 for TR0 and 47 ± 13 ind m2 for TR+) (Figure 2). However, ANOVA results showed no significant effect on nematode density, with p-value of 0.611. Additionally, the interaction of earthworm inoculation and cropping system had no influence on earthworm (p-value = 0.740) and nematode (p-value = 0.391) densities.

3.2. Soil Physico-Chemical Properties

3.2.1. Soil Properties at the Beginning of the Experiment

ANOVA showed that plot position had a significant effect on initial soil properties, despite soil homogenization before cultivation (Table 1). Plots without earthworm inoculation (lower parts of the four blocks) had significant higher total C (+17.5%), total N (+17.9%), and Pi (+146%) contents, compared to plots where earthworms were inoculated. ANOVA also showed an effect of blocks, block 1 showing significantly lower N contents and stocks compared to other blocks (data not shown).

3.2.2. Soil Properties at the End of the 3rd Cropping Season

At the end of the 3rd cropping season, we found that soil total C, NH4+, and Pi contents in EW0 were higher compared to values measured in EW+, with an average increase of 3%, 7%, and 29%, respectively (Table 2). However, according to ANCOVA results, the difference between earthworm treatments was not significant (Table 3). Additionally, total C, NH4+, and Pi stocks were also not affected by cropping system and initial soil characteristics. A significant block effect was found on total C content and stock. Regarding soil physical properties, the earthworm inoculation had no significant effect on bulk density (p-value = 0.214). Conversely, ANOVA showed that soil macroaggregates were significantly affected by earthworm inoculation (p-value < 0.001) (Table 3). In EW+, the proportion of soil macroaggregates to the total soil was 28%, while it was 17% in EW0 (Table 2). The cropping system had no significant effect on soil macroaggregates and bulk density, with p-values of 0.384 and 0.702, respectively (Table 3).

3.2.3. Changes in Soil Properties with Time

Earthworm inoculation had no significant effect on the difference between C content at the 3rd cropping season and the initial C content (∆C content) (p-value = 0.135, Table 4 and Table 5). The same pattern was observed for ∆C stock with a p-value of 0.218. C stock changes were significantly modified by the cropping system, with a higher C stock increase for conservation agriculture than for tillage systems: ∆C stocks were 13.33 ± 1.42 Mg ha−1 in CA, 9.06 ± 1.29 Mg ha−1 in TR+ and 9.63 ± 0.79 Mg ha−1 in TR0 (p-value = 0.056). With regards to Pi stock changes, neither earthworm inoculation nor cropping system had a significant effect, with p-values of 0.659 and 0.120, respectively.

3.3. Plant Parameters

3.3.1. Rice Growth at the End of the 3rd Cropping Season

The effect of earthworm inoculation on rice height depends on the physiological stages: At the beginning of the tillering stage, there was no significant effect (p-value = 0.306), while at the end of the tillering stage and at the flowering stage, rice height was significantly higher in the presence of earthworms (p-value = 0.028 and p-value = 0.008, respectively). At the end of the flowering stage, rice height was 10 cm higher with earthworms than without (61.1 ± 1.6 cm and 51.3 ± 2.5 cm, respectively) (Figure 3). Similarly, rice shoot biomass at the harvest significantly increased by 26% in treatments where earthworms were inoculated (p-value = 0.060). The shoot:root ratio was also significantly affected by earthworm inoculation (p-value = 0.091, 2.62 ± 0.64 in EW+ and 1.65 ± 0.59 in EW0 treatments) but not root biomass (p-value = 0.984). We observed a significant effect of the cropping system on rice height at the beginning of the tillering stage (p-value = 0.028) and non-significant effects at later stages (p-value = 0.351 at the end of the tillering stage, p-value = 0.403 at the flowering stage). At the beginning of tillering, rice height was lower in CA (16.3 ± 0.9 cm) than in both TR treatments (18.9 ± 0.4 cm for TR0 and 18.0 ± 0.5 cm for TR+). At the harvest, no significant effect of cropping system was found on shoot rice biomass (p-value = 0.300), but it was significant for root biomass (p-value = 0.094) with lower values in CA compared to tillage systems (Figure 3). Likewise, the shoot:root ratio was not significantly influenced by cropping system (p-value = 0.352).

3.3.2. Rice Grain Yield and Nutrient Acquisition at the end of the 3rd Cropping Season

Statistical analyses showed that earthworm inoculation significantly affected the rice grain yield at the end of the third cropping season (p-value of 0.040, Table 6). We found that earthworm presence had a positive effect on rice yield (+44.6%). In EW+ treatments, rice grain yield was 1.20 ± 0.13 Mg ha−1 while it was 0.83 ± 0.09 Mg ha−1 in EW0. Similarly, rice nutrition was enhanced in treatments where earthworms were inoculated. Grain N and P amounts were respectively, significantly, increased by 43% and 28% in EW+ treatments compared to EW0 treatments (Table 6).
Regarding cropping system, rice grain yield tended to be lower in CA (0.88 ± 0.15 Mg ha−1) than in tillage systems (1.18 ± 0.12 Mg ha−1 for TR0 and 1.04 ± 0.17 Mg ha−1 for TR+), but the difference was not significant (p-value = 0.281). Additionally, the interaction of earthworm inoculation and cropping system had no significant effect on any of the plant variables (Table 6).

4. Discussion

4.1. Success of Earthworm Inoculation

In our experiment, we did not eliminate native earthworms from the experimental field before earthworm inoculation. The earthworm density measured at the end of the experiment in the treatment without earthworm inoculation thus characterizes the earthworm population developing naturally in crop conditions. This density of 32 ind m−2 is similar to other observations from upland agricultural fields in the Highlands of Madagascar, which also showed the dominance of the species P. corethrurus [31,39]. This value is half the density recommended by [25,40] (30 g m−2, i.e., 60 ind m−2) to observe a significant positive effect of earthworms on plant growth. With earthworm inoculation, we reached a density of 91 ind m−2, which is above this recommended value. There is no evidence regarding if this value results only from the last inoculation (Year 4), or from the inoculations of the four years of the experiment. The relatively high presence of juveniles in treatments EW+ may indicate that earthworms survive and develop well in this type of environment, especially when they are introduced at relatively high densities. These observations indicate evidence of earthworm activities throughout the experiment and prove the success of earthworm inoculation in upland agricultural field conditions in the Highlands of Madagascar. Survival and reproduction of introduced earthworms reinforce the interest of this inoculation technique. Conditions were thus reached to observe a positive effect of earthworms on the different plant and soil properties measured. Complementary studies should be realized to analyze whether inoculation should be practiced each year or whether initial inoculation during the first years result in the installation of a stable and abundant earthworm population. It was also interested to notice that earthworm density doubled in CA compared to tillage systems. This confirms the interest of such cropping systems without soil tillage, alternative to traditional ones, in increasing soil biota and especially earthworms [39,41,42,43]. The positive effect of the CA system is imputed to the reduction of soil tillage and the presence of plant cover, bringing organic matter and buffering climatic stresses.

4.2. Effects on Soil and Plant Properties

In our experiment, soil properties at the beginning of the experiment were variable, especially between subplots where earthworms were introduced and subplots without earthworm inoculation. In that sense, it seems more relevant and interesting to focus on the changes in soil properties with time in the different plots.
We measured that C stock changes after 3 years (∆C stocks) were slightly more important in EW+ than in EW0, but this was not significant. This may indicate that our experiment was too short to observe a positive effect of earthworms on C storage as it has been observed in other experiments [44].
The inoculation of earthworms resulted in significant effects on plant properties at the third year of the experiment. Earthworms increased plant growth (height at the end of the tillering stage and at the flowering stage). The absence of effect on rice height at the beginning of the tillering stage may be explained by the short span (only one week) between earthworm inoculation and this measure of height. This highlights the importance of inoculating earthworms rapidly after sowing to observe rapid and positive effects on plant growth. The positive effect of earthworms was measured for shoot plant biomass (+26% in EW+ compared to EW0) and rice grain yield (+44.6%) as well. This fits with conclusions given by [12] in a meta-analysis, who gives an average increase in aboveground biomass of 23% and in yield of 25%. This increase is particularly visible for rice, which aboveground biomass strongly increased in the presence of earthworms compared to treatments without earthworms: +35% in Peru [21] and +42% in the review article by [40]. This positive effect of earthworms on plant growth and yield can be explained by different mechanisms [45]: increase in mineralization of organic matter [8] and higher availability of nutrients, especially N and P [2,34]; modification of soil porosity and aggregation improving water retention [46]; stimulation of microbial activity and production of auxin-like substances [47]; regulation of pests and pathogens [7,48]; stimulation of microbial plant symbiotes [45]; horizontal transfer of fertility through earthworm biomass decomposition [34]. In our experiment, we observed higher grain N and P amounts in the presence of earthworms, indicating higher availability of nutrients with earthworms. This is also confirmed by higher values of shoot:root ratio in the presence of earthworms, suggesting that plants are less constrained by nutrient acquisition and can invest less energy in roots and more energy in the aerial parts, following the theory of ‘functional equilibrium’ [49]. This biomass allocation towards aerial parts was observed in different studies (for instance [50,51]).
In contrast, cropping systems had no significant effects, neither on soil properties nor on plant growth, nutrition, or yield. We only observed a higher increase in C storage (∆C storage) in CA compared to tillage treatments, which was already described for Ferralsols of the Highlands of Madagascar [52]. These authors also observed an increase in P stock under CA aged 8 years old, which is longer than our experiment. This also confirms that it takes some years before the positive effects of CA on C and nutrient stocks become significant [53,54]. It was also noticed that the amount of crop residues in CA is important for defining the amplitude of CA effects on soil and plant properties. For instance, [55] noticed a significant positive effect of CA on weed pressure for a residue amount of 10 Mg ha−1. In our experiment, the residue amount was lower than this value (<5 mg ha−1).

5. Conclusions

In the Highlands of Madagascar, the inoculation of earthworms appears to be an efficient technology to improve both crop productivity (yield) and sustainability; sustainability seen here as the improvement of ecological functions and services such as increase in soil biota, C storage, and soil aggregation. The average surface of farms in this region is 0.9 ha, only 22% of the surface being used for upland cultivation, indicating small surfaces for upland fields. Earthworms of the species P. corethrurus, well known by farmers and locally called ‘kanka-mena’ (literally red worms) can be found easily, at high densities, in different places in the environment: around lowland fields, under trees, or under perennial grasses. It does not take much time to collect thousands of earthworms and to inoculate them in fields, at the beginning of the rainy season. This seems feasible with regards to the very small size of rain-fed rice plots in the family farms in the Highlands of Madagascar.
Inoculated earthworms survived at least during the whole rainy season and especially in CA situations characterized by no-tillage and rotation of rice with maize+dolichos. Inoculation resulted in positive effects on different functions and services, especially soil aggregation, rice nutrition, and rice yield (+44.6%). Studies giving such results in field conditions are relatively rare. Our study brings evidence that intensifying earthworm activity in cropped soils may increase upland rice yield and other agrosystem services such as C storage and erosion control.

Author Contributions

Conceptualization, O.R. and E.B.; methodology, O.R., E.B., M.R. (Malalatiana Razafindrakoto), T.R. and M.R. (Michel Rabenarivo); validation, T.R., L.B., T.B. and J.T.; formal analysis, O.R., E.B. and M.R. (Malalatiana Razafindrakoto); investigation, O.R., E.B. and M.R. (Malalatiana Razafindrakoto); data curation, O.R., E.B. and M.R. (Malalatiana Razafindrakoto); writing—original draft preparation, O.R. and E.B.; writing—review and editing, all authors; supervision, E.B.; project administration, E.B.; funding acquisition, E.B. All authors have read and agreed to the published version of the manuscript.

Funding

This study was partly funded by the project SECuRE (Soil ecological function restoration to enhance agrosystem services in rainfed rice cropping systems in agroecological transition), supported by the Agropolis Foundation under the reference ID 1605-007 through the “Investissements d’Avenir” program (Labex Agro: ANR-10-LABX-0001-01) under the framework I-SITE MUSE (ANR-16-IDEX-0006). The IRD (French Institute of Research for the Development) also contributed to the funding via the ERA-Net ARD II call (2012), CAMES Project (Impacts of Conservation Agriculture on Macrofauna diversity and related Ecosystem Services for improved farmer’s cropping systems and livelihoods in Highlands of Madagascar).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

We are thankful to Modeste Rakotondramanana, Marie-Paule Razafimanantsoa, and Koloina Rahajaharilaza for technical help.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Zhang, W.; Hendrix, P.F.; Dame, L.E.; Burke, R.A.; Wu, J.; Neher, D.A.; Li, J.; Shao, Y.; Fu, S. Earthworms Facilitate Carbon Sequestration through Unequal Amplification of Carbon Stabilization Compared with Mineralization. Nat. Commun. 2013, 4, 2576. [Google Scholar] [CrossRef] [PubMed]
  2. Van Groenigen, J.W.; Van Groenigen, K.J.; Koopmans, G.F.; Stokkermans, L.; Vos, H.M.; Lubbers, I.M. How Fertile Are Earthworm Casts? a Meta-Analysis. Geoderma 2019, 338, 525–535. [Google Scholar] [CrossRef]
  3. Bityutskii, N.P.; Kaidun, P.; Yakkonen, K. Earthworms Can Increase Mobility and Bioavailability of Silicon in Soil. Soil Biol. Biochem. 2016, 99, 47–53. [Google Scholar] [CrossRef]
  4. Blanchart, E.; Albrecht, A.; Brown, G.G.; Decaëns, T.; Duboisset, A.; Lavelle, P.; Mariani, L.; Roose, E. Effects of Tropical en-Dogeic Earthworms on Soil Erosion: A Review. Agric. Ecosyst. Environ. 2004, 104, 303–315. [Google Scholar] [CrossRef]
  5. Blanchart, E.; Marilleau, N.; Chotte, J.-L.; Drogoul, A.; Perrier, E.; Cambier, C. SWORM: An Agent-Based Model to Simulate the Effect of Earthworms on Soil Structure. Eur. J. Soil Sci. 2008, 60, 13–21. [Google Scholar] [CrossRef]
  6. Bertrand, M.; Blouin, M.; Barot, S.; Charlier, A.; Marchand, D.; Roger-Estrade, J. Biocontrol of Eyespot Disease on Two Winter Wheat Cultivars by an Anecic Earthworm (Lumbricus terrestris). Appl. Soil Ecol. 2015, 96, 33–41. [Google Scholar] [CrossRef] [Green Version]
  7. Blanchart, E.; Ratsiatosika, O.; Raveloson, H.; Razafimbelo, T.; Razafindrakoto, M.; Sester, M.; Becquer, T.; Bernard, L.; Trap, J. Nitrogen Supply Reduces the Earthworm-Silicon Control on Rice Blast Disease in a Ferralsol. Appl. Soil Ecol. 2020, 145, 103341. [Google Scholar] [CrossRef]
  8. Bernard, L.; Chapuis-Lardy, L.; Razafimbelo, T.; Razafindrakoto, M.; Pablo, A.-L.; Legname, E.; Poulain, J.; Brüls, T.; O’Donohue, M.; Brauman, A.; et al. Endogeic Earthworms Shape Bacterial Functional Communities and Affect Organic Matter Mineralization in a Tropical Soil. ISME J. 2011, 6, 213–222. [Google Scholar] [CrossRef] [Green Version]
  9. Butenschoen, O.; Poll, C.; Langel, R.; Kandeler, E.; Marhan, S.; Scheu, S. Endogeic Earthworms Alter Carbon Translocation by Fungi at the Soil–Litter Interface. Soil Biol. Biochem. 2007, 39, 2854–2864. [Google Scholar] [CrossRef]
  10. Villenave, C.; Rabary, B.; Kichenin, E.; Djigal, D.; Blanchart, E. Earthworms and Plant Residues Modify Nematodes in Tropical Cropping Soils (Madagascar): A Mesocosm Experiment. Appl. Environ. Soil Sci. 2010, 2010, 1–7. [Google Scholar] [CrossRef] [Green Version]
  11. Eisenhauer, N. The Action of an Animal Ecosystem Engineer: Identification of the Main Mechanisms of Earthworm Impacts on Soil Microar-Thropods. Pedobiologia 2010, 53, 343–352. [Google Scholar] [CrossRef]
  12. Van Groenigen, J.W.; Lubbers, I.M.; Vos, H.M.J.; Brown, G.G.; De Deyn, G.B.; Van Groenigen, K.J. Earthworms Increase Plant Production: A Meta-Analysis. Sci. Rep. 2015, 4, 6365. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Lavelle, P.; Barois, I.; Martin, A.; Zaidi, Z.; Schaefer, R. Management of Earthworm Populations in Agro-Ecosystems: A Possible Way to Maintain Soil Quality? In Ecology of Arable Land—Perspectives and Challenges; Springer Science and Business Media LLC: Berlin/Heidelberg, Germany, 1989; pp. 109–122. [Google Scholar]
  14. Jouquet, P.; Blanchart, E.; Capowiez, Y. Utilization of Earthworms and Termites for the Restoration of Ecosystem Functioning. Appl. Soil Ecol. 2014, 73, 34–40. [Google Scholar] [CrossRef]
  15. Wubs, E.R.J.; Van Der Putten, W.H.; Bosch, M.; Bezemer, T.M. Soil Inoculation Steers Restoration of Terrestrial Ecosystems. Nat. Plants 2016, 2, 16107. [Google Scholar] [CrossRef] [Green Version]
  16. Brun, J.J.; Cluzeau, D.; Tréhen, P.; Bouché, M.B. Biostimulation: Perspectives et Limites de L’Amélioration Biologique Des Sols par Stimulation ou Introduction d’Espèces Lombriciennes. Rev. Ecol. Biol. Sols 1987, 24, 685–701. [Google Scholar]
  17. Butt, K.R. Earthworms in Soil Restoration: Lessons Learned from United Kingdom Case Studies of Land Reclamation. Restor. Ecol. 2008, 16, 637–641. [Google Scholar] [CrossRef]
  18. Butt, K.R.; Frederickson, J.; Morris, R.M. The Earthworm Inoculation Unit Technique: An Integrated System for Cultivation and Soil-Inoculation of Earthworms. Soil Biol. Biochem. 1997, 29, 251–257. [Google Scholar] [CrossRef]
  19. Lavelle, P.; Barros, E.; Blanchart, E.; Brown, G.G.; Desjardins, T.; Mariani, L.; Rossi, J.-P. SOM Management in the Tropics: Why Feeding the Soil Macrofauna? Nutr. Cycl. Agroecosyst. 2001, 61, 53–61. [Google Scholar] [CrossRef]
  20. Butt, K.R. Inoculation of Earthworms into Reclaimed Soils: The UK Experience. Land Degrad. Dev. 1999, 10, 565–575. [Google Scholar] [CrossRef]
  21. Pashanasi, B.; Lavelle, P.; Alegre, J.; Charpentier, F. Effect of the Endogeic Earthworm Pontoscolex Corethrurus on Soil Chemical Characteristics and Plant Growth in a Low-Input Tropical Agroecosystem. Soil Biol. Biochem. 1996, 28, 801–810. [Google Scholar]
  22. Senapati, B.; Lavelle, P.; Giri, S.; Pashanasi, B.; Alegre, J.; Decaëns, T.; Jimenez, J.J.; Albrecht, A.; Blanchart, E.; Mahieu, M.; et al. In-Soil Earthworm Technologies for Tropical Agroecosystems. In Earthworm Management in Tropical Agroecosystems; Lavelle, P., Brussaard, L., Hendrix, P., Eds.; CABI Publishing: Wallingford, UK, 1999; pp. 199–238. [Google Scholar]
  23. Blanchart, E.; Albrecht, A.; Chevallier, T.; Hartmann, C. The Respective Roles of Biota (Roots and Earthworms) in the Restoration of Physical Properties in Vertisol under a Digitaria De-cumbens Pasture (Martinique). Agric. Ecosyst. Environ. 2004, 103, 343–355. [Google Scholar] [CrossRef]
  24. Senapati, B.K.; Lavelle, P.; Panigrahi, P.K.; Giri, S.; Brown, G.G. Restoring Soil Fertility and Enhancing Productivity in Indian Tea Plantations with Earthworms and Organic Fertilizers. In Proceedings of the International Technical Workshop on Biological Management of Soil Ecosystems for Sustainable Agriculture, Londrina, Brazil, 24–27 June 2002; FAO (Food and Agriculture Organization): Rome, Italy, 2003; Volume 1, pp. 172–190. [Google Scholar]
  25. Spain, A.; Lavelle, P.; Mariotti, A. Stimulation of Plant Growth by Tropical Earthworms. Soil Biol. Biochem. 1992, 24, 1629–1633. [Google Scholar] [CrossRef]
  26. Raminoarison, M.; Razafimbelo, T.; Rakotoson, T.; Becquer, T.; Blanchart, E.; Trap, J. Multiple-Nutrient Limitation of Upland Rainfed Rice in Ferralsols: A Greenhouse Nutrient-Omission Trial. J. Plant Nutr. 2019, 43, 270–284. [Google Scholar] [CrossRef]
  27. Razafindrakoto, M.; Csuzdi, C.S.; Blanchart, E. New and Little-Known Giant Earthworms from Madagascar (Oligochaeta: Ky-notidae). Afr. Invertebr. 2011, 52, 285–294. [Google Scholar] [CrossRef]
  28. Csuzdi, C.; Razafindrakoto, M.; Blanchart, E. New and Little-Known Earthworm Species from Central Madagascar (Oligo-chaeta: Kynotidae). Zootaxa 2012, 3578, 36–42. [Google Scholar] [CrossRef] [Green Version]
  29. Razafindrakoto, M.; Csuzdi, C.; James, S.; Blanchart, E. New Earthworms from Madagascar with Key to the Kynotus Species (Oligochaeta: Kynotidae). Zool. Anz. A J. Comp. Zool. 2017, 268, 126–135. [Google Scholar] [CrossRef]
  30. Hong, Y.; Razafindrakoto, M.; Blanchart, E.; Szederjesi, T.; Csuzdi, C. New Acanthodrilid Species from Madagascar (Clitellata, Acanthodrilidae). Acta Zool. Acad. Sci. Hung. 2019, 65, 215–233. [Google Scholar] [CrossRef]
  31. Rakotomanga, D.; Blanchart, E.; Rabary, B.; Randriamanantsoa, R.; Razafindrakoto, M.; Autfray, P. Diversité de la Macrofaune Des Sols Cultivés Sur Les Hautes Terres de Madagascar. Biotechnol. Agron. Soc. Environ. 2016, 20, 495–507. [Google Scholar]
  32. Taheri, S.; Pelosi, C.; Dupont, L. Harmful or Useful? A Case Study of the Exotic Peregrine Earthworm Morphospecies Pontoscolex Corethrurus. Soil Biol. Biochem. 2018, 116, 277–289. [Google Scholar] [CrossRef]
  33. Coq, S.; Barthès, B.G.; Oliver, R.; Rabary, B.; Blanchart, E. Earthworm Activity Affects Soil Aggregation and Organic Matter Dynamics According to the Quality and Localization of Crop Residues—an Experimental Study (Madagascar). Soil Biol. Biochem. 2007, 39, 2119–2128. [Google Scholar] [CrossRef]
  34. Trap, J.; Blanchart, E.; Ratsiatosika, O.; Razafindrakoto, M.; Becquer, T.; Andriamananjara, A.; Morel, C. Effects of the Earthworm Pontoscolex Corethrurus on Rice P Nutrition and Plant-Available Soil P in a Tropical Ferralsol. Appl. Soil Ecol. 2021, 160, 103867. [Google Scholar] [CrossRef]
  35. Fonte, S.J.; Quintero, D.C.; Velásquez, E.; Lavelle, P. Interactive Effects of Plants and Earthworms on the Physical Stabilization of Soil Organic Matter in Aggregates. Plant Soil 2012, 359, 205–214. [Google Scholar] [CrossRef]
  36. André, J.; Villenave, C.; Charpentier, F.; Bernier, N.; Lavelle, P.; Ponge, J.; Derouard, L.; Pashanasi, B.; Rossi, J.-P. Effects of Earthworms on Soil Organic Matter and Nutrient Dynamics at a Landscape Scale over Decades. In Earthworm Ecology; Informa UK Limited: London, UK, 2004; pp. 145–160. [Google Scholar]
  37. Hooper, D.J. Extraction of Free-Living Stages from Soil. In Laboratory Methods for Work with Plant and Soil Nematodes; Southey, J.F.E., Ed.; Her Majesty’s Stationary Office: London, UK, 1986; pp. 5–30. [Google Scholar]
  38. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2014. [Google Scholar]
  39. Rabary, B. Impact du Semis Direct Sous Couverture Végétale Sur la Macrofaune ET la Microflore des Sols Ferrallitiques D’Andranomanelatra, Antsirabe, Hautes Terres Malgaches. Ph.D. Thesis, University of Antananarivo, Antananarivo, Madagascar, 2011; p. 211. [Google Scholar]
  40. Brown, G.G.; Pashanasi, B.; Villenave, C.; Patron, J.C.; Senapati, B.K.; Giri, S.; Barois, I.; Lavelle, P.; Blanchart, E.; Blakemore, R.J.; et al. Effects of Earthworms on Plant Production in the Tropics. in Earthworm Management. In Tropical Agroecosystems; Lavelle, P., Brussaard, L., Hendrix, P., Eds.; CABI Publishing: Wallingford, UK, 1999; pp. 87–148. [Google Scholar]
  41. Blanchart, E.; Bernoux, M.; Sarda, X.; Siqueira Neto, M.; Cerri, C.C.; Piccolo, M.; Douzet, J.M.; Scopel, E.; Feller, C. Effect of Direct Seeding Mulch-Based Systems on Soil Carbon Storage and Macrofauna in Central Brazil. Agric. Conspec. Sci. 2007, 72, 81–87. [Google Scholar]
  42. Eriksen-Hamel, N.S.; Speratti, A.B.; Wang, W.; Légère, A.; Madramootoo, C.A. Earthworm Populations and Growth Rates Related to Long-Term Crop Residue and Tillage Management. Soil Tillage Res. 2009, 104, 311–316. [Google Scholar] [CrossRef]
  43. Henneron, L.; Bernard, L.; Hedde, M.; Pelosi, C.; Villenave, C.; Chenu, C.; Bertrand, M.; Girardin, C.; Blanchart, E. Fourteen Years of Evidence for Positive Effects of Conservation Agriculture and Organic Farming on Soil Life. Agron. Sustain. Dev. 2014, 35, 169–181. [Google Scholar] [CrossRef] [Green Version]
  44. Lavelle, P.; Spain, A.V.; Fonte, S.; Bedano, J.C.; Blanchart, E.; Galindo, V.; Grimaldi, M.; Jimenez, J.J.; Velasquez, E.; Zangerlé, A. Soil Aggregation, Ecosystem Engineers and the C Cycle. Acta Oecol. 2020, in press. [Google Scholar] [CrossRef]
  45. Brown, G.G.; Edwards, C.A.; Brussaard, L. How Earthworms Affect Plant Growth: Burrowing into the Mechanisms. In Earth-Worm Ecology, 2nd ed.; Edwards, C.A., Ed.; CRC Press: Boca Raton, FL, USA, 2004; Chapter 2; p. 441. [Google Scholar]
  46. Blanchart, E.; Albrecht, A.; Alegre, J.; Duboisset, A.; Villenave, C.; Pashanasi, B.; Lavelle, P.; Brussaard, L. Effects of Earthworms on Soil Structure and Physical Properties. In Earthworm Management in Tropical Agroecosystems; Lavelle, P., Brussaard, L., Eds.; CABI Publishing: Wallingford, UK, 1999; pp. 149–172. [Google Scholar]
  47. Muscolo, A.; Bovalo, F.; Gionfriddo, F.; Nardi, S. Earthworms Humic Matter Produces Auxin-Like Effects on Daucus Carota Cell Growth and Nitrate Metabolism. Soil Biol. Biochem. 1999, 31, 1303–1311. [Google Scholar] [CrossRef]
  48. Blouin, M. Chemical Communication: An Evidence for Co-evolution between Plants and Soil Organisms. Appl. Soil Ecol. 2018, 123, 409–415. [Google Scholar] [CrossRef]
  49. Brouwer, R. Functional Equilibrium: Sense or Nonsense? Neth. J. Agric. Sci. 1983, 31, 335–348. [Google Scholar] [CrossRef]
  50. Kreuzer, K.; Bonkowski, M.; Langel, R.; Scheu, S. Decomposer Animals (Lumbricidae, Collembola) and Organic Matter Distri-Bution Affect the Performance of Lolium Perenne (Poaceae) and Trifolium Repens (Fabaceae). Soil Biol. Biochem. 2004, 36, 2005–2011. [Google Scholar] [CrossRef]
  51. Coulis, M.; Bernard, L.; Gerard, F.; Hinsinger, P.; Claude, P.; Villeneuve, M.; Blanchart, E. Endogeic Earthworms Modify Soil Phosphorus, Plant Growth and Interactions in a Legume–Cereal Intercrop. Plant Soil 2014, 379, 149–160. [Google Scholar] [CrossRef]
  52. Razafimbelo, T.M.; Albrecht, A.; Oliver, R.; Chevallier, T.; Chapuis-Lardy, L.; Feller, C. Aggregate Associated-C and Physical Protection in a Tropical Clayey Soil under Malagasy Conventional and No-Tillage Systems. Soil Tillage Res. 2008, 98, 140–149. [Google Scholar] [CrossRef]
  53. Parr, J.; Papendick, R.; Hornick, S.; Meyer, R. Soil Quality: Attributes and Relationship to Alternative and Sustainable Agriculture. Am. J. Altern. Agric. 1992, 7, 5–11. [Google Scholar] [CrossRef]
  54. Altieri, M.A. Agroecology: The Science of Sustainable Agriculture; Westview Press: Boulder, CO, USA, 1995; p. 182. [Google Scholar]
  55. Ranaivoson, L.; Naudin, K.; Ripoche, A.; Affholder, F.; Rabeharisoa, R.L.; Corbeels, M. Agro-Ecological Functions of Crop Residues under Conservation Agriculture. A Review. Agron. Sustain. Dev. 2017, 37, 26. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Experimental design. Treatments are abbreviated as follows: EW0: without earthworm inoculation, EW+: with earthworm inoculation, TR0: traditional tillage without rice residues restitution, TR+: traditional tillage with rice residues restitution, CA: conservation agriculture.
Figure 1. Experimental design. Treatments are abbreviated as follows: EW0: without earthworm inoculation, EW+: with earthworm inoculation, TR0: traditional tillage without rice residues restitution, TR+: traditional tillage with rice residues restitution, CA: conservation agriculture.
Agriculture 11 00060 g001
Figure 2. Earthworm density (ind m−2) at the end of the rainy period of the 4th cropping season according to treatments (earthworm inoculation, cropping system) (a) and proportion of juveniles to total earthworms (b). Error bars indicate standard errors. Letters indicate significant differences between EW0 and EW+ for different cropping systems according to Student t-test. See Figure 1 for legend.
Figure 2. Earthworm density (ind m−2) at the end of the rainy period of the 4th cropping season according to treatments (earthworm inoculation, cropping system) (a) and proportion of juveniles to total earthworms (b). Error bars indicate standard errors. Letters indicate significant differences between EW0 and EW+ for different cropping systems according to Student t-test. See Figure 1 for legend.
Agriculture 11 00060 g002
Figure 3. Plant properties at the 3rd cropping season according to treatments (earthworm inoculation, cropping system). Error bars indicate standard errors. Letters indicate significant differences between EW0 and EW+ for different cropping systems according to Student t-test. (A) Plant height, (B) Plant biomass, (C) Plant production and nutrition. Legend: see Figure 1. MAS = month after sowing.
Figure 3. Plant properties at the 3rd cropping season according to treatments (earthworm inoculation, cropping system). Error bars indicate standard errors. Letters indicate significant differences between EW0 and EW+ for different cropping systems according to Student t-test. (A) Plant height, (B) Plant biomass, (C) Plant production and nutrition. Legend: see Figure 1. MAS = month after sowing.
Agriculture 11 00060 g003
Table 1. Mean ± SE and one-way ANOVA statistics for soil variables measured at the beginning of the experiment (2013). Treatments are abbreviated as follows: EW0: treatments where earthworms will not be inoculated, EW+: treatments where earthworms will be inoculated. Significant p-values are in bold.
Table 1. Mean ± SE and one-way ANOVA statistics for soil variables measured at the beginning of the experiment (2013). Treatments are abbreviated as follows: EW0: treatments where earthworms will not be inoculated, EW+: treatments where earthworms will be inoculated. Significant p-values are in bold.
UnitsTreatmentsOne-Way ANOVA
EW0EW+Plot Position
Total C contentg kg−116.8 ± 0.614.3 ± 0.60.005
Total N contentg kg−11.84 ± 0.091.56 ± 0.05<0.001
Pi contentmg kg−110.65 ± 0.834.33 ± 0.55<0.001
Total C stockMg ha−119.3 ± 0.816.7 ± 0.70.028
Total N stockMg ha−12.11 ± 0.111.83 ± 0.080.011
Pi stockkg ha−112.1 ± 0.95.1 ± 0.7<0.001
Bulk densityg cm−31.14 ± 0.041.18 ± 0.040.058
Table 2. Mean ± SE for soil variables measured at the end of the 3rd cropping season (2016). Treatments are abbreviated as follows: TR0: traditional tillage without rice residues restitution, TR+: traditional tillage with rice residues restitution, CA: conservation agriculture system; EW0: without earthworm inoculation, EW+: with earthworm inoculation.
Table 2. Mean ± SE for soil variables measured at the end of the 3rd cropping season (2016). Treatments are abbreviated as follows: TR0: traditional tillage without rice residues restitution, TR+: traditional tillage with rice residues restitution, CA: conservation agriculture system; EW0: without earthworm inoculation, EW+: with earthworm inoculation.
Soil VariablesUnitsTreatments
TR0TR+CA
EW0EW+EW0EW+EW0EW+
Chemical soil variables
Total C content
g kg−123.7 ± 1.122.6 ± 0.624.3 ± 1.022.7 ± 1.423.6 ± 0.622.2 ± 1.0
NH4+ contentg kg−14.18 ± 0.834.60 ± 0.514.41 ± 0.434.14 ± 0.674.38 ± 0.545.07 ± 0.28
Pi contentmg kg−111.86 ± 1.905.70 ± 1.9313.36 ± 1.568.99 ± 2.6414.38 ± 0.909.73 ± 2.02
Total C stockMg ha−129.7 ± 1.828.1 ± 1.030.0 ± 1.428.6 ± 1.929.3 ± 1.328.1 ± 1.6
Total NH4+ stockMg ha−15.21 ± 0.995.70 ± 0.595.44 ± 0.515.20 ± 0.845.41 ± 0.646.42 ± 0.41
Pi stockkg ha−114.9 ± 2.57.1 ± 2.516.5 ± 2.011.3 ± 3.317.8 ± 0.912.4 ± 2.7
Physical soil variables
Bulk density
g cm−31.25 ± 0.021.24 ± 0.011.23 ± 0.011.26 ± 0.011.24 ± 0.031.27 ± 0.02
Macroaggregates%17.3 ± 0.730.0 ± 2.318.5 ± 0.824.2 ± 2.016.3 ± 1.530.2 ± 2.3
Table 3. p-Values of main effects (earthworm inoculation, cropping system, block) and their interaction and covariates (initial soil properties) on soil physico-chemical variables. ANCOVA and ANOVA were processed respectively for soil chemical and physical variables. Significant p-values are in bold.
Table 3. p-Values of main effects (earthworm inoculation, cropping system, block) and their interaction and covariates (initial soil properties) on soil physico-chemical variables. ANCOVA and ANOVA were processed respectively for soil chemical and physical variables. Significant p-values are in bold.
Factors Chemical Soil VariablesPhysical Soil Variables
Total C
Content
NH4+PiTotal CNH4+PiSoilBulk
ContentContentStockStockStockMacroaggregatesDensity
Main effects
Earthworm inoculation (E)0.3970.2160.2760.7080.2810.225<0.0010.214
Cropping system (S)0.8400.6430.6320.8270.4660.7040.3840.702
Block0.0290.0660.2550.0310.1280.1660.4900.077
Interaction
E*S0.8780.6200.9660.9630.4330.9540.0630.423
Covariates
Initial total C content0.9990.2290.801
Initial total N content0.2100.3250.417
Initial Pi content0.2890.7290.276
Initial total C stock 0.7560.1410.393
Initial total N stock 0.3360.3780.278
Initial Pi stock 0.2400.6020.148
Table 4. Mean ± SE of changes of soil properties with time (2013 to 2016). Legend: see Table 2.
Table 4. Mean ± SE of changes of soil properties with time (2013 to 2016). Legend: see Table 2.
Soil VariablesUnitsTreatments
TR0TR+CA
EW0EW+EW0EW+EW0EW+
∆C contentg kg−16.02 ± 0.767.48 ± 0.654.92 ± 1.307.40 ± 1.718.55 ± 1.499.77 ± 1.70
∆Pi contentmg kg−10.04 ± 0.711.29 ± 1.882.94 ± 1.353.63 ± 2.634.01 ± 1.285.66 ± 2.40
∆C stockMg ha−18.95 ± 1.4310.30 ± 0.747.78 ± 1.1910.35 ± 2.3012.54 ± 2.0214.12 ± 2.23
∆Pi stockkg ha−11.12 ± 1.291.91 ± 2.514.54 ± 1.664.92 ± 3.386.23 ± 1.313.20 ± 0.28
Table 5. Three-way ANOVA statistics (p-values) of changes of soil properties with time (2013 to 2016). Legend: see Table 2. Significant p-values are in bold.
Table 5. Three-way ANOVA statistics (p-values) of changes of soil properties with time (2013 to 2016). Legend: see Table 2. Significant p-values are in bold.
Soil VariablesUnitsThree-Way ANOVA
Earthworm Inoculation (E)Cropping System (S)BlockE*S
∆C contentg kg−10.1350.0890.3970.882
∆Pi contentmg kg−10.4650.1320.8340.970
∆C stockMg ha−10.2180.0560.4150.934
∆Pi stockkg ha−10.6590.1200.7090.971
Table 6. ANOVA p-values of main effects (earthworm inoculation, cropping system, block) and their interaction on plant variables at the 3rd cropping season (2016–2017). E*S means the interaction between Earthworm inoculation (E) and Cropping system (S).
Table 6. ANOVA p-values of main effects (earthworm inoculation, cropping system, block) and their interaction on plant variables at the 3rd cropping season (2016–2017). E*S means the interaction between Earthworm inoculation (E) and Cropping system (S).
Plant VariablesThree-Way ANOVA
Main EffectsInteraction
Earthworm Inoculation (E)Cropping System (S)BlockE*S
Plant growth
Height 1MAS0.3060.0280.0860.771
Height 2MAS0.0280.3510.7130.375
Height 3MAS0.0080.4030.9830.477
Shoot biomass0.1260.0840.2680.280
Grain biomass0.0470.3590.6310.813
Shoot biomass0.0600.3000.5350.584
Root biomass0.9840.0940.7610.828
Total biomass0.1790.1200.8470.673
Shoot:root ratio0.0910.3520.2620.997
Plant nutrition
Shoot N amount0.6800.3100.1840.314
Grain N amount0.0370.4210.3450.843
Aboveground N amount0.0630.4860.2710.714
Shoot P amount0.1400.5470.3110.671
Grain P amount0.0560.3920.0650.848
Aboveground P amount0.1670.6370.0700.781
Plant production
Grain yield0.0400.2810.6680.802
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ratsiatosika, O.; Razafindrakoto, M.; Razafimbelo, T.; Rabenarivo, M.; Becquer, T.; Bernard, L.; Trap, J.; Blanchart, E. Earthworm Inoculation Improves Upland Rice Crop Yield and Other Agrosystem Services in Madagascar. Agriculture 2021, 11, 60. https://0-doi-org.brum.beds.ac.uk/10.3390/agriculture11010060

AMA Style

Ratsiatosika O, Razafindrakoto M, Razafimbelo T, Rabenarivo M, Becquer T, Bernard L, Trap J, Blanchart E. Earthworm Inoculation Improves Upland Rice Crop Yield and Other Agrosystem Services in Madagascar. Agriculture. 2021; 11(1):60. https://0-doi-org.brum.beds.ac.uk/10.3390/agriculture11010060

Chicago/Turabian Style

Ratsiatosika, Onja, Malalatiana Razafindrakoto, Tantely Razafimbelo, Michel Rabenarivo, Thierry Becquer, Laetitia Bernard, Jean Trap, and Eric Blanchart. 2021. "Earthworm Inoculation Improves Upland Rice Crop Yield and Other Agrosystem Services in Madagascar" Agriculture 11, no. 1: 60. https://0-doi-org.brum.beds.ac.uk/10.3390/agriculture11010060

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop