Prediction of concentrated flow width in ephemeral gully channels

Empirical prediction equations of the form W = aQ^b have been reported for rills and rivers, but not for ephemeral gullies. In this study six experimental data sets are used to establish a relationship between channel width (W, m) and flow discharge (Q, m³ s^?1) for ephemeral gullies formed on cropland. The resulting regression equation (W = 2·51 Q^0·412; R² = 0·72; n = 67) predicts observed channel width reasonably well. Owing to logistic limitations related to the respective experimental set ups, only relatively small runoff discharges (i.e. Q < 0·02 m³s^?1) were covered. Using field data, where measured ephemeral gully channel width was attributed to a calculated peak runoff discharge on sealed cropland, the application field of the regression equation was extended towards larger discharges (i.e. 5 × 10^?4m³s^?1 < Q < 0·1 m³s^?1). Comparing W–Q relationships for concentrated flow channels revealed that the discharge exponent (b) varies from 0·3 for rills over 0·4 for gullies to 0·5 for rivers. This shift in b may be the result of: (i) differences in flow shear stress distribution over the wetted perimeter between rills, gullies and rivers, (ii) a decrease in probability of a channel formed in soil material with uniform erosion resistance from rills over gullies to rivers and (iii) a decrease in average surface slope from rills over gullies to rivers. The proposed W–Q equation for ephemeral gullies is valid for (sealed) cropland with no significant change in erosion resistance with depth. Two examples illustrate limitations of the W–Q approach. In a first example, vertical erosion is hindered by a frozen subsoil. The second example relates to a typical summer situation where the soil moisture profile of an agricultural field makes the top 0·02 m five times more erodible than the underlying soil material. For both cases observed W values are larger than those predicted by the established channel width equation for concentrated flow on cropland. For the frozen soils the equation W = 3·17 Q^0·368 (R² = 0·78; n = 617) was established, but for the summer soils no equation could be established. Copyright © 2002 John Wiley & Sons, Ltd.     

Hydraulic characteristics of a small Coastal Plain stream of the southeastern United States: effects of hydrology and season

Various physical and biological properties affect solute transport patterns in streams. We measured hydraulic characteristics of Payne Creek, a low‐gradient upper Coastal Plain stream, using tracer experiments and parameter estimation with OTIS‐P (one‐dimensional transport with inflow and storage with parameter optimization). The primary objective of this study was to estimate the effects of varying discharge, season, and litter accumulation on hydraulic parameters. Channel area A ranged from 0·081 to 0·371 m² and transient storage area A_s ranged from 0·027 to 0·111 m². Dispersion D ranged from 1·5 to 11·1 m² min⁻¹ and exchange coefficient α ranged from 0·009 to 0·038 min⁻¹. Channel area and dispersion were positively correlated to discharge Q, whereas storage area and exchange coefficient were not. Relative storage size A_s/A ranged from 0·17 to 0·59, and was higher during fall than other seasons under a similar Q. The fraction of median travel time due to transient storage ranged from 8·8 to 34·5% and was significantly correlated with Q through a negative power function. Both metrics indicated that transient storage was a significant component affecting solute transport in Payne Creek, especially during the fall. Comparison between the measured channel area A_c and A suggested that surface storage was dominant in Payne Creek. During fall, accumulation of leaf litter resulted in larger A and A_s and lower velocity and D than during other seasons with similar discharge. Seasonal changes in discharge and organic matter accumulation, and dynamic channel morphology affected the magnitude of transient storage and overall hydraulic characteristics of Payne Creek. Copyright © 2005 John Wiley & Sons, Ltd.     

Predicting unit plot soil loss in Sicily,south Italy

