Abstract

Our objective is to study interaction between physical and chemical properties of soils and their earthworm community characteristics in different areas irrigated by wastewaters and well waters. The fields have different topography and agricultural practices conditions and are located in two regions of Batna department (Eastern Algeria). Both regions are characterized by a semiarid climate with cold winters and Calcisol soils. Nine fields were subject of this study. Three of these fields are located in Ouled Si Slimane region whose irrigation is effectuated by natural waters of Kochbi effluent. The other six fields are located at edges of Wed El Gourzi, effluent from Batna city, and partially treated through water treatment station. The best rates of water saturation and infiltration as well as abundance of earthworms were recorded at sites characterized by irrigation with wastewaters downstream of El Gourzi effluent. PCA characterizes two major groups: a group of hydrodynamic infiltration parameters and structural index stability of soil, explained by fields irrigated with wastewaters downstream of El Gourzi effluent. This group includes chemical characteristics: pH and electric conductivity. The second group is the characteristics of earthworms and includes organic matter content, active limestone levels, and Shannon Biodiversity Index.

1. Introduction

Treated or untreated urban wastewater has been used commonly for agricultural irrigation in arid and semiarid regions of the world. According to the estimations, at least 20 million hectares of agricultural land worldwide is irrigated with treated and untreated wastewaters [1]. In areas suffering from water shortage, such as arid and semiarid areas, where insufficiency and water scarcity inevitably pose problems, especially from the point of view of meeting quantitative needs of irrigation in agriculture, the wastewater recycling as a complementary water resource is interesting for cultivated soils to solve the problem posed in agriculture by water resources insufficiency in these zones [2]. Its use has increased recently because there are inadequate freshwater resources. The population suffering from water scarcity is presently 11% of the total worldwide population; it is estimated that the population with adequate water will be 38% in 2025 [3].

Reuse of raw or treated agricultural wastewater is a widespread practice downstream of urban centers and mostly in areas affected by water shortages [4]. Although proscribed, this practice is appreciated by farmers because wastewater is a regularly available and abundant resource and contains the fertilizer needed to improve soil properties and yield crops [5, 6].

Urban wastewater contains higher levels of organic matter, nutrients, and pollutants (heavy metals and suspended solids) compared to fresh water. Although wastewater application provides positive effects on soil properties and crop productivity because of its high organic matter and macro- and micronutrient contents, the pollutants in wastewater may cause some problems to soil and crops [7]. Nevertheless, wastewater compounds in particular affect soil porosity and hence hydrological properties [8]. Nadav et al. [9] indicated that the physicochemical properties of soils were altered with treated wastewater irrigation, because long-term wastewater application caused the accumulation of organic matter in soil. High organic matter in wastewater is cement for the improvement of soil aggregates. Therefore, lower bulk density and higher infiltration and water retention have been obtained under the wastewater irrigation conditions. However, suspended solids in wastewater negatively affect the soil porosity.

Earthworm’s activities of a mechanical or metabolic nature allow transforming organic matter, recycling nutrients, and transferring matter and energy between different horizons. Their importance in mixing soil horizons has been underlined by many authors [1013]. Earthworms have a very important role for soil fertility; they essentially favor the formation of clay-humic complexes; anecic worms together feed on organic and mineral matter, increase stability of aggregates, and make it possible to unpack compacted soils [14]. Above all, they create loose structure elements in soil, possessing high porosity, which increases useful reserves and allows good aeration; this structure is very favorable to germination and nodulation of clover and increases yield of certain crops [15]. Roots use earthworm galleries to grow [14]. Earthworms also modify the microbial communities through digestion, stimulation, and dispersion in casts. Consequently, changes in activities of earthworm communities, as a result of soil management practices, can also be used as indicators of soil fertility and quality [16]. Despite the importance of this fauna and work done on these organisms, it is not important in Algeria [17]. In Algeria, work on biodiversity of earthworms is still insufficient. This very varied biogeographic space in terms of climate, soil, and vegetation from the littoral to the desert could reveal a great diversity of earthworms with certainly species very adapted to the drought. Studies on this subject are difficult: on one hand, identification and classification of these organisms remain difficult due to lack of qualified taxonomists [18] and, on the other hand, study of earthworms is not obvious due to several constraints related to nature of soil and the complexity of these organisms [19].

