Soil macropore dynamics affected by tillage and irrigation for a silty loam alluvial soil in southern Portugal

Soil tillage can have a significant effect on soil porosity and water infiltration. This study reports field measurements of near saturated hydraulic conductivity in an undisturbed soil under two tillage treatments, conventional tillage (CT) and minimum tillage (MT). The objective was to determine effective macro and mesoporosities, porosity dynamics during the irrigation season, and their contribution to water flow. Field observations were performed during the 1998 maize (Zea mays L.) cropping season in an Eutric Fluvisol with a silty loam texture, located in the Sorraia River Watershed in the south of Portugal. Infiltration measurements were done with a tension infiltrometer. At each location an infiltration sequence was performed corresponding to water tensions (φ) of 0, 3, 6 and 15 cm. Five sets of infiltration measurements were taken in both treatments in the top soil layer between May and September. One set of measurements was done at the depth of 30 cm at the bottom of the plowed layer in the CT plot. After 5 years of continuous tillage treatments the results show that regardless of the tillage treatment, saturated conductivity values K(φ₀) were several times larger than near saturation conductivity K(φ₃). This indicates that subsurface networks of water conducting soil pores can exist in both CT and MT maize production systems. In CT, the moldboard plow created macro and mesoporosity in the top soil layer while breaking pore continuity at 30 cm depth. This porosity was partially disrupted by the first irrigation, resulting in a significant decrease of 45% in the macropore contribution to flow. Later in the season, the irrigation effect was overlaid by the root development effect creating new channels or continuity between existing pores. In MT macroporosity contribution to flow did not show significant differences in time, representing 85% of the total flow. In both the treatments, macropores were the main contributing pores to the total flow, in spite of the very low macroporosity volumes.     

Dynamics of soil hydraulic properties during fallow as affected by tillage

There is limited information on the effects of tillage practices on soil hydraulic properties, especially changes with time. The objective of this study was to evaluate on a long-term field experiment the influence of conventional tillage (CT), reduced tillage (RT) and no-tillage (NT) on the dynamics of soil hydraulic properties over 3 consecutive 16–18 month fallow periods. Surface measurements of soil dry bulk density (ρ_b), soil hydraulic conductivity (K(ψ)) at −14, −4, −1 and 0 cm pressure heads using a tension disc infiltrometer, and derived hydraulic parameters (pore size, number of pores per unit of area and water-transmission porosity) calculated using the Poiseuille's Law were taken on four different dates over the fallow period, namely, before and immediately after primary tillage, after post-tillage rains and at the end of fallow. Under consolidated structured soil conditions, NT plots presented the most compacted topsoil layer when compared with CT and RT. Soil hydraulic conductivity under NT was, for the entire range of pressure head applied, significantly lower (P < 0.05) than that measured for CT and RT. However, NT showed the largest mean macropore size (0.99, 0.95 and 2.08 mm for CT, RT and NT, respectively; P < 0.05) but the significantly lowest number of water-conducting pores per unit area (74.1, 118.5 and 1.4 macropores per m² for CT, RT and NT, respectively; P < 0.05). Overall, water flow was mainly regulated by macropores even though they represented a small fraction of total soil porosity. No significant differences in hydraulic properties were found between CT and RT. In the short term, tillage operations significantly increased K (P < 0.05) for the entire range of pressure head applied, which was likely a result of an increase in water-conducting mesopores despite a decrease in estimated mesopore diameter. Soil reconsolidation following post-tillage rains reduced K at a rate that increased with the intensity of the rainfall events.     

Using the simplified falling head technique to detect temporal changes in field-saturated hydraulic conductivity at the surface of a sandy loam soil

Determining temporal changes in field-saturated hydraulic conductivity (K_fs) is important for understanding and modeling hydrological phenomena at the field scale. Little is known about temporal variability of K_fs values measured at permanent sampling points. In this investigation, the simplified falling head (SFH) technique was used for an approximately 2-year period to determine temporal changes in K_fs at 11 permanent sampling points established at the surface of a sandy loam soil. Additional K_fs measurements were obtained by the single-ring pressure infiltrometer (PI) technique to also compare the SFH and PI techniques. The lowest mean values of K_fs, M(K_fs), were detected in December and January (20.5 ≤ M(K_fs) ≤ 146.2 mm h⁻¹), whereas higher results (190.5 ≤ M(K_fs) ≤ 951.9 mm h⁻¹) were obtained in the other months of the year. The K_fs values were higher and less variable in the dry soil (θ_i ≤ 0.21 m³ m⁻³, M(K_fs) = 340.6 mm h⁻¹, CV(K_fs) = 106%) than in the wet one (θ_i > 0.21 m³ m⁻³, M(K_fs) = 78.4 mm h⁻¹, CV(K_fs) = 185%). Both wet and dry soil were less conductive at the end of the study period than at the beginning one but a more appreciable change was detected for the dry soil (K_fs decreasing by 83.4%) than for the wet one (K_fs decreasing by 63.0%). The simple SFH technique yielded K_fs results similar to the more laborious and time-consuming PI technique (i.e., mean values differing at the most by a factor of two). It was concluded that (i) the soil water content was an important factor affecting the K_fs results obtained in a relatively coarse-textured soil, (ii) the impact of time from the beginning of the experiment on the saturated hydraulic conductivity was larger for a repeated sampling of dry soil than of wet soil and (iii) the SFH technique yielded reliable K_fs results in a relatively short period of time without the need for extensive instrumentation or analytical methodology.     