Predicting soil loss is necessary to establish soil conservation measures. Variability of soil and hydrological parameters complicates mathematical simulation of soil erosion processes. Methods for predicting unit plot soil loss in Sicily were developed by using 5 years of data from replicated plots. At first, the variability of the soil water content, runoff, and unit plot soil loss values collected at fixed dates or after an erosive event was investigated. The applicability of the Universal Soil Loss Equation (USLE) was then tested. Finally, a method to predict event soil loss was developed. Measurement variability decreased as the mean increased above a threshold value but it was low also for low values of the measured variable. The mean soil loss predicted by the USLE was lower than the measured value by 48%. The annual values of the soil erodibility factor varied by seven times whereas the mean monthly values varied between 1% and 244% of the mean annual value. The event unit plot soil loss was directly proportional to an erosivity index equal to , being Q_RR_e the runoff ratio times the single storm erosion index. It was concluded that a relatively low number of replicates of the variable of interest may be collected to estimate the mean for both high and particularly low values of the variable. The USLE with the mean annual soil erodibility factor may be applied to estimate the order of magnitude of the mean soil loss but it is not usable to estimate soil loss at shorter temporal scales. The relationship for estimating the event soil loss is a modified version of the USLE‐M, given that it includes an exponent for the Q_RR_e term. Copyright © 2007 John Wiley & Sons, Ltd.     

Plot‐scale measurement of soil erosion at the experimental area of Sparacia (southern Italy)

Obtaining good quality soil loss data from plots requires knowledge of the factors that affect natural and measurement data variability and of the erosion processes that occur on plots of different sizes. Data variability was investigated in southern Italy by collecting runoff and soil loss from four universal soil‐loss equation (USLE) plots of 176 m², 20 ‘large’ microplots (0·16 m²) and 40 ‘small’ microplots (0·04 m²). For the four most erosive events (event erosivity index, R_e ≥ 139 MJ mm ha^?1 h^?1), mean soil loss from the USLE plots was significantly correlated with R_e. Variability of soil loss measurements from microplots was five to ten times greater than that of runoff measurements. Doubling the linear size of the microplots reduced mean runoff and soil loss measurements by a factor of 2·6–2·8 and increased data variability. Using sieved soil instead of natural soil increased runoff and soil loss by a factor of 1·3–1·5. Interrill erosion was a minor part (0·1–7·1%) of rill plus interrill erosion. The developed analysis showed that the USLE scheme was usable to predict mean soil loss at plot scale in Mediterranean areas. A microplot of 0·04 m² could be used in practice to obtain field measurements of interrill soil erodibility in areas having steep slopes. Copyright © 2003 John Wiley & Sons, Ltd.     

Average sand particle trajectory examined by the Raindrop Detachment and Wind‐driven Transport (RD‐WDT) process

Recent studies of soil loss by the integrated action of raindrop impact and wind transport have demonstrated the significance of this mechanism. This paper presents data obtained during wind‐tunnel experiments examining the ‘Raindrop Detachment and Wind‐driven Transport’ (RD‐WDT) process to investigate average sand particle trajectory and the spatial extent at which the process operates. In the experimental design, at the same time as the horizontal wind velocities of 6·4, 10, and 12 m s^–1 passed through the tunnel, rainfall was simulated falling on very well sorted dune sand. The aspect and slope of the sand bed was varied to reproduce both windward (Ww) and leeward (Lw) slopes of 4º and 9º with respect to the prevailing wind direction. The average sand particle trajectories by the RD‐WDT process ( ) were estimated by a mass‐distribution function, which was integrated over a 7‐m uniform slope segment. The results showed that depended statistically upon the wind shear velocity (u*), and the effect of the slope gradient (θ) was insignificant on . This was different from that of the windless rain process ( ), ‘Raindrop Detachment and Splash‐driven Transport’ (RD‐ST), the spatial range of which relies strongly on θ. Additionally, was approximately 2·27 ± 2·2 times greater than the average path of a typical saltating sand particle of the rainless wind ( ), ‘Wind Erosion Saltation Transport’ (WE‐ST). Copyright © 2009 John Wiley & Sons, Ltd.     

FAST HANKEL TRANSFORMS*

Inspired by the linear filter method introduced by D. P. Ghosh in 1970 we have developed a general theory for numerical evaluation of integrals of the Hankel type: Replacing the usual sine interpolating function by sinsh (x) =a· sin (ρx)/sinh (aρx), where the smoothness parameter a is chosen to be “small”, we obtain explicit series expansions for the sinsh-response or filter function H*. If the input function f(λ exp (iω)) is known to be analytic in the region o < λ < ∞, |ω|≤ω₀ of the complex plane, we can show that the absolute error on the output function is less than (K(ω₀)/r) · exp (?ρω₀/Δ), Δ being the logarthmic sampling distance. Due to the explicit expansions of H* the tails of the infinite summation ((m?n)Δ) can be handled analytically. Since the only restriction on the order is ν > ? 1, the Fourier transform is a special case of the theory, ν=± 1/2 giving the sine- and cosine transform, respectively. In theoretical model calculations the present method is considerably more efficient than the Fast Fourier Transform (FFT).     