In this complex set of contexts, we have found it interesting to contribute to studying the interaction between physical and physicochemical properties of soil and the properties of earthworm community in El Madher region, where soils are cultivated and irrigated by the partially treated waters of El Gourzi effluent and are in incessant extension. We have chosen Ouled Si Slimane area as a control field, which is irrigated by water from natural sources. El Madher and Ouled Si Slimane areas are located in Batna department, characterized by a semiarid climate with cold winters and Calcisol soils according to WRB in 1999.

2. Materials and Methods

2.1. Presentation of El Madher Region

El Madher town is located in the northeast part of Batna department at a distance of 23 km from Batna town (Figure 1). The geographical coordinates of this region are as follows: latitude: 35°37′52′′; Longitude: 6°22′09′′. It is at an altitude of 1100 meters.

The geomorphology of this region is related to brittle tectonics, where a collapse zone corresponds to a plain and a high zone corresponds to the mounts: Dj. Bou Arif (1744 m), Koudiat Tfouda in Dj. Tafraout (1080 m), and Dj. Sarif (1744 m). The hydrographic network is mainly represented by Gourzi effluent or El Madher river, which springs from Batna watershed and empties into the Chott Gadine, with a multitude of Chaabas and Talwegs.

Precipitation is very low distributed over a period of 30 years (1985–2015), with total annual average precipitation of 319,28 mm/year. The average monthly temperature range is from 5.23°C in January to 25.66°C in July. The Gaussen Ombrothermic diagram shows a relatively long dry period, spreading from May to October. The study area is classified according to Emberger’s Climogram in semiarid bioclimatic zone with cool winter and hot summer.

2.2. Presentation of Ouled Si Slimane Region

Ouled Si Slimane town is located in the western part of Batna department at a distance of 90 km (Figure 1). The chief town of Daïra is limited by the following Lambert coordinates: 35°37′N for latitude North and 5°38′ for longitude East, with an average altitude of 760 m.

The main geological structure of Ouled Si Slimane region is divided into three strata as follows:(i)The mountainous region consists mainly of calcareous rocks of the Cretaceous.(ii)The plains consist mainly of calcareous marls of lower Miocene.(iii)There are also slopes and Wadis (rivers or effluents). Those are alluviums, regs, and terraces of the Quaternary constituted principally of clays, silt, and pebbles.Hydrography is represented by an effluents network all oriented from north to south and from north to east of the region. This network is practically seasonal, dry in summer and torrential in winter. It includes Wed Laïoune and Wed Kochbi effluents.

Precipitation is very low distributed over a period of 25 years, with total annual average precipitation of 318,12 mm/year. The Gaussen Ombrothermic diagram shows a relatively long dry period, spreading from May to October. The study area is classified according to Emberger’s Climogram in semiarid bioclimatic zone with cold winter and hot summer.

2.3. Selection and Description of Sampling Sites

We selected in our first study area six different sites with different peculiarities of topography (upstream and downstream), vegetation, agricultural practice, and irrigation by wastewaters from El Gourzi effluent for more than 10 years. These sites are distributed along the effluent upstream at an altitude of 924 m to downstream at an altitude of 868 m for a distance of about 25 km. El Gourzi effluent is the main collector of sewage network of the city of Batna as well as rainwater. It is an open effluent that crosses this town, with a flow that varies according to seasons. Before leaving the city, it passes through industrial zone, where it collects, in addition to urban waste, all industrial waste [20]. Only a part of these discharges (15,000–2,000 m3 per day) are treated in treatment station before joining the effluent [21].

Maalem et Ghanem [22] studied the water quality of El Gourzi effluent at 5 upstream to downstream points; he showed that these wastewaters are characterized by higher values ​​of pH, EC, MES, BOD5, and COD upstream than downstream (Table 1). This was also proven in the study of Tamrabet et al. [20]. According to the same study [20] (Table 2), this effluent is generally more charged upstream than downstream by nitrates, ammonium, orthophosphates, iron, copper, zinc, manganese, fecal coliforms, and fecal Streptococcus. The process of evolution from upstream to downstream direction acts by current force on water flow itself on suspended matter flow, mineral and organic matter, and that of living organisms [23]. In conditions of our study area, problems of phytotoxicity by Trace Element Metals do not arise for the moment because of their low content on one side and physicochemical characteristics of soil on the other side. The soil of El Madher region is of clay loam nature; according to Marschner [24], these conditions favor the precipitation of TEM. However, Alouni [25] mentions that the epidemiological risks of pathogens such as Salmonella and Vibrio cholerae are zero. Only the risks related to fecal coliforms and fecal streptococci are latent in wastewaters.

