Next Article in Journal
Allelopathic Potential of Plant Aqueous Mixtures on Euphorbia heterophylla
Previous Article in Journal
Diversity and Management Strategies of Plant Parasitic Nematodes in Moroccan Organic Farming and Their Relationship with Soil Physico-Chemical Properties
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Tillage Systems and Cropping Patterns on Soil Physical Properties in Mozambique

1
Mozambique Agricultural Research Institute (IIAM), Maputo 3658, Mozambique
2
Department of Soil, Crop, and Climate Sciences, University of the Free State, Bloemfontein 9300, South Africa
*
Author to whom correspondence should be addressed.
Submission received: 26 August 2020 / Revised: 23 September 2020 / Accepted: 27 September 2020 / Published: 30 September 2020

Abstract

:
Conservation agriculture (CA) practices are advocated to reduce soil degradation, resulting in more sustainable food production as compared to conventional tillage (CT). In this study, the short-term effects of two tillage systems in combination with cropping patterns on selected soil physical parameters on four experimental sites in Mozambique were studied. The study sites differ according to their climatic conditions, soil types, and crop adaptation. Tillage systems evaluated were CA and CT, while the cropping pattern had four levels of sole cropping and three levels of intercropping. In general, soil physical properties showed significant changes due to the tillage systems, but the cropping pattern and their interaction with tillage systems did not yield significant impacts on the soil physical properties. CA increased bulk density, penetration resistance, and saturated hydraulic conductivity as compared to CT. A significant difference due to the tillage system was observed across the four sites, and in general, evaporation was higher in CT compared to CA. The presence of crop residues in CA contributed to lower evaporation. Thus, in the short term, CA practices could be a sustainable option to conserve soil water through higher infiltration and less evaporation.

1. Introduction

In sub-Saharan Africa (SSA), unsustainable farming practices, such as excessive soil tillage with no soil cover, have led to a decrease in soil fertility, accelerated erosion, and degradation of arable land. In Mozambique, around 68% of the population live in rural areas, and up to 90% are smallholder farmers [1]. The most common tillage system practiced by smallholder farmers in Mozambique is where hand hoes are used to till the soil, with crop residue collection and burning during land preparation. This is referred to as conventional tillage (CT) throughout this study. The use of inappropriate tillage systems (i.e., CT) has led to the deterioration of soil quality through decay in organic matter, disruption of soil structure, and decrease in aggregate stability. This resulted in lower macroporosity, thus promoting soil erosion, organic matter loss, and fertility depletion [2,3,4,5]. The development of sustainable management practices to reduce the physical degradation of soil should, therefore, be a high priority. Conservation agriculture (CA) practices are proposed to minimize soil degradation and improve soil quality. CA is defined as a tillage system where at least 30% of crop residues are left on the field surface, and it is a conservation practice to reduce soil erosion, improve water use efficiency, and regulate upper soil temperature [6,7]. CA integrates three principles, namely minimal soil disturbance, crop residue retention, and crop diversification through rotation or intercropping [6,8,9]. These CA practices are advocated to increase the microporosity in relation to macroporosity, consequently restoring soil aggregate stability [5,8,10]. Fortunately, to reduce soil degradation there is increased adoption of CA practices by smallholder farmers, as this is an attractive solution.
Crop production management practices cause changes in the physical attributes of soil, mainly due to changes in soil structure that disrupt the distribution and continuity of pores. Bulk density (Db), penetration resistance (PR), saturated hydraulic conductivity (Ks), available soil water, evaporation, and aggregate stability are the most vulnerable parameters used to assess the impact of tillage systems on soil physical quality (e.g., Zarnoza et al. [7], Busari and Salako [11]). The impact of different tillage operations on the physical properties is, however, controversial. The aim of CT is generally to loosen soil, resulting in higher macroporosity, and thus lowered bulk density and penetration resistance even over the short term [12]. Intensive use of CT practices can, however, promote subsurface compaction through the reduction of pore sizes, which coincides with an increase in Db and decrease in Ks. On the contrary, short-term impacts of CA can also lead to higher Db and PR values and a reduction in Ks under CA when compared to CT systems, mainly due to the absence or reduction of soil disturbance [2,13,14]. In the long term, the accumulation of organic matter under CA practices is expected to decrease the Db and PR and increase Ks, all of which are associated with increased macroaggregation in the top soil [14]. Furthermore, the presence of crop residues in CA practices will promote surface infiltration rates, reduce the impact of raindrops (sealing), and reduce evaporation [14,15].
The consortia between cereals and legumes under CA practices allow for higher organic matter accumulation. Biomass productions from cereal and legume associations are a critical source of organic carbon input that improves physical attributes [16] and decreases soil degradation. Organic matter maintains aggregate stability and improves soil physical quality by improving water storage through higher infiltration and lower runoff and evaporation. Thus, soil loss by water erosion is minimized, while the lower bulk density and penetration resistance promote root growth and development [3,17,18].
Although the literature suggests that physical attributes (e.g., Db, PR, and Ks) are improved under CA practices over the long term, it is often the short term impacts and benefits that promote adoption. These impacts must, therefore, be understood and optimized. In Mozambique, research on the individual effects of tillage systems and cropping patterns on soil physical properties has been conducted [19,20,21], but the information is scanty on the combined effects of tillage systems and cropping pattern on soil physical properties. Therefore, this study is conducted at four sites over two growing seasons to evaluate the impact of tillage systems and cropping patterns and their interactions on selected soil physical properties. This study hypothesizes that the tillage system and cropping pattern have significant impacts on selected soil physical properties, such as bulk density, penetration resistance, hydraulic conductivity, and evaporation in the short term.

2. Materials and Methods

2.1. Site Description

This study was conducted at four research stations of the Agricultural Research Institute of Mozambique (IIAM). These are the Nhacoongo station in the Inhambane province in the southern region, Gurúè station in the Zambézia province in the central region, Mutuali station in the Nampula province in the northern region, and Lichinga station in the Niassa province in the northern region of Mozambique (Figure 1). Some characteristics of the research stations are given in Table 1.

