Prediction of the soil structures produced by tillage

Data are presented for the amount of clods >50 mm produced when five different soils were tilled at a range of different, naturally occurring water contents. The optimum water content for soil tillage is defined as that at which the amount of clods produced is minimum. The amount clods produced at this optimum water content is shown to be linearly and negatively correlated with the value of Dexter's index S of soil physical quality. This results in a rational model for soil tillage that enables predictions to be made for all different soils and conditions. Pedo-transfer functions can be used to estimate the input parameters for the model for cases, for which measured values are not available. It is concluded that for soils with good physical condition (i.e. S > 0.035), no clods >50 mm are produced during tillage.     

Predicting the saturated hydraulic conductivity of compacted subsoils using the non-similar media concept

It is widely recognized that saturated hydraulic conductivity is dominated by the micromorphology of soil pores rather than by the merely total porosity or dry bulk density. Nevertheless, some researchers are reporting that the decrease in saturated hydraulic conductivity of subsoil is simply associated with the decrease in soil porosity or increase in dry bulk density. Based on these understandings in published papers and on our preliminary field investigation, we assumed that micromorphology of soil pores in topsoils is subjected to be destroyed with continuous disturbance by frequent tillage while subsoils tend to be compacted without serious changes of micromorphology of soil pores. Thus, we focused on finding the dependence of saturated hydraulic conductivity on dry bulk density by separating the soils into tilled layer and compacted layer. The objective of this study was to describe the relationship between saturated hydraulic conductivity and dry bulk density using a theoretical model, the non-similar media concept (NSMC) model, capable of predicting saturated hydraulic conductivities of soils with different values of dry bulk densities. The study area was located near the Tone River in Saitama Prefecture, Japan, where the soils were classified into Haplic Brown Lowland Soils according to the Classification of Cultivated Soils in Japan (Eutric Fluvisol according to FAO/UNESCO). Two sites, where the topsoils were seasonally tilled while the subsoils were sustained as it is, and another site where the topsoil was seasonally tilled, too, but extra deep tillage (1 m tillage depth) had been done, were chosen for the measurements. The saturated hydraulic conductivities and dry bulk densities of undisturbed soil cores from different depths were measured in the laboratory. The NSMC model was carefully applied only when the soil textures were the same among samples. The well-known conventional equations formulated by Kozeny–Carman and by Campbell, were used to compare the applicabilities with the NSMC model. The NSMC model succeeded in predicting the saturated hydraulic conductivities in the compacted subsoils. On the other hand, the NSMC model was not applicable to the tilled topsoils and to the deeply tilled subsoil. The saturated hydraulic conductivity of tilled topsoils and deeply tilled subsoil was always lower than that of compacted subsoils at the same dry bulk densities. The Kozeny–Carman and Campbell equations both failed in the prediction of saturated hydraulic conductivity in subsoil. It was concluded that the saturated hydraulic conductivity of subsoils under compaction without extreme disturbance is well related with its dry bulk density by the NSMC model.     

Quantitative analysis of earthworm burrow systems with respect to biological soil-structure regeneration after soil compaction

 On arable land, tilled with conventional tillage (CT) and conservation tillage (CS) respectively, plots were compacted by wheeling them 6 times with a 5 Mg wheel load in spring 1995. Immediately after compaction, undisturbed soil monoliths were excavated from the compacted and uncompacted plots. The monoliths were defaunated and inoculated with either Lumbricus terrestris or Aporrectodea caliginosa. One monolith from each plot remained uninoculated as a control. After 6 months the monoliths were defaunated again and then scanned with X-ray helical computed tomography. The data were transformed, the void systems inside the monoliths were reconstructed and visualised, and the parameters total void length, total void volume, tortuosity and continuity were quantified. The parameters' values were generally lower in the controls than in the inoculated monoliths. Differences in burrow construction could be explained by the different life strategies of the two earthworm species. Changes in burrow morphology due to tillage system and soil compaction were minor. Only the continuity of the burrow systems clearly changed: decreasing for L. terrestris and increasing for A. caliginosa. This can be explained by a change in the earthworms' burrowing activity to minimise energy expenditure in compacted soil. By extrapolating field data, we concluded that earthworms have great potential for biologically regenerating the soil structure after a single compaction event. Due to higher earthworm abundances in soil managed by CS the regeneration of the soil structure is assumed to be better in these plots than those tilled by CT. Received: 17 December 1997     

