Reducing runoff by managing crop location at the catchment level,considering agronomic constraints at farm level

Runoff and erosion cause frequent damage through muddy floods in the loess belt of Northern Europe. One possibility for reducing damage is to lower runoff on agricultural land by spatially alternating different crops at the catchment level. But crop location results from decisions taken at the farm level. This study aimed to assess the existing leeway to modify crop location in the farms of a catchment, in order to reduce runoff at the catchment's outlet. The case study was the Bourville catchment (1086 ha), cultivated by 28 farmers and located in Pays de Caux, France. First, crop location rules in the 14 main farms of the catchment were analysed on the basis of surveys carried out with farmers, distinguishing spatial constraints from temporal ones. These rules made it possible to simulate crop location on each farm territory for the 2001–2002 crop year. Each field of the catchment was classified depending on whether one or several crops could be sown, taking into account both field history and farmer decision rules. Then two extreme scenarios of crop location in the Bourville catchment were built. Runoff simulation at the outlet with the STREAM model showed that runoff could be reduced while sticking to current farmer decision rules in terms of crop location. Depending on rainfall event characteristics, runoff reduction varied between 13·5 per cent and 4·5 per cent. Copyright © 2006 John Wiley & Sons, Ltd.     

Remote-sensing data as an alternative input for the ‘STREAM’ runoff model

Water erosion of cropland constitutes an issue for natural environments along runoff flowpaths due to property damage by soil-laden water and the associated transfer of nutrients and pesticides. In the Pays de Caux region of northwestern France, the silty soils with crusting properties induce a high risk of runoff and erosion. Changes in agricultural practices, land use and landscape patterns appear to have increased the occurrence of erosion and mud flows over the past few decades.A runoff and erosion model called STREAM, applicable to single rainfall events at catchment scale, has been developed to simulate the impacts of land-use modifications. The model takes into account processes that degrade surface states when calculating infiltration rates, as well as agricultural aspects when computing the runoff circulation network. STREAM is based on an expert-system approach that focuses on the dominant processes whilst having only a few input parameters: three of these are used to determine the runoff circulation network, and the other four to calculate infiltration rates. Input nevertheless requires field observations, which restricts application of the model to small catchments.Satellite data covering large areas is considered as an alternative input for such a model, the main objectives being to adapt STREAM accordingly, and to compare the obtained results with field data. In view of previous work involving the extraction and validation of roughness indices using RADARSAT data, this study is based on RADARSAT and LANDSAT TM data collected during the winter of 1998.After adaptation to receive remote-sensing data, the resulting STREAM-TED model requires less input, namely (1) slope and orientation, (2) land-use classification from optical remote-sensing data, (3) roughness indices from radar remote-sensing data, and (4) previous rainfall.Runoff volumes at a gauged catchment outlet (Bourville in Upper Normandy, France) are simulated by four successive versions of the model ranging from the original STREAM to the adapted STREAM-TED. Predictions of the four versions are compared, and performance of the successive simulations is assessed in relation to measured values and according to five statistical indices.Predictions of runoff volume at the catchment outlet using STREAM-TED are similar to those using the original STREAM model, but with a tendency towards overestimation. The final STREAM-TED version is capable of identifying areas sensitive to runoff within a group of catchments and could be used as a planning decision tool in the implementation of conservation practices.     

Simulation of Nitrogen and Phosphorus Losses in Loess Landforms of Northern Iran

Population growth, urban expansion and intensive agriculture and thus increased use of fertilizers aimed at increasing the production capacity cause extensive loss of nutrients such as nitrogen and phosphorus and lead to reduced quality of soil and water. Therefore, identification of nutrients in the soil and their potential are essential. The aim of this study was to evaluate the capability of the SWAT model in simulating runoff, sediment, and nitrogen and phosphorus losses in Tamer catchment. Runoff and sediment measured at Tamar gauging station were used to calibrate and validate the model. Simulated values were generally consistent with the data observed during calibration and validation period (0.6 < R² and 0.5 < NS). In the case of nitrogen loss, the model performed an almost good simulation (0.6 < R² and 0.47 < NS), but phosphorus simulation yielded better results (0.76 < R² and 0.66 < NS). The results showed that cultivated lands had higher loss of nitrogen and phosphorus than other types of land use. Among the various forms of nitrogen and phosphorus, the loss of organic nitrogen and nitrate and soluble phosphorus and mineral phosphorus attached to the sediments showed the greatest sensitivity to the type of land use. Results also showed that the average nutrient loss caused by erosion in this catchment, was 6.99 kg/ha for nitrogen, 0.35 kg/ha for nitrate, 1.3 kg/ha for organic phosphorus, 0.015 kg/ha for soluble phosphorus, and 0.45 kg/ha for mineral phosphorus.     