2.2. Experimental Design and Layout

The trials were conducted over two cropping seasons, i.e., 2016–2017 to 2017–2018, to evaluate the response on soil physical properties (Table 2). At each location, a field experiment was set up as a two factor split-plot experiment in a randomized complete block design and replicated four times. The first order factor was tillage systems, i.e., two levels—conventional tillage (CT) and conservation agriculture (CA), with 26 m × 37.5 m (975 m, including 2 m buffer zones between blocks and 1 m buffer zones between plots) main plots, and the second order factor was cropping patterns (seven levels; see Table 2) with 5.0 m × 4.5 m (22.5 m2, excluding any buffer zones) subplots (see Figure 2). The four sites had the same experimental design, and the treatments were applied to the same plots during the two years of study. For the CA system, glyphosate was applied at a rate of 2.5 L ha−1 (480 g L−1 of isopropylamine salt active ingredient) two weeks before planting. Under the CA system, the soil was disturbed only with a dibble stick planter when the seeds were planted, while under the CT system, soil was disturbed with a hand hoe to a depth of approximately 15 cm for seedbed preparation and planting. For the CT system, crop residues were removed, and for the CA system, the surface area was covered with at least 30% of the previous season harvested crop residues (2–3 t ha−1). No fertilizers were applied, since this is not common practice amongst the smallholders.
The sowing and harvesting dates were respectively in November and May every year for all crops, except for pigeon pea, due to a longer growing period, which was harvested in July of each year. Maize and pigeon pea seed were sown in rows at the recommended spacing of 75 cm between rows and 30 cm within rows, while the other legumes varied from 50 cm × 30 cm for common beans, cowpea, and groundnut, and 50 cm × 10 cm for soybean. The arrangement of row intercropping was 1:1 (one row of maize alternated by one row of legume). Maize and legumes were sown simultaneously as a sole crop or intercropped in alternative rows. Two seeds were sown per hill and later thinned to one per hill; gapping was done five days after emergence. Postemergence weeds were controlled with hand hoeing.

2.3. Measurement of Soil Properties

Soil physical properties were measured at the end of the second growing season (July 2018). Undisturbed core samples were collected from randomly located points in all plots for bulk density determination. Disturbed samples were collected at the same locations to measure particle size distribution. Saturated hydraulic conductivity, penetration resistance, and evaporation were determined in situ using methods discussed below.

2.3.1. Particle Size Distribution and Gravel

Disturbed composite soil samples of 50 g were from topsoil horizon of selected plots, oven-dried at 105 °C, and passed through a 2 mm sieve. The hydrometer method was used for the determination of particle size distribution [27]: sand (2.0–0.05 mm); silt (0.05–0.002 mm), and clay (<0.002 mm). This method is based on the contrasting settling velocities of the particle sizes within a water column. The textural class was deducted using the United States Department of Agriculture (USDA) textural triangle [27]. The gravel retained on the 2 mm sieve was used to estimate the percentage gravel in a composite soil sample.

2.3.2. Bulk Density

Undisturbed core samples were collected at the end of the second growing season (July 2018) from four randomly located points in all subplots. Cores with a 2.5 cm radius and 5 cm in length were used for sampling at Lichinga, Mutuali, and Gurúè, while slightly larger cores (3.5 cm radius and 12 cm length) were used in Nhacoongo. Bulk density was determined using the undisturbed core sampling method of Okalebo et al. [28]. Bulk density, Db (g·cm−3), was calculated as the relationship between dry soil mass and total ring volume using the equation below:
b   ( g · c m 3 ) = M V
where:
Db = Bulk density (g·cm−3);
M = Dry weight of the core sample (g);
V = Volume of the soil core (cm3).

2.3.3. Penetration Resistance

At the end of the second cropping season (July 2018), penetration resistance of the soil of each plot at all four sites was measured using a manual cone penetrometer (hand penetrometer Eijkelkamp 06.01.14). This instrument measures the soil resistance to penetration in Newton (N). Measurements were taken at 10 randomly selected points per plot within each plot to a depth of 15 cm, and the average value was calculated. The penetration resistance measurements were made by the same person in all four sites to ensure that the pressures were applied consistently on the penetrometer from one place to another using the equation below.
P R ( k P a ) =   F A
where:
PR = penetration resistance (N·cm−2 then converted to kPa);
F = normal force (N);
A = base area of the cone (cm2).

2.3.4. Saturated Hydraulic Conductivity

Saturated hydraulic conductivity was determined using the double-ring falling head infiltrometer method. The double ring infiltrometer method consisted of a 10 cm inner ring diameter and a 30 cm outer ring diameter inserted approximately 10 cm into the soil at two randomly selected points in selected plots at the experimental sites. Both the outer and inner rings were then saturated with water. The water level in the outer ring was kept above that of the inner ring to avoid lateral drainage from the inner ring. The time for water to infiltrate 2 cm was then recorded (from 10 to 8 cm above the soil surface). The infiltration test was repeated until the infiltration time became nearly constant (or less than a 5% decline from the previous infiltration time). The saturated hydraulic conductivity (Ks) was then calculated using the modified Bouwer and Rice [29] equation below.
K s = L t × ln ( h 0 + L ) ( h 1 + L )
where:
Ks = saturated hydraulic conductivity;
L = thickness of horizon (L);
t = time until constant infiltration rate was obtained (T);
h 0 and h 1 = head of water above the surface before and at the start of the test and after t , respectively.

2.3.5. Evaporation

Evaporation (E) was measured using microlysimeters (MLs) following the second cropping season (August, September, October, and November of 2018) according to the method proposed by Flumignam et al. [30]. The MLs (Figure 3B) consist of rigid polyvinyl chloride (PVC) pipes that were 15 cm long with 10 cm inner diameter. The wall thickness of the pipes was 0.25 cm. One of the sides was sharpened to ease penetration into the soil. A cap was inserted at the bottom to prevent deep drainage (Figure 3A). A slightly larger envelope was then inserted into the soil used for hosting the MLs (Figure 3C). The bottom of the MLs was sealed and raised 0.5 cm to eliminate the effects of drainage and runoff.
The MLs were saturated with water (approximately 100 mL), and the wet soil was weighed on a 2 g resolution electronic scale. Water evaporation was measured as the reduction in weight from the wet soil in the MLs on a daily basis for a week. This was repeated four times aligned to the site visits. The difference in weight between measurements was attributed to evaporation and was calculated using the equation below.
E   ( m m ) =   Δ M M L s A M L s
where:
E = water evaporation from the MLs (mm);
ΔMMLs = MLs variation in mass (kg);
AMLs = MLs surface area with a fixed value of 0.0785 m2.

