A new least-squares minimization approach to depth and shape determination from magnetic data

We have developed a least‐squares minimization approach to depth determination using numerical second horizontal derivative anomalies obtained from magnetic data with filters of successive window lengths (graticule spacings). The problem of depth determination from second‐derivative magnetic anomalies has been transformed into finding a solution to a non‐linear equation of the form, f(z) = 0. Formulae have been derived for a sphere, a horizontal cylinder, a dike and a geological contact. Procedures are also formulated to estimate the magnetic angle and the amplitude coefficient. We have also developed a simple method to define simultaneously the shape (shape factor) and the depth of a buried structure from magnetic data. The method is based on computing the variance of depths determined from all second‐derivative anomaly profiles using the above method. The variance is considered a criterion for determining the correct shape and depth of the buried structure. When the correct shape factor is used, the variance of depths is less than the variances computed using incorrect shape factors. The method is applied to synthetic data with and without random errors, complicated regionals, and interference from neighbouring magnetic rocks. Finally, the method is tested on a field example from India. In all the cases examined, the depth and the shape parameters are found to be in good agreement with the actual parameters.     

High‐accuracy practical spline‐based 3D and 2D integral transformations in potential‐field geophysics

Potential, potential field and potential‐field gradient data are supplemental to each other for resolving sources of interest in both exploration and solid Earth studies. We propose flexible high‐accuracy practical techniques to perform 3D and 2D integral transformations from potential field components to potential and from potential‐field gradient components to potential field components in the space domain using cubic B‐splines. The spline techniques are applicable to either uniform or non‐uniform rectangular grids for the 3D case, and applicable to either regular or irregular grids for the 2D case. The spline‐based indefinite integrations can be computed at any point in the computational domain. In our synthetic 3D gravity and magnetic transformation examples, we show that the spline techniques are substantially more accurate than the Fourier transform techniques, and demonstrate that harmonicity is confirmed substantially better for the spline method than the Fourier transform method and that spline‐based integration and differentiation are invertible. The cost of the increase in accuracy is an increase in computing time. Our real data examples of 3D transformations show that the spline‐based results agree substantially better or better with the observed data than do the Fourier‐based results. The spline techniques would therefore be very useful for data quality control through comparisons of the computed and observed components. If certain desired components of the potential field or gradient data are not measured, they can be obtained using the spline‐based transformations as alternatives to the Fourier transform techniques.     

Combining 3D Hydraulic Tomography with Tracer Tests for Improved Transport Characterization

Hydraulic tomography (HT) is a method for resolving the spatial distribution of hydraulic parameters to some extent, but many details important for solute transport usually remain unresolved. We present a methodology to improve solute transport predictions by combining data from HT with the breakthrough curve (BTC) of a single forced‐gradient tracer test. We estimated the three dimensional (3D) hydraulic‐conductivity field in an alluvial aquifer by inverting tomographic pumping tests performed at the Hydrogeological Research Site Lauswiesen close to Tübingen, Germany, using a regularized pilot‐point method. We compared the estimated parameter field to available profiles of hydraulic‐conductivity variations from direct‐push injection logging (DPIL), and validated the hydraulic‐conductivity field with hydraulic‐head measurements of tests not used in the inversion. After validation, spatially uniform parameters for dual‐domain transport were estimated by fitting tracer data collected during a forced‐gradient tracer test. The dual‐domain assumption was used to parameterize effects of the unresolved heterogeneity of the aquifer and deemed necessary to fit the shape of the BTC using reasonable parameter values. The estimated hydraulic‐conductivity field and transport parameters were subsequently used to successfully predict a second independent tracer test. Our work provides an efficient and practical approach to predict solute transport in heterogeneous aquifers without performing elaborate field tracer tests with a tomographic layout.     

2D and 3D potential-field upward continuation using splines

