Hydrological processes of a closed catchment‐lake system in a semi‐arid Mediterranean environment

Data collected in 4 years of field observations were used in conjunction with continuous simulation models to study, at the small‐basin scale, the water balance of a closed catchment‐lake system in a semi‐arid Mediterranean environment. The open water evaporation was computed with the Penman equation, using the data set collected in the middle of the lake. The surface runoff was partly measured at the main tributary and partly simulated using a distributed, catchment, hydrological model, calibrated with the observed discharge. The simplified structure of the developed modelling mainly concerns soil moisture dynamics and bedrock hydraulics, whereas the flow components are physically based. The calibration produced high efficiency coefficients and showed that surface runoff is greatly affected by soil water percolation into fractured bedrock. The bedrock reduces the storm‐flow peaks and the interflow and has important multi‐year effects on the annual runoff coefficients. The net subsurface outflow from the lake was calculated as the residual of the lake water balance. It was almost constant in the dry seasons and increased in the wet seasons, because of the moistening of the unsaturated soil. During the years of observation, rainfall 30% higher than average caused abundant runoff and a continuous rise in the lake water levels. The analysis allows to predict that, in years with lower than the average rainfall, runoff will be drastically reduced and will not be able to compensate for negative balance between precipitation and lake evaporation. Such highly unsteady situations, with great fluctuations in lake levels, are typical of closed catchment‐lake systems in the semi‐arid Mediterranean environment. Copyright © 2012 John Wiley & Sons, Ltd.     

A hillslope hydrology approach for catchment‐scale slope stability analysis

Regional analysis of slope stability is often constrained by availability of data. Model requirements for input data cannot be met at the desired spatial resolution because data are either site‐speci?c or non‐existent. Faced with these dif?culties it has often been the practice to assume that certain parameters are uniform throughout the area of interest. An alternative approach proposed here allows a more detailed discrimination of slope stability conditions. Based on the principles of hillslope hydrology, hydrologic information can be generated at suf?cient resolution to allow higher resolution slope stability analysis. Measurements from an instrumented network in a small area have been used to establish index‐based models for topographic and climate‐related controls of piezometric response. The ability to relate groundwater levels to rainfall and topographic parameters provides a means of up‐scaling to larger catchments and ultimately the opportunity to generate a catchment‐wide prediction of the distribution, magnitude and frequency of rainstorm‐generated groundwater levels. The example provided in this study uses the topography index of TOPMODEL in GIS to predict the spatial patterns of groundwater elevation for seasonal soil moisture conditions and given rainfall inputs. This allows modelling of catchment‐wide response of soil water to rainstorms with different return periods (representing different magnitudes), and is an essential prerequisite for a probabilistic regional slope stability analysis. Copyright © 2004 John Wiley & Sons, Ltd.     

Scale relationships in hillslope runoff and erosion

Eight runoff plots, located within a small catchment within the Walnut Gulch Experimental Watershed, southern Arizona, were constructed to test the argument that sediment yield (kg m^?2) decreases as plot length increases. The plots ranged in length from 2 m to 27·78 m. Runoff and sediment loss from these plots were obtained for ten natural storm events. The pattern of sediment yield from these plots conforms to the case in which sediment yield first increases as plot length increases, but then subsequently decreases. Data from the present experiment indicate that maximum sediment yield would occur from a plot 7 m long. Analysis of both runoff and sediment yield from the plots indicates that the relationship of sediment yield to plot length derives both from the limited travel distance of individual entrained particles and from a decline in runoff coefficient as plot length increases. Particle‐size analysis of eroded sediment confirms the role of travel distance in controlling sediment yield. Whether in response to the finite travel distance of entrained particles or the relationship of runoff coefficient to plot length, the experiment clearly demonstrates that the erosion rates for hillslopes and catchments cannot be simply extrapolated from plot measurements, and that alternative methods for estimating large‐area erosion rates are required. Copyright © 2006 John Wiley & Sons, Ltd.     

Modelling catchment‐scale erosion patterns in the East African Highlands