2.4. Data Analysis

Data collected were tested for normality with the Shapiro–Wilk test and for homogeneity with the Levene test before analysis of variance (ANOVA) with XLSTAT v 17.3 software for Excel [31]. Treatment means were separated by Tukey’s multicomparison test, and differences were reported at the 5% probability level. The interaction between the tillage system and cropping pattern was given separately for those attributes that had a significant influence, and these were described accordingly.

3. Results

3.1. Particle Size Distribution and Gravel Content

The particle size distribution and resulting textural classes of the soil samples from the Nhacoongo, Mutuali, Lichinga, and Gurúè experimental sites are presented in Table 3. Soil texture is unlikely to be influenced by the treatments. It is, however, important to consider the texture when the impact on other physical properties are assessed. It is also important to consider the texture if recommendations and extrapolations are made to other sites. Soils from Nhacoongo are loamy sand, Mutuali and Gurúè are sandy loam, and Lichinga is sandy clay loam. In general, the highest gravel content was observed in Gurúè (5.00%), followed by Lichinga (3.28%), Mutuali (1.20%), and Nhacoongo (0.50%).

3.2. Bulk Density

The Db was influenced significantly by the tillage system at Nhacoongo, Mutuali, and Lichinga but not at Gurúè (Table 4). However, at all four sites, CA was associated with increases in Db of 1.29% for Nhacoongo, 2.63% for Mutuali, 7.50% for Lichinga, and 1.92% for Gurúè in comparison to under CT. On the other hand, the cropping pattern had no significant influence on Db at all four sites varying from 1.03 to 1.57 g·cm−3. The interaction between the tillage system and cropping pattern significantly influenced Db only in Nhacoongo and Lichinga but not in Mutuali and Gurúè.

3.3. Penetration Resistance

The effect of tillage systems on PR was significant only at Nhacoongo and Gurúè (Table 4). However, PR values across the four study sites were consistently higher under the CA system (Nhacoongo = 90.4%, Mutuali = 4.5%, Lichinga = 3.1%, and Gurúè = 12.2%) compared to the CT system. Cropping patterns, however, had a significant influence on PR only in Mutuali (Figure 4). The interaction effects between tillage systems and the cropping patterns had no significant influence on PR in Mutuali, Lichinga, and Gurúè. Overall, the interaction produced a higher PR under CA than under CT.

3.4. Saturated Hydraulic Conductivity

The Ks was significantly affected (p ≤ 0.05) by the tillage system in Nhacoongo, Lichinga, and Gurúè but not in Mutuali (Table 4). Generally, average Ks values were higher under the CA system (Nhacoongo = 14.7%, Mutuali = 13.9%, Lichinga = 65.5% and Gurúè = 74.2%) than under the CT system. The current study did not show any significant cropping pattern on Ks across the four sites ranging from 335 to 773 mm h−1. This was also the case with the interaction of tillage systems by cropping patterns.

3.5. Evaporation

A significant effect of tillage system on evaporation was obtained in some months at Nhacoongo, Mutuali, Lichinga, and Gurúè (Figure 5a–d). Significant differences in evaporation were at Nhacoongo and Gurúè, recorded in September, October, and November; at Lichinga in August, September, and November; and at Mutuali only in November. In most cases, evaporation was higher under CT than CA. Across the four sites, lower evaporation was recorded in the sandy loam soils of Gurúè (CA = 3.42 and CT = 3.61 mm) and Mutuali (CA = 6.22 and CT = 6.57 mm) than in the loamy sand soil of Nhacoongo (CA = 8.60 and CT = 11.27 mm) and the sandy clay loam soil of Lichinga (CA = 11.10 and CT = 10.98 mm). Unexpectedly, in August, higher evaporation values were obtained under CA than CT in Nhacoongo, Mutuali, and Gurúè.

4. Discussion

4.1. Bulk Density

Measurements made in the four experimental sites showed significant differences in Db as a result of tillage system. The CA system increased Db, resulting in increased compaction and less aeration compared to the CT system. Several studies (Fabrizzi et al. [2], Rashidi and Keshavarzpour [12], Malecka et al. [13], and Dembele [32]) reported greater Db values under CA than CT. Management practices adopted for CT (ploughing) tend to disturb the soil to a 30 cm depth, creating more pore space. Furthermore, the absence of tillage (CA system) increases Db because the soil becomes more compacted after only two cropping seasons. In the long term, decomposition of crop residues and organic matter accumulation might result in a decrease in Db [33,34]. It is clear that the expected long-term benefit associated with the higher organic matter has not yet been obtained. Across the four sites, the highest Db was recorded in the loamy sand soil of Nhacoongo (CA = 1.57 and CT = 1.55 g·cm−3), followed by sandy clay loam soil of Lichinga (CA = 1.29 and CT = 1.20 g·cm−3), sandy loam soil of Mutuali (CA = 1.17 and CT=1.14 g·cm−3), and the sandy loam soil of Gurúè (CA = 1.06 and CT = 1.04 g·cm−3). This was not unexpected and can be attributed to the fact that soil texture (loamy sand with >80% sand) differed from that in Mutuali and Gurúè (sandy loam with less than 75% sand). Sandy loam soils tend to have higher Db than more clayey variants [35]. At Lichinga, the soil was sandy clay loam with an average of 45% clay. Lower clay contents contribute to higher compaction, which explains the high Db in Nhacoongo (clay content < 9%) compared to other sites (clay content > 14%). Increased Db in sandy soil as compared to clayey soils was also recorded [36]. This can be ascribed firstly to the relatively high clay contents, which are not prone to compaction, and secondly to the short study period.
Cropping pattern (CP) had no significant influence on Db at all four sites. It was expected that plots with legumes and with more crop residues accumulation would have a positive effect on soil physical quality by improving the structure and decreased Db. However, no clear trends were observed between different cropping patterns. Similarly, Oliveira et al. [37] found no significant difference between beans and millet intercropping effects on Db.
A significant interaction effect between the tillage systems and CP on Db was recorded in Nhacoongo and Lichinga but not in Mutuali and Gurúè. The significant differences could be attributed to single outliers, and no clear trend was observed. Several other researchers have reported similar contradictory results. Nabanita [38] and Alam et al. [39] found insignificant interaction effects on Db between tillage systems and cropping patterns, whereas Rusu et al. [40] obtained significant interaction effects. The overall trend of tillage systems and cropping pattern interaction effects on Db under CA was higher than that under CT.

