Reconsolidation of the soil surface after tillage discontinuity, with and without cultivation, related to erosion and its prediction with RUSLE

Site-specific regression coefficient values are essential for erosion prediction with empirical models. With the objective to investigate the surface-soilconsolidation factor, Cf, linked to the RUSLE's prior-land-use subfactor, PLU, an erosion experiment using simulated rainfall on a 0.075 m m-1 slope, sandy loam Paleudult soil, was conducted at the Agriculture Experimental Station of the Federal University of Rio Grande do Sul (EEA/UFRGS), in Eldorado do Sul, State of Rio Grande do Sul, Brazil. Firstly, a row-cropped area was excluded from cultivation (March 1995), the existing crop residue removed from the field, and the soil kept clean-tilled the rest of the year (to get a degraded soil condition for the intended purpose of this research). The soil was then conventional-tilled for the last time (except for a standard plot which was kept continuously cleantilled for comparison purposes), in January 1996, and the following treatments were established and evaluated for soil reconsolidation and soil erosion until May 1998, on duplicated 3.5 x 11.0 m erosion plots: (a) fresh-tilled soil, continuously in clean-tilled fallow (unit plot); (b) reconsolidating soil without cultivation; and (c) reconsolidating soil with cultivation (a crop sequence of three corn- and two black oats cycles, continuously in no-till, removing the crop residues after each harvest for rainfall application and redistributing them on the site after that). Simulated rainfall was applied with a Swanson's type, rotating-boom rainfall simulator, at 63.5 mm h-1 intensity and 90 min duration, six times during the two-and-half years of experimental period (at the beginning of the study and after each crop harvest, with the soil in the unit plot being retilled before each rainfall test). The soil-surface-consolidation factor, Cf, was calculated by dividing soil loss values from the reconsolidating soil treatments by the average value from the fresh-tilled soil treatment (unit plot). Non-linear regression was used to fit the Cf = e b.t model through the calculated Cf-data, where t is time in days since last tillage. Values for b were -0.0020 for the reconsolidating soil without cultivation and -0.0031 for the one with cultivation, yielding Cf-values equal to 0.16 and 0.06, respectively, after two-and-half years of tillage discontinuation, compared to 1.0 for fresh-tilled soil. These estimated Cf-values correspond, respectively, to soil loss reductions of 84 and 94 %, in relation to soil loss from the fresh-tilled soil, showing that the soil surface reconsolidated intenser with cultivation than without it. Two distinct treatmentinherent soil surface conditions probably influenced the rapid decay-rate of Cf values in this study, but, as a matter of a fact, they were part of the real environmental field conditions. Cf-factor curves presented in this paper are therefore useful for predicting erosion with RUSLE, but their application is restricted to situations where both soil type and particular soil surface condition are similar to the ones investigate in this study.

Relations between soil surface roughness, tortuosity, tillage treatments, rainfall intensity and soil and water losses from a red yellow latosol

The soil surface roughness increases water retention and infiltration, reduces the runoff volume and speed and influences soil losses by water erosion. Similarly to other parameters, soil roughness is affected by the tillage system and rainfall volume. Based on these assumptions, the main purpose of this study was to evaluate the effect of tillage treatments on soil surface roughness (RR) and tortuosity (T) and to investigate the relationship with soil and water losses in a series of simulated rainfall events. The field study was carried out at the experimental station of EMBRAPA Southeastern Cattle Research Center in São Carlos (Fazenda Canchim), in São Paulo State, Brazil. Experimental plots of 33 m² were treated with two tillage practices in three replications, consisting of: untilled (no-tillage) soil (NTS) and conventionally tilled (plowing plus double disking) soil (CTS). Three successive simulated rain tests were applied in 24 h intervals. The three tests consisted of a first rain of 30 mm/h, a second of 30 mm/h and a third rain of 70 mm/h. Immediately after tilling and each rain simulation test, the surface roughness was measured, using a laser profile meter. The tillage treatments induced significant changes in soil surface roughness and tortuosity, demonstrating the importance of the tillage system for the physical surface conditions, favoring water retention and infiltration in the soil. The increase in surface roughness by the tillage treatments was considerably greater than its reduction by rain action. The surface roughness and tortuosity had more influence on the soil volume lost by surface runoff than in the conventional treatment. Possibly, other variables influenced soil and water losses from the no-tillage treatments, e.g., soil type, declivity, slope length, among others not analyzed in this study.