Three other sites in second studied area are irrigated with healthy potable and unpolluted waters from Kochbi and differ only in conditions of vegetation and cultural practices. El Kochbi is the effluent from natural water source of Kochbi. These sites are distributed in two localities of Guerza and Teniet El Abed at an altitude oscillating between 780 and 849 m for a distance of about 20 km.

The sites are all characterized by fine granulometric size textures (Table 3). The sites are characterized by their cultural practices (cereal cultivation and arboriculture) and tillage; irrigation takes place according to climatic demands of crops. The input of fertilizers is made as basic fertilizer every two years generally with plowing, in form of NPK. The sites chosen as control are the borders of effluent under natural vegetation and without any practice; they are influenced by their proximity to effluents and are wet by capillary ascension. We carried out 5 sampling points of soil and earthworms in each of the sites. In total, 45 points were studied (Table 4) in 2014, 2015, and 2016 in spring period of March, April, and May. This period is the most suitable for removing earthworms because the soil is moist, and spring temperatures favor worm activity.

2.4. Methods Sampling

The method used is manual sorting [1012]. It is a physical method of extracting earthworms. Each sampling point consists of a single sample using a shovel with a volume of soil of 30 × 30 × 30 cm. Worms were collected at the same time as soil and stored in vials with 4% formalin. Soil sampling was carried out after sorting the earthworms and deposited in labeled and numbered polyethylene bags and taken to the laboratory.

2.5. Methods of Soil Physical Properties Analysis

Determination of different particle size soil fractions was performed after organic matter destruction and particles dispersion; the clay, the fine loam, and the coarse silt content were estimated by Robinson pipette, fine sand, and coarse sand by sieving.

Soil moisture at time of earthworms’ removal was determined after passing in oven set at 105°C of a soil moist weight for 24 hours. It is estimated by the following formula:% is the soil moisture at time of earthworms removal; is soil moist weight; is soil dry weight.

Soil moisture at maximum retention capacity (Cr%) was estimated by submersion of soil samples by water until saturation. Steaming is at 105°C.

The method measuring soil bulk density (Da, g/cm3) consisted in using metal cylinders of known volume (126,60 cm3). Porosity (P%) was calculated by bulk density measurements and soil real density of 2.65 g/cm3.

Hénin’s structural instability index (IS) according to Mussy and Soutter [26] is by definition proportional to soil sensitivity phenomena of bursting on one side and of swelling dispersion on the other. Three samples of 10 g of a soil sample of 300 g are taken without the destruction of the organic matter, previously sieved to dryness at 2 mm. Two of them are, respectively, treated with 5 cm3 of alcohol at 95° and 5 cm3 of benzene. After 5 minutes of imbibition, the three samples are suddenly immersed in water and then after half an hour of contact sieved to 0.2 mm, which makes it possible to determine the terms AgA, AgB, and AgE. IS was expressed by the following relation:where AgA, AgB, and AgE are fractions determined by weight which are collected after sieving with a tamis of 0,2 mm. These fractions remain after immersion in, respectively, alcohol, benzene, and water. Sg is a fraction of sand included in AgA, AgB, and AgE fractions.

The filtration rate K (h) at saturation was estimated from Darcy’s sense by the method of a soil maintained in a permanent water flow. K (h) was carried out by open PVC cylindrical tube at both ends, having a diameter of 9 cm and a height of 40 cm, and is gently pressed into soil layer and carefully removed without being disturbed. A piece of tulle is put on using an elastic ring at the bottom of the tube. Water soil infiltrating in a beaker is collected every 5 minutes in a test tube to determine volume. This happens until reaching a constant volume during 1 hour. K (h) sat. (mm/h) is determined for the water heights of 0.05, 0.3, 0.6, and 1 kPa and calculated by (3) of Darcy law cited by Mussy and Soutter [26].where is height of soil column by cm; is percolated water for 1 hour of infiltration by cm3; is height of water column by cm; is internal section of tube by cm2.

2.6. Methods of Soils Chemical Properties Analysis

The pH and electrical conductivity (EC) were measured by direct reading on a pH meter and conductivity meter in a suspension with a soil/water ratio of 1 : 2.5 and 1 : 5, respectively [27].

Determination of organic matter content (OM%) is performed using the Walkley and Black method, which is based on the oxidation of carbon with potassium dichromate in strongly acid media [28].

Total CaCO3% is determined by volumetric method of the “Bernard Calcimeter,” while active CaCO3% is determined by contact with a specific extraction reagent, “ammonium oxalate,” at 0.2 N.