Cotton management in a compacted subsurface microirrigated coastal plain soil of the southeastern US

A loamy sand Acrisol (Aquic Hapludult) that had been microirrigated for 6 years became so severely compacted that it had root limiting values of soil cone index in the Ap horizon and a genetic hardpan below it. Deep and surface tillage systems were evaluated for their ability to alleviate compaction. Deep tillage included subsoiling or none. Both deep tillage treatments were also surface tilled by disking, chiseling, or not tilling. Subsoiling was located in row or between rows to avoid microirrigation tubes (laterals) that were buried under every other mid row or every row. Cotton (Gossypium hirsutum) was planted in 0.96-m wide rows. Cotton yield was improved by irrigation from 485 to 1022 kg ha⁻¹ because both 2001 and 2002 were dry years. Tillage loosened the soil by an average of 0.5–1.3 MPa; but compacted zones remained outside tilled areas. Subsoiling improved yield by 131 kg ha⁻¹ when performed in row where laterals were placed in the mid rows; but subsoiling did not improve yield when it was performed in mid rows. For subsurface irrigation management in these soils, the treatment with laterals buried under every other mid row was able to accommodate in-row subsoiling which improved yield; and this treatment was just as productive as and had been shown to be less expensive to install than burying laterals under every row.     

A simplified modelling approach for quantifying tillage effects on soil carbon stocks

Soil tillage has been shown to affect long‐term changes in soil organic carbon (SOC) content in a number of field experiments. This paper presents a simplified approach for including effects of tillage in models of soil C turnover in the tilled‐soil layer. We used an existing soil organic matter (SOM) model (CN‐SIM) with standard SOC data for a homogeneous tilled layer from four long‐term field experiments with conventionally tilled (CT) and no‐till (NT) treatments. The SOM model was tested on data from long‐term (>10 years) field trials differing in climatic conditions, soil properties, residue management and crop rotations in Australia, Brazil, the USA and Switzerland. The C input for the treatments was estimated using data on crop rotation and residue management. The SOM model was applied for both CT and NT trials without recalibration, but incorporated a ‘tillage factor’ (TF) to scale all decomposition and maintenance parameters in the model. An initial value of TF = 0.57 (parameter uncertainty, PU = 0.15) for NT (with TF set to 1.0 for CT) was used on the basis of a previous study with observations of soil CO₂ respiration. The simulated and observed changes in SOC were then compared using slopes of linear regressions of SOC changes over time. Results showed that the SOM model captured observed changes in SOC content from differences in rotations, N application and crop residue management for conventional tillage. On the basis of SOC change data a mean TF of 0.48 (standard deviation, SD = 0.12) was estimated for NT. The results indicate that (i) the estimated uncertainty of tillage effects on SOC turnover may be smaller than previously thought and (ii) simple scaling of SOM model parameters may be sufficient to capture the effects of soil tillage on SOM turnover in the tilled layer. Scenario analyses showed that the average extra C input needed to compensate for soil tillage was 762 (SD = 351) kg C ha⁻¹ year⁻¹. Climatic conditions (temperature and precipitation) also affected how much extra C was needed, with substantially larger inputs being required for wetter and warmer climates.     

Modeling tillage effects on soil physical properties

Tillage refers to the manipulation of soil by an implement powered by humans, animals or machines. Tillage operation generally create two zones: (1) a zone where soil has been fractured and then turned over leading to rough surface conditions; and (2) a zone where soil has been compacted by the weight of the machinery. Thus, modeling tillage effects on soil physical properties involves two separate approaches depending upon the zone under consideration.

Modeling tillage systems offers an opportunity to: (1) synthesize the extensive experimental data in the literature; (2) develop tools for site specific management recommendations; and (3) identify areas of research where additional information is needed. Modeling tillage systems involves modeling the soil physical, chemical and biological properties and processes and then linking them with crop growth models to simulate crop yields or environmental impacts. This paper reviews models for predicting tillage effects on state soil physical properties. Specifically, we reviewed models which predict bulk density, surface microrelief, aerodynamic roughness length, water retention characteristics, hydraulic conductivity function, thermal conductivity, volumetric heat capacity and gas diffusion coefficient. Since most of the existing models for predicting soil physical properties are developed for untilled soils, the paper outlines procedures to adapt these models to fractured and compacted zones in tilled soils. The paper also identifies specific assumptions that need both laboratory and field testing.     