Prompt location of areas exposed to high erosion is of the utmost importance for soil and water conservation planning. Erosion models can be useful tools to locate sources of sediment and areas of deposition within a catchment, but the reliability of model predictions of spatial patterns of erosion at catchment scale has seldom been validated against observations. This study aimed to evaluate the performance of a simple empirical model (Morgan, Morgan and Finney model, MMF) in predicting spatial patterns of erosion at two small catchments in the East African Highlands: Kwalei (Tanzania) and Gikuuri (Kenya). Erosion maps predicted by the MMF model were compared with erosion maps obtained by direct survey. In Kwalei, erosion features were especially frequent in fields of annual crops. In Gikuuri, slope was the critical erosion factor, with estimated erosion rates >10 kg m^?2 a^?1 on slopes >18 per cent. Predicted erosion rates were mainly transport‐limited and ranged from <0·01 to 13·50 kg m^?2 a^?1 in Kwalei and 9·29 kg m^?2 a^?1 in Gikuuri. The performance of the MMF model in predicting the spatial patterns of erosion was acceptable in Kwalei, but poor in Gikuuri. However, by excluding the elements at the valley bottoms in Gikuuri Catchment, the performance of the model improved dramatically. The spatial pattern of erosion predicted by the MMF model was driven by the accumulation of surface runoff, which did not consider the possibility of re‐infiltration along the slope. As a result, the MMF erosion patterns predicted by the model increased invariably from the ridges to the valley bottoms, hampering the model suitability for locating areas subjected to high and very high erosion. It is concluded that the model predictions could be substantially improved by introducing a more realistic hydrological component for the prediction of surface runoff along the hill‐slope. Copyright © 2005 John Wiley & Sons, Ltd.     

A simplified model for estimating field‐scale surface runoff hydrographs

The problem of obtaining field‐scale surface response to rainfall events is complicated by the spatial variability of infiltration characteristics of the soil and rainfall. In this paper, we develop and test a simplified model for generating surface runoff over fields with spatial variation in both rainfall rate and saturated hydraulic conductivities. The model is able to represent the effects of local variation in infiltration, as well as the run‐on effect that controls infiltration of excess water from saturated upstream areas. The effective rainfall excess is routed to the slope outlet using a simplified solution of the kinematic wave approximation. Model results are compared to averaged hydrographs from numerically‐intensive Monte–Carlo simulations for observed and design rainfall events and soil patterns that are typical of Central Italy. The simplified model is found to yield satisfactory results at a relatively small computational expense. A proposal to include a simple channel routing scheme is also presented as a prelude to extend this conceptualization to watershed scales. Copyright © 2006 John Wiley & Sons, Ltd.     

Sensitivity of runoff and soil erosion to climate change in two Mediterranean watersheds. Part I: model parameterization and evaluation

Climate change is expected to effect storm runoff and erosion processes in Mediterranean watersheds at multiple spatial scales. Models are typically applied to estimate these impacts; however, the scarcity of spatially distributed data for parameterization, calibration and validation often prevents application of these models, particularly for larger catchments. This report, the first part of a two‐part article, presents an application and evaluation of the MEFIDIS model for two Mediterranean meso‐scale watersheds (115 and 290 km²) in a data‐scarce environment. A multi‐scale assessment method was used that combines quantitative validation and qualitative evaluation, consisting of three steps: (1) calibration at the small (field) scale using results from rainfall simulation experiments; (2) calibration and validation for catchment‐scale results while changing catchment‐scale parameters only (channel roughness and a parameter controlling the distribution of saturated areas); and (3) qualitative evaluation of within‐watershed erosion processes using empirical estimates of sediment delivery ratio and gully location. The results indicate that calibrating MEFIDIS at the field scale can provide reasonable results for catchment runoff and sediment export and for within‐watershed erosion processes. Copyright © 2009 John Wiley & Sons, Ltd.     

Capabilities and limitations of detailed hillslope hydrological modelling