The effect of conservation tillage on runoff erosivity and soil erodibility during concentrated flow

Crop residues in conservation tillage systems are known to cause both a reduction in the erosive runoff power and an increase in the topsoil erosion resistance. In this study, the relative importance of both mechanisms in reducing soil loss by concentrated flow erosion is examined. Therefore, a method to calculate the effective flow shear stress responsible for soil detachment in the presence of a residue cover is applied. The determination of effective flow shear stress is based on the recalculation of the hydraulic radius for residue treatments. The method was tested in a laboratory flume by comparing soil detachment rates of identical pairs of soil samples that only differ in the presence or absence of crop residues. This shear stress partitioning approach and a soil detachment correction were then applied to a dataset of soil detachment measurements on undisturbed topsoil samples from a no‐till field plot on a loess‐derived soil, sampled during one growing season. Results indicate that only a small fraction (10% on average) of the difference in soil detachment rate between conventional and conservation tillage can be attributed to the dissipation of shear forces on the residues. The remaining decrease in soil detachment during concentrated runoff after a two‐year application of conservation tillage can be explained by the increased dry bulk density and root and crop residue content in the topsoil that reduces soil erodibility. After correcting for the presence of residues, the temporal variability in soil detachment rates (D_r) during concentrated flow for a given flow shear stress (τ) for both treatments can be predicted fairly well (R² = 0·87) from dry soil bulk density (DBD, representing consolidation effects), soil moisture content (SMC, representing antecedent rainfall conditions), the dry mass of organic material (OM, representing root growth and residue decomposition) and saturated soil shear strength σ_{s, sat} using an equation of the form: This study is the first to show that the effect of conservation tillage on soil detachment rates is a result of soil property modifications affecting soil erodibility, rather than a result of the surface residue decreasing flow erosivity. Copyright © 2007 John Wiley & Sons, Ltd.     

Interrill erosion on cultivated Greek soils: modelling sediment delivery

For interrill erosion, raindrop‐induced detachment and transport of sediment by rainfall‐disturbed sheet flow are the predominant processes, while detachment by sheet flow and transport by raindrop impact are negligible. In general, interrill subprocesses are inter‐actively affected by rainfall, soil and surface properties. The objective of this work was to study the relationships among interrill runoff and sediment loss and some selected para‐meters, for cultivated soils in central Greece, and also the development of a formula for predicting single storm sediment delivery. Runoff and soil loss measurement field experiments have been conducted for a 3·5‐year period, under natural storms. The soils studied were developed on Tertiary calcareous materials and Quaternary alluvial deposits and were textured from sandy loam to clay. The second group of soils showed greater susceptibility to sealing and erosion than the first group. Single storm sediment loss was mainly affected by rain and runoff erosivity, being significantly correlated with rain kinetic energy (r = 0·64***), its maximum 30‐minute intensity (r = 0·64***) and runoff amount (r = 0·56***). Runoff had the greatest correlation with rain kinetic energy (r = 0·64***). A complementary effect on soil loss was detected between rain kinetic energy and its maximum 30‐minute intensity. The same was true for rain kinetic energy and topsoil aggregate instability, on surface seal formation and thus on infiltration characteristics and overland flow rate. Empirical analysis showed that the following formula can be used for the successful prediction of sediment delivery (D_i): D_i = 0·638βEI₃₀tan(θ) (R² = 0·893***), where β is a topsoil aggregate instability index, E the rain kinetic energy, I₃₀ the maximum 30‐minute rain intensity and θ the slope angle. It describes soil erodibility using a topsoil aggregate instability index, which can be determined easily by a simple laboratory technique, and runoff through the product of this index and rain kinetic energy. Copyright © 2006 John Wiley & Sons, Ltd.     

Rotatory vibration of a thick spherical shell of isotropic non-homogeneous elastic material