Tillage management changes size-distribution of aggregates and macro-structure of soils used for irrigated row-crops

The Tatura system for the preparation of seed-beds for irrigated annual row-crops is described, where the soil is tilled when wet and friable and so requires few passes with implements to become suitable for crops, and where seeds are sown into wet soil. In soil prepared by the Tatura system, the percentage of aggregates < 0.5 mm diameter (as measured by dry-sieving) in the seed-bed was about half that found in commercially prepared seed-beds which were tilled up to 50 times when dry. With the Tatura system, the wetter the soil (up to 22% water content) when tilled, or the more passes (up to 4) of the implement at a water content of 22%, the less dust (< 0.5 mm diameter) and/or fewer clods (> 20 mm diameter) were formed.

The macro-structure of the surface layer of soil tilled at different water contents by the Tatura system was also quantified statistically by the method of wax-impregnation. The macro-structures were compared at the 10 mm, 20 mm, 40 mm, 60 mm and 80 mm depths in beds of soil prepared for irrigated annual row-crops by a system which has been described previously. Within each treatment (21.7%; 19.0%; 11.6% water content at 0–100 mm depth at tillage), the linear porosity and mean pore-size each tended to decrease with depth to 40 mm, with no further change or slight decrease to 80 mm depth. In all treatments, the mean aggregate-size tended to increase with depth from 10 mm depth to 80 mm depth. The sizes of pores and aggregates varied across each bed and possibly depended on the position of tines within the bed at each pass at tillage. Water content at tillage led to small differences in structure of the beds of soil. Soil tilled at a water content slightly above the Casagrande Plastic Limit generally had slightly larger pores and aggregates than soil tilled at lower water contents.     

Experimental evidence for the role of earthworms in compacted soil regeneration based on field observations and results from a semi-field experiment

In laboratory experiments, earthworms are often observed to burrow through compacted soil layers, leading to the general assumption that these animals play a significant role in regenerating compacted soils in agricultural plots. To demonstrate this role under field conditions, the abundance of earthworm macropores inside compacted zones was estimated on plots under reduced (RT) or conventional tillage (CT). Then, different types of compacted zones typically found in CT (plough pan and compacted clods) and RT plots (compacted volume under wheel tracks) were experimentally simulated in wooden boxes, buried in the field and inoculated with different earthworm species. After 6 weeks of incubation, the number of macropores inside the compacted zones was examined. Field observations showed that approximately 10% and 30% of the compacted zones were colonised by at least one macropore in CT and RT plots, respectively. A significantly greater number of anecics was found in RT plots, but we could not conclude that this ecological type of earthworm plays a more major role in the regeneration process in these plots since there were fewer compacted zones and these covered a smaller area in CT. The semi-field experiment provided evidence that earthworm-mediated regeneration of compacted zones is possible and its nature varies between ecological types of earthworm. Lumbricus terrestris, which makes individual burrows that are vertical and deep, was the main species to cross through the plough pan. The other three earthworm species (Aporrectodea giardi, A. caliginosa and A. rosea) did burrow inside the other types of compacted zones (“wheel tracks” and “compacted clods”). In every case, however, macropore density was far greater in non-compacted zones, illustrating that avoidance of compacted soil by earthworms is important and should be taken into account when extrapolating results from laboratory studies.     

Cotton management in a compacted subsurface microirrigated coastal plain soil of the southeastern US

A loamy sand Acrisol (Aquic Hapludult) that had been microirrigated for 6 years became so severely compacted that it had root limiting values of soil cone index in the Ap horizon and a genetic hardpan below it. Deep and surface tillage systems were evaluated for their ability to alleviate compaction. Deep tillage included subsoiling or none. Both deep tillage treatments were also surface tilled by disking, chiseling, or not tilling. Subsoiling was located in row or between rows to avoid microirrigation tubes (laterals) that were buried under every other mid row or every row. Cotton (Gossypium hirsutum) was planted in 0.96-m wide rows. Cotton yield was improved by irrigation from 485 to 1022 kg ha⁻¹ because both 2001 and 2002 were dry years. Tillage loosened the soil by an average of 0.5–1.3 MPa; but compacted zones remained outside tilled areas. Subsoiling improved yield by 131 kg ha⁻¹ when performed in row where laterals were placed in the mid rows; but subsoiling did not improve yield when it was performed in mid rows. For subsurface irrigation management in these soils, the treatment with laterals buried under every other mid row was able to accommodate in-row subsoiling which improved yield; and this treatment was just as productive as and had been shown to be less expensive to install than burying laterals under every row.     

