Soil hydraulic conductivities of latosols under pasture,forest and teak in Rondonia,Brazil

We investigated the changes of saturated hydraulic conductivity, K_sat , with depth of latosols developed on Precambrian basement rocks under primary rainforest, pasture and teak. In all cases, K_sat decreased with depth, with most of the decrease occurring between the surface and a depth of 30 cm. In conjunction with prevailing rainfall intensities and frequencies, this anisotropy supports a pronounced lateral component of hillslope flow paths, and also of overland flow under pasture. Our results are at variance with data from other latosols where K_sat tends to increase with depth, and hence suggest that considerable restraint is needed in generalization and extrapolation until results from a co‐ordinated effort at hydrology‐oriented data collection become available. Copyright © 1999 John Wiley & Sons, Ltd.     

Spatial variability of hydraulic conductivity and bulk density along a blanket peatland hillslope

This article presents the results of a field investigation of saturated hydraulic conductivity K_sat and bulk density (ρ_bd) in an Atlantic blanket bog in the southwest of Ireland. Starting at a peatland stream and moving along an uphill transect toward the peatland interior, ρ_bd and K_sat were examined at regular intervals. Saturated horizontal hydraulic conductivity (Kh_sat) and vertical (Kv_sat) was estimated at two depths: 10–20 and 30–40 cm below the peat surface, whereas ρ_bd was estimated for the full profile. We consider two separate zones, one a riparian zone extending 10 m from the stream and a second zone in the bog interior. We found that the K_sat was higher (~10^–5 m s^–1) in the bog interior than that in the riparian zone (~10^–6 m s^–1), whereas the converse applied to bulk density, with lowest density (~0.055 g cm^–3) at the interior and highest (~0.11 g cm^–3) at the riparian zone. In general, we found Kh_sat to be approximately twice the Kv_sat. These results support the idea that the lower K_sat at the margins control the hydrology of blanket peatlands. It is therefore important that the spatial variation of these two key properties be accommodated in hydrological models if the correct rainfall runoff characteristics are to be correctly modelled. Stream flow analysis over 3 years at the peatland catchment outlet showed that the stream runoff was composed of 8% base flow and 92% flood flow, suggesting that this blanket peatland is a source rather than a sink for floodwaters. Copyright © 2011 John Wiley & Sons, Ltd.     

The soil hydrologic response to forest regrowth: a case study from southwestern Amazonia

As a large and dynamic land‐use category, tropical secondary forests may affect climate, soils, and hydrology in a manner different from primary forests or agricultural areas. We investigated the saturated hydraulic conductivity K_sat of a Kandiudult under different land uses in Rondonia, Brazil. We measured K_sat at four depths (12·5, 20, 30 and 50 cm) under (a) primary forest, (b) a former banana–cacao plantation (SF1), and (c) an abandoned pasture (SF2). At 12·5 cm, all three land uses differ significantly (α = 0·1), but not at the 20 and 30 cm depths. At 50 cm, K_sat was significantly greater in the former pasture than in other land uses. Lateral subsurface flow is expected during intense rainfall (about 30 times per year) at 30 cm depth in SF1 and at 50 cm depth in the forest, whereas the relatively low permeability at shallow 12·5 cm in the SF2 may result not only in lateral subsurface flow, but also saturation overland flow. For modelling purposes, recovering systems seem to have K_sat values distinct from primary forest at shallow depths, whereas at deeper layers (>20 cm) they may be considered similar to forests. Copyright © 2002 John Wiley & Sons, Ltd.     

Hillslope hydrology in tropical rainforest steeplands in Brunei