Hillslope hydrological modelling is considered to be of great importance for the understanding and quantification of hydrological processes in hilly or mountainous landscapes. In recent years a few comprehensive hydrological models have been developed at the hillslope scale which have resulted in an advanced representation of hillslope hydrological processes (including their interactions), and in some operational applications, such as in runoff and erosion studies at the field scale or lateral flow simulation in environmental and geotechnical engineering. An overview of the objectives of hillslope hydrological modelling is given, followed by a brief introduction of an exemplary comprehensive hillslope model, which stimulates a series of hydrological processes such as interception, evapotranspiration, infiltration into the soil matrix and into macropores, lateral and vertical subsurface soil water flow both in the matrix and preferential flow paths, surface runoff and channel discharge. Several examples of this model are presented and discussed in order to determine the model's capabilities and limitations. Finally, conclusions about the limitations of detailed hillslope modelling are drawn and an outlook on the future prospects of hydrological models on the hillslope scale is given.The model presented performed reasonable calculations of Hortonian surface runoff and subsequent erosion processes, given detailed information of initial soil water content and soil hydraulic conditions. The vertical and lateral soil moisture dynamics were also represented quite well. However, the given examples of model applications show that quite detailed climatic and soil data are required to obtain satisfactory results. The limitations of detailed hillslope hydrological modelling arise from different points: difficulties in the representations of certain processes (e.g. surface crusting, unsaturated–saturated soil moisture flow, macropore flow), problems of small‐scale variability, a general scarcity of detailed soil data, incomplete process parametrization and problems with the interdependent linkage of several hillslopes and channel–hillslope interactions. Copyright © 1999 John Wiley & Sons, Ltd.     

Quantifying rainfall controls on catchment‐scale landslide erosion in Taiwan

Landslide erosion is a dominant hillslope process and the main source of stream sediment in tropical, tectonically active mountain belts. In this study, we quantified landslide erosion triggered by 24 rainfall events from 2001 to 2009 in three mountainous watersheds in Taiwan and investigated relationships between landslide erosion and rainfall variables. The results show positive power‐law relations between landslide erosion and rainfall intensity and cumulative rainfall, with scaling exponents ranging from 2·94 to 5·03. Additionally, landslide erosion caused by Typhoon Morakot is of comparable magnitude to landslide erosion caused by the Chi‐Chi Earthquake (M_W = 7·6) or 22–24 years of basin‐averaged erosion. Comparison of the three watersheds indicates that deeper landslides that mobilize soil and bedrock are triggered by long‐duration rainfall, whereas shallow landslides are triggered by short‐duration rainfall. These results suggest that rainfall intensity and watershed characteristics are important controls on rainfall‐triggered landslide erosion and that severe typhoons, like high‐magnitude earthquakes, can generate high rates of landslide erosion in Taiwan. Copyright © 2012 John Wiley & Sons, Ltd.     

Exploring effects of rainfall intensity and duration on soil erosion at the catchment scale using openLISEM: Prado catchment,SE Spain

In semi‐arid areas, high‐intensity rainfall events are often held responsible for the main part of soil erosion. Long‐term landscape evolution models usually use average annual rainfall as input, making the evaluation of single events impossible. Event‐based soil erosion models are better suited for this purpose but cannot be used to simulate longer timescales and are usually applied to plots or small catchments. In this study, the openLISEM event‐based erosion model was applied to the medium‐sized (～50 km²) Prado catchment in SE Spain. Our aim was to (i) test the model's performance for medium‐sized catchments, (ii) test the ability to simulate four selected typical Mediterranean rainfall events of different magnitude and (iii) explore the relative contribution of these different storms to soil erosion using scenarios of future climate variability. Results show that because of large differences in the hydrologic response between storms of different magnitudes, each event needed to be calibrated separately. The relation between rainfall event characteristics and the calibration factors might help in determining optimal calibration values if event characteristics are known. Calibration of the model features some drawbacks for large catchments due to spatial variability in K_sat values. Scenario calculations show that although ～50% of soil erosion occurs as a result of high frequency, low‐intensity rainfall events, large‐magnitude, low‐frequency events potentially contribute significantly to total soil erosion. The results illustrate the need to incorporate temporal variability in rainfall magnitude–frequency distributions in landscape evolution models. Copyright © 2011 John Wiley & Sons, Ltd.     

A catchment scale evaluation of the SIBERIA and CAESAR landscape evolution models

Landscape evolution models provide a way to determine erosion rates and landscape stability over times scales from tens to thousands of years. The SIBERIA and CAESAR landscape evolution models both have the capability to simulate catchment–wide erosion and deposition over these time scales. They are both cellular, operate over a digital elevation model of the landscape, and represent fluvial and slope processes. However, they were initially developed to solve research questions at different time and space scales and subsequently the perspective, detail and process representation vary considerably between the models. Notably, CAESAR simulates individual events with a greater emphasis on fluvial processes whereas SIBERIA averages erosion rates across annual time scales. This paper describes how both models are applied to Tin Camp Creek, Northern Territory, Australia, where soil erosion rates have been closely monitored over the last 10 years. Results simulating 10 000 years of erosion are similar, yet also pick up subtle differences that indicate the relative strengths and weaknesses of the two models. The results from both the SIBERIA and CAESAR models compare well with independent field data determined for the site over different time scales. Representative hillslope cross‐sections are very similar between the models. Geomorphologically there was little difference between the modelled catchments after 1000 years but significant differences were revealed at longer simulation times. Importantly, both models show that they are sensitive to input parameters and that hydrology and erosion parameter derivation has long‐term implications for sediment transport prediction. Therefore selection of input parameters is critical. This study also provides a good example of how different models may be better suited to different applications or research questions. Copyright © 2010 John Wiley & Sons, Ltd and Commonwealth of Australia     