2.7. Methods Studying Earthworms

The determination of taxa concerns only adult worms. In this study, we tried to use much more external characters of worms collected using internal anatomical features for certain taxa to determine species in question. The external and internal characters used are those explained by Bouché and Bachelier [10, 11]. The key chosen for taxa determination is modified by Blakemore in 2007 [29]. For classification of ecological categories, we used that of Bouché [10, 30] which was used by Bazri [31] in his study of the Northern Algerian earthworm population.

Abundance is expressed as point abundance, that is, total number of earthworms present in a sampling point, and is expressed as number of individuals per m2.

Earthworm biomass is expressed as population biomass taken and expressed by individual biomass (weight of each worm) and worm point biomass (total weight of worms taken at a sampling point).

To estimate diversity of earthworm community for each site, we used Shannon diversity index derived from a function established by Shannon and Wiener which became the Shannon diversity index.

Value of Shannon diversity index varies between 0.5 and 4.5 bits and is given by the following formula [32, 33]:where (see (5)) is number of species individuals in relation to total number of individuals identified :This index is independent of sample size and takes into account distribution of individuals number per species [34].

3. Results and Discussion

3.1. Soil Compartment
3.1.1. Physical Properties of Soil

The means of soil moisture (H%) at moment of sampling of earthworm are in Table 5. These close rates of soil moisture among the nine sites seem to be more favorable to the biological activity of earthworms. These moisture levels are due to spring irrigation or precipitation. Average of highest water retention rate (Cr%) was recorded at S1, S4, and S7 at border of effluents. Overall, mean water retention rates in study sites appear to be closely related to fine particle contents of study fields, especially at the edges of effluents downstream, where they are deposited after training by the flow of water. Meanwhile, Mojid and Wyseure [35] determined enhanced water retention capacity values in the wastewater irrigation conditions compared to fresh water irrigation.

Porosity is highest in S2 characterized by wastewaters downstream and by cereal crops, while it was lowest at S7 at edges of natural water Kochbi effluent (Figure 2). Effectively, Gharaibeh et al. [36] noted that irrigation treated for long periods with wastewaters resulted in slight decrease in bulk density compared to the control. Many researches obtained lower bulk density or higher porosity values under wastewater irrigation [35, 37, 38]. Moreover, elevation of dissolved organic matter in treated wastewaters is reported to cause/increase clay dispersion by reducing the differential viscosity of clay suspensions and decreasing the attraction forces between clay particles [39]. Wang et al. [40] indicated that wastewater irrigation caused a slight increase in soil compaction. Particularly high-suspended solid concentration in wastewater may increase the bulk density, while lower concentrations may not significantly affect it [41]. These values reflect, on one hand, the plowing practices in the study sites and, on the other hand, a less favorable effect of the accumulation of organic matter in the fields located upstream of the effluent.

Structural stability is better in S6 characterized by irrigation with wastewaters upstream and alternations of vegetable crops and cereal crops (Figure 3). It is low in S8 with irrigation by natural waters and arboriculture. Indeed, Kirkham [42] notes that inputs of sludge in successive 4 years have raised organic matter content of soil of the first 15 cm from 1.2 to 2.4%. Input of organic matter and calcium ions play a favorable role in cements stability, which causes improvement of soil particles. This same author underlined an improvement of unfavorable structures for clayey soils following application of sludge. Also, Gharaibeh et al. [36] noticed that aggregate stability (AS%) results revealed that irrigation with treated wastewaters significantly increased the percentage of stable soil aggregates compared to the control. Miller and Kemper [43] observed increases in water-stable aggregates at least for one growing season following alfalfa incorporation, after which there was a decrease. They attributed this increase to the production of cementing substances through microbial activity by fungal and actinomycete mycelia which provides substrates to stabilize soil aggregates. However, Vogeler [38] attributed higher AS value to the higher total carbon content or to the composition of soil organic matter in areas irrigated for long periods with treated wastewaters compared with the control block.

The best rates of filtration at saturation Ks (h) for 0,05 kPa suction were recorded in S2 and S3 sites characterized by wastewaters downstream of El Gourzi effluent (Figure 4). But globally for this pressure, Ks is more important in the irrigated wastewaters. For the other pressures, difference is less remarkable. Rapid infiltration capacity of a soil is estimated by macroporosity. Rabbi et al. [44] mention that percentages of clay, bulk density, and porosity have a strong influence on Ks. These results are in agreement with Bardhan et al. [45] who reported a decrease in the infiltration rate due to clogging of soil pores by suspended materials present in the treated wastewaters. Previous studies with similar soil texture (clay) of the current study reached the same conclusion that application of wastewater resulted in pore clogging, which leads to reduced soil porosity and subsequently a decrease in the soil infiltrability [46].