Many remaining areas of tropical rainforest in south‐east Asia are located on landscapes dominated by deep valleys and very steep slopes. Now that logging activities are extending into these steeplands, it is essential to understand how the natural rainforest system behaves if any kind of realistic assessment of the effects of such disturbance is to be made. This paper examines the hydrological behaviour of an undisturbed rainforest system on steep topography in the Temburong District of Brunei, north‐west Borneo. The physical and hydrological properties of the regolith material are generally typical of tropical residual soils. The regolith has a clay texture and a low dry bulk density beneath a superficial litter/organic horizon. The infiltration capacity of the surface soil was several hundred mm h⁻¹. That of the exposed mineral subsoil was an order of magnitude less, similar to the saturated hydraulic conductivity (K_sat) of around 180 mm h⁻¹ at a depth of 150 cm. There was no indication that K_sat reduced with depth except very near the bedrock interface. Soil tensions were measured using a two‐dimensional array of tensiometers on a 30° slope. During dry season conditions, infiltrating rain‐water contributes to soil moisture, and drying of the soil is dominated by transpiration losses. During wet season conditions, perched water tables quickly develop during heavy rainfall, giving rise to the rapid production of return flow in ephemeral channels. No infiltration excess or saturation overland flow was observed on hillslopes away from channel margins. Subsurface storm flow combined with return flow produce stream flow hydrographs with high peak discharges and very short lag times. Storm event runoff coefficients are estimated to be as high as 40%. It is concluded that the most distinctive feature of the hydrology of this ‘steepland rainforest’ is the extremely ‘flashy’ nature of the catchment runoff regime produced by the combination of thin but very permeable regolith on steep slopes. Copyright © 2000 John Wiley & Sons, Ltd.     

Determining hydraulic properties of a loam soil by alternative infiltrometer techniques

Testing infiltrometer techniques to determine soil hydraulic properties is necessary for specific soils. For a loam soil, the water retention and hydraulic conductivity predicted by the BEST (Beerkan Estimation of Soil Transfer parameters) procedure of soil hydraulic characterization was compared with data collected by more standard laboratory and field techniques. Six infiltrometer techniques were also compared in terms of saturated soil hydraulic conductivity, K_s. BEST yielded water retention values statistically similar to those obtained in the laboratory and K_s values practically coinciding with those determined in the field with the pressure infiltrometer (PI). The unsaturated soil hydraulic conductivity measured with the tension infiltrometer (TI) was reproduced satisfactorily by BEST only close to saturation. BEST, the PI, one‐potential experiments with both the TI and the mini disk infiltrometer (MDI), the simplified falling head (SFH) technique and the bottomless bucket (BB) method yielded statistically similar estimates of K_s, differing at the most by a factor of three. Smaller values were obtained with longer and more soil‐disturbing infiltration runs. Any of the tested infiltration techniques appears usable to obtain the order of magnitude of K_s at the field site, but the BEST, BB and PI data appear more appropriate to characterize the soil at some stage during a rainfall event. Additional investigations on both similar and different soils would allow development of more general procedures to apply infiltrometer techniques for soil hydraulic characterization. Copyright © 2015 John Wiley & Sons, Ltd.     

Geomorphic controls on the physical and hydrologic properties of soils in a valley floor

The distribution of soil hydraulic and physical properties strongly influences runoff processes in landscapes. Although much work has been done to quantify and predict the properties of hillslope soils, far less is known about the distribution of soil properties in valley floors. A technique that links the estimation and distribution of soil hydraulic properties in valleys, with easily identified geomorphic features, was developed along a 2 km length of a valley at Brooks Creek in New South Wales, Australia. Soil physical and hydraulic property data were collected across a set of floodplain and fan features within the valley and analysed statistically to determine if soil properties varied significantly between geomorphic features and stratigraphic layers. The results show that the depth‐averaged saturated hydraulic conductivity, K_s, of the soil varies significantly with landform: fan units have K_g values that are twice that of floodplains and colluvial toeslope deposits have K_s values four times higher than floodplains. Given the notorious variability of K_s values in space, the strong statistical separation of soil properties by landform, backed up by strong separation of soil particle size by landform, suggests a way forward in understanding the distribution of soil properties in valleys and their influence on catchment hydrology. Copyright © 2000 John Wiley & Sons, Ltd.     

Identification of optimal soil hydraulic functions and parameters for predicting soil moisture