The dominant upward‐continuation technique used in the potential‐field geophysics industry is the fast Fourier transform (FFT) technique. However, the spline‐based upward‐continuation technique presented in this paper has some advantages over the FFT technique. The spline technique can be used to carry out level‐to‐uneven surface 2D and 3D potential‐field upward continuation. An example of level‐to‐uneven surface upward continuation of 3D magnetic data using the spline technique is shown, and it is evident that the continued anomalies are very close to the theoretical values. The spacing can be irregular. Synthetic examples using the spline technique to continue noise‐contaminated gravity and magnetic data upward to an altitude of 15 km on irregular grids are shown. Gaussian noise with a zero mean and a standard deviation of 1% does not cause much error and can readily be tolerated. Through comparison with the FFT technique, it is found that for low‐altitude gravity and magnetic upward continuation, both the FFT technique and the spline technique are suitable; for high‐altitude upward continuation, the FFT technique is inaccurate, whereas the spline technique works very well. Also, upward continuation by the spline technique has a smaller edge effect than upward continuation by the FFT technique. The spline‐based upward continuation technique works fairly well even when the periphery of a grid is not quiet: it is rather robust in general. A real example shows that the spline technique can be employed to perform upward continuation of total‐field magnetic data and to de‐emphasize near‐surface noise.     

An optimal Savitzky–Golay derivative filter with geophysical applications: an example of self-potential data

We propose a strategy in designing an optimal set of filter parameters, such as the order of interpolating polynomial and the filter length for a Savitzky–Golay derivative filter. The proposed strategy is based on the ‘principle of parsimony’ while satisfying the optimality conditions. The optimality conditions are based on the Durbin–Watson lag-1 test statistic and the Derringer–Suich desirability function. While the former checks for an appropriate data fitting, the latter, on the other hand, ensures minimal shape distortion of the reconstructed response. The proposed strategy of designing filter parameters is developed and validated through numerical experiments using Gaussian pulse as a test function which is contaminated with additive white Gaussian noise. In the numerical tests, the polynomial orders used were 3, 5 and 7, but the filter length for each polynomial was varying until the optimality conditions were satisfied. The Savitzky–Golay derivative filtering is used in obtaining the robust reconstruction of noisy geophysical anomaly and the robust estimation of its first- and second-order derivatives. We validated the proposed technique on the published self-potential anomaly data using a data-based interpretation technique where the reconstructed anomaly and its first- and second-order derivatives were used in estimating model parameters. The data-based interpretation using the proposed technique of Savitzky–Golay derivative filtering provides a close agreement with the published results.     

基于深度卷积神经网络的地震数据局部信噪比估计

估计地震数据的信噪比对于地震数据的处理和解释具有重要作用.以往估计地震数据信噪比的方法都需要分离数据中的有效信号和噪声,然后再估计相应的信噪比.这些估计方法的精度严重依赖信号估计方法或噪声压制方法的有效性,往往存在偏差.本文提出一种估计地震数据局部信噪比的深度卷积神经网络模型,通过迭代训练优化参数,构建从含噪地震数据到其信噪比的特征映射.然后使用该神经网络完成信噪比的推理预测,不需要分离地震数据中的有效信号和噪声.模拟数据和实际资料的处理结果都表明,本文的方法可以准确而高效地估计局部地震数据的信噪比,为地震数据质量的定量评价提供依据.     

Surface nuclear magnetic resonance signals recovery by integration of a non‐linear decomposition method with statistical analysis

Presence of noise in the acquisition of surface nuclear magnetic resonance data is inevitable. There are various types of noise, including Gaussian noise, spiky events, and harmonic noise that affect the signal quality of surface nuclear magnetic resonance measurements. In this paper, we describe an application of a two‐step noise suppression approach based on a non‐linear adaptive decomposition technique called complete ensemble empirical mode decomposition in conjunction with a statistical optimization process for enhancing the signal‐to‐noise ratio of the surface nuclear magnetic resonance signal. The filtering procedure starts with applying the complete ensemble empirical mode decomposition method to decompose the noisy surface nuclear magnetic resonance signal into a finite number of intrinsic mode functions. Afterwards, a threshold region based on de‐trended fluctuation analysis is defined to identify the noisy intrinsic mode functions, and then the no‐noise intrinsic mode functions are used to recover the partially de‐noised signal. In the second stage, we applied a statistical method based on the variance criterion to the signal obtained from the initial phase to mitigate the remaining noise. To demonstrate the functionality of the proposed strategy, the method was evaluated on an added‐noise synthetic surface nuclear magnetic resonance signal and on field data. The results show that the proposed procedure allows us to improve the signal‐to‐noise ratio significantly and, consequently, extract the signal parameters (i.e., and V₀) from noisy surface nuclear magnetic resonance data efficiently.     

Implementation of online model updating in hybrid simulation