Soil fertility, nutrition and yield of maize and barley with gypsum application on soil surface in no-till

Annual crop yield and nutrition have shown differentiated responses to modifications in soil chemical properties brought about by gypsum application. The aim of this study was to evaluate the effect of gypsum application rates on the chemical properties of a Latossolo Bruno (Clayey Oxisol), as well as on the nutrition and yield of a maize-barley succession under no-till. The experiment was set up in November 2009 in Guarapuava, Parana, Brazil, applying gypsum rates of 0.0, 1.5, 3.0, 4.5, and 6.0 Mg ha-1 to the soil surface upon sowing maize, with crop succession of barley. Gypsum application decreased the levels of Al3+ and Mg2+ in the 0.0-0.1 m layer and increased soil pH in the layers from 0.2-0.6 m depth. Gypsum application has increased the levels of Ca2+ in all soil layers up to 0.6 m, and the levels of S-SO4(2-) up to 0.8 m. In both crops, the leaf concentrations of Ca and S were increased while Mg concentrations have decreased as a function of gypsum rates. There was also an effect of gypsum rates on grain yield, with a quadratic response of maize and a linear increase for barley. Yield increases were up to 11 and 12 % in relation to control for the maximum technical efficiency (MTE) rates of 3.8 and 6.0 Mg ha-1 of gypsum, respectively. Gypsum application improved soil fertility in the profile, especially in the subsurface, as well as plant nutrition, increasing the yields of maize and barley.

Characterization of the soil fertility and root system of restinga forests

The Restinga vegetation consists of a mosaic of plant communities, which are defined by the characteristics of the substrates, resulting from the type and age of the depositional processes. This mosaic complex of vegetation types comprises restinga forest in advanced (high restinga) and medium regeneration stages (low restinga), each with particular differentiating vegetation characteristics. The climate along the coast is tropical (Köppen). Of all ecosystems of the Atlantic Forest, Restinga is the most fragile and susceptible to anthropic disturbances. Plants respond to soil characteristics with physiological and morphological modifications, resulting in changes in the architecture (spatial configuration) of the root system. The purpose of this study was to characterize the soil fertility of high and low restinga forests, by chemical and physical parameters, and its relation to the root system distribution in the soil profile. Four locations were studied: (1) Ilha Anchieta State Park, Ubatuba; (2) two Ecological Stations of Jureia-Itatins and of Chauás, in the municipality of Iguape; (3) Vila de Pedrinhas in the municipality of Ilha Comprida; and (4) Ilha do Cardoso State Park, Cananeia. The soil fertility (chemical and physical properties) was analyzed in the layers 0-5, 0-10, 0-20, 20-40 and 40-60 cm. In addition, the distribution of the root system in the soil profile was evaluated, using digital images and the Spring program. It was concluded that the root system of all vegetation types studied is restricted to the surface layers, 0-10 and 10-20 cm, but occupies mainly the 0-10 cm layer (70 %); that soil fertility is low in all environments studied, with base saturation values below 16 %, since most exchange sites are occupied by aluminum; and that restinga vegetation is edaphic.

Variability of apparently homogeneous soilscapes in São Paulo state, Brazil: I. spatial analysis

The spatial variability of strongly weathered soils under sugarcane and soybean/wheat rotation was quantitatively assessed on 33 fields in two regions in São Paulo State, Brazil: Araras (15 fields with sugarcane) and Assis (11 fields with sugarcane and seven fields with soybean/wheat rotation). Statistical methods used were: nested analysis of variance (for 11 fields), semivariance analysis and analysis of variance within and between fields. Spatial levels from 50 m to several km were analyzed. Results are discussed with reference to a previously published study carried out in the surroundings of Passo Fundo (RS). Similar variability patterns were found for clay content, organic C content and cation exchange capacity. The fields studied are quite homogeneous with respect to these relatively stable soil characteristics. Spatial variability of other characteristics (resin extractable P, pH, base- and Al-saturation and also soil colour), varies with region and, or land use management. Soil management for sugarcane seems to have induced modifications to greater depths than for soybean/wheat rotation. Surface layers of soils under soybean/wheat present relatively little variation, apparently as a result of very intensive soil management. The major part of within-field variation occurs at short distances (< 50 m) in all study areas. Hence, little extra information would be gained by increasing sampling density from, say, 1/km² to 1/50 m². For many purposes, the soils in the study regions can be mapped with the same observation density, but residual variance will not be the same in all areas. Bulk sampling may help to reveal spatial patterns between 50 and 1.000 m.