Abstract

The accuracy of six combined methods formed by three commonly-used soil hydraulic functions and two methods to determine soil hydraulic parameters based on a soil hydraulic parameter look-up table and soil pedotransfer functions was examined for simulating soil moisture. A novel data analysis and modelling approach was used that eliminated the effects of evapotranspiration so that specific sources of error among the six combined methods could be identified and quantified. By comparing simulated and observed soil moisture at six sites of the USDA Soil Climate Analysis Network, we identified the optimal soil hydraulic functions and parameters for predicting soil moisture. Through sensitivity tests, we also showed that adjusting only the soil saturated hydraulic conductivity, K_s , is insufficient for representing important effects of macropores on soil hydraulic conductivity. Our analysis illustrates that, in general, soil hydraulic conductivity is less sensitive to K_s than to the soil pore-size distribution parameter.

Editor D. Koutsoyiannis; Associate editor D. Hughes

Citation Pan, F., McKane, R.B. and Stieglitz, M., 2012. Identification of optimal soil hydraulic functions and parameters for predicting soil moisture. Hydrological Sciences Journal, 57 (4), 723–737.     

Meta‐analysis of field‐saturated hydraulic conductivity recovery following wildland fire: Applications for hydrologic model parameterization and resilience assessment

Hydrologic recovery after wildfire is critical for restoring the ecosystem services of protecting of human lives and infrastructure from hazards and delivering water supply of sufficient quality and quantity. Recovery of soil‐hydraulic properties, such as field‐saturated hydraulic conductivity (K_fs), is a key factor for assessing the duration of watershed‐scale flash flood and debris flow risks after wildfire. Despite the crucial role of K_fs in parameterizing numerical hydrologic models to predict the magnitude of postwildfire run‐off and erosion, existing quantitative relations to predict K_fs recovery with time since wildfire are lacking. Here, we conduct meta‐analyses of 5 datasets from the literature that measure or estimate K_fs with time since wildfire for longer than 3‐year duration. The meta‐analyses focus on fitting 2 quantitative relations (linear and non‐linear logistic) to explain trends in K_fs temporal recovery. The 2 relations adequately described temporal recovery except for 1 site where macropore flow dominated infiltration and K_fs recovery. This work also suggests that K_fs can have low hydrologic resistance (large postfire changes), and moderate to high hydrologic stability (recovery time relative to disturbance recurrence interval) and resilience (recovery of hydrologic function and provision of ecosystem services). Future K_fs relations could more explicitly incorporate processes such as soil‐water repellency, ground cover and soil structure regeneration, macropore recovery, and vegetation regrowth.     

DHI evaluation by combining rock physics simulation and statistical techniques for fluid identification of Cambrian-to-Cretaceous clastic reservoirs in Pakistan

The use of seismic direct hydrocarbon indicators is very common in exploration and reservoir development to minimise exploration risk and to optimise the location of production wells. DHIs can be enhanced using AVO methods to calculate seismic attributes that approximate relative elastic properties. In this study, we analyse the sensitivity to pore fluid changes of a range of elastic properties by combining rock physics studies and statistical techniques and determine which provide the best basis for DHIs. Gassmann fluid substitution is applied to the well log data and various elastic properties are evaluated by measuring the degree of separation that they achieve between gas sands and wet sands. The method has been applied successfully to well log data from proven reservoirs in three different siliciclastic environments of Cambrian, Jurassic, and Cretaceous ages. We have quantified the sensitivity of various elastic properties such as acoustic and extended elastic (EEI) impedances, elastic moduli (K _sat and K _sat–μ), lambda–mu–rho method (λρ and μρ), P-to-S-wave velocity ratio (V _P/V _S), and Poisson’s ratio (σ) at fully gas/water saturation scenarios. The results are strongly dependent on the local geological settings and our modeling demonstrates that for Cambrian and Cretaceous reservoirs, K _sat–μ, EEI, V _P/V _S, and σ are more sensitive to pore fluids (gas/water). For the Jurassic reservoir, the sensitivity of all elastic and seismic properties to pore fluid reduces due to high overburden pressure and the resultant low porosity. Fluid indicators are evaluated using two metrics: a fluid indicator coefficient based on a Gaussian model and an overlap coefficient which makes no assumptions about a distribution model. This study will provide a potential way to identify gas sand zones in future exploration.     