TILTH MELLOWING

Effects of weathering action, mainly wetting and drying cycles, on the strength of the clods produced by tillage are studied. Experiments were carried out on sandy loam soils at two sites in South Australia, and on silt loam and clay soils at Wye College, England. It is found that tillage increases the amplitude of soil water content fluctuations. These bigger soil water content fluctuations resulted in a decrease in the clod strength and this in turn modified the size distribution of the clods produced by tillage in the South Australian soils. The decrease in clod strength, as measured by the drop shatter test, was followed by an increase in the proportion of the smaller aggregate size fraction produced by a second implement pass. It is suggested that, for soils in which the increase in the soil water content fluctuations after the first tillage implement pass decreases clod strength, a further implement pass should be delayed for several days. By doing this, the soil can be tilled with minimum energy and cost to produce a good seed bed.     

The structure of two alluvial soils in Italy after 10 years of conventional and minimum tillage

Micro and macroporosity, pore shape and size distribution, aggregate stability, saturated hydraulic conductivity and crop yield were analysed in alluvial silty loam (Fluventic Eutrochrept) and clay soils (Vertic Eutrochrept) following long-term minimum and conventional tillage. The soil structure attributes were evaluated by characterizing porosity by means of image analysis of soil thin sections prepared from undisturbed soil samples.

The interaggregate microporosity, measured by mercury intrusion porosimetry, increased in the minimally tilled soils, with a particular increase in the storage pores (0.5–50 μm). The amount of elongated transmission pores (50–500 μm) also increased in the minimally tilled soils. The resulting soil structure was more open and more homogeneous, thus allowing better water movement, as confirmed by the greater hydraulic conductivity of the minimally tilled soils. The aggregate stability was less in the conventionally tilled soils and this resulted in a greater tendency to form surface crusts and compacted structure, compared with the minimally tilled soils. The latter tillage practice seemed to maintain, in the long-term, better soil structure conditions and, therefore, maintain favourable conditions for plant growth. In the silt loam, the crop yield did not differ significantly between the two tillage systems, while in the clay soil it decreased in the minimum tilled soil because of problems of seed bed preparation at the higher surface layer water content.     

Effects of autumn tillage and residue management on soil respiration in a long‐term field experiment in Sweden

Several previous field studies in temperate regions have shown decreased soil respiration after conventional tillage compared with reduced or no‐tillage treatments. Whether this decrease is due to differences in plant residue distribution or changes in soil structure following tillage remains an open question. This study investigated (1) the effects of residue management and incorporation depth on soil respiration and (2) biological activity in different post‐tillage aggregates representing the actual size and distribution of aggregates observed in the tilled layer. The study was conducted within a long‐term tillage experiment on a clay soil (Eutric Cambisol) in Uppsala, Sweden. After 38 y, four replicate plots in two long‐term treatments (moldboard plowing (MP) and shallow tillage (ST)) were split into three subplots. These were then used for a short‐term trial in which crop residues were either removed, left on the surface or incorporated to about 6 cm depth (ST) or at 20 cm depth (MP). Soil respiration, soil temperature, and water content were monitored during a 10‐d period after tillage treatment. Respiration from aggregates of different sizes produced by ST and MP was also measured at constant water potential and temperature in the laboratory. The results showed that MP decreased short‐term soil respiration compared with ST or no tillage. Small aggregates (< 16 mm) were biologically most active, irrespective of tillage method, but due to their low proportion of total soil mass they contributed < 1.5% to total respiration from the tilled layer. Differences in respiration between tillage treatments were found to be attributable to indirect effects on soil moisture and temperature profiles and the depth distribution of crop residues, rather than to physical disturbance of the soil.     

含水量和容重对旱地耕层土壤热导率的影响及预测