Feasibility of using FGD gypsum to conserve water and reduce erosion from an agricultural soil in Georgia

Crop production in Georgia and the Southeastern U.S. can be limited by water. Highly-weathered, drought-prone soils are susceptible to runoff and erosion. Rainfall patterns generate runoff producing storms followed by extended periods of drought during the crop growing season. Thus, supplemental irrigation is often needed to sustain profitable crop production. Increased water retention and soil conservation would efficiently improve water use and reduce irrigation amounts/costs and sedimentation, and sustain productive farm land, thus improving producer's profit margin. Soil amendments, such as flue gas desulfurization (FGD) gypsum, have been shown to retain rainfall and/or irrigation water through increased infiltration while decreasing runoff (R) and sediment (E). Objectives were to quantify rainfall partitioning and sediment delivery improvements with surface applied FGD gypsum from an Ultisol managed to conventional till (CT) and to assess the feasibility of using FGD gypsum on agricultural land in southern Georgia. A field study (Faceville loamy sand, Typic Kandiudult) was established (2006, 2007) near Dawson, GA managed to CT, irrigated cotton (Gossypium hirsutum L.). FGD gypsum application rates evaluated were 0, 1.1, 2.2, 4.5, and 9 Mg ha^− 1. Gypsum treatments and simulated rainfall (50 mm h^− 1 for 1 h) were applied to 2-m wide × 3-m long field plots (n = 3). Runoff and E were measured from each 6-m² plot (slope = 1%). FGD gypsum plots averaged 26% more infiltration (INF), 40% less R, 58% less E, 27% lower maximum R rates (R_max), and 2 times lower maximum E rates (E_max) than control plots. Values of INF and water for crop use increased, and R, E, R_max, and E_max decreased as FGD gypsum application rate increased. Values of INF, R, E, R_max, and E_max for 9 Mg ha^− 1 plots were as much as 17% greater, 35% less, 1.9 times less, 35% less, and 1.9 times less than those from other FGD gypsum plots, respectively; and 40% greater, 40% less, 2.2 times less, 52% less, and 2.9 times less than those from control plots, respectively. Applying FGD gypsum to agricultural lands is a cost-effective management practice for producers in Georgia that beneficially impacts natural resource conservation, producer profit margins, and environmental quality. Agriculture in the Southeast provides a viable market for the electric power industry to convert disposal costs of FGD gypsum into a profitable commodity.     

Residue cover and rainfall intensity effects on runoff soil organic carbon losses

Soil cover and rainfall intensity (RI) are recognized to have severe impacts on soil erosion and an interaction exists between them. This study investigates the effect of rainfall intensity (RI) and soil surface cover on losses of sediment and the selective enrichment of soil organic carbon (SOC) in the sediment by surface runoff. A field rainfall simulator was used in the laboratory to produce 90 min rainfall events of three rainfall intensities (65, 85 and 105 mm h^− 1) and four cover percentages (0%, 25%, 50% and 75%) on soil material at 9% slope. A strong negative exponential relation was observed between cover percentage and RI on sediment loss under 85 and 105 mm h^− 1 of rain, while under RI of 65 mm h^− 1, the highest sediment loss was observed under 25% cover. Overall, higher RI and lower cover produced higher sediment and consequently higher nutrient loss, but resulted in a lower SOC enrichment ratio (ER_SOC) in the sediment. The amount of runoff sediment rather than the ER_SOC in the sediment was the determinant factor for the amount of nutrients lost. The values of ER_SOC were high and positively correlated with the ER values of particles smaller than 20 µm (p < 0.01). Although the sediment contained substantially more fine fractions (fine silt and clay, < 20 µm), the original soil and runoff sediment were still of the same texture class, i.e. silt clay loam.     

Accuracy of methods for field assessment of rill and ephemeral gully erosion