Sensitivity of soil parameters in unsaturated zone modelling and the relation between effective,laboratory and in situ estimates

Simulation of soil moisture content requires effective soil hydraulic parameters that are valid at the modelling scale. This study investigates how these parameters can be estimated by inverse modelling using soil moisture measurements at 25 locations at three different depths (at the surface, at 30 and 60 cm depth) on an 80 by 20 m hillslope. The study presents two global sensitivity analyses to investigate the sensitivity in simulated soil moisture content of the different hydraulic parameters used in a one‐dimensional unsaturated zone model based on Richards' equation. For estimation of the effective parameters the shuffled complex evolution algorithm is applied. These estimated parameters are compared to their measured laboratory and in situ equivalents. Soil hydraulic functions were estimated in the laboratory on 100 cm³ undisturbed soil cores collected at 115 locations situated in two horizons in three profile pits along the hillslope. Furthermore, in situ field saturated hydraulic conductivity was estimated at 120 locations using single‐ring pressure infiltrometer measurements. The sensitivity analysis of 13 soil physical parameters (saturated hydraulic conductivity (K_s), saturated moisture content (θ_s), residual moisture content (θ_r), inverse of the air‐entry value (α), van Genuchten shape parameter (n), Averjanov shape parameter (N) for both horizons, and depth (d) from surface to B horizon) in a two‐layer single column model showed that the parameter N is the least sensitive parameter. K_s of both horizons, θ_s of the A horizon and d were found to be the most sensitive parameters. Distributions over all locations of the effective parameters and the distributions of the estimated soil physical parameters from the undisturbed soil samples and the single‐ring pressure infiltrometer estimates were found significantly different at a 5% level for all parameters except for α of the A horizon and K_s and θ_s of the B horizon. Different reasons are discussed to explain these large differences. Copyright © 2004 John Wiley & Sons, Ltd.     

Synthesis of soil‐hydraulic properties and infiltration timescales in wildfire‐affected soils

We collected soil‐hydraulic property data from the literature for wildfire‐affected soils, ash, and unburned soils. These data were used to calculate metrics and timescales of hydrologic response related to infiltration and surface runoff generation. Sorptivity (S) and wetting front potential (Ψ_f) were significantly different (lower) in burned soils compared with unburned soils, whereas field‐saturated hydraulic conductivity (K_fs) was not significantly different. The magnitude and duration of the influence of capillarity during infiltration was greatly reduced in burned soils, causing faster ponding times in response to rainfall. Ash had large values of S and K_fs but moderate values of Ψ_f, compared with unburned and burned soils, indicating ash has long ponding times in response to rainfall. The ratio of S²/K_fs was nearly constant (~100 mm) for unburned soils but more variable in burned soils, suggesting that unburned soils have a balance between gravity and capillarity contributions to infiltration that may depend on soil organic matter, whereas in burned soils the gravity contribution to infiltration is greater. Changes in S and K_fs in burned soils act synergistically to reduce infiltration and accelerate and amplify surface runoff generation. Synthesis of these findings identifies three key areas for future research. First, short timescales of capillary influences on infiltration indicate the need for better measurements of infiltration at times less than 1 min to accurately characterize S in burned soils. Second, using parameter values, such as Ψ_f, from unburned areas could produce substantial errors in hydrologic modeling when used without adjustment for wildfire effects, causing parameter compensation and resulting underestimation of K_fs. Third, more thorough measurement campaigns that capture soil‐structural changes, organic matter impacts, quantitative water repellency trends, and soil‐water content along with soil‐hydraulic properties could drive the development of better techniques for numerically simulating infiltration in burned areas.     