4.2. Penetration Resistance

PR values across the four study sites were consistently higher under CA systems compared to the CT system. Therefore, the CA system offered more resistance to root development than CT. This concurs with the results of Rashidi and Keshavarzpour [12] and Çelik [41]. The higher PR under CA can be due to almost no soil disturbance, which causes compaction and a decrease of macropore volume on the surface. Ploughing loosens the soil, leading to a lower PR. Across the four sites, the highest PR was recorded in Gurúè (CA = 2692 and CT = 2399 kPa) followed by Mutuali (CA = 1989 and CT = 1903 kPa) and Lichinga (CA = 1861 and CT = 1805 kPa), and unexpectedly, the loamy sand soil of Nhacoongo (CA = 1055 and CT = 554 kPa), despite that it recorded higher Db (CA = 1.57 and CT = 1.55 g·cm−3) than Mutuali (CA = 1.17 and CT = 1.14 g·cm−3), Lichinga (CA = 1.29 and CT = 1.20 g·cm−3) and Gurúè (CA = 1.06 and CT = 1.04 g·cm−3). Sandy soils (Gurúè and Mutuali) are more prone to compaction and surface crusting than clayey soils (Lichinga), thus explaining the difference in results obtained at the various sites.
The cropping pattern influenced PR significantly in Mutuali only. The interaction effects between tillage systems and CP had no significant influence on PR in Mutuali, Lichinga, and Gurúè. Across the four sites, higher PR values were observed under the CA system with CP interaction than under the CT system with CP interaction. According to some literature, the interaction between tillage systems and CP significantly influences PR [42]. At the same time, other researchers reported insignificant interaction effects, especially with short term studies [36].

4.3. Saturated Hydraulic Conductivity

In general, the hydraulic conductivity on all the sites was high (>300 mm h−1), which exceeds typical rainfall intensity. The Ks values were higher under CA systems than CT. Across the four sites, higher Ks values were recorded in the loamy sand soil of Nhacoongo (CA = 743 and CT = 648 mm h−1), while the lowest Ks values were recorded in the sandy loam soil of Mutuali (CA = 517 and CT = 454 mm h−1). Singh et al. [43] emphasized that the lower Ks values recorded under CT could be ascribed to the destruction of aggregates and reduction of macroporosity in ploughed soils. The positive impact of CA on Ks was somewhat surprising, as the general theory is that an increase in compaction and bulk density is typically associated with lower conductivity. The accumulation of organic matter associated with reduced physical aeration and accumulation of crop residues could be attributed to the increased conductivity. Organic matter and continuous root channels will increase the macroporosity under the CA, thereby increasing Ks. These results are consistent with those reported by Mcgarry et al. [44] and Kumar et al. [45], who found an increase in Ks under CA, probably due to the retention of crop residues that modified soil structure, hence aggregate stability. Although our study sites had high inherent Ks, the short term impact of CA in increasing Ks can be beneficial in other areas with low infiltration rates.

4.4. Evaporation

The tillage system had a significant influence on evaporation in some months at Nhacoongo, Mutuali, Lichinga, and Gurúè: Significant differences in evaporation were recorded at Nhacoongo and Gurúè in September, October, and November; at Lichinga in August, September, and November; and at Mutuali only in November. In most instances, evaporation was higher under the CT system than CA. These results are in agreement with those of Dembele [32], Singh et al. [46], and Thierfelder et al. [47], who reported improved water use efficiency due to the application of CA instead of CT. The higher evaporation under the CT system reflects the disadvantages of ploughing, which increases soil aeration and decay of soil aggregates, leading to enhanced water loss by evaporation. Conservation agriculture systems are effective in reducing evaporation and increase water holding capacity due to the presence of crop residues, which act as a protective skin for the soil by intercepting a great part of the incident solar radiation, creating thus a wetter microclimate at the soil surface, which enhances crop water availability. These findings are in agreement with those of Hobbs et al. [8], Alam et al. [39], and van Donk et al. [48], who reported that soils dry slowly on the surface due to the presence of crop residues under the CA system. The reduced evaporation is ascribed to shading, which moderates surface soil temperature and wind effects. Crop residues in CA systems also prevent the impact of direct rainwater on the soil, thus avoiding the formation of surface sealing, in turn maintaining water conductivity. Crop residues can reduce evaporation by 10–20% [49]. However, contradictory results were noted in Nhacoongo and Gurúè in August, Mutuali in August and October, and Lichinga in August, with higher evaporation recorded under CA than CT. The inconsistent trends recorded in August can be attributed to the fact that crop residues delay heat loss from the surface during winter, and this could promote warming during summer, leading to increased evaporation [50,51]. The inconsistency in the results could also be attributed to the use of microlysimeters, since they might not contain sufficient soil water to maintain realistic evaporation.

5. Conclusions