Hybrid simulation combines numerical and experimental methods for cost‐effective, large‐scale testing of structures under simulated earthquake loading. Structural system level response can be obtained by expressing the equation of motion for the combined experimental and numerical substructures, and solved using time‐stepping integration similar to pure numerical simulations. It is often assumed that a reliable model exists for the numerical substructures while the experimental substructures correspond to parts of the structure that are difficult to model. A wealth of data becomes available during the simulation from the measured experiment response that can be used to improve upon the numerical models, particularly if a component with similar structural configuration and material properties is being tested and subjected to a comparable load pattern. To take advantage of experimental measurements, a new hybrid test framework is proposed with an updating scheme to update the initial modeling parameters of the numerical model based on the instantaneously‐measured response of the experimental substructures as the test progresses. Numerical simulations are first conducted to evaluate key algorithms for the selection and calibration of modeling parameters that can be updated. The framework is then expanded to conduct actual hybrid simulations of a structural frame model including a physical substructure in the laboratory and a numerical substructure that is updated during the tests. The effectiveness of the proposed framework is demonstrated for a simple frame structure but is extendable to more complex structural behavior and models. Copyright © 2013 John Wiley & Sons, Ltd.     

Investigating the transport dynamics and the properties of bedload material with a hydro‐acoustic measuring system

This article deals with the following two questions. Are acoustic measurements in running waters appropriate for a highly resolved investigation of the bedload transport? Which characterizations of the bedload regarding mass and shape are possible via the acoustic signals? The signals were recorded by means of data recorders (Tascam Inc. DAP1 Portable Data Recorder) and hydrophones (International Transducer Corp. ITC‐4001 A). The ITC‐4001 is a shallow water omnidirectional transducer containing a flexural disc transducer utilizing Channelite‐5400 ceramics mounted in a rugged corrosion‐resistant housing. These hydrophones were screwed onto the bottom side of stainless steel plates, serving as a contact surface for the bedload in motion above them. After more than 100 series of tests in the laboratory, which indicated the basic relations between the dimension, shape and weight of the bedload and the resulting signal, field tests of the measuring system were conducted. By artificially produced flood waves in the small brooks Riverisbach, Olewiger Bach and by a winter flood wave in the River Moselle, it is possible to elaborate similar structures of the signal course of the bedload movement. The highest transport rates can be observed at the beginning of the increasing limbs and behind the peaks of the waves. At the beginning of the waves, the increasing transport power of the water and the loose material can be considered as the cause for this result. The high stream velocity behind the wave peaks explains the increase in the bedload transport so that material from the channel beds is unfastened and will be mobilized. The characterization of the bedload regarding the shape and mass is still limited regarding the field measurements and could be solved only for homogeneous grain sizes and single stones under laboratory conditions. Copyright © 2007 John Wiley & Sons, Ltd.     

基于隐马尔可夫模型平滑估计的随机噪声压制方法

以地震勘探记录去噪为目标,本文提出了一种隐马尔可夫模型平滑估计方法.它是在基本隐马尔可夫模型滤波基础之上,运用信号检测环节将带噪信号段和无信号段加以区分,构建带噪地震记录的状态转移模型,在贝叶斯框架下,利用平滑密度函数进行状态估计,从而达到压制噪声的目的.数值模拟表明,无论对信噪比还是均方误差,隐马尔可夫模型平滑估计处理后的重构信号优于常规的维纳滤波所恢复信号.我们可以期待这种方法会成为实际地震记录噪声压制的有效手段.     

Assessing burn severity and comparing soil water repellency,Hayman Fire,Colorado

An important element of evaluating a large wildfire is to assess its effects on the soil in order to predict the potential watershed response. After the 55 000 ha Hayman Fire on the Colorado Front Range, 24 soil and vegetation variables were measured to determine the key variables that could be used for a rapid field assessment of burn severity. The percentage of exposed mineral soil and litter cover proved to be the best predictors of burn severity in this environment. Two burn severity classifications, one from a statistical classification tree and the other a Burned Area Emergency Response (BAER) burn severity map, were compared with measured ‘ground truth’ burn severity at 183 plots and were 56% and 69% accurate, respectively. This study also compared water repellency measurements made with the water drop penetration time (WDPT) test and a mini‐disk infiltrometer (MDI) test. At the soil surface, the moderate and highly burned sites had the strongest water repellency, yet were not significantly different from each other. Areas burned at moderate severity had 1·5 times more plots that were strongly water repellent at the surface than the areas burned at high severity. However, the high severity plots most likely had a deeper water repellent layer that was not detected with our surface tests. The WDPT and MDI values had an overall correlation of r = ?0·64(p < 0·0001) and appeared to be compatible methods for assessing soil water repellency in the field. Both tests represent point measurements of a soil characteristic that has large spatial variability; hence, results from both tests reflect that variability, accounting for much of the remaining variance. The MDI is easier to use, takes about 1 min to assess a strongly water repellent soil and provides two indicators of water repellency: the time to start of infiltration and a relative infiltration rate. Copyright © 2005 John Wiley & Sons, Ltd.     