Spatial and temporal variations in soil erosion and deposition due to land‐levelling in a semi‐arid area of Basilicata (Southern Italy)

This paper examines the changes from 1955 to 2002 in soil erosion and deposition due to changes in land‐use patterns in the semi‐arid territory of Craco, which is characterized by landsliding and badland erosion. The area underwent continuous degradation during the last century due not only to its lithological vulnerability but also to the anthropic pressure favoured by the introduction of Common Agricultural Policy (CAP) measures, which has led to the reclamation of scrub lands and badlands for durum wheat cultivation. Our analysis integrates the Unit Stream Power Erosion Deposition (USPED) model with a geographic information system (GIS) to quantify erosion risk and predict deposition patterns. Soil data, land use inventory, digital elevation data and climatic atlases were used as resource data sets to generate USPED factor values. The obtained results correlate well with field measured erosion data by other researchers. In the investigated 47 years, stable areas decreased by about 280 ha (3·8% of the total surface area), largely attributable to the increase of the low and moderate erosion intensity without significant change in sedimentation. Results from this study have implications related to understanding the geomorphic response of sites that were abandoned following remodelling due to the application of the F measure of Regulation CEE 2078/92. The average annual erosion rates estimated for abandoned and remodelled sites are respectively 15·99 and 10·64 t ha^?1, meaning that the total amount of erosion in 20 years could be estimated at around 100 t ha^?1. Copyright © 2007 John Wiley & Sons, Ltd.     

Application of a data‐based mechanistic modelling (DBM) approach for predicting runoff generation in semi‐arid regions

This paper addresses the application of a data‐based mechanistic (DBM) modelling approach using transfer function models (TFMs) with non‐linear rainfall filtering to predict runoff generation from a semi‐arid catchment (795 km²) in Tanzania. With DBM modelling, time series of rainfall and streamflow were allowed to suggest an appropriate model structure compatible with the data available. The model structures were evaluated by looking at how well the model fitted the data, and how well the parameters of the model were estimated. The results indicated that a parallel model structure is appropriate with a proportion of the runoff being routed through a fast flow pathway and the remainder through a slow flow pathway. Finally, the study employed a Generalized Likelihood Uncertainty Estimation (GLUE) methodology to evaluate the parameter sensitivity and predictive uncertainty based on the feasible parameter ranges chosen from the initial analysis of recession curves and calibration of the TFM. Results showed that parameters that control the slow flow pathway are relatively more sensitive than those that control the fast flow pathway of the hydrograph. Within the GLUE framework, it was found that multiple acceptable parameter sets give a range of predictions. This was found to be an advantage, since it allows the possibility of assessing the uncertainty in predictions as conditioned on the calibration data and then using that uncertainty as part of the decision‐making process arising from any rainfall‐runoff modelling project. Copyright © 2001 John Wiley & Sons, Ltd.     

Precipitation‐runoff processes in Shimen hillslope micro‐catchment of Taihang Mountain,north China