Summary Rotatory vibrations of a thick spherical shell of isotropic non-homogeneous material with rigidity and density given by (i) = ₀ r ^-2 withQ =Q ₀ r ^-2 e ^2mr and (ii) = ₀ r ^m with =Q ₀ r ⁿ have been discussed and the frequency equation is derived with numerical enumeration of frequency in each case.     

A one‐dimensional model for simulating armouring and erosion on hillslopes: 2. Long term erosion and armouring predictions for two contrasting mine spoils

This paper investigates the dynamics of soil armouring as a result of fluvial erosion for a non‐cohesive sandy gravel spoil from the Ranger Mine, Australia, and a cohesive silt loam spoil from the Northparkes Mine, Australia, using a model for hillslope soil armouring. These long term predictions concentrate on the temporal and spatial changes of the spoil grading and erosion over 100–200 years for the flat cap regions (1–2%) and steep batter edges (10–30%) typically encountered on waste rock dumps. The existence of a significant rock fragment fraction in the Ranger spoil means that it armours readily, while Northparkes does not. For Ranger the waste rock showed reductions in (1) cumulative erosion of up to 81% from that obtained by extrapolating the initial erosion rate out 100 years and (2) the erosion/year by more than 10‐fold. For Northparkes reductions were less marked, with the maximum reduction in erosion/year being 37% after 200 years. For Ranger the reductions were greatest and fastest for intermediate gradient hillslopes. For the steepest hillslopes the armouring decreased because the flow shear stresses were large enough to mobilize all material in the armour layer. Model uncertainty was assessed with probabilistic confidence limits demonstrating that these erodibility reductions were statistically significant. A commonly used hillslope erosion model (sediment flux = β₁ discharge m₁ slope n₁) was fitted to these predictions. The erodibility, β₁, and m₁ decreased with time, which was consistent with our physical intuition about armouring. At Ranger the parameter m₁ asymptoted to 1·5–1·6 while at Northparkes it asymptoted to 1·2–1·3. At Ranger transient spatial trends in armouring led to a short term (50–200 years in the future) reduction in n₁, to below zero under certain circumstances, recovering to an asymptote of about 0·5–1. At Northparkes n₁ asymptoted to about 0·6, with no negative transients predicted. The m₁ and n₁ parameters predicted for Ranger were shown to be consistent with field data from a 10‐year‐old armoured hillslope and consistent with published relationships between erodibility and rock content for natural hillslopes. Copyright © 2007 John Wiley & Sons, Ltd.     

Runoff dependent erosivity and slope length factors suitable for modelling annual erosion using the Universal Soil Loss Equation

Despite revisions and refinements, the Revised Universal Soil Loss Equation (RUSLE), which is the revised version of the Universal Soil Loss Equation (USLE), over predicts small annual soil losses and under predicts large annual soil losses. To some large extent, this results from the equation over estimating small event soil losses and under estimating large event soil losses. Replacing the USLE/RUSLE event erosivity index (EI₃₀) by the product of EI₃₀ and the runoff ratio (Q_R) significantly reduces the errors in estimating event erosion when runoff is measured, but the USLE‐M, the USLE variant that uses the Q_REI₃₀ index, requires crop and support practice factors that differ from those used in the RUSLE. The theory which enables the Q_REI₃₀ index to be used in association with the RUSLE crop and support practice factors is presented. In addition, the USLE/RUSLE approach was developed for conditions where runoff is produced uniformly over a hill slope. A runoff dependent slope length factor that takes account of runoff variations over a hill slope is presented and demonstrated for the situation where runoff from a low runoff producing area passes onto an area where runoff is produced more readily. Copyright © 2007 John Wiley & Sons, Ltd.     

Applying the QREI30 index within the USLE modelling environment