Variation in soil saturated hydraulic conductivity along the hillslope of collapsing granite gullies

Saturated hydraulic conductivity (K_s) affects the soil hydrological process and is influenced by many factors that exhibit strong spatial variations. To accurately measure K_s and its scale, spatial variability and relationship with collapsing gullies, we analysed four double-ring infiltrometer diameters in three soil layers during in situ experiments designed to measure K_s in two typical collapsing gullies (three slope sites) in Tongcheng County of China. The results showed that K_s increased with increasing inner ring diameter, but no significant difference existed between inner diameters of 30 and 40 cm. The K_s in red soil layers was higher than that in sandy soil layers, the transition layers had the lowest values. K_s also varied with slope position, gradually decreasing with distance from the gully head. The suggestion is that the spatial variation in K_s is affected not only by the intrinsic soil properties but also by the interaction with the collapsing gully.     

Effect of anisotropy on solute transport in degraded fen peat soils

Peat soils are heterogeneous, anisotropic porous media. Compared to mineral soils, there is still limited understanding of physical and solute transport properties of fen peat soils. In this study, we aimed to explore the effect of soil anisotropy on solute transport in degraded fen peat. Undisturbed soil cores, taken in vertical and horizontal direction, were collected from one drained and one restored fen peatland both in a comparable state of soil degradation. Saturated hydraulic conductivity (K _s) and chemical properties of peat were determined for all soil cores. Miscible displacement experiments were conducted under saturated steady state conditions using potassium bromide as a conservative tracer. The results showed that (1) the K _s in vertical direction (K _sv) was significantly higher than that in horizontal direction (K_sh), indicating that K _s of degraded fen peat behaves anisotropically; (2) pronounced preferential flow occurred in vertical direction with a higher immobile water fraction and a higher pore water velocity; (3) the 5% arrival time (a proxy for the strength of preferential flow) was affected by soil anisotropy as well as study site. A strong correlation was found between 5% arrival time and dispersivity, K _s and mobile water fraction; (4) phosphate release was observed from drained peat only. The impact of soil heterogeneity on phosphate leaching was more pronounced than soil anisotropy. The soil core with the strongest preferential flow released the highest amount of phosphate. We conclude that soil anisotropy is crucial in peatland hydrology but additional research is required to fully understand anisotropy effects on solute transport.     

Using a transient infiltrometric technique for intensively sampling field‐saturated hydraulic conductivity of a clay soil in two runoff plots

The point measurement of soil properties allows to explain and simulate plot scale hydrological processes. An intensive sampling was carried out at the surface of an unsaturated clay soil to measure, on two adjacent plots of 4 × 11 m² and two different dates (May 2007 and February–March 2008), dry soil bulk density, ρ_b, and antecedent soil water content, θ_i, at 88 points. Field‐saturated soil hydraulic conductivity, K_fs, was also measured at 176 points by the transient Simplified Falling Head technique to determine the soil water permeability characteristics at the beginning of a possible rainfall event yielding measurable runoff. The ρ_b values did not differ significantly between the two dates, but wetter soil conditions (by 31%) and lower conductivities (1.95 times) were detected on the second date as compared with the first one. Significantly higher (by a factor of 1.8) K_fs values were obtained with the 0.30‐m‐diameter ring compared with the 0.15‐m‐diameter ring. A high K_fs (> 100 mm h^?1) was generally obtained for low θ_i values (< 0.3 m³m^?3), whereas a high θ_i yielded an increased percentage of low K_fs data (1–100 mm h^?1). The median of K_fs for each plot/sampling date combination was not lower than 600 mm h^?1, and rainfall intensities rarely exceeded 100 mm h^?1 at the site. The occurrence of runoff at the base of the plot needs a substantial reduction of the surface soil permeability characteristics during the event, probably promoted by a higher water content than the one of this investigation (saturation degree = 0.44–0.62) and some soil compaction due to rainfall impact. An intensive soil sampling reduces the risk of an erroneous interpretation of hydrological processes. In an unstable clay soil, changes in K_fs during the event seem to have a noticeable effect on runoff generation, and they should be considered for modeling hydrological processes. Copyright © 2012 John Wiley & Sons, Ltd.     