Using Cs and Pbex for quantifying soil organic carbon redistribution affected by intensive tillage on steep slopes

Studying on spatial and temporal variation in soil organic carbon (SOC) is of great importance because of global environmental concerns. Tillage-induced soil erosion is one of the major processes affecting the redistribution of SOC in fields. However, few direct measurements have been made to investigate the dynamic process of SOC under intensive tillage in the field. Our objective was to test the potential of ¹³⁷Cs and ²¹⁰Pb_ex for directly assessing SOC redistribution on sloping land as affected by tillage. Fifty plowing operations were conducted over a 5-day period using a donkey-drawn moldboard plow on a steep backslope of the Chinese Loess Plateau. Profile variations of SOC, ¹³⁷Cs and ²¹⁰Pb_ex concentrations were measured in the upper, middle and lower positions of the control plot and the plot plowed 50 times. ¹³⁷Cs concentration did not show variations in the upper 0–30 cm of soil whereas ²¹⁰Pb_ex showed a linear decrease (P < 0.05) with soil depth in the upper and middle positions, and an exponential decrease (P < 0.01) at the lower position of the control plot. The amounts of SOC, ¹³⁷Cs and ²¹⁰Pb_ex of sampling soil profiles increased in the following order: lower > middle > upper positions on the control plot. Intensive tillage resulted in a decrease of SOC amounts by 35% in the upper and by 44% in the middle positions for the soil layers of 0–45 cm, and an increase by 21% in the complete soil profile (0–100 cm) at the lower position as compared with control plot. Coefficients of variation (CVs) of SOC in soil profile decreased by 18.2% in the upper, 12.8% in the middle, and 30.9% in the lower slope positions whereas CVs of ¹³⁷Cs and ²¹⁰Pb_ex decreased more than 31% for all slope positions after 50 tillage events. ¹³⁷Cs and ²¹⁰Pb_ex in soil profile were significantly linearly correlated with SOC with R² of 0.81 and 0.86 (P < 0.01) on the control plot, and with R² of 0.90 and 0.86 (P < 0.01) on the treatment plot. Our results evidenced that ³⁷Cs and ²¹⁰Pb_ex, and SOC moved on the sloping land by the same physical mechanism during tillage operations, indicating that fallout ¹³⁷Cs and ²¹⁰Pb_ex could be used directly for quantifying dynamic SOC redistribution as affected by tillage erosion.     

Effect of tillage and mode of straw mulch application on soil erosion in the submontaneous tract of Punjab, India

The submontaneous tract of Punjab comprising 10% of the state, is prone to soil erosion by water. Soils of the area are coarse in texture, low in organic matter and poor in fertility. High intensity rains during the monsoon season result in fertile topsoil removal. There is an urgent need to control soil erosion in this region so as to improve soil productivity. A field study was conducted to estimate the effect of tillage and different modes of mulch application on soil erosion losses. Treatments comprised two levels of tillage, viz. minimum (T_m) and conventional (T_c) in the main plots and five modes of straw mulch application, viz. mulch spread over whole plot (M_w), mulch spread on lower one-third of plot (M_1/3), mulch applied in strips (M_s), vertical mulching (M_v) and unmulched control (M_o), in subplots in a replicated split plot design. Rate of mulch application was 6 t ha⁻¹ in all modes. Compared with M_o, M_w reduced runoff by 33%. Runoff and soil loss were 5 and 40% higher under T_c than under T_m. Though other modes of straw mulch application (M_1/3, M_s and M_v) controlled soil loss better than M_o, their effectiveness was less than M_w. T_m was more effective in conserving soil moisture than T_c. Compared with M_o, M_w had 3–7% higher soil moisture content in the 0–30 cm soil depth under T_m. Minimum soil temperature of the surface layer was 1.4–2.4 °C lower under M_w than under M_o. Straw mulching reduced maximum soil temperature and helped in conserving soil moisture. Minimum tillage coupled with M_w was highly effective in reducing soil erosion losses, decreasing soil temperature and increasing moisture content by providing maximum surface cover.     

Comparison of concentrated-flow detachment equations for low shear stress