Range imaging: a new method for high‐resolution topographic measurements in small‐ and medium‐scale field sites

Topographic measurements are essential for the study of earth surface processes. Three‐dimensional data have been conventionally obtained through terrestrial laser scanning or photogrammetric methods. However, particularly in steep and rough terrain, high‐resolution field measurements remain challenging and often require new creative approaches. In this paper, range imaging is evaluated as an alternative method for obtaining surface data in such complex environments. Range imaging is an emerging time‐of‐flight technology, using phase shift measurements on a multi‐pixel sensor to generate a distance image of a surface. Its suitability for field measurements has yet not been tested. We found ambient light and surface reflectivity to be the main factors affecting error in distance measurements. Low‐reflectivity surfaces and strong illumination contrasts under direct exposure to sunlight lead to noisy distance measurements. However, regardless of lighting conditions, the accuracy of range imaging was markedly improved by averaging multiple images of the same scene. For medium ambient lighting (shade) and a light‐coloured surface the measurement uncertainty was approximately 9 mm. To further test the suitability of range imaging for field applications we measured a reach of a steep mountain stream with a horizontal resolution of approximately 1 cm (in the focal plane of the camera), allowing for the interpolation of a digital elevation model on a 2 cm grid. Comparison with an elevation model obtained from terrestrial laser scanning for the same site revealed that both models show similar degrees of topographic detail. Despite limitations in measurement range and accuracy, particularly at bright ambient lighting, range imaging offers three‐dimensional data in real time and video mode without the need of post‐processing. Therefore, range imaging is a useful complement or alternative to existing methods for high‐resolution measurements in small‐ to medium‐scale field sites. Copyright © 2012 John Wiley & Sons, Ltd.     

Reduction of the number of knots in the approximation by smoothed cubic splines

Summary A method of approximation, minimizing the number of knots and satisfying the smoothing and fitting properties of the cubic spline approximation, is suggested. Two modifications of the method are presented. The first approach is specified to provide an adequate simple approximation to the digitalized smooth curve (such as the velocity-depth function in seismology), while the second is more general and may be used to approximate geophysical measurements, which may include errors.     

Cross-borehole flowmeter tests for transient heads in heterogeneous aquifers

Cross-borehole flowmeter tests have been proposed as an efficient method to investigate preferential flowpaths in heterogeneous aquifers, which is a major task in the characterization of fractured aquifers. Cross-borehole flowmeter tests are based on the idea that changing the pumping conditions in a given aquifer will modify the hydraulic head distribution in large-scale flowpaths, producing measurable changes in the vertical flow profiles in observation boreholes. However, inversion of flow measurements to derive flowpath geometry and connectivity and to characterize their hydraulic properties is still a subject of research. In this study, we propose a framework for cross-borehole flowmeter test interpretation that is based on a two-scale conceptual model: discrete fractures at the borehole scale and zones of interconnected fractures at the aquifer scale. We propose that the two problems may be solved independently. The first inverse problem consists of estimating the hydraulic head variations that drive the transient borehole flow observed in the cross-borehole flowmeter experiments. The second inverse problem is related to estimating the geometry and hydraulic properties of large-scale flowpaths in the region between pumping and observation wells that are compatible with the head variations deduced from the first problem. To solve the borehole-scale problem, we treat the transient flow data as a series of quasi-steady flow conditions and solve for the hydraulic head changes in individual fractures required to produce these data. The consistency of the method is verified using field experiments performed in a fractured-rock aquifer.     

On the acceleration‐based adaptive inverse control of shaking tables