This study revealed some significant differences related to the tillage system across the four experimental sites in Mozambique after two years of employing different cultivation practices. Higher values of Db, PR, and Ks were recorded under the CA system when compared to the CT system. A change in soil physical properties usually requires a significant accumulation of organic matter, which requires a long time to achieve. The increased compaction, as indicated by Db and PR, is ascribed to the absence of tillage of the soil, whereas the increase in conductivity related to increased macroporosity is associated with the accumulation of organic matter through reduced tillage and crop residues. Evaporation at the soil surface is higher in the CT system than in the CA system. The decline in evaporation under CA, when compared to CT, is a short term benefit of CA that should be explored in the future to promote CA adoption under smallholder farmers. In the present study, tillage methods that significantly affected the soil physical properties due to decreased evaporation by reducing water loss were noted. Therefore, smallholder farmers should be encouraged to adopt CA practices to increase water availability in the soil during the crop growth period. Two years of experimentations were not enough to develop clear trends of sole cropping compared to intercropping. Thus, there is a need to conduct long-term studies to quantify cropping pattern effects on soil physical properties.

Author Contributions

Conceptualization, O.C.; J.v.T.; C.D.P. and G.C.; Field and lab analysis, O.C. and J.v.T.; writing—original draft preparation, O.C.; writing—review and editing, J.v.T.; C.D.P. and G.C.; project administration, O.C.; funding acquisition, O.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Agricultural Productivity Program for Southern Africa (APPSA) subproject financed by the World Bank. Any opinions, findings, conclusions, or recommendations expressed in this publication are those of the authors and do not necessarily reflect the view of the donors mentioned previously.

Acknowledgments