Precipitation runoff is a critical hillslope hydrological process for downslope streamflow and piedmont/floodplain recharge. Shimen hillslope micro‐catchment is strategically located in the central foothill region of Taihang Mountains, where runoff is crucial for water availability in the piedmont corridors and floodplains of north China. This study analyzes precipitation‐runoff processes in the Shimen hillslope micro‐catchment for 2006–2008 using locally designed runoff collection systems. The study shows that slope length is a critical factor, next only to precipitation, in terms of runoff yield. Regression analysis also shows that runoff is related positively to precipitation, and negatively to slope length. Soil mantle in the study area is generally thin and is therefore not as critical a runoff factor as slope length. The study shows a significant difference between overland and subsurface runoff. However, that between the 0–10 and 10–20 cm subsurfaces is insignificant. Runoff hardly occurs under light rains (<10 mm), but is clearly noticeable under moderate‐to‐rainstorm events. In the hillslope catchment, vertical infiltration (accounting for 42–84% of the precipitation) dominates runoff processes in subsurface soils and weathered granite gneiss bedrock. A weak lateral flow (at even the soil/bedrock interface) and the generally small runoff suggest strong infiltration loss via deep percolation. This is critical for groundwater recharge in the downslope piedmont corridors and floodplains. This may enhance water availability, ease water shortage, avert further environmental degradation, and reduce the risk of drought/flood in the event of extreme weather conditions in the catchment and the wider north China Plain. Copyright © 2011 John Wiley & Sons, Ltd.     

Improving the predictions of a MIKE SHE catchment‐scale application by using a multi‐criteria approach

This paper suggests a multi‐criteria protocol for appropriately evaluating the predictions of hydrologic models during calibration and evaluation stages. The protocol includes different statistical, analytical and visual criteria such as analysis of peak and low flows, cumulative volumes, extreme value statistics, performance statistics, etc. Furthermore, the protocol assesses the physical consistency of model predictions by filtering the total observed hydrograph into different flow‐components (baseflow, interflow and overland flow) and using these filtered data in the calibration and evaluation processes. Based on the distributed modelling of a medium size catchment, it is shown that application of the suggested protocol, and in particular the use of the filtered flow‐components in model calibration, enhances the physical consistency of model predictions, adding considerable value to the calibration process. Copyright © 2007 John Wiley & Sons, Ltd.     

Implications of model uncertainty for the mapping of hillslope‐scale soil erosion predictions

An integrated modelling approach (MIRSED) which utilizes the process‐based soil erosion model WEPP (Water Erosion Prediction Project) is presented for the assessment of hillslope‐scale soil erosion at five sites throughout England and Wales. The methodology draws upon previous uncertainty analysis of the WEPP hillslope soil erosion model by the authors to qualify model results within an uncertainty framework. A method for incorporating model uncertainty from a range of sources is discussed as a first step towards using and learning from results produced through the GLUE (Generalized Likelihood Uncertainty Estimation) technique. Results are presented and compared to available observed data, which illustrate that levels of uncertainty are significant and must be taken into account if a meaningful understanding of output from models such as WEPP is to be achieved. Furthermore, the collection of quality, observed data is underlined for two reasons: as an essential tool in the development of soil erosion modelling and also to allow further constraint of model uncertainty. Copyright © 2001 John Wiley & Sons, Ltd.     

Estimating runoff from a glacierized catchment using natural tracers in the semi‐arid Andes cordillera

This paper presents a methodology for hydrograph separation in mountain watersheds, which aims at identifying flow sources among ungauged headwater sub‐catchments through a combination of observed streamflow and data on natural tracers including isotope and dissolved solids. Daily summer and bi‐daily spring season water samples obtained at the outlet of the Juncal River Basin in the Andes of Central Chile were analysed for all major ions as well as stable water isotopes, δ¹⁸O and δD. Additionally, various samples from rain, snow, surface streams and exfiltrating subsurface water (springs) were sampled throughout the catchment. A principal component analysis was performed in order to address cross‐correlation in the tracer dataset, reduce the dimensionality of the problem and uncover patterns of variability. Potential sources were identified in a two‐component U‐space that explains 94% of the observed tracer variability at the catchment outlet. Hydrograph separation was performed through an Informative‐Bayesian model. Our results indicate that the Juncal Norte Glacier headwater sub‐catchment contributed at least 50% of summer flows at the Juncal River Basin outlet during the 2011–2012 water year (a hydrologically dry period in the Region), even though it accounts for only 27% of the basin area. Our study confirms the value of combining solute and isotope information for estimating source contributions in complex hydrologic systems, and provides insights regarding experimental design in high‐elevation semi‐arid catchments. The findings of this study can be useful for evaluating modelling studies of the hydrological consequences of the rapid decrease in glacier cover observed in this region, by providing insights into the origin of river water in basins with little hydrometeorological information. Copyright © 2016 John Wiley & Sons, Ltd.     

Comparison of two stream gauging systems for measuring runoff and sediment yield for a semi‐arid watershed