Accurate reproduction of time series with diverse frequency characteristics is a central issue in structural testing. This is true not only for simple experimental tests performed by reaction walls or shaking tables but also for more sophisticated ones, such as hybrid testing. Especially in the latter case, where actual feedback from an ongoing test is used in the calculation of the next excitation value, any possible mismatch may be fatal for both the validity of the test and the safety. The objective of this study is to propose a framework for the adaptive inverse control of shaking tables, which succeeds in this matching to a certain degree. By formulating a critical set of design specifications that correspond to safety, implementation, robustness and ease of use, the conducted research results in a design that is based on a modified version of the filtered‐X algorithm with very competitive features. These are the following: (i) default operation in hard real‐time and acceleration mode; (ii) very low hardware requirements; (iii) effective cancelation of the shaking table's dynamics; and (iv) robustness against specimen dynamics. For its practical evaluation, the method is applied to shaking table waveform replication tests under the installation of an approximately linear specimen of sufficiently high mass and complex geometry. The results are promising and suggest further research toward this field, especially in conjunction with hybrid testing, as the method retains certain global applicability attributes and it can be easily extended to other transfer systems, apart from shaking tables. Copyright © 2014 John Wiley & Sons, Ltd.     

Enforcing smoothness and assessing uncertainty in non-linear one-dimensional prestack seismic inversion

Estimation of elastic properties of rock formations from surface seismic amplitude measurements remains a subject of interest for the exploration and development of hydrocarbon reservoirs. This paper develops a global inversion technique to estimate and appraise 1D distributions of compressional‐wave velocity, shear‐wave velocity and bulk density, from normal‐moveout‐corrected PP prestack surface seismic amplitude measurements. Specific objectives are: (a) to evaluate the efficiency of the minimization algorithm (b) to appraise the impact of various data misfit functions, and (c) to assess the effect of the degree and type of smoothness criterion enforced by the inversion. Numerical experiments show that very fast simulated annealing is the most efficient minimization technique among alternative approaches considered for global inversion. It is also found that an adequate choice of data misfit function is necessary for a reliable and efficient match of noisy and sparse seismic amplitude measurements. Several procedures are considered to enforce smoothness of the estimated 1D distributions of elastic parameters, including predefined quadratic measures of length, flatness and roughness. Based on the general analysis of global inversion techniques, we introduce a new stochastic inversion algorithm that initializes the search for the minimum with constrained random distributions of elastic parameters and enforces predefined autocorrelation functions (semivariograms). This strategy readily lends itself to the assessment of model uncertainty. The new global inversion algorithm is successfully tested on noisy synthetic amplitude data. Moreover, we present a feasibility analysis of the resolution and uncertainty of prestack seismic amplitude data to infer 1D distributions of elastic parameters measured with wireline logs in the deepwater Gulf of Mexico. The new global inversion algorithm is computationally more efficient than the alternative global inversion procedures considered here.     

Toward a full multiscale approach to interpret potential fields

The way potential fields convey source information depends on the scale at which the field is analysed. In this sense a multiscale analysis is a useful method to study potential fields particularly when the main field contributions are caused by sources with different depths and extents. Our multiscale approach is built with a stable transformation, such as depth from extreme points. Its stability results from mixing, in a single operator, the wavenumber low‐pass behaviour of the upward continuation transformation of the field with the enhancement high‐pass properties of n‐order derivative transformations. So, the complex reciprocal interference of several field components may be efficiently faced at several scales of the analysis and the depth to the sources may be estimated together with the homogeneity degrees of the field. In order to estimate the source boundaries we use another multiscale method, the multiscale derivative analysis, which utilizes a generalized concept of horizontal derivative and produces a set of boundary maps at different scales. We show through synthetic examples and application to the gravity field of Southern Italy that this multiscale behaviour makes this technique quite different from other source boundary estimators. The main result obtained by integrating multiscale derivative analysis with depth from extreme points is the retrieval of rather effective information of the field sources (horizontal boundaries, depth, structural index). This interpretative approach has been used along a specific transect for the analysis of the Bouguer anomaly field of Southern Apennines. It was set at such scales, so to emphasize either regional or local features along the transect. Two different classes of sources were individuated. The first one includes a broad, deep source with lateral size of 45∼50 km, at a depth of 13 km and having a 0.5 structural index. The second class includes several narrower sources located at shallowest depths, ranging from 3–6 km, with lateral size not larger than 5 km and structural indexes ranging from 1–1.5. Within a large‐scale geological framework, these results could help to outline the mean structural features at crustal depths.     

Modelling mid‐channel bars in meandering channels