Soil fertility, organic carbon and fractions of the organic matter at different distances from eucalyptus stumps

Knowledge on variations in vertical, horizontal and temporal characteristics of the soil chemical properties under eucalyptus stumps left in the soil is of fundamental importance for the management of subsequent crops. The objective of this work was to evaluate the effect of eucalyptus stumps (ES) left after cutting on the spatial variability of chemical characteristics in a dystrophic Yellow Argisol in the eastern coastal plain region of Brazil. For this purpose, ES left for 31 and 54 months were selected in two experimental areas with similar characteristics, to assess the decomposition effects of the stumps on soil chemical attributes. Soil samples were collected directly around these ES, and at distances of 30, 60, 90, 120 and 150 cm away from them, in the layers 0-10, 10-20 and 20-40 cm along the row of ES, which is in-between the rows of eucalyptus trees of a new plantation, grown at a spacing of 3 x 3 m. The soil was sampled in five replications in plots of 900 m² each and the samples analyzed for pH, available P and K (Mehlich-1), exchangeable Al, Ca and Mg, total organic carbon (TOC) and C content in humic substances (HS) and in the free light fraction. The pH values and P, K, Ca2+, Mg2+ and Al3+ contents varied between the soil layers with increasing distance from the 31 and 54-monthold stumps. The highest pH, P, K, Ca2+ and Mg2+ values and the lowest Al3+ content were found in the surface soil layer. The TOC of the various fractions of soil organic matter decreased with increasing distance from the 31 and 54-month-old ES in the 0-10 and 10-20 cm layers, indicating that the root (and stump) cycling and rhizodeposition contribute to maintain soil organic matter. The C contents of the free light fraction, of the HS and TOC fractions were higher in the topsoil layer under the ES left for 31 months due to the higher clay levels of this layer, than in those found under the 54-month-old stumps. However, highest C levels of the different fractions of soil organic matter in the topsoil layer reflect the deposition and maintenance of forest residues on the soil surface, mainly after forest harvest.

Temporal variability of soil water content under different surface conditions in the semiarid region of the Pernambuco state

Rainfall in the semiarid region of Pernambuco is characterized by irregular distribution in time and space, which significantly hinders the rainfed agriculture in the region. This work aims to evaluate the temporal profile of soil moisture in the semiarid region of the Pernambuco State (Brazil) and the effect of different soil surface conditions on soil water content variation and the yield of rainfed beans. To monitor soil water content, five plots 4.5 m wide by 11 m long were installed in a Yellow Argisol (Ultisol). The following treatments were adopted in the experimental plots: natural vegetation, bean intercropped with cactus, beans planted down the slope, beans planted along contour lines with mulch and rock barriers, and bare soil. In each plot, eight PVC access tubes were installed for monitoring the soil water content profile at depths of 0.20 and 0.40 m using a neutron probe device. The surface condition significantly influenced the soil water content variation, both in the dry and rainy seasons. The use of mulch, associated with rock barriers, provided higher soil water content levels than the other treatments and increased the rainfed beans production.

New methods to quantify NH3 volatilization from fertilized surface soil with urea

Gaseous N losses from soil are considerable, resulting mostly from ammonia volatilization linked to agricultural activities such as pasture fertilization. The use of simple and accessible measurement methods of such losses is fundamental in the evaluation of the N cycle in agricultural systems. The purpose of this study was to evaluate quantification methods of NH3 volatilization from fertilized surface soil with urea, with minimal influence on the volatilization processes. The greenhouse experiment was arranged in a completely randomized design with 13 treatments and five replications, with the following treatments: (1) Polyurethane foam (density 20 kg m-3) with phosphoric acid solution absorber (foam absorber), installed 1, 5, 10 and 20 cm above the soil surface; (2) Paper filter with sulfuric acid solution absorber (paper absorber, 1, 5, 10 and 20 cm above the soil surface); (3) Sulfuric acid solution absorber (1, 5 and 10 cm above the soil surface); (4) Semi-open static collector; (5) 15N balance (control). The foam absorber placed 1 cm above the soil surface estimated the real daily rate of loss and accumulated loss of NH3N and proved efficient in capturing NH3 volatized from urea-treated soil. The estimates based on acid absorbers 1, 5 and 10 cm above the soil surface and paper absorbers 1 and 5 cm above the soil surface were only realistic for accumulated N-NH3 losses. Foam absorbers can be indicated to quantify accumulated and daily rates of NH3 volatilization losses similarly to an open static chamber, making calibration equations or correction factors unnecessary.