The authors acknowledged the staff of Nhacoongo, Mutuali, Lichinga, and Gurúè research stations for their assistance in the course of the experiments and data collection.

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. AGRA (Alliance for a Green Revolution in Africa). Africa Agriculture Status Report: Youth in Agriculture in Sub-Saharan Africa; Alliance for a Green Revolution in Africa (AGRA): Nairobi, Kenya, 2015. [Google Scholar]
  2. Fabrizzi, K.P.; Garcia, F.O.; Costa, J.L.; Picone, L.I. Soil water dynamics, physical properties and corn and wheat responses to a minimum and no-tillage system in the southern Pampas of Argentina. Soil Tillage Res. 2005, 81, 57–69. [Google Scholar] [CrossRef]
  3. Bolliger, A.; Magid, J.; Amado, T.J.C.; Neto, F.S.; Ribeiro, M.F.S.; Calegari, A.; Ralisch, R.; De Neergaard, A. Taking stock of the Brazilian “Zero-till revolution”: A review of landmark research and farmers’ practice. Adv. Agron. 2006, 91, 47–110. [Google Scholar]
  4. Wang, X.B.; Cai, D.X.; Hoogmoed, W.B.; Oenema, O.; Perdok, U.D. Developments in conservation tillage in rainfed regions of North China. Soil Tillage Res. 2007, 93, 239–250. [Google Scholar] [CrossRef]
  5. Alam, M.K.; Salahin, N.; Begum, R.; Hasanuzzaman, M.; Islam, M.; Rahman, M. Patterns of change in soil organic matter, physical properties and crop productivity under tillage practices and cropping systems in Bangladesh. J. Agric. Sci. 2017, 155, 216–238. [Google Scholar] [CrossRef]
  6. Mupangwa, W.; Twomlow, S.; Walker, S.; Hove, L. Effect of minimum tillage and mulching on maize (Zea mays L.) yield and water content of clay and sandy soils. Phys. Chem. Earth 2007, 32, 1127–1134. [Google Scholar] [CrossRef]
  7. Busari, M.A.; Salako, F.K. Soil hydraulic properties and maize root growth after application of poultry manure under different tillage systems in Abeokuta, Southern Nigeria. Arch. Agron. Soil Sci. 2015, 61, 223–237. [Google Scholar] [CrossRef]
  8. Hobbs, P.R.; Sayre, K.D.; Gupta, R. The role of conservation agriculture in sustainable agriculture. Philos. Trans. R. Soc. Ser. B Biol. Sci. 2007, 363, 543–555. [Google Scholar] [CrossRef]
  9. Thierfelder, C.; Matemba-Mutasa, R.; Runinamhodzi, L. Yield response of maize (Zea mays L.) to conservation agriculture cropping system in Southern Africa. Soil Tillage Res. 2014, 146, 230–242. [Google Scholar] [CrossRef]
  10. Bottinelli, N.; Angers, D.A.; Hallaire, V.; Michot, D.; Le Guillou, C.; Cluzeau, D.; Heddadj, D.; Menasseri-Aubry, S. Tillage and fertilization practices affect soil aggregate stability in a Humic Cambisol of Northwest France. Soil Tillage Res. 2017, 170, 14–17. [Google Scholar] [CrossRef]
  11. Zarnoza, R.; Acosta, J.A.; Bastida, F.; Dominguez, S.G.; Toledo, D.M.; Faz, A. Identification of sensitive indicators to access the interrelationship between soil quality, management practices, and human health. Soil 2015, 1, 173–185. [Google Scholar] [CrossRef] [Green Version]
  12. Rashidi, M.; Keshavarzpour, F. Effect of different tillage methods on soil physical properties and crop yield of melon (Cucumis melo). J. Agric. Biol. Sci. 2008, 3, 41–46. [Google Scholar]
  13. Malecka, I.; Blecharczyk, A.; Sawinska, Z.; Dobrzeniecki, T. The effect of various long term tillage systems on soil properties and spring barley yield. Turkish J. Agric. For. 2012, 36, 217–226. [Google Scholar]
  14. Singh, R.; Serawat, M.; Singh, A. Effect of tillage and crop residue management on soil physical properties. J. Soil Salin. Water Qual. 2018, 10, 200–206. [Google Scholar]
  15. Shukla, M.K.; Lal, R.B.; Owens, L.; Unkefer, P. Land use management impacts on structure and infiltration characteristics of soil in the north Appalachian region of Ohio. Soil Sci. 2003, 168, 167–177. [Google Scholar] [CrossRef]
  16. Jat, M.L.; Gathala, M.K.; Saharawat, Y.S.; Tetarwal, J.P.; Gupta, R.; Yadvinder, S. Double no till and permanent raised beds in maize-wheat rotation of north-western Indo-gangetic plains of India: Effects on crop yields water productivity, profitability and soil physical properties. Field Crops Res. 2013, 149, 291–299. [Google Scholar] [CrossRef]
  17. Strudley, M.W.; Green, T.R.; Ascough II, J.C. Tillage effects on soil hydraulic properties in space and time: State of the science. Soil Tillage Res. 2008, 99, 4–48. [Google Scholar] [CrossRef]
  18. Thierfelder, C.; Wall, P.C. Effects of conservation agriculture techniques on infiltration and soil water content in Zambia and Zimbabwe. Soil Tillage Res. 2009, 105, 217–227. [Google Scholar] [CrossRef]
  19. Tsujimoto, Y.; Pedro, J.A.; Boina, G.; Murracama, M.V.; Ito, O.; Tobita, S.; Oya, T.; Cuambe, C.E.; Martinho, C. Performance of maize-soybean intercropping under various N application rates and soil moisture conditions in Northern Mozambique. Plant Prod. Sci. 2015, 18, 365–376. [Google Scholar] [CrossRef]
  20. Chilundo, M. Effects of Irrigation and Fertilizer Management on Water and Nitrogen Use Efficiency in Maize on a Semi-Arid Loamy Sandy Soil. Ph.D. Thesis, Swedish University of Agricultural Sciences, Uppsala, Swedish, 2017; p. 84. [Google Scholar]
  21. Nascente, A.S.; Dambiro, J.; Cosntantico, C. Effects of grain-producing cover crops on rice grain yield in Cabo Delgado. Rev. Ceres 2017, 64, 607–615. [Google Scholar] [CrossRef] [Green Version]
  22. MAE (Ministério da Administração Estatal). Perfil do Distrito de Inharrime, Província de Inhambane; Ministério da Administração Estatal: Maputo, Mozambique, 2005; pp. 1–53. [Google Scholar]
  23. Yamada, R.; Koide, J.; Oya, T. Farm Management in Northern Mozambique; Yokendo: Sendai, Japan, 2017; pp. 1–82. [Google Scholar]
  24. MAE (Ministério da Administração Estatal). Perfil do Distrito de Lichinga, Província de Niassa; Ministério da Administração Estatal: Maputo, Mozambique, 2005; pp. 1–54. [Google Scholar]
  25. MAE (Ministério da Administração Estatal). Perfil do Distrito de Gurúè, Província da Zambézia; Ministério da Administração Estatal: Maputo, Mozambique, 2005; pp. 1–59. [Google Scholar]
  26. FAO (Food and Agriculture Organization of the United Nations). FAO/UNESCO Soil map of the World; FAO: Rome, Italy, 2016; Available online: Http://www.fao.org/fileadmin/user_upload/soils/docs/Soil_map_FAOUNESCO/acrobat/Africa_VI.pdf (accessed on 29 September 2020).
  27. United States Department of Agriculture (USDA). Soil Mechanics Level I. Module 3—USDA Textural Soil Classification; Soil Conservation Service, United States Department of Agriculture: Washington, DC, USA, 1987; pp. 1–48. [Google Scholar]
  28. Okalebo, J.R.; Gathua, K.W.; Woomer, P.L. Laboratory Methods of Soil and Plant Analysis: A Working Manual; TSBF: Nairobi, Kenya, 2002; pp. 1–131. [Google Scholar]
  29. Bouwer, H.; Rice, R. A slug test for determining hydraulic conductivity of unconfined aquifers with completely or partially penetrating wells. Water Resour. Res. 1976, 12, 423–428. [Google Scholar] [CrossRef] [Green Version]
  30. Flumignam, D.L.; Faria, R.T.; Lena, B.P. Test of a microlysimeter for measurement of soil evaporation. Eng. Agríc. 2012, 32, 80–90. [Google Scholar] [CrossRef] [Green Version]
  31. Addinsoft. XLSTAT 17.3 Statistical: Data Analysis and Statistic Software for Microsoft Excel; Addinsoft: Paris, France, 2019. [Google Scholar]
  32. Dembele, A. The Impact of Different Tillage and Soil Amendments on Soil Moisture Storage, Erosion and the Yield of Maize in the Moist Semi-Deciduous Forest Zone Of Ghana. Master’s Thesis, Faculty of Agriculture—Kwame Nkrumah University of Science and Technology, Kumasi, Ghana, 2015; pp. 1–125. [Google Scholar]
  33. Alam, M.K.; Salahin, N. Changes in soil physical properties and crop productivity as influenced by different tillage depths and cropping patterns. Bangladesh J. Agric. Res. 2013, 38, 289–299. [Google Scholar] [CrossRef]
  34. Barut, Z.B.; Celik, I. Tillage effects on some soil physical properties in a semi-arid Mediterranean region of Turkey. Chem. Eng. Trans. 2017, 58, 217–222. [Google Scholar]
  35. Landon, J.R. Booker Tropical Soil Manual. A Handbook for Soil Survey and Agricultural Land Evaluation in the Tropics and Subtropics; Longman: New York, NY, USA, 1991; p. 474. [Google Scholar]
  36. Malviya, S. Effect of Conservation Agricultural Practices on Selected Soil Physical Properties and Carbon Pools in Black Soils of Central India. Master’s Thesis, College of Agriculture—Department of Soil Science and Agricultural Chemistry, Jabalpur, India, 2014. [Google Scholar]
  37. Oliveira, F.C.C.; Ferreira, G.W.D.; Souza, J.L.S.; Vieira, M.E.O.; Pedrotti, A. Soil physical properties and soil organic carbon in northeast Brazil: Long-term tillage systems effects. Sci. Agric. Soil Plant Nut. 2020, 77, 1–6. [Google Scholar] [CrossRef] [Green Version]
  38. Nabanita, B. Impacts of Maize Based Conservation Agriculture Production System on Soil Properties in North Central Plateau Zone of Odisha. Master’s Thesis, College of Agriculture and Technology Bhubaneswar: Odisha, India, 2012. [Google Scholar]
  39. Alam, M.K.; Salahin, N.; Rashid, M.H.; Salam, M.A. Effects of different tillage practices and cropping patterns on soil physical properties and crop productivity. J. Trop. Resour. Sustain. Sci. 2013, 1, 51–61. [Google Scholar]
  40. Rusu, T.; Moraru, P.I.; Ranta, O.; Drocas, I.; Bogdan, I.; Pop, A.I.; Sopterean, M.L. No tillage and minimum tillage their impact on soil compaction, water dynamics, soil temperature and production on wheat, maize and soybean crop. Bull. UASVM Agric. 2011, 68, 318–323. [Google Scholar]
  41. Çelik, I. Effects of tillage methods on penetration resistance, bulk density and saturated hydraulic conductivity in a clayey soil conditions. J. Agric. Sci. 2011, 17, 143–156. [Google Scholar]
  42. Grant, C.A.; Lafond, G.P. The effects of tillage systems and crop sequences on soil bulk density and penetration resistance on a clay soil in southern Saskatchewan. Can. J. Soil Sci. 1993, 73, 223–232. [Google Scholar] [CrossRef]
  43. Singh, Y.; Bharadwaj, A.K.; Singh, S.P.; Singh, R.K.; Chaudhary, D.C.; Saxena, A.; Singh, V.; Singh, S.P.; Kuamr, A. Effect of rice (Oryza sativa) establishment methods, tillage practices in wheat (Triticum aestivum) and fertilization on soil physical properties and rice-wheat system productivity on silty clay Mollisol of Uttaranchal. Indian J. Agric. Sci. 2002, 72, 200–205. [Google Scholar]
  44. Mcgarry, D.; Bridge, B.J.; Radford, B.J. Contrasting soil physical properties after zero and traditional tillage of an alluvial soil in the semiarid subtropics. Soil Tillage Res. 2000, 53, 105–115. [Google Scholar] [CrossRef]
  45. Kumar, V.; Kumar, M.; Singh, S.K.; Jat, R.K.; Chandra, S.K. An appraisal of conservation agriculture on soil properties and crop—A review. Int. J. Pure Appl. Biosci. 2017, 5, 558–565. [Google Scholar] [CrossRef]
  46. Singh, B.; Humphreys, E.; Eberbach, P.L.; Katupitiya, A.; Singh, Y.; Kukal, S.S. Growth, yield and water productivity of zero till wheat as affected by rice straw mulch and irrigation schedule. Field Crops Res. 2011, 121, 209–225. [Google Scholar] [CrossRef]
  47. Thierfelder, C.; Chivenge, P.; Mupangwa, W.; Rosenstock, T.S.; Lamanna, C.; Eyre, J.X. How Climate-smart is conservation agriculture (CA)—It’s potential to deliver on adaptation, mitigation and productivity on smallholder farms in Southern Africa. Food Secur. 2017, 9, 537–560. [Google Scholar] [CrossRef] [Green Version]
  48. Van Donk, S.J.; Martin, D.L.; Irmak, S. Crop residues cover effects on evaporation, soil water, and yield of deficit-irrigated corn in West-central Nebraska. Trans. ASABE 2010, 53, 1787–1797. [Google Scholar] [CrossRef]
  49. Scopel, E.; Douzet, J.M.; Silva, F.A.M.; Cardoso, A.; Fideling, A.; Bernoux, M. Impactos do sistema de plantio direto com cobertura vegetal (SPDCV) sobre a dinâmica da água, do nitrogênio mineral e do carbono do solo Cerrado brasileiro. Cad. Ciência Tecnol. 2005, 22, 169–183. [Google Scholar]
  50. Horton, R.; Bristow, K.L.; Kluitenberg, G.J.; Sauer, T.J. Crop residue effects on surface radiation and energy balance—Review. Theor. Appl. Clim. 1996, 54, 27–37. [Google Scholar] [CrossRef] [Green Version]
  51. Sharratt, B. Freezing and Thawing- Cycles. In Encyclopedia of Soils in the Environment; Academic Press: New York, NY, USA, 2005. [Google Scholar]