土壤热导率是研究地表能量平衡和土壤水热运移过程中的一个基础参数。受土壤耕作、干湿交替和根系生长等过程的影响,耕层土壤的含水率和结构呈现较强的变异特征,而目前缺乏关于定量分析耕层土壤热导率变异特征的研究。该研究利用田间定位试验,采用热脉冲技术测定了含水率和容重变化条件下耕层土壤热导率的变异特征,并利用传递函数模型对耕层土壤热导率进行了预测。结果表明:含水率和容重是影响耕层土壤热导率变异的主要因子,而耕作强度和干湿交替是这种变异的关键驱动力;与翻耕和旋耕处理相比,免耕处理提高了土壤容重和含水率,从而增大了土壤热导率;在干湿交替作用下,翻耕后土壤容重逐步增加,耕层热导率也呈现上升趋势,波动幅度与含水率的变化相关。基于含水率、容重和质地信息,土壤热导率传递函数模型可以给出可靠的田间土壤热导率估计值,其均方根误差和平均偏差分别为0.09和-0.01 W/(m·K);考虑耕层土壤容重的动态信息,可以提高该模型预测土壤热导率的准确性。     

Assessment of soil erosion indicators for maize-based agro-ecosystems in Kenya

Measuring soil loss is costly, must cover a range of field situations, is not standardized, and is season dependent. In addition, use of sparse soil loss data (from other studies) compromises the integrity of many erosion models. Easily assessable soil erosion indicators to monitor the cumulative effect of erosion between tillage/weeding and harvesting called eroding clods, flow surfaces, pre-rills, and rills were surveyed directly after the 1995 rainy season in the Taita Taveta district of Kenya, to assessed the utility of each indicator. Their incidences were modeled using CPA. In the area, 70 maize plots in 11 map units, having considerable variation in altitude, land cover, rainfall, and geomorphology, were surveyed. Soil loss was considered variable between plots due to differences in surface soil, land cover, infrastructure (trash lines, grass strips, and Fanya-juu), crop management, slope, and map unit. The eroding clods indicator proved of little significance because the initial clods cover was unknown; the indicator probably relates better to soil erodibility then to soil loss. Flow surfaces, formed during erosive showers, were less present on fields with a higher ground and canopy cover, if the area of eroding clods was high, and if the topsoil had no loam which reduces chances of sealing; no impact of infrastructure, tillage, and weeding were detected. Fewer pre-rills were present where the fraction of groundcover was high, where Pigeon Peas were not grown (they cause micro-relief and concentrated flows), where weeding ended late (time effect), where more flow surfaces occurred, where Fanya-juu was constructed (less steep slopes), where the top-soil contained little sand (less sediment entrainment), and where maize was intercropped with vegetables (positive canopy cover effects). The model was not map unit specific and had an Adjusted R² of 67%. The log-linear relationship indicates that combined positive conditions exponentially reduce the occurrence of pre-rills. The “pre-rill” indicator related best to management-affected site conditions and seems to reflect best the cumulative effects of soil loss over time. Rills were found at 18 sites located in drier areas on sandy–clay soils. The model suggested more rills if the topsoil contains no silt; this makes the soil susceptible to compaction, peptisation when wet, and rill formation.     

Least limiting water range for different soil management practices in dryland farming in Iran

Integrated evaluation of soil physical properties using the least limiting water range (LLWR) approach may allow a better knowledge of soil water availability. We determined the LLWR for four tillage practices consisted of conventional tillage (CT), reduced tillage (RT), no-tillage (NT) and fallow no-tillage (NT_f). In addition, LLWR was determined for abandoned soils (i.e. control), compacted soils, ploughed compacted soils and abandoned soils with super absorbent polymers (SAPs) application. Soil water retention, penetration resistance (PR), air-filled porosity and bulk density were determined for the 0–5 and 0–25-cm depths. Mean LLWR (0.07–0.08 cm³ cm^?3) was lower in compacted soils than the soils under CT, NT, NT_f, RT, tilled, abandoned and SAP practices but it was not different among tillage practices. The values of LLWR were 0.12 cm³ cm^?3 for NT and CT. LLWR for tilled plots (0.12 cm³ cm^?3) became greater than compacted soils by 1.3 times. Analysis of the lower and upper limits of the LLWR further indicated that PR was the only limiting factor for soil water content, but aeration was not a limiting factor. The LLWR was more dependent on soil water content at permanent wilting point and at PR.     

Reassessment of tillage erosion rates by manual tillage on steep slopes in northern Thailand