An analytical modelling framework is proposed to reproduce the frequently observed but poorly studied occurrence of mid‐channel bars in meandering channels. Mid‐channel bars occur in meanders and may characterize transitional morphologies between pure meandering and braided rivers. Based on existing field and experimental observations, we propose that two different mechanisms can generate central topographical patterns in meanders. A former mechanism (‘width‐forced’) is related to spatial width oscillations which determine a laterally symmetrical bed shear stress pattern that promotes mid‐channel bars. A second mechanism (‘curvature‐forced’) can take place also in curvilinear equiwidth streams since also longitudinal variations of channel curvature can produce laterally symmetrical alterations of the sediment transport capacity. A perturbation approach is employed to model both mechanisms within a common framework, allowing reproduction, at least qualitatively, of several observed features. While width‐forced mid‐channel bars are a symmetric linear altimetric response, to reproduce curvature‐forced mid‐channel bars requires modelling nonlinear flow‐bed topography interactions at the second order of the perturbation expansion. Hypotheses on how these mechanisms operate are further discussed through an application to field cases. The amplitude of the nonlinear response can be relevant compared to that of the point bar in equiwidth meanders and the location of mid‐channel bars seldom coincides with bend apexes, mainly depending upon the intrinsic meander wavelength. Central bars tend to symmetrically divert the flow against the two banks, a process which is proposed as a possible cause of cross‐sectional overwidening, along with the asymmetry between the rates of bank erosion and of the opposite bank accretion. The outcomes of this first modelling step on the subject allow discussion of the mutual feedback processes that characterize interactions between mid‐channel bars and width variations in river meanders. Copyright © 2010 John Wiley & Sons, Ltd.     

The application of spatial derivatives to non‐potential field data interpretation

Source/body edge detection is a common feature in the processing and interpretation of potential field data sets. A wide range of spatial derivatives is available to enhance the information contained in the basic data. Here the ability of these procedures to assist with the mapping interpretation of non‐potential field data is considered. The study uses airborne electromagnetic (conductivity) data but also provides a general context for other conductivity/resistivity data, provided the non‐potential field nature of active and thus spatially‐focused, measurements is acknowledged. The study discusses and demonstrates the application of a range of common spatial derivative procedures, including the analytic signal and upward continuation, to both magnetic and conductivity data. The ability of the tilt derivative to provide enhanced mapping of conductivity data is considered in detail. Tilt and its associated functions are formed by taking combinations of vertical and horizontal derivatives of the data set. Theoretical forward modelling studies are first carried out to assess the performance of the tilt derivative in relation to the detection and definition of concealed conductivity structure. The tilt derivative embodies automatic gain control that normalizes the detection and definition of both weak and strong conductivity gradients across an appropriate subsurface depth range. The use of high‐order spatial derivatives inevitably results in a degree of noise (cultural perturbation) amplification that is survey and technique specific. Both of these aspects are considered using practical case studies of jointly obtained magnetic and conductivity data at a variety of spatial scales.     

Water condition control of in situ soil water repellency: an observational study from a hillslope in a Japanese humid‐temperate forest

To evaluate the relationship between the moisture conditions and the water repellency of soil surfaces in situ, we periodically conducted field surveys for more than a year in a humid‐temperate forest in Japan. Measurements were made in four plots with varying soil physicochemical properties and under different topographic conditions across a hillslope. Each plot contained permanent quadrats with measurement points in a grid pattern. At each point, we measured the volumetric water content at 0‐ to 5‐cm depths and the water repellency at soil surfaces approximately twice a month. The repeated measurements enabled us to estimate the critical water content (CWC) below which soils repelled water at each point. We defined the representative CWC (RCWC) of a plot as the median of all CWCs in a plot and estimated the representative critical water potential (RCWP) on the basis of the RCWC using the water retention curve. The RCWC values differed among plots, but the corresponding RCWP values were similar (pF = 3.5–3.9). The relationship of the areal fraction showing water repellency against soil water potentials was similar across plots, but the relationship differed among plots against the soil moisture content. These results suggest that soil water potential is more indicative of the spatial occurrence of water repellency than moisture content on a hillslope where soil physicochemical properties vary. Plots located on ridge crests frequently exhibited lower water potentials and showed a higher areal fraction of water repellency, implying a greater chance of generating surface runoff by rainfall events. Copyright © 2011 John Wiley & Sons, Ltd.