Figure 1. Locations of the Lichanga (a), Gurúè (b), Mutuali (c), and Nhacoongo (d) research stations in Mozambique where the study was conducted.
Figure 1. Locations of the Lichanga (a), Gurúè (b), Mutuali (c), and Nhacoongo (d) research stations in Mozambique where the study was conducted.
Agriculture 10 00448 g001
Figure 2. Experimental layout for both conservation agriculture or conventional tillage treatments (buffers between plots and blocks are 1 and 2 m, respectively).
Figure 2. Experimental layout for both conservation agriculture or conventional tillage treatments (buffers between plots and blocks are 1 and 2 m, respectively).
Agriculture 10 00448 g002
Figure 3. Components of a microlysimeter system for measurement of evaporation: (A) Cap to avoid drainage, (B) microlysimeter, (C) outer envelope, and (D) installed microlysimeter. Source: Flumignam et al. [30].
Figure 3. Components of a microlysimeter system for measurement of evaporation: (A) Cap to avoid drainage, (B) microlysimeter, (C) outer envelope, and (D) installed microlysimeter. Source: Flumignam et al. [30].
Agriculture 10 00448 g003
Figure 4. Cropping pattern in Mutuali; Maize, sole maize, Ppea, sole pigeon pea; Sbean, sole soybean; Cpea, sole cowpea; M_Ppea, maize–pigeon pea intercropping; M_Cpea, maize–cowpea intercropping; M_Sbean, maize–soybean intercropping. Means followed by different letters within a column are significantly different using Tukey’s test at p ≤ 0.05.
Figure 4. Cropping pattern in Mutuali; Maize, sole maize, Ppea, sole pigeon pea; Sbean, sole soybean; Cpea, sole cowpea; M_Ppea, maize–pigeon pea intercropping; M_Cpea, maize–cowpea intercropping; M_Sbean, maize–soybean intercropping. Means followed by different letters within a column are significantly different using Tukey’s test at p ≤ 0.05.
Agriculture 10 00448 g004
Figure 5. Evaporation at Nhacoongo (a), Mutuali (b), Lichinga (c), and Gurúè (d) experimental sites. Bars represent means values and error bars one standard error; bars followed by different letters are significantly different using Tukey’s test at p ≤ 0.05.
Figure 5. Evaporation at Nhacoongo (a), Mutuali (b), Lichinga (c), and Gurúè (d) experimental sites. Bars represent means values and error bars one standard error; bars followed by different letters are significantly different using Tukey’s test at p ≤ 0.05.
Agriculture 10 00448 g005
Table 1. Selected characteristics of the Nhacoongo, Mutuali, Lichinga, and Gurúè research stations [22,23,24,25] where the experiments were conducted.
Table 1. Selected characteristics of the Nhacoongo, Mutuali, Lichinga, and Gurúè research stations [22,23,24,25] where the experiments were conducted.
Research StationAgroecological ZoneProvinceDistrictLocationAnnual Rainfall (mm)Temperature °CAltitude (M.A.S.L.)Soil Type *Textural Class **
NhacoongoR2InhambaneInharrime24°19′49.00″ S1000–120018–3368Luvic ArenosolsLoamy sand
35°12′55.00″ E
MutualiR7NampulaMalema14°52′14.02″ S300–130015–36.6574Orthic FerralsolsSandy loam
37°00′15.98″ E
LichingaR10NiassaLichinga13°18′46.01″ S1200–140016.1–32.91396Orthic FerralsolsSandy clay loam
35°14′26.02″ E
GurúèR10ZambéziaGurúè19°09′05.00″ S1800–200015–23678Humic NitisolsSandy loam
36°42′43.90″ E
* Based on FAO soil classification (FAO, [26]). ** Based on USDA textural soil classification (USDA, [27]).
Table 2. Details of tillage systems and cropping pattern at each of the four study sites.
Table 2. Details of tillage systems and cropping pattern at each of the four study sites.
Tillage System: Conservation Agriculture (CA)/Conventional Tillage (CT)
NhacoongoMutualiLichingaGurúè
Cropping pattern at each study siteSole maize (Matuba) *Sole maize (Matuba)Sole maize (Matuba)Sole maize (Matuba)
Sole groundnut (JL24)Sole pigeon pea (00554)Sole pigeon pea (00554)Sole cowpea (IT16)
Sole pigeon pea (00554)Sole cowpea (IT16)Sole soybean (TGX 1937-1F)Sole soybean (TGX 1937-1F)
Sole cowpea (IT16)Sole soybean (TGX 1937-1F)Sole common beans (Nua 45)Sole common beans (Nua 45)
Maize–groundnut intercroppingMaize–pigeon pea intercroppingMaize–pigeon pea intercroppingMaize–cowpea intercropping
Maize–pigeon pea intercroppingMaize–cowpea intercroppingMaize–soybean intercroppingMaize–soybean intercropping
Maize–cowpea intercropMaize–soybean intercropMaize–common beans intercroppingMaize–common beans intercropping
* Crops planted: maize (Zea mays L.), cowpea (Vigna unguiculata L.), groundnut (Arachis hypogaea L.), pigeon pea (Cajanus cajan L.), soybean (Glycine max L.), and common beans (Phaseolus vulgaris L.).
Table 3. Average particle size distribution and gravel content of the soil at the four experimental sites.
Table 3. Average particle size distribution and gravel content of the soil at the four experimental sites.
LocationClay (%)Silt (%)Sand (%)Textural Class *Gravel (%)
Nhacoongo8.85.385.9loamy sand0.5
Mutuali14.314.471.3sandy loam1.2
Lichinga31.513.554.7sandy clay loam3.3
Gurúè12.08.879.2sandy loam5.0
* Based on USDA soil textural classification (USDA, [26]).
Table 4. Effect of tillage systems on soil physical properties at Nhacoongo, Mutuali, Lichinga, and Gurúè experimental sites.
Table 4. Effect of tillage systems on soil physical properties at Nhacoongo, Mutuali, Lichinga, and Gurúè experimental sites.
Db (g·cm−3)PR (kPa)Ks (mm h−1)
Tillage SystemNhacoongoMutualiLichingaGurúèNhacoongoMutualiLichingaGurúèNhacoongoMutualiLichingaGurúè
CA1.57 a1.17 a1.29 a1.06 a1055 a1989 a1861 a2692 a743 a517 a619 a669 a
CT1.55 b1.14b1.20 b1.04 a554 b1903 a1805 a2399 b648 b454 a374 b384 b
L.S.D.0.010.030.040.036113722016289196137228
CPnsnsnsnsns*nsnsnsnsnsns
CP x TS*ns*nsnsnsnsnsnsnsnsns
CV (%)2.477.748.926.9114.0513.0822.2111.8223.7537.3544.6353.07
CA, conservation agriculture; CT, conventional tillage; CP, cropping pattern; TS, tillage system; Db, bulk density; PR, penetration resistance; Ks, saturated hydraulic conductivity; L.S.D., least significant difference at p ≤ 0.05; ns, not significantly different at p ≤ 0.05; *, significantly different at p ≤ 0.05. Means followed by different letters within a column are significantly different at p ≤ 0.05.

Share and Cite

MDPI and ACS Style

Chichongue, O.; van Tol, J.; Ceronio, G.; Du Preez, C. Effects of Tillage Systems and Cropping Patterns on Soil Physical Properties in Mozambique. Agriculture 2020, 10, 448. https://0-doi-org.brum.beds.ac.uk/10.3390/agriculture10100448

AMA Style

Chichongue O, van Tol J, Ceronio G, Du Preez C. Effects of Tillage Systems and Cropping Patterns on Soil Physical Properties in Mozambique. Agriculture. 2020; 10(10):448. https://0-doi-org.brum.beds.ac.uk/10.3390/agriculture10100448

Chicago/Turabian Style

Chichongue, Oscar, Johan van Tol, Gert Ceronio, and Chris Du Preez. 2020. "Effects of Tillage Systems and Cropping Patterns on Soil Physical Properties in Mozambique" Agriculture 10, no. 10: 448. https://0-doi-org.brum.beds.ac.uk/10.3390/agriculture10100448

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