Our ability to understand erosion processes in semi‐arid ecosystems depends on establishing relationships between rainfall and runoff. This requires collection of extensive and accurate hydrologic and sediment data sets. A supercritical flume with a total load traversing slot sediment sampler used on several sites at the Walnut Gulch Experimental Watershed (WGEW) near Tombstone, AZ has proven to be a reliable way to measure flow and sediment discharge from small watersheds. However, it requires installation of a costly structure that is only suitable for relatively small flows. A more commonly used method based on ease of installation and expense is the pump sampler. One example of this is a set of instrumentation developed by the Australian Commonwealth Scientific and Industrial Research Organization (CSIRO), in which the pump sediment sampler is part of an in‐channel, fully automated system for measuring water velocity, depth, turbidity and collecting runoff samples. A 3.7 ha arid watershed at WGEW was instrumented with both systems and hydrologic and sediment data were collected and compared during a 2 year period. Total sediment yield for the entire period measured by the CSIRO pump sampler (11.6 t ha^‐1) was similar to that by traversing slot sampler (11.5 t ha^‐1). The pump sampler accurately estimated the amount of fine (< 0.5 mm) sediment fractions exported, but consistently underestimated the coarse (>0.5 mm) sediment fractions. Median sediment diameter of samples collected by traversing slot and pump sampler were 0.32 and 0.22 mm, respectively. This study outlines the benefits and limitations of the pump sampler based system for monitoring sediment concentration and yield in high‐energy headwater catchments, and makes recommendations for improvement of its performance. Copyright © 2012 John Wiley & Sons, Ltd.     

Effects of hillslope topography on runoff response in a small catchment in the Fudoji Experimental Watershed,central Japan

We investigated the role of different hillslope units with different topographic characteristics on runoff generation processes based on field observations at two types of hillslopes (0·1 ha): a valley‐head (a convergent hillslope) and a side slope (a planar hillslope), as well as at three small catchments having two types of slopes with different drainage areas ranging from 1·9 to 49·7 ha in the Tanakami Mountains, central Japan. We found that the contribution of the hillslope unit type to small catchment runoff varied with the magnitude of rainfall. When the total amount of rainfall for a single storm event was < 35 mm, runoff in the small catchment was predominantly generated from the side slope. As the amount of rainfall increased (>35 mm), the valley‐head also began to contribute to the catchment runoff, adding to runoff from the side slope. Although the direct runoff from the valley‐head was greater than that from the side slope, the contribution from the side slope was quantitatively greater than that from the valley‐head due to the proportionally larger area occupied by the side slope in the small catchment. The storm runoff responses of the small catchments reflected the change in the runoff components of each hillslope unit as the amount of rainfall increased and rainfall patterns changed. However, similar runoff responses were found for the small catchments with different areas. The similarity of the runoff responses is attributable to overlay effects of different hillslope units and the similar composition ratios of the valley‐head and side slope in the catchments. This study suggests that the relative roles of the valley‐head and side slope are important in runoff generation and solute transport as the catchment size increases from a hillslope/headwater to a small catchment. Copyright © 2011 John Wiley & Sons, Ltd.     

Determining soil erosion rates on semi‐arid watersheds using radioisotope‐derived sedimentation chronology

This study investigates erosion dynamics of the past 90 years in three small semi‐arid watersheds with histories of grazing and vegetation change. Activity of ¹³⁷Cs and excess ²¹⁰Pb from 18 cores collected from sedimentation ponds were measured using a gamma spectrometer. The sediment was dated using a constant rate of supply (CRS) model. This study represents the first time that reservoir sediment accumulation rates determined from fallout isotopes have been verified by direct volumetric measurements of aggradation based on topographic surveys. Measured sedimentation in the ponds ranged between 1.9 and 2.3 cm y^?1, representing average sediment delivery rates from the watersheds of between 0.6 and 2.0 t ha^?1 y^?1. These sediment delivery rates were in agreement with those established by other methods for similar catchments in the region. Past variations in sedimentation rates were identified and correlated with recorded history of anthropogenic disturbance. ¹³⁷Cs and ²¹⁰Pb methods are suitable for use in arid environments and can complement each other to increase reliability of erosion rate estimates. The abundance of stock ponds in southwestern USA presents an opportunity to quantify historic erosion and sediment transfer dynamics in areas that have not been well studied or instrumented. Published 2016. This article is a U.S. Government work and is in the public domain in the USA