To properly assess soil erosion in agricultural areas, it is necessary to determine precisely the volume of ephemeral gullies and rills in the field by using direct measurement procedures. However, little information is available on the accuracy of the different methods used. The main purpose of this paper is to provide information for a suitable assessment of rill and ephemeral gully erosion with such direct measurement methods. To achieve this objective: a) the measurement errors associated to three methods used for field assessment of channel cross sectional areas are explored; b) the influence of the number of cross sections used per unit channel length on the assessment accuracy, is analysed and; c) the effect of the channel size and shape on measurement errors is examined. The three methods considered to determine the cross sectional areas were: the micro-topographic profile meter (1); the detailed measurement of section characteristic lengths with a tape (2); and the measurement of cross section width and depth with a tape (3). Five reaches of different ephemeral gully types 14.0 or 30.0 m long and a set of six 20.4 to 29.4 m long rill reaches were selected. On each gully reach, the cross sectional areas were measured using the three above mentioned methods, with a separation (s) between cross sections of 1 m. For rills, the cross sectional areas were measured with methods 1 and 3, with s = 2 m. Then, the corresponding total erosion volumes were computed. The volume calculated with method 1 with s = 1 m for gullies and s = 2 m for rills was taken as the reference method. For each channel, and for each one of the possible combinations of s and measurement method (m), the relative measurement error and the absolute value of the relative measurement error (E_sm^r and |E_sm^r|), defined with respect to the reference one, was calculated. |E_sm^r| much higher than 10% were obtained very easily, even for small s values and for apparently quasi prismatic channels. Channel size and shape had a great influence on measurement errors. In fact, the selection of the more suitable method for a certain gully shape and size seemed to be much more important than s, at least when s < 10 m. Method 1 always provided the most precise measurements, and its results were the less dependent on s. However, s must be < 5 m to guarantee an error smaller than 10%. Method 2 is not recommended, because it is difficult, time consuming and can lead to large errors. Method 3 seems to be enough for small, wide and shallow gullies, and for small rills, but only if s is shorter than 5 m. Results obtained after the analysis of rill measurement errors were similar to those of gullies. The analysis of E_sm^r and |E_sm^r| when calculating channel volumes using a unique representative cross section highlighted the importance of correctly selecting the adequate cross section. Due to the high error values that this method can entail, it is not considered as advisable whenever accurate erosion measurements are pursued.     

Development and application of a simple hydrologic model simulation for a Brazilian headwater basin

Physically based hydrologic models for watersheds are important tools to support water resources management and predict hydrologic impacts produced by land-use change. Grande River Basin is located in southern Minas Gerais State, and the Grande River is the main tributary of Basin which has 2080 km² draining into the Camargos Hydropower Plant Reservoir (CEMIG — “Minas Gerais State Energy Company”). The objectives of this work were: 1) to create a semi-physically based hydrologic model in semi-distributed to sub-basins approach and based on GIS and Remote Sensing tools and, 2) to simulate the hydrologic responses of the Grande River Basin, thus creating an important tool for management and planning of water resources for region. The hydrologic model is based on the SCS Curve Number (SCS-CN) and MGB/IPH models, and structured into three hydrologic components: estimation of the flow components (quick runoff, hortonian and base flows), propagation into the respective soil reservoirs (surface, sub-surface and shallow saturated zone) and propagation into the channels. Precipitation and discharge data sets were obtained from the Brazilian National Water Agency (ANA). Reference evapotranspiration (ETo) data were obtained from the Brazilian National Meteorological Institute (INMET). In order to estimate actual evapotranspiration, crop coefficient, soil moisture and satellite image interpretation of actual land-use were applied. The long-term hydrologic data series were structured for period between 1990 and 2003. The calibration and validation process was carried out by evaluating the behavior of the Nash–Sutcliffe Coefficient (C_NS), obtained from three different combinations of calibration and validation years. This allowed us to evaluate the model performance to simulate years in which El Niño (EN) and La Niña (LN) events were registered (1997–1998 and 1999–2000, respectively). The combinations of calibration and validation years were: the first 7 years to calibrate and remaining 6 years to validate; the first 9 years to calibrate and remaining 4 years to validate; and first 11 years to calibrate and the last 2 years to validate. The statistical precision showed that the model was able to simulate the hydrologic impacts, including years of EN and LN events, with C_NS scores greater than 0.70 in both situations. The evaluation of the C_NS scores showed small variation in the coefficient as the years of validation decreased. In addition, the model was also able to simulate the hydrologic impacts of land-use change in the Grande River Basin, based on the C_NS scores of 0.80 for different combinations of validation periods. The hydrologic impacts in Grande River Basin produced from grassland area converted to eucalyptus under three specific scenarios were evaluated, which predicted annual runoff mean reductions of up to 17.8%, due to an increase in evapotranspiration rate for the eucalyptus plantation.     

Development and application of a simple hydrologic model simulation for a Brazilian headwater basin