Meanwhile, Morel et al. [47] showed that permeability of a soil increases significantly in plots enriched with sludge, as organic matter richness of sludge improves water balance and increases reserve of useful water; this favors processes of stabilization of soil aggregates and in particular gives soil a better permeability in relation to a more stable structure. There is a close relationship between the pore size distribution and soil water content due to the fact that macropores control the aeration and drainage, mesopores control the water conductivity, micropores control the water retention and the most available water for plants [48].

3.1.2. Chemical Soil Properties

The pH (Table 5) in study fields shows moderate alkalinity in sites irrigated with natural waters, as well as fields irrigated by wastewater upstream. The pH levels are more alkaline in sites downstream of wastewater effluent than in the sites with upstream wastewaters. Standard deviations between pH values of sampling points of sites often irrigated by upstream wastewater are much greater. Schipper et al. [15] indicate that soil pH increases as a result of a long period of irrigation with wastewater. They attribute this increase to chemical composition of cations effluent such as Na, Ca, and Mg. The pH of soil irrigated with wastewater decreases following oxidation of organic compounds and nitrification of ammonium [4951].

Electrical conductivity values are the lowest at sampling points of S5, characterized by irrigation by wastewater upstream of El Gourzi effluent. Standard deviation values are much higher at sampling points of S1, as previously described. Contrary to the results mentioned in the work of Maalem and Ghanem [22] on the values of El Gourzi effluent’s wastewaters, which are higher upstream than downstream, the soil has higher CE values downstream than upstream. This could be due to the organic matter, which acts as a buffer. In a general way, according to standards presented by DIAEA/DRHA/SEEN in El Oumlouki et al. [52], EC of various treatments soils of our experimentation remains weak.

The organic matter is lowest at S8 with well water irrigation and arboriculture (Figure 5). It is highest at S4 upstream of El Gourzi effluent border. Indeed, Tamrabet et al. [6] reported that a comparison of control soil averages and irrigated soils by wastewater shows that irrigation with this wastewater has a very significant effect on OM improvement. However, many other studies showed increase in the OM% within creased period of treated wastewaters irrigation [38, 53, 54]. This is the case in our study. Moreover, high clay content of soil may physically protect the OM from the decomposition.

The average percentage rate of total CaCO3% is lowest at S3 downstream of El Gourzi effluent. It is highest at Kochbi natural water effluent border. These values are generally strongly to very strongly calcareous according to grading scale proposed by GEPPA in Baize [27]. The limestone levels appear to be related to soil pH.

3.2. Earthworm Compartment (Table 6)

The total abundance and biomass average of earthworms at different sampling points of study sites are the lowest in S9 with irrigation by natural waters of El Kochbi effluent and cereal crop. They were the highest in S5 with irrigation by El Gourzi wastewater upstream and cereal growing (Figure 6). This abundance was represented by 54% of adults and 46% of juvenile worms (Figure 7). Our results explain a negative effect of plowing on the abundance and biomass of earthworms. Contrariwise, the irrigation with wastewaters improves significantly their abundance and biomass. In Bazri et al.’s [17] study, the density of earthworm in eastern Algeria from the coast to the desert was, respectively, 6.00 ± 1.41 to 29.60 ± 11.83 individuals/m2 and 0.28 ± 0.39 to 13.13 ± 7.94 g/m2. These results are almost similar to those of Omodeo and Martinucci [55] in Northern Algeria who found earthworm densities ranging from 11.0 to 12.7 individuals/m2 and biomass ranging from 1.25 to 3.0 g/m2. Edwards and Bohlen [56] explained that soils with low organic matter contents usually do not support high densities of earthworms. The highest values are usually found in fertilized pastures and the lowest ones are in acid or arid soils [12, 57].

Taxonomic study was carried out on a total of 594 adult worms. The taxonomic key and nomenclature quoted by Rougerie et al. [18] allowed classification of 7 species of Lumbricidae: Aporrectodea trapezoides (Dugés, 1828), Aporrectodea rosea (Savigny, 1826), Aporrectodea caliginosa (Savigny, 1826), Allolobophora molleri (Rosa, 1889), Octodrilus complanatus (Dugés, 1828), Eiseniella tetraedra tetraedra (Savigny, 1826), and Proctodrilus antipai antipai (Michaelsen, 1891). Percentage of overall specific abundance shows dominance of Aporrectodea trapezoides, followed by Aporrectodea rosea (Figure 8).