Spatial patterns of preconsolidation pressure and soil moisture along transects in two directions under coffee

Information on the spatial structure of soil physical and structural properties is needed to evaluate the soil quality. The purpose of this study was to investigate the spatial behavior of preconsolidation pressure and soil moisture in six transects, three selected along and three across coffee rows, at three different sites under different tillage management systems. The study was carried out on a farm, in Patrocinio, state of Minas Gerais, in the Southeast of Brazil (18 º 59 ' 15 '' S; 46 º 56 ' 47 '' W; 934 m asl). The soil type is a typic dystrophic Red Latosol (Acrustox) and consists of 780 g kg-1 clay; 110 g kg-1 silt and 110 g kg-1 sand, with an average slope of 3 %. Undisturbed soil cores were sampled at a depth of 0.10-0.13 m, at three different points within the coffee plantation: (a) from under the wheel track, where equipment used in farm operations passes; (b) in - between tracks and (c) under the coffee canopy. Six linear transects were established in the experimental area: three transects along and three across the coffee rows. This way, 161 samples were collected in the transect across the coffee rows, from the three locations, while 117 samples were collected in the direction along the row. The shortest sampling distance in the transect across the row was 4 m, and 0.5 m for the transect along the row. No clear patterns of the preconsolidation pressure values were observed in the 200 m transect. The results of the semivariograms for both variables indicated a high nugget value and short range for the studied parameters of all transects. A cyclic pattern of the parameters was observed for the across-rows transect. An inverse relationship between preconsolidation pressure and soil moisture was clearly observed in the samples from under the track, in both directions.

Impacts of land leveling on lowland soil physical properties

The practice of land leveling alters the soil surface to create a uniform slope to improve land conditions for the application of all agricultural practices. The aims of this study were to evaluate the impacts of land leveling through the magnitudes, variances and spatial distributions of selected soil physical properties of a lowland area in the State of Rio Grande do Sul, Brazil; the relationships between the magnitude of cuts and/or fills and soil physical properties after the leveling process; and evaluation of the effect of leveling on the spatial distribution of the top of the B horizon in relation to the soil surface. In the 0-0.20 m layer, a 100-point geo-referenced grid covering two taxonomic soil classes was used in assessment of the following soil properties: soil particle density (Pd) and bulk density (Bd); total porosity (Tp), macroporosity (Macro) and microporosity (Micro); available water capacity (AWC); sand, silt, clay, and dispersed clay in water (Disp clay) contents; electrical conductivity (EC); and weighted average diameter of aggregates (WAD). Soil depth to the top of the B horizon was also measured before leveling. The overall effect of leveling on selected soil physical properties was evaluated by paired "t" tests. The effect on the variability of each property was evaluated through the homogeneity of variance test. The thematic maps constructed by kriging or by the inverse of the square of the distances were visually analyzed to evaluate the effect of leveling on the spatial distribution of the properties and of the top of the B horizon in relation to the soil surface. Linear regression models were fitted with the aim of evaluating the relationship between soil properties and the magnitude of cuts and fills. Leveling altered the mean value of several soil properties and the agronomic effect was negative. The mean values of Bd and Disp clay increased and Tp, Macro and Micro, WAD, AWC and EC decreased. Spatial distributions of all soil physical properties changed as a result of leveling and its effect on all soil physical properties occurred in the whole area and not specifically in the cutting or filling areas. In future designs of leveling, we recommend overlaying a cut/fill map on the map of soil depth to the top of the B horizon in order to minimize areas with shallow surface soil after leveling.

EFFECT OF SOIL TILLAGE AND PLANT RESIDUE ON SURFACE ROUGHNESS OF AN OXISOL UNDER SIMULATED RAIN