The USLE/RUSLE model was designed to predict long‐term (~20 years) average annual soil loss by accounting for the effects of climate, soil, topography and crops. The USLE/RUSLE model operates mathematically in two steps. The first step involves the prediction of soil loss from the ‘unit’ plot, a bare fallow area 22.1 m long on a 9% slope gradient with cultivation up and down the slope. Appropriate values of the factors accounting for slope length, gradient, crops and crop management and soil conservation practice are then used to adjust that soil loss to predict soil loss from areas that have conditions that are different from the unit plot. Replacing EI₃₀, the USLE/RUSLE event erosivity index, by the product of the runoff ratio (Q_R) and EI₃₀, can enhance the capacity of the model to predict short‐term soil loss from the unit plot if appropriate data on runoff is available. Replacing the EI₃₀ index by another index has consequences on other factors in the model. The USLE/RUSLE soil erodibility factor cannot be used when the erosivity factor is based on Q_REI₃₀. Also, the USLE/RUSLE factors for slope length, slope gradient crops and crop management, and soil conservation practice cannot be used when runoff from other than the unit plot is used to calculate Q_R. Here, equations are provided to convert the USLE/RUSLE factors to values suitable for use when the erosivity factor is based on the Q_REI₃₀ index under these circumstances. At some geographic locations, non linear relationships exist between soil loss from bare fallow areas and the Q_REI₃₀ index. The effect of this on the slope length factor associated with the Q_REI₃₀ index is demonstrated using data from runoff and soil loss plots located at the Sparacia site, Sicily. Copyright © 2012 John Wiley & Sons, Ltd.     

Sediment delivery from hillslopes and the Universal Soil Loss Equation: some perceptions and misconceptions

The Universal Soil Loss Equation (USLE) or the revised USLE (RUSLE) are often used together with sediment delivery ratios in order to predict sediment delivery from hillslopes. In using sediment delivery ratios for this purpose, it is assumed that the sediment delivery ratio for a given hillslope does not vary with the amount of erosion occurring in the upslope area. This assumption is false. There is a perception that hillslope erosion is calculated on the basis that hillslopes are, in effect, simply divided into 22·1 m long segments. This perception fails to recognize the fact the inclusion of the 22·1 m length in the calculation has no physical significance but simply produces a value of 1·0 for the slope length factor when slopes have a length equal to that of the unit plot. There is a perception that the slope length factor is inappropriate because not all the dislodged sediment is discharged. This perception fails to recognize that the USLE and the RUSLE actually predict sediment yield from planar surfaces, not the total amount of soil material dislocated and removed some distance by erosion within an area. The application of the USLE/RUSLE to hillslopes also needs to take into account the fact that runoff may not be generated uniformly over that hillslope. This can be achieved by an equation for the slope length factor that takes account of spatial variations in upslope runoff on soil loss from a segment or grid cell. Several alternatives to the USLE event erosivity index have been proposed in order to predict event erosion better than can be achieved using the EI₃₀ index. Most ignore the consequences of changing the event erosivity index on the values for the soil, crop and soil conservation protection factors because there is a misconception that these factors are independent of one another. Copyright © 2007 John Wiley & Sons, Ltd.     

Using palm‐mat geotextiles on an arable soil for water erosion control in the UK

To date, most studies of the effectiveness of geotextiles on soil erosion rates and processes have been conducted in laboratory experiments for less than 1 h. Hence, at Hilton (52°33′ N, 2°19′ W), UK, the effectiveness of employing palm‐mat geotextiles for soil erosion control under field conditions on arable loamy sands was investigated. Geotextile‐mats constructed from Borassus aethiopum (Borassus palm of West Africa) and Mauritia flexuosa (Buriti palm of South America) leaves are termed Borassus mats and Buriti mats, respectively. Duplicate runoff plots (10 m × 1 m on a 15° slope) had five treatments (bare, permanent grass, Borassus total plot cover, Borassus buffer strip and Buriti buffer strip). Borassus covered plots had about 72% ground cover and to differentiate between this treatment and Borassus buffer strips, the former treatment is termed Borassus completely‐covered. Runoff and eroded soil were collected from each bounded plot in a concrete gutter, leading to a receptacle. Results from 08/01/2007–23/01/2009 (total precipitation = 1776·5 mm; n = 53 time intervals) show that using Borassus buffer strips (area coverage ～10%) on bare soil decreased runoff volume by about 71% (P > 0·05) and soil erosion by 92% (P < 0·001). Bare plots had nearly 29·1 L m^?2 runoff and 2·36 kg m^?2 soil erosion during that period. Borassus buffer strip, Buriti buffer strip and Borassus completely‐covered plots had similar effects in decreasing runoff volume and soil erosion. Runoff volumes largely explain the variability in soil erosion rates. Although buffer strips of Borassus mats were as effective as whole plot cover of the same mats, the longevity of Borassus mats was nearly twice that of Buriti mats. Thus, use of Borassus mats as buffer strips on bare plots is highly effective for soil erosion control. The mechanisms explaining the effectiveness of buffer strips require further studies under varied pedo‐climatic conditions. Copyright © 2011 John Wiley & Sons, Ltd.     

