Estimating field‐saturated soil hydraulic conductivity by a simplified Beerkan infiltration experiment

Field‐saturated soil hydraulic conductivity, K_fs, is highly variable. Therefore, interpreting and simulating hydrological processes, such as rainfall excess generation, need a large number of K_fs data even at the plot scale. Simple and reasonably rapid experiments should be carried out in the field. In this investigation, a simple infiltration experiment with a ring inserted shortly into the soil and the estimation of the so‐called α* parameter allowed to obtain an approximate measurement of K_fs. The theoretical approach was tested with reference to 149 sampling points established on Burundian soils. The estimated K_fs with the value of first approximation of α* for most agricultural field soils (α* = 0.012 mm^?1) differed by a practically negligible maximum factor of two from the saturated conductivity obtained by the complete Beerkan Estimation of Soil Transfer parameters (BEST) procedure for soil hydraulic characterization. The measured infiltration curve contained the necessary information to obtain a site‐specific prediction of α*. The empirically derived α* relationship gave similar results for K_fs (mean = 0.085 mm s^?1; coefficient of variation (CV) = 71%) to those obtained with BEST (mean = 0.086 mm s^?1; CV = 67%), and it was also successfully tested with reference to a few Sicilian sampling points, since it yielded a mean and a CV of K_fs (0.0094 mm s^?1 and 102%, respectively) close to the values obtained with BEST (mean = 0.0092 mm s^?1; CV = 113%). The developed method appears attractive due to the extreme simplicity of the experiment. Copyright © 2012 John Wiley & Sons, Ltd.     

Soil water content and salinity determination using different dielectric methods in saline gypsiferous soil / Détermination de la teneur en eau et de la salinité de sols salins gypseux à l'aide de différentes méthodes diélectriques

Abstract

Measurements of dielectric permittivity and electrical conductivity were taken in a saline gypsiferous soil collected from southern Tunisia. Both time domain reflectometry (TDR) and the new WET sensor based on frequency domain reflectometry (FDR) were used. Seven different moistening solutions were used with electrical conductivities of 0.0053–14 dS m^?1. Different models for describing the observed relationships between dielectric permittivity (K _a ) and water content (θ), and bulk electrical conductivity (EC _a ) and pore water electrical conductivity (EC _p ) were tested and evaluated. The commonly used K _a –θ models by Topp et al. (1980) and Ledieu et al. (1986) cannot be recommended for the WET sensor. With these models, the RMSE and the mean relative error of the predicted θ were about 0.04 m³ m^?3 and 19% for TDR and 0.08 m³ m^?3 and 54% for WET sensor measurements, respectively. Using the Hilhorst (2000) model for EC _p predictions, the RMSE was 1.16 dS m^?1 and 4.15 dS m^?1 using TDR and the WET sensor, respectively. The WET sensor could give similar accuracy to TDR if calibrated values of the soil parameter were used instead of standard values.     

Persistence of road runoff generation in a logged catchment in Peninsular Malaysia

Measurements of saturated hydraulic conductivity (K_s) and diagnostic model simulations show that all types of logging road/trail in the 14·4 ha Bukit Tarek Experimental Catchment 3 (BTEC3) generate substantial Horton overland flow (HOF) during most storms, regardless of design and level of trafficking. Near‐surface K_s(0–0·05 m) on the main logging road, skid trails and newly constructed logging terraces was less than 1, 2 and 34 mm h^?1, respectively. Near‐surface K_s on an abandoned skid trail in an adjacent basin was higher (62 mm h^?1), owing to the development of a thin organic‐rich layer on the running surface over the past 40 years. Saturated hydraulic conductivity measured at 0·25 m below the surface of all roads was not different (all <6 mm h^?1) and corresponded to the K_s of the adjacent hillslope subsoil, as most roads were excavated into the regolith more than 0·5–1 m. After 40 years, only limited recovery in near‐surface K_s occurred on the abandoned skid trail. This road generated HOF after the storage capacity of the upper near‐surface layer was exceeded during events larger than about 20 mm. Thus, excavation into low‐K_s substrate had a greater influence on the persistence of surface runoff production than did surface compaction by machinery during construction and subsequent use during logging operations. Overland flow on BTEC3 roads was also augmented by the interception of shallow subsurface flow traveling along the soil–saprolite/bedrock interface and return flow emerging from the cutbank through shallow biogenic pipes. The most feasible strategy for reducing long‐term road‐related impacts in BTEC3 is limiting the depth of excavation and designing a more efficient road network, including minimizing the length and connectivity of roads and skid trails. Copyright © 2007 John Wiley & Sons, Ltd.     

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.     

Infiltration and runoff generation processes in fire‐affected soils

Post‐wildfire runoff was investigated by combining field measurements and modelling of infiltration into fire‐affected soils to predict time‐to‐start of runoff and peak runoff rate at the plot scale (1 m²). Time series of soil‐water content, rainfall and runoff were measured on a hillslope burned by the 2010 Fourmile Canyon Fire west of Boulder, Colorado during cyclonic and convective rainstorms in the spring and summer of 2011. Some of the field measurements and measured soil physical properties were used to calibrate a one‐dimensional post‐wildfire numerical model, which was then used as a ‘virtual instrument’ to provide estimates of the saturated hydraulic conductivity and high‐resolution (1 mm) estimates of the soil‐water profile and water fluxes within the unsaturated zone. Field and model estimates of the wetting‐front depth indicated that post‐wildfire infiltration was on average confined to shallow depths less than 30 mm. Model estimates of the effective saturated hydraulic conductivity, K_s, near the soil surface ranged from 0.1 to 5.2 mm h^?1. Because of the relatively small values of K_s, the time‐to‐start of runoff (measured from the start of rainfall), t_p, was found to depend only on the initial soil‐water saturation deficit (predicted by the model) and a measured characteristic of the rainfall profile (referred to as the average rainfall acceleration, equal to the initial rate of change in rainfall intensity). An analytical model was developed from the combined results and explained 92–97% of the variance of t_p, and the numerical infiltration model explained 74–91% of the variance of the peak runoff rates. These results are from one burned site, but they strongly suggest that t_p in fire‐affected soils (which often have low values of K_s) is probably controlled more by the storm profile and the initial soil‐water saturation deficit than by soil hydraulic properties. Published 2013. This article is a U.S. Government work and is in the public domain in the USA.     

滆湖水体光学性质初步研究

基于2009年7月至2010年6月滆湖全湖15个采样点的水体光学参数及相关水质理化因子数据,分析滆湖水体周年光合有效辐射(PAR)衰减特性,以期为滆湖沉水植物生态修复提供相关水体光学资料.结果表明,滆湖水体PAR衰减系数(Kd)周年变化范围为1.32～17.42 m-1.秋季Kd相对最小,平均值为2.35 m-1,变化范围为1.32 ～3.70 m-1;夏季Kd相对最大,平均值为6.23 m-1,变化范围为3.68～17.42 m-1.春、秋、冬季,滆湖水体真光层平均深度均满足沉水植物的生长需求,而在夏季滆湖水体真光层平均深度仅为0.84m,小于全湖平均水深(1.20 m),因此夏季PAR是限制沉水植物恢复的因子之一.滆湖水体Kd与透明度(SD)在秋、冬季的关系为:Kd =2.089 +0.705/SD.叶绿素a浓度和悬浮物浓度是影响滆湖水体Kd的重要因子之一.     

Multi‐decadal water table manipulation alters peatland hydraulic structure and moisture retention

A peatland complex disturbed by berm construction in the 1950s was used to examine the long‐term impact of water table (WT) manipulation on peatland hydraulic properties and moisture retention at three adjacent sites with increasing depth to WT (WET, INTermediate reference and DRY). Saturated hydraulic conductivity (K_s) was found to decrease with depth by several orders of magnitude over a depth of 1–1.5 m at all sites. The depth dependence of WT response to rainfall was similar across sites: WT response increased from 1 : 1 at the surface, to 5 : 1 at 50 cm depth. While surface specific yield (S_y) values were similar across all sites, it decreased with depth at a rate of 0.014 cm^?1 in hollows and 0.007 cm^?1 in hummocks. Bulk density (ρ_b) exhibited similar depth‐dependent trends as S_y and explains a high amount of variance (r² > 0.69) in moisture retention across a range of pore water pressures (?15 to ?500 cm H₂O). Because of higher ρ_b, hollow peat had greater moisture retention, where site effects were minimal. However, the estimated residual water content for surface Sphagnum samples, while on average lower in hummocks (0.082 m³ m^?3) versus hollows (0.087 m³ m^?3), increased from WET (0.058 m³ m^?3) to INT (0.088 m³ m^?3) to DRY (0.108 m³ m^?3) which has important implications for moisture stress under conditions of persistent WT drawdown. Given the potential importance of microtopographic succession for altering peatland hydraulic structure, our findings point to the need for a better understanding of what controls the relative height and proportional coverage of hummocks in relation to long‐term disturbance‐response dynamics. Copyright © 2014 John Wiley & Sons, Ltd.     

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.     

Predicting physical equations of soil detachment by simulated concentrated flow in Ultisols (subtropical China)

Soil detachment in concentrated flow is due to the dislodging of soil particles from the soil matrix by surface runoff. Both aggregate stability and shear strength of the topsoil reflect the erosion resistance of soil to concentrated runoff, and are important input parameters in predicting soil detachment models. This study was conducted to develop a formula to predict soil detachment rate in concentrated flow by using the aggregate stability index (As), root density (R_d) and saturated soil strength (σ_s) in the subtropical Ultisols region of China. The detachment rates of undisturbed topsoil samples collected from eight cultivated soil plots were measured in a 3.8 m long, 0.2 m wide hydraulic flume under five different flow shear stresses (τ = 4.54, 9.38, 15.01, 17.49 and 22.54 Pa). The results indicated that the stability index (As) was well related with soil detachment rate, particularly for results obtained with high flow shear stress (22.54 Pa), and the stability index (As) has a good linear relationship with concentrated flow erodibility factors (K_c). There was a positive linear relationship between saturated soil strength (σ_s) and critical flow shear stress (τ_c) for different soils. A significant negative exponential relationship between erodibility factors (K_c) and root density (R_d) was detected. This study yielded two prediction equations that allowed comparison of their efficiency in assessing soil detachment rate in concentrated flow. The equation including the root density (R_d) may have a better correlation coefficient (R² = 0.95). It was concluded that the formula based on the stability index (As), saturated soil strength (σ_s) and root density (R_d) has the potential to improve methodology for assessing soil detachment rate in concentrated flow for the subtropical Chinese Ultisols. Copyright © 2012 John Wiley & Sons, Ltd.     

The thermal boundary-layer interpretation of D″ and its role as a plume source

The “anomalous” layer in the lowermost mantle, identified as D″ in the notation of K.E. Bullen, appears in the PREM Earth model as a 150 km-thick zone in which the gradient of incompressibility with pressure, $\text{d}\text{K}\text{d}\text{P}$, is almost 1.6, instead of 3.2 as in the overlying mantle. Since PREM shows no accompanying change in density or density gradient, we identify D″ as a thermal boundary layer and not as a chemically distinct zone. The anomaly in $\text{d}\text{K}\text{d}\text{P}$ is related to the temperature gradient by the temperature dependence of K_s, for which we present a thermodynamic identity in terms of accessible quantities. This gives the numerical result (?K_s/?T)_P=?1.6×10⁷ Pa K^?1 for D″ material. The corresponding temperature increment over the D″ range is 840 K. Such a layer cannot be a static feature, but must be maintained by a downward motion of the lower mantle toward the core-mantle boundary with a strong horizontal flow near the base of D″. Assuming a core heat flux of 1.6 × 10¹² W, the downward speed is 0.07 mm y^?1 and the temperature profile in D″, scaled to match PREM data, is approximately exponential with a scale height of 73 km. The inferred thermal conductivity is 1.2 W m^?1 K^?1. Using these values we develop a new analytical model of D″ which is dynamically and thermally consistent. In this model, the lower-mantle material is heated and softened as it moves down into D″ where the strong temperature dependence of viscosity concentrates the horizontal flow in a layer ～ 12 km thick and similarly ensures its removal via narrow plumes.     

Effects of biological soil crusts on soil detachment process by overland flow in the Loess Plateau of China

Biological soil crusts (BSCs) cover up to 60 to 70% of the soil surface in grasslands after the ‘Grain for Green’ project was implemented in 1999 to rehabilitate the Loess Plateau. However, few studies exist that quantify the effects of BSCs on the soil detachment process by overland flow in the Loess Plateau. This study investigated the potential effects of BSCs on the soil detachment capacity (D_c), and soil resistance to flowing water erosion reflected by rill erodibility and critical shear stress. Two dominant BSC types that developed in the Loess Plateau (the later successional moss and the early successional cyanobacteria mixed with moss) were tested against natural soil samples collected from two abandoned farmland areas. The samples were subjected to flow scouring under six different shear stresses ranging from 7.15 to 24.08 Pa. The results showed that D_c decreased significantly with crust coverage under both moss and mixed crusts. The mean D_c of bare soil (0.823 kg m^?2 s^?1) was 2.9 to 48.4 times greater than those of moss covered soil (0.017–0.284 kg m^?2 s^?1), while it (3.142 kg m^?2 s^?1) was 4.9 to 149.6 times greater than those of mixed covered soil (0.021–0.641 kg m^?2 s^?1). The relative detachment rate of BSCs compared with bare soils decreased exponentially with increasing BSC coverage for both types of BSCs. The D_c value can be simulated by flow shear stress, cohesion, and BSC coverage using a power function (NSE ≥ 0.59). Rill erodibility also decreased with coverage of both crust types. Rill erodibility of bare soil was 3 to 74 times greater than those of moss covered soil and was 2 to 165 times greater than those of mixed covered soil. Rill erodibility could also be estimated by BSC coverage in the Loess Plateau (NSE ≥ 0.91). The effect of crust coverage on critical shear stress was not significant. Copyright © 2015 John Wiley & Sons, Ltd.     

Measuring streambank erosion due to ground water seepage: correlation to bank pore water pressure,precipitation and stream stage

Limited information exists on one of the mechanisms governing sediment input to streams: streambank erosion by ground water seepage. The objective of this research was to demonstrate the importance of streambank composition and stratigraphy in controlling seepage flow and to quantify correlation of seepage flow/erosion with precipitation, stream stage and soil pore water pressure. The streambank site was located in Northern Mississippi in the Goodwin Creek watershed. Soil samples from layers on the streambank face suggested less than an order of magnitude difference in vertical hydraulic conductivity (K_s) with depth, but differences between lateral K_s of a concretion layer and the vertical K_s of the underlying layers contributed to the propensity for lateral flow. Goodwin Creek seeps were not similar to other seeps reported in the literature, in that eroded sediment originated from layers underneath the primary seepage layer. Subsurface flow and sediment load, quantified using 50 cm wide collection pans, were dependent on the type of seep: intermittent low‐flow (LF) seeps (flow rates typically less than 0·05 L min^?1), persistent high‐flow (HF) seeps (average flow rate of 0·39 L min^?1) and buried seeps, which eroded unconsolidated bank material from previous bank failures. The timing of LF seeps correlated to river stage and precipitation. The HF seeps at Goodwin Creek began after rainfall events resulted in the adjacent streambank reaching near saturation (i.e. soil pore water pressures greater than ?5 kPa). Seep discharge from HF seeps reached a maximum of 1·0 L min^?1 and sediment concentrations commonly approached 100 g L^?1. Buried seeps were intermittent but exhibited the most significant erosion rates (738 g min^?1) and sediment concentrations (989 g L^?1). In cases where perched water table conditions exist and persistent HF seeps occur, seepage erosion and bank collapse of streambank sediment may be significant. Copyright © 2007 John Wiley & Sons, Ltd.     

Spectral attenuation of solar radiation in Patagonian fjord and coastal waters and implications for algal photobiology

The spectral attenuation of solar irradiation was measured during summer in two types of coastal waters in southern Chile, a north Patagonian fjord (Seno Reloncaví) and open coast (Valdivia). In order to relate the light availability with the light requirements of upper subtidal seaweeds, the saturating irradiance for photosynthesis (E_k) from P–I curves was measured. In addition the UV risk was assessed. Based on the z_1% of PAR, the lower limit of the euphotic zone in the studied systems averaged 21 m (K_d 0.24 m^?1) in Seno Reloncaví and 18 m (K_d 0.27 m^?1) in the coast of Valdivia. Photosynthesis of the studied seaweeds was saturated at markedly lower irradiances than found in their natural depths at the time of the study. Solar radiation penetrating into these depths at both locations largely supports the light requirements for the photosynthesis of subtidal species: 50–160 μmol m^?2 s^?1 for seaweeds from Seno Reloncaví (7 m tidal range) and 20–115 μmol m^?2 s^?1 for Valdivia assemblages (2 m tidal range). Optimal light conditions to saturate photosynthesis (E_k) were present at 10–16 m water depth. The attenuation of solar irradiation did not vary significantly between the fjord and coastal sites of this study. However, the underwater light climates to which seaweeds are exposed in these sites vary significantly because of the stronger influence of tidal range affecting the fjord system as compared with the open coastal site. The patterns of UV-B penetration in these water bodies suggest that seaweeds living in upper littoral zones such as the intertidal and shallow subtidal (<3 m) may be at risk.     

Occurrence of water ponding on soil surfaces depending on infiltration rates on Mongolian rangeland

The occurrence of water ponding on soil surfaces during and after heavy rainfall produces surface run‐off or surface water accumulation in low‐lying areas, which might reduce the water supply to soils and result in a reduction of the soil water that plants can use, especially in arid climates. On Mongolian rangeland, we observed ponded water on the surface of a specific soil condition subjected to a heavy rainfall of 30 mm/hr. By contrast, ponded water was not observed for the same type of soil where livestock grazing had been removed for 6–8 years via a fence or for nearby soil containing less clay. We measured the infiltration rate (the saturated hydraulic conductivity of the surface soil, K_s) of the three sites by applying ponded water on the soil surface (an intake rate test). The results showed that K_s in the rangeland was lower than the rainfall intensity in the site where water ponded on the soil surface; however, K_s of the soil inside of the fence has recovered to 3 times that of the soil outside of the fence to exceed the rainfall intensity. Heavy rainfall that exceeds the infiltration rate occurs several times a year at the livestock grazing site where we observed ponded water. Slight water repellency of the soil reduces rain infiltration to increase the possibility of surface ponding for the soil.     

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.     

A dripper-TDR method for in situ determination of hydraulic conductivity and chemical transport properties of surface soils

Field determined hydraulic and chemical transport properties can be useful for the protection of groundwater resources from land-applied chemicals. Most field methods to determine flow and transport parameters are either time or energy consuming and/or they provide a single measurement for a given time period. In this study, we present a dripper-TDR field method that allows measurement of hydraulic conductivity and chemical transport parameters at multiple field locations within a short time period. Specifically, the dripper-TDR determines saturated hydraulic conductivity (K_s), macroscopic capillary length (λ_c), immobile water fraction (θ_im/θ), mass exchange coefficient (α) and dispersion coefficient (D_m). Multiple dripper lines were positioned over five crop rows in a field. Background and step solutions were applied through drippers to determine surface hydraulic conductivity parameters at 44 locations and surface transport properties at 38 locations. The hydraulic conductivity parameters (K_s, λ_c) were determined by application of three discharge rates from the drippers and measurements of the resultant steady-state flux densities at the soil surface beneath each dripper. Time domain reflectometry (TDR) was used to measure the bulk electrical conductivity of the soil during steady infiltration of a salt solution. Breakthrough curves (BTCs) for all sites were determined from the TDR measurements. The K_s and λ_c values were found to be lognormally distributed with average values of 31.4 cm h⁻¹ and 6.0 cm, respectively. BTC analysis produced chemical properties, θ_im/θ, α, and D_m with average values of 0.23, 0.0036 h⁻¹, and 1220 cm² h⁻¹, respectively. The estimated values of the flow and transport parameters were found to be within the ranges of values reported by previous studies conducted at nearby field locations. The dripper TDR method is a rapid and useful technique for in situ measurements of hydraulic conductivity and solute transport properties. The measurements reported in this study give clear evidence to the occurrence of non-equilibrium water and chemical movement in surface soil. The method allows for quantification of non-equilibrium model parameters and preferential flow. Quantifying the parameters is a necessary step toward determining the influences of surface properties on infiltration, runoff, and vadose zone transport.     

Pan evaporation modelling and changing attribution analysis on the Tibetan Plateau (1970–2012)

Pan evaporation (E_p) is an important indicator of water and energy and the decline of E_p has been reported in many regions over the last decades. The climate and E_p are dependent on each other. In this study, the temporal trends of E_p and main E_p drivers, namely mean air temperature (T_a), wind speed (u), global solar radiation (R_s), net long‐wave radiation(R_nl) and vapour pressure deficit (D) from 1970 to 2012, were calculated on the basis of 26 meteorological stations on the Tibetan Plateau. The arithmetic average of E_p from 26 stations decreased with the rate of ?11.91 mm a^?2; the trends of R_s, R_nl, T_a, u and D were ?1.434 w m^?2 decade^?1, 0.2511 w m^?2 decade^?1, 0.3590°C decade^?1, ?0.2376 m s^?1 decade^?1 and 9.523 Pa decade^?1, respectively. The diffuse irradiance is an essential parameter to model E_p and quantify the contribution of climatic factors to changing E_p. 60 724 observations of R_s and diffuse solar irradiance (R_d) from seven of the 26 stations were used to develop the correlation between the diffuse fraction (R_d/R_s), and the clearness index (R_s/R_o). On the basis of the estimation of the diffuse component of R_s and climatic data, we modified the PenPan model to estimate Chinese micro‐pan evaporation (E_p) and assess the attribution of E_p dynamics using partial derivatives. The results showed that there was a good agreement between the observed and calculated daily E_p values. The observed decrease in E_p was mostly due to declining wind speed (?13.7 mm a^?2) with some contributions from decreasing solar irradiance (?3.1 mm a^?2); and the increase of temperature had a large positive effect (4.55 mm a^?2) in total whilst the increase of R_nl had insignificant effect (0.35 mm a^?2) on E_p rates. The change of E_p is the net result of all the climatic variables. Copyright © 2014 John Wiley & Sons, Ltd.     

The effect of interstitial gaseous pressure on the thermal conductivity of a simulated Apollo 12 lunar soil sample

The thermal conductivity of a simulated Apollo 12 lunar soil sample was measured with a needle probe under vacuum. The result showed that the sample, with bulk densities of 1.70–1.85 g cm^?3 held in a vertical cylinder (2.54 cm in diameter and 6.99 cm long) has a thermal conductivity ranging from 8.8 to 10.9 mW m^?1 K^?1. This is comparable to the lunar regolith's thermal conductivity as determined in situ. Besides the dense packing of the soil particles, an enhanced intergranular thermal contact, due to the self-compression of the sample, is necessary to raise the sample's thermal conductivity from the level of loose soil (< 5 mW m^?1 K^?1) to that of the lunar regolith deeper than 35 cm (～ 10 mW m^?1 K^?1). A model of the lunar regolith, a thin layer of loose soil resting on a compacted self-compressed substratum, is consistent with the lunar regolith's surface structure as deduced from an observation of the lunar surface's brightness temperature. Martian regolith surface structure is similar, except that its surface layer may be missing in places because of aeolian activity. Measurements of thermal conductivity under simulated martian surface conditions showed that the thermal properties of loose and compacted soils agreed with the two peak values of the martian surface's thermal inertia as observed from “Viking” orbiters, suggesting that drifted loose soil and exposed compacted soil are responsible for the bimodal distribution of the martian surface's thermal inertia near zero elevation. For compacted soil exposed to the martian surface to have the same thermal conductivity as that buried under the surface layer, a cohesion of the soil particles must be assumed.     

青藏高原典型湖库光学吸收特性与光合有效辐射衰减系数初步研究

青藏高原地区分布的湖泊数量众多、面积较大、分布范围广泛.受制于恶劣的自然条件,对该地区湖泊的光学吸收特性以及光合有效辐射衰减系数(Kd(PAR))的研究鲜有成果.本文依据2014和2015年间采集的13个典型高海拔、湖泊面积较大的湖库的现场实测数据和实验室测定数据,分析了采样湖库各个采样点的Kd(PAR)特征以及有色可溶性有机物(CDOM)、藻类颗粒物吸收及非色素颗粒物吸收特性,计算并分析了Kd(PAR)与透明度以及光学活性物质的关系.研究结果表明:青藏高原地区各湖库平均各项颗粒物吸收系数均较低,总颗粒物吸收系数在400~700 nm波段内不超过0.14 m~(-1)、CDOM吸收系数在355 nm波长处最高,为1.23 m~(-1)、最低接近于0、藻类颗粒物吸收特性不明显;实验数据完整的巴木错、格仁错和班公错的主导吸收组分各异,其中巴木错为CDOM吸收主导,格仁错与班公错为非色素颗粒物吸收主导;青藏高原采样湖泊总体Kd(PAR)平均值较小,仅为0.26 m~(-1),样点最大值出现在可鲁克湖(1.17 m~(-1)),最小值出现在普莫雍错(0.10 m~(-1));在采样湖泊中Kd(PAR)与透明度呈显著相关;Kd(PAR)与CDOM的相关性最强,叶绿素a浓度次之,与总悬浮颗粒物浓度的相关性最不显著.     

Factors controlling the chemical evolution of travertine‐depositing rivers of the Barkly karst,northern Australia

Groundwaters feeding travertine‐depositing rivers of the northeastern segment of the Barkly karst (NW Queensland, Australia) are of comparable chemical composition, allowing a detailed investigation of how the rate of downstream chemical evolution varies from river to river. The discharge, pH, temperature, conductivity and major‐ion concentrations of five rivers were determined by standard field and laboratory techniques. The results show that each river experiences similar patterns of downstream chemical evolution, with CO₂ outgassing driving the waters to high levels of calcite supersaturation, which in turn leads to widespread calcium carbonate deposition. However, the rate at which the waters evolve, measured as the loss of CaCO₃ per kilometre, varies from river to river, and depends primarily upon discharge at the time of sampling and stream gradient. For example, Louie Creek (Q = 0·11 m³ s^?1) and Carl Creek (Q = 0·50 m³ s^?1) have identical stream gradients, but the loss of CaCO₃ per kilometre for Louie Creek is twice that of Carl Creek. The Gregory River (Q = 3·07 m³ s^?1), O'Shanassy River (Q = 0·57 m³ s^?1) and Lawn Hill Creek (Q = 0·72 m³ s^?1) have very similar gradients, but the rate of hydrochemical evolution of the Gregory River is significantly less than either of the other two systems. The results have major implications for travertine deposition: the stream reach required for waters to evolve to critical levels of calcite supersaturation will, all others things being equal, increase with increasing discharge, and the length of reach over which travertine is deposited will also increase with increasing discharge. This implies that fossil travertine deposits preserved well downstream of modern deposition limits are likely to have been formed under higher discharge regimes. Copyright © 2002 John Wiley & Sons, Ltd.