Surface roughness of the soil is formed by mechanical tillage and is also influenced by the kind and amount of plant residue, among other factors. Its persistence over time mainly depends on the fundamental characteristics of rain and soil type. However, few studies have been developed to evaluate these factors in Latossolos (Oxisols). In this study, we evaluated the effect of soil tillage and of amounts of plant residue on surface roughness of an Oxisol under simulated rain. Treatments consisted of the combination of the tillage systems of no-tillage (NT), conventional tillage (CT), and minimum tillage (MT) with rates of plant residue of 0, 1, and 2 Mg ha-1 of oats (Avena strigosa Schreb) and 0, 3, and 6 Mg ha-1 of maize (Zea mays L.). Seven simulated rains were applied on each experimental plot, with intensity of 60±2 mm h-1 and duration of 1 h at weekly intervals. The values of the random roughness index ranged from 2.94 to 17.71 mm in oats, and from 5.91 to 20.37 mm in maize, showing that CT and MT are effective in increasing soil surface roughness. It was seen that soil tillage operations carried out with the chisel plow and the leveling disk harrow are more effective in increasing soil roughness than those carried out with the heavy disk harrow and leveling disk harrow. The roughness index of the soil surface decreases exponentially with the increase in the rainfall volume applied under conditions of no tillage without soil cover, conventional tillage, and minimum tillage. The oat and maize crop residue present on the soil surface is effective in maintaining the roughness of the soil surface under no-tillage.

SPATIAL VARIABILITY OF SOIL PROPERTIES IN AN AGRARIAN REFORM SETTLEMENT

ABSTRACT The study of soil chemical and physical properties variability is important for suitable management practices. The aim of this study was to evaluate the spatial variability of soil properties in the Malhada do Meio settlement to subsidize soil use planning. The settlement is located in Chapadinha, MA, Brazil, and has an area of 630.86 ha. The vegetation is seasonal submontane deciduous forest and steppe savanna. The geology is formed of sandstones and siltstones of theItapecuru Formation and by colluvial and alluvial deposits. The relief consists of hills with rounded and flat tops with an average altitude of 67 m, and frequently covered over by ferruginous duricrusts. A total of 183 georeferenced soil samples were collected at the depth of 0.00-0.20 m inPlintossolos, Neossolo andGleissolo. The following chemical variables were analyzed: pH(CaCl2), H+Al, Al, SB, V, CEC, P, K, OM, Ca, Mg, SiO2, Al2O3, and Fe2O3; along with particle size variables: clay, silt, and sand. Descriptive statistical and geostatistical analyses were carried out. The coefficient of variation (CV) was high for most of the variables, with the exception of pH with a low CV, and of sand with a medium CV. The models fitted to the experimental semivariograms of these variables were the exponential and the spherical. The range values were from 999 m to 3,690 m. For the variables pH(CaCl2), SB, and clay, there are three specific areas for land use planning. The central part of the area (zone III), where thePlintossolos Pétricos and Neossolos Flúvicos occur, is the most suitable for crops due to higher macronutrient content, organic matter and pH. Zones I and II are indicated for environmental preservation.

Surface and subsurface decomposition of a desiccated grass pasture biomass related to erosion and its prediction with RUSLE

Erosion is deleterious because it reduces the soil's productivity capacity for growing crops and causes sedimentation and water pollution problems. Surface and buried crop residue, as well as live and dead plant roots, play an important role in erosion control. An efficient way to assess the effectiveness of such materials in erosion reduction is by means of decomposition constants as used within the Revised Universal Soil Loss Equation - RUSLE's prior-land-use subfactor - PLU. This was investigated using simulated rainfall on a 0.12 m m-1 slope, sandy loam Paleudult soil, at the Agriculture Experimental Station of the Federal University of Rio Grande do Sul, in Eldorado do Sul, State of Rio Grande do Sul, Brazil. The study area had been covered by native grass pasture for about fifteen years. By the middle of March 1996, the sod was mechanically mowed and the crop residue removed from the field. Late in April 1996, the sod was chemically desiccated with herbicide and, about one month later, the following treatments were established and evaluated for sod biomass decomposition and soil erosion, from June 1996 to May 1998, on duplicated 3.5 x 11.0 m erosion plots: (a) and (b) soil without tillage, with surface residue and dead roots; (c) soil without tillage, with dead roots only; (d) soil tilled conventionally every two-and-half months, with dead roots plus incorporated residue; and (e) soil tilled conventionally every six months, with dead roots plus incorporated residue. Simulated rainfall was applied with a rotating-boom rainfall simulator, at an intensity of 63.5 mm h-1 for 90 min, eight to nine times during the experimental period (about every two-and-half months). Surface and subsurface sod biomass amounts were measured before each rainfall test along with the erosion measurements of runoff rate, sediment concentration in runoff, soil loss rate, and total soil loss. Non-linear regression analysis was performed using an exponential and a power model. Surface sod biomass decomposition was better depicted by the exponential model, while subsurface sod biomass was by the power model. Subsurface sod biomass decomposed faster and more than surface sod biomass, with dead roots in untilled soil without residue on the surface decomposing more than dead roots in untilled soil with surface residue. Tillage type and frequency did not appreciably influence subsurface sod biomass decomposition. Soil loss rates increased greatly with both surface sod biomass decomposition and decomposition of subsurface sod biomass in the conventionally tilled soil, but they were minimally affected by subsurface sod biomass decomposition in the untilled soil. Runoff rates were little affected by the studied treatments. Dead roots plus incorporated residues were effective in reducing erosion in the conventionally tilled soil, while consolidation of the soil surface was important in no-till. The residual effect of the turned soil on erosion diminished gradually with time and ceased after two years.