Several equations exist to describe the relationship between concentrated-flow detachment and shear stress (τ). However, the advantages and disadvantages of these equations for specific circumstances remain unclear. This study examines the performance of linear and power equations with and without a critical shear stress (τ_c) term for concentrated-flow detachment at low shear stress. Equations were fit to data collected from experiments on five midwestern US soils using flume experimental data at low shear stress levels. Field experimental data were also available for these soils. The linear equation was simple to use and parameter values were easily estimated with linear regression. However, significant lack of fit was found when the linear equation was applied to data collected from low to medium shear stress levels. The value of soil erodibility (K) for a soil varied by a factor of 3 and critical shear stress (τ_c) varied by a factor of 2.5. The linear equation prediction underestimated detachment (D) by 25% at high shear stress and overestimated detachment by 30% at low shear stress. In contrast, the power equations gave more stable erodibility parameters because these equations reduced the systematic nature of the observation residuals found with the linear equation. Correlation between rill detachment D and τ was generally lower with the linear compared to the power equations for conditions tested. For higher shear stresses and longer slopes, the linear equation may be acceptable where field experiments show a linear trend. It is suggested that τ_c only be used when it has a value significantly different from zero.     

Wind tunnel simulation of ridge-tillage effects on soil erosion from cropland

In the arid and semi-arid regions, ridge tillage was often used as an alternative practice for wind erosion control on the croplands without sufficient crop residues left during the fallow period. Through wind tunnel experiments, wind erosion rate and vertical mass flux profile of blown sand under the simulated conditions of ridge tillage and flat tillage were studied in 15, 10, 10, 5, 3 min exposures at the wind velocities of 8, 10, 15, 20, 24 m s⁻¹, respectively. The results for the soil tested indicate that the mean rate of wind erosion under flat tillage was 129.89 g m⁻² min⁻¹, while that under ridge tillage were 20–60% less. Under ridge tillage with different structures, average wind erosion rate had a positive correlation with the spacing between adjacent ridges. For the same ridge height, average wind erosion rate decreased with increasing ratio between the height of ridge and the width of furrow. For the same ratio between the height of ridge and the width of furrow, average wind erosion rate increased with increasing height of ridge. Power function relationships were found between wind erosion rate and wind velocity on all the simulated tillage conditions. A wind velocity of 15 m s⁻¹ was the critical velocity, above which wind erosion rate increased rapidly for the soil and simulated tillage conditions tested. Compared with flat tillage, ridge tillage remarkably decreased wind erosion rates when wind velocities were beyond 15 m s⁻¹. Under ridge tillage, the total mass of sand transported at a height of 0–20 cm above soil surface (Q_0–20), and the fraction of that travelling at a height of 0–4 cm (Q_0–4/Q_0–20), were less man mat under flat tillage. For the same ridge height, Q_0–4/Q_0–20 increased with increasing ratio between the height of ridge and the width of furrow. For the same ratio between the height of ridge and the width of furrow, Q_0–4/Q_0–20 decreased with increasing height of the ridge. Sand transport rate under flat tillage decreased with increasing height by a negative exponential function, while negative linear functions were found under ridge tillage. Thus ridge tillage decreased the rate of wind erosion and sand transportation near soil surface, reduced the loss of soil nutrient caused by wind erosion and plant damage caused by blown sand abrasion, which make it an effective agricultural technology for wind erosion control in the arid and semi-arid regions.     

Soil quality indicators under continuous cropping systems in the Argentinean Pampas

Over the last two decades, soil cultivation practices in the southern Argentinean Pampas have been changing from a 7 year cash-crop production system alternated with 2–3 years under pasture, to a continuous cropping system. A better understanding of the impact of the period of time a field has been under continuous cropping on a broad spectrum of soil properties related to soil quality is needed to target for sustainable cropping systems. The objectives of this study were to: (i) assess the relationship between physical and chemical soil parameters related to soil quality and (ii) identify soil quality indicators sensitive to soil changes under continuous cropping systems in the Argentinean Pampas.

Correlation analysis of the 29 soil attributes representing soil physical and chemical properties (independent variables) and years of continuous cropping (dependent variable) resulted in a significant correlation (p < 0.05) in 78 of the 420 soil attribute pairs. We detected a clear relationship between hydraulic conductivity at tension h (K_h) and structural porosity (ρ_e); ρ_e being a simple tool for monitoring soil hydraulic conditions.

Soil tillage practice (till or no-till) affected most of the soil parameters measured in our study. It was not possible to find only one indicator related to the years under continuous cropping regardless of the cultivation practice. We observed a significant relationship between years under continuous cropping and K_h under no-till (NT) and wheat fallow (p < 0.001, R² = 0.70). Under these conditions, K₋₄₀ diminished as the number of years under continuous cropping increased.