Physically based hydrologic models for watersheds are important tools to support water resources management and predict hydrologic impacts produced by land-use change. Grande River Basin is located in southern Minas Gerais State, and the Grande River is the main tributary of Basin which has 2080 km² draining into the Camargos Hydropower Plant Reservoir (CEMIG — “Minas Gerais State Energy Company”). The objectives of this work were: 1) to create a semi-physically based hydrologic model in semi-distributed to sub-basins approach and based on GIS and Remote Sensing tools and, 2) to simulate the hydrologic responses of the Grande River Basin, thus creating an important tool for management and planning of water resources for region. The hydrologic model is based on the SCS Curve Number (SCS-CN) and MGB/IPH models, and structured into three hydrologic components: estimation of the flow components (quick runoff, hortonian and base flows), propagation into the respective soil reservoirs (surface, sub-surface and shallow saturated zone) and propagation into the channels. Precipitation and discharge data sets were obtained from the Brazilian National Water Agency (ANA). Reference evapotranspiration (ETo) data were obtained from the Brazilian National Meteorological Institute (INMET). In order to estimate actual evapotranspiration, crop coefficient, soil moisture and satellite image interpretation of actual land-use were applied. The long-term hydrologic data series were structured for period between 1990 and 2003. The calibration and validation process was carried out by evaluating the behavior of the Nash–Sutcliffe Coefficient (C_NS), obtained from three different combinations of calibration and validation years. This allowed us to evaluate the model performance to simulate years in which El Niño (EN) and La Niña (LN) events were registered (1997–1998 and 1999–2000, respectively). The combinations of calibration and validation years were: the first 7 years to calibrate and remaining 6 years to validate; the first 9 years to calibrate and remaining 4 years to validate; and first 11 years to calibrate and the last 2 years to validate. The statistical precision showed that the model was able to simulate the hydrologic impacts, including years of EN and LN events, with C_NS scores greater than 0.70 in both situations. The evaluation of the C_NS scores showed small variation in the coefficient as the years of validation decreased. In addition, the model was also able to simulate the hydrologic impacts of land-use change in the Grande River Basin, based on the C_NS scores of 0.80 for different combinations of validation periods. The hydrologic impacts in Grande River Basin produced from grassland area converted to eucalyptus under three specific scenarios were evaluated, which predicted annual runoff mean reductions of up to 17.8%, due to an increase in evapotranspiration rate for the eucalyptus plantation.     

Predicting runoff from semi-arid hillslopes as source areas for water harvesting in the Sierra de Gador, southeast Spain

The effectiveness of water harvesting systems collecting surface runoff form rangeland hillslopes in semi-arid regions is difficult to predict, since the hydrological response at the outlet depends on the heterogeneity of hydrological processes. The lack of continuous runoff pathways, due to the irregular spatial patterns of soil properties and the variety of antecedent soil moisture conditions directly influence runoff generation and control discharge into the water harvesting cisterns. The aim of this paper is to evaluate the effectiveness of semi-arid hillslopes in generating runoff for water harvesting systems. Runoff was estimated by the STREAM expert-based model which was applied to three semi-arid hillslopes (0.4 to 6 ha). On the one hand the STREAM model rules were adapted to the regional conditions i.e. an antecedent precipitation index was adjusted to local soil moisture conditions and the rainfall duration was defined as the total rainfall event quantity and the effective rainfall duration (P_tot/t_eff). On the other hand, the distribution of rock outcrop and vegetation cover along the slopes was used to define homogeneous hydrological units. Final infiltration capacities were attributed to these hydrological units based on values found in the literature. The prediction performances are acceptable for the three water harvesting systems with an RMSE of 13.9 m³. It was shown that the rainfall/runoff model was more sensitive to the duration of the storm than to the antecedent soil moisture conditions. The use of a unique set of hydrological parameters for the three water harvesting systems on representative hillslopes allows the runoff prediction from any rangeland hillslope within the same region. Furthermore, the spatial patterns of soil surface characteristics are crucial for collecting runoff at the outlet of the system. Model runs demonstrated that degradation of vegetation and sealing of very small areas within flow paths can lead to an increase of annual runoff by as much as a factor two.     

Effect of snow amount on runoff,soil loss and suspended sediment during periods of snowmelt in southern West Siberia

Surface runoff, soil loss, suspended sediment concentration (SSC), texture of eroded soils and suspended sediment were determined on slightly eroded chernozems (mouldboard fall-ploughed) during years with different amounts of snow in three areas of southern West Siberia (Predsalairye, Priobye and Kuznetsk hollow). These areas have different geomorphological and climatic characteristics and soils. Observations were made from 1969 to 2007. The soil loss during very low-snow and low-snow years did not exceed 2 t ha^− 1. After winters with normal amounts of snow, the runoff led to slight soil loss (2–5 t ha^− 1). Soil losses in high-snow and very high-snow years varied from slight to severe (4.8–15.8 t ha^− 1) depending on studied area. The main sediment exported during intensive snowmelt and the 1 mm of runoff transported from 35 to 150 kg ha^− 1 of soil material. The removal of soil particles < 0.01 mm (especially clay) prevailed during the initial and final stages of snowmelt. Clay removal by meltwater from the ploughed layer in high-snow and very high-snow years varied from 3300 to 4200 kg ha^− 1 and, in the initial and final stages of snowmelt clay removal, accounted for 1260–1,500 kg ha^− 1. Among the three studied regions, Predsalairye had decreased soil erosion resistance and was the area with the greatest danger of erosion.     