Ap. trapezoides is present in greater quantity (13.4 individuals/m2) in S7 characterized by nonpractical plowing and natural vegetation at Kochbi edge where there is healthy natural water (Figure 9). However, Ap. rosea is present in greater amounts (13.0 and 12.0 individuals/m2) in S4 and S5 at edges of effluent and irrigated by wastewaters upstream. Pr. antipai is the most dominant species in S6, characterized by wastewater irrigation upstream of the effluent and arboriculture. On biomass, Ap. trapezoides is highest (12.48 g/m2) in S7 at edges of natural water effluent. Ap. rosea is present on a biomass of 11.7 g/m2 in S5 irrigated by wastewaters and characterized by cereal cultivation. These results allow us to deduce that wastewaters have intense effect on development of earthworm’s biodiversity. Species Ap. rosea was the most concerned by this improvement.

Bazri et al. [17] noted in their study that it is curious and hard to explain that Ap. trapezoides, the most common and dominant species, was not previously quoted by other authors working on Algerian lumbricoid fauna. Possible explanations include identification problems in this complex group of species, recent introduction and posterior expansion, or older introduction as suggested by the large geographic distribution of the species. In the semiarid regions, the taxon A. molleri is important; it locates preferentially in wet points (notably at the edge of effluents). In mountains of semiarid regions, the species O. complanatus dominates. As for Ap. rosea, it is the only one observed in arid bioclimatic stage at points where there is sufficient water.

Among the 7 species inventoried (Figure 10), the endogeic species accounted for 62% (Ap. rosea, A. molleri, P. antipai, and Ap. caliginosa), epigeic species for 3% (E. tetraedra), and anecic species for 22% (O. complanatus and Ap. trapezoides); the latter may be considered as anecic, endogeic, or endoanecic because it varies according to the strain.

Bouché [10] separated earthworms into three categories, based on morphological and behavioral characteristics. Epigeic species are consumer litter living and feeding on or surface soil. Anecic earthworms live in permanent vertical burrows within the soil and may emerge to feed on surface layer; endogeic species live in temporary horizontal burrows and feed on the soil. This species is geophagous, since it gains its nutrients by eating the soil and the green morph is characterized by Bouché [10] as more epigeic.

Endogeic abundance is much higher in S4 at edges upstream of the El Gourzi effluent, which is without cultural practices (Figure 11). However, S7 at edges of El Kochbi effluent from natural waters, without cultural practices, is better represented by an anecic population. The epigeic population is important only on S4 and S6 upstream of the wastewater effluent. Endogeic earthworms are a major component of soil fauna communities in most natural ecosystems of the humid tropics [58].

The Shannon diversity index of different study sites is the lowest at S7 corresponding to edges of Kochbi effluent. It is the highest at S4 characterized by wastewaters upstream of El Gourzi effluent and any cultural practices. This result allows us to deduce that wastewater loaded with organic matter promotes considerably the inheritance of the population of earthworms.

3.3. ANOVA

Statistical analysis with the Newman-Keuls multiple comparison posttest at ≤ 0,05% (Table 5) revealed a highly significant effect of different irrigation practices and cultural practices on soil’s physical parameters (Is, Cr, P%, K 0.05, K 0.3, K 0.6, and K 1). Water retention, infiltration, and hydraulic conductivity are among the soil’s hydraulic properties affected by soil porosity and pore size distribution [40, 59]. Spatial variability of soil’s physical properties within or among agricultural fields is inherent in nature due to geologic and pedologic soil forming factors, but some of the variability may be induced by tillage and other management practices. These factors interact with each other across spatial and temporal scales and are further modified locally by erosion and deposition processes [60]. In the regions with a cool climate, soils are exposed to freeze-thaw cycles, especially in the spring period. Aggregation and therefore soil structure may be either positively or negatively affected by freeze-thaw cycles [61]. Therefore, the impacts of wastewater irrigation on main soil properties in agricultural areas under cool climate conditions may be different.

ANOVA reveals highlight at very highly significant effect of different irrigation practices and cultural practices on studied chemical parameters of soil (MO, CE, CaCO3 tot., CaCO3 act., and pH). Effectively, the soil’s physical properties are associated with nutrient applied and environmental soil science availability, solute and pollutant movement, microbial activity, and soil organic matter stabilization [62].