The change in mean weight diameter (CMWD) was the only physical parameter related to the number of years under continuous cropping, explaining 36% of the variability in the number of years under continuous cropping (p < 0.001) The combination of three soil quality indicators (CMWD, partial R² = 0.38; slope of the soil water retention curve at its inflexion point (S), partial R² = 0.14 and cation exchange capacity (CEC), partial R² = 0.13) was able to explain, in part, the years under continuous cropping (R² = 0.65; p value > 0.001), a measure related to soil quality.     

Estimating the saturated hydraulic conductivity in a spatially variable soil with different permeameters: a stochastic Kozeny–Carman relation

The spatial variability of the saturated hydraulic conductivity (K_s) of a greenhouse banana plantation volcanic soil was investigated with three different permeameters: (a) the Philip-Dunne field permeameter, an easy to implement and low cost device; (b) the Guelph field permeameter; (c) the constant head laboratory permeameter. K_s was measured on a 14×5 array of 2.5 m×5 m rectangles at 0.15 m depth using the above three methods. K_s differences obtained with the different permeameters are explained in terms of flow dimensionality and elementary volume explored by the three methods. A sinusoidal spatial variation of K_s was coincident with the underlying alignment of banana plants on the field. This was explained in terms of soil disturbances, such as soil compaction, originated by management practices and tillage. Soil salinity showed some coincidence in space with the hydraulic conductivity, because of the irrigation system distribution, but a causal relationship between the two is however difficult to support. To discard the possibility of an artefact, the original 70 point mesh was doubled by intercalation of a second 14×5 grid, such that the laboratory K_s was finally determined on a 140 points 2.5 m×2.5 m square grid. Far from diluting such anisotropy this was further strengthened after inclusion of the new 70 points. The porosity (φ) determined on the same laboratory cores shows a similar sigmoid trend, thus pointing towards a plausible explanation for such variability. A power-law relationship was found between saturated hydraulic conductivity and porosity, K_sφⁿ (r²=0.38), as stated by the Kozeny–Carman relation. A statistical reformulation of the Kozeny–Carman relation is proposed that both improved its predictability potential and allows comparisons between different representative volumes, or K_s data sets with different origin. Although the two-field methods: Guelph and Philip-Dunne, also follow a similar alignment trend, this is not so evident, suggesting that additional factors affect K_s measured in the field. Finally, geostatistical techniques such as cross correlograms estimation are used to further investigate this spatial dependence.     

Soil physical responses to novel rice cultural practices in the rice–wheat system: Comparative evidence from a swelling soil in Nepal

Soil puddling in advance of rice (Oryza sativa L.) transplanting disperses surface aggregates and generates compaction at depth. As a management scheme for rice, puddling is typically considered advantageous for maximizing resource availability and yield. However, some experimental findings suggest a conflict between edaphic conditions created by this establishment technique and the performance of subsequent non-rice crops like wheat (Triticum aestivum L.). At a site in the mid-hills region of Nepal on a silt loam soil with vertic characteristics, we compared the impact of six rice tillage (surface tillage—T₁, shank subsoiler—T₂, shank subsoiler + moldboard plough—T₃) and establishment (soil puddling + transplanting—TPR, direct seeding—DSR) combinations on soil physical properties over two cycles of the rice–wheat rotation. For the rice season, 0–20 cm saturated hydraulic conductivity (K_sat) in the DSR plots was 2.6 and 4.3 times higher than their TPR counterparts in the first (Y₁) and second (Y₂) years, respectively (TPR-Y₁ = 93 mm day⁻¹, DSR-Y₁ = 241 mm day⁻¹, TPR-Y₂ = 133 mm day⁻¹, DSR-Y₂ = 582 mm day⁻¹), whereas tillage method did not significantly influence K_sat in this soil layer. The impact of rice establishment method was reflected in higher TPR bulk densities in the 5–10 (DSR = 1.19 g cm⁻³, TPR = 1.24 g cm⁻³) and 10–15 cm (DSR = 1.24 g cm⁻³, TPR = 1.29 g cm⁻³) depth increments in the wet season. Although none of the treatments significantly influenced the position or thickness of the plough sole, penetration resistance profiles suggest that vertical fractures with reduced soil strength were created within the pan region by deep tillage (T₂ and T₃), although these features were not associated with higher hydraulic conductivities from 20 to 50 cm. As the soils dried at the end of the rice season, crack propagation in the deep tilled plots (T₂ and T₃) was more pervasive. During the wheat season, comparable bulk density profiles and soil moisture retention characteristics across the treatments suggest that many of the edaphic changes induced by contrasting rice tillage and establishment practices did not persist in the self-mulching, vertic soils at our site. Conversely, significant increases in K_sat among the DSR plots from Y₁ to Y₂ (Y₁ = 241 mm day⁻¹, Y₂ = 582 mm day⁻¹) imply a temporal element to soil structural regeneration with adoption of direct seeding.