Concentrated flow erosion rates reduced through biological geotextiles

Soil erosion by concentrated flow can cause serious environmental damage. Erosion‐control geotextiles have considerable potential for reducing concentrated flow erosion. However, limited data are available on the erosion‐reducing potential of geotextiles. In this study, the effectiveness of three biological geotextiles in reducing soil losses during concentrated flow is investigated. Hereto, runoff was simulated in a concentrated flow flume, filled with an erodible sandy loam on three slope gradients (13·5, 27·0 and 41·5%). Treatments included three biological geotextiles (borassus, buriti and bamboo) and one bare soil surface. Darcy–Weisbach friction coefficients ranged from 0·01 to 2·84. The highest values are observed for borassus covered soil surfaces, followed by buriti, bamboo and bare soil, respectively. The friction coefficients are linearly correlated with geotextile thickness. For the specific experimental conditions of this study, borassus geotextiles reduced soil detachment rate on average to 56%, buriti geotextiles to 59% and bamboo geotextiles to 66% of the soil detachment rate for bare soil surfaces. Total flow shear stress was the hydraulic parameter best predicting soil detachment rate for bare and geotextile covered surfaces (R² = 0·75–0·84, p <0·001, n = 12–15). The highest resistance against soil detachment was observed for the borassus covered soil surfaces, followed by buriti, bamboo and bare soil surfaces, respectively. Overall, biological geotextiles are less effective in controlling concentrated flow erosion compared with interrill erosion. Copyright © 2009 John Wiley & Sons, Ltd.     

Soil erosion by surface water flow on a stony,semiarid hillslope

Soil erosion on hillslopes occurs by processes of soil splash from raindrop impacts and sediment entrainment by surface water flows. This study investigates the process of soil erosion by surface water flow on a stony soil in a semiarid environment. A field experimental method was developed whereby erosion by concentrated flow could be measured in predefined flow areas without disturbing the soil surface. The method allowed for measurements in this study of flow erosion at a much wider range of slopes (2·6 to 30·1 per cent) and unit discharge rates (0·0007 to 0·007 m2 s−1) than have been previously feasible. Flow velocities were correlated to discharge and hydraulic radius, but not to slope. The lack of correlation between velocity and slope might have been due to the greater rock cover on the steeper slopes which caused the surface to be hydraulically rougher and thus counteract the expected effect of slope on flow velocity. The detachment data illustrated limitations in applying a linear hydraulic shear stress model over the entire range of the data collected. Flow detachment rates were better correlated to a power function of either shear stress (r2 = 0·51) or stream power (r2 = 0·59). Published in 1999 by John Wiley & Sons, Ltd.     

Magnitude–intensity and intensity–attenuation relationships for atlas region and Algerian earthquakes

This paper presents the results of an investigation of the magnitude–intensity and intensity–attenuation relationships for earthquakes in the Atlas block and Algeria using macroseismic data. This work is based on a selected sample of isoseismal maps from 32 events which were recently revised. Surface-wave magnitudes, M_s, are recalculated using the Prague formula and range from 4·2 to 7·45. Because the Atlas mountains block is in a collision zone, earthquakes occur in general within a layer 15 km deep. Expressions of general form for the magnitude–intensity and intensity–attenuation correlations are adopted and are, respectively, and where R² = d² + h², d the source distance in km, h the focal depth in km, M_s the revised surface-wave magnitude, M_sc the predicted surface-wave magnitude, I_i the intensity at isoseismal i, I the predicted intensity, σ the standard deviation and P is zero for 50-percentile values and one for 84-percentile, and the coefficients A's and B's are determined by regression analysis. The results of this study show that the intensity–attenuation models are adequate to predict quite well the die-out of intensity with distance in the Atlas zone and coastal Algeria; it is also found that magnitude can be predicted accurately by calibrating isoseismal radii against revised instrumental surface-wave magnitude. Such magnitude–intensity relationships may be used to evaluate the magnitude of historical earthquakes in the region under survey, with no instrumental data, for which isoseismal radii and intensities are available.     