Nitrogen application timing and soil inorganic nitrogen dynamics under no-till oat/maize sequential cropping

The timing of N application to maize is a key factor to be considered in no-till oat/maize sequential cropping. This study aimed to evaluate the influence of pre-planting, planting and sidedress N application on oat residue decomposition, on soil N immobilisation and remineralisation and on N uptake by maize plants in no-till oat/maize sequential cropping. Undisturbed soil cores of 10 and 20 cm diameter were collected from the 0-15 cm layer of a no-till Red Latossol, when the oat cover crop was in the milk-grain stage. Two greenhouse experiments were conducted simultaneously. Experiment A, established in the 10 cm diameter cores and without plant cultivation, was used to asses N dynamics in soil and oat residues. Experiment B, established in the 20 cm diameter cores and with maize cultivation, was used to assess plant growth and N uptake. An amount of 6.0 Mg ha-1 dry matter of oat residues was spread on the surface of the cores. A rate of 90 kg N ha-1 applied as ammonium sulphate in both experiments was split in pre-planting, planting and sidedress applications as follows: (a) 00-00-00 (control), (b) 90-00-00 (pre-planting application, 20 days before planting), (c) 00-90-00 (planting application), (d) 00-30-60 (split in a planting and a sidedress application 31 days after emergence), (e) 00-00-00* (control, without oat residue) and (f) 90-00-00* (pre-planting application, without oat residue). The N concentration and N content in oat residues were not affected during decomposition by N fertilisation. Most of the fertiliser NH4+-N was converted into NO3--N within 20 days after application. A significant decrease in NO3--N contents in the 0-4 cm layer was observed in all treatments between 40 and 60 days after the oat residue placement on the soil surface, suggesting the occurrence of N immobilisation in this period. Considering that most of the inorganic N was converted into NO3- and that no immobilisation of the pre planting fertiliser N occurred at the time of its application, it was possible to conclude that pre-planting applied N was prone to losses by leaching. On the other hand, with split N applications, maize plants showed N deficiency symptoms before sidedress application. Two indications for fertiliser-N management in no-till oat/maize sequential cropping could be suggested: (a) in case of split application, the sidedress should be earlier than 30 days after emergence, and (b) if integral application is preferred to save field operations, this should be done at planting.

Numerical analysis of the impact of charcoal production on soil hydrological behavior, runoff response and erosion susceptibility

The impact of charcoal production on soil hydraulic properties, runoff response and erosion susceptibility were studied in both field and simulation experiments. Core and composite samples, from 12 randomly selected sites within the catchment of Kotokosu were taken from the 0-10 cm layer of a charcoal site soil (CSS) and adjacent field soils (AFS). These samples were used to determine saturated hydraulic conductivity (Ksat), bulk density, total porosity, soil texture and color. Infiltration, surface albedo and soil surface temperature were also measured in both CSS and AFS. Measured properties were used as entries in a rainfall runoff simulation experiment on a smooth (5 % slope) plot of 25 x 25 m grids with 10 cm resolutions. Typical rainfall intensities of the study watershed (high, moderate and low) were applied to five different combinations of Ks distributions that could be expected in this landscape. The results showed significantly (p < 0.01) higher flow characteristics of the soil under charcoal kilns (increase of 88 %). Infiltration was enhanced and runoff volume reduced significantly. The results showed runoff reduction of about 37 and 18 %, and runoff coefficient ranging from 0.47-0.75 and 0.04-0.39 or simulation based on high (200 mm h-1) and moderate (100 mm h-1) rainfall events over the CSS and AFS areas, respectively. Other potential impacts of charcoal production on watershed hydrology were described. The results presented, together with watershed measurements, when available, are expected to enhance understanding of the hydrological responses of ecosystems to indiscriminate charcoal production and related activities in this region.