Development of probability distributions for urban hydrologic model parameters and a Monte Carlo analysis of model sensitivity

This paper proposes an approach to estimating the uncertainty related to EPA Storm Water Management Model model parameters, percentage routed (P_R) and saturated hydraulic conductivity (K_sat), which are used to calculate stormwater runoff volumes. The methodology proposed in this paper addresses uncertainty through the development of probability distributions for urban hydrologic parameters through extensive calibration to observed flow data in the Philadelphia collection system. The established probability distributions are then applied to the Philadelphia Southeast district model through a Monte Carlo approach to estimate the uncertainty in prediction of combined sewer overflow volumes as related to hydrologic model parameter estimation. Understanding urban hydrology is critical to defining urban water resource problems. A variety of land use types within Philadelphia coupled with a history of cut and fill have resulted in a patchwork of urban fill and native soils. The complexity of urban hydrology can make model parameter estimation and defining model uncertainty a difficult task. The development of probability distributions for hydrologic parameters applied through Monte Carlo simulations provided a significant improvement in estimating model uncertainty over traditional model sensitivity analysis. Copyright © 2013 John Wiley & Sons, Ltd.     

A comparison between the single ring pressure infiltrometer and simplified falling head techniques

Testing the relative performances of the single ring pressure infiltrometer (PI) and simplified falling head (SFH) techniques to determine the field saturated soil hydraulic conductivity, K_fs, at the near point scale may help to better establish the usability of these techniques for interpreting and simulating hydrological processes. A sampling of 10 Sicilian sites showed that the measured K_fs was generally higher with the SFH technique than the PI one, with statistically significant differences by a factor varying from 3 to 192, depending on the site. A short experiment with the SFH technique yielded higher K_fs values because a longer experiment with the PI probably promoted short‐term swelling phenomena reducing macroporosity. Moreover, the PI device likely altered the infiltration surface at the beginning of the run, particularly in the less stable soils, where soil particle arrangement may be expected to vary upon wetting. This interpretation was supported by a soil structure stability index, SSI, and also by the hydraulic conductivity data obtained with the tension infiltrometer, i.e. with a practically negligible disturbance of the sampled soil surface. In particular, a statistically significant, increasing relationship with SSI and an unsaturated conductivity greater than the saturated one were only detected for the K_fs data obtained with the PI. The SFH and PI techniques should be expected to yield more similar results in relatively rigid porous media (low percentages of fine particles and structurally stable soils) than in soils that modify appreciably their particle arrangement upon wetting. The simultaneous use of the two techniques may allow to improve K_fs determination in soils that change their hydrodynamic behaviour during a runoff producing rainfall event. Copyright © 2013 John Wiley & Sons, Ltd.     

Hydraulic properties of peat soils along a bulk density gradient—A meta study

Our understanding of hydraulic properties of peat soils is limited compared with that of mineral substrates. In this study, we aimed to deduce possible alterations of hydraulic properties of peat soils following degradation resulting from peat drainage and aeration. A data set of peat hydraulic properties (188 soil water retention curves [SWRCs], 71 unsaturated hydraulic conductivity curves [UHCs], and 256 saturated hydraulic conductivity [K_s] values) was assembled from the literature; the obtained data originated from peat samples with an organic matter (OM) content ranging from 23 to 97 wt% (weight percent; and according variation in bulk density) representing various degrees of peat degradation. The Mualem‐van Genuchten model was employed to describe the SWRCs and UHCs. The results show that the hydraulic parameters of peat soils vary over a wide range confirming the pronounced diversity of peat. Peat decomposition significantly modifies all hydraulic parameters. A bulk density of approximately 0.2 g cm^?3 was identified as a critical threshold point; above and below this value, macroporosity and hydraulic parameters follow different functions with bulk density. Pedotransfer functions based on physical peat properties (e.g., bulk density and soil depth) separately computed for bog and fen peat have significantly lower mean square errors than functions obtained from the complete data set, which indicates that not only the status of peat decomposition but also the peat‐forming plants have a large effect on hydraulic properties. The SWRCs of samples with a bulk density of less than 0.2 g cm^?3 could be grouped into two to five classes for each peat type (botanical composition). The remaining SWRCs originating from samples with a bulk density of larger than 0.2 g cm^?3 could be classified into one group. The Mualem‐van Genuchten parameter values of α can be used to estimate K_s if no K_s data are available. In conclusion, the derived pedotransfer functions provide a solid instrument to derive hydraulic parameter values from easily measurable quantities; however, additional research is required to reduce uncertainty.     