On parameters of earthworm characteristics (Table 6), ANOVA revealed a significant effect of study sites differences except for Ap. caliginosa and E. tetraedra tetraedra species’ abundance. The significant effect of sites has been well revealed on endogeic and anecic abundance. However, there is no significant effect of different irrigation and cropping practices on epigeic abundance.

3.4. PCA

In order to study interaction between physicochemical soil parameters and earthworm communities of studied sites, Principal Component Analysis (PCA) was used. This procedure makes it possible to group or distribute the sampling sites around principal axes in function of the physicochemical and earthworm parameters, thus facilitating observation of possible links between variables and places where they are most represented (Figure 12).

Contribution of principal axes to total variation is 41.90% for axis 1 and 19.96% for axis 2, which makes a total of 61.86%, which is well acceptable. On variables distribution graph, both axes, 1 and 2, contrast two groups of variables: biological parameters of earthworms (total biomass, total abundance, and adult worms’ abundance) and physical infiltration parameters (infiltration at saturation, K 0.05 kPa, K 0.3 kPa, K 0.6 kPa, and K 1 kPa). The first group of biological variables is represented by S5 site, characterized by irrigations with wastewaters upstream of El Gourzi effluent and by cereal crop (wheat/alfalfa). The group of infiltration characteristics is represented by S2 and S3 sites. These sites are irrigated with downstream wastewaters and characterized by cereal and arboreal crops. This group includes physicochemical characteristics: pH and EC. Axis 1 is intensely correlated with structural stability index and water retention capacity. It divides the characteristics of earthworms studied into two subgroups. One of them is well represented by S4 and S6 with irrigation by wastewater at upstream. This group includes OM and active CaCO3 as well as Shannon Biodiversity Index.

Francis and Fraser [63] and Capowiez et al. [64] reported that galleries excavated by terrestrial bioturbation activity contributed to water transfers. Also, Bottinelli et al. [65] and Pérès et al. [66] reported that earthworms contribute to improving soil porosity. Indeed, the PCA groups the S4, S5, and S6 sites with abundance and biomass of earthworms and soil porosity. For his part, Supersperg [67] found that annual applications of liquid sludge cause an increase in compactness of a heavy soil by clogging and a decrease in pores volume.

By this PCA, it can be deduced that very high organic charge of partially treated waters has negative repercussions on hydraulic conductivity if this load exceeds certain doses, because the PCA rallied the hydraulic conductivity at saturation by different pressures applied to S1, S2, and S3 sites downstream of El Gourzi effluent, as well as S8 and S9 sites irrigated with natural waters of El Kochbi effluent. This is contradictory to the results of Minhas and Samra [68] and Ababsa et al. [69] who have worked on the effect of earthworms on hydraulic conductivity in soils irrigated by wastewaters. Meanwhile, our findings concur with the works of some authors, Wang et al., Viviani and Iovino, and Molahoseini [40, 46, 70], who have to work on the same axis of research.

This opposition of results can be explained by more or less rapid duration of organic matter evolution of wastewaters. This change in state of OM from pores plugging form to soil stabilizing and structuring form is closely related to soil physicochemical parameters (pH, active CaCO3%, and EC). Indeed, Callot and Dupuis [71] assume the role of active CaCO3% in protection and conservation of organic matter in carbonate soils. Active CaCO3% thus acts as a protector against biological degradation of organic matter of these soils to make it a stable phase not sensitive to decomposition [72]. All the soils being studied are carbonate soils; their active CaCO3% should play an important role in humification process by rapidly sequestering lignified organic compounds as inherited or residual humus [28].

4. Conclusion

Our objective was to study interaction between earthworm community characteristics and physicochemical properties of soils characterized by irrigation with wastewaters and well waters in conditions of different topography (upstream and downstream) and cropping practices (natural vegetation, cereal cultivation, and arboriculture) in two regions of Batna department.

The results show that the maximum water retention capacity is better explained by the fine texture of the soil at the effluent border. The porosity and rates of infiltration at saturation are higher in sites characterized by irrigation with wastewaters downstream of effluent. This was explained by clogging of soil pores by suspended materials present in the wastewaters. However, structural stability is better in site characterized by irrigation with wastewaters upstream of El Gourzi effluent and cereal crop. This is due, on one hand, to input of organic matter and calcium ions which play the role of cements and cause improvement of particles soil aggregation according to Kirkham [42]. On the other hand, at cereal crop, that effect increases production of cementing substances through microbial activity by fungal and actinomycete mycelia that provide substrates to stabilize soil aggregates according to Miller and Kemper [43].

