Simulations of water movement and solute transport through different soil texture configurations under negative‐pressure irrigation

This study examined the effects of different soil texture configurations on water movement and solute transport to provide a reliable scientific basis for the application of negative‐pressure irrigation (NPI) technology. HYDRUS‐2D was used to analyse water movement and solute transport under NPI. The main results are as follows: (a) HYDRUS‐2D can be used to simulate water movement and solute transport under NPI, as there was good agreement between the simulated and measured values for water contents, NaCl concentrations, cumulative water infiltration, and wetting distances in the horizontal and vertical directions; the Nash–Sutcliffe efficiency coefficients were in the range of 0.94–0.97. (b) Layered soils have obvious effects on water movement under NPI. With the emitter position in the loam layer, when a coarse texture of loamy sand was present below the loam layer (namely, L‐LS), irrigation water accumulated in the topsoil, and this led to an increase in evaporation compared with the homogeneous loam profile. However, fine texture silty loam or silty clay loam layers beneath the loam layer (namely, L‐SiL or L‐SiCL, respectively) was more conducive to water infiltration into the lower layer, and this increased the amount of water infiltration and simultaneously reduced the surface evaporation effectively. (c) Layered soils have obvious effects on solute transport under NPI, and salt accumulation will readily occur in the clay‐rich soil layer at the interface. The maximum soil salt accumulation of L‐LS occurred above the soil interface between the two soil layers with a value of 21.80 g/kg; however, for L‐SiCL and L‐SiL, the maximum salt accumulation occurred below the soil interface between the two soil layers, with values of 23.80 g/kg and 20.08 g/kg, respectively. (d) Interlayered soils showed remarkable changes in the water infiltration characteristics and salt‐leaching intensities under NPI, and the properties for the soil profile with a silty loam interlayer were better than those for the soil profile with a silty clay loam interlayer. The soil profile with a loamy sand interlayer had the lowest amount of water infiltration, which resulted in reductions of the salt‐leaching intensities. Thus, NPI is clearly not suitable for loamy sand soil. Overall, the results demonstrated that soil texture configurations affected water movement and solute transport under NPI. Therefore, careful consideration should be given to the use of NPI to achieve target soil water and solution conditions and reduce water loss.     

Modelling water infiltration in cracked paddy field soil

Periodic paddy field flooding is a major source of groundwater recharge. Many paddy fields thus are used as groundwater recharge ponds after harvesting the first crop of the summer. Following rice harvesting, paddy field surfaces may crack into fissures as a result of drainage and exposure to sunlight. Field observation indicates that applying precipitation to the paddy field can increase the rate of infiltration. To quantitatively evaluate the amount of infiltration in a cracked paddy field, this study sets up a simple soil crack model to simulate the field infiltration process. A three‐dimensional groundwater model FEMWATER is adopted to simulate water movement in the paddy field subjected to various crack conditions. Using the field and laboratory data of irrigation water requirements, soil physical properties, hydraulic conductivities and soil profiles obtained from Ten‐Chung, FEMWATER simulates the water movement in the dry cracked paddy. Simulation results show that if the cracks develop extensively and penetrate the ploughed soil, the infiltration rate may increase significantly. The infiltration fluxes of crack with depths of 80, 60 and 27·5 cm are 18·77, 14·50 and 8·06 times higher than that of 20 cm, respectively. The simulation results of cracks with 80 cm depth correlated closely with field observations. The results of the study elucidate the processes of unsaturated water movement in a dry cracked paddy field. Copyright © 2004 John Wiley & Sons, Ltd.     

Characterising deep vadose zone water movement and solute transport under typical irrigated cropland in the North China Plain

Understanding the dynamics and mechanisms of soil water movement and solute transport is essential for accurately estimating recharge rates and evaluating the impacts of agricultural activities on groundwater resources. In a thick vadose zone (0–15 m) under irrigated cropland in the piedmont region of the North China Plain, soil water content, matric potential, and solute concentrations were measured. Based on these data, the dynamics of soil water and solutes were analysed to investigate the mechanisms of soil water and solute transport. The study showed that the 0–15‐m vadose zone can be divided into three layers: an infiltration and evaporation layer (0–2 m), an unsteady infiltration layer (2–6 m), and a quasi‐steady infiltration layer (6–15 m). The chloride, nitrate, and sulphate concentrations all showed greater variations in the upper soil layer (0–1 m) compared to values in the deep vadose zone (below 2 m). The average concentrations of these three anions in the deep vadose zone varied insignificantly with depth and approached values of 125, 242, and 116 mg/L. The accumulated chloride, sulphate, and nitrate were 2,179 ± 113, 1,760 ± 383, and 4,074 ± 421 kg/ha, respectively. The soil water potential and solute concentrations indicated that uniform flow and preferential flow both occurred in the deep vadose zone, and uniform flow was the dominant mechanism of soil water movement in this study. The piston‐like flow velocity of solute transport was 1.14 m per year, and the average value of calculated leached nitrate nitrogen was 107 kg/ha?year below the root zone. The results can be used to better understand recharge processes and improve groundwater resources management.     

LISEM: A SINGLE-EVENT PHYSICALLY BASED HYDROLOGICAL AND SOIL EROSION MODEL FOR DRAINAGE BASINS. I: THEORY,INPUT AND OUTPUT

A new physically based hydrological and soil erosion model has been developed, which can be used for planning and conservation purposes: the LImburg Soil Erosion Model (LISEM). The LISEM model is one of the first examples of a physically based model that is completely incorporated in a raster Geographical Information System. This incorporation facilitates easy application in larger catchments, improves the user friendliness by avoiding conversion routines and allows remotely sensed data to be used. Processes incorporated in the model are rainfall, interception, surface storage in micro-depressions, infiltration and vertical movement of water in the soil, overland flow, channel flow, detachment by rainfall and throughfall, detachment by overland flow and transport capacity of the flow. Special attention has been given to the influence of tractor wheelings, small roads and surface sealing. Vertical movement of water in the soil is simulated using the Richard's equation. Optionally, the user can choose the Holtan or the Green–Ampt infiltration model. For the distribution flow routing, a four-point finite-difference solution of the kinematic wave is used together with Manning's equation.     

Modified models for better prediction of infiltration rates in trapezoidal permeable stormwater channels

ABSTRACT

In stormwater management, it is important to accurately quantify the infiltration rates to solve urban runoff-related problems. This study proposes a method to improve estimates of the infiltration rate in permeable stormwater channels. As part of the analysis, five infiltration models were evaluated: the Kostiakov, Horton, modified Kostiakov, Philip and SCS (Soil Conservation Service) models. Infiltration tests with various initial water levels were performed on channel models with differing base width and side slopes. The results show that the addition of three parameters that describe the trapezoidal cross-sectional area, i.e. the depth, side slope and base width, in the infiltration models yielded better estimates of the infiltration rate. A comparison of the infiltration capacity values obtained from the models after the three parameters were added with those that were experimentally obtained, shows that the improved modified Kostiakov model is the most suitable model to predict infiltration rates in trapezoidal permeable stormwater channels.     

Surfactant and irrigation effects on wettable soils: runoff,erosion, and water retention responses

Surfactants are chemical compounds that can change the contact angle of a water drop on solid surfaces and are commonly used to increase infiltration into water repellent soil. Since production fields with water repellent soil often contain areas of wettable soil, surfactants applied to such fields worldwide will likely be applied to wettable soil, with unknown consequences for irrigation‐induced erosion, runoff, or soil water relations. We evaluated surfactant and simulated sprinkler irrigation effects on these responses for three wettable, Pacific Northwest soils, Latahco and Rad silt loams, and Quincy sand. Along with an untreated control, we studied three surfactants: an alkyl polyglycoside (APG) in solution at a concentration of 18 g active ingredient (AI) kg^?1, a block copolymer at 26 g kg^?1, and a blend of the two at 43 g kg^?1. From 2005 to 2009 in the laboratory, each surfactant was sprayed at a rate of 46·8 l ha^?1 onto each soil packed by tamping into 1·2‐ by 1·5‐m steel boxes. Thereafter, each treated soil was irrigated twice at 88 mm h^?1 with surfactant‐free well water. After each irrigation, runoff and sediment loss were measured and soil samples were collected. While measured properties differed among soils and irrigations, surfactants had no effect on runoff, sediment loss, splash loss, or tension infiltration, compared to the control. Across all soils, however, the APG increased volumetric water contents by about 3% (significant at p≤0·08) at matric potentials from 0 to ? 20 kPa compared to the control. With a decrease in the liquid–solid contact angle on treated soil surfaces, surfactant‐free water appeared able to enter, and be retained in pores with diameters ≥ 15 µm. All told, surfactants applied at economic rates to these wettable Pacific Northwest soils posed little risk of increasing either runoff or erosion or harming soil water relations. Moreover, by increasing water retention at high potentials, surfactants applied to wettable soils may allow water containing pesticides or other agricultural chemicals to better penetrate soil pores, thereby increasing the efficacy of the co‐applied materials. Copyright © 2010 John Wiley & Sons, Ltd.     

Soil macroporosity and infiltration characteristics of a forest podzol

Spatial distribution of soil macroporosity was determined for a forest podzol from tension infiltrometer measurements at the soil surface. Surface‐derived macroporosity values were compared with point infiltration characteristics obtained from soil water content and soil water chemistry measurements during an experimental irrigation, and with parameters of a kinematic wave model applied to soil water content data. Macroporosity estimated by the tension infiltrometer ranged from 0·00087 to 0·0219% of soil volume, and infiltration at these two sites was dominated by propagation of a well‐defined wetting front through the soil profile and bypass flow via soil macropores, respectively. Infiltration at sites with intermediate macroporosities reflected a combination of these two processes, although results were inconclusive at one site owing to lateral flow at the base of the soil profile. There was no agreement between macroporosities estimated by the tension infiltrometer and the kinematic wave model. The maximum soil conductance parameter within the profile at a site, however, was related directly to the surface‐derived macroporosity. The partial agreement between surface‐derived macroporosity estimates and point infiltration characteristics shown here supports the use of tension infiltrometry as a rapid, non‐destructive method of assessing spatial variations in the relative contribution of macropore flow to the infiltration process. Copyright © 2000 John Wiley & Sons, Ltd.     

Modelling soil solute release and transport in run‐off on a loessial slope with and without surface stones

Stone covers on loessial slopes can increase the time of infiltration by slowing the velocity of the overland flow, which reduces the transport of solutes, but few mechanistic models have been tested under water‐scouring conditions. We carried out field experiments to test a previously proposed, physically based model of water and solute transport. The area of soil infiltration was calculated from the uncovered surface area, and Richards' equation and the kinematic wave equation were used to describe water infiltration and flow along slopes with stone covers. The transport of chemicals into the run‐off from the surface soil, presumably by diffusion, and their movement in the soil profile could be described by the convection–diffusion equations of the model. The simulated and measured data correlated well. The stones on the soil surface reduced the area available for infiltration but increased the Manning coefficient, eventually leading to increased water infiltration and decreased solute loss with run‐off. Our results indicated that the traditional model of water movement and solute migration could be used to simulate water transport and solute migration for stone‐covered soil on loessial slopes.     

A soil‐water‐balance approach to quantify groundwater recharge from irrigated cropland in the North China Plain

Rapidly depleting unconfined aquifers are the primary source of water for irrigation on the North China Plain. Yet, despite its critical importance, groundwater recharge to the Plain remains an enigma. We introduce a one‐dimensional soil‐water‐balance model to estimate precipitation‐ and irrigation‐generated areal recharge from commonly available crop and soil characteristics and climate data. To limit input data needs and to simplify calculations, the model assumes that water flows vertically downward under a unit gradient; infiltration and evapotranspiration are separate, sequential processes; evapotranspiration is allocated to evaporation and transpiration as a function of leaf‐area index and is limited by soil‐moisture content; and evaporation and transpiration are distributed through the soil profile as exponential functions of soil and root depth, respectively. For calibration, model‐calculated water contents of 11 soil‐depth intervals from 0 to 200 cm were compared with measured water contents of loam soil at four sites in Luancheng County, Hebei Province, over 3 years (1998–2001). Each 50‐m² site was identically cropped with winter wheat and summer maize, but received a different irrigation treatment. Average root mean‐squared error between measured and model‐calculated water content of the top 180 cm was 4·2 cm, or 9·3% of average total water content. In addition, model‐calculated evapotranspiration compared well with that measured by a large‐scale lysimeter. To test the model, 12 additional sites were simulated successfully. Model results demonstrate that drainage from the soil profile is not a constant fraction of precipitation and irrigation inputs, but rather the fraction increases as the inputs increase. Because this drainage recharges the underlying aquifer, improving irrigation efficiency by reducing seepage will not reverse water‐table declines. Copyright © 2003 John Wiley & Sons, Ltd.     

Application of a water infiltration model for simulating water repellency of humus soil

One-dimensional infiltration experiments were conducted using hydrophilic and water-repellent soils from the Guishui River Basin to study the effects of soil water repellency on cumulative infiltration (CI) and the infiltration rate (IR). The test results show that, for the hydrophilic soil (HS) sample, CI increases monotonously with time and IR decreases monotonously. For the water-repellent soil (W-RS), however, the following characteristics were observed: (a) There is an inflection point in CI and a sudden increase in IR. Larger values of the initial soil water content produce an earlier and more significant inflection point in CI, and a larger peak value of IR. (b) The post-peak stable IR is greater than the pre-peak value, ignoring the beginning of rapid infiltration, and the overall IR presents a single peak. The applicability of various water infiltration models was analysed for the two soil types. Numerical analysis suggests the following conclusions: (a) For both HS and W-RS, the Kostiakov function, Gamma function, and Beta function (BF) models exhibit good applicability. (b) For W-RS, the Gauss function model not only reflects the monotonous decrease in IR, but also produces a steady IR in the initial stage, a gradual increase before the peak value, and a gradual decrease after the peak value. Similarly, the BF model reflects the monotonous decrease in IR. A piecewise BF reproduces the U-shaped change in rapid infiltration before the inflection point, as well as the gradual increase and right-skewed distribution curve of W-RS infiltration before and after the inflection point. The BF model achieves the best simulation accuracy and has the widest applicability.     

Using bromide data tracer and HYDRUS-1D to estimate groundwater recharge and evapotranspiration under film-mulched drip irrigation in an arid Inland Basin,Northwest China

Accurate estimation of groundwater recharge (GR) and evapotranspiration (ET) are essential for sustainable management of groundwater resources, especially in arid and semi-arid regions. In the Manas River Basin (MRB), water shortage is the main factor restricting sustainable development of irrigated agriculture, which relies heavily on groundwater. Film-mulched drip irrigation significantly changes the pattern and dominant processes of water flow in the unsaturated zone, which increases the difficulty of GR and ET estimation. To better estimate GR and ET under film-mulched drip irrigation in the MRB, bromide tracer tests and soil lithologic investigation were conducted at 12 representative sites. A one-dimensional variably saturated flow model (HYDRUS-1D) was calibrated at each site using soil evaporation data inferred from the bromide tracer tests. The results showed that average annual soil evaporation in uncultivated lands calculated from bromide trace tests was 25.55 mm. The annual GR ranged from 5.5 to 37.0 mm under film-mulched drip irrigation. The annual ET ranged from 507.0 to 747.1 mm, with soil evaporation between 35.7 and 117.0 mm and transpiration between 460.9 and 642.3 mm. Soil evaporation represented 7% to 16% of the total ET and more than 70% of precipitation and irrigation water was used by cotton plants. Spatial variations of soil lithology, water table depth and initial soil water content led to the spatial differences of GR and ET in the MRB. Our study indicated that bromide tracer tests are useful for inferring ET in the arid and semi-arid oases. The combination of bromide tracer tests and HYDRUS-1D enhances reliability for estimation of GR and ET under film-mulched drip irrigation in the MRB and shows promise for other similar arid inland basins around the world.     

Comparison of the effects of litter covering and incorporation on infiltration and soil erosion under simulated rainfall

Plant litter can either cover on soil surface or be incorporated into top-soil layer in natural ecosystems. Their effects on infiltration and soil erosion are likely quite different. This study was performed to compare the effects of litter covering on soil surface and being incorporated into top-soil layer on infiltration and soil erosion under simulated rainfall. Four litter types (needle-leaf, broad-leaf, brush, and herb) were collected from fields and applied to cover on soil surface or to be incorporated into top-soil layer (5 cm) at the same rate (0.2 kg/m²). The simulated rainfalls (40 and 80 mm/hr) were run at two slope angles (10° and 20°). The results showed that the mean infiltration rate of litter covering treatment was 1.4 times as great as that of litter incorporated. Litter covering enhanced infiltration via protecting surface from soil sealing. Whereas, litter incorporation affected infiltration by its water repellency. Soil erosion of litter incorporated treatment was 5.4 times as large as that of litter covered treatment, which was attributed to the changes in surface litter coverage and soil erosion resistance. Litter type affected soil erosion through the variations in litter coverage and litter morphology. For litter covering treatment, litter coverage can explain the major variance of soil loss on the slopes. Whereas, for litter incorporated treatment, both the influences of litter coverage and litter length on soil erosion resistance were considered necessary to well explain the variance of soil loss. The results also showed that the benefits of litter to control soil erosion declined with rainfall intensity and slope gradient for both covering and incorporated treatments. The results of this study are helpful to understand the mechanisms of litter influencing hydrological and erosion processes on hillslopes.     

Responses of soil water dynamic processes and groundwater recharge to irrigation intensity and antecedent moisture in the vadose zone

In this study, a field experiment was conducted to investigate the soil water dynamics and water percolation through the deep vadose zone. A calibrated HYDRUS‐1D model was used to simulate the process of soil water movement and the water budget. Based on the measured volumetric soil water contents, the model was well calibrated and validated. Then, we conducted scenario analyses to determine the combined effects of irrigation amount (IA), antecedent soil moisture (AM), crop evapotranspiration, and deep percolation (DP) in an irrigation event. Four IAs (5, 10, 15, and 20 cm) and three AM conditions (AM‐1, AM‐2, and AM‐3) were controlled in the scenario analyses. The results indicate that according to the S_e's (effective saturation) values status and the observed or simulated depth, there could be different conclusions on the influence of DP. Under different IAs in dry (AM‐1) and medium (AM‐2) AM status, DP changed slightly; it was 0.39 and 2.47 cm in AM‐1 and 0.40 and 2.48 cm in AM‐2 for the summer maize and winter wheat crop, respectively; the AM had a crucial contribution to DP. While under the condition of wet AM (AM‐3) or small observation depth, the water inputs could have a significant effect on DP. According to increasing irrigation intensity, the higher values of S_e (>0.6) in the whole profile were only displayed between 70 and 300 cm at AM‐1, 70–500 cm at AM‐2, and 70‐below 600 cm at AM‐3, which were gradually extended and moved down with increasing AM. Hence, the IA significantly affected the water percolation at a depth of 200 cm, whereas there was a weak influence at 600 cm except in AM‐3. Furthermore, in the higher values of the S_e (>0.65) domain, the correlation between IA and DP was an exponential function and significantly under P < 0.05. In addition, DP began to occur when the soil water content was equal to or greater than 0.75 times that of the field water capacity or the S_e > 0.65. When the coarse silt layer became embedded in the silt clay soil profile, it lagged the process of water transport but did not affect permeability in the end.     

The effect of geometrical parameters on the discharge capacity of meandering compound channels

A number of methods and formulae has been proposed in the literature to estimate the discharge capacity of compound channels. When the main channel has a meandering pattern, a reduction in the conveyance capacity for a given stage is observed, which is due to the energy dissipations caused by the development of strong secondary currents and to the decrease of the main channel bed slope with respect to the valley bed slope. The discharges in meandering compound channels are usually assessed applying, with some adjustments, the same methods used in the straight compound channels. Specifically, the sinuosity of the main channel is frequently introduced to account for its meandering pattern, although some methods use different geometric parameters.In this paper the stage—discharge curves for several compound channels having identical cross-sectional area, roughness and bed slope but different planimetric patterns are numerically calculated and compared, in order to identify which geometric parameter should be efficaciously used in empirical formulae to account for meandering patterns. The simulations are carried out using a 3D finite-volume model that solves the RANS equations using a k-ε turbulence model. The numerical code is validated against experimental data collected in both straight and meandering compound channels.The numerical results show that the sinuosity is the main parameter to be accounted for in empirical formulae to assess the conveyance capacity of meandering compound channels. Comparison of the stage—discharge curves in the meandering compound channels with that obtained in a straight channel having identical cross-sectional area clearly shows the reduction of discharge due to the presence of bends in the main channel. The effect of other geometric parameters, such as the meander-belt width and the mean curvature radius, results very weak.     

OPTIMUM CHANNEL GEOMETRY AND MINIMUM ENERGY DISSIPATION RATE

1. INTRODUCTIONA natural river adjusts its roughness, channel geometry, slope I length, and Pattern to changing hydrologic, geologic, and manmade conditions. However, some sort of mean condition or quasi-equilibriumcondition appears to exist for natural streams. The earliest attempts to describe the quasi--equlibrium geometry of rivers were the empirical approaches, such as the regime approach proposed by Lacey (1929) andLeoPOld et al. (1953 ). Later, theoretical approaches based on the…     

REGIME MORPHOLOGY OF EQUILIBRIUM ALLUVIAL CHANNELS

IINTRODUCTIONThequestfordeterminingthedesigncharacteristicsofregimechannelshasbeengoingonforalongtime.Peoplehavebeenexcavatingnewormodifyingexistingchannelstousethemforirrigation,watersupply,navigation,floodcontrolandotherpurposes.Recently,archeologistsdiscoveredwhatiscurrentlybelievedtobetheoldestman(madecanalsystem.ItwasfoundintheareawhereMesopotamiausedtoexistanditisdatedbacktoabollt4,000BC.Ifachanne]isnotproperlydesigned,erosionofitsbanksordepositionofsedimentwithinitscross-sectionw…     

A study of soil water movement combining soil water potential with stable isotopes at two sites of shallow groundwater areas in the North China Plain

In the shallow groundwater areas of the North China Plain (NCP), precipitation infiltration and evapotranspiration in the vertical direction are the main processes of the water cycle, in which the unsaturated zone plays an important role in the transformation process between precipitation and groundwater. In this paper, two typical sites in Cangzhou (CZ) and Hengshui (HS) of Hebei province with shallow water tables were selected to analyse the relationship among precipitation, soil water and groundwater. At each site, precipitation, soil water at depths 10, 20, 30, 50, 70, 100, 150, 200, 300 cm, and groundwater were sampled to analyse the stable isotope compositions of hydrogen and oxygen. The soil water potentials at the corresponding depths were observed. Although the climates at the two sites are similar, there are some differences in the infiltration process, soil water movement and groundwater recharge sources. Evaporation occurred at the upper depths, which led to the decrease of soil potential and the enrichment of heavy isotopes. At the CZ site, precipitation infiltrated with piston mode, and an obvious mixture effect existed during the infiltration process. Preferential flow may exist in the soil above 100 cm depth. However, at the HS site soil water moved in piston mode, and groundwater was mainly recharged by precipitation. When precipitation recharged the groundwater it experienced a strong evaporation effect. The results of the soil water movement mechanism provides the transformation relationship among precipitation, soil water and groundwater in the middle and eastern NCP. Copyright © 2009 John Wiley & Sons, Ltd.     

Responses of soil moisture at different topographic positions to rainfall events along a steep subtropical forested hillslope

Understanding the dynamic response of soil moisture to rainfall is crucial for describing hydrological processes at the hillslope scale. However, because of sparse monitoring coupled with the complexity of water movement and steep topography, the findings of rainfall-related soil moisture dynamics have not always been consistent, indicating a need for systematic investigations of soil moisture dynamics and infiltration patterns following rainfall inputs at multiple topographic positions along a hillslope. This study aimed to examine the nature of these responses by characterizing and quantifying the response amplitude, rate and time for 37 large rainfall events at 25 combinations of topographic positions and soil depths along a steep forested hillslope. Our results showed that soil moisture responses under different rainfall patterns could be attributed to one or the other rainfall characteristics, such as rainfall intensity and amount. However, soil moisture dynamics at different hillslope positions after rainfall varied widely due to the controls of soil properties, topography, and non-equilibrium flow. Preferential flow was more evident under dry initial soil conditions than under wet initial soil conditions. Findings of this study reveal that the dynamic response patterns of soil moisture to rainfall do not always follow topographic controls, which can improve our understanding of water cycling related to the infiltration process at the hillslope scale, and support water resources management in subtropical mountain ecosystems.     

STUDY OF INFILTRATION INTO A HETEROGENEOUS SOIL USING MAGNETIC RESONANCE IMAGING

In this paper, the feasibility of using magnetic resonance imaging (MRI) to study water infiltration into a heterogeneous soil is examined, together with its difficulties and limitations. MRI studies of ponded water infiltration into an undisturbed soil core show that the combination of one- and two-dimensional imaging techniques provides a visual and non-destructive means of monitoring the temporal changes of soil water content and the moisture profile, and the movement of the wetting front. Two-dimensional images show air entrapment in repetitive ponded infiltration experiments. During the early stages of infiltration, one-dimensional images of soil moisture profiles clearly indicate preferential flow phenomena. The observed advance of wetting fronts can be described by a linear relationship between the square root of infiltration time (√t) and the distance of the wetting front from the soil surface. Similarly, the cumulative infiltration is also directly proportional to √t. Furthermore, from the MRI infiltration moisture profiles, it is possible to estimate the parameters that feature in infiltration equations. © 1997 by John Wiley & Sons, Ltd.     

Effects of frozen soil and snow cover on cold‐season soil water dynamics in Tokachi,Japan

Despite the potential impact of winter soil water movements in cold regions, relatively few field studies have investigated cold‐season hydrological processes that occur before spring‐onset of snowmelt infiltration. The contribution of soil water fluxes in winter to the annual water balance was evaluated over 5 years of field observations at an agricultural field in Tokachi, Hokkaido, Japan. In two of the winters, soil frost reached a maximum depth of 0·2 m (‘frozen’ winters), whereas soil frost was mostly absent during the remaining three winters (‘unfrozen’ winters). Significant infiltration of winter snowmelt water, to a depth exceeding 1·0 m, occurred during both frozen and unfrozen winters. Such infiltration ranged between 126 and 255 mm, representing 28–51% of total annual soil water fluxes. During frozen winters, a substantial quantity of water (ca 40 mm) was drawn from deeper layers into the 0–0·2 m topsoil layer when this froze. Under such conditions, the progression and regression of the freezing front, regulated by the thickness of snow cover, controlled the quantity of soil water flux below the frozen layer. During unfrozen winters, 13–62 mm of water infiltrated to a depth of 0·2 m, before the spring snowmelt. These results indicate the importance of correctly evaluating winter soil water movement in cold regions. Copyright © 2010 John Wiley & Sons, Ltd.