Integral energy of conventional available water, least limiting water range and integral water capacity for better characterization of water availability and soil physical quality

Different approaches have been proposed for quantification of soil water availability for plants but mostly they do not fully describe how water is released from the soil to be absorbed by the plant roots. A new concept of integral energy (E_I) was suggested by Minasny and McBratney (Minasny, B., McBratney, A.B. 2003. Integral energy as a measure of soil-water availability. Plant and Soil 249, 253-262) to quantify the energy required for plants to take up a unit mass of soil water over a defined water content range. This study was conducted to explore the E_I concept in association with other new approaches for soil water availability including the least limiting water range (LLWR) and the integral water capacity (IWC) besides conventional plant available water (PAW). We also examined the relationship between E_I and Dexter's index of soil physical quality (S-value). Twelve agricultural soils were selected from different regions in Hamadan province, western Iran. Soil water retention and penetration resistance, Q, were measured on undisturbed samples taken from the 5-10 cm layer. The PAW, LLWR and IWC were calculated with two matric suctions (h) of 100 and 330 hPa for field capacity (FC), and then the E_I values were calculated for PAW, LLWR and IWC. There were significant differences (P < 0.01) between the E_I values calculated for PAW₁₀₀, PAW₃₃₀, LLWR₁₀₀, LLWR₃₃₀ and IWC. The highest (319.0 J kg⁻¹) and the lowest (160.7 J kg⁻¹) means of E_I were found for the E_I(IWC) and E_I(PAW330), respectively. The E_I values calculated for PAW₁₀₀, LLWR₁₀₀ and LLWR₃₃₀ were 225.6, 177.9 and 254.1 J kg⁻¹, respectively. The mean value of E_I(PAW330) was almost twice as large as the mean of E_I(IWC) showing that IWC is mostly located at lower h values when compared with PAW₃₃₀. Significant relationships were obtained between E_I(IWC) and h at Q = 1.5 MPa, and E_I(LLWR100) or E_I(LLWR330) and h at Q = 2 MPa indicating strong dependency of E_I on soil strength in the dry range. We did not find significant relationships between E_I(PAW100) or E_I(PAW330) and bulk density (ρ_b) or relative ρ_b (ρ_b-rel). However, E_I(LLWR100) or E_I(LLWR330) was negatively and significantly affected by ρ_b and ρ_b-rel. Both E_I(PAW100) and E_I(PAW330) increased with increasing clay content showing that a plant must use more energy to absorb a unit mass of PAW from a clay soil than from a sandy soil. High negative correlations were found between E_I(PAW100) or E_I(PAW330) and the shape parameter (n) of the van Genuchten function showing that soils with steep water retention curves (coarse-textured or well-structured) will have lower E_I(PAW). Negative and significant relations between E_I(PAW100) or E_I(PAW330) and S were obtained showing the possibility of using S to predict the energy that must be used by plants to take up a unit mass of water in the PAW range. Our findings show that E_I can be used as an index of soil physical quality in addition to the PAW, LLWR, IWC and S approaches.     

Modelling of snowmelt erosion and sediment yield in a small low-mountain catchment in Germany

Temporal variability and spatial heterogeneity of surface runoff generation triggers the dynamics of source areas of sediment and sediment-associated nutrient transport. Reliable modelling of hydrological special situations i.e. snowmelt is of high importance for the quality of erosion and sediment yield modelling. Data from the research catchment Schäfertal demonstrate the individuality of snowmelt events in terms of runoff coefficient and delivery ratio. This 1.44 km² low mountain catchment is characterised by a high portion of arable land with a winter grain/winter rape crop rotation. The integrated winter erosion and nutrient load model (IWAN) considers these dynamic aspects by coupling a hydrological model with a sediment load model. Cell size of this raster-based approach is 10 × 10 m². Additionally, snowmelt rill erosion is simulated with a newly developed physically based model that is firstly applied on a catchment scale. A sensitivity analysis of this model system component demonstrates the plausibility of the model approach and the overall robustness of the model system IWAN. The results of the long-term hydrological modelling from 1991 to 2003 are reliable and form the basis for the simulation of six snowmelt events which were observed in the Schäfertal catchment. The estimated total runoff volumes for these events match the observations well. The modelled overland runoff coefficients vary from 0.001 to 0.72. The mean values of cell erosion, which were modelled with one set of parameters for all six events range from 0.0006 to 0.96 t ha^− 1. The total modelled erosion for the events with unfrozen soil and low amount of surface runoff is of a factor 50 below those with partly frozen soil. In addition to these distinctions, the major differences are caused by flow accumulation in shallow depressions in variable parts of the catchment. However, the validation of these results on the single event scale is restricted due to limited spatial data. Total simulated sediment yield at the catchment outlet was as high as 13.84 t which underestimates the observed values, with the exception of one event. Oversimplification of the modelled channel processes may be a reason. The temporal variability and spatial heterogeneity of the surface roughness parameter, which was identified to be sensitive, also causes uncertainty in the parameter estimation. Despite these findings, the model system IWAN was applied successfully on the catchment scale and the simulated results are reliable.     