Changing land-use practices in northern Thailand have increased tillage intensity. This study re-assesses the rate of tillage erosion by manual hoeing on steep slopes (17–82%) in northern Thailand. Previously collected soil translocation data during an on-farm tillage erosion experiment and additionally collected data during an on-farm tillage erosion survey have been analysed whereby a new calculation method (i.e. trapezoid tillage step) has been used. A comparison with previously collected data indicates that the trapezoid tillage step method and the tracer method are the most reliable methods to assess downslope translocation by manual tillage. Based on newly acquired understanding of the processes involved, soil fluxes by tillage erosion are quantified by linear functions for different slope gradient classes rather than one single diffusion-type equation for the whole slope range. For slope gradients smaller than 3%, soil fluxes are close to zero as farmers do not have a preferred tillage direction. For slope gradients between 3% and 70%, soil is tilled only in the downslope direction and soil fluxes range between 16 and 67 kg m⁻¹ tillage pass⁻¹. On slopes with gradients in excess of 70%, the angle of repose for soil clods is often exceeded resulting in a sliding down of the complete tilled top layer. These data are used to assess the soil flux for complete cropping cycles for the most dominant cropping systems in the highlands of northern Thailand: i.e. upland rice, maize, (soy) beans, cabbage and ginger. The on-site effects of tillage erosion will be very pronounced if parcels are short with respect to their slope length, cultivated for upland rice or cabbage, or when weed pressure is high. Tillage erosion results in a tillage step with low soil fertility and low infiltration capacity. Solutions to reduce tillage erosion intensity depend on the degree that tillage intensity can be reduced. This might happen by an improved weed management or by changing landuse to perrenial cropping. Other strategies are concentrating nutrients on the truncated hillslope sections and retaining soil on the field by vegetative buffers.     

Spatial variability of short-term effect of tillage on soil penetration resistance

The effect of tillage on soil properties varies within field due to spatial variability of soils. Mapping changes of soil penetration resistance (PR) would be useful to understand and assess tillage practices to alleviate soil compaction. The objectives were to model the short-term effect of tillage on PR and its spatial structure, and to delineate homogeneous zones based on soil response in a Typic Argiudoll previously managed under no-till. A grid sampling for PR and soil water content (SWC) were performed before and after chiselling. Spatial analysis was performed on the effect of tillage on PR data by 10 cm layers and homogeneous zones were delineated by k-means cluster analysis. The effect of tillage was ?0.33 MPa in 10–20 and 20–30 cm layer. No differences of PR were found at 0–10 cm. Short range (5–7 m) spatial structure on the horizontal plane was observed in all layers. Only 45% of the field showed a marked effect of tillage on PR. Mapping the effect of tillage on PR would be a useful approach for evaluating the global and local response of soil to tillage, as well as for delineating of areas within field for site-specific tillage practices.     

Electrical resistivity tomography to detect the effects of tillage in a soil with a variable rock fragment content

Electrical resistivity tomography (ERT) is a promising non‐destructive tool to characterize agricultural soils where management effects are superimposed on natural variability. The aim of our study was to test whether ERT was capable of detecting stones and tillage effects in a soil with a variable rock fragment content. Field experiments were conducted by performing a set of three two‐dimensional (2D) resistivity tomographies across two management systems (tillage/no tillage) replicated twice on each transect, using dipole‐dipole configuration and 0.25‐m inter‐electrode spacing. Soil texture, bulk density and water content were measured destructively. Greater average electrical resistivity (ER) was found in tilled plots, with maximum values of up to 1700 Ohm m. However, when the spatial correlation structure was considered in a mixed‐effects model, no significant difference in ER was found between tilled and untilled plots. Empirical semivariograms showed less spatial continuity and more noise in tilled plots. Resistivity was strongly correlated with rock fragment content (r = 0.68), with greater average values in ploughed plots, which may possibly be linked to kinetic sieving after ploughing. ERT was able to identify the position of gravel lenses and was also sensitive to the presence of clay (r = ?0.45): a linear trend in resistivity across the field (r = 0.80) was consistent with a decreasing clay content (r = ?0.68). Resistivity was correlated with rock fragments, clay and an interaction variable (water × rock fragments). There was a poor fit for the tilled plot where resistivity peaks could be linked to the presence of voids, but their detection would have required a resolution greater than that which we adopted.