Testing alternative erosivity indices to predict event soil loss from bare plots in Southern Italy

Methods for predicting unit plot soil loss for the ‘Sparacia’ Sicilian (Southern Italy) site were developed using 316 simultaneous measurements of runoff and soil loss from individual bare plots varying in length from 11 to 44 m. The event unit plot soil loss was directly proportional to an erosivity index equal to (Q_REI₃₀)^1·47, being Q_REI₃₀ the runoff ratio (Q_R) times the single storm erosion index (EI₃₀). The developed relationship represents a modified version of the USLE‐M, and therefore it was named USLE‐MM. By the USLE‐MM, a constant erodibility coefficient was deduced for plots of different lengths, suggesting that in this case the calculated erodibility factor is representative of an intrinsic soil property. Testing the USLE‐M and USLE‐MM schemes for other soils and developing simple procedures for estimating the plot runoff ratio has practical importance to develop a simple method to predict soil loss from bare plots at the erosive event temporal scale. Copyright © 2009 John Wiley & Sons, Ltd.     

Groundwater quality in Mesogea basin in eastern Attica (Greece)

Studies on the hydrogeological conditions of the Mesogea basin in east Attica reveal that the aquifers developed on the post‐alpine formations at the inner part of the coastal brackish zone exhibit positive hydraulic head. These Neogene and Quaternary deposits present high salt concentrations. Selected points were sampled (total 85: 51 wells and 34 boreholes) in order to obtain hydrogeological and hydrochemical data for a better understanding of the structure, operation and dynamics of the aquifer of the area. Statistical methods, R‐mode factor analysis and scatter‐plot diagrams were used for the hydrochemical analysis and presentation of the data. The groundwater resources are relatively weak and there is significant quality degradation due to the geological structure of the greater area, as well as the bad management of the aquifer and anthropogenic activities. Groundwater is characterized by high salt concentrations. Electrical conductivity values range between 260 and 6970 µS cm^?1. High salt concentrations at the coastal aquifers are due to sea intrusion, whereas they are attributed to the dissolution of minerals of the geological environment in the inland area. The groundwaters of the study area can be classified into five water types: Ca–HCO₃, Mg–HCO₃, Na–HCO₃, Na–Cl and Mg–Cl. They are saturated in dolomite and calcite, whereas they are unsaturated in anhydrite. High ion concentrations, e.g. ] (0‐221 mg l^?1), ] (0·01‐1·88 mg l^?1), ] (0·01‐6·75 mg l^?1), as well as high heavy metals concentrations are attributed to anthropogenic impacts. Copyright © 2006 John Wiley & Sons, Ltd.     

The spatial variation of weathering and soil depth across a Triassic sandstone outcrop

In this study, we used an archive of borehole logs from the British Geological Survey to collect information on the spatial structure of weathering that extends from the surface to competent bedrock across the Triassic Sherwood Sandstone Group outcrop (750 km²), in the East Midlands, UK. The borehole logs were used to estimate the thickness of the soil (n = 280) and soil and saprolite (S&S) to competent rock (n = 500). The weathering profile of the sandstone consisted of soil (median thickness ～ 1·5 m) overlying a transition zone of compacted and weakly cemented weathered sandstone saprolite over bedrock. Topographic analysis using a NEXTMAP 5 m × 5 m digital elevation model (DEM) revealed no significant relationships between slope properties (relief, flow length, flow accumulation or slope angle) and soil or S&S thickness. A weak, but statistically significant correlation was found between the thickness of the soil and S&S (r_s = 0·25, p < 0·001, n = 192). The variation in soil thickness may be related to changes in current and historic and land‐use, variation in sandstone properties and the influence of glacial/peri‐glacial processes. The thickness of the saprolite was more variable towards the southern part of the study area, where it increased to a maximum 40 m. We hypothesize and provide evidence that the greater weathering thickness is related to the occurrence of increased faulting in this part of the study region, allowing increased access to meteoric waters. A possible source of increased water supply is meltwater from Quaternary ice sheets; the overburden of ice may have increased sub‐glacial pore water pressure, with the fractures and faults acting as a drainage system for the removal of dissolved weathering products. British Geological Survey © NERC 2010