Spatial and temporal variability in CH4 and N2O fluxes from a Scottish ombrotrophic peatland: Implications for modelling and up-scaling

Peatlands typically exhibit significant spatial heterogeneity which can lead to large uncertainties when catchment scale greenhouse gas fluxes are extrapolated from chamber measurements (generally <1 m²). Here we examined the underlying environmental and vegetation characteristics which led to within-site variability in both CH₄ and N₂O emissions and the importance of such variability in up-scaling. We also consider within-site variation in the controls of temporal dynamics. Net annual emissions (and coefficients of variation) for CH₄ and N₂O were 1.06 kg ha⁻¹ y⁻¹ (300%) and 0.02 kg ha⁻¹ y⁻¹ (410%), respectively. The riparian zone was a significant CH₄ hotspot contributing ∼12% of the total catchment emissions whilst covering only ∼0.5% of the catchment area. In contrast to many other studies we found smaller CH₄ emissions and greater uptake in chambers containing either sedges or rushes. We also found clear differences in the drivers of temporal CH₄ dynamics across the site, e.g. water table was important only in chambers which did not contain aerenchymous plants. We suggest that depending on the heterogeneity of the site, flux models could be improved by incorporating a number of spatially distinct sub-models, rather than a single model parameterized using whole-catchment averages.     

Harvest intensity of aromatic shrubs vs. soil erosion: An equilibrium for sustainable agriculture (SE Spain)

In the Mediterranean Europe, where rainfall is scarce and irregular but often of high intensity, wild shrubs protect the soil against erosivity of raindrops. Moreover, some of these plants are the economic income for local farmers. Particularly in SE Spain, soil erosion is a core factor in environmental degradation attributed primarily to the cultivation practices and human pressure on the land. Over a four-year period, soil erosion and runoff were monitored in erosion plots on a mountainside, comparing four harvest intensities of four aromatic shrubs (Lavandula lanata L., Santolina rosmarinifolia L. Origanum bastetanum, and Salvia lavandulifolia V.): 0% (HI-0), 25% (HI-25), 50% (HI-50), and 75% (HI-75). Also, the fresh biomass and essential-oil content were quantified in each treatment. The erosion plots were located in Lanjarón (Granada, SE Spain) on the southern flank of the Sierra Nevada Mountains, on a 20% slope, and of 96 m² in area. The average soil loss for HI-0, HI-25, HI-50, and HI-75 during the study period was 144.6, 187.2, 256.0, and 356.0 kg ha^− 1, respectively, and runoff 2.6, 3.2, 3.4, and 4.7 mm, respectively. The lowest average soil erosion and runoff rates for the study period were recorded with plant cover of S. lavandulifolia V. 67.6 kg ha^− 1 and 1.3 mm, respectively. Since no significant differences were found between HI-25 and HI-50 for soil erosion and runoff, and harvest and distillation of wild-aromatic plants currently persists as an important economic activity in mountainous areas of the study zone, we recommend a rational harvest (HI-50), leaving the 50% of the plant biomass in the field (especially for sage and lavender) to avoid the soil degradation. In this sense, the harvest of 50% of fresh herb of sage, santolina, lavender, and oregano produced reasonable essential-oil yield of 12.7, 14.0, 19.7, and 18.3 L ha^− 1, respectively. The inappropriate harvest of aromatic plants and the intensity farming systems of mountain areas endanger land conservation, and there is an urgent need to implement appropriate land management which has a large-scale perspective but acts at the local level.     

Soil biochemical activity and phosphorus transformations and losses from acidified forest soils