The pH levels are more alkaline in downstream sites than upstream sites of wastewater effluent. The pH of soils irrigated with wastewater decreases following oxidation of organic compounds and nitrification of ammonium. Electrical conductivity values are the lowest at sampling points of site, characterized by irrigation with wastewater upstream of El Gourzi effluent and cereal crop. This could probably be due to the level of organic matter which acts as a buffer. The organic matter is highest at site upstream of El Gourzi effluent border characterized by natural vegetation, which is certainly explained by an important load of organic matter from the wastewater upstream of the effluent, and high clay content of soil may physically protect the OM from the decomposition. The average percentage rate of total CaCO3% is higher at Kochbi natural water effluent border, while it showed low values at downstream sites. However, the active CaCO3% rate is the lowest at downstream site of El Gourzi effluent border. But it presents the highest values at S6 site irrigated by wastewaters upstream of effluent. This result appears to be related to soil pH.

The total abundance and biomass average of earthworms at different sampling points of study sites are the lowest at site irrigated by natural waters of El Kochbi effluent which is cultivated by cereal. On the contrary, they are the highest at upstream site irrigated by El Gourzi wastewater and cereal growing. This result explains that abundance and biomass of earthworms are related to wastewaters irrigation. Earthworms are much more present where irrigation water is heavily loaded with organic materials.

Taxonomic study allowed classification of 7 species of Lumbricidae. Percentage of overall specific abundance shows dominance of Aporrectodea trapezoides, followed by Aporrectodea rosea. The result shows that wastewaters have intense effect on development of earthworm’s biodiversity. Species Ap. rosea was the most concerned by this improvement. Among the 7 species inventoried, the endogeic earthworms represent the majority. Endogeic abundance is much higher in site at edges upstream of the El Gourzi effluent, which is without cultural practices. However, site at edges of El Kochbi effluent from natural waters, without cultural practices, is better represented by an anecic population. The Shannon diversity index of different study sites shows that wastewater promotes considerably the inheritance of the population of earthworms.

ANOVA revealed a highly significant effect of different irrigation types as well as different cultivation practices on all physicochemical properties of soil and on endogeic and anecic population. This was not the case for the epigeic population.

The PCA allowed characterizing two major groups. The first group is defined by sites that are explained by the hydrodynamic infiltration. These sites are irrigated by wastewater; they are located downstream of effluent (lowest load of organic matter) and characterized by cereal and arboreal crops. This group includes other chemical characteristics: pH and EC. The second group formed by upstream sites irrigated by wastewaters encompasses the characteristics of earthworms explained by highest load of organic matter. It includes chemical characteristics: OM and active limestone. It also includes the Shannon Biodiversity Index. We explained this by active CaCO3%, which acts as a protector against biological degradation of organic matter of these soils to make it a stable phase not sensitive to decomposition.

Overall, the results of the PCA show us that the effect of wastewater is positive on improvement of abundance and biomass and biodiversity of earthworms population. On the contrary, the wastewaters are unfavorable for the infiltration physical parameters of water into soil. This will surely have long-term repercussions on the physical and biological parameters of soil and will induce a new organization of macro- and mesoporosity due to improvement of biological parameters of earthworms.

Finally, we propose complementing these results in the field and in the laboratory to see the evolution of infiltration parameters of soil water in parallel with the dynamics of earthworm community in soil irrigated by wastewaters. The biodiversity of earthworms may induce an improvement in infiltration and reduce the problems of clogging pores and solve other problems that wastewaters can introduce into soils.

Abbreviations

Cr:Water retention capacity
Is:Structural instability index
Da:Bulk density
:Porosity
 (h) sat.:Saturated hydraulic conductivity at potential difference
CE:Electric conductivity
MO:Organic matter
TEM:Trace Element Metals
MES:Suspended matter
BOD5:5-Day biological oxygen demand
COD:Chemical oxygen demand
Abd:Abundance
Juv.:Juvenile
Adt.:Adult
Biom.:Biomass
Endo.:Endogeic
Epig.:Epigeic
Anec.:Anecic
Ap t: Aporrectodea trapezoides
Ap r: Aporrectodea rosea
Ap c: Aporrectodea caliginosa
Am: Allolobophora molleri
P ant: Proctodrilus antipai antipai
Oc: Octodrilus complanatus
E t: Eiseniella tetraedra tetraedra
SDI:Shannon Biodiversity Index
Dj:Djebel or mountain
Wed:effluent.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.