Improved understanding of spatio‐temporal controls on regional scale groundwater flooding using hydrograph analysis and impulse response functions

Controls on the spatio‐temporal extent of groundwater flooding are poorly understood, despite the long duration of groundwater flood events and distinct social and economic impacts. We developed a novel approach using statistical analysis of groundwater level hydrographs and impulse response functions (IRFs) and applied it to the 2013/2014 Chalk groundwater flooding in the English Lowlands. We proposed a standardized index of groundwater flooding which we calculated for monthly groundwater levels for 26 boreholes in the Chalk. We grouped these standardized series using k‐means cluster analysis and cross‐correlated the cluster centroids with the Standardized Precipitation Index accumulated over time intervals between 1 and 60 months. This analysis reveals 2 spatially coherent groups of standardized hydrographs that responded to precipitation over different timescales. We estimated IRF models of the groundwater level response to effective precipitation for 3 boreholes in each group. The IRF models corroborate the Standardized Precipitation Index analysis showing different response functions between the groups. We applied identical effective precipitation inputs to each of the IRF models and observed differences between the hydrographs from each group. It is suggested this is due to the hydrogeological properties of the Chalk and of overlying relatively low permeability superficial deposits (recent unconsolidated sediments overlying the bedrock, such as clays and tills), which are extensive over 1 of the groups. The overarching controls on groundwater flood response are concluded to be a complex combination of antecedent conditions, rainfall, and catchment hydrogeological properties. These controls should be taken into consideration when anticipating and managing future groundwater flood events. The approach presented is generic and parsimonious and can be easily applied where sufficient groundwater level and rainfall data are available.     

Comprehensive evaluation of regional flood frequency analysis by L- and LH-moments. I. A re-visit to regional homogeneity

As part I of a sequence of two papers, previously developed L-moments by Hosking (J R Stat Soc Ser B Methodol 52(2):105–124, 1990), and the LH-moments by Wang (Water Resour Res 33(12):2841–2848, 1997) are re-visited. New relationships are developed for regional homogeneity analysis by the LH-moments, and further establishment of regional homogeneity is investigated. Previous works of Hosking (J R Stat Soc Ser B Methodol 52(2):105–124, 1990) and Wang (Water Resour Res 33(12):2841–2848, 1997) on L-moments and LH-moments for the generalized extreme value (GEV) distribution are extended to the generalized Pareto (GPA) and the generalized logistic (GLO) distributions. The Karkhe watershed, located in western Iran is used as a case study area. Regional homogeneity was investigated by first assuming the entire study area as one regional cluster. Then the entire study area was designated “homogeneous” by the L-moments (L); and was designated “heterogeneous” by all four levels of the LH-moments (L₁ to L₄). The k-means method was used to investigate the case of two regional clusters. All levels of the L- and LH-moments designated the upper watershed (region A), “homogeneous”, and the lower watershed (region B) “possibly-homogeneous”. The L₃ level of the GPA and the L₄ level of the GLO were selected for regions A and B, respectively. Wang (Water Resour Res 33(12):2841–2848, 1997) identified a reversing trend in improved performance of the GEV distribution at the LH-moments level of L₃ (during the goodness-of-fit test). Similar results were also obtained in this research for the GEV distribution. However, for the case of the GPA distribution the reversing trend started at L₄ for region A; and at L₂ for region B. As for the case of the GLO, an improved performance was observed for all levels (moving from L to L₄); for both regions.     

Role of reservoir regulation and groundwater feedback in a simulated ground-soil-vegetation continuum: A long-term regional scale analysis

The regional terrestrial water cycle is strongly altered by human activities. Among them, reservoir regulation is a way to spatially and temporally allocate water resources in a basin for multi-purposes. However, it is still not sufficiently understood how reservoir regulation modifies the regional terrestrial- and subsequently, the atmospheric water cycle. To address this question, the representation of reservoir regulation into the terrestrial component of fully coupled regional Earth system models is required. In this study, an existing process-based reservoir network module is implemented into NOAH-HMS, that is, the terrestrial component of an atmospheric–hydrologic modelling system, namely, the WRF-HMS. It allows to quantitatively differentiate role of reservoir regulation and of groundwater feedback in a simulated ground-soil-vegetation continuum. Our study focuses on the Poyang Lake basin, where the largest freshwater lake of China and reservoirs of different sizes are located. As compared to streamflow observations, the newly extended NOAH-HMS slightly improves the streamflow and streamflow duration curves simulation for the Poyang Lake basin for the period 1979–1986. The inclusion of reservoir regulation leads to major changes in the simulated groundwater recharges and evaporation from reservoirs at local scale, but has minor effects on the simulated soil moisture and surface runoff at basin scale. The performed groundwater feedback sensitivity analysis shows that the strength of the groundwater feedback is not altered by the consideration of reservoir regulation. Furthermore, both reservoir regulation and groundwater feedback modify the partitioning of the simulated evapotranspiration, thus affecting the atmospheric water cycle in the Poyang Lake region. This finding motivates future research with our extended fully coupled atmospheric–hydrologic modelling system by the community.     

Assessment of chemical water types and their spatial variation using multi-stage cluster analysis, Queensland, Australia

A multivariate assessment has been adapted to the classification of a large, irregular dataset of approximately 34,000 surface water samples accumulated over more than 30 years. A two-stage K-means clustering method was designed to analyse chemical data in the form of percentages of major ions (Na, Mg, Ca, Cl, HCO₃ and SO₄); the first stage of clustering produced 347 groups, which were then re-clustered to generate the final nine water types. The analysis enabled the definition of provinces of water composition and highlighted natural processes influencing the surface water chemistry. On a statewide basis, sodium is the dominant cation and around 50% at all stream flows, while proportions of calcium and magnesium are about equal. Chloride and bicarbonate constitute the bulk of anions present, while sulfate occurs occasionally and tends to be localised. On a global basis, Queensland surface waters are relatively high in sodium, chloride and magnesium, and low in calcium and sulfate. It was also found that the geographical location has a greater impact on major ion ratios than does the stage of stream flow.

The regional chemical trends are consistent with geology and climate. Streams in northeast Queensland, with short, steep catchments and high rainfall, yield low salinity, sodium-dominated water; this is also the case for sandy southern coastal catchments. Both also reflect an oceanic influence. The proportions of sodium and chloride decrease westward; streams draining the western side of the Great Dividing Range or flowing into the Gulf of Carpentaria have low salinity but relatively hard water. Streams in western Queensland are higher in calcium and bicarbonate. In the large catchments flowing from Queensland into central Australia, the water composition is highly variable, commonly with elevated sulfate. Also in Queensland, there are several other clearly definable water provinces such as the high magnesium waters of basaltic areas.

The findings of this study confirm that the application of such analytical methods can provide a useful assessment of controls over water composition and support management at regional level; the approach used is shown and are applicable to large, disparate datasets.     

Sensitivity analysis in groundwater vulnerability assessment based on GIS in the Mahdia-Ksour Essaf aquifer,Tunisia: a validation study

Abstract

The assessment of groundwater vulnerability to pollution has proved to be an effective tool for water resource management, especially in arid and semi-arid regions like Mahdia and Ksour Essaf. The main objective of this study is to assess the aquifer vulnerability by applying the DRASTIC method as well as using sensitivity analysis to evaluate the effect of each DRASTIC parameter on the final vulnerability map. An additional objective is to demonstrate the role of the GIS techniques in the vulnerability assessment. The DRASTIC method assigns a high vulnerability to the coast of the Mahdia-Ksour Essaf. The lowest values are observed in the southern part of the study area. A sensitivity analysis applied in this study suggests that net recharge, aquifer media and depth of groundwater are the key factors determining vulnerability. The model is validated with groundwater quality data and the results have shown strong relationships between modified DRASTIC Vulnerability Index and nitrate and chloride concentrations.

Citation Saidi, S., Bouri, S. & Ben Dhia, H. (2011) Sensitivity analysis in groundwater vulnerability assessment based on GIS in the Mahdia-Ksour Essaf aquifer, Tunisia: a validation study. Hydrol. Sci. J. 56(2), 288–304.     

T‐NET,a dynamic model for simulating daily stream temperature at the regional scale based on a network topology

Currently, the distribution areas of aquatic species are studied by using air temperature as a proxy of water temperature, which is not available at a regional scale. To simulate water temperature at a regional scale, a physically based model using the equilibrium temperature concept and including upstream‐downstream propagation of the thermal signal is proposed. This model, called Temperature‐NETwork (T‐NET), is based on a hydrographical network topology and was tested at the Loire basin scale (10⁵ km²). The T‐NET model obtained a mean root mean square error of 1.6 °C at a daily time step on the basis of 128 water temperature stations (2008–2012). The model obtained excellent performance at stations located on small and medium rivers (distance from headwater <100 km) that are strongly influenced by headwater conditions (median root mean square error of 1.8 °C). The shading factor and the headwater temperature were the most important variables on the mean simulated temperature, while the river discharge influenced the daily temperature variation and diurnal amplitude. The T‐NET model simulates specific events, such as temperature of the Loire during the floods of June 1992 and the thermal regime response of streams during the heatwave of August 2003, much more efficiently than a simple point‐scale heat balance model. The T‐NET model is very consistent at a regional scale and could easily be transposed to changing forcing conditions and to other catchments. Copyright © 2016 John Wiley & Sons, Ltd.     

核环吊缩尺模型抗震性能动力时程分析

为研究核环吊在地震作用下的瞬态响应机理,完成了核环吊缩尺模型的抗震实验,并采用ABAQUS时程分析法研究了缩尺模型的动力响应。首先介绍了核环吊缩尺模型实验的设计方法和准则;其次通过与谱分析的对比分析说明了进行动力时程分析的必要性;最后介绍了ABAQUS中时程分析法的具体方法和步骤。仿真与实验结果的对比分析表明:ABAQUS的Slide-Plane连接单元可准确地模拟大车与小车不跳轨情况下的轮轨约束,验证了关键参数设置方法的准确性,为今后的核环吊产品的抗震设计提供了借鉴和参考。     

Impacts of post‐glacial rebound on landslide spatial distribution at a regional scale in northern Iceland (Skagafjörður)

In de‐glaciated areas, para‐glaciation (i.e. the conditioning of landscapes by prior glaciation) has often been considered a major predisposing factor in landslide occurrence; its consequences have been particularly well identified at a fine scale (especially on bedrock jointing). Hitherto, the relative impacts of para‐glaciation on hillslope dynamics at a regional scale had nevertheless not been quantified statistically. We examine Skagafjörður area (northern Iceland) where landslides are widespread (at least 108 were mapped in an area of c. 3000 km²). We compare the role of para‐glaciation (debuttressing, influence of post‐glacial rebound) with that of classic factors (topography, lithology, etc.) in landslide occurrence and location, using a spatial analysis based on a chi‐square test. On the one hand, the results highlight that landslides are over‐represented in areas where post‐glacial rebound was at its maximum, with a stronger concentration of landslides in the northern part of the fjord. On the other hand, the distribution of landslides did not show any clear relationship with the pattern of glacial debuttressing. Tschuprow coefficient highlights that the influence of post‐glacial rebound on landslide location is higher than the combined influence of slope gradient, curvature or geological structure. This result is supported by our initial evidence for the timing of landslides in the area: most landslides occurred during the first half of the Holocene, and a period of hillslope instability was initiated when the post‐glacial uplift was at its maximum. Finally, the mechanisms that link post‐glacial rebound and landsliding as well as the geomorphic impacts of landslides, are discussed. Copyright © 2013 John Wiley & Sons, Ltd.     

Using multivariate statistical analysis of groundwater major cation and trace element concentrations to evaluate groundwater flow in a regional aquifer

Groundwater samples were collected from 11 springs in Ash Meadows National Wildlife Refuge in southern Nevada and seven springs from Death Valley National Park in eastern California. Concentrations of the major cations (Ca, Mg, Na and K) and 45 trace elements were determined in these groundwater samples. The resultant data were subjected to evaluation via the multivariate statistical technique principal components analysis (PCA), to investigate the chemical relationships between the Ash Meadows and Death Valley spring waters, to evaluate whether the results of the PCA support those of previous hydrogeological and isotopic studies and to determine if PCA can be used to help delineate potential groundwater flow patterns based on the chemical compositions of groundwaters. The results of the PCA indicated that groundwaters from the regional Paleozoic carbonate aquifers (all of the Ash Meadows springs and four springs from the Furnace Creek region of Death Valley) exhibited strong statistical associations, whereas other Death Valley groundwaters were chemically different. The results of the PCA support earlier studies, where potentiometric head levels, δ¹⁸O and δD, geological relationships and rare earth element data were used to evaluate groundwater flow, which suggest groundwater flows from Ash Meadows to the Furnace Creek springs in Death Valley. The PCA suggests that Furnace Creek groundwaters are moderately concentrated Ash Meadows groundwater, reflecting longer aquifer residence times for the Furnace Creek groundwaters. Moreover, PCA indicates that groundwater may flow from springs in the region surrounding Scotty's Castle in Death Valley National Park, to a spring discharging on the valley floor. The study indicates that PCA may provide rapid and relatively cost‐effective methods to assess possible groundwater flow regimes in systems that have not been previously investigated. Copyright © 1999 John Wiley & Sons, Ltd.     

Exploring landscape and geologic controls on spatial patterning of streambank groundwater discharge in a mixed land use watershed

Preferential groundwater discharge features along stream corridors are ecologically important at local and stream network scales, yet we lack quantification of the multiscale controls on the spatial patterning of groundwater discharge. Here we identify physical attributes that best explain variation in the presence and lateral extent of preferential groundwater discharges along two 5th order streams, the Housatonic and Farmington Rivers, and 32 1st to 4th order reaches across the Farmington River network. We mapped locations of preferential groundwater discharge exposed along streambanks using handheld thermal infrared cameras paired with high-resolution topographic and land use land cover datasets, surficial soil characteristic maps, and depth-to-bedrock geophysical measurements. The unconfined Housatonic River, MA, USA (12 km) had fewer discharge locations and less lateral extent (41 discharge locations with 38 m of active discharge/km of river) compared to the partially confined Farmington River, CT, USA (26 km; 169 discharge locations with 129 m of active discharge/km of river). Using a moving window analysis, we found along both rivers that discharge was more likely to occur where bank slopes were steeper, floodplain extent was narrower, and degree of confinement was higher. Along the Farmington River, groundwater discharge was more likely to occur where saturated hydraulic conductivity was higher and depth-to-bedrock was shallower. Among the 32 stream reaches surveyed (33.2 km of total stream length) within the Farmington River watershed, preferential discharge was observed in all but two stream reaches, varied from 0 to 25% of lateral extent along stream banks (mean = 6%), and was more likely to occur where stream reach slopes were steep, saturated hydraulic conductivity was high, and watershed urbanization was low. Our results show that, though both surface (e.g., topographic, land use land cover) and subsurface (e.g., soil characteristics, bedrock depth) factors control the prevalence of streambank preferential groundwater discharge, the dominant controls vary across valley settings and stream sizes.     

Evidence of short-term groundwater recharge signal propagation from the Andes to the central Atacama Desert: a singular spectrum analysis approach

The effect of pressure-driven groundwater recharge signal propagation in the Andean-Atacamenian environment is investigated by assessing a record of 15 years of water table fluctuations of an unconfined–confined aquifer system. Based on a singular spectrum analysis of water table time series, it is shown that, in the given case, groundwater levels in the central Atacama Desert are hydraulically controlled by two distant recharge areas associated with the Andes. The maximum observed range of the pressure signal propagation is ~50 km over an elevation difference of more than 3000 m at a lag of ~25 months. Several findings indicate that an often-cited study misinterpreted a water level rise at the same site as an in-situ alluvial fan recharge. Thus, the effect’s impact on groundwater dynamics in complex aquifer systems can easily be overlooked. Singular spectrum analysis could be of use to investigate pressure effects at hydrologically comparable sites.     

Impacts of nutrient management and decrease in paddy field area on groundwater nitrate concentration: a case study at the Nasunogahara alluvial fan,Tochigi Prefecture,Japan

In order to help evaluate the trends in the NO₃-N concentration in groundwater with a view to preventing further degradation in water quality in the future, a distributed groundwater quality model was constructed for the Nasunogahara basin. The best fit for the groundwater table elevations by the flow component of the model was achieved with average mean absolute errors (MAEs) of 0·92 m for the calibration period and 0·83 m for the validation period. Moreover, the best fit for the NO₃-N concentration by the water quality component was achieved with average mean relative errors (MREs) of 29·8% for the calibration period and 30·3% for the validation period. After developing a robust model, various change scenarios were tested; specifically, the effects of effluent load control and a decrease in paddy field area on the NO₃-N concentration in groundwater were predicted. The most intensively farmed area contributed about 40% of the total effluent load because of livestock farming in the basin. When the effluent load from this area was decreased by 50%, the average NO₃-N concentrations at sites S1, S2 and S3 were reduced by about 15%; however, the average concentrations at S4 and S5 were reduced by only 1%. Furthermore, when the total effluent load from the concentrated livestock area was removed completely, the average groundwater NO₃-N concentrations at S1, S2 and S3 were reduced by about 30% as compared with the original calculated results. In contrast, decreasing the area of the paddy fields in the basin did not greatly influence the groundwater NO₃-N concentration. In the case of a 70% reduction in paddy field area, average NO₃-N concentrations increased by about 7% at S1, S2 and S3. Copyright © 2008 John Wiley & Sons, Ltd.     

A numerical coupling scheme for nonlinear time history analysis of buildings on a regional scale considering site‐city interaction effects

Seismic damage simulation of buildings on a regional scale is important for loss estimation and disaster mitigation of cities. However, the interaction among densely distributed buildings in a city and the site, ie, the “site‐city interaction (SCI) effects,” is often neglected in most regional simulations. Yet, many studies have found that the SCI effects are very important in regional simulations containing a large number of tall buildings and underground structures. Therefore, this work proposed a numerical coupling scheme for nonlinear time history analysis of buildings on a regional scale considering the SCI effects. In this study, multiple‐degree‐of‐freedom models are used to represent different buildings above the ground, while an open source spectral element program, SPEED, is used for simulating wave propagation in underlying soil layers. The proposed numerical scheme is firstly validated through a shaking table test. Then, a detailed discussion on the SCI effects in a 3D basin is performed. Finally, a nonlinear time history analysis of buildings on a regional scale is performed using the Tsinghua University campus in Beijing as a case study. The Tsinghua University campus case results show that the SCI effects will reduce the seismic responses of most buildings. However, some buildings will suffer much more severe damage when the SCI effects are considered, which may depend on the input motions, site characteristics, and building configurations.     

Using geostatistics and spectral analysis to study spatial patterns in the topography of southeastern Washington State,U.S.A.

Topography for four areas in the Palouse region of southeastern Washington State having different patterns and encompassing areas of about 900 ha each were quantitatively compared and described using two-dimensional semi-variograms and periodograms. The four areas studied were from the Garfield, La Crosse, Thornton, and Wilcox quadrangle topographic maps. Semivariance of elevation residuals were modelled using a combination of spherical, periodic, or linear semivariogram models. The range of the spherical component was interpreted as a relatively short-range scale of correlation which was not periodic. For each of the study areas a model was developed to describe the variation in range with orientation. Values for the range from this model reached maximum values of from about 700 m to 800 m at an orientation of from 35° to 55° (approximately northeast) in each study area. This orientation was interpreted as the dominant direction of non-periodic small-scale landscape features. The wavelength and amplitude of the periodic semivariogram component, which were highly correlated, were interpreted in terms of parallel northeast-trending ridges having relatively long scales of periodic correlation varying in size from about 1350 m to 2100 m. This attempt to identify the dominant orientation of periodic landscape features using models for the wavelength and amplitude was, however, not completely definitive. Two-dimensional spectral analysis provided significantly more detail concerning orientation and wavelengths of the periodic topographic patterns than the semivariogram analysis. In the Garfield study area, spectral analysis identified north-trending ridge systems separated by a wavelength of 1494 m and northeast-trending systems with a wavelength of 747 m. In the La Crosse study area, both north and northeast-trending patterns were identified having periodic spacings of 980 and 735 m. North-trending ridges separated by wavelengths of 996 m and northeast-trending ridges separated by wavelengths of 747 or 996 m were the predominant periodic features of topography in the Thornton study area. In the Wilcox study area, northeast-trending ridge systems separated by wavelengths of 373, 498, or 996 m were detected. A comparison of the results from geostatistical and spectral analysis of these complex topographic surfaces shows that each approach had significant strengths and weaknesses. Two-dimensional analysis with semivariograms was the only method which could be used for identifying the correlation scale and orientation of relatively small non-periodic landscape features. Two-dimensional spectral analysis accurately identified the predominant orientation of relatively large periodic features in topography, whereas semivariogram analysis was somewhat inconclusive. Also, semivariograms were generally unable to detect the presence of multiple or harmonic periodicities operating at different wavelengths along a given orientation.     

Stochastic analysis of immiscible displacement of the fluids with arbitrary viscosities and its dependence on support scale of hydrological data

Stochastic analysis is commonly used to address uncertainty in the modeling of flow and transport in porous media. In the stochastic approach, the properties of porous media are treated as random functions with statistics obtained from field measurements. Several studies indicate that hydrological properties depend on the scale of measurements or support scales, but most stochastic analysis does not address the effects of support scale on stochastic predictions of subsurface processes. In this work we propose a new approach to study the scale dependence of stochastic predictions. We present a stochastic analysis of immiscible fluid–fluid displacement in randomly heterogeneous porous media. While existing solutions are applicable only to systems in which the viscosity of one phase is negligible compare with the viscosity of the other (water–air systems for example), our solutions can be applied to the immiscible displacement of fluids having arbitrarily viscosities such as NAPL–water and water–oil. Treating intrinsic permeability as a random field with statistics dependant on the permeability support scale (scale of measurements) we obtained, for one-dimensional systems, analytical solutions for the first moments characterizing unbiased predictions (estimates) of system variables, such as the pressure and fluid–fluid interface position, and we also obtained second moments, which characterize the uncertainties associated with such predictions. Next we obtained empirically scale dependent exponential correlation function of the intrinsic permeability that allowed us to study solutions of stochastic equations as a function of the support scale. We found that the first and second moments converge to asymptotic values as the support scale decreases. In our examples, the statistical moments reached asymptotic values for support scale that were approximately 1/10000 of the flow domain size. We show that analytical moment solutions compare well with the results of Monte Carlo simulations for moderately heterogeneous porous media, and that they can be used to study the effects of heterogeneity on the dynamics and stability of immiscible flow.     

Functional evaluation of a simplified scaling method for assessing the spatial variability of soil hydraulic properties at the hillslope scale

Abstract

Mapping soil hydraulic parameters with traditional scaling methods that use laboratory-determined hydraulic characteristics (the LAB method) is not always feasible as it involves expensive, time-consuming and sophisticated measurements on soil samples collected in several locations of the study area. An alternative scaling method (the AP method) has been recently proposed to indirectly retrieve the soil hydraulic properties following the Arya-Paris physico-empirical pedotransfer function, which makes use of particle-size distribution and bulk density values. In this synthetic study we verify the performance of the AP method from a functional perspective, by evaluating the differences in the simulated soil water budget through a Monte Carlo approach. Notwithstanding that the AP method can provide soil hydraulic property patterns with faster experimental procedures and minor costs, we observe significant bias in the predicted spatially-averaged soil water budget due to a poor parametric calibration of the AP method and an imprecise identification of the spatial correlation structure of the AP-estimated scaling factors.

Citation Nasta, P., Romano, N., and Chirico, G.B., 2013. Functional evaluation of a simplified scaling method for assessing the spatial variability of soil hydraulic properties at the hillslope scale. Hydrological Sciences Journal, 58 (5), 1059–1071.     

Estimation of soil loss at a regional scale based on plot measurements—some critical considerations

The paper is prompted by apparent deficiencies in the design of plot studies in regional erosion surveys. The principal shortcomings of observational erosion research have been poor sampling design and inadequate analyses of data. The paper identifies various sources of bias which must be taken into account before plot data can be extrapolated to land units in a regional survey. Judging from soil loss data of a case-study in the Ardèche rangelands one may conclude that even accurate plot measurements can still be rather a rough basis for regional erosion assessment. Finally, the paper highlights strategies that might be used to improve erosion sampling.     

Landscape‐scale controls on the spatial distribution of caesium 137: A study based on an airborne geophysical survey across Northern Ireland

The spatial distribution of ¹³⁷Cs across the landscape and the processes controlling its redistribution are of interest because (i) ¹³⁷Cs has been widely used to quantify the movement of soil and sediments and (ii) substantial fallout of ¹³⁷Cs after the Chernobyl accident has led to contamination of foodstuffs in some places. A high‐resolution airborne geophysical radiometric survey of Northern Ireland has provided an opportunity to study the distribution and possible redistribution of ¹³⁷Cs. The ¹³⁷Cs activity (recorded at 1·2 million points) is distributed in a series of bands oriented approximately 160° and 115° clockwise from north. Geostatistical analysis of the data shows a strong, short‐range structure (correlation ranges between 0·6 and 8 km) in ¹³⁷Cs activity across the vast majority of the region; the spatial distribution shows association with a published, coarse‐scale depositional pattern of ¹³⁷Cs from Chernobyl. Two indices of land form derived from a digital elevation model, namely compound topographic index and the length–slope factor of the Revised Universal Soil Loss Equation, account for only 3% of the variance in ¹³⁷Cs activity. In contrast, soil type and land cover in combination (including their interaction) account for 20% of the variance. In areas that received moderate fallout from Chernobyl, soil type alone accounts for a substantial proportion of the spatially correlated ¹³⁷Cs activity. We attribute this to each soil type having a fairly uniform radiocaesium interception potential that differs from those of other soil types and that this potential controls the vertical migration of ¹³⁷Cs. Over the granitic Mourne Mountains there is a strong spatial cross‐correlation between ¹³⁷Cs activity and airborne estimates of soil potassium, suggesting that the latter provides a measure of the soil's radiocaesium interception potential; this is probably dominated by the quantity of the mineral illite. Copyright © 2010 John Wiley & Sons, Ltd. and British Geological Survey     

Assessing the impact of the hydraulic properties of a crusted soil on overland flow modelling at the field scale

Soil surface crusts are widely reported to favour Hortonian runoff, but are not explicitly represented in most rainfall‐runoff models. The aim of this paper is to assess the impact of soil surface crusts on infiltration and runoff modelling at two spatial scales, i.e. the local scale and the plot scale. At the local scale, two separate single ring infiltration experiments are undertaken. The first is performed on the undisturbed soil, whereas the second is done after removal of the soil surface crust. The HYDRUS 2D two‐dimensional vertical infiltration model is then used in an inverse modelling approach, first to estimate the soil hydraulic properties of the crust and the subsoil, and then the effective hydraulic properties of the soil represented as a single uniform layer. The results show that the crust hydraulic conductivity is 10 times lower than that of the subsoil, thus illustrating the limiting role the crust has on infiltration. Moving up to the plot scale, a rainfall‐runoff model coupling the Richards equation to a transfer function is used to simulate Hortonian overland flow hydrographs. The previously calculated hydraulic properties are used, and a comparison is undertaken between a single‐layer and a double‐layer representation of the crusted soil. The results of the rainfall‐runoff model show that the soil hydraulic properties calculated at the local scale give acceptable results when used to model runoff at the plot scale directly, without any numerical calibration. Also, at the plot scale, no clear improvement of the results can be seen when using a double‐layer representation of the soil in comparison with a single homogeneous layer. This is due to the hydrological characteristics of Hortonian runoff, which is triggered by a rainfall intensity exceeding the saturated hydraulic conductivity of the soil surface. Consequently, the rainfall‐runoff model is more sensitive to rainfall than to the subsoil's hydrodynamic properties. Therefore, the use of a double‐layer soil model to represent runoff on a crusted soil does not seem necessary, as the increase of precision in the soil discretization is not justified by a better performance of the model. Copyright © 2005 John Wiley & Sons, Ltd.