Three Bohemian Forest catchments, Plešné, ?erné and ?ertovo, were studied. These catchments have similar climatic conditions, relief and vegetation, but differ in their bedrock composition. The granitic bedrock in the Plešné catchment was more susceptible to phosphorus (P) leaching under acid conditions than was the mica schist bedrock in the other catchments. The goal of this study was to determine if higher P leaching from the Plešné catchment was associated with differences in microbial P transformations and enzymatic P hydrolysis. Phosphorus and nitrogen contents in soil microbial biomass (P_MB, N_MB; chloroform fumigation), C mineralisation rate (C_min; CO₂ production by GC) and phosphatase activity (MUF-phosphate), were measured in three successive years. Phosphatase activity, P_MB, and C_min were used to characterise the enzymatic hydrolysis of organic P, microbial P accumulation, and microbial mineralisation rates of organic compounds, respectively. Soil chemical properties were characterised by C, N and P content, pH, and by oxalate-extractable P, Fe and Al. Spatial variability in N_MB, P_MB, C_min and phosphatase activity within the catchment was higher (coefficient of variation, CV<50%) than their temporal variability (CV<30%). Multivariate analysis revealed a significant soil layer effect but not that of catchment. When soil layers were evaluated separately, a difference between the Plešné and ?erné or ?ertovo catchments was found in litter and mineral layers, even though the variability within one catchment was high. Within soil profile, phosphatase activity was positively correlated with C_tot, N_MB and C_min (r²=0.89-0.92) being very correlated with P_MB (r²=0.99). Phosphatase activity was higher in the litter (14.0 nmol g⁻¹ h⁻¹) and humus (8.65 nmol g⁻¹ h⁻¹) layers of Plešné than in the same layers of the ?erné (9.65 and 6.40 nmol g⁻¹ h⁻¹) and ?ertovo (12.8 and 6.0 nmol g⁻¹ h⁻¹) soils. Similarly, P_MB in the litter and humus layers of Plešné soil (161 and 93 μg g⁻¹) was higher than P_MB of the same layers of the ?erné (120 and 66 μg g⁻¹) and ?ertovo (148 and 89 μg g⁻¹) soils. High MUFP hydrolysis rate: C_min molar ratio (0.16-1.17 M of P per 1 M of respired C) indicated that potential enzymatic P hydrolysis exceeded estimated microbial P demand (0.034 M of P per 1 M of respired C) in all catchments. The results suggest that higher microbial P transformations and enzymatic P hydrolysis could contribute to enhanced P leaching from the Plešné catchment, which could be enhanced by the lower Fe content in the soil of this catchment as compared to the ?erné and ?ertovo catchments.     

Surface runoff phosphorus (P) loss in relation to phosphatase activity and soil P fractions in Florida sandy soils under citrus production

Phosphorus losses by surface runoff from agricultural lands have been of public concern due to increasing P contamination to surface waters. Five representative commercial citrus groves (C1-C5) located in South Florida were studied to evaluate the relationships between P fractions in soils, surface runoff P, and soil phosphatase activity. A modified Hedley P sequential fractionation procedure was employed to fractionate soil P. Soil P consisted of mainly organically- and Ca/Mg-bound P fractions. The organically-bound P (biological P, sum of organic P in the water, NaHCO₃ and NaOH extracts) was dominant in the acidic sandy soils from the C2 and C3 sites (18% and 24% of total soil P), whereas the Ca/Mg-bound P (HCl-extractable P) accounted for 45-60% of soil total P in the neutral and alkaline soils (C1, C4 and C5 soils). Plant-available P (sum of water and NaHCO₃ extractable P fractions) ranged from 27 to 61 mg P kg⁻¹ and decreased in the order of C3>C4>C1>C2>C5. The mean total P concentrations (TP) in surface runoff water samples ranged from 0.51 to 2.64 mg L⁻¹. Total P, total dissolved P (TDP), and PO₄³⁻-P in surface runoff were significantly correlated with soil biological P and plant-available P forms (p<0.01), suggesting that surface runoff P was directly derived from soil available P pools, including H₂O- and NaHCO₃- extractable inorganic P, water-soluble organic P, and NaHCO₃- and NaOH-extractable organic P fractions, which are readily mineralized by soil microorganisms and/or enzyme mediated processes. Soil neutral (55-190 mg phenol kg⁻¹ 3 h⁻¹) and natural (measured at soil pH) phosphatase activities (77-295 mg phenol kg⁻¹ 3 h⁻¹) were related to TP, TDP, and PO₄³⁻-P in surface runoff, and plant-available P and biological P forms in soils. These results indicate that there is a potential relationship between soil P availability and phosphatase activities, relating to P loss by surface runoff. Therefore, the neutral and natural phosphatase activities, especially the natural phosphatase activity, may serve as an index of surface runoff P loss potential and soil P availability.     

Using Cs measurements to validate the application of the AGNPS and ANSWERS erosion and sediment yield models in two small Devon catchments

Distributed erosion and sediment yield models are being increasingly used for predicting soil erosion and sediment yields in agricultural catchments. In most applications, validation of such models has commonly been restricted to comparison of the predicted and measured sediment output from a catchment, because spatially distributed information on rates and patterns of soil redistribution within the catchment has been lacking. However, such spatially distributed data are needed for rigorous model testing, in order to validate the internal functioning of a model and its applicability at different spatial scales. The study reported in this paper uses two approaches to test the performance of the agricultural non-point source pollution (AGNPS) and areal non-point source watershed environmental response simulation (ANSWERS) erosion and sediment yield models in two small catchments in Devon, UK. These involve, firstly, comparison of observed and predicted runoff and sediment output data for individual storm events monitored at the basin outlets and, secondly, information on the spatial pattern of soil redistribution within the catchments derived from ¹³⁷Cs measurements. The results obtained indicate that catchment outputs simulated by both models are reasonably consistent with the recorded values, although the AGNPS model appears to provide closer agreement between observed and predicted values. However, the spatial patterns of soil redistribution and the sediment delivery ratios predicted for the two catchments by the AGNPS and ANSWERS models differ significantly. Comparison of the catchment sediment delivery ratios and the pattern of soil redistribution in individual fields predicted by the models with equivalent information derived from ¹³⁷Cs measurements indicates that the AGNPS model provides more meaningful predictions of erosion and sediment yield under UK conditions than the ANSWERS model and emphasises the importance of using information on both catchment output and sediment redistribution within the catchment for model validation.     

Effects of 30-year-old Aleppo pine plantations on runoff,soil erosion,and plant diversity in a semi-arid landscape in south eastern Spain

Forest management policies in Mediterranean areas have traditionally encouraged land cover changes, with the establishment of tree cover (Aleppo pine) in natural or degraded ecosystems for soil conservation purposes: to reduce soil erosion and to increase the vegetation structure. In order to evaluate the usefulness of these management policies on reduced erosion in semi-arid landscapes, we compared 5 vegetation cover types (bare soil, dry grassland, shrublands, afforested dry grasslands and afforested thorn shrublands), monitored in 15 hydrological plots (8 × 2 m), in the Ventós catchment (Alicante, SE Spain), over 4 years (1996 to 1999). Each cover type represented a different dominant patch of the vegetation mosaic on the north-facing slopes of this catchment. The results showed that runoff coefficients of vegetated plots were less than 1% of the precipitation volume; whereas runoff in denuded areas was nearly 4%. Soil losses in vegetation plots averaged 0.04 Mg ha^− 1 year^− 1 and increased 40-fold in open-land plots. The evaluation of these forest management policies, in contrast with the natural vegetation communities, suggests that: (1) thorn shrublands and dry grassland communities with vegetation cover could control runoff and sediment yield as effectively as Aleppo pine afforestation in these communities, and (2) afforestation with a pine stratum improved the stand's vertical structure resulting in pluri-stratified communities, but reduced the species richness and plant diversity in the understorey of the plantations.     

Effectiveness of geotextiles in reducing runoff and soil loss: A synthesis

Despite geotextiles having potential for soil conservation, limited scientific data are available to assess the effects of geotextiles in reducing runoff and water erosion. Hence, the objective of this review is to analyse the effects of plot length (L) and other possible affecting factors [cover percentage (C, %), slope gradient (S), rainfall duration (D), rainfall intensity (I), sand, silt and clay contents, soil organic matter (SOM) content and geotextile type (natural or synthetic)] on the effectiveness of geotextiles in reducing soil and water loss, based on reported experimental data. From linear regressions, C (%) and soil sand, silt and clay contents are found to be the most important variables in reducing SLR (ratio of soil loss in bare plots to that in geotextile treated plots) for splash, C (%) for interrill and D (min) for rill and interrill erosion processes, respectively. Soil clay and silt contents and D are key variables in decreasing RR (ratio of runoff from bare plots to that from geotextile treated plots) for interrill, and clay content for rill and interrill erosion processes, respectively. The linear relationship between mean b-value (geotextile effectiveness factor in reducing soil loss) and L of all studies was not significant (P > 0.05). The same is true for the relationship between L and SLR, and L and RR. However, when L is added to an equation as an interaction term with C (%), geotextile cover is significantly (P < 0.05) more effective in reducing SLR on shorter plots than longer ones for both interrill and rill and interrill erosion processes. Buffer strip plots (area coverage ∼ 10%) with Borassus and Buriti mats have the highest b-values.