Seasonal variations in infiltrability of moss‐dominated biocrusts on aeolian sand and loess soil in the Chinese Loess Plateau

Biocrust effects on soil infiltration have attracted increasing attention in dryland ecosystems, but their seasonal variations in infiltrability have not yet been well understood. On the Chinese Loess Plateau, soil infiltrability indicated by saturated hydraulic conductivity (K_s) of biocrusts and bare soil, both on aeolian sand and loess soil, was determined by disc infiltrometer in late spring (SPR), midsummer (SUM), and early fall (FAL). Then their correlations with soil biological and physiochemical properties and water repellency index (RI) were analysed. The results showed that the biocrusts significantly decreased K_s both on sand during SPR, SUM, and FAL (by 43%, 66%, and 35%, respectively; P < .05) and on loess (by 42%, 92%, and 10%, respectively; P <.05). As compared with the bare soil, the decreased K_s in the biocrusted surfaces was mostly attributed to the microorganism biomass and also to the increasing content of fine particles and organic matter. Most importantly, both the biocrusts and bare soil exhibited significant (F ≥ 11.89, P ≤ .003) seasonal variations in K_s, but their patterns were quite different. Specifically, the K_s of bare soil gradually decreased from SPR to SUM (32% and 42% for sand and loess, respectively) and FAL (29% and 39%); the K_s of biocrusts also decreased from SPR to SUM (59% and 92%) but then increased in FAL (36% and 588%). Whereas the seasonal variations in K_s of the biocrusts were closely correlated with the seasonal variations in RI, the RI values were not high enough to point at hydrophobicity. Instead of that, the seasonal variations of K_s were principally explained by the changes in the crust biomass and possibly by the microbial exopolysaccharides. We conclude that the biocrusts significantly decreased soil infiltrability and exhibited a different seasonal variation pattern, which should be carefully considered in future analyses of hydropedological processes.     

Changes in Entrapped Gas Content and Hydraulic Conductivity with Pressure

Water table fluctuations continuously introduce entrapped air bubbles into the otherwise saturated capillary fringe and groundwater zone, which reduces the effective (quasi‐saturated) hydraulic conductivity, K_quasi, thus impacting groundwater flow, aquifer recharge and solute and contaminant transport. These entrapped gases will be susceptible to compression or expansion with changes in water pressure, as would be expected with water table (and barometric pressure) fluctuations. Here we undertake laboratory experiments using sand‐packed columns to quantify the effect of water table changes of up to 250 cm on the entrapped gas content and the quasi‐saturated hydraulic conductivity, and discuss our ability to account for these mechanisms in ground water models. Initial entrapped air contents ranged between 0.080 and 0.158, with a corresponding K_quasi ranging between 2 and 6 times lower compared to the K_s value. The application of 250 cm of water pressure caused an 18% to 26% reduction in the entrapped air content, resulting in an increase in K_quasi by 1.16 to 1.57 times compared to its initial (0 cm water pressure) value. The change in entrapped air content measured at pressure step intervals of 50 cm, was essentially linear, and could be modeled according to the ideal gas law. Meanwhile, the changes in K_quasi with compression–expansion of the bubbles because of pressure changes could be adequately captured with several current